USGS ScienceSearch

SEARCH · USGS Science

Results for “Environmental Science and Technology - Water”

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,665 records · Page 71Linked to original sources

Unprecedented burning in tropical peatlands during the 20th century compared to the previous two millennia

Tropical peatland wildfire incidence has risen in recent decades, driven by drainage for land use and intensified by severe droughts with global climate change. These disturbances have altered vegetation structure, disrupted ecosystem functioning, and increased carbon emissions, particularly in Southeast Asia. However, the long-term history and characteristics of wildfires in tropical peatlands remain largely unknown. Here, we compiled fifty-eight macro-charcoal records from peatlands across the tropics, ranging from lowland forested to montane peatlands, to assess millennia-scale changes and controlling factors of tropical peatland burning. We divided the datasets into four main sub-regions: Neotropical, Afrotropical, Indomalayan and Australasian ecoregions to explore regional variability. Tropical peatlands had high burning levels between 0 and 850 ce , followed by a relatively low and stable period until a marked increase during the 20th century. The general trend in tropical peatland burning follows changes in global temperature, and climate variables that control the length and severity of drought events have a notable influence on peat burning before 1900 ce . During the 20th century, regional differences were observed, with declining fire trends in the Neotropical and Afrotropical regions and increasing fire trends in the Indomalayan and Australasian regions. This difference is likely attributable to human activities, and such intervention is also evident in palm swamps and hardwood swamps under similar wet, weakly seasonal climates. With the increase in anthropogenic pressures on peatlands and greater climate variability, future wildfires in peatlands are likely to become more frequent and widespread across all tropical ecoregions. Conservation and sustainable land-use practices could be used to mitigate and control peatland burning and protect these carbon-rich sinks.

Global Change Biology

Leveraging wildfire to augment forest management and amplify forest resilience

Successive catastrophic wildfire seasons in western North America have escalated the urgency around reducing fire risk to communities and ecosystems. In historically frequent-fire forests, fuel buildup as a result of fire exclusion is contributing to increased fire severity. The probability of high-severity fire can be reduced by active forest management that reduces fuels, prompting federal and state agencies to commit significant resources to increase the pace and scale of fuel reduction treatments. However, lower severity areas of wildfires also have the potential to act as “treatments,” and even catastrophic fires with large areas of high severity can still have substantial areas of lower severity fire that may be improving forest conditions locally. We quantified active management and wildfire severity across yellow pine and mixed conifer (YPMC) forests in the Sierra Nevada of California over a 22-year period (2001–2022). We did not detect increases in the area treated through time, but the area of beneficial wildfire (low to moderate severity) increased substantially, exceeding active treatment area in 8 of 22 years. Overall, beneficial wildfire treated ~17% more area than all treatments combined, and roughly four times more area than fire-related treatments alone. We then used disturbance history to evaluate resistance to high-severity wildfire and forest loss across the YPMC range. Of the 2.3 million ha YPMC of forests in 2001, 20% lost mature forests due to high-severity fire by 2022, which is nearly half of all YPMC area burned. Most of the landscape (47%) remains at risk of high-severity fire because it had no restorative disturbances, but 33% of the study area has some level of resistance to high-severity wildfire. In these areas, resistance will need to be enhanced and maintained over time via active management or managed wildfire, but these treatment needs will likely outpace capacity even under optimistic implementation scenarios. Given limited resources for implementing active management and the likelihood of a more fiery future, incorporating beneficial wildfire into landscape-level treatment planning has the potential to amplify the impact of active management treatments.

California, Nevada

Tectonic controls on volcanism and associated hydrothermal activity in a sediment-dominated mid-ocean ridge; Escanaba Trough

