USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research- Earth Surface”

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 361 records · Page 20Linked to original sources
USGS2000

Crustal deformation associated with glacial fluctuations in the eastern Chugach Mountains, Alaska

The changes of the solid Earth in south central Alaska in response to two major glacial fluctuations on different temporal and spatial scales have been estimated and we evaluated their influence on the stress state and ongoing tectonic deformation of the region. During the recent (1993–1995) Bering Glacier surge, a large transfer of ice from the Bagley Ice Field to the Bering Glacier terminus region occurred. We estimated the elastic displacement of the solid Earth due to ice mass redistribution from Global Positioning System (GPS) measurements at sites near the surging glacier. We can account for these displacements by transfer of an ice volume of about 14 km 3 from the surge reservoir area to the terminus region. We examined the background seismicity (M L ≥ 2.5) before, during, and after the surge. We found that the occurrence of small earthquakes (M L ≤ 4.0) in the surge reservoir region increased during the surge time interval possibly in response to a decrease in ice mass. This suggests that a small decrease in the vertical stress, σ 3 , could be enough to modulate the occurrence of small, shallow earthquakes in this dominantly thrust fault setting. During this century the southern Alaska coastal glaciers have been undergoing an overall decrease in volume. Based on our compilation of changes in the extent and thickness of the coastal glaciers between the Malaspina and Bering, we calculated surface displacements due to the Earth's viscoelastic response to annual thinning and to the cumulative retreat over the last 100 years. The uplift of the region due to an average annual thinning rate of 1–6 m/yr in the ablation region is 1–12 mm/yr. For our reference model with a viscosity of 5×10 19 Pa s for depths between ≈ 40 and 200 km the total viscoelastic response due to the retreat over the last century may be as much as a couple of meters within the coastal ablation zone near Icy Bay. The maximum decrease in σv between 0 and 10 km was ≈ 1.0 MPa, which is significant in relation to the stress drops in recent earthquakes (≈ 2 to 10 MPa) but small in relation to the estimated tectonic stress magnitude. Therefore the occurrence of an earthquake such as the St. Elias (1979, M S = 7.2) may have been advanced in time; however, most of the ongoing stress accumulation would be primarily due to tectonic forces.

Journal of Geophysical Research B: Solid Earth↗
USGS2014

Eruptions at Lone Star geyser, Yellowstone National Park, USA: 2. Constraints on subsurface dynamics

We use seismic, tilt, lidar, thermal, and gravity data from 32 consecutive eruption cycles of Lone Star geyser in Yellowstone National Park to identify key subsurface processes throughout the geyser's eruption cycle. Previously, we described measurements and analyses associated with the geyser's erupting jet dynamics. Here we show that seismicity is dominated by hydrothermal tremor (~5–40 Hz) attributed to the nucleation and/or collapse of vapor bubbles. Water discharge during eruption preplay triggers high-amplitude tremor pulses from a back azimuth aligned with the geyser cone, but during the rest of the eruption cycle it is shifted to the east-northeast. Moreover, ~4 min period ground surface displacements recur every 26 ± 8 min and are uncorrelated with the eruption cycle. Based on these observations, we conclude that (1) the dynamical behavior of the geyser is controlled by the thermo-mechanical coupling between the geyser conduit and a laterally offset reservoir periodically filled with a highly compressible two-phase mixture, (2) liquid and steam slugs periodically ascend into the shallow crust near the geyser system inducing detectable deformation, (3) eruptions occur when the pressure decrease associated with overflow from geyser conduit during preplay triggers an unstable feedback between vapor generation (cavitation) and mass discharge, and (4) flow choking at a constriction in the conduit arrests the runaway process and increases the saturated vapor pressure in the reservoir by a factor of ~10 during eruptions.

Wyoming↗
USGS2015

Holocene geologic slip rate for the Banning strand of the southern San Andreas Fault, southern California

