USGS Science⌕ Search

SEARCH · USGS Science

Results for “Near Surface Geophysics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47Linked to original sources

Thermal infrared spectroscopy of experimentally shocked anorthosite and pyroxenite: Implications for remote sensing of Mars

The feldspar and pyroxene mineralogies on Mars revealed by the Thermal Emission Spectrometer (TES) on Mars Global Surveyor likely record a variety of shock effects, as suggested by petrologic analyses of the Martian meteorites and the abundance of impact craters on the planet's surface. To study the effects of shock pressures on thermal infrared spectra of these minerals, we performed shock recovery experiments on orthopyroxenite and anorthosite samples from the Stillwater Complex (Montana) over peak pressures from 17 to 63 GPa. We acquired emissivity and hemispherical reflectance spectra (350–1400 cm −1 ; ∼7–29 μm) of both coherent chips and fine-grained powders of shocked and unshocked samples. These spectra are more directly comparable to remotely sensed data of Mars (e.g., TES) than previously acquired absorption or transmission spectra of shocked minerals. The spectra of experimentally shocked feldspar show systematic changes with increasing pressure due to depolymerization of the silica tetrahedra. For the spectra of chips, this includes the disappearance of small bands in the 500–650 cm −1 region and a strong band at 1115 cm −1 , and changes in positions of a strong band near 940 cm −1 and the Christiansen feature near 1250 cm −1 . Spectra of the shocked powders show the gradual disappearance of a transparency feature near 830 cm −1 . Fewer changes are observed in the pyroxene spectra at pressures as high as 63 GPa. Spectra of experimentally shocked minerals will help identify more precisely the mineralogy of rocks and soils not only from TES but also from Mars instruments such as miniTES and THEMIS.

Journal of Geophysical Research E: Planets↗

Rheologic and structural controls on the deformation of Okmok volcano, Alaska: FEMs, InSAR, and ambient noise tomography

Interferometric synthetic aperture radar (InSAR) data indicate that the caldera of Okmok volcano, Alaska, subsided more than a meter during its eruption in 1997. The large deformation suggests a relatively shallow magma reservoir beneath Okmok. Seismic tomography using ambient ocean noise reveals two low‐velocity zones (LVZs). The shallow LVZ corresponds to a region of weak, fluid‐saturated materials within the caldera and extends from the caldera surface to a depth of 2 km. The deep LVZ clearly indicates the presence of the magma reservoir beneath Okmok that is significantly deeper (>4 km depth) compared to previous geodetic‐based estimates (3 km depth). The deep LVZ associated with the magma reservoir suggests magma remains in a molten state between eruptions. We construct finite element models (FEMs) to simulate deformation caused by mass extraction from a magma reservoir that is surrounded by a viscoelastic rind of country rock embedded in an elastic domain that is partitioned to account for the weak caldera materials observed with tomography. This configuration allows us to reduce the estimated magma reservoir depressurization to within lithostatic constraints, while simultaneously maintaining the magnitude of deformation required to predict the InSAR data. More precisely, the InSAR data are best predicted by an FEM simulating a rind viscosity of 7.5 × 10 16 Pa s and a mass flux of −4.2 × 10 9 kg/d from the magma reservoir. The shallow weak layer within the caldera provides a coeruption stress regime and neutral buoyancy horizon that support lateral magma propagation from the central magma reservoir to extrusion near the rim of the caldera.

Alaska↗

Crustal structure of east central Oregon: Relation between Newberry Volcano and regional crustal structure

A 180-km-long seismic refraction transect from the eastern High Cascades, across Newberry Volcano, to the eastern High Lava Plains is used to investigate the subvolcanic crustal and upper mantle velocity structure there. Near-surface volcanic flows and sedimentary debris (1.6-4.7 km/s), ranging from 3 to 5 km in thickness, overlie subvolcanic Basin and Range structures. East and west of Newberry Volcano, the subvolcanic basement (5.6 km/s) has been downwarped, producing 5-km-deep basins. The midcrust (8- to 28-km depth) is characterized by velocities ranging from 6.1 to 6.5 km/s and varies laterally in thicknesses. The lower crust is characterized by an unusually high velocity (about 7.4 km/s), and its geometry mirrors the subvolcanic basement geometry. The Moho is located at a depth of 37 km and represents a transition to an upper mantle velocity of 8.1 km/s. The shallow subsurface (1.2 km) beneath Newberry Volcano is characterized by high-velocity 5.6 km/s, versus 4.1 km/s for the surrounding area) intrusions and appears to be located on a basement high. Beneath the seismic refraction array at Newberry Volcano, an absence of low-velocity anomalies suggests that large silicic magma chambers do not exist in the upper crust, but apparent high attenuation of the seismic wave field may be consistent with either partial melts in small volumes, elevated crustal temperatures, and/or poor geophone-recording site coupling.

Oregon↗

