USGS ScienceSearch

SEARCH · USGS Science

Results for “Frontiers in Enviornmental Science”

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

Seismic velocity changes from repetitive seismicity at Mauna Loa prior to and during its 2022 eruption

Mauna Loa’s short-lived eruption from late November to early December 2022 marked the culmination of nearly a decade of elevated seismic activity and geodetic inflation. The volcano has been monitored by a network of permanent, short period and broadband seismometers. I used the continuous waveform data from that network starting in 2012 to generate a catalog of seismicity that enhances the US Geological Survey Hawaiian Volcano Observatory’s public seismic catalog with four times the number of earthquakes, which were then grouped by waveform similarity. Analysis of subtle delays in the timing of arrivals of scattered waves between pairs of earthquakes in this catalog yields a history of small changes in the shallow seismic velocity structure of the volcano. Seismic velocities have been shown at other volcanoes to change during unrest and eruption. My results show a decrease in seismic velocity centered on the summit beginning in September 2022, corresponding to the onset of a vigorous precursory swarm of seismic activity and shallow inflation. During the eruption itself, I observe large changes due likely to dike opening along the northeast rift zone and deflation of the summit reservoir. However, seismic velocity changes associated with non-volcanic sources such as ground shaking from large earthquakes and meteorological influences at seasonal and diurnal time scales are also observed, and these dominate the velocity changes prior to the eruption. Proper accounting of these effects will be a requirement for use in real-time monitoring, and this work serves as a starting point in that endeavor for Mauna Loa.

Hawaii

Improved prediction of postfire debris flows through rainfall anomaly maps

Predicting where runoff-generated debris flows might occur during rainfall on steep, recently burned terrain is challenging. Studies of mass-movement processes in unburned areas indicate that event locations are well-predicted by rainfall anomaly, R* , in which peak observed rainfall is normalized by local rainfall climatology. Here, we use remote and field methods to map debris flows triggered within the 2020 Dolan Fire burn area in coastal California, demonstrate that a short-duration R* metric predicts debris-flow occurrence more effectively than absolute peak intensity or longer-duration rainfall metrics, and show that incorporating an R* criterion into an existing debris-flow likelihood model can reduce false positive predictions and improve accuracy. We test R * at three other climatically distinct fires in California, demonstrating its utility for mapping likely debris-flow locations in different climates. We also consider how R* can benefit postfire debris-flow prediction given recent increases in climatological variability within individual burn perimeters.

Callifornia

The 3D Elevation Program—Supporting the economy of Massachusetts

Introduction Massachusetts extends from the mountains of the Appalachian system in the west of the State to the sandy beaches and rocky shorelines of the Atlantic coast in the east. Inland topographic data support a wide range of important activities, including geologic mapping, transportation planning, forest and wildlife management, quantifying ecological services, water supply protection, commonwealth-wide infrastructure planning, local site planning, and flood-plain management. Nearshore bathymetry can be used to support coastal portions of the Commonwealth by addressing the combined threats of ocean warming, strong storm surge, and rising sea levels. The maintenance and (or) expansion of Massachusetts ports (for instance, Boston, New Bedford) and Cape Cod sediment management depends upon the accurate mapping of bathymetry and the frequent influx of sediment and redeposition. Critical applications that address the broad range of requirements depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features. The 3D Elevation Program (3DEP) is managed by the U.S. Geological Survey (USGS) in partnership with Federal, State, Tribal, U.S. territorial, and local agencies to acquire consistent lidar coverage at quality level 2 or better to meet the many needs of the Nation and Massachusetts. The status of available and in-progress 3DEP baseline lidar data in Massachusetts is shown in figure 1. 3DEP baseline lidar data include quality level 2 or better, 1-meter or better digital elevation models, and lidar point clouds, and must meet the Lidar Base Specification version 1.2 ( https://www.usgs.gov/3dep/lidarspec ) or newer requirements. The National Enhanced Elevation Assessment identified user requirements and conservatively estimated that availability of lidar data would result in at least $1.23 million in new benefits annually to Massachusetts. The top 10 Massachusetts business uses for 3D elevation data, which are based on the estimated annual conservative benefits of 3DEP, are shown in table 2.

Massachusetts

Haunted Summerville: Ghostly lights or earthquake lights?

Among the colorful local lore in the Charleston, South Carolina, area, are a number of ghost stories, shared not only over campfires but also in published books. Among the most well-known of the stories is the tale of the Summerville Light. Local lore holds that a strange light sometimes seen in a remote area is a lantern carried by the ghost of a local woman who once waited hours for her husband, who turned out to have been decapitated earlier that day in a train accident (DePoppe, 2023). Extant sources suggest the ghost stories began to circulate in the 1950s to 1960s. So pervasive was the lore that (Old) Sheep Island Road became known among local residents as Light Road, with a local stretch of road known today as Old Light Road. Reviewing the location where the lights appear as well as the nature of accounts, I suggest that many if not all of the anecdotal observations can be most readily attributed to natural phenomena, including earthquake lights from earthquakes that were too small to be felt. Accounts of lights near Summerville cluster in proximity to the generally accepted epicenter of the 1886 Charleston, South Carolina, earthquake, where foreshocks to the 1886 mainshock were apparently concentrated, and within a few kilometers of three M3.5 – 4.4 earthquakes in 1959 and 1960.

South Carolina

Flood-inundation maps of the Blue River and major tributaries in and near Kansas City, Missouri, 2023–25

Digital flood-inundation maps for 35.5 miles of the Blue River, in and near Kansas City, Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the City of Kansas City. Maps were also developed for a combined 7.3 miles of Brush Creek (4.4 miles) and Indian Creek (2.9 miles), two primary tributaries to the Blue River in the study reach. The flood-inundation maps, available through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding, depth of water above first-floor elevations in selected buildings, water velocity, water-surface elevations, and time of travel corresponding to selected water levels or gage heights (hereafter referred to as “stages”) at nine reference USGS streamgages—five on the Blue River (USGS station numbers 06893590, 06893578, 06893500, 06893150, and 06893100), two on Brush Creek (USGS station numbers 06893562 and 06893557), one on Indian Creek (USGS station number 06893390), and one on the Missouri River (USGS station number 06893000). Near-real-time stages at these streamgages may be obtained from USGS Water Data for the Nation at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at seven of these sites (USGS station numbers 06893000, 06893590, 06893578, 06893500, 06893150, 06893557, and 06893390). Flood profiles were computed for eight map reaches by means of two-dimensional hydraulic models. The models were calibrated using the USGS-developed stage–streamflow relations at each reference streamgage. Two map reaches on the lower Blue River also include the effects of backwater from the Missouri River on flood inundation extent and hydraulic characteristics. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “Action Stage” or near bankfull to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root-mean-square error) to delineate the area flooded at each water level and the associated hydraulic characteristics. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that could be critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Kansas, Missouri

Central Valley Hydrologic Model version 2 (CVHM2): Decision support tool for groundwater and land subsidence management

The San Joaquin Valley (SJV) of California is one of the world’s most productive agricultural regions. Reliance on groundwater has led to some of the greatest rates of human-induced land subsidence in the world in the 20th century, as well as more recently. The United States Geological Survey (USGS) has recently developed an integrated surface–subsurface hydrologic model, the Central Valley Hydrologic Model 2 (CVHM2), that represents the major components of the hydrologic system of California’s Central Valley. In this study, CVHM2 was applied as a decision support tool while simulating various management strategies to mitigate the land subsidence caused by the extraction of groundwater. CVHM2 was extended through to 2073 and applied to simulate management scenarios in terms of three primary drivers and their impact on subsidence along the Delta–Mendota Canal (DMC), a critical piece of infrastructure in the western SJV. The drivers considered were agricultural water demands, managed aquifer recharge (MAR), and changes in future climate. The results show that future subsidence is most sensitive to water demands, second most sensitive to future changes in climate, and relatively insensitive to MAR when it is applied as a surface application in the western SJV. However, we demonstrate via proof-of-concept scenarios that the MAR is capable of arresting subsidence when implemented via injection below the Corcoran Clay Member of the Tulare Formation instead of as a surface application. We also examine the uncertainty that is the result of climate variability and how to use the tool to identify the most appropriate strategies to constrain future subsidence to acceptable levels.

California

Seismic techniques and suggested instrumentation to monitor volcanoes

Introduction Changes in the pressure or location of magma can stress or break surrounding rocks and trigger flow of nearby waters and gases, causing seismic signals, such as discrete earthquakes and tremor. These phenomena are types of seismic unrest that commonly precede eruption and can be used to forecast volcanic activity. Mass movements at the surface, including avalanches, debris flows, and lahars, may also generate seismic signals that are specifically addressed in chapter H , this volume (Thelen and others, 2024). Our focus in this chapter is to determine the levels of instrumentation recommended to produce high-quality, well-constrained seismic observations important for early warning of impending eruptions, detecting changes in ongoing eruptions, and characterizing other hazardous volcanic events. There are emerging techniques and new types of instrumentation, such as distributed acoustic sensing or rotational seismometers, that we do not consider here. These types of instrumentation show promise for monitoring but still require maturation before being considered more generally in volcano monitoring. Most of the capabilities mentioned below are universal for all types of volcanic systems, although some are best applied to stratovolcanoes with an apical single vent. In some settings, such as calderas or shield volcanoes, we must broaden coverage to include multiple possible storage regions or vent locations. As an example, Thelen (2014) discretized the long rift zones of shield volcanoes in Hawaiʻi as a set of evenly spaced “vents.” In this construct, each vent comes with recommendations, and several thousand network configurations were simulated to assess the effect on network quality levels and to determine the most efficient network design. The same process could be applied in a caldera setting or a volcanic field, where an evenly spaced grid of potential vents is considered. Localized recommendations for each unique system are beyond the scope of this report and left up to local experts to assess based on the conditions, restrictions, and requirements of each volcano.

Scientific Investigations Report

Ordovician stratigraphy, structure, and karst of the Falling Spring Valley, Alleghany County, Virginia, USA

This one-day trip highlights new findings on a preliminary bedrock geologic map that shows results from ongoing geologic mapping in the Falling Spring Valley of Alleghany County, Virginia, USA, which is the southern end of the larger Warm Springs Valley, an elongated anticlinal valley rimmed by Ordovician and Silurian siliciclastic rocks, and which is famous for its thermal springs. This mapping includes stratigraphic, structural, and karst field and lab research focused on the Ordovician strata exposed in the area, the oldest of which is the dolomitic upper part of the Beekmantown Formation (Lower Ordovician, Darriwilian), and the youngest of which is the Juniata Formation (Upper Ordovician, Katian), a sequence of siliciclastic redbeds. Warm Springs Valley is the location of the only known caves in the eastern United States—three at present—with thermal waters flowing in some of their passages. Stops on the trip will highlight key details from mapping efforts, primarily within the structurally deformed Ordovician carbonate sequence that is exposed in the core and limbs of the anticline, as well as the associated karst features that are developed in those carbonate rocks, including results of recent dye traces and water temperature monitoring that have improved our understanding of the karst hydrogeologic systems developed in these strata.

Virginia

Channel morphology and large wood control postfire debris-flow erosion and deposition

Runoff-generated debris flows are a known response to wildfire, and accurately predicting the volume of these debris flows is important for estimating the magnitude of downstream hazards. Prior data collection efforts have focused on debris-flow volume measurements at catchment outlets, but few studies have considered how erosion and deposition modulate the volume of sediment arriving at catchment outlets. This study takes advantage of a high-resolution dataset to document the factors that control the total debris-flow volume reaching the catchment outlet during a fatal postfire debris flow. Using pre- and post-event airborne lidar, satellite imagery and field mapping, we found that a postfire debris flow in the Black Hollow catchment in northern Colorado eroded 136,000 ± 30,000 m 3 and redeposited 27,000 ± 7,500 m 3 in the main channel. Most of the in-channel deposition (52% by volume) occurred where a confined channel reach transitioned to an unconfined channel reach downstream, allowing the flow to widen and deposit material. Wood jams played multiple roles in the debris-flow dynamics, both nucleating deposition (25% of the deposit volume was stored behind wood jams) and exacerbating erosion (50% of the total erosion occurred downstream from a wood dam break). The remaining deposition occurred due to spatial changes in channel slope as well as deposition observed at newly formed channel bars. Using these data in this study, we identified topographic and vegetation metrics that can be used (pre-event) to anticipate where deposition may occur in channels prior to a debris flow.

Colorado

Conceptual and numerical groundwater flow model of the Iowa River alluvial aquifer near Tama County, Iowa, 1980 through 2022

The Iowa River alluvial aquifer is an important source of water on the Meskwaki Settlement in Tama County, Iowa, which is land owned by the Sac & Fox Tribe of the Mississippi in Iowa (commonly known as the Meskwaki Nation). The U.S. Geological Survey constructed a groundwater flow model, including a conceptual and numerical model, of the Iowa River alluvial aquifer and underlying hydrogeologic units near the Meskwaki Settlement in Tama County, Iowa, for the period of January 1980–August 2022 to estimate the fraction of water pumped from the Iowa River alluvial aquifer by Meskwaki Settlement wells that is derived from streamflow depletion in the Iowa River and its tributaries. Streamflow depletion is a reduction in streamflow caused by groundwater pumping and includes the interception by groundwater production wells of water that otherwise would have been discharged to streams (called “captured groundwater discharge”) and induced infiltration of streamflow to the production wells. Calibrated model runs were performed with no simulated pumping and simulated pumping only at Meskwaki Settlement wells, and the change in simulated flow rates between the groundwater system and streams for the two model runs represents the amount of streamflow depletion in the Iowa River and tributary streams resulting from pumping at the Meskwaki Settlement wells. Streamflow depletion in the Iowa River and its tributaries as a percentage of simulated pumping at the Meskwaki Settlement wells was calculated by dividing this difference by the total simulated pumping rate for the Meskwaki Settlement wells. The model results demonstrate that the mean monthly streamflow depletion, including induced infiltration and captured discharge, in the Iowa River and its tributary streams as a percentage of mean monthly pumping at the Meskwaki Settlement wells was 97.4 percent and ranged from 65.4 to 112 percent. Of the total streamflow depletion, mean monthly induced recharge was 20.9 percent and ranged from 4.9 to 37.2 percent. Mean monthly captured discharge was 76.5 percent and ranged from 57.1 to 97.1 percent. These results indicate that most of the water pumped from the Meskwaki Settlement wells is the result of streamflow depletion, in the form of both induced infiltration and captured discharge.

Iowa

Lunar analog study using portable gamma-ray neutron detector: Radiochemical mapping of silicic-to-basaltic volcanic terrains in the San Francisco Volcanic Field, Arizona

Studying lunar silicic volcanism provides key insights into the Moon’s crustal evolution, magmatic processes, and volcanic history. The Lunar-VISE (Lunar Vulkan Imaging and Spectroscopy Explorer) mission will investigate the Gruithuisen domes, a unique lunar region hypothesized to have formed through silicic volcanism. Using instruments on a Firefly Aerospace lander and a Honeybee Robotics rover, Lunar-VISE will analyze mineralogy, geochemistry, and surface properties to determine the origin and evolution of the domes, with a gamma-ray and neutron spectrometer (LV-GRNS) among its payload instruments. In preparation for this mission, we conducted preliminary fieldwork using a handheld gamma-ray neutron detector with NaI(Tl) and Cs LiYCl :Ce (CLYC) scintillators. This study focuses on various rhyolitic and basaltic volcanic centers in the San Francisco Volcanic Field (SFVF) near Flagstaff, Arizona—specifically Sugarloaf Peak, Bonito Lava Flow, and Robinson Mountain, a region containing several well-characterized lunar analog sites. The SFVF was selected for this study due to its broad compositional diversity spanning basaltic to rhyolitic compositions, its well-preserved volcanic morphologies, and its extensive use in previous NASA field campaigns and astronaut training exercises, making it an ideal terrestrial laboratory for testing planetary exploration techniques. We measured natural radioactivity, specifically from potassium (K), thorium (Th), and uranium (U) and other elements in their decay chains, which serve as diagnostic tracers of magmatic differentiation and crustal evolution processes, to assess geochemical variability across compositionally diverse terrains. The detector’s sensitivity was assessed across varying concentration levels. Regionally averaged concentrations of radioisotopes were determined by gamma-ray spectroscopy at selected sites. Our measurements reveal significant compositional variations between sites, with Sugarloaf Peak (rhyolitic, silicic dome) exhibiting the highest average radioisotope concentrations (K: 3.29 ± 0.24 wt%, U: 14.35 ± 1.81 ppm, Th: 27.14 ± 2.43 ppm), while Bonito Lava Flow (K: 1.13 ± 0.08 wt%, U: 3.86 ± 0.52 ppm, Th: 7.52 ± 1.02 ppm), and Robinson Mountain (K: 1.46 ± 0.16 wt%, U: 5.59 ± 1.06 ppm, Th: 11.75 ± 1.98 ppm) show lower concentrations aligned with basaltic compositions. Gamma-ray fluxes were elevated by a factor of approximately three to four at Sugarloaf Peak relative to nearby basaltic terrains, consistent with expected geochemical differentiation patterns. Furthermore, at meter scales, proximity to geological features significantly affects measurements. Th concentrations adjacent to a cliff face at Sugarloaf’s base were 28% higher than values measured 3-5 meters away from the same feature, demonstrating localized compositional heterogeneity. Analysis of station-by-station measurements reveals that K, U, and Th concentrations show an elevation trend at Robinson Mountain and generally higher concentrations near Sugarloaf Peak’s summit, suggesting progressive magmatic differentiation and/or the presence of more evolved lithologies at higher elevations. By mapping these elements on Earth using GRNS, we aim to optimize measurement techniques for Lunar-VISE and similar rover-borne missions by demonstrating the applicability of portable GRNS instruments for planetary surface exploration in an analog environment.

Arizona

Effects of stochastically-simulated near-fault ground motions on soil liquefaction

The scarcity of historically recorded near-fault ground motions poses a challenge to systematically understanding the influence of near-fault effects on various types of seismic demands for engineering purposes. In particular, the current state of knowledge of the influence of ground-shaking intensity on soil liquefaction and its consequences does not specifically account for the effects of near-fault ground motion characteristics. In this study, the influence of near-fault ground motions on liquefaction triggering and lateral spreading are investigated using non-linear modeling of a hypothetical liquefiable soil column in the finite-element computational platform OpenSees subjected to simulated ground motion time series that represent strong earthquake shaking in the near field. The simulated ground motion time series and resulting datasets are based on a parametric stochastic model and are developed for a range of source and path parameters to represent a realistic variability of ground motion characteristics. Dependencies between ground motion intensity measures (IMs) and liquefaction demand parameters are investigated for near-fault pulse and nonpulse-like ground motion sets. Evolutionary IMs, such as cumulative absolute velocity (CAV) and the time-varying magnitude-adjusted peak ground acceleration (PGAM), are considered in developing liquefaction triggering probability density functions. Post-liquefaction triggering responses such as lateral spreading displacements are examined in relation to PGA M and CAV. The ground motion simulations are validated by comparing their liquefaction-capacity PGA M fragilities and post-triggering CAV vulnerability relationships to historical records from the 1994 Northridge earthquake in California, USA. Finally, a path forward for future studies that includes finding systematic differences in the IM-liquefaction demand relationships between near-fault and far-field stochastic ground motion sets is outlined.

Conference Paper

Rainfall thresholds for postfire debris-flow initiation vary with short-duration rainfall climatology

The size, frequency, and geographic scope of severe wildfires are expanding across the globe, including in the Western United States. Recently burned steeplands have an increased likelihood of debris flows, which pose hazards to downstream communities. The conditions for postfire debris-flow initiation are commonly expressed as rainfall intensity-duration thresholds, which can be estimated given sufficient observational history. However, the spread of wildfire across diverse climates poses a challenge for accurate threshold prediction in areas with limited observations. Studies of mass-movement processes in unburned areas indicate that thresholds vary with local climate, such that higher rainfall rates are required for initiation in climates characterized by frequent intense rainfall. Here, we use three independent methods to test whether initiation of postfire runoff-generated debris flows across the Western United States varies similarly with climate. Through the compilation of observed thresholds at various fires, analysis of the spatial density of observed debris flows, and quantification of feature importance at different spatial scales, we show that postfire debris-flow initiation thresholds vary systematically with short-duration rainfall-intensity climatology. The predictive power of climatological data sets that are readily available before a fire occurs offers a much-needed tool for hazard management in regions that are facing increased wildfire activity, have sparse observational history, and/or have limited resources for field-based hazard assessment. Furthermore, if the observed variation in thresholds reflects long-term adjustment of the landscape to local climate, rapid shifts in rainfall intensity related to climate change will likely induce spatially variable shifts in postfire debris-flow likelihood.

Arizona, California, Colorado, Nevada, New Mexico,

Niobium-rich minerals from the Sheep Creek carbonatite-related deposits, Montana, USA

The Sheep Creek area in southwest Montana contains a number of small, tabular, carbonatite-related deposits that are rich in Nb- and REE-minerals. Two of the “vein-dikes” were mined at a small scale for Nb in the late 1950s. The most abundant Nb-bearing minerals are columbite-(Fe), nioboaeschynite-(Ce), baotite, and two varieties of Ti-Nb-Fe-oxide. An early oxide with Ti/Nb (atomic) near unity later broke down to a mixture of columbite + TiO2 with distinctly lower Nb content. Zircon intergrown with the Nb-REE minerals was dated by U-Pb methods to 116±2 Ma. This age is in agreement with previous dating (Th-Pb on monazite, Re-Os on molybdenite), and indicates that the Sheep Creek deposits are early Cretaceous in age.

Montana

Intense alteration on early Mars revealed by high-aluminum rocks at Jezero Crater

The NASA Perseverance rover discovered light-toned float rocks scattered across the surface of Jezero crater that are particularly rich in alumina ( ~ 35 wt% Al 2 O 3 ) and depleted in other major elements (except silica). These unique float rocks have heterogeneous mineralogy ranging from kaolinite/halloysite-bearing in hydrated samples, to spinel-bearing in dehydrated samples also containing a dehydrated Al-rich phase. Here we describe SuperCam and Mastcam-Z observations of the float rocks, including the first in situ identification of kaolinite or halloysite on another planet, and dehydrated phases including spinel and apparent partially dehydroxylated kaolinite. The presence of spinel in these samples is likely detrital in origin, surviving kaolinitization, pointing to an ultramafic origin. However, the association of low hydration with increased Al 2 O 3 abundances suggests heating-induced dehydration which could have occurred during the lithification or impact excavation of these rocks. Given the orbital context of kaolinite-bearing megabreccia in the Jezero crater rim, we propose an origin for these rocks involving intense aqueous alteration of the parent material, followed by dehydration/lithification potentially through impact processes, and dispersion into Jezero crater through flood or impact-related processes.

Communications Earth & Environment

Young explosive eruptions from the Clear Lake volcanic field

The Clear Lake volcanic field is the northernmost and youngest field in a chain of volcanic fields in and near the California Coast Range mountains. For 2 million years, numerous eruptions have happened around (and through) Clear Lake. The most recent period of activity in the Clear Lake volcanic field probably started around 40,000 years ago and was mainly explosive eruptions concentrated on and near faults in and around the lake. The combination of hot magma and groundwater created explosive releases of hot steam. The craters, called maars, threw out volcanic ash, pumice, and lava fragments, which fell back to the ground and draped the surrounding landscape. The deposits are visible in many places around the southeast end of the lake.

California

Comparing subduction ground-motion models to observations for Cascadia

We evaluate Cascadia subduction ground-motion models (GMMs), considered for the 2023 US National Seismic Hazard Model (NSHM) update, by comparing observations to model predictions. The observations comprise regional recordings from intraslab earthquakes, including contributions from 2021 and 2022 events in southern Cascadia and global records from interface earthquakes. Since the 2018 NSHM update, new GMMs for Cascadia have been published by the Next Generation Attenuation (NGA)-Subduction Project that require independent evaluation. In the regional intraslab comparisons, we highlight a characteristic frequency dependence for Cascadia data, with short periods having lower ground motions and longer periods being comparable to other subduction zones. We evaluate differences in northern and southern Cascadia and find that the NGA-Subduction GMMs developed using southern Cascadia data perform better in this region than the model that did not consider these data. We compare ground-motion variability in Cascadia with the NGA-Subduction model predictions and find differences at short periods ( T = 0.1 s) due to the use of global versus regional data in the development of these models. Moreover, the within-event component of aleatory variability from the GMMs overpredicts the standard deviation of Cascadia recordings at very short periods ( T < 0.05 s). Using global interface earthquakes as a proxy to evaluate the Cascadia GMMs, we find long-period overprediction from a simulation-based GMM and some of the empirical GMMs. When comparing recent observations, we find a similar misfit to GMMs and the 2010 and 2022 Ferndale earthquakes. Finally, we observe different basin amplification factors arising in different subsets of the data, which indicate that differences in basin factors between empirical GMMs could arise from the data selection choices by the developers. As part of evaluating the regional basin terms, we apply basin amplification factors from the magnitude 9 Cascadia earthquake simulations to the empirical GMMs for interface earthquakes. The comparisons presented in this study indicate that the NGA-Subduction GMMs for Cascadia perform well relative to observations and older subduction GMMs.

British Columbia, California, Oregon, Washington

Multi-scale geophysical mapping of the brine and bedrock surfaces along the Dolores River, Paradox Valley, Colorado, December 2023

Total dissolved solids derived from salt dome–sourced brine in the underlying alluvial aquifer substantially increase with distance in the reach of the Dolores River that passes through Paradox Valley in southwestern Colorado. The area has been the site of salinity control operations since the 1990s to reduce salt loading to the downstream Colorado River. Previous airborne and ground/water-based electromagnetic (EM) geophysical data have successfully mapped the top of the brine surface, albeit with relatively coarse near-surface resolution and limited spatial coverage. This present December 2023 study used ground-based high-resolution EM and passive seismic (horizontal-to-vertical spectral ratio, HVSR) tools to map in detail the depth and thickness of the brine zone in the alluvial aquifer (top of the brine down to bedrock contact) in areas immediately surrounding the Dolores River where previous airborne EM (AEM) results indicated brine within 10 m of land surface. Results indicate the deepest bedrock is generally associated with the shallowest brine and local depressions in the collapse breccia (caprock to the Paradox Formation salt) may facilitate vertical migration of brine into the alluvial aquifer. Additionally, the ground-based EM mapping corroborated general patterns in depth to brine that were observed in previous AEM results while also revealing additional detail, including suspected focused brine discharge zones to the Dolores River. A river-based survey wherein EM data, channel depth and river water–specific conductance information were collected augmented these findings and indicated specific areas on both the western and eastern sides of the river where focused brine discharge may occur. This study comprises a large-scale, ground- and water-based geophysical mapping effort, including hundreds of HVSR soundings and 100s of kilometres of EM data, which were successfully translated into digital brine and bedrock surfaces that could be incorporated into groundwater modelling efforts, future well siting or other decision-making.

Colorado