USGS Science⌕ Search

SEARCH · USGS Science

Results for “Article”

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.

979 records · Page 26Linked to original sources
USGS2024

Tungsten skarn quantitative mineral resource assessment and gold, rare earth elements, graphite, and uranium qualitative assessments of the Kuldjuktau and Auminzatau Ranges, in the central Kyzylkum region, Uzbekistan

A new quantitative mineral resource assessment for tungsten skarn was conducted for the Auminzatau and Kuldjuktau mountain ranges in Central Uzbekistan, along with qualitative assessments of orogenic gold, rare earth elements (REEs), amorphous graphite, and uranium. By integrating a variety of geological, geochemical, geophysical, and remote sensing data sets, estimates of undiscovered tungsten skarn deposits in permissive tracts are combined with grade and tonnage distributions of known deposits to generate probabilistic estimates of undiscovered resources. Undiscovered deposits in Auminzatau are estimated to contain median resources of 98 thousand metric tons (kt) of WO 3 with a 70 percent (%) probability of at least 28 kt and a 10% probability of at least 468 kt, of which 16 kt to 293 kt may be economic to extract. In Kuldjuktau, the undiscovered deposits are estimated to contain median resources of 27 kt of WO 3 with a 60% probability of at least 12 kt and a 10% probability of at least 208 kt, of which 5 kt to 132 kt may be economic to extract. Our results suggest that the Auminzatau–Kuldjuktau Mountains area is highly prospective for additional discovery of significant Au and U resources and has low prospectivity for discovery of significant REE and graphite resources.

Minerals↗
USGS2025

Latitudinal gradients of snow contamination in the Rocky Mountains associated with anthropogenic sources

Seasonal snow is an important source of drinking water and recreation, and for agriculture in the Rocky Mountain region. Monitoring snow-water quality can inform on the effects to the albedo and energy balance of the snowpack, and the sources of natural and anthropogenic aerosol and gases. This study analyzed metals in the seasonal snowpack from water year (WY) 2018 for 49 sites. Calcium, lanthanum, and cerium concentrations support the importance of mineral dust to the southern Rocky Mountains. Mercury (Hg), zinc (Zn), and cadmium (Cd) concentrations showed a similar spatial pattern to mineral dust, whereas antimony (Sb) concentrations were highest in the northern Rocky Mountains. To assess the relative contributions from dust versus anthropogenic contaminant sources, enrichment factors (EF) were calculated, with values above 10 indicating anthropogenic contamination. For Cd, Hg, Sb, and Zn, EF values exceeded 10 at northern sites. These observations were compared to spatial trends of EF values of Hg from WY2009 to WY2018, regional monitoring networks, and back trajectory analyses. The agreement between these datasets revealed temporally consistent contaminant sources and/or transport processes to the northern Rocky Mountains snowpack. Sources include current and historical mining and smelting in the region. Strategies to limit the emissions of these metals to the Northern Rockies could benefit from focusing on remediation of contaminated sites, and continued monitoring and mitigation of active mining and smelting.

Colorado , Idaho, Montana, New Mexico, Utah, Wyomi↗
USGS2025

A new groundwater energy transport model for the MODFLOW hydrologic simulator

Heat transport in the subsurface is an important aspect of research related to the effects of a warming climate on ecological services (i.e., cold-water refugia); the development of geothermal resources for energy banking schemes (i.e., aquifer thermal energy storage [ATES]); and the effects of temperature on other aspects of groundwater quality, such as nutrient cycling. Historically, simulation of heat transport using the MODFLOW groundwater simulator and related codes was performed by scaling the input parameters of a solute-transport model to emulate heat transport. However, that approach required additional pre- and post-processing of input and output and could not account for the variation in effective thermal storage and transport properties during transient, unsaturated flow, for example. True heat-transport capabilities in the context of MODFLOW were first introduced in a variant called USG-Transport. More recently, a new groundwater energy-transport (GWE) model type has been added to MODFLOW 6, the core version of the MODFLOW hydrologic simulator. GWE supports the simulation of heat transport on structured or unstructured grids as well as within and between features of advanced packages that represent streams, lakes, multi-aquifer wells, and the unsaturated zone. GWE is integrated within MODFLOW 6 and is accessible through the FloPy Python package and the MODFLOW 6 application programming interface (API). An example simulation demonstrates conduction between grid cells through both the water and the solid aquifer material, including thermal bleeding from saturated overburden cells into a groundwater flow field.