Assessment of Mars Pathfinder landing site predictions

Remote sensing data at scales of kilometers and an Earth analog were used to accurately predict the characteristics of the Mars Pathfinder landing site at a scale of meters. The surface surrounding the Mars Pathfinder lander in Ares Vallis appears consistent with orbital interpretations, namely, that it would be a rocky plain composed of materials deposited by catastrophic floods. The surface and observed maximum clast size appears similar to predictions based on an analogous surface of the Ephrata Fan in the Channeled Scabland of Washington state. The elevation of the site measured by relatively small footprint delay-Doppler radar is within 100 m of that determined by two-way ranging and Doppler tracking of the spacecraft. The nearly equal elevations of the Mars Pathfinder and Viking Lander 1 sites allowed a prediction of the atmospheric conditions with altitude (pressure, temperature, and winds) that were well within the entry, descent, and landing design margins. High-resolution (~38 m/pixel) Viking Orbiter 1 images showed a sparsely cratered surface with small knobs with relatively low slopes, consistent with observations of these features from the lander. Measured rock abundance is within 10% of that expected from Viking orbiter thermal observations and models. The fractional area covered by large, potentially hazardous rocks observed is similar to that estimated from model rock distributions based on data from the Viking landing sites, Earth analog sites, and total rock abundance. The bulk and fine-component thermal inertias measured from orbit are similar to those calculated from the observed rock size-frequency distribution. A simple radar echo model based on the reflectivity of the soil (estimated from its bulk density), and the measured fraction of area covered by rocks was used to approximate the quasi-specular and diffuse components of the Earth-based radar echos. Color and albedo orbiter data were used to predict the relatively dust free or unweathered surface around the Pathfinder lander compared to the Viking landing sites. Comparisons with the experiences of selecting the Viking landing sites demonstrate the enormous benefit the Viking data and its analyses and models had on the successful predictions of the Pathfinder site. The Pathfinder experience demonstrates that, in certain locations, geologic processes observed in orbiter data can be used to infer surface characteristics where those processes dominate over other processes affecting the Martian surface layer. Copyright 1999 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Mars Exploration Program 2007 Phoenix landing site selection and characteristics

To ensure a successful touchdown and subsequent surface operations, the Mars Exploration Program 2007 Phoenix Lander must land within 65° to 72° north latitude, at an elevation less than −3.5 km. The landing site must have relatively low wind velocities and rock and slope distributions similar to or more benign than those found at the Viking Lander 2 site. Also, the site must have a soil cover of at least several centimeters over ice or icy soil to meet science objectives of evaluating the environmental and habitability implications of past and current near‐polar environments. The most challenging aspects of site selection were the extensive rock fields associated with crater rims and ejecta deposits and the centers of polygons associated with patterned ground. An extensive acquisition campaign of Odyssey Thermal Emission Imaging Spectrometer predawn thermal IR images, together with ∼0.31 m/pixel Mars Reconnaissance Orbiter High Resolution Imaging Science Experiment images was implemented to find regions with acceptable rock populations and to support Monte Carlo landing simulations. The chosen site is located at 68.16° north latitude, 233.35° east longitude (areocentric), within a ∼50 km wide (N‐S) by ∼300 km long (E‐W) valley of relatively rock‐free plains. Surfaces within the eastern portion of the valley are differentially eroded ejecta deposits from the relatively recent ∼10‐km‐wide Heimdall crater and have fewer rocks than plains on the western portion of the valley. All surfaces exhibit polygonal ground, which is associated with fracture of icy soils, and are predicted to have only several centimeters of poorly sorted basaltic sand and dust over icy soil deposits.

Journal of Geophysical Research E: Planets↗

Relative age of interior layered deposits in southwest Candor Chasma based on high-resolution structural mapping

High‐resolution topography generated from stereo HiRISE (High‐Resolution Imaging Science Experiment) imagery reveals the meter‐scale structure of interior layered deposits (ILD) in southwest Candor Chasma. This study seeks to determine the age of the local ILD relative to any normal faults that can be attributed to chasma formation. The study area is located near the contact of these ILD and the wall rock and is in an area where chasma‐forming normal faults have been proposed. We find that while normal faults are found in the study area, these faults are not sufficiently large nor appropriately located or oriented to accommodate the roughly northeast‐southwest extension that is required for normal faults that can be attributed to chasma formation. Additionally, bedding exposed in the local ILD generally dips toward the center of Candor Chasma, consistent with sediment deposition in a preexisting basin. Further, pit craters of Tithonia Catena, presumed to predate or be contemporaneous with the formation of west Candor Chasma, do not cut into the ILD within the study area. These independent lines of evidence support a postchasma age for the ILD exposed within the study area. Chasma‐related normal faults may exist within these ILD at depth but are not exposed at the surface. Approximately 2 km of conformable stratigraphy is exposed in the study area, and therefore at least several kilometers of the local ILD were deposited subsequent to any chasma‐related normal faulting that may have occurred in this part of Candor Chasma.

