USGS Science⌕ Search

FIND YOUR NEXT DISCOVERY

Results for “Earthquake Spectra”

Original records, connected by a shared subject.

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

266 records · Page 2Linked to original sources
USGS2025

Earthquake stress-drop values delineate spatial variations in maximum shear stress in the Japanese forearc lithosphere

Earthquake stress drop (Δσ) may increase with depth and stress in the brittle lithosphere. However, the range of uncertainty in Δσ and the lack of constraints on absolute stress make it difficult to establish whether they are correlated. Here, we investigate Δσ dependence on depth and maximum shear stress ( τ max ) based on ~11 years of seismicity in the northeastern Japanese forearc following the 2011 Tohoku-Oki megathrust earthquake. We interpret Δσ estimates computed using both individual spectra and spectral-ratio methods and find that Δσ exhibits a clear depth dependence within the seismically active upper ~60 km of the forearc lithosphere ( ~ 0.8 MPa per 10 km). We further compare Δσ values with quantitative τ max estimates from finite-element models of force balance. We find that median Δσ values increase with τ max in the brittle forearc lithosphere and that earthquake stress release is proportional to τ max . The dependence of Δσ on τ max explains the apparent depth dependence of Δσ and suggests that average Δσ values provide a relative measure of the stress at failure. In the northeastern Japanese forearc, Δσ values remained roughly constant in the decade following the Tohoku-Oki earthquake, suggesting negligible changes in failure stress in the forearc since the mainshock.

Communications Earth and Environment↗
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

Magnitude, depth and methodological variations of spectral stress drop within the SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest Earthquake Sequence

We present the first ensemble analysis of the 56 different sets of results submitted to the ongoing Community Stress Drop Validation Study using the 2019 Ridgecrest, California, earthquake sequence. Different assumptions and methods result in different estimation of the source contribution to recorded seismograms, and hence to the source parameters (principally corner frequency, f c ⁠ , spectral stress drop, Δσ, and seismic moment, M 0 ⁠ ) obtained from modeling calculated source spectra. For earthquakes smaller than magnitude (M) 2.5 there is negligible correlation between the f c values obtained by different studies, implying that no present method is reliable using available data. For larger magnitude events, correlation between f c measurements of different studies, within even a small M range is always higher than spectral ⁠Δσ , because the f c measurements simply reflect the underlying physical decrease in f c with increasing M. We model the observed trends of submitted f c with both magnitude and depth. Most methods report an increase in spectral Δσ with M, although a magnitude‐invariant spectral Δσ is within the confidence limits. The depth dependence is smaller and depends on whether a study allows attenuation to vary with source depth; a combination of depth‐dependent attenuation correction, and depth‐dependent shear‐wave velocity can compensate for reported depth trends. We model the submitted values to remove differing M and depth variation to investigate the relative interevent variability. We find consistent relative variation between individual events, and also lower relative spectral Δσ in the northwest of the aftershock sequence, and higher on the cross fault and in the region of main fault intersection. This large‐scale comparison implies that absolute spectral Δσ estimates are dependent on the methods used; studies of different regions or using different methods should not be directly compared and improved constraints on path and site corrections are needed to resolve these absolute spectral Δσ differences.

California↗
USGS2025

Earthquake magnitude and source parameter estimation with a distributed acoustic sensing dataset in the Gorda subduction zone

Distributed acoustic sensing (DAS) systems offer a cost‐effective way to create large‐scale strainmeter arrays for seismological applications using fiber‐optic cables. DAS‐based strain measurements are known to be influenced by various factors, bringing into question their general reliability for accurate earthquake characterization. A 15‐km‐long DAS deployment in northern California was operational within 3 days of the 2022 M w 6.4 Ferndale earthquake and ran continuously throughout the aftershock sequence. We utilize these aftershock data to validate DAS‐based strain measurements in two ways. We first test the accuracy of DAS‐based magnitude estimates from peak dynamic strains by comparing them with magnitude and attenuation scaling relations derived independently from traditional borehole strainmeter (BSM) data. We demonstrate that DAS‐based magnitudes are comparable to BSM‐based magnitudes when corrections for variations in site response along the fiber‐optic cable are properly made. Magnitude errors are spatially correlated, potentially because of factors such as finite‐fault effects (e.g., stress drop) or more complex, unmodeled path attenuation or because of wave propagation effects in heterogeneous media. We then apply more advanced source characterization methodology to the DAS data using a time‐domain empirical Green’s function (EGF) deconvolution approach to measure details of the moment rate history. The EGF approach using DAS data depends on careful treatment of distorting factors such as anthropogenic sources of noise and optical phase wrapping but successfully isolates source spectra for moderate‐magnitude earthquakes: source spectral ratios obtained from DAS data, broadband seismometer data, and BSM data in the same region show consistent results, revealing differences in directivity and spectral shape among earthquakes. Although further research is needed to refine source‐time‐function estimation techniques for DAS data, particularly for larger magnitude events, these case studies demonstrate the clear potential of DAS for earthquake source characterization.