Northwest directed slip from the southern San Andreas Fault is transferred to the Mission Creek, Banning, and Garnet Hill fault strands in the northwestern Coachella Valley. How slip is partitioned between these three faults is critical to southern California seismic hazard estimates but is poorly understood. In this paper, we report the first slip rate measured for the Banning fault strand. We constrain the depositional age of an alluvial fan offset 25 ± 5 m from its source by the Banning strand to between 5.1 ± 0.4 ka (95% confidence interval (CI)) and 6.4 + 3.7/−2.1 ka (95% CI) using U-series dating of pedogenic carbonate clast coatings and 10 Be cosmogenic nuclide exposure dating of surface clasts. We calculate a Holocene geologic slip rate for the Banning strand of 3.9 + 2.3/−1.6 mm/yr (median, 95% CI) to 4.9 + 1.0/−0.9 mm/yr (median, 95% CI). This rate represents only 25–35% of the total slip accommodated by this section of the southern San Andreas Fault, suggesting a model in which slip is less concentrated on the Banning strand than previously thought. In rejecting the possibility that the Banning strand is the dominant structure, our results highlight an even greater need for slip rate and paleoseismic measurements along faults in the northwestern Coachella Valley in order to test the validity of current earthquake hazard models. In addition, our comparison of ages measured with U-series and 10 Be exposure dating demonstrates the importance of using multiple geochronometers when estimating the depositional age of alluvial landforms.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Acceleration spectra for subduction zone earthquakes

We estimate the source spectra of shallow earthquakes from digital recordings of teleseismic P wave groups, that is, P + p P + s P , by making frequency dependent corrections for the attenuation and for the interference of the free surface. The correction for the interference of the free surface assumes that the earthquake radiates energy from a range of depths. We apply this spectral analysis to a set of 12 subduction zone earthquakes which range in size from M S = 6.2 to 8.1, obtaining corrected P wave acceleration spectra on the frequency band from 0.01 to 2.0 Hz. Seismic moment estimates from surface waves and normal modes are used to extend these P wave spectra to the frequency band from 0.001 to 0.01 Hz. The acceleration spectra of moderate subduction zone earthquakes, that is, earthquakes whose seismic moments are less than 10 27 dyn cm, exhibit ω-square or Brune-type spectra, while the acceleration spectra of large subduction zone earthquakes, that is, earthquakes whose seismic moments are greater than 10 27 dyn cm, exhibit intermediate slopes where ü(ω) ∝ ω 5/4 for frequencies from 0.005 to 0.05 Hz. For this set of earthquakes, spectral shape appears to be a discontinuous function of seismic moment. Using reasonable assumptions for the phase characteristics, we transform the spectral shape observed for large earthquakes into the time domain to fit Ekström's (1987) moment rate functions for the M S = 8.1 Michoacan earthquake of September 19, 1985, and the M S = 7.6 Michoacan aftershock of September 21, 1985.

Journal of Geophysical Research Solid Earth↗
USGS2009

Moon meteoritic seismic hum: Steady state prediction

We use three different statistical models describing the frequency of meteoroid impacts on Earth to estimate the seismic background noise due to impacts on the lunar surface. Because of diffraction, seismic events on the Moon are typically characterized by long codas, lasting 1 h or more. We find that the small but frequent impacts generate seismic signals whose codas overlap in time, resulting in a permanent seismic noise that we term the "lunar hum" by analogy with the Earth's continuous seismic background seismic hum. We find that the Apollo era impact detection rates and amplitudes are well explained by a model that parameterizes (1) the net seismic impulse due to the impactor and resulting ejecta and (2) the effects of diffraction and attenuation. The formulation permits the calculation of a composite waveform at any point on the Moon due to simulated impacts at any epicentral distance. The root-mean-square amplitude of this waveform yields a background noise level that is about 100 times lower than the resolution of the Apollo long-period seismometers. At 2 s periods, this noise level is more than 1000 times lower than the low noise model prediction for Earth's microseismic noise. Sufficiently sensitive seismometers will allow the future detection of several impacts per day at body wave frequencies. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2001

Vertical tectonics in northern Escanaba Trough as recorded by thick late Quaternary turbidites

Escanaba Trough, the southernmost segment of the Gorda Ridge, is filled by as much as 500 m of late Quaternary turbidite and hemipelagic sediment. Coring at Deep Sea Drilling Project Site 35 and Ocean Drilling Program (ODP) Sites 1037 and 1038 together with 4.5-kHz deep-tow and 3.5-kHz surface-ship seismic reflection profiles enable a distinct pattern of reflections to be mapped throughout Escanaba Trough in the upper part of this sediment fill. The uppermost 80 m of turbidite sediment, which includes at least 11 turbidity current events, were deposited in <7 kyr. Nine of these turbidites are found throughout Escanaba Trough at water depths >3200 m. The turbidity currents were trapped upon entering Escanaba Trough, resulting in all of the sediment in suspension in the flows being deposited. The thickness of the turbidite layers reflects both the flow thickness and the vertical grain concentration within the flow that deposited the layer. Variations in the turbidite thickness with respect to water depth can be used to estimate the degree of relative vertical movement within the floor of Escanaba Trough. In the area of hydrothermal activity near ODP Site 1038, uplift of as much as 140 m has occurred over the past 8 kyr.