Journal of Geophysical Research E: Planets↗

Remote compositional analysis of lunar olivine-rich lithologies with Moon Mineralogy Mapper (M3) spectra

A systematic approach for deconvolving remotely sensed lunar olivine‐rich visible to near‐infrared (VNIR) reflectance spectra with the Modified Gaussian Model (MGM) is evaluated with Chandrayaan‐1 Moon Mineralogy Mapper (M 3 ) spectra. Whereas earlier studies of laboratory reflectance spectra focused only on complications due to chromite inclusions in lunar olivines, we develop a systematic approach for addressing (through continuum removal) the prominent continuum slopes common to remotely sensed reflectance spectra of planetary surfaces. We have validated our continuum removal on a suite of laboratory reflectance spectra. Suites of olivine‐dominated reflectance spectra from a small crater near Mare Moscoviense, the Copernicus central peak, Aristarchus, and the crater Marius in the Marius Hills were analyzed. Spectral diversity was detected in visual evaluation of the spectra and was quantified using the MGM. The MGM‐derived band positions are used to estimate the olivine's composition in a relative sense. Spectra of olivines from Moscoviense exhibit diversity in their absorption features, and this diversity suggests some variation in olivine Fe/Mg content. Olivines from Copernicus are observed to be spectrally homogeneous and thus are predicted to be more compositionally homogeneous than those at Moscoviense but are of broadly similar composition to the Moscoviense olivines. Olivines from Aristarchus and Marius exhibit clear spectral differences from those at Moscoviense and Copernicus but also exhibit features that suggest contributions from other phases. If the various precautions discussed here are weighed carefully, the methods presented here can be used to make general predictions of absolute olivine composition (Fe/Mg content).

Journal of Geophysical Research E: Planets↗

Variations in fault slip and strain accumulation at Parkfield, California: Initial results using two-color geodimeter measurements 1984-1988

Repeated length measurements of several geodetic baselines near Parkfield, California, have revealed significant variations in the local rates of shallow fault slip and strain. This network of baselines, surveyed several times each week, straddles the San Andreas fault in the transition zone between the creeping section to the northwest and the locked section to the southeast. The length measurements, characterized by a precision approaching 0.1 ppm, reveal large fluctuations in the rates of baseline extension. Principal mode analysis of the length change data indicates that the two largest components of the signal are (1) secular extensions and contractions consistent with surface slip on the main strand of the San Andreas fault, and (2) a large seasonal oscillation with no obvious spatial coherence. On most of these baselines, the second component appears to be in phase with seasonal rainfall. When data from the baselines with the largest amplitude of the seasonal signal are excluded, the remaining data can be modeled in terms of both spatial and temporal variations in surface slip, variations in the components of the changes in uniform strain, and the possible displacement of the central monument in this radial network. In parameterizing this model, the spatial variation of slip beneath the near surface is reflected by changes in shear strain. Although the computed secular shear is highly dependent upon the specified parameterization of surface slip, the data are consistent with the hypothesis that slip at intermediate depths lags behind the surface slip rate. However, the range in models that fit the data does not necessarily imply that there is a deficit in slip at depth relative to the surface. Comparison of the inferred values of surface slip from the model with the observed fault slip measured by very short baseline creep meters indicates close agreement in secular rates, but the short-term variations observed with the creep meters are either highly attenuated or nonexistent in the modeled slip since the modeled slip is a spatial average which smooths out possible short-wavelength variations in the surface slip for which the creep instruments are most sensitive. An interesting conclusion from the two-color data is that surface slip on the San Andreas fault appears to be spread over a 2-km-wide zone on the south flank of Middle Mountain but is confined to a very narrow zone to the south as the fault passes through the center of the network. This conclusion is dependent upon the assumption that a few critical monuments are stable and track tectonic displacements in the long term. Finally, the largest observed strain change is an extensional strain coincident with the Kettleman Hills earthquake M 5.5 in August 1985.

Journal of Geophysical Research Solid Earth↗

Rupture process of a multiple main shock sequence: analysis of teleseismic, local and field observations of the Tennant Creek, Australia, earthquakes of January 22, 1988