California↗
USGS2026

Efficient physics‐informed ground‐motion simulations with reduced‐order models: CyberShake implications and high‐resolution site terms for southern San Andreas fault earthquakes

Recent advances in Probabilistic Seismic Hazard Analysis (PSHA) leverage physics‐based ground‐motion simulations to estimate seismic hazard, such as the CyberShake project. However, computational costs quickly escalate when performing PSHA for numerous faults or sites and can become prohibitively expensive. To reduce computational demands, CyberShake uses reciprocity and interpolates physics‐informed corrections from simulations conducted at fewer locations, but the accuracy of these interpolations remains poorly quantified. To quantify the interpolation accuracy, we derive high‐resolution, frequency‐dependent site terms for southern California and compare them with interpolated site terms using the CyberShake approach. We accomplish this by performing a set of earthquake point‐source simulations distributed along the nonplanar fault geometry for the southern San Andreas fault (SSAF) extending from Bombay Beach to Lake Hughes. Using SeisSol, we simulate three minutes of viscoelastic seismic wave propagation for these sources and store the horizontal‐component Green’s functions for 480,000 sites. We then use a scientific machine learning approach based on interpolated proper orthogonal decomposition to construct an accurate reduced‐order model of the Green’s functions to efficiently predict effective amplitude spectra (EAS) for finite‐source rupture models of SSAF earthquakes. Using minimum curvature interpolation with tension, as used in CyberShake, we compare the interpolated site terms against our high‐resolution site terms. We identify local discrepancies with EAS differing by up to a factor of approximately three. Furthermore, we identify locations where unexpectedly high or low ground motions are missed when using the interpolated dataset for these earthquakes. We estimate that our approach may be used within CyberShake to reduce the time‐to‐solution by a factor of 336 for the entire earthquake rupture forecast. Our analysis of physics‐based site terms provides more insight into the seismic hazard due to SSAF ruptures and guides future developments by combining high‐performance computing and reduced‐order modeling techniques for PSHA.

California↗
USGS2026

Remote compositional analyses of space-weathered lunar maria

Visible-to-shortwave infrared (VSWIR) reflectance spectroscopy has revolutionized our understanding of planetary surface compositions. However, space-weathering processes on airless bodies complicate quantitative compositional analyses. Here, we present a framework to isolate the signatures of space weathering in VSWIR spectra of lunar maria by leveraging radiative transfer modeling under the assumptions that (i) a space-weathered target can be expressed as a mixture of fresh and fully space-weathered components and (ii) remaining signatures can be modeled by including agglutinates as an end-member component. We first validate this approach against laboratory spectra of space-weathered Apollo mare soils of known mineral compositions using a probabilistic Markov Chain Monte Carlo implementation of the Hapke radiative transfer model. Second, we illustrate how this approach can be applied to orbital Moon Mineralogy Mapper data. The proposed space-weathering correction workflow for lunar maria could be expanded to other lunar lithologies and applied to existing and future data sets.

Planetary Science Journal↗
USGS2025

High-pass corner frequency selection and review tool for use in ground-motion processing