Groundwater↗
USGS2025

A review of post-wildfire adaptations of surface-water-quality models: Synthesis, gaps, and opportunities

As wildfires increasingly affect water-supply watersheds, the demand for models to predict water-quality responses is increasing. This work reviews and synthesizes existing post-wildfire applications of water-quality models in the context of geographic and ecohydrological distribution, hydrologic and water-quality response process representation, model parameterization, model and input data scales, model calibration data availability, as well as calibration and performance evaluation approaches. Emphasis is placed on models that simulate water-quality output, rather than sediment and erosional response as the primary focus. Here, identified gaps and opportunities to advance the post-wildfire application of water-quality models include: 1. applying models in under-represented geographic and ecohydrologic regions, 2. simulating multiple streamflow generation mechanisms, including groundwater, with an emphasis on shifting dominant flow pathways as the landscape recovers following wildfire, 3. adding studies that include the simulation of metals, 4. incorporating more biogeochemical and in-stream processes to model applications, 5. applying finer spatial and temporal resolution of precipitation data input as well as finer spatial resolution hydrologic response units, 6. implementing fully distributed grid or element models or finer resolution response units to capture burn severity heterogeneity, 7. collecting enhanced water-quality data for model calibration and validation, 8. conducting model-intercomparison studies, and 9. developing model parameter value guidance in post-wildfire applications. These identified gaps and opportunities may assist users in deciding on key processes and approaches to consider in modeling post-wildfire water-quality conditions.

Science of the Total Environment↗
USGS2026

Dynamic drainage reorganization in Eastern Tibet: Insights from the Yangtze River first bend

The modern drainage network of eastern Tibet is widely believed to have developed through a series of river capture and flow reversal events; however, the timing and mechanisms driving this reorganization remain contentious. Among these events, the river capture that formed the First Bend of the Yangtze River (YFB) stands out as both iconic and particularly debated. Here we present sedimentary provenance data from the Late Miocene–Quaternary Dali Basin, located south of the YFB, which indicate that a southward-flowing Jinsha River (i.e., the present-day upper Yangtze River) sourced sediment to the Dali basin at ∼7.4–6.4 Ma in a drainage configuration different from that of today. Because this interval postdates the initial establishment of a near-modern Jinsha River system prior to the Miocene, our results imply at least two discrete fluvial reorganizations occurred at the YFB—one preceding ∼7.4 Ma and another following ∼6.4 Ma. By integrating these findings with landscape evolution modeling, we infer that the initiation of rapid uplift of the Yulong-Haba Mountains and the Diancang Shan may have been responsible for these drainage reorganizations. These results underscore that Cenozoic drainage systems on the eastern Tibetan Plateau have evolved dynamically on a short timescale of ∼10 5 –10 6 -year, rather than remaining in a long-term stationary configuration on ∼10 7 -year timescales.

eastern Tibetan Plateau, first bend of the Yangtze↗
USGS2026

Alternative approaches to dryland reclamation enhance vegetation cover and soil stability at a former uranium mine

Introduction Conventional reclamation methods to recover dryland ecosystems after mineral extraction often have low success. Alternative reclamation techniques may help overcome the many challenges to establishing persistent native vegetation and stable soils. Objectives Compare conventional reclamation approaches that focus primarily on soil stabilization to alternative approaches that may improve native ecosystem restoration at the Kanab North Mine in northern Arizona. Methods In a fully crossed experiment, conventional seed mix (native and non-native species) and seeding methods (drill-seeding) were compared to alternative practices, including a native-only seed mix, connectivity modifiers (ConMods) with hand seeding, and a biocrust inoculation. Restoration metrics like plant establishment, species composition, soil surface cover and stability, biocrust establishment, and near-surface wind-driven sediment were monitored over the first 5 years post-reclamation. Results While all restoration metrics improved with time, compared to traditional methods, ConMods increased total, graminoid, and native plant species cover, and the native-only seed mix increased woody, native, and perennial species cover. Biocrust inoculation increased late-successional biocrust development, but biocrusts only persisted over 5 years in combination with ConMods. Soil stability increased in all treatments, but both ConMods and biocrust inoculation increased stability compared to traditional practices. Conclusion Overall, alternative approaches to reclamation tested here (ConMods, native seed, and biocrust inoculation) were more effective than traditional practices at increasing total plant cover, facilitating biocrust establishment, and promoting soil stability after uranium mine reclamation.

