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

Geodetic evidence for en echelon dike emplacement and concurrent slow slip during the June 2007 intrusion and eruption at Kilauea volcano, Hawaii

A series of complex events at Kīlauea Volcano, Hawaii, 17 June to 19 June 2007, began with an intrusion in the upper east rift zone (ERZ) and culminated with a small eruption (1500 m 3 ). Surface deformation due to the intrusion was recorded in unprecedented detail by Global Positioning System (GPS) and tilt networks as well as interferometric synthetic aperture radar (InSAR) data acquired by the ENVISAT and ALOS satellites. A joint nonlinear inversion of GPS, tilt, and InSAR data yields a deflationary source beneath the summit caldera and an ENE-striking uniform-opening dislocation with ~2 m opening, a dip of ∼80° to the south, and extending from the surface to ~2 km depth. This simple model reasonably fits the overall pattern of deformation but significantly misfits data near the western end of an inferred dike-like source. Three more complex dike models are tested that allow for distributed opening including (1) a dike that follows the surface trace of the active rift zone, (2) a dike that follows the symmetry axis of InSAR deformation, and (3) two en echelon dike segments beneath mapped surface cracks and newly formed steaming areas. The en echelon dike model best fits near-field GPS and tilt data. Maximum opening of 2.4 m occurred on the eastern segment beneath the eruptive vent. Although this model represents the best fit to the ERZ data, it still fails to explain data from a coastal tiltmeter and GPS sites on Kīlauea's southwestern flank. The southwest flank GPS sites and the coastal tiltmeter exhibit deformation consistent with observations of previous slow slip events beneath Kīlauea's south flank, but inconsistent with observations of previous intrusions. Slow slip events at Kīlauea and elsewhere are thought to occur in a transition zone between locked and stably sliding zones of a fault. An inversion including slip on a basal decollement improves fit to these data and suggests a maximum of ~15 cm of seaward fault motion, comparable to previous slow-slip events.

Hawai'i↗
USGS2020

Methods for rapidly estimating velocity precision from GNSS time series in the presence of temporal correlation: A new method and comparison of existing methods

Time series of position estimates from Global Navigational Satellite System (GNSS) are used to measure the velocities of points on the surface of the Earth. Along with the velocity estimates, a measure of the precision is needed to assess the quality of the velocity measurement. Here, I evaluate rate uncertainties provided by four different methods that have been applied to geodetic time series. The most rigorous approach uses a data covariance that incorporates a variety of noise processes relevant to geodetic time series but is computationally demanding. Two other approaches are efficient algorithms and are used widely, but both can provide less rigorous estimates of the rate uncertainty. I propose and evaluate a fourth method, which provides estimates of rate uncertainty closer to the rigorous approach but is significantly less computationally demanding. I have evaluated all three methods against the more rigorous method using both simulations and time series from 190 GNSS sites. For data best characterized as having a flicker type noise process, one of the widely used methods overestimates the uncertainty by up to a factor of 2, while the other widely used method underestimates the uncertainty by less than a factor of 2. For a random-walk process, both methods underestimate the rate uncertainty by a factor of 3 to 5.

Journal of Geophysical Research↗
USGS2012

Cooling rate and thermal structure determined from progressive magnetization of the dacite dome at Mount St. Helens, Washington

Our study of a magnetic anomaly associated with the recently active dacite dome at Mount St. Helens suggests that the dome consists of a hot, nonmagnetized core surrounded by a cool, magnetized carapace and flanking talus. The talus does not contribute to the anomaly because its constituent blocks are randomly oriented. Temporal changes in the magnetic anomaly indicate that the magnetized carapace thickened at an average rate of 0.03±0.01 m/d from 1984 to 1986. Petrographic and rock magnetic properties of dome samples indicate that the dominant process responsible for these changes is magnetization of extensively oxidized rock at progressively deeper levels within the dome as the rock cools through its blocking temperature, rather than subsequent changes in magnetization caused by further oxidation. Newly extruded material cools rapidly for a short period as heat is conducted outward in response to convective heat loss from its surface. The cooling rate gradually declines for several weeks, and thereafter the material cools at a relatively constant rate by convective heat loss from its interior along fractures that propagate inward. The rate of internal convective heat loss through fractures varies with rainfall, snowmelt, and large-scale fracturing during subsequent eruptive episodes. In accordance with a model for solidification of the 1959 lava lake at Kilauea Iki, Hawaii, we picture the dome's magnetized carapace as being a two-phase, porous, convective zone separated from the nonmagnetized core of the dome by a thin, single-phase conductive zone. As a consequence of the heat balance between the conductive and convective zones, the blocking-temperature isotherm migrates inward at a relatively constant rate. If the dome remains inactive, the time scale for its complete magnetization is estimated to be 18–36 years, a forecast which can be refined by shallow drilling into the dome and by continuing studies of its growing magnetic anomaly.

Journal of Geophysical Research Solid Earth↗
USGS2012

Comparison of ground motion from tremors and explosions in deep gold mines