Raw seismological waveform data contain noise from the instrument’s surroundings and the instrument itself that can dominate recordings at low and high frequencies. To use these data in ground‐motion modeling, the effects of noise on the signals must be reduced and the signals’ usable frequency range identified. We present automated procedures to efficiently reduce low‐frequency noise that are implemented in the software package gmprocess. These procedures check for, and as needed remove, low‐frequency artifacts in the displacement record using polynomial fits, which can be used in combination with existing signal‐to‐noise ratio (SNR)‐based corner‐frequency selection procedures. The automated selections are then efficiently verified and refined using a graphical user interface (GUI) that plots relevant ground‐motion time series and spectra and tracks modifications to signal processing parameters. We demonstrate these procedures using recordings from the 2020 M 5.1 Sparta, North Carolina, and the 2013 M 4.7 southern Ontario earthquakes. Data processed with the SNR‐only and polynomial criteria for these events contain displacement artifacts in 37% and 23% of processed traces, respectively. Records with remaining artifacts are corrected manually using the GUI. These processing steps illustrate the workflow for efficient data processing with quality control.

Seismological Research Letters↗
USGS2025

Site response and wave propagation effects in the eastern United States

Fourier amplitude spectra from regional earthquakes in the eastern United States are used in a parametric inversion for source, path, and site effects. Five earthquakes are selected for analysis during the installation of the United States National Seismic Network (US), Earthscope’s USArray Transportable Array (TA), and other temporary arrays to maximize station coverage. A global search algorithm is used to solve for site response from 0.1 to 15 Hz, corner frequency, geometrical spreading ( r - γ ), and frequency dependent anelastic attenuation in the form Q(f) = Q o f α . Tradeoff between moment and geometric spreading is handled by fixing the moment. The tradeoff between corner frequency and Q(f) is solved by selecting the value of corner frequency that minimizes an objective function defined over all stations. Values of site response and attenuation parameters show a strong spatial correlation with the physiographic provinces of the eastern United States. Site response for the Atlantic Coastal Plain is consistent with previous work using spectral ratios relative to a reference site, defined by strong resonance peaks correlated with the thickness of sediments. Site response for the other physiographic provinces is markedly different from the coastal plain, with a lack of distinct resonance peaks and a broad moderate high at frequences from 0.1 to 0.5 Hz consistent with the hard-rock geology of the regions. Like site response, Q(f) has a strong correlation with physiographic province, showing lower values on the coastal plain and higher values inland. Geometric spreading exponent, γ, decreases with increasing hypocenter distance from just above 1 at a few tens of kilometers to 0.9 at 500 km. The limited range in geometric spreading values is attributed to starting the Fourier transform window at the S ‐wave arrival for all distances and averaging over multiple wave types.

eastern United States↗
USGS2026

A methods framework for evaluating measurement consistency across spectrometers for multispectral uncrewed aerial system vegetation mapping applications

The U.S. Geological Survey collects remote sensing data to support national scientific assessments of natural resources, hazards, and landscape change. Spectrometers and spectroradiometers are essential for gathering point-based spectral measurements used in applications such as uncrewed aerial systems (UAS) multispectral image calibration, validation, and analysis. Evaluating how different instruments perform in laboratory and field environments helps determine whether they provide consistent, interoperable measurements. Such verification can expand access to spectral ground data during UAS operations by allowing scientists to use alternative instruments when budgets, logistics, or field conditions limit options. We propose and test a methodological framework for evaluating spectrometers for measurement consistency during UAS multispectral vegetation mapping applications. There are three central evaluation components to the framework: laboratory, field, and relative to UAS multispectral imagery. By evaluating the instruments in both relatively controlled and uncontrolled environments, we thoroughly examine measurement consistency and when/why measurements may differ. We opportunistically selected two instruments for a case study in a coastal marsh setting: a compact laboratory spectrometer we modified for field use and a field-ready spectroradiometer. The instruments produced consistent measurements in both environments. We found differences between the field spectra and UAS spectra that likely reflect the perspectives of ground vs. aerial data and indicate that further radiometric calibration may be needed.

Massachusetts↗
USGS2026

Apatite (U-Th)/He thermochronology from Marsh Creek anticline reconciles Cenozoic and Holocene strain patterns and elucidates the Cenozoic canning displacement zone of Arctic Alaska