Arizona↗
USGS2026

Modeling legacy nitrogen transport under instantaneous, steady-state, and transient groundwater flow conditions

In hydrologic settings where groundwater discharge contributes substantially to surface waters, legacy nitrogen in groundwater can confound surface water nitrogen loads estimated exclusively from current terrestrial sources. Additionally, legacy nitrogen in groundwater can contribute to lagged responses to nitrogen management efforts. Some methods of estimating groundwater contributions to surface water nitrogen loads account for legacy nitrogen, while others do not. The resulting differences are rarely quantified. We used a numerical modeling framework to compare three methods of estimating time-varying annual groundwater nitrogen loads to surface water receptors on eastern Long Island, New York. The instantaneous load method used steady-state contributing areas and includes no temporal groundwater lag. The second method used numerical simulations of nitrogen loads under steady-state flow, which captures groundwater transport lags but omits the annual variability in transient hydrologic stresses. The third method numerically simulated both transient groundwater flow and nitrogen transport to explicitly capture the effects of legacy nitrogen in groundwater. Depending on antecedent nitrogen and hydrologic conditions, historical nitrogen loads estimated from the numerical simulations were sometimes similar (<10% difference) and other times substantially different (±100%) from the instantaneous load estimates. Additionally, simulated future surface water nitrogen loads responded asymptotically over several decades following reductions in terrestrial nitrogen sources, further highlighting the effect of groundwater transport lag times. The comparison of the three methods, quantification of historical interannual variability, and prediction of lagged responses to nitrogen source reductions provide important context for decision makers using estimated groundwater nitrogen loads to help evaluate nitrogen management efficacy.

New York↗
USGS2025

Review: The hydrogeology of critical mineral resources relevant to the energy transition

Attaining the goals of the international treaty on climate change (the Paris Agreement) will greatly increase the demand for the critical minerals required to implement clean-energy technologies. This poses both challenges and opportunities to the hydrogeologic community from several perspectives. Here, important insights that the hydrogeological sciences have to offer for mineral exploration, mineral production, and addressing environmental issues related to mining and mine decommissioning are summarized. This study focuses on copper, cobalt, lithium, and rare earths to represent the broad spectrum of critical minerals and illustrate their relevance by referring to projected demands and production rates. The current understanding of the hydrogeologic processes that form major deposits of these minerals are then summarized. Ore is defined as the naturally occurring material from which minerals of economic value can be extracted, where most ore deposits are the products of complex hydrogeologic couplings between fluid flow, heat transport, solute transport, chemical reactions, and mechanical deformation. Exploration models for the discovery of deeper, hidden deposits are potentially informed by hydrogeologic theory and hydrogeochemical processes. Hydrogeologic understanding and methods are also essential to production and recovery. Longstanding challenges are mine dewatering and (conversely) mine water supply, as well as mineral-extraction practices such as spoil heap leaching and in situ mining. New challenges arise from element extraction from subsurface brines. Finally, the quantity of water use and potential environmental impacts of mining on water quality are at the core of ‘social license’: the approval and acceptance of society to mining activities.

Hydrogeology Journal↗
USGS2025

Statistical approaches for modeling correlated grade and tonnage distributions and applications for mineral resource assessments

Correlations between grade and tonnage exist in mineral resource data compiled from published reports, but they are not always addressed during quantitative assessment of undiscovered mineral resources. Failure to account for correlated grade and tonnage distributions can result in geologically unrealistic assessment results. Current software tools simulate univariate ore tonnage and multivariate resource grades of undiscovered deposits independently. As a result, analysts are forced to rely on ad-hoc solutions to minimize the correlation issues by: 1) creating subsets of data with restricted criteria; 2) truncating grade and tonnage distributions; and 3) testing model robustness using exploratory data analysis. While these methods represent pragmatic solutions, the statistical solutions presented here provide additional options to address real correlations in grade and tonnage data used for mineral resource assessments. We present a modified version of the MapMark4 package in R that introduces two alternatives for modeling grade and tonnage distributions, consisting of a multivariate solution that accounts for correlations between ore tonnage and metal grades and an empirical solution that utilizes simple random sampling with replacement to reproduce coupled grades and tonnages from the input data. We present simulations for contained ore and metal for three case studies representing tungsten skarn, komatiite-hosted nickel, and sediment-hosted carbonate amagmatic zinc-lead (Mississippi Valley-type) deposits. Employing the methods presented here yields quantitative mineral resource assessment results that more closely reflect the empirical distributions of grades and tonnages observed in nature and expands the applicability of these tools for ongoing critical mineral resource assessments.