California↗
USGS2012

Modeling crater topography and albedo from monoscopic Viking orbiter images: 1. Methodology

A new photoclinometric technique for extraction of topographic data from single planetary images is presented that overcomes many previous limitations of photoclinometry. The procedure fully compensates for oblique viewing geometry prevalent in spacecraft images. Albedo variations have been one of the most serious obstacles in the application of photoclinometry to planetary surfaces. This problem is overcome in the topographic solution by simultaneously utilizing brightness data from a pair of profiles; both segments are assumed to have the same topographic and albedo variations along their lengths. Profile directions are chosen where the orientation of downslope or upslope is obvious, thus resolving a major ambiguity in photoclinometry. This requirement is particularly easy to satisfy for craters and not very difficult for many irregular features. An additional procedure is presented that eliminates even the requirement of topographic symmetry along the pair of profiles. If two profiles have the same relief but their shapes are very different, another method can be used in an iterative process to derive topographic profiles; however, this procedure does, require that the albedo not vary along the profiles. Test results indicate that both procedures have an accuracy and precision of approximately 2° for slopes of typical bowl-shaped craters, which translates to approximately 5% for depths.

Journal of Geophysical Research Solid Earth↗
USGS2007

Generation of fine hydromagmatic ash by growth and disintegration of glassy rinds

The deposits of mafic hydromagmatic eruptions are more fine grained and variable in vesicularity than dry magmatic deposits. Blocky, equant shapes of many hydromagmatic clasts also contrast with droplet, thread, and bubble wall morphology of dry magmatic fragments. Small (<∼180 μ m), blocky hydromagmatic pyroclasts have traditionally been interpreted to result from discrete vapor explosions, although such explosions tend to occur only under certain conditions. This paper considers a process of hydromagmatic ash formation that involves repeated growth and disintegration of glassy rinds on pyroclast surfaces as they deform within turbulent flows. This process, termed “turbulent shedding”, may occur during the expansion phase of vapor explosions or during turbulent but nonexplosive mixing of magma with water, steam, or water sprays. The occurrence of turbulent shedding and the resulting fragment sizes depend on the timescale for rind growth and the timescale between disturbances that remove or disintegrate glassy rinds. Turbulent shedding is directly observable in some small littoral jets at Kilauea. Calculations suggest that, in the presence of liquid water or water sprays, glassy rinds having a thickness of microns to millimeters should form in milliseconds to seconds. This is similar to the timescale between turbulent velocity fluctuations that can shred lava globules and remove such rinds. The fraction of a deposit consisting of fine ash should increase with the duration of this process: Large‐scale Surtseyan jets generate hundreds or thousands of shedding events; bubble bursts or tephra jets at Kilauea's coast may produce only a few.

Journal of Geophysical Research B: Solid Earth↗
USGS2025

Advancing current understanding of Martian impact-generated hydrothermal systems through novel coupled modeling: Insights from Gale, Jezero, and other craters

Impact-generated hydrothermal systems are considered potentially habitable environments on Mars, Earth, and other planetary bodies for microbial life. However, there is an ongoing debate regarding what geological features on Mars provide definitive evidence for such systems. Although earlier studies have modeled hydrothermal processes in Martian craters, they often lacked integration with shock physics hydrocodes to constrain initial impact conditions. The importance of this two-code coupling was demonstrated by successfully replicating alteration signatures in the Earth's Haughton impact structure. In this study, we use a similar two-code approach, combining the iSALE hydrocode with the HYDROTHERM hydrothermal model to simulate the full evolution of impact-generated hydrothermal systems. We apply this method to craters the size of Jezero (∼50 km) and Gale (∼154 km) in diameter. Although Jezero's interior is largely buried, our results align with hypothesized hydrothermal vents and alteration minerals near central uplifts in similarly sized exposed craters, such as Toro and Auki. Furthermore, our models correspond to alteration patterns observed by the Curiosity in the lower layers of Mount Sharp, which may represent remnants of impact-driven hydrothermal activity. A key finding is that these systems may persist much longer than previously estimated. Our simulations suggest that a Jezero-sized system could remain habitable for thermophiles for approximately 720,000 years, whereas a Gale-sized system could persist for nearly 2 million years. Additionally, simulations under unsaturated crustal conditions reveal that air-dominated near-surface layers can suppress vertical fluid flow, enabling deep subsurface alteration without producing detectable mineral signatures at the surface.