The northeastern Brooks Range is a vast fold-thrust belt that records the northernmost expression of Cenozoic deformation in Alaska, likely with complex dynamic linkages to distant tectonic elements. The principal focus of this manuscript is the Marsh Creek anticline (MCA), which is a large and recently active structure within the broader northeastern Brooks Range. The MCA is a complex, polyphase structure mostly buried beneath the Arctic coastal plain. We present apatite (U-Th)/He (AHe) dates from six outcrop samples, all near the structural crest of the anticline, as well as supporting detrital zircon U/Pb (ZUPb) age spectra from five outcrop samples. The AHe sample transect spans ~3 km of structural relief. The four shallowest AHe samples are not reset. However, the two deepest samples exhibit low intrasample variability among replicates and strong date versus grain radius correlations, suggesting post-depositional burial heating to within the AHe partial retention zone. Inverse thermal history modeling of these data resolves a protracted phase of exhumational cooling that began around 33–28 Ma. Detrital ZUPb age spectra from Paleogene strata support an interpretation of diverse inherited thermal histories recorded by the AHe data, and they are also consistent with extant reconstructions of Paleogene sediment dispersal. Placed into regional context, the results constrain a phase of increasing structural relief in the northeastern Brooks Range fold-thrust belt that began by the Eocene. The relief developed in a zone characterized by diminishing structural relief to the west called the Canning displacement zone. Although this zone was active since the Eocene, several structures including the studied part of the MCA were rapidly exhumed in the Oligocene. Lastly, we argue for kinematic compatibility between Cenozoic contraction and active strike-slip in the northeastern Brooks Range, with the temporal change in strain style likely reflecting the increase in structural relief and topography.

Alaska↗
USGS2008

Geology and geochemistry of Jurassic plutonic rocks, Baboquivari Mountains, south-central Arizona

Among the plutons of the Jurassic magmatic arc segment in southern Arizona, northern Sonora, and southern California, two distinctive rock types dominate, occurring together in numerous mountain ranges: porphyritic, titanite-bearing hornblende-biotite granodiorite; and biotite leucogranite. These characteristic rocks are particularly well preserved and well exposed around Kitt Peak, in the Baboquivari Mountains, southern Arizona. The Middle to early Late Jurassic Kitt Peak Plutonic Suite (KPPS) comprises three units: Aguirre Peak Quartz Diorite (APQD;<170 Ma, >=165 Ma, U-Pb), chiefly hornblende mesodiorite; Kitt Peak Granodiorite (KPGD; 165 Ma), hornblende-biotite granodiorite and monzogranite; and Pavo Kug Granite (PKG; 159 Ma), equigranular biotite leucogranite. This range also hosts a fourth, regionally unusual, Jurassic plutonic unit: the Baboquivari Peak Perthite Granite, related to post-arc, Late Jurassic crustal extension. Binary variation diagrams, incompatible element patterns, and REE spectra indicate APQD and KPGD are consanguineous. APQD basaltic mesodiorite evidently approximates the parental magma for the APQD–KPGD series. More evolved rocks of this series evidently formed by combinations of fractional crystallization and mixing with or assimilation of crustal material. PKG is less directly related to APQD and KPGD than these two are to one another.

Arizona, Sonora↗
USGS2024

Effective site coefficients for the 2024 International Building Code (IBC)