On January 22, 1988, three large intraplate earthquakes (with M S 6.3, 6.4, and 6.7) occurred within a 12-hour period near Tennant Creek, Australia. These earthquakes, which occurred over a small interval of time and within a small volume of space, present a unique opportunity to study the rupture process of the class of intraplate earthquakes that occur as multiple main shocks. Broadband displacement and velocity records of body waves from teleseismically recorded data are analyzed to determine source mechanisms, depths, and complexity of rupture of each of the three main shocks. Hypocenters of an additional 150 foreshocks and aftershocks constrained by local arrival time data and field observations of surface rupture are used to complement the source characteristics of the main shocks in order to derive as complete a description of the rupture process as possible. The interpretation of the combined data sets suggests that the overall rupture process involved unusually complicated stress release. As locations of the main shock hypocenters progressively moved from west to east, we infer that the first and third main shocks, denoted as MS 1 and MS 3 , produced the southeast-northwest trending scarps observed at the western end (the Kunayungku fault) and at the eastern end (the east end of the Lake Surprise fault), respectively, of the rupture zone. The epicenter of the only immediate foreshock was located in the gap between these two fault scarps. MS 1 nucleated near this epicenter and ruptured upward and to the northwest from a depth of 6.5 km. MS 3 ruptured predominantly to the SE at a depth of 4.5 km. The second main shock, MS 2 , is inferred to have produced the deformation of the southwest trending central scarp segment (the western end of the Lake Surprise fault). From the sense of thrusting seen at the surface and from the distribution of aftershock hypocenters, the south dipping nodal planes derived from waveform modeling are identified as the fault planes for earthquakes MS1 and MS 3 . In contrast, the dip of the central fault scarp is reversed relative to the dips of the western and eastern fault scarps. The rupture process Of MS 2 turns out to be commensurately complex and sufficiently explains the geological complexity. MS 2 consisted of three subevents. The southeast dipping nodal plane of the first two subevents is coplanar with a southeast dipping plane implied by locations of aftershocks which did not break the surface. Choice of the north dipping plane as the rupture plane of the third subevent, consistent with the surface deformation and coplanar with a second plane delineated by aftershocks, would imply conjugate faulting. The majority of the aftershocks are concentrated near the edges of the fault planes, and there is an absence of activity in the center of the planes. The areas of absent activity may represent the failed asperities of the main shocks in which substantial stress relief occurred. The rupture process of each main shock is characterized by the rapid release of energy followed by a much slower release of moment and by aftershock zones whose dimensions exceed the inferred dimensions of the rupture. These characteristics suggest that substantial slow slip occurred on each of the three fault interfaces that was not accompanied by major energy release. The first main shock nucleated at the deepest part of its aftershock zone and ruptured upward. In contrast, MS 2 and MS 3 nucleated at depths that were one third to one half of the maximum depth of their aftershock zones. This variation of focal depth and the strong increase of moment and radiated energy with each main shock imply that lateral variations of strength were more important than vertical gradients of shear stress in controlling the progression of rupture.

Journal of Geophysical Research Solid Earth↗

Lunar prospector epithermal neutrons from impact craters and landing sites: Implications for surface maturity and hydrogen distribution

Initial studies of neutron spectrometer data returned by Lunar Prospector concentrated on the discovery of enhanced hydrogen abundances near both lunar poles. However, the nonpolar data exhibit intriguing patterns that appear spatially correlated with surface features such as young impact craters (e.g., Tycho). Such immature crater materials may have low hydrogen contents because of their relative lack of exposure to solar wind-implanted volatiles. We tested this hypothesis by comparing epithermal* neutron counts (i.e., epithermal −0.057 × thermal neutrons) for Copernican-age craters classified as relatively young, intermediate, and old (as determined by previous studies of Clementine optical maturity variations). The epithermal* counts of the crater and continuous ejecta regions suggest that the youngest impact materials are relatively devoid of hydrogen in the upper 1 m of regolith. We also show that the mean hydrogen contents measured in Apollo and Luna landing site samples are only moderately well correlated to the epithermal* neutron counts at the landing sites, likely owing to the effects of rare earth elements. These results suggest that further work is required to define better how hydrogen distribution can be revealed by epithermal neutrons in order to understand more fully the nature and sources (e.g., solar wind, meteorite impacts) of volatiles in the lunar regolith.

Journal of Geophysical Research E: Planets↗

Geophysical well log analysis of fractured granitic rocks at Atikokan, Ontario, Canada