Seismic body waves, from tamped chemical explosions, two with yields of 50 and one of 150 kg, were compared with corresponding data from three mining-induced tremors with a view to testing methods of discriminating between the two types of events. Detonated at depths of about 2 km, all three explosions generated P waves for which the low-frequency spectral asymptotes agree well with corresponding results scaled down from nuclear shots at the Nevada Test Site. For the two smaller explosions, recorded underground at hypocentral distances ranging from 234 to 871 m, the P wave corner frequencies, defined by the intersections of the low- and high-frequency asymptotes, are also in fine agreement with scaled-down results from the Nevada Test Site; for the 150-kg explosion, which was recorded only at the surface, it appears that 2.39 km of upward propagation caused the corner frequency to be at least a factor of 4 lower than anticipated due to attenuation. All three explosions generated S waves that appear to be a consequence of deviatoric stress release in the immediate environs of the explosions. The three tremors analyzed here generated P and S waves whose spectra and source parameters agree well with standard earthquake source models. With regard to discrimination based on P wave spectra, we conclude that for events of fixed low-frequency spectral asymptotes, the explosions typically have higher corner frequencies than tremors or earthquakes, although counterexamples certainly exist. Interestingly, the 150-kg explosion was identified as such on the basis of P and S wave polarities that are incompatible with the normally expected double-couple source model; instead, these initial motions are consistent with an explosion in conjunction with normal faulting. The body wave spectra of this explosion and those of a nearby tremor, however, were indistinguishable.

Journal of Geophysical Research Solid Earth↗
USGS2026

Role of volcano-tectonic interactions during early-phase magma-assisted continental rifting

We study volcanic-tectonic interactions during the early phase of rifting in the Natron Rift, a magmatic segment of the East African Rift System. We use the open-source crustal deformation modeling code, PyLith, to model deflation during the 2007–2008 explosive eruptions and inflation during 2016–2023. We compare predictions of surface displacements resulting from magma reservoir inflation from 3D models with and without topography and homogeneous and heterogeneous material properties against Global Navigation Satellite System displacements from 2016 to 2023. We prefer the model with heterogeneous material properties and topography for fitting the observed displacements. Including topography improves the fit to GNSS data and modulates stress changes from shallow magma sources. Our numerical modeling of Coulomb stress changes suggests that deflation of a magma reservoir under the active Ol Doinyo Lengai volcano could trigger oblique normal slip on the Natron Fault with Coulomb stress increases exceeding ∼0.1 MPa on the volcano-adjacent fault segment. Fault slip on the Natron Fault is less likely during inflation of the magma reservoir because normal compressive stress increases on the fault. We also demonstrate that Coulomb stress changes are sensitive to the inclusion of topography in our model, particularly because the depth of the magma reservoir is similar to the amplitude of the topographic relief.

Ol Doinyo Lenga↗
USGS2003

Oceans on Mars: An assessment of the observational evidence and possible fate

If the large Late Hesperian outflow channels were eroded by extensive floods, as appears likely, then large bodies of water must have once occupied the northern plains during that period. Previous estimates of the sizes of bodies of water in the northern lowlands range up to 3 ?? 108 km3. Several contacts have been previously mapped around the edges of the northern plains and interpreted to be shorelines remaining from these former standing bodies of water. We examine the elevations and geologic relations along these contacts in detail and find little support for their interpretation as shorelines. Some contacts are clearly of volcanic origin, and all have significant variations in elevation. Better support for the former presence of water over large parts of the northern plains is provided by the Vastitas Borealis Formation (VBF). Most of the post-Noachian fill within the northern basin is ridged plains of Lower Hesperian age, interpreted to be volcanic in origin. Overlying the ridged plains is the VBF, a thin veneer of material of Upper Hesperian age. The VBF may have been deposited from large floods. Support for this interpretation is the similarity in age between the outflow channels and the VBF, the presence of the VBF at the lower ends of the outflow channels, and identification of numerous features in the outcrop areas of the VBF that are suggestive of basal melting of an ice sheet. To cover all the area over which the VBF is exposed would require ???2.3 ?? 107 km3 of water. Spread over the entire surface of Mars, this volume is equal to a global layer (global equivalent layer, or GEL) ???156 m deep. We find no support for the larger estimates of ocean volumes that range up to 3 ?? 108 km3 and which imply comparable amounts of water per unit area as are currently present on the surface of the Earth. Under present climatic conditions on Mars an ocean would freeze in a geologically short time period (???104 years), then would sublimate away at rates strongly dependent on the presence or absence of debris on the ice surface. The present VBF is interpreted as a sublimation residue from the ponded outflow channel effluents. The fate of a volume of water thought to have been emplaced by the outflow channels (???2.3 ?? 107 km3) is largely accounted for by the presence of other existing reservoirs on the planet. An approximately 20-30 m GEL of water is estimated to be in the present polar caps, and a 50 m GEL may have escaped to space since the Hesperian, leaving ???80 m GEL unaccounted for. This amount may be partly trapped in other volatile-rich deposits on the surface, and a significant amount could have reentered the groundwater system by south polar basal melting and been progressively cold-trapped at the base of a growing cryosphere. On the basis of our assessment of the Hesperian-aged deposits, we predict that testing of the Clifford and Parker [2001] hypothesis that a Noachian-aged ocean covered up to one third of the surface of Mars will be made very difficult by the enhanced degradation rates in the Noachian and subsequent geological events in the northern lowlands.