Mid-ocean ridges, the Earth's most extensive volcanic system, exhibit unique characteristics in sediment-dominant environments. Thick sediment cover insulates the crust and channels fluid along pathways that can lead to the formation of distinct crustal alteration patterns, exceptionally large mineral deposits, and specialized chemosynthetic ecosystems. This study presents an interdisciplinary investigation into the tectonics of the Escanaba Trough, a heavily sedimented axial valley at the southern Gorda Ridge in the Northeast Pacific Ocean. A primary challenge in such environments is overcoming the masking effect of thick sediments on basement structures that control magmatic and hydrothermal activity. We address this by employing three-dimensional (3D) magnetic modeling of high-resolution near-seafloor magnetic data collected by an autonomous underwater vehicle (AUV). The 2022 surveys with AUV Sentry provided data for 3D magnetic susceptibility models, refining our understanding of the geometry of sub-sediment laccoliths/saucer-shaped sills and hydrothermal alteration. In conjunction with a new 1:100,000 scale lithostratigraphic map, we outline the tectonic controls on the emplacement of Escanaba Trough's three main volcanic centers, characterize the geometry of its spreading segments, and provide volumetric data on the distribution of sub-sediment volcanism in the southern Gorda Ridge.

Escanaba Trough, Pacific Ocean

Mesoproterozoic to Paleozoic tectonics, Pleistocene landforms, and Holocene seismicity in the Blue Ridge: Results from integrated studies of the 9 August 2020, Mw 5.1 earthquake area near Sparta, North Carolina, USA

This field trip examines the results of integrated geologic studies of the 9 August 2020, M w 5.1 earthquake near Sparta, North Carolina, USA. The earthquake generated ~4 km of coseismic surface rupture of the Little River fault and uplifted a surface area of ~11 km 2 . The Little River fault is a thrust fault oriented 110–130°/45–70°SW, and mapped fault segments are en echelon with scarp heights from <5–30 cm. The epicenter is in polydeformed rocks of the Ashe and Alligator Back Metamorphic Suites in the eastern Blue Ridge. Bedrock structure formed during multiple Paleozoic orogenies; the regional foliation strikes NE-SW and dips SE (mean orientation 063°/52°SE). Mapping identified late Paleozoic veins and shear zones, a regional joint set striking 330–340° and 250–240°, and brittle faults that cut the Paleozoic foliation. Brittle faults oriented similar to the Little River fault are mapped up to 4 km along strike from the coseismic rupture along Bledsoe Creek valley, and the combined length of the Little River fault system is ~8 km. Paleoseismic trenches across the Little River fault corroborate the reactivation of an older fault by the 2020 earthquake and reveal two events during late Pleistocene (<50 ka). Surficial mapping identified several terrace deposits, including a deposit along Bledsoe Creek that yielded a 26 Al/ 10 Be isochron burial age of 0.46 ± 0.13 Ma and overlies a brittle fault, thus constraining the timing of movement of the fault at that location. Paleoliquefaction studies document soft-sediment deformation features in alluvium that may represent paleoseismic events. Collectively, these results highlight long-lived paleoseismicity of the Blue Ridge and that the 9 August 2020 earthquake reactivated an older, suitably oriented brittle fault in the bedrock. The Little River fault is an example of a previously unknown but active fault lying outside of known seismic zones with demonstrated recurrence of paleo-ruptures, raising questions about the assumption that damaging earthquakes are limited to areas of ongoing background seismicity, which is counter to seismic hazard assessments in the eastern United States. Bedrock mapping separates eastern Blue Ridge lithostratigraphy of the Lynchburg Group and Ashe and Alligator Back Metamorphic Suites into separate fault-bound packages juxtaposed over various 1.3–1.0 Ga basement rocks of the northern French Broad massif by the Gossan Lead fault.

North Carolina

Shallow storage of the explosive Earthquake Flat Pyroclastics magma body, Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand: Evidence from phase-equilibria experiments

Rhyolitic tuffs range widely in their crystal contents from nearly aphyric to crystal-rich, and their crystal cargoes inform concepts of upper crustal magma reservoirs. The Earthquake Flat pyroclastics (Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand) are 10 km 3 of rhyolitic tuffs with abundant (~ 40 vol.%) plagioclase and quartz, minor biotite, hornblende, and orthopyroxene, and accessory Fe-Ti oxides, apatite, and zircon, set in high-silica rhyolitic glass. Major minerals form large, euhedral phenocrysts and abundant glomerocrysts with few disequilibrium textures excepting some faintly resorbed quartz. Plagioclase phenocrysts have thick rims of nearly constant composition near An 30 , and hornblende is weakly zoned or unzoned. The abundant and texturally complex mineral assemblage contrasts with the nearby (~ 25 km), nearly synchronous, but more voluminous and crystal-moderate rhyolite tuffs from Rotoiti caldera. New H 2 O-saturated phase-equilibria results on the erupted Earthquake Flat melt (glass) determine its co-saturation with the partial phenocryst assemblage of plagioclase, quartz, biotite, and Fe-Ti oxides at: 140 MPa, 755 ºC. These closely approximate the conditions of the pre-eruptive magma body assuming it was saturated with nearly pure H 2 O and at an f O 2 of ~ Ni–NiO. Absence of hornblende and orthopyroxene from the synthesized assemblages may result from those minerals being in a peritectic reaction relation with melt to produce biotite, so they would not grow from the liquid used as starting material. Experimental results on Rotoiti rhyolite (Nicholls et al. 1992) show that the two bodies resided at similar pressures, temperatures, and f O 2 s. Lower crystal abundance of the Rotoiti tuffs may result from slight compositional differences. We interpret that the Earthquake Flat pyroclastics were sourced from the crystal-rich periphery of a mushy reservoir system with the Rotoiti occupying a more melt-rich central location. Uncertain is whether this was a single intrusion zoned continuously in crystallinity, or discrete adjacent intrusions, but our results illustrate and quantify complexities of magma storage across relatively short distances.

North Island, Taupo Volcanic Zone

Predictive regressive models of recent marsh sediment thickness improve the quantification of coastal marsh sediment budgets

Coastal marsh wetlands experience variations in vertical gains and losses through time, which have allowed them to infill relict topography and record variations in drivers. The stratigraphic unit associated with the development of the marsh also reflects the long-term importance of key ecosystem services supplied by the marsh environment, including carbon storage and storm mitigation. Mapping these coastal wetland sediments and the marsh unit thickness is challenging as traditional coastal geophysical tools are not easily deployable (acoustic methods) or are unreliable in saline-soil environments (e.g., ground-penetrating radar), leaving core-based methods the most viable mapping method. In the present study, we utilized prior information on the geologic architecture of the region to select spatial and physical metrics that likely persisted throughout evolution of the marsh during the late Holocene. We then assessed the individual and collective power of these metrics to predict marsh thickness observed from cores. Employing regressive predictive models powered by these data, we improve the quantification of marsh thickness for a coastal fringing marsh within the Grand Bay estuary in Mississippi and Alabama (USA). The information gained from this approach yields improved estimates of the carbon stocks in this environment. Additionally, the stored sediment masses reflect the past, and potential future, persistence of the Grand Bay marsh under historical and present marsh-estuarine sediment exchange fluxes. Such improvements to both the sediment budget of recent marsh stratigraphic units and the spatial extent provide new resources for comparison with large-scale landscape models, the latter of which may be used, when validated, to predict future change and ecosystem transformations.

Alabama. Mississippi

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

Extending the Boore and Abrahamson (2023) modified square-root-impedance method for the development of site amplifications consistent with the full-resonance approach to a range of VS30 values

The square-root-impedance (SRI) method is commonly used to approximate the seismic site amplifications computed using the full-resonance (FR) method for gradient shear-wave velocity ( V S ) profiles that are smoothly varying with depth. The SRI site amplifications have been observed to systematically underpredict the FR site amplifications by a ratio of FR/SRI amplifications around 1.05 to 1.3 across a wide frequency range (Boore, 2013). Recently, Boore and Abrahamson (2023; hereafter, BA23) related this difference in the SRI and FR methods to differences in the exponent η of the ratio of seismic impedances between the two methods. They proposed the implementation of a modified frequency-dependent η in the SRI method to improve its match to the FR site amplifications. This modified η was derived using only five V S profiles. We investigate the performance of the BA23 η for a wide range of realistic gradient V S profiles with V S30 ranging from 180 to 1500 m/s. These gradient V S profiles are constructed using two power-law functions of depth and are constrained by the assigned VS30 value, the depth and velocity of the half-space, and depths to shear-wave velocity horizons of 1.0 and 2.5 km/s ( Z 1.0 and Z 2.5 ) based on western United States sites. Despite observing a V S30 dependence of η, we find that the BA23 η generally works reasonably well for the range of V S profiles analyzed. Using the VS30 -dependent η derived in this study results in improvements in matching the FR site amplification compared to using the BA23 η. These improvements are more pronounced for the soft-site conditions and become modest to negligible for the stiff site conditions

Bulletin of the Seismological Society of America

Mitigating climate change by abating coal mine methane: A critical review of status and opportunities

Methane has a short atmospheric lifetime compared to carbon dioxide (CO 2 ), ∼decade versus ∼centuries, but it has a much higher global warming potential (GWP), highlighting how reducing methane emissions can slow the rate of climate change. When considering the contribution of greenhouse gas (GHG) emissions to current global warming (2010–2019) relative to the industrial revolution (1850–1900) levels, methane contributes 0.5 °C or ∼ a third of the total. The most recent post-2023 global estimates of methane emissions by bottom-up (BU) and top-down (TD) approaches for the coal mining sector are in the range of ∼41 ± 3 Tg yr −1 and 33 ± 5 Tg yr −1 , respectively. This divergence, notwithstanding overlapping confidence intervals, is a result of differences between applied TD global inversion models and BU emission inventories. Further research can help to better refine emissions from the various contributing coal mine methane (CMM) emissions sources. The coal mining sector accounts for over 10 % of global anthropogenic methane emissions. The contribution of CMM emissions to the global budget have increased since 2000, although upward and downward regional trends have been observed.

International Journal of Coal Geology

Preliminary field report of landslide hazards following Hurricane Helene

Executive Summary This report reflects our knowledge regarding the widespread landslide activity associated with Hurricane Helene observed during the U.S. Geological Survey’s (USGS) mission assignment to North Carolina in October 2024. The material in this report was originally prepared for the Federal Emergency Management Agency under mission assignment DR-4827-NC. The data and commentary in this report are reflective of a report provided to the Federal Emergency Management Agency (FEMA) on October 18, 2024, as well as information provided in briefings at the Buncombe County Emergency Operations Center. The report has been modified for public dissemination. This assessment was based on systematic visual examination and mapping of landslide locations from aerial and satellite imagery, visual and photographic observations from low-level helicopter overflights and conversations with local landslide experts from the North Carolina Geological Survey and Appalachian Landslide Consultants PLLC, and more than 50 years of combined landslide hazard professional experience of the mission-assigned field team. No systematic field investigations were done by the USGS. While responding to the event, the USGS did not identify any landslides that posed an immediate major threat to recovery personnel in parts of nine counties in North Carolina (Avery, Buncombe, Henderson, McDowell, Mitchell, Polk, Rutherford, Watauga, and Yancey); however, threats from renewed landslide activity may remain heightened in localized areas for months or even years. Known areas of the most abundant landslide occurrence include Bat Cave, Lake Lure, Chimney Rock, Swannanoa, Black Mountain, Fairview, steep areas in Asheville, and the Blue Ridge Parkway. The USGS shared detailed locations of known landslides with the Emergency Operations Centers. The thousands of landslide scars on hillsides and landslide deposits on flatter ground may present some threat to recovery activities. Soil and rocks will continue to erode from newly exposed landslide scars and may pose a threat to people and infrastructure who are immediately nearby. In general, the steeper and taller the landslide scar, the greater the potential threat. This threat is heightened during periods of rainfall and increases with the duration and intensity of rainstorms. Very heavy rainfall, or repeated rainfall events during short periods, could also initiate new landslides on steep slopes. Excavation of landslide deposits, particularly excavation of those deposits directly adjacent to steep slopes, may also pose a threat to nearby people and equipment. An interagency collaborative mapping effort led by the USGS that informed this assessment identified 1,155 landslide locations by the October 2024 briefings, but that number increased to 2,217 in a final reviewed version of the locations published in January 2025. Locations were mapped from satellite imagery, fixed-wing and helicopter surveys, media and social media, and field reports in the 3 weeks following the passage of the remnants of Hurricane Helene. USGS products outlined in this report are publicly available and include geotagged photographs from aerial reconnaissance, hazard models, an interactive view of mapped landslide locations, and landslide safety and education resources.

North Carolina, South Carolina, Tennessee, Virgini

Mafic alkaline magmatism and rare earth element mineralization in the Mojave Desert, California: The Bobcat Hills connection to Mountain Pass

Occurrences of alkaline and carbonatite rocks with high concentrations of rare earth elements (REE) are a defining feature of Precambrian geology in the Mojave Desert of southeastern California. The most economically important occurrence is the carbonatite stock at Mountain Pass, which constitutes the largest REE deposit in the United States. A central scientific goal is to understand the genesis of the carbonatite ore body in the context of widespread REE-rich igneous activity. A swarm of mafic alkaline (shonkinite) dikes has been mapped and sampled at Bobcat Hills, 65 km southeast of the Mountain Pass mine. Whole-rock geochemistry and zircon geochronology demonstrate a clear affinity to the ca. 1.4 Ga Mountain Pass intrusive system. Bobcat Hills dikes have comparably high REE concentrations (La ∼1,000× chondritic) and an error-weighted mean 207 Pb/ 206 Pb zircon crystallization age of 1,426 ± 2 Ma (2 σ ). Unlike the alkaline intrusions at Mountain Pass, which have abundant inherited zircon from Paleoproterozoic basement rocks and crustally influenced oxygen isotope compositions (δ 18 O zircon = 6.5–7.5‰), the Bobcat Hills dikes lack any evidence of crustal assimilation and have oxygen isotope values that overlap a mantle range (Bobcat Hills average δ 18 O zircon = 5.6 ± 0.3‰). The dikes were a high-temperature, early center of mafic alkaline magmatism in the Mojave Desert that serve as a snapshot of melt generation from a spatially extensive, metasomatized mantle source. We propose that modification of the crust over many tens of Myr at Mountain Pass created an environment that favored crustal assimilation and enabled ascent of late-stage, REE-rich carbonatite magmas.

California

A methodology to estimate CO2 and energy gas storage resources in depleted conventional gas reservoirs

Depleted hydrocarbon reservoirs are subsurface geological structures capable of sequestering vast quantities of carbon dioxide (CO 2 ) as well as storing other energy gases for later usage, such as natural gas, and potentially hydrogen (H 2 ). Here we outline a methodology to quantify multi-gas storage resources in depleted conventional gas reservoirs for usage in assessments by the United States Geological Survey (USGS) at the scale of sedimentary basins. The methodology consists first of quantifying accessible pore volume in a depleted reservoir for natural gas storage using up to three equations. Input data are derived from commonly reported or estimated reservoir parameters and natural gas production volumes, and equations may be combined in linear models to improve pore volume estimates. Storage estimates from these equations are tested and validated for 31 reservoirs in the Michigan Basin Province, USA that were previously converted to underground gas storage facilities and have known (federally reported) natural gas storage capacities. Secondly, natural gas storage capacities can be transformed via fluid substitution calculations to estimate the storage resources for non-native fluids, applied here for, CO 2 , H 2 , and methane-H 2 blends, accounting for molecule-specific deviations from ideal gas behavior at reservoir pressures and temperatures as well as differing storage efficiencies. Importantly, the storage of non-native fluids may not be appropriate in all depleted gas reservoir settings due to potential risks like leakage, in particular in the case of H 2 storage, requiring additional knowledge of caprock sealing capacity. Given this caveat, we demonstrate the fluid substitution method for natural gas reservoirs of the Northern Niagaran Reef and Southern Niagaran Reef USGS plays in the Michigan Basin Province, as these trends of Silurian pinnacle reefs are capped with tight-sealing evaporite facies. The deterministic equations outlined from this methodology can be incorporated into future probabilistic USGS gas storage assessments for CO 2 , H 2 , and natural gas resources in the United States.

Michigan

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

Summary of results from electromagnetic and galvanic soundings on Kilauea Iki lava lake, Hawaii

The purpose of this report is to summarize various electrical sounding studies made on the Kilauea Iki lava lake and to present some of the preliminary interpretations that resulted from the measurements. The following discussion is not intended as a complete interpretation, but is intended to summarize the major points made in an invited oral presentation at the American Geophysical Union fall meeting held in San Francisco, Calif. on December 6, 1976. This paper was part of a group of papers dealing with various geophysical studies recently made on the Kilauea Iki lava lake. A general plan map of Kilauea Iki lava lake is shown in figure 1. Electromagnetic (EM) loop-loop soundings, D.C. resistivity soundings, and VLF apparent resistivity measurements were made on the lave lake surface. All of these surveys were made to determine the geoelectric structure of the lava lake and to determine parameters of the molten basaltic lens that exist in the partially solidified lava lake.

Hawaii

The 2023 Alaska National Seismic Hazard Model

US Geological Survey (USGS) National Seismic Hazard Models (NSHMs) are used extensively for seismic design regulations in the United States and earthquake scenario development, as well as risk assessment and mitigation for both buildings and infrastructure. This 2023 update of the long-term, time-independent Alaska NSHM includes substantial changes to both the earthquake rupture forecast (ERF) and ground motion models (GMMs). The ERF includes numerous additions to the finite-fault model, considers two deformation models, and introduces updated declustering and smoothing algorithms in the gridded background seismicity model. For the Alaska–Aleutian subduction zone, megathrust earthquakes occur on an updated structural and segmentation model, and the moment magnitude (M) 8+ rupture and rate model include a logic tree branch that considers slip rates derived from geodetic models of interface coupling. The megathrust model considers multiple models of down-dip width, and magnitudes are computed using newly developed scaling relations. For subduction intraslab events and subduction interface events with M < 7, the 2023 update uses a smoothed seismicity model with rupture depths derived from Slab2. The 2023 model updates GMMs in all tectonic settings using the recently published Next Generation Attenuation Subduction (NGA-Sub) GMMs for subduction interface and intraslab events, and the NGA-West2 GMMs for active crustal settings. Collectively, additions and updates to the Alaska NSHM result in hazard increases across most of south-central Alaska relative to the previous model, published in 2007. These changes are primarily due to the adoption of updated rate models for the large-magnitude interface events and the NGA-Sub GMMs that have much higher aleatory variability (sigma), consistent with global observations, and that include models of epistemic uncertainty.

Alaska

Relationship between peak and cumulative ground motions from 49 Mw 3-6 earthquakes in the San Francisco Bay Area

We examine the relationship between peak ground velocity (PGV) and cumulative absolute displacement (CAD) for 49 M w 3 – 6 earthquakes in the San Francisco Bay Area (SFBA) and gain insight into the spatiotemporal partitioning of seismic energy in ground motion records with respect to source, path, and site effects. PGV and CAD are positively correlated, but there can be large deviations from the average trend. For example, ground motion records with either very long duration resonance or short pulse-like motions will have higher or lower CAD, respectively, but could have very similar PGV. We perform principal component analysis (PCA) on PGV-CAD for >7000 records in the SFBA with the goal of investigating what influences positive or negative anomalies in cumulative motions. PCA rotates the PGV-CAD datapoints into two principal components, where the one with the larger variance, which we call the “primary intensity” represents mostly the distance-dependence of ground motion amplitudes. The other principal component, which we call the “excess motion”, represents the deviation from cumulative motions that would be typical for a ground motion record with a given PGV. The excess motion will be positive in the case of records with long duration ringing and will be negative for short duration pulse-like ground motions. We find that excess motion is generally positive at sites in sedimentary basins and in soft sediments around the SF Bay. Excess motion is generally negative in the very near field, as well as at sites on hard bedrock. We discuss the findings here in terms of implications for seismic hazard applications and other wave propagation phenomena.

California

Eruptive history within the vicinity of Al Madīnah in northern Harrat Rahat, Kingdom of Saudi Arabia

The northernmost part of the Harrat Rahat volcanic field contains early Pleistocene to Holocene mafic eruptive products within the vicinity of the city of Al Madīnah, Kingdom of Saudi Arabia. A detailed geologic investigation into the eruptive history of a 570 square kilometer (km 2 ) area covering Al Madīnah and the surrounding area has yielded 33 mapped Quaternary volcanic units consisting of lava flows, scoria cones, and shield volcanoes. These eruptive products consist of continental, intraplate alkalic and minor transitional basalts, hawaiites, and a single mugearite that were emplaced from at least 1,014±14 thousand years ago (ka) to a single Holocene eruption in 1256 C.E. Lava flows are generally 10 to 15 kilometers (km) long (but can reach 23 km long), 1 to 3 km wide, and at least 10 meters thick. Most of the mapped units erupted episodically between 400 and 340 ka and 180 and 100 ka. Despite small individual volumes (less than 1 cubic kilometers dense rock equivalent), each unit represents eruption of a distinct magma batch that was strongly influenced by clinopyroxene, olivine, and plagioclase fractionation. Some of these units are interpreted to have undergone magma mixing pre- and (or) syneruptively. Integrating eruption ages, geochemistry, and paleomagnetic data yields evidence that some eruptions were temporally and (or) spatially clustered. Aligned scoria cones and elongate vent edifices were constructed atop fissure vent systems that reflect the local stress field, which controls dike ascent through the middle and upper crust.

Professional Paper

HyFlood: A surrogate-model-based framework for compound coastal flooding

Compound coastal flooding is a major threat to low-lying coastal regions and is expected to intensify under future climate change projections. However, modeling the joint interaction of waves, storm surge, tides, and rainfall remains computationally demanding, limiting the development of fast and reliable forecast tools. Here we present HyFlood, a hybrid statistical-numerical downscaling framework capable of computing and mapping high-resolution compound flood hazards while substantially reducing the computational cost compared with fully process-based hydrodynamic modeling. HyFlood combines statistical sampling and selection algorithms with a cascade of reduced-complexity surrogate models that emulate nearshore wave transformation, surf-zone hydrodynamics, and coastal, fluvial, and pluvial flooding. The surrogate models employ machine-learning and regression algorithms applied to a low-dimensional representation of the flooding outputs, obtained through statistical dimensionality reduction. The framework is demonstrated in southern O'ahu, Hawai'i, a region exposed to elevated sea levels driven by tides, waves, and storm surge along with frequent precipitation-driven flash flooding. Validation of the surrogates against the physics-based model outputs demonstrates that HyFlood accurately reproduces daily maxima of spatially distributed flooding depths. This hybrid approach offers a scalable and efficient tool to better quantify how changes in flooding drivers translate into hazard and impact assessments, and to support compound-flood risk assessments and climate-change adaptation planning.

Hawaii