Two boreholes, drilled to approximate depths of 750 and 1,260 m in a granitic intrusion located near Atikokan, Ontario, were studied by obtaining a full suite of conventional borehole geophysical logs. In addition, selected intervals in these boreholes were logged with a borehole acoustic televiewer that produces a high-resolution image of the borehole wall, an acoustic waveform-logging system using 34-kiloHertz magnetostrictive and 5-kiloHertz sparker sources, and a highly sensitive heat-pulse flowmeter. Emphasis was on identifying and characterizing fracture zones that represent groundwater conduits in deeper portions of the granite, and on characterizing the properties of the largest intervals of unfractured granite. Major fracture zones were indicated by correlating geophysical log anomalies detected on the suite of conventional logs (unpublished data from Atomic Energy of Canada). However, several other anomalies, were identified as mafic intrusions of approximately the same thickness as major fracture zones. Geophysical log anomalies were compared for all major fracture zones that could serve as significant groundwater conduits, and fracture zone permeability is estimated on the basis of acoustic tube-wave attenuation in these intervals. Acoustic televiewer logs obtained at depths below 1,000 m in the deeper well indicate that most of the few fractures identified on core at these depths do not remain open enough under in situ conditions to produce detectable anomalies in acoustic refraction. Flowmeter data indicate that some groundwater circulation occurs in the upper portion of both boreholes. Water in the shallower of the two holes was observed to flow at 2.0 L/min; most of this flow entered the borehole at a depth < 25 m, and no flow occurred below a depth of 100 m. Downflow at rates < 0.5 L/min was determined to enter the deeper borehole within 20 m of the surface, and to exist at various fractures down to a depth of 250 m. (Author 's abstract)

Water-Resources Investigations Report↗

Implications of low-temperature cooling history on a transect across the Colorado Plateau-Basin and Range boundary, west central Arizona

Fission track ages of apatite and zircon from metamorphic, plutonic, and sedimentary rocks along a 80-km transect across the Colorado Plateau-Basin and Range boundary in west central Arizona show differences in the low-temperature cooling histories between the provinces. The transect extends from Cypress Mountain in the Colorado Plateau transition zone to the eastern Buckskin Mountains in the Basin and Range. Along the northeast margin of the Basin and Range province, metamorphic rocks exposed in the footwall of a major detachment fault system yield zircon and apatite fission track ages of 16–10 Ma. These ages are similar to K-Ar fusion ages of biotite and age minima of K-feldspar on 40 Ar/ 39 Ar age spectra and collectively indicate rapid cooling. One K-feldspar age spectrum has an age maximum of about 22.8 Ma, an age minimum of 11.8 Ma, and a spectrum whose shape is suggestive of reheating, possibly in middle Miocene time. The heating event was probably related to hydrothermal activity during emplacement of Cu and Mn deposits in and above the detachment fault zone. Effects of this heating are only locally detected in rocks above the detachment fault. In the Poachie Range fission track ages of apatite and zircon are 60–50 and 80–70 Ma, respectively. The disparity between the apatite and zircon ages indicates that the rocks cooled slowly in Late Cretaceous and early Tertiary time, probably due to gradual uplift and erosion. Total uplift and denudation in the area of the Poachie Range since Cretaceous time is 6 km or more. To the northeast of the range, fission track ages of apatite and zircon increase and diverge, indicating that apparent uplift decreased in that direction. The apatite ages from the Poachie Range are concordant with early Tertiary hornblende ages determined in other studies of lower plate rocks near the southwest end of the transect. The ages represent cooling of crystalline rocks after Cretaceous regional metamorphism and magmatism. Near Bagdad, 20 km northeast of the Poachie Range, 2 km or less of erosion has occurred since intrusion of high-level plutons and dikes and caldera formation in Late Cretaceous time. Remnants of an erosion surface that developed in middle Tertiary time are preserved in the transition zone. Volcanic and sedimentary rocks at least as old as early Miocene were deposited on the erosion surface and filled valleys cut into it. Dissection of these deposits began about 8 Ma. We interpret these data combined with those from other studies to indicate that in Cretaceous time southward thrusting and later extensive magmatism in the middle crust led to thickening and heating of the crust. The Cretaceous igneous rocks at Bagdad are high-level manifestations of this magmatism. Uplift and slow cooling occurred in Late Cretaceous and early Tertiary time. In late Oligocene and early Miocene time during northeast-southwest extension, middle crustal rocks moved southwest put from beneath the southwest margin of the transition zone. Tectonic denudation rapidly exposed the crust that had been brought up from a depth of 10 km or more and rapidly cooled in the eastern Buckskin and Harcuvar mountains. Middle Miocene reheating occurred locally in the lower plate, along the detachment, and in nearby parts of the upper plate.

Journal of Geophysical Research↗

Reconstructing the dynamics of the highly similar May 2016 and June 2019 Iliamna Volcano, Alaska ice–rock avalanches from seismoacoustic data

Surficial mass wasting events are a hazard worldwide. Seismic and acoustic signals from these often remote processes, combined with other geophysical observations, can provide key information for monitoring and rapid response efforts and enhance our understanding of event dynamics. Here, we present seismoacoustic data and analyses for two very large ice–rock avalanches occurring on Iliamna Volcano, Alaska (USA), on 22 May 2016 and 21 June 2019. Iliamna is a glacier-mantled stratovolcano located in the Cook Inlet, ∼200 km from Anchorage, Alaska. The volcano experiences massive, quasi-annual slope failures due to glacial instabilities and hydrothermal alteration of volcanic rocks near its summit. The May 2016 and June 2019 avalanches were particularly large and generated energetic seismic and infrasound signals which were recorded at numerous stations at ranges from ∼9 to over 600 km. Both avalanches initiated in the same location near the head of Iliamna's east-facing Red Glacier, and their ∼8 km long runout shapes are nearly identical. This repeatability – which is rare for large and rapid mass movements – provides an excellent opportunity for comparison and validation of seismoacoustic source characteristics. For both events, we invert long-period (15–80 s) seismic signals to obtain a force-time representation of the source. We model the avalanche as a sliding block which exerts a spatially static point force on the Earth. We use this force-time function to derive constraints on avalanche acceleration, velocity, and directionality, which are compatible with satellite imagery and observed terrain features. Our inversion results suggest that the avalanches reached speeds exceeding 70 m s−1, consistent with numerical modeling from previous Iliamna studies. We lack sufficient local infrasound data to test an acoustic source model for these processes. However, the acoustic data suggest that infrasound from these avalanches is produced after the mass movement regime transitions from cohesive block-type failure to granular and turbulent flow – little to no infrasound is generated by the initial failure. At Iliamna, synthesis of advanced numerical flow models and more detailed ground observations combined with increased geophysical station coverage could yield significant gains in our understanding of these events.

Alaska↗

Ancient winds, waves, and atmosphere in Gale Crater, Mars, inferred from sedimentary structures and wave modeling

Wave modeling and analysis of sedimentary structures were used to evaluate whether four examples of symmetrical, reversing, or straight-crested bedforms in Gale crater sandstones are preserved wave ripples; deposition by waves would demonstrate that the lake was not covered by ice at that time. Wave modeling indicates that regardless of atmospheric density, winds that exceeded the threshold of aeolian sand transport could have generated waves capable of producing nearshore wave ripples in most grain sizes of sand. Reversing 3-m-wavelength bedforms in the Kimberley formation are interpreted not as wave ripples but rather as large aeolian ripples that formed in an atmosphere approximately as thin as at present. These exhumed bedforms define many of the ridges at outcrops that appear striated in satellite images. At Kimberley these bedforms demonstrably underlie and therefore predate subaqueous beds, suggesting that a thin atmosphere existed at least temporarily before subaqueous deposition ceased in the crater. The other three candidate wave ripples (Square Top, Hunda, and Voe) are consistent with modeled waves, but other origins cannot be excluded. The predominance of flat-laminated (non-rippled) beds in the lacustrine Murray formation suggests that some aspect of the lake was not conducive to formation or preservation of recognizable wave ripples. Water depths may generally have been too deep, lakebed sediment may have been too fine-grained, the lake may have been smaller than modeled, or the lake may have been covered by ice. Plain Language Summary Wave modeling and analysis of sedimentary structures were used to evaluate whether ancient lake deposits in Gale crater contain ripples formed by waves on the surface of the lake. Deposition by waves would show that the lake was not covered by ice at that time. Modeling shows that regardless of atmospheric density, winds capable of moving sand on land would generally have been strong enough to form waves that would produce ripples near shore. Large bedforms in the Kimberley formation are interpreted as ripples formed by the wind in an atmosphere similar to that of Mars today. These bedforms underlie and are older than other beds deposited in water, thereby showing that a thin atmosphere existed at least temporarily before deposition in water ceased in the crater. Three other candidate wave ripples are consistent with modeled waves, but other origins are possible. Thick sequences of sedimentary rock in Gale crater are flat-laminated rather than rippled, suggesting that some aspect of the lake was not favorable for their formation or preservation. Much of the lake may have been too deep or ice-covered, or the lake may have been smaller than modeled or had sediment too fine to form easily observed ripples.

Journal of Geophysical Research: Planets↗

Internal structure of the Sierra Nevada batholith based on specific gravity and gravity measurements

About 6,000 specific‐gravity (SG) measurements of samples collected from nearly 200 granitic plutons comprising the central Sierra Nevada batholith yield a SG contour map across the batholith from 36.25° to 38° north latitude. With notable exceptions, SG decreases from values generally greater than 2.7 in the west to less than 2.6 over a few small areas of high‐silica, high‐potassium granites near the east edge. A good correlation between measured SG and analyzed weight percent SiO 2 enables estimation of average silica variations across the batholith. The average SG is 2.69 corresponding to an average of 68 wt. % SiO 2 for the 18,000 km² central part of the batholith. A 1‐km gridded version of the SG measurements has been used to generate a series of synthetic gravity maps, assuming that the SG of rocks at the surface extends unchanged to various depths. The synthetic map computed for a depth of 10 km shows the best correspondence with the isostatic residual gravity map indicating that variations in the observed gravity residuals are largely caused by SG variations of the plutonic rocks exposed at the surface that apparently extend downward to an average depth of about 10 km. Although the 10‐km synthetic gravity map gives the best overall fit to the observed gravity data, comparison of individual anomalies indicates that the bottoms of the plutons as defined by SG variations at the surface are generally shallower along the west edge of the Sierra Nevada (7±2 km) and deeper in the younger and more felsic eastern part (12±3 km). These depths do not necessarily represent a distinct base of the Sierra Nevada batholith. They may indicate the depth below which density homogenization occurs, either by igneous, or possibly, structural processes.

California↗

Imaging the mantle beneath Iceland using integrated seismological techniques

Using a combination of body wave and surface wave data sets to reveal the mantle plume and plume head, this study presents a tomographic image of the mantle structure beneath Iceland to 400 km depth. Data comes primarily from the PASSCAL-HOTSPOT deployment of 30 broadband instruments over a period of 2 years, and is supplemented by data from the SIL and ICEMELT networks. Three sets of relative teleseismic body wave arrival times are generated through cross correlation: S and SKS arrivals at 0.03–0.1 Hz, and P and PKIKP arrivals at 0.03–0.1 and 0.8–2.0 Hz. Prior to inversion the crustal portion of the travel time anomalies is removed using the crustal model ICECRTb. This step has a significant effect on the mantle velocity variations imaged down to a depth of ∼250 km. Inversion of relative arrival times only provides information on lateral velocity variations. Surface waves are therefore used to provide absolute velocity information for the uppermost mantle beneath Iceland. The average wave number for the Love wave fundamental mode at 0.020 and 0.024 Hz is measured and used to invert for the average S velocity. Combination of the body wave and surface wave information reveals a predominantly horizontal low-velocity anomaly extending from the Moho down to ∼250 km depth, interpreted as a plume head. Below the plume head a near-cylindrical low-velocity anomaly with a radius of ∼100 km and peak V P and V S anomalies of −2% and −4%, respectively, extends down to the maximum depth of resolution at 400 km. Within the plume head, in the uppermost mantle above the core of the plume, there is a relatively high velocity with a maximum V P and V S anomaly of +2%. This high-velocity anomaly may be the result of the extreme degree of melt extraction necessary to generate the thick (46 km) crust in central Iceland. Comparison of the plume volumetric flux implied by our images, the crustal generation rate, and the degree of melting suggested by rare earth element inversions, suggests that (1) mantle material must be flowing horizontally away from the plume core faster than the overlying lithosphere and (2) the bulk of the plume material does not participate in melting beneath Iceland.

Journal of Geophysical Research B: Solid Earth↗

A high-resolution seismic reflection/refraction study of the Chugach-Peninsular terrane boundary, southern Alaska

We present results from a high-resolution seismic refraction analysis of the shallow (approximately 2 km) crustal structure along the 107-km-long Trans-Alaska Crustal Transect Chugach reflection line in southern Alaska and a comparison with laboratory measurements of field samples. The refraction analysis includes the two-dimensional interpretation of several thousand first- and secondary-arrival travel times digitized from 1024-channel split-spread common shot gathers. The velocity model derived from this analysis better defines the location and geometry of terrane boundaries than does the normal incidence reflection section and agrees well with surface mapping of lithologies. Furthermore, the model predicts travel times within 100 ms of the reflection times recorded from the base of the Quaternary on the Chugach reflection section. Thicknesses of Quaternary deposits, with velocities between 1.1 and 2.0 km/s, correlate inversely with the quantity of observed lower crustal reflections on the Chugach section, suggesting that the presence or absence of these sediments in sufficient thickness exerted primary control on the quality of the deeper portion of the section. There is a significant velocity contrast between crystalline rocks across the Border Ranges fault (5.0 versus 5.6 km/s), the major contact between the Chugach and Peninsular terranes, in agreement with laboratory measurements of field specimens. In the Peninsular terrane the modeling indicates that an unnamed fault delimiting the southern flank of the Copper River Basin dips steeply northward at 50° and has about 1300 m of vertical offset. Laboratory measurements document a maximum velocity anisotropy of 20% for phyllitic schists of the Valdez Group in the Chugach terrane. In agreement with the observed E-W strike and near-vertical dip of the Valdez Group, we determined a significant (14%) velocity anisotropy for ray paths oriented N-S versus NE-SW.

Journal of Geophysical Research Solid Earth↗

Geology and water resources of Owens Valley, California

Owens Valley, a long, narrow valley located along the east flank of the Sierra Nevada in east-central California, is the main source of water for the city of Los Angeles. The city diverts most of the surface water in the valley into the Owens River-Los Angeles Aqueduct system, which transports the water more than 200 miles south to areas of distribution and use. Additionally, ground water is pumped or flows from wells to supplement the surface-water diversions to the river-aqueduct system. Pumpage from wells needed to supplement water export has increased since 1970, when a second aqueduct was put into service, and local concerns have been expressed that the increased pumpage may have had a detrimental effect on the environment and the indigenous alkaline scrub and meadow plant communities in the valley. The scrub and meadow communities depend on soil moisture derived from precipitation and the unconfined part of a multilayered aquifer system. This report, which describes the hydrogeology of the aquifer system and the water resources of the valley, is one in a series designed to (1) evaluate the effects that ground-water pumping has on scrub and meadow communities and (2) appraise alternative strategies to mitigate any adverse effects caused by pumping. Two principal topographic features are the surface expression of the geologic framework--the high, prominent mountains on the east and west sides of the valley and the long, narrow intermountain valley floor. The mountains are composed of sedimentary, granitic, and metamorphic rocks, mantled in part by volcanic rocks as well as by glacial, talus, and fluvial deposits. The valley floor is underlain by valley fill that consists of unconsolidated to moderately consolidated alluvial fan, glacial and talus, and fluvial and lacustrine deposits. The valley fill also includes interlayered recent volcanic flows and pyroclastic rocks. The bedrock surface beneath the valley fill is a narrow, steep-sided graben that is structurally separated into the Bishop Basin to the north and the Owens Lake Basin to the south. These two structural basins are separated by (1) a bedrock high that is the upper bedrock block of an east-west normal fault, (2) a horst block of bedrock (the Poverty Hills), and (3) Quaternary basalt flows and cinder cones that intercalate and intrude the sedimentary deposits of the valley fill. The resulting structural separation of the basins allowed separate development of fluvial and lacustrine depositional systems in each basin. Nearly all the ground water in Owens Valley flows through and is stored in the saturated valley fill. The bedrock, which surrounds and underlies the valley fill, is virtually impermeable. Three hydrogeologic units compose the valley-fill aquifer system, a defined subdivision of the ground-water system, and a fourth represents the valley fill below the aquifer system and above the bedrock. The aquifer system is divided into horizontal hydrogeologic units on the basis of either (1) uniform hydrologic characteristics of a specific lithologic layer or (2) distribution of the vertical hydraulic head. Hydrogeologic unit 1 is the upper unit and represents the unconfined part of the system, hydrogeologic unit 2 represents the confining unit (or units), and hydrogeologic unit 3 represents the confined part of the aquifer system. Hydrogeologic unit 4 represents the deep part of the ground-water system and lies below the aquifer system. Hydrogeologic unit 4 transmits or stores much less water than hydrogeologic unit 3 and represents either a moderately consolidated valley fill or a geologic unit in the valley fill defined on the basis of geophysical data. Nearly all the recharge to the aquifer system is from infiltration of runoff from snowmelt and rainfall on the Sierra Nevada. In contrast, little recharge occurs to the system by runoff from the White and Inyo Mountains or from direct precipitation on the valley floor. Ground water flows from the margins of the valley towards the center of the valley; the ground water then flows south to the terminus of the system at Owens (dry) Lake. Ground water flows south from Bishop Basin to Owens Lake Basin through the narrows that constrict the flow opposite Poverty Hills. The aquifer system in the northern half of Owens Lake Basin is divided into east and west halves by the barrier effect caused by the Owens Valley fault. Discharge from the aquifer system is primarily by pumpage and evapotranspiration, and to a lesser extent by flowing wells, springs, underflow, and leakage to the Owens River-Los Angeles aqueduct system. Withdrawals from pumped or flowing wells is the largest component of discharge and accounts for about 50 percent of the outflow from the system. Transpiration by scrub and meadow plant communities, and to a lesser extent by irrigated alfalfa pasture, accounts for about 40 percent of the system's discharge. Natural hydraulic conductivity ranges from less than 400 to about 12,000 feet per day in the basalt flows, the more permeable material in the aquifer system. Where the basalts are fractured by explosives and drilling techniques, actual transmissivities can be greater than 1,000,000 feet squared per day. Hydraulic conductivities in sedimentary deposits of the aquifer system range from less than a few feet per day in lacustrine clays to more than 300 feet per day in gravel stringers and beach deposits in the transition zone between alluvial fan deposits and fluvial and lacustrine deposits. Degree of confinement in the aquifer system generally increases to the south and east in both the Bishop and Owens Lake Basins. The vertical hydraulic gradient across hydrogeologic unit 2 and confining beds in hydrogeologic units 1 and 3 is a function of (1) the asymmetric recharge and hydraulic head created by the dominant recharge from Sierra Nevada runoff and (2) the areal extent and thickness of the confining beds. Although most of the pumpage is from hydrogeologic unit 3, some coincident drawdown has been recorded in nonpumped wells that tap unit 1. Drawdown in hydrogeologic unit 1 is a function of changes in (1) lateral flow through hydrogeologic unit 1, (2) upward flow of ground water through the confining beds, (3) downward leakage of water from hydrogeologic unit 1 to unit 3 through wells, (4) direct withdrawal from well intervals open to hydrogeologic unit 1, and (5) increased evapotranspiration. The water in the aquifer system is generally of excellent quality for public supply and irrigation, with the exception of water stored in thick sequences of lacustrine silts and clays near Owens Lake. The water is principally a calcium bicarbonate type and dissolved-solids concentrations range from approximately 104 to 325 milligrams per liter. Water in the lacustrine sediments of Owens (dry) Lake is a sodium bicarbonate type and dissolved-solids concentrations are about 5,400 milligrams per liter.

California↗