Search USGSSearch

SEARCH · Search USGS

Results for “Remote Sensing Applications: Society and Environment”

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.

1,658 records · Page 7Linked to original sources

Turbidite correlation for paleoseismology

Marine turbidite paleoseismology relies on the assumption of synchronous triggering of turbidity currents by earthquake shaking to infer rupture extent and recurrence. Such inference commonly depends on age dating and correlation of the physical stratigraphy of deposits carried by turbidity currents (i.e., turbidites) across great distances. Along the Cascadia subduction zone, which lies offshore the Pacific Northwest, USA, turbidite facies in core photographs, X-ray computed tomography images, and magnetic susceptibility (MS) data exhibit differences in character over relatively short distances, which implies that not all deposits can be correlated with confidence. Thus, subjective correlation based on expected similarity over great distances and weak age constraints does not independently support paleoseismic models. We present a new method for correlating turbidites along the Cascadia margin that can yield a more objective and repeatable stratigraphic framework to underpin earthquake recurrence. We use dynamic time warping to correlate MS logs and measure correlation coefficients of core pairs to evaluate correlation strength. We then compare these measures to a distribution of correlation coefficients of randomly generated turbidite sequences and find that only a small number of core pairs can be correlated more confidently than randomly stacked turbidites. This methodology promises a more robust correlation strategy for future stratigraphic studies.

Oregon, Washington

Testing characteristic magnitude distributions in modern PSHA models

The characteristic magnitude distribution hypothesis predicts a higher rate of large earthquakes than a Gutenberg–Richter extrapolation of the small‐earthquake rate would imply. Characteristic magnitude distributions have been commonly applied to faults in probabilistic seismic hazard analysis (PSHA), and in modern models they can emerge from the way short‐term seismicity constraints are combined with long‐term geologic and geodetic constraints. We test the characteristic magnitude distribution hypothesis by comparing the fault‐based magnitude distributions from the 2023 update to the National Seismic Hazard Model (NSHM23) in the Western United States with observed seismicity over the past 93 yr. We find that observed magnitude distributions fall outside the model‐predicted confidence bounds in regions where NSHM23 produces characteristic magnitude distributions: in these regions, the model predicts higher rates of large earthquakes than are observed. An analysis of the earlier California model (Uniform California Earthquake Rupture Forecast, version 3) also reveals discrepancies between the modeled and observed magnitude distributions. In addition, we find that observed magnitude distributions near modeled faults are not significantly different from those in background regions. These results challenge the prevalence of characteristic magnitude distributions in fault‐based seismic hazard models and call for a reassessment of how disparate data sets are integrated in PSHA.

western United States

Use of precariously balanced rocks to constrain postglacial earthquake magnitudes in New England, United States

We use the age and fragility of 13 precariously balanced rocks (PBRs) in Massachusetts, New Hampshire, and southwestern Maine to evaluate the estimated probabilistic seismic hazard (PSH) from the U.S. Geological Survey’s National Seismic Hazard Model (NSHM), and also develop a map of maximum postglacial earthquake magnitudes in the region. The PBRs are erratic boulders left behind by the retreat of the Laurentide ice sheet between 13 and 18 ka, and published ice retreat maps are used to estimate their fragility ages. The fragilities and ages of the rocks appear to be consistent with the NSHM, in that the estimated earthquake magnitudes and resulting ground motions that would topple them are only expected on timescales longer than their fragility ages. This result, on the North American passive margin, contrasts with PBR studies at active plate boundaries, where the PBR constraints imply large reductions to estimated PSH. It may be that the PBRs are not old enough in the eastern United States to reveal the same PSH modeling issues as those revealed in the plate boundary settings. The study provides constraints on maximum magnitudes of postglacial earthquakes across much of Massachusetts, New Hampshire, and southwestern Maine, with the constraints limited to the ages of the PBRs.

Maine. Massachusetts, New Hampshire

Temporal and spatial changes in seismic attenuation associated with inferred fluid migration in the 2016 central Apennines earthquake sequence

Prior work suggests that high‐frequency seismic attenuation acts as a highly sensitive proxy for crustal permeability and fluid mobility in fractured media. We test the hypothesis that the fault system responsible for the 2016–2017 Amatrice–Visso–Norcia–Capitignano sequence acted as an impermeable seal, compartmentalizing pressurized fluids until dynamic rupture triggered widespread fluid diffusion. By tracking across the sequence the spatiotemporal evolution of the S ‐wave anelastic attenuation parameter, we identify large, positive low‐frequency attenuation anomalies emerging within the hanging wall following the Amatrice mainshock and strictly preceding subsequent large ruptures. Conversely, we observe weaker, negative anomalies in the footwall, anticorrelated in time with those of the hanging wall, revealing a massive asymmetry in fluid redistribution and permeability evolution across the fault system. Furthermore, aftershock migration rates reveal distinct linear alignments in a distance‐reduced time space, allowing us to explicitly track and quantify episodes of lateral and upward fluid migration. These physically consistent patterns suggest that stress‐driven fluid diffusion directly weakens adjacent fault patches, dictating the spatiotemporal migration of seismicity. We conclude that near‐real‐time monitoring of seismic attenuation may help detect fluid redistribution in active fault systems and may provide useful information for time‐dependent seismic hazard assessment.

central Apennines

High‐resolution surface deformation and slip distribution observations for the 2023 Kahramanmaraş, Türkiye, earthquake sequence help constrain the rupture process

Splay, or branch, faults are a common geometric feature of earthquake surface ruptures and may provide constraints on the rupture behavior of an earthquake. The 2023 M w 7.8 Pazarcık and M w 7.5 Elbistan, Türkiye, earthquakes are examples of ruptures with multiple small splays, and the Pazarcık earthquake nucleated on a splay fault, the Narlı fault, before rupturing bilaterally on the East Anatolian fault (EAF). Here, we present 3‐m‐resolution surface displacement from subpixel correlation of Planet Dove optical images for the entirety of both ruptures with corresponding surface slip distributions. For a 30‐km‐long study region spanning the Narlı‐EAF intersection, we compare surface slip derived from five data sets with different resolutions (on‐the‐ground, WorldView, Planet Dove, Sentinel‐2, and Sentinel‐1) to elucidate complementary information. In addition, we integrate information from the surface expression of faulting with published dynamic rupture simulations and rupture process studies to constrain a rupture evolution for the Pazarcık earthquake that is consistent across data sets. This work highlights the complementary nature of disparate surface slip data sets and the role that high‐resolution surface displacement information, including from fault splays, can play in constraining nonunique rupture models and refining understanding of the earthquake rupture process.

Bulletin of the Seismological Society of America

Characterizing directivity in small (M 2.4-5) aftershocks of the Ridgecrest sequence

Directivity, or the focusing of energy along the direction of an earthquake rupture, is a common property of earthquakes of all sizes and can cause increased hazard due to azimuthally dependent ground‐motion amplification. For small earthquakes, the effects of directivity are generally less pronounced due to reduced rupture size, yet the directivity in small events can bias source property estimates and provide important insights into general regional faulting patterns. However, due to observational limitations, directivity is usually only measured and modeled for large events. As such, many studies of small earthquakes either ignore directivity altogether or assume a constant rupture direction for all events in a cluster. In our study, we apply a refined directivity fitting method constrained with two separate methods of source deconvolution to the dataset of aftershocks of the 2019 Ridgecrest earthquakes, which contain a large number of well‐recorded small‐to‐mid sized earthquakes occurring in close proximity to each other. The revealed directivity of 100+ small (M 2.4–5) earthquakes is highly heterogeneous and primarily oblique to and away from the main fault strike, suggesting a complex postseismic stress redistribution. In addition, the energy focusing effect of directivity appears to bias the selection of high‐quality data from stations in the direction of rupture, leading to average stress‐drop increases of 50% if directivity is not accounted for.

California

Late Jurassic igneous rocks in south-central Arizona and north-central Sonora: Magmatic accompaniment of crustal extension

In Middle Jurassic time, the region of south-central Arizona and north-central Sonora was part of a continental magmatic arc. In early Late Jurassic time, this magmatic arc gradually gave way to a extensional regime. Late Jurassic rifting was accompanied by magmatism, typically small in volume relative to the preceding arc phase, and with mildly alkaline compositions. These igneous rocks, the 158- to 146-Ma Ko Vaya Suite, make up much or most of the Comobabi, Artesa, Quijotoa, and Brownell Mountains and Sierra del Cobre. The bimodal Ko Vaya Suite constitutes volcanic to shallow plutonic complexes comprising chiefly trachyandesite and volcanic wacke; monzodiorite; rhyolite porphyry; distinctive, compositionally and texturally heterogeneous, quartzpoor granite (the Ko Vaya Granite); A-type perthite granite and quartz syenite; aphyric intrusive rhyolite; and locally abundant hematite veins. Alteration, particularly potassic alteration, is widespread. Minimally altered samples of granite and quartz syenite have the composition of trachyte or trachydacite and alkali rhyolite. The Ko Vaya Suite is nonconformably overlain by the latest Jurassic Sand Wells Formation, approximately correlative with the Glance Conglomerate, the basal unit of the Bisbee Group. The Ko Vaya Suite is in some ways analogous to the coeval Independence dike swarm of eastern California, but less diffuse. We infer that the Ko Vaya Suite represents especially intense or focused Late Jurassic extensional magmatism.

Arizona, Sonora

Top-down targeted network analysis of critical mineral commodities applied to international geochemistry database

The global demand for critical mineral commodities is rapidly increasing, making domestic production an important factor in supporting the economy and national security. Large scale, publicly available geochemical databases allow for the application of data informatics methods to interrogate critical mineral commodities data for correlations in deposit formation and distribution, particularly for identifying enrichment of multiple critical mineral commodities at the same deposit. In this study, we applied network analysis to the Critical Minerals Mapping Initiative (CMMI) ore geochemistry (Critical Minerals in Ores, CMiO) database to identify the high concentration (defined as 10× bulk crustal abundance) co-occurrence of different critical mineral commodities across a mineral system hierarchy from deposit environments to individual deposits. Identifying patterns or unique outliers in enrichment in network communities will allow for the location of secondary critical mineral commodity resources from under-utilized deposits. We find trends in the enrichment of critical mineral commodities in network-communities between the elements praseodymium (Pr), neodymium (Nd), terbium (Tb), and dysprosium (Dy) across multiple CMiO database deposit environments and groups down to specific deposit types and sites. A separate trend in network community deposition is observed as well between iridium (Ir) and platinum (Pt) in deposit environments, groups, types, and sites. Network analysis focused on critical minerals in magmatic-hydrothermal deposits identified multiple deposit sites from different deposit types within the CMiO database with concentrations of Dy, Nd, Tb, Pr, Ir, and Pt that are at least ten times greater than the crustal average. This approach can be applied to any target element(s) or deposit(s) of interest, allowing broad investigation of co-enriched critical mineral commodities.

Journal of Geochemical Exploration

Uranium—Deposits, production and resources, market dynamics, and supply chain risks

Introduction Interest in nuclear power for the generation of electricity has risen with the increase in the need for more diverse baseload power, enhanced energy security, and the development of new technologies, such as small modular reactors (SMRs), which could provide power for remote areas, industrial applications, and artificial intelligence (AI) data centers. In 2024, the U.S. Department of Energy received $2.7 billion in congressional funding to bolster the domestic uranium production and nuclear fuel supply chain and address reliance on imports from foreign suppliers. In 2025, the U.S. Government issued several Executive and Secretary’s orders aimed at revitalizing the U.S. nuclear sector. If SMRs are to be as widely deployed in the United States and worldwide as envisioned, demand for uranium (nuclear reactor fuel) will likely increase. After the Fukushima nuclear accident in 2011, the market spot price of uranium began a decline, followed by a decrease in U.S. and global uranium exploration and mine development expenditures that led to a uranium supply deficit until 2020, when prices started to recover, prompting a resurgence in uranium exploration and development. In January of 2024, the uranium spot price rose to a 17-year high $106 (U.S. dollars) per pound of U3O8 (triuranium oxide, commonly known as “yellowcake”), which is expected to increase uranium exploration, mine development, and uranium production domestically and worldwide.

Fact Sheet

‘The fish that stop’: Drivers of historical decline for Pacific cod and implications for modern management in an era of rapidly changing climate

n the Gulf of Alaska, a series of marine heat waves depleted Pacific cod ( Gadus macrocephalus ) biomass to the lowest abundance ever recorded and led to the fishery’s closure in 2020. Although the fishery has been productive for decades, this collapse may have historical precedents. Traditional knowledge holders refer to cod as ‘the fish that stop’, and there is a suggested period of decline in the 1930s. Here we conduct a catch reconstruction of the early commercial fishery (1864–1950), confirming a rapid catch decline in the 1920s and 1930s. Next, we evaluate evidence for possible drivers. We document changes to demand and technology that contributed to declining catch. However, we also find both qualitative and quantitative evidence of depletion, suggesting catch declines were not driven entirely by social factors. Overfishing may have contributed to localized catch declines as evidenced by declining catch rates in heavily fished localities. We also find evidence for climate as a driver of regional decline, with the period of catch decline characterized by up to 2°C higher temperatures as compared to the earlier period of high fisheries production. Our analysis underscores the importance of understanding long-term drivers of fisheries productivity and the value of linking fisheries and climate histories.

Alaska

Constraining source and path effects of large magnitude earthquakes using ground motion simulations

The purpose of this study is to use ground‐motion simulations to investigate ways in which source and path effects for large‐magnitude earthquakes can be represented in nonergodic ground‐motion models (GMMs). To achieve this, we designed a ground‐motion study in the San Francisco Bay Area that includes earthquakes with a broad range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of kinematic simulations (magnitude 4–7), we then develop a nonergodic GMM with the simulated data. We find that trends in the within‐site residuals are affected significantly by the earthquake radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within‐site residuals among multiple source realizations. Finally, after removing the source effects from the within‐site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site as the large events, whereas the second approach considers all small events on the fault plane. The results indicate that it is difficult to satisfactorily approximate the path effects of large events with those of small events using either approach, at least in the case of simulations.

California

Near-fault amplification and ground motion variability during the 2019 Ridgecrest, California sequence

We estimate ground-motion variability near the 2019 M 7.1 Ridgecrest earthquake sequence. Accurate seismic hazard estimation requires understanding ground-motion spatial correlations, yet many studies lack the dense station coverage needed to resolve small-scale variability. The 2019 M 7.1 Ridgecrest earthquake sequence presents a unique opportunity to examine ground motions and their spatial correlations at a range of interstation distances. The permanent seismic network was augmented with hundreds of temporary stations including several fault-crossing nodal arrays. We compute the event ( δE i ) and within-event ( δW ij ) residuals from the observed peak ground velocity and peak ground acceleration data to isolate potential sources of ground-motion variability. We then compare δW ij between station pairs that record an event to understand the semivariance of the ground motion versus interstation distance. By fitting an exponential model to the semivariances, we determine a correlation range of 25 km for the Ridgecrest region. Although the exponential model fits the broad-scale increase of semivariance with interstation distance, we also observe smaller-scale trends. We find that ground motions are less correlated for station pairs that are near or across faults that ruptured during the 2019 Ridgecrest sequence. We also find large, positive median δW ij with relative values 2–3 times larger than nearby stations for individual stations’ near-fault traces. Near-fault amplification and greater ground-motion variability can delineate fault zones and may locally increase the seismic hazard.

California

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface

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

Ground motions from finite faults with unknown geometry

Earthquake emergency response and loss modeling require accurate models of ground‐motion intensity metrics (IMs), such as those produced by ShakeMap. In the critical minutes and hours after an earthquake, however, knowledge of the earthquake source geometry limits the accuracy of these models, especially in areas where seismic instrumentation is insufficient to constrain the shaking intensity pattern. Previous attempts to overcome this limitation have centered on computing the mean distance metrics from a suite of assumed source geometries. That approach is inadequate with modern ground‐motion models (GMMs) because they typically include terms for the hanging wall effect, depth to top of rupture, and other nonlinear terms. Such approaches also propagate the distance uncertainty through the GMMs and may not properly allocate the additional uncertainty into the between‐ and within‐event components. Here we present Finite‐Fault SIMulation with Mixed Effects Regression (FFSIMMER), a method to address these limitations by computing IM distributions from a suite of rupture realizations that fully accounts for all geometry‐related GMM terms and allows for the separation of the between‐ and within‐event components of the total uncertainty. The accuracy of this method is limited by the available constraints on the source parameters. To address this, FFSIMMER can incorporate constraints based on knowledge of the source, such as the approximate strike and dip based on knowledge of the plate tectonics and focal mechanism of the earthquake. We demonstrate the improvement in accuracy that can be achieved with the 2025 M 7.7 Mandalay, Burma (Myanmar), earthquake by incorporating constraints on the range of plausible source parameters.

Bulletin of the Seismological Society of America

Generalized Bancroft algorithm for locating earthquakes with P- and S-wave arrival times

Because of similarities between locating an earthquake with seismic stations and locating a Global Positioning System (GPS) receiver from satellites, the Bancroft algorithm developed for GPS processing can be used to locate earthquakes. Such an approach to earthquake location differs from the conventional method of choosing an initial or trial solution and then iteratively improving the solution until convergence. The Bancroft algorithm has the advantage of being a direct, noniterative solution but with the disadvantage of only being able to accommodate a homogeneous velocity model. An additional limitation of the standard Bancroft algorithm is that it considers arrival times in a medium with a single propagation velocity. This poses no problem for GPS processing because electromagnetic waves travel at the speed of light; however, for seismic waves it means the algorithm can be applied to collections of either P ‐ or S ‐wave arrival times. Here, I show how the Bancroft algorithm can be generalized to handle both P ‐ and S ‐wave arrival‐time measurements simultaneously. I also show how to accommodate depth‐varying P ‐ and S ‐wave velocity models. I apply the generalized Bancroft algorithm to microearthquakes beneath Tanaga Volcano in Alaska and compare standard locations from the widely used HYPOINVERSE location code to Bancroft locations and to the output of HYPOINVERSE when setting the trial location to the Bancroft location. I find the Bancroft locations outperform the results from the other methods for shallow earthquakes near sea level, where a quantity known as the geometric dilution of precision is large and linearized approaches such as HYPOINVERSE are expected to struggle.

Alaska

Site response in the Walnut Creek–Concord region of the San Francisco Bay, California: Ground motion amplification in a fault-bounded basin

Thirty‐seven portable accelerometers were deployed in the eastern San Francisco Bay communities of Walnut Creek and Concord to study site response in a fault‐bounded, urban, sedimentary basin. Local earthquakes were recorded for a period of two years from 2017 to 2019 resulting in 101 well‐recorded events. Site response is estimated by two methods: the reference site spectral ratio method and a source‐site spectral inversion method. The reference site spectral ratio method allows investigation of the variability of site amplification with source azimuth and frequency. The source‐site spectral inversion method yields the best least‐squares fit to site response for a database of ground‐motion records. Both methods show substantial amplification in the Walnut Creek–Concord basin below 2 Hz indicating strong surface‐wave development. Greater amplification is seen for sources aligned along the long axis of the basin. Inversion using close‐in sources at short distances yields lower amplification at longer periods than the entire data set due to reduced surface‐wave generation for steeper angles of incidence. Inversion of site response spectra for shallow shear‐wave velocity using a global search algorithm yields V S30 values consistent with generalized mapping results based on geology and topography but with greater variability due to local site variations. 3D finite‐element modeling shows greater amplification in the Walnut Creek–Concord basin with a basin‐edge effect likely contributing to higher ground motions. Topography is also seen to lead to increased scattering and shadowing effects.

California

High frequency and region-scale simulations of large (Mw7+) earthquakes on the southern Whidbey Island fault, Washington, USA

We simulate ground shaking in western Washington State from hypothetical M w 7.0–7.5 earthquakes on the southern Whidbey Island fault (SWIF). Ground motions are modeled considering kinematic source distributions on a complex fault plane, a 3D seismic velocity model, and region‐specific soil velocity models. We run simulations with varying model resolutions, including regional‐scale simulations with a maximum‐modeled frequency of ∼1 Hz and local‐scale simulations with a maximum‐modeled frequency of ∼2.5 Hz. Additional local‐scale simulations are run considering high‐resolution surface topography. We explore how source parameters (i.e., magnitude, hypocenter location, and dip direction) and 3D velocity structure impact peak shaking intensity and its variability. In particular, we find that earthquakes on the SWIF would likely produce strong shaking throughout the populated Puget Lowland, including in the cities of Everett, Seattle, Bellevue, and Tacoma, Washington. Simulated short‐period ( T ≤ 2 s) spectral accelerations are strong throughout the Puget Lowland, and long‐period shaking ( T ≥ 5 s) is strong in the deep regional sedimentary basins, especially the Everett and Seattle basins. Source parameters strongly influence intra‐ and interevent variability in response, primarily through changes in source and site geometry, as well as rupture directivity. We also note a potential coupling between rupture directivity and basin effects, wherein directivity pulses are seemingly guided through the region’s deep, interconnected sedimentary basins. Overall, this work highlights the impacts of 3D source, path, and site effects on seismic hazard in the U.S. Pacific Northwest and substantially expands the catalog of simulated ground motions for Puget Sound area crustal faults.

Washington