Journal of Geophysical Research - Planets↗
USGS2012

Analysis of ERS 1 synthetic aperture radar data of frozen lakes in northern Montana and implications for climate studies

Lakes that freeze each winter are good indicators of regional climate change if key parameters, such as freeze-up and breakup date and maximum ice thickness, are measured over a decade-scale time frame. Synthetic aperture radar (SAR) satellite data have proven to be especially useful for measurement of climatologically significant parameters characteristic of frozen lakes. In this paper, five lakes in Glacier National Park, Montana, have been studied both in the field and using Earth Remote-Sensing Satellite (ERS) 1 SAR data during the 1992-1993 winter. The lakes are characterized by clear ice, sometimes with tubular or rounded bubbles, and often with a layer of snow ice on top of the clear ice. They are also often snow covered. Freeze-up is detected quite easily using ERS 1 SAR data as soon as a thin layer of ice forms. The effect of snow ice on the backscatter is thought to be significant but is, as yet, undetermined. On the five lakes studied, relative backscatter was found to increase with ice thickness until a maximum was reached in February. Breakup, an often ill-defined occurrence, is difficult to detect because surface water causes the SAR signal to be absorbed, thus masking the ice below. Comparison of the bubble structure of thaw lakes in northern Alaska with lakes in northern Montana has shown that the ice structure is quite different, and this difference may contribute to differential SAR signature evolution in the lakes of the two areas.

Journal of Geophysical Research↗
USGS2023

Strength recovery in quartzite is controlled by changes in friction in experiments at hydrothermal conditions up to 200°C

The rate of fault zone restrengthening between earthquakes can be influenced by both frictional and cohesive healing processes. Friction is dependent on effective normal stress while cohesion is independent of normal stress, potentially explaining—in part—the lack of depth dependence of earthquake stress drops. Although amenable to laboratory testing, few studies have systematically addressed the normal stress dependence of restrengthening rate. This is partially due to difficulty in separating relative contributions of friction and cohesion in recovery of fault strength. We present results from a series of slide-hold-slide tests on thin layers (≤10 𝜇m) of ultrafine quartz gouge that develop during shearing of initially bare-surface quartzite. Tests were conducted at 10 MPa constant pore pressure, 20–200 MPa constant effective normal stress, and temperatures of 22°–200°C. Restrengthening, defined as the difference between peak shear stress measured after resumption of sliding and steady-state sliding shear stress, increases with the log of hold duration. The 200°C healing rate, 0.014 per e-fold increase in time, is comparable to that determined from seismological observations along the Calaveras Fault, California. Construction of Mohr-Coulomb failure envelopes shows that changes in cohesion are small (<1 MPa) and independent of hold durations to 10 5 s, indicating that the increased strength is due to changes in the friction coefficient. These experimental results are inconsistent with the hypothesis that cohesive healing explains the depth independence of earthquake stress drop, but higher temperatures, longer time-scales, and more complex mineralogy could facilitate cohesive healing in natural fault systems.

Journal of Geophysical Research - Solid Earth↗
USGS2012

Strain accumulation rates in the western United States between 1970 and 1978

The rate of dilatation and the rate and direction of shear have been determined from trilateration data for 23 Geodolite networks in the western United States. Sixteen nets are located along the San Andreas fault system between Point Reyes, California, and the United States‐Mexico border. Other locations are across the Garlock fault in California; across Puget Sound near Seattle, Washington; near Hanford in eastern Washington; near Hebgen Lake in Montana; across the Wasatch fault at Ogden, Utah; across the Rio Grande rift at Socorro, New Mexico; and Dixie Valley in Nevada; and at the northern end of Owens Valley on the California‐Nevada border. Implicit in the treatment are the assumptions that the strain was accumulating at a constant rate over the time period (within the interval 1970–1978) and over the local area (usually about 50‐km diameter) covered by the surveys. Of the nets located away from the San Andreas fault, only Ogden and Hebgen show significant strain accumulation. At Ogden the deformation is principally an east‐west compression of 0.23±0.05 μstrain/yr and at Hebgen Lake a northeast‐southwest extension of 0.17±0.03 μstrain/yr. Along the San Andreas fault system the rate of shear is 0.2 to 0.4μ/yr. The direction of shear agrees very well with the surface strike of nearby faults. This agreement is maintained even in regions like the ‘big bend,’ where both the fault strike and the observed shear direction are more westerly than they are elsewhere. Shear strain in northern California appears to be concentrated more closely on the faults, whereas in southern California the strain is a broader, smoother feature. In the San Francisco Bay area the strain data indicate slip at depth on both the San Andreas and the Calaveras faults. In addition to the observed shear the nets in California indicate a negative dilatation (areal decrease) of about 0.2 μstrain/yr. This dilatation is unexplained, but the following sources appear unlikely: (1) systematic survey error, (2) an association with the southern California uplift, (3) an association with the big bend in the San Andreas fault in Southern California, or (4) the result of the superposition of a uniaxial strain on the Pacific‐North American plate boundary shear.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Inducing in situ, nonlinear soil response applying an active source

[1] It is well known that soil sites have a profound effect on ground motion during large earthquakes. The complex structure of soil deposits and the highly nonlinear constitutive behavior of soils largely control nonlinear site response at soil sites. Measurements of nonlinear soil response under natural conditions are critical to advancing our understanding of soil behavior during earthquakes. Many factors limit the use of earthquake observations to estimate nonlinear site response such that quantitative characterization of nonlinear behavior relies almost exclusively on laboratory experiments and modeling of wave propagation. Here we introduce a new method for in situ characterization of the nonlinear behavior of a natural soil formation using measurements obtained immediately adjacent to a large vibrator source. To our knowledge, we are the first group to propose and test such an approach. Employing a large, surface vibrator as a source, we measure the nonlinear behavior of the soil by incrementally increasing the source amplitude over a range of frequencies and monitoring changes in the output spectra. We apply a homodyne algorithm for measuring spectral amplitudes, which provides robust signal-to-noise ratios at the frequencies of interest. Spectral ratios are computed between the receivers and the source as well as receiver pairs located in an array adjacent to the source, providing the means to separate source and near-source nonlinearity from pervasive nonlinearity in the soil column. We find clear evidence of nonlinearity in significant decreases in the frequency of peak spectral ratios, corresponding to material softening with amplitude, observed across the array as the source amplitude is increased. The observed peak shifts are consistent with laboratory measurements of soil nonlinearity. Our results provide constraints for future numerical modeling studies of strong ground motion during earthquakes.

Journal of Geophysical Research B: Solid Earth↗
USGS2023

Spatial and temporal analysis of geologic slip rates, Cucamonga Fault, California, USA: Implications for along-strike applications and multi-fault rupture

To constrain fault processes and hazard, fault slip rates may be extrapolated over different fault lengths or time intervals. Here, we investigate slip rates for the Cucamonga Fault (CF). The CF is located at the junction of the Transverse Range fault system with the San Andreas and San Jacinto Faults, and it is hypothesized to connect with these faults, promoting the propagation of large, multi-fault earthquakes. Previous work has shown that CF displacements on late Quaternary alluvial fan surfaces are highly variable along strike. We present two new 10 Be surface exposure ages from depth profiles on the alluvial fans. Slip rates are consistent with a rate of 1.4 ± 0.3 m/kyr over time intervals of ∼20, ∼30, and ∼40 kyr. If the CF participates in multi-fault ruptures, then these earthquakes occur either rarely or with sufficient regularity to maintain apparently steady rates over multiple intervals. We also explore along-strike fault displacement variability using a calibrated morphological model. The model successfully reproduces scarp profiles and indicates that fault displacement variability can be explained in part by scarp age but not uplift rate. We infer that both erosion by ephemeral gullying and distributed deformation contribute to fault displacement variability, although both are difficult to detect confidently without excavations across the scarp. These investigations show that better characterization of cumulative-slip variability along strike may improve accuracy and precision of slip rates. Slip rates that do not consider epistemic uncertainties may not be suitable for extrapolation over longer fault sections.

California↗
USGS1999

P wave crustal velocity structure in the greater Mount Rainier area from local earthquake tomography

We present results from a local earthquake tomographic imaging experiment in the greater Mount Rainier area. We inverted P wave arrival times from local earthquakes recorded at permanent and temporary Pacific Northwest Seismograph Network seismographs between 1980 and 1996. We used a method similar to that described by Lees and Crosson [1989], modified to incorporate the parameter separation method for decoupling the hypocenter and velocity problems. In the upper 7 km of the resulting model there is good correlation between velocity anomalies and surface geology. Many focal mechanisms within the St. Helens seismic zone have nodal planes parallel to the epicentral trend as well as to a north-south trending low-velocity trough, leading us to speculate that the trough represents a zone of structural weakness in which a moderate ( M 6.5–7.0) earthquake could occur. In contrast, the western Rainier seismic zone does not correlate in any simple way with anomaly patterns or focal mechanism fault planes, leading us to infer that it is less likely to experience a moderate earthquake. A ∼10 km-wide low-velocity anomaly occurs 5 to 18 km beneath the summit of Mount Rainier, which we interpret to be a signal of a region composed of hot, fractured rock with possible small amounts of melt or fluid. No systematic velocity pattern is observed in association with the southern Washington Cascades conductor. A midcrustal anomaly parallels the Olympic-Wallowa lineament as well as several other geophysical trends, indicating that it may play an important role in regional tectonics.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Composition and variation of noise recorded at the Yellowknife Seismic Array, 1991-2007

We analyze seismic noise recorded on the 18 short-period, vertical component seismometers of the Yellowknife Seismic Array (YKA). YKA has an aperture of 23 km and is sited on cratonic lithosphere in an area with low cultural noise. These properties make it ideal for studying natural seismic noise at periods of 1-3 s. We calculated frequency-wave number spectra in this band for over 6,000 time windows that were extracted once per day for 17 years (1991-2007). Slowness analysis reveals a rich variety of seismic phases originating from distinct source regions: R g waves from the Great Slave Lake; L g waves from the Atlantic, Pacific, and Arctic Oceans; and teleseismic P waves from the north Pacific and equatorial mid-Atlantic regions. The surface wave energy is generated along coastlines, while the body wave energy is generated at least in part in deep-water, pelagic regions. Surface waves tend to dominate at the longer periods and, just as in earthquake seismograms, L g is the most prominent arrival. Although the periods we study are slightly shorter than the classic double-frequency microseismic band of 4-10 s, the noise at YKA has clear seasonal behavior that is consistent with the ocean wave climate in the Northern Hemisphere. The temporal variation of most of the noise sources can be well fit using just two Fourier components: yearly and biyearly terms that combine to give a fast rise in microseismic power from mid-June through mid-October, followed by a gradual decline. The exception is the R g energy from the Great Slave Lake, which shows a sharp drop in noise power over a 2-week period in November as the lake freezes. The L g noise from the east has a small but statistically significant positive slope, perhaps implying increased ocean wave activity in the North Atlantic over the last 17 years. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

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↗
USGS2017

Dynamic rupture modeling of the M7.2 2010 El Mayor-Cucapah earthquake: Comparison with a geodetic model

The 2010 M w 7.2 El Mayor-Cucapah earthquake is the largest event recorded in the broader Southern California-Baja California region in the last 18 years. Here we try to analyze primary features of this type of event by using dynamic rupture simulations based on a multifault interface and later compare our results with space geodetic models. Our results show that starting from homogeneous prestress conditions, slip heterogeneity can be achieved as a result of variable dip angle along strike and the modulation imposed by step over segments. We also considered effects from a topographic free surface and find that although this does not produce significant first-order effects for this earthquake, even a low topographic dome such as the Cucapah range can affect the rupture front pattern and fault slip rate. Finally, we inverted available interferometric synthetic aperture radar data, using the same geometry as the dynamic rupture model, and retrieved the space geodetic slip distribution that serves to constrain the dynamic rupture models. The one to one comparison of the final fault slip pattern generated with dynamic rupture models and the space geodetic inversion show good agreement. Our results lead us to the following conclusion: in a possible multifault rupture scenario, and if we have first-order geometry constraints, dynamic rupture models can be very efficient in predicting large-scale slip heterogeneities that are important for the correct assessment of seismic hazard and the magnitude of future events. Our work contributes to understanding the complex nature of multifault systems.

California↗