Journal of Geophysical Research E: Planets↗
USGS2024

Limited preservation of strike-slip surface displacement in the geomorphic record

Offset geomorphic markers are commonly used to interpret slip history of strike-slip faults and have played an important role in forming earthquake recurrence models. These data sets are typically analyzed using cumulative probability methods to interpret average amounts of slip in past earthquakes. However, interpretation of the geomorphic record to infer surface slip history is complicated by slip variability, measurement uncertainty, and modification of offset features in the landscape. To investigate how well geomorphic data record surface slip, we use offset measurements from recent strike-slip surface ruptures ( n = 39), faults with geomorphic evidence of multiple strike-slip earthquakes ( n = 29), and synthetic slip distributions with added noise ( n > 10,000) to examine the constraints of the geomorphic record and the underlying assumptions of the cumulative offset probability distribution analysis method. We find that the geomorphic record is unlikely to resolve more than two paleo-slip distributions, except in specific cases with low slip variability, high slip-per-event, and semiarid climate. In cases where site-specific conditions allow for interpretation of more than two earthquakes, lateral extrapolation along a fault is not straightforward because on-fault displacement and distributed deformation may be spatially variable in each earthquake. We also find that average slip in modern earthquakes is adequately recovered by probability methods, but the reported prevalence of strike-slip faults with characteristic slip history is not supported by geomorphic data. We also propose updated methods to interpret slip history and construct uncertainty bounds for paleo-slip distributions.

Journal of Geophysical Research: Solid Earth↗
USGS2006

Erosion rates at the Mars Exploration Rover landing sites and long-term climate change on Mars

Erosion rates derived from the Gusev cratered plains and the erosion of weak sulfates by saltating sand at Meridiani Planum are so slow that they argue that the present dry and desiccating environment has persisted since the Early Hesperian. In contrast, sedimentary rocks at Meridiani formed in the presence of groundwater and occasional surface water, and many Columbia Hills rocks at Gusev underwent aqueous alteration during the Late Noachian, approximately coeval with a wide variety of geomorphic indicators that indicate a wetter and likely warmer environment. Two-toned rocks, elevated ventifacts, and perched and undercut rocks indicate localized deflation of the Gusev plains and deposition of an equivalent amount of sediment into craters to form hollows, suggesting average erosion rates of ???0.03 nm/yr. Erosion of Hesperian craters, modification of Late Amazonian craters, and the concentration of hematite concretions in the soils of Meridiani yield slightly higher average erosion rates of 1-10 nm/yr in the Amazonian. These erosion rates are 2-5 orders of magnitude lower than the slowest continental denudation rates on Earth, indicating that liquid water was not an active erosional agent. Erosion rates for Meridiani just before deposition of the sulfate-rich sediments and other eroded Noachian areas are comparable with slow denudation rates on Earth that are dominated by liquid water. Available data suggest the climate change at the landing sites from wet and likely warm to dry and desiccating occurred sometime between the Late Noachian and the beginning of the Late Hesperian (3.7-3.5 Ga). Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2022

Three Mw ≥ 4.7 earthquakes within the Changning (China) shale gas field ruptured shallow faults intersecting with hydraulic fracturing wells

From 2017 to 2019, three destructive earthquakes (27 January 2017 Mw 4.7, 16 December 2018 Mw 5.2, and 3 January 2019 Mw 4.8) occurred in the Changning shale gas field in the southwest Sichuan Basin, China. Previous seismological studies attributed these events to hydraulic fracturing (HF), but were unable to identify the causative seismogenic faults and their slip behaviors. Here, we use Sentinel-1 synthetic aperture radar data to measure surface deformation triggered by the three events and conduct geodetic inversions to characterize their rupture models. The resulting coseismic interferograms show prominent surface deformation with the maximum line-of-sight displacements of up to 4 cm. The inversion results show that all three earthquakes mainly ruptured sedimentary formations above the shale gas bed, in the upper 3 km of the crust, with slip magnitudes ranging from 8.5 to 15 cm, and stress drops ranging from ∼1.8 to ∼3.3 MPa. Their source faults intersect with horizontal HF wells, but do not root in the crystalline basement. Combined with the reported difficulty of increasing HF operation pressures prior to the three events, we argue that they were most likely induced by direct injection of pressurized fluids into the fault zones. Crustal deformation patterns inferred from regional topography and GPS velocities highlight that the Changning field is located within a triple junction region near the southeastern margin of the Tibetan Plateau with large deformation gradients; such conditions are not only favorable to the development of critically stressed faults, but also facilitate the occurrence of at least moderate magnitude earthquakes.

Changning shale gas field, Sichuan Basin↗