Applied Computing and Geosciences↗
USGS2025

Irrigated agriculture influences selenium levels in an endangered marsh bird

Selenium bioaccumulation in aquatic food webs poses risks to wildlife, particularly in wetlands receiving irrigation runoff. The Salton Sea, California’s largest lake, is primarily sustained by agricultural drainage. This drainage creates wetland habitat along the lakeshore that many bird species depend on, including the federally endangered Yuma Ridgway’s rail ( Rallus obsoletus yumanensis ). However, these marshes may pose an ecological trap – attracting rails despite high selenium exposure. We captured rails during the 2020–2023 breeding seasons and compared rail selenium levels within three types of marshes (fed with irrigation runoff, Colorado River water, or groundwater). We collected blood, breast feathers, and head feathers of rails in all three water sources for selenium comparisons. We tagged adult rails with GPS transmitters to locate nests and foraging locations where we collected eggshells, unhatched eggs, and prey. We assessed selenium exposure by collecting multiple prey species commonly eaten by rails in all three water sources. Selenium concentrations varied among sampling locations. Selenium concentrations in most sample types were predominately influenced by water source and marsh inflow velocity (sometimes in combination with marsh size). Distance to inflow, however, did not influence selenium concentrations in any sample type. Selenium concentrations were highest in agricultural-fed marshes compared to river-fed and spring-fed marshes. Increased marsh inflow velocities resulted in lower selenium concentrations. Given the risk of an ecological trap, our results suggest that supplementing wetlands with Colorado River water could mitigate selenium bioaccumulation in Yuma Ridgway’s rails.

California↗
USGS2026

Corundum discovered by SuperCam and the Perseverance rover at Jezero crater, Mars

Mars is primarily composed of mafic mineral assemblages and their alteration products, but small, scattered rocks strewn across the landscape offer clues to greater petrological diversity. While traversing the Jezero crater rim, the Perseverance rover encountered several plagioclase-rich light-toned float rocks. SuperCam identified the distinctive signature of corundum (α-Al 2 O 3 ) in these rocks using time-resolved luminescence spectroscopy. Two strong peaks (692.7 and 694.1 nm) with millisecond lifetimes, and additional supporting lines, are consistent with Cr 3+ substitution for Al 3+ in corundum. Corundum forms in Al-rich, Si-depleted environments through magmatic or metamorphic processes. Given the rocks' small size, association with plagioclase, and location on the crater rim, we interpret the most plausible formation scenario to be impact induced metamorphism at the interface of a felsic and a mafic/ultramafic member with the likely action of fluids at some stage, although other possibilities are not excluded.

Geophysical Research Letters↗
USGS2025

A high-resolution 3-D P-wave velocity structure of the south-central Cascadia subduction zone from wide-angle shore-crossing seismic refraction data

This study addresses a significant gap in understanding the features of the south-central Cascadia subduction zone, a region characterized by complex geologic, tectonic, and seismic transitions both offshore and onshore. Unlike other segments along this margin, this area lacks a 3-D velocity model to delineate its structural and geological features on a fine scale. To address this void, we developed a high-resolution 3-D P-wave velocity model using active source seismic data from ship-borne seismic shots recorded on temporary and permanent onshore seismic stations and ocean-bottom seismometers. Our model shows velocity variations across the region with distinct velocity-depth profiles for the Siletz, Franciscan, and Klamath terranes in the overlying plate. We identified seaward dipping high-velocity static backstops associated with the Siletz and Klamath terranes, situated near the shoreline and further inland, respectively. Regions of reduced crustal velocity are associated with crustal faults. Moreover, there is significant along-strike depth variation in the subducting slab, which is about 4 km deeper near the thick, dense Siletz terrane and becomes shallower near the predominantly less-dense Franciscan terrane. This highlights a sudden tectonic and geologic transition at the southern boundary of the Siletz terrane. Our velocity model also indicates slightly increased hydration, though still minimal, in both the oceanic crust and the upper mantle of the subducting plate compared to other parts of the margin.

British Columbia, California, Oregon, Washington↗
USGS2024

Onset of aftershocks: Constraints on the Rate-and-State model

Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R (t)=K(c+t) −p ⁠ , with parameters K , c , and p . The rate‐and‐state (RS) model, which is currently the most popular physics‐based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c ‐value that depends on the size of the coseismic stress change. Because the mainshock‐induced stresses strongly vary in space, the c ‐value should vary accordingly. Short‐time aftershock incompleteness (STAI) in earthquake catalogs has prevented a detailed test of this prediction so far, but the newly developed a ‐positive method for reconstructing the true earthquake rate now allows its testing. Using previously published slip models, we calculate the coseismic stress changes for the six largest mainshocks in Southern California in recent decades and estimate the maximum shear as a scalar proxy of the coseismic stress tensor. Aftershock rates reconstructed for events in different stress ranges show that the rates follow a power law with p = 1 independent of stress with no clear sign of a c ‐value. The onset of the power‐law decay is abrupt and more delayed in areas with smaller stress changes. The observations do not necessarily contradict the RS model, as STAI limits the resolution for early aftershocks, and the RS model can reproduce the observations for specific Aσ values. However, the observations lead to strong constraints, namely Aσ <10 kPa and a power‐law decay of the background rate with distance to the fault, with exponent 2.7.

Seismological Research Letters↗
USGS2024

Testing rate‐and‐state predictions of aftershock decay with distance

We analyze aftershocks of the 2019 M 7.1 Ridgecrest mainshock and isolated M 5–6 mainshocks in southern California to test predictions made by the rate‐and‐state friction model of Dieterich (1994) . Rate‐and‐state friction predicts that the seismicity rate after a stress step follows Omori decay, where the Omori c ‐value, which is the saturation in aftershock rate observed at small times, is larger for smaller stress steps. Put in the context of an aftershock sequence, this predicts that the Omori c ‐value will be systematically larger at greater distances from the mainshock. To our knowledge, this predicted effect has not been observed. In part this may be because the Omori c ‐value is difficult to measure because it often reflects short‐term catalog incompleteness rather than a true saturation in aftershock rate. We explore the dependence of the Omori c ‐value on the distance to the mainshock by applying the “ a ‐positive” method ( van der Elst and Page, 2023 ). This method is insensitive to short‐term aftershock incompleteness and allows resolution of the true aftershock rate deep into the mainshock coda. For aftershocks of the Ridgecrest mainshock and stacked M 5–6 mainshocks, we observe systematic differences in early aftershock rates, relative to mainshock distance, consistent with the predictions of rate‐and‐state friction. Furthermore, for the larger Ridgecrest dataset, we observe that aftershocks nearer to the mainshock start earlier, and we resolve a flattening of the Omori curve consistent with a larger Omori c ‐value for the farthest aftershocks, as predicted by Dieterich (1994) .

Seismological Research Letters↗
USGS2025

Characterizing the scale of regional landslide triggering from storm hydrometeorology

Rainfall strongly affects landslide triggering; however, understanding how storm characteristics relate to the severity of landslides at the regional scale has thus far remained unclear, despite the societal benefits that would result from defining this relationship. As mapped landslide inventories typically cover a small region relative to a storm system, here we develop a dimensionless index for landslide-inducing rainfall, A * , based on extremes of modeled soil water relative to its local climatology. We calibrate A * using four landslide inventories, comprising over 11 000 individual landslides over four unique storm events, and find that a common threshold can be applied to estimate regional shallow-landslide-triggering potential across diverse climatic regimes in California (USA). We then use the spatial distribution of A * , along with topography, to calculate the landslide potential area (LPA) for nine landslide-inducing storm events over the past 20 years, and we test whether atmospheric metrics describing the strength of landfalling storms, such as integrated water vapor transport, correlate with the magnitude of hazardous landslide-inducing rainfall. We find that although the events with the largest LPA do occur during exceptional atmospheric river (AR) storms, the strength of landfalling atmospheric rivers does not scale neatly with landslide potential area, and even exceptionally strong ARs may yield minimal landslide impacts. Other factors, such as antecedent soil moisture driven by storm frequency and mesoscale precipitation features within storms, are instead more likely to dictate the patterns of landslide-generating rainfall throughout the state.

California↗
USGS2026

Surface-wave relocation and characterization of the 2024–2025 dyking episode along the Fentale–Dofen segment of the Ethiopian rift

From late-December 2024 to mid-March 2025, a 50-km-long dyke intrusion triggered over 300 earthquakes (magnitude 4 to 5.9) between Fentale and Dofen volcanoes along the Northern Main Ethiopian Rift. Dyke intrusions periodically occur along the Fentale–Dofen magmatic segment and are an expression of ongoing rift extension. Preliminary analyses using interferometric synthetic aperture radar revealed extensive ground deformation (up to 60 cm), which closely matched the temporal and spatial evolution of surface manifestations and earthquake locations from global catalogues. While global catalogues are critical for real-time monitoring, the precision of locations in remote and or sparsely instrumented regions can be low. In this investigation, we present surface-wave relocation results of the dyking episode that began near Fentale volcano in December 2024. We estimate relative locations using differential traveltimes measured from regional-to-teleseismic distance surface-wave observations of earthquakes reported by the U.S. Geological Survey. Relative relocations reduce the initial region of diffuse seismicity to a 50-km-long narrow band bounding the strike of surface manifestations and the zone of maximum surface deformation. We demonstrate the precision of surface-wave relocations over incremental time periods, capturing the progression of dyking from seismic onset through seismic migration and caldera subsidence. Results showcase the utility of surface-wave relocations in the characterization of dyking episodes and provide complementary insights into the current understanding of the Fentale–Dofen volcanic plumbing system.

Geophysical Journal International↗
USGS2026

Moment magnitude for small earthquakes in the Delaware basin of west Texas and southeast New Mexico, USA

The Delaware Basin region of west Texas and southeast New Mexico has become one of the most prolific regions of seismic activity in the continental United States due to widespread hydraulic fracturing and wastewater disposal injection. In response to the increased number of earthquakes in this region, rapid and accurate characterization of earthquake sources is necessary to understand the evolution of seismic activity and level of seismic hazard associated with these earthquakes. This study re-evaluates earthquake magnitudes, estimating moment magnitude (MW) for small earthquakes in the Delaware Basin using 1) moment-rate spectra derived from S-wave coda envelopes, and 2) a relative magnitude method that relies exclusively on the ratio of waveform amplitudes between highly correlated waveform pairs. The coda-envelope method produces accurate M W estimates for small earthquakes ( M 1.5 – 3) that are consistent with independent, waveform modeled moment magnitudes for events with M W > 3 . Using the relative amplitudes method to extend these M W magnitudes to many other events, we successfully provide relative moment magnitude ( M W,rel ) values for 81% of the Texas Seismological Network catalog in the Delaware Basin region, and 45% of the USGS Induced Seismicity Project’s catalog of events in southeast New Mexico. The adoption and integration of the calibrated M W,rel method with current magnitude estimation methods offers valuable insights into the relationships between local and moment magnitude and will contribute to improved characterization of widespread induced seismicity.

New Mexico, Texas↗
USGS2025

Earthquake stress drop, source spectra from laboratory experiments, and the apparent unsteadiness of frictional slip

Over a range of co-seismic total energy release, high sample rate slip, stress, and slip velocity measurements were made during frictional failure to determine seismic source properties. The near-field measurements were used to estimate stress drops using spectral approaches that are analogous to the standard seismological moment-corner determined Δσ and acceleration amplitude Δσ arms stress drops. Rapid frictional sliding produces slip velocity amplitude spectra that vary with frequency -1 so long as the energy dissipated by friction approaches the total energy release. Under these conditions Δσ and Δσ arms are similar in size. Since these are nearly direct measurements of fault slip velocity and because they produce flat acceleration spectra between the corner and an apparent fmax, friction in these experiments seems to produce dynamic acceleration and deceleration at all resolved frequencies. The physical origin of this unsteady motion is not known exactly but likely reflects rapid dynamic variations in fault strength or applied stress. On the other hand, when the total energy stored prior to failure is artificially made somewhat larger relative to that which can be dissipated by friction the velocity spectra vary with frequency -2 immediately below the corner. At these conditions Δσ and Δσ arms diverge. The excess stored energy is preferentially partitioned into radiated energy upon failure. Collectively the experiments imply that to produce unsteady, white noise accelerations that are observed for natural earthquakes, requires that at least 95% of the energy released does not escape the source region to be radiated to the far-field.

Bulletin of the Seismological Society of America↗