The U.S. National Seismic Hazard Models (NSHMs), developed by the U.S. Geological Survey (USGS), have long been the scientific foundation of seismic design guidelines and have been used to compute design ground motions for construction of new buildings and retrofit of existing buildings in the United States and its territories. The 2018 NSHM is adopted by the 2024 International Building Code (IBC). Prior to the 2018 NSHM update, hazard calculations were required only at one reference site condition defined by V S30 =760 m/s (where V S30 is the time-averaged shear wave velocity from the surface to a depth of 30 m) and three periods (peak ground acceleration, PGA, and pseudo spectral accelerations at periods of 0.2 and 1 s, S s and S 1 ). Site coefficients, F PGA , F a , and F v , were then defined by the Building Seismic Safety Council (BSSC) Provisions Update Committee (PUC) in the site-specific procedures of National Earthquake Hazard Reduction Program (NEHRP) Recommended Seismic Provisions to calculate ground motions for other site classes with different V S30 values at the given periods. The design ground motions at other periods were then estimated using a generic spectral shape that was also defined by the BSSC PUC in NEHRP provisions. In recent years, the engineering community has realized there were deficiencies with the F a and F v site coefficients and the generic spectral shape. To avoid potentially dangerous underestimations of design ground motions for long period structures on soft site conditions in high seismicity regions, the BSSC PUC recommended the use of multi-period response spectra (MPRS) in 2017. As a result, the USGS produced multi-period and multi-V S30 response spectral values in the 2018 NSHM for calculations of design ground motions and the site coefficients F a and F v were eliminated from the 2020 NEHRP Provisions. As these site coefficients were widely used inside and outside of the United States, in this study we back-calculate the “effective” site coefficients F a,eff , and F v,eff by comparing MPRS for various site classes with the MPRS for the reference site condition, and discuss the changes that are observed in the 2024 IBC compared to its previous version in 2021. The effective site coefficients are presented for test site locations and their dependence on various factors including period, ground motion intensity, and regional models are discussed. Ratio maps between the new effective site coefficients and the old ones are then presented for soft site classes and for short and long periods. For soft site classes at short periods, the new effective site coefficients are lower than the old site coefficients for high seismicity regions and higher for low seismicity regions. As it was expected, for soft site classes at long periods and high seismicity regions, the new effective site coefficients are much larger than the old site coefficients without imposing the 50% increase as a penalty that was suggested in the 2021 IBC, whereas they could be much smaller if the 50% increase would have been imposed particularly around New Madrid and Charleston high seismicity regions. For low seismicity regions, the long period effective site coefficients can be smaller or larger by 20% compared to the 2021 IBC coefficients.

Conference Paper↗
USGS2025

Using peak geometry and shifts in the x-ray spectrum of carbon from electron probe microanalysis to determine thermal maturity of organic matter

During the burial of mudstones, the associated organic matter undergoes gradual thermal maturation, a key process that can influence the reactivity of organic matter during catagenesis, the formation of hydrocarbon deposits and the chemical weathering of mudstones. Conventional methods for assessing the thermal maturity of organic matter often fail to reflect the geochemical heterogeneity between individual organic phases in mudstone samples. Here, we report an alternative, non-destructive, surficial and micro-scale (analytical spot size of ~ 300 nm with about 4 μm diffusion depth for micrometre-size organic grains) method to evaluate the thermal maturity of organic matter in mudstones using the carbon K α X-ray spectrum measured by field emission-electron probe microanalyser (FE-EPMA). Using this method, we observed correlations between parameter values derived from FE-EPMA spectra, including the peak position, the peak area and the intra-sample heterogeneity of these measurements, and independently measured vitrinite/solid bitumen reflectance for a suite of mudstones, representing different age, geological context and burial depth. With the increased values in peak area and position, we identified an increase in the carbon mass fraction of organic matter and the mean nominal oxidation state of carbon approaching zero. These trends, which are consistent with aromatisation and graphitisation, provide the rationale for using FE-EPMA to estimate the thermal maturity of organic matter. To explore some of these trends in more detail, we employed time-of-flight secondary ionisation mass spectrometry, X-ray photoelectron spectroscopy and optical reflectance measurements on a subset of samples.

Geostandards and Geoanalytical Research↗
USGS2024

An orbital comparison of a late mantling unit on Aeolis Mons with other erosion-resistant strata explored by MSL in Gale crater, Mars

Previous investigations along Curiosity 's traverse in Gale crater have explored the relationship between orbital and in situ observations. This work aims to better understand the geologic environment of units only observable from orbit and compare them to the properties of units examined by Curiosity . Here, we map an erosion-resistant dark-toned mantling unit that overlies the modern topography of Aeolis Mons (informally known as Mt. Sharp) and compare this unit to two other previously mapped dark-toned resistant units, the marker band and the mound skirting unit (MSU), that have been inferred to represent different geologic environments (lacustrine and aeolian, respectively). Visible to short wave infrared spectra from the Compact Reconnaissance Imaging Spectrometer for Mars and visual images from the High Resolution Imaging Science Experiment and Context Cameras aboard the Mars Reconnaissance Orbiter are used for this comparison. Spectral data suggest a mafic composition with minor alteration, although the composition varies more with location around Mt. Sharp rather than between units. Morphologically, the mantling unit has strong similarities to the marker band based on their consistent low-albedo, erosion-resistance, and smooth appearance, contrasting with the highly variable surface texture of the MSU. We hypothesize that all three units had a similar sediment source but experienced aqueous alteration at different times: early ubiquitous cementation in a surface aqueous environment in the mantling unit and marker band versus patchy late diagenesis in the MSU. If true, these results suggest that water activity continued within the Gale crater long after the erosion of Mt. Sharp.

JGR Planets↗
USGS2025

Pre-Acadian tectonics of the eastern Orange-Milford Belt, south-central Connecticut

This excursion presents a reinterpretation of mapping and new analytical data from the eastern Orange-Milford belt (OMB) in south-central Connecticut. The OMB is a fault-bound terrane of argillites and mafic rocks of anomalously low metamorphic grade—and of poorly constrained ages and tectonic affinity—wedged between kyanite/sillimanite-grade peri- Laurentian rocks to the west and anatectic peri-Gondwanan rocks to the east (Fig. 1A). Our data demonstrate that Ordovician(?) igneous and sedimentary rocks of the OMB were variably metamorphosed in the Ordovician and Silurian but escaped regionally pervasive, high-grade Devonian and later metamorphism. Previous interpretations (Fritts 1963a, 1965a, 1965b; Burger, 1967; Burger and others, 1968; Rodgers, 1985) described these rocks as a conformable sequence of low-grade, Ordovician to Devonian metasediments and metavolcanics. Our results reveal that the “metavolcanics” are not extrusive rocks but rather slivers of lower oceanic crust with complicated high- and low-grade metamorphic fabrics, intruded by a swarm of Silurian sheeted basalt dikes, and in fault contact with the surrounding metasediments. The purpose of this trip is to show evidence of early Paleozoic (Taconic) deformation and metamorphism preserved in rocks of the eastern OMB. These rocks remained shallow, cool, and sufficiently dry during the regionally dominant Acadian and Alleghanian orogenies to have avoided significant overprinting. As such, these rocks serve as windows into a geologic history otherwise unavailable between anatectic rocks of the peri-Gondwanan Bronson Hill, Avalon, and Gander terranes east of the Hartford basin and sillimanite-grade rocks of the peri-Laurentian Hartland and gneiss dome belts west of the OMB. We present major and trace element geochemistry including rare-earth element patterns for all mafic units in the eastern OMB as well as 40 Ar/ 39 Ar age spectra of amphibole, muscovite, and K-feldspar from rocks of the Maltby Lakes complex (of Deasy and others, 2017), Savin Schist, and Wepawaug Schist. Our evidence demonstrates that the units of the OMB have been assembled by faulting or intrusion, and that no stratigraphic relationships exist between the argillaceous schists and the metaigneous rocks.

Connecticut↗
USGS2025

Confocal laser-scanning microscopy (CLSM)-based thermal maturity of Tasmanites and progress in standardization of fluorescence microspectrometry

Evaluation of thermal maturity in vitrinite-free or vitrinite-deficient sediments via fluorescence microspectrometry can provide relevant information related to petroleum exploration and thermal history assessment. However, variation in spectral fluorescence properties of alginite macerals with increasing thermal maturity is largely underexplored. Here, authors of this study have applied confocal laser-scanning microscopy (CLSM) in conjunction with fluorescence microspectrometry to a maturity series of marine Upper Devonian Tasmanites algae from the Ohio Shale (Huron Member) and a single sample from the Marcellus Formation of the Appalachian Basin. Spectral fluorescence properties of Tasmanites were evaluated in relation to orientation, measurement location, and the number of measurements per sample, and were compared to published literature. Emission spectra of Tasmanites from continuous wave laser excitation (405 nm) were acquired from sections perpendicular and parallel to bedding and at different positions within individual Tasmanites bodies. The results showed a progressive red-shift in emission maxima (λ max ) in a large sample sized maturity series ( N = 19), e.g., 493 to 578 nm for the perpendicular section at middle position. Further, blue-shifted apex and mineral-adjacent positions within sections perpendicular to bedding were observed, with the latter being reported here for the first time. While blue-shift at apex positions can be attributed to mechanical deformation-induced reorientation of photoselected fluorophores, the blue-shifted mineral-adjacent positions could result from strain loading and development of a plastic deformation region at the mineral contact zone with Tasmanites . A decrease in standard deviation with increasing number of measured emission maxima is well-observed, and 15 to 20 individual measurements per sample appears sufficient for low standard deviation and coefficient of variance. CLSM-derived thermal maturity parameters indicated that a moderate positive correlation of red/green quotient (Q 650/500 ; R 2 = 0.67) with solid bitumen reflectance (BR o in %) exists. For reproducible results, the determination of λ max and Q 650/500 should be conducted exclusively in the middle position at perpendicular and parallel sections of the polished whole-rock pellets, where the lowest standard deviation in measurement was observed. These results strengthen the suitability and relevance of the CLSM technique in thermal maturity studies of dispersed organic matter (DOM) and contribute to the standardization of fluorescence microspectrometry methods in organic petrology investigation.

Kentucky, Ohio, Virginia, West Virginia↗
USGS2024

Constraining the earthquake recording threshold of intraslab earthquakes with turbidites in southcentral Alaska’s lakes and fjords

Strong ground motion from intraslab earthquakes, which do not produce primary paleoseismic evidence, may initiate gravity-driven turbidity flows in subaqueous basins. The resulting deposits (turbidites) can provide a paleoseismic proxy if the conditions that initiate these flows are known. To better constrain the initiating conditions, we use two recent intraslab earthquakes in southcentral Alaska, the M w 7.1 30 November 2018 Anchorage earthquake and the M w 7.1 24 January 2016 Iniskin earthquake, as calibration events. Through a multilake investigation, we document the occurrence, or the absence, of earthquake-generated turbidity flows from these two events. Both earthquakes are recorded by centimeter-scale turbidites that can be differentiated from climatically generated deposits, as well as other seismic sources based on deposit thickness, sedimentological properties, and deposit age. We show that a Modified Mercalli Intensity (MMI) of ∼V–V1/2 is the minimum shaking intensity required to generate localized sediment remobilization from deltaic slopes, and an MMI of ∼V1/2 is required to produce a deposit of sufficient thickness that a seismic origin can be confidently assigned. The documentation of seismically generated deposits in quick succession (∼2 years) with diagnostic features highlights the utility of using recent earthquakes as calibration events to investigate the subaqueous response to strong ground motion.

Alaska↗
USGS2026

Stress states on the eve of past earthquakes inform earthquake rupture through fault complexity along the San Andreas and San Jacinto faults

Estimating the evolving state of stress along active fault systems can provide insight into the conditions that generated past ground‐rupturing earthquakes and influenced their ability to propagate through areas of geometric complexity, such as fault branches and stepovers. We use quasi‐static forward numerical models that incorporate the 3D complex configuration of active faults in southern California to estimate shear tractions on the geometrically complex southern San Andreas and San Jacinto faults from 1000 to 1900 C.E. These tractions include interseismic accumulation of traction due to tectonic loading, viscoelastic relaxation of shear stress within the upper crust between earthquakes, and effects of other earthquakes on the fault network. We simulate ground‐rupturing earthquakes based on the along‐strike earthquake extents modeled by Scharer and Yule (2020) , assuming that stress drop is complete in each earthquake. We use Monte Carlo simulations to estimate uncertainty in evolving shear tractions due to uncertainties in earthquake timing and in upper‐crustal viscosity. Pre‐earthquake shear tractions typically do not exceed ∼2 MPa. Although ruptures with length <200 km have pre‐earthquake shear tractions that range from near zero to ∼1.75 MPa, these tractions are not less than ∼0.4 MPa for earthquakes with rupture length >200 km. Earthquakes with long (>200 km) ruptures occur only in the single‐stranded part of the system, whereas those with short (<125 km) rupture length and high pre‐earthquake shear traction occur near fault stepovers and branches. This suggests that high accumulated shear traction encourages longer rupture propagation, but may not be sufficient to overcome geometric complexities. This modeling approach informs our understanding of rupture propagation and provides estimates of fault shear tractions that are unavailable from direct measurements.

Callifornia↗