USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Volcanology and Geothermal Research”

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

Hydrological control on soil redox condition and carbon loss of coastal wetland under sea-kevel rise

Coastal wetlands are critical carbon sinks with their biogeochemical and ecological functioning shaped by dynamic hydrological conditions that are increasingly influenced by climate change. A key unresolved question is how hydrologic flow, vegetation response, and rising sea levels interact to regulate soil redox condition and carbon loss in coastal wetlands. Using a field-tested hydrological–biogeochemical–ecological modeling framework, we reveal how the interplay between terrestrial groundwater discharge and tidal fluctuations generates complex groundwater flow patterns at the terrestrial–aquatic interface, and how these patterns modulate soil redox conditions, in turn influencing soil organic matter decomposition and carbon loss. Notably, rising sea levels suppress soil CO 2 emissions while reducing lateral dissolved carbon losses, thereby enhancing litter carbon sequestration under anoxic conditions. As vegetation responds to sea-level rise and carbon inputs diminish, litter carbon subsequently declines. These findings underscore a critical hydrological control on carbon cycling, advancing our understanding of coastal ecosystem resilience in a warming world.

Louisiana

Numerical model of the groundwater-flow system near the southeastern part of Puget Sound, Washington

Groundwater flow in the active model area (AMA) was simulated using a groundwater-flow model. A steady-state model version of the model simulates equilibrium conditions, and a transient model version simulates monthly variability. The model corresponds to the physical and temporal dimensions of the conceptual model and groundwater budget. The steady-state model version represents average conditions for an 11-year period (January 1, 2005–December 31, 2015), and the transient model represents monthly hydrologic variability within that period. The 13-layer model was constructed using MODFLOW-NWT with a uniformly spaced grid consisting of 416 rows, 433 columns, and cells with a horizontal dimension of 500 feet (ft) on a side. The model was calibrated to measured values of water levels in wells and lakes and estimated base flow for selected streamflow measurement stations, commonly referred to as streamgages. Model calibration was accomplished using a combination of manual and automatic methods, including the Model-Independent Parameter Estimation (PEST) program that adjusted model input parameters with the aim of minimizing the difference between estimated and model-simulated values of hydraulic head and base flow. Model boundary conditions consist of all simulated groundwater inflow to and outflow from the AMA. For example, a stream reach that simulates a gain from or loss to groundwater is a boundary condition that allows water to exit or enter, respectively, the groundwater system. Other boundary conditions include springs, seeps, precipitation recharge, groundwater exchange with lakes and Puget Sound, and groundwater pumping. A comparison of the estimated groundwater budget to that simulated by the steady-state model version indicates that the relative percentages of total inflow or total outflow for six major categories of boundary conditions are similar for the two budgets. The model was used to simulate three suites of scenarios of potential drought and water-use changes. Scenario 1 suite consisted of the steady-state model version that was run with 0, 15, 20, and 25 percent reduction of precipitation recharge to assess the corresponding reductions in base flow with decreasing recharge. The last simulation for the scenario 1 suite consisted of the transient model version simulating 3 years of consecutive seasonal drought, defined by the months of May through September, to assess the corresponding base-flow reductions. Scenario 2 suite consisted of the steady-state model version with all simulated groundwater use removed, compared with a simulation that includes current groundwater use to evaluate changes to potentiometric surfaces and base flows. Scenario 3 suite consisted of a transient model version of the model that simulated pumping increases for four different categories of water-supply wells (compared to no pumping increases) to evaluate resulting reductions in base flow. Although, these scenarios provide examples of model applications and useful insights, many other scenarios could be simulated. A description of how to download the model is described in the body of this report. Uncertainty is associated with most model inputs. Groundwater levels, lake levels, and land-surface altitudes are relatively certain; other model inputs are far less certain, including precipitation recharge, base flow, hydraulic properties, water use, and the three-dimensional structure of subsurface hydrogeologic units. Models are useful not because of high levels of accuracy of all model inputs, but because they combine the best information and estimates available, thereby providing the best predictions available related to physical processes. The model described in this report simulates groundwater flow on a regional scale, which has inherent limitations for simulating hydrologic scenarios at local scales. Model structures and inputs were generalized to be consistent with this regional scale. For example, the actual groundwater system has much greater heterogeneity of hydraulic conductivity than is possible within the model’s degrees of freedom. Variations in hydraulic gradients over distances less than 500 ft cannot be simulated. The distances between model features, such as a pumping well and a stream, must be placed at 500-ft intervals and are co-located if both features are within the same model cell.

Washington

A regional model comparison between MODPATH and MT3D of groundwater travel time distributions

Groundwater quality changes in wells and streams lag behind changes to land use due to groundwater travel times. Two contaminant transport methods were compared to assess differences in their simulated travel time distributions (TTDs) to streams and wells in the Wisconsin Central Sands. MODPATH simulates advective groundwater flow with particle tracking, while MT3D simulates age-mass using a finite difference solution without dispersion to allow for direct comparison of the two methods. MODPATH appropriately simulates groundwater TTDs from the water table to surface discharge but is subject to inaccuracies at weak-sink well cells due to the flow-model grid discretization and imprecise location of well discharge within well cells. MT3D better represents weak-sink well cells since it removes mass in proportion to the prescribed pumping rate, although travel time within well cells is neglected. Conversely, MT3D's treatment of surface water boundary cells is not as accurate as MODPATH because mass should be removed from the water table rather than the full cell volume. MT3D simulations of TTDs can also be confounded by the instantaneous vertical distribution of mass introduced throughout recharge cells instead of at the water table, which initiates mass along deeper flow paths. We evaluated 9 MODPATH and 13 MT3D implementations, generating differences in median travel times of up to 18 years. Both methods have strengths and weaknesses, with MT3D better representing weak-sink well cell behavior and MODPATH better representing surficial recharge and discharge. The effect of these characteristics on simulated TTDs, along with ideas for ameliorating method weaknesses, is discussed.

Wisconsin

Integrity of production wells and confining unit at the Naval Weapons Industrial Reserve Plant, Dallas, Texas, 1995

Ground water in the shallow alluvial aquifer is contaminated at the Naval Weapons Industrial Reserve Plant, Dallas, Texas. Five production wells at the site are cased through the alluvial aquifer and underlying units and are screened in either the Paluxy or Twin Mountains aquifers. Three abandoned wells, originally completed in the Twin Mountains aquifer but filled with drilling mud in 1958, also penetrate the alluvial aquifer. The Paluxy and Twin Mountains aquifers are used for drinking-water supplies in and around the Dallas-Fort Worth metroplex. Trichloroethylene and its degradation products, dichloroethylene and vinyl chloride, and the metal chromium previously have been detected in the shallow alluvial aquifer. Current (1995) analyses of water-quality samples taken from the static water column of the five production wells and one of the abandoned wells indicate no trichloroethylene, dichloroethylene, or vinyl chloride in the water column of these wells. Chromium was detected in all samples, but concentrations were less than the practical quantitation limit, which is the regulatory action level for this site. The results of borehole geophysical log analysis indicate that two of the production wells could have vertically connected intervals where cement bonding in the well annulus is poor. The other production wells have overall good bonding. Temperature logs do not indicate flow behind casing except in the screened interval of one well. Geophysical logs show the Eagle Ford Shale ranges from 147 to 185 feet thick at the site. The Eagle Ford Shale has low permeability and a high plasticity index. These physical characteristics make the Eagle Ford Shale an excellent confining unit.

Texas

Status and performance of the ShakeAlert® earthquake early warning system: 2019-2023

The U.S. Geological Survey (USGS)‐operated ShakeAlert® system is the United States West Coast earthquake early warning system ( Given et al. , 2018 ). In this study we detail ShakeAlert’s performance during some of the largest events seen by the system thus far. Statewide public alerting using ShakeAlert messages was authorized in California in October 2019. Over the next few years, public alerts were expanded into Oregon and then into Washington ( U.S. Geological Survey, 2024 ). ShakeAlert source results are routinely compared to the USGS Comprehensive Catalog (ComCat; Guy et al. , 2015 ; U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which contains the earthquake location and magnitude determined using complete waveform data. M 4.5 and larger is the threshold used for public alerting and was deliberately set below the level where damage is likely to compensate for cases where the system underestimates the magnitude. Between 17 October 2019 and 1 September 2023, the ShakeAlert system created 95 events with maximum magnitude estimates of M ≥4.5, the public alerting threshold. 94 of the 95 events were due to real earthquakes. Seven were categorized “false” per ShakeAlert’s internal definition that there was no matching catalog event within 100 km and 30 s of origin time; however, all but one of these were real earthquakes that were poorly located, primarily because they were at the edges of the seismic network. Three detected events were labeled “missed” because they were very poorly located (>100 km location error). In addition, the system did not produce solutions for four ComCat events M ≥4.5 ( U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which were all at the edge of the alerting and network boundaries. The ShakeAlert system has accurately detected the majority of earthquakes that have occurred within the operational region since completing the public rollout, and alerts from the system have been delivered to millions of cell phone users throughout the West Coast.

California, Oregon, Washington

Status and understanding of groundwater quality in the Mojave Basin Domestic-Supply Aquifer study unit, 2018—California GAMA Priority Basin Project

Groundwater quality in the western part of the Mojave Desert in San Bernardino County, California, was investigated in 2018 as part of the California State Water Resources Control Board Groundwater Ambient Monitoring and Assessment Program Priority Basin Project. The Mojave Basin Domestic-Supply Aquifer study unit (MOBS) region was divided into two study areas—floodplain and regional—to assess differences between the two major aquifers used for drinking water supply in the area. This assessment characterized the quality of ambient groundwater and not the quality of treated drinking water. The study included three components: (1) a status assessment, which characterized the quality of groundwater resources used for domestic drinking-water supply in the floodplain and regional study areas; (2) a brief understanding assessment, which evaluated factors that could potentially affect the quality of groundwater used by domestic wells in the region; and (3) a comparative assessment between the groundwater resources used by domestic wells and public-supply wells in the two study areas. The domestic-well assessment was based on data collected by the U.S. Geological Survey from 48 domestic wells in January–May 2018. The public-supply assessment was based on data for samples from 322 public-supply wells in 2008–18, either collected by the U.S. Geological Survey or compiled from the California State Water Resources Control Boards Division of Drinking Water publicly available database. Concentrations of water-quality constituents in ambient groundwater were compared to regulatory and non-regulatory benchmarks typically used by the State of California and Federal agencies as health-based or aesthetic standards for public drinking water. Relative concentrations, defined as the measured concentration divided by the benchmark concentration, were classified as high (greater than 1.0), moderate (greater than 0.5 for inorganic constituents or 0.1 for organic and special-interest constituents, and not high), or low (concentrations lower than moderate). The floodplain and regional study areas were divided into 15 and 35 grid cells, respectively, and grid-based methods were used to compute the areal proportions of the two study areas with high, moderate, or low relative concentrations of individual constituents and classes of constituents. For the domestic-supply assessment, one or more inorganic constituents with health-based benchmarks were detected at high relative concentrations in 58 percent of the regional study area and 13 percent of the floodplain study area. The inorganic constituents with health-based benchmarks detected at high relative concentrations in the regional study area were arsenic, chromium and hexavalent chromium, fluoride, adjusted gross alpha particle activity, uranium, molybdenum, strontium, and nitrate; only arsenic was detected at high relative concentrations in the floodplain study area. One or more inorganic constituents with secondary maximum contaminant level benchmarks were detected at high concentrations in 15 and 6.7 percent of the regional and floodplain study areas, respectively. The constituents detected at high relative concentrations in the regional study area were total dissolved solids, chloride, sulfate, and iron; only total dissolved solids and sulfate were detected at high relative concentrations in the floodplain study area. Organic constituents were not detected at moderate or high relative concentrations in either the regional or floodplain study areas. Volatile organic compounds were detected at low relative concentrations in 21 and 27 percent of the regional and floodplain study areas, respectively, and pesticides were detected at low relative concentrations in 9.1 and 20 percent of the regional and floodplain study areas, respectively. The only individual organic constituent detected in more than 10 percent of either study area was the trihalomethane trichloromethane. Total coliform bacteria were detected in 15 and 27 percent of the grid wells in the regional and floodplain study areas, respectively. The greater prevalence of high relative concentrations of many inorganic constituents in the regional study area compared to the floodplain area likely indicates the greater diversity of geologic material at depth in aquifer material and generally finer-grained alluvium compared to the floodplain study area combined with generally older groundwater that has had more contact time with aquifer materials. In general, trace element concentrations (1) increased with increasing groundwater age, (2) increased with distance from recharge sources in the mountains, and (3) increased with closer proximity to some types of geological units. In general, groundwater from domestic wells in the floodplain study area is young, with most samples containing a component of modern groundwater based on tritium and unadjusted carbon-14 activities, whereas groundwater from domestic wells in the regional study area generally is old, with most samples having unadjusted carbon-14 ages of 5,000–40,000 years. Public-supply wells in MOBS generally were deeper than domestic wells and presumably are in contact with older, more weathered alluvium that may have more mobile trace elements, such as arsenic or uranium. However, only 26 percent of the public-supply regional study area had high relative concentrations of inorganic constituents, compared to 58 percent for the domestic regional study area. The percentages of the public-supply and domestic floodplain study areas with high relative concentrations of inorganic constituents were 11 and 13 percent, respectively. The ages of groundwater used by public-supply and domestic wells in each study area were similar, which was not expected given the greater depth of the public-supply wells. Three potential factors may contribute to these results: (1) greater spatial footprint of domestic well network, which may result in domestic wells pumping groundwater from fractured bedrock or mineralized areas not used by public-supply wells; (2) greater pumping rates in public-supply wells, resulting in more water being withdrawn from coarse-grained, heterogeneous alluvium than finer-grained layers, which may have higher concentrations of (or more mobile) inorganic constituents; and (3) a greater degree of well management with public-supply wells, which may include pausing use of or decommissioning wells if treating or blending water is not feasible to lower constituent concentrations.

California

Evidence of nitrate attenuation in intertidal and subtidal groundwater in a subterranean estuary at a Cape Cod embayment, East Falmouth, Massachusetts, 2015–16

Nitrogen dynamics in intertidal and nearshore subtidal groundwater (subterranean estuary) adjacent to the Seacoast Shores peninsula, Falmouth, Massachusetts, were investigated during 2015–16 by the U.S. Geological Survey. The peninsula is a densely populated residential area with septic systems and cesspools that are substantial sources of nitrogen to groundwater. The study area is in the Eel River, an estuarine saltwater embayment connected to the ocean adjacent to the western shore of the peninsula, that was the subject of an earlier study by Colman and others (2018, https://doi.org/10.3133/sir20185095 ) on nitrogen transport and transformations in groundwater between onshore and offshore locations. The previous study documented the distribution of nitrate concentrations and nitrate attenuation reactions in fresh groundwater beneath the peninsula and the estuary. The current study extended those observations with more detailed sampling and analysis of shallow groundwater from wells near discharge sites beneath the estuary. The current field investigation included sampling of existing wells and installation and sampling of clusters of wells and temporary sampling points in the subterranean estuary, including (1) shallow transects 0.3 to 1.2 meters (m) deep extending from 1 to 13.5 m offshore and (2) deeper wells (from 1.83 to 4.88 m deep) extending from 4.3 to 14.3 m offshore. Measurements of hydraulic-head gradients 2–5 m below the sediment/water interface in the intertidal and nearshore subtidal zones indicated that groundwater flow generally was upwards (towards the estuary) under all tide conditions in October 2016. The magnitude of the gradient was greatest during low tide conditions, indicating that groundwater discharge likely decreased during high tides. Measurements of specific conductance in shallow groundwater in the subterranean estuary in three transects perpendicular to shore were consistent with the existence of saltwater flow cells (infiltration of overlying saline water, mixing with fresh groundwater, and discharge to the overlying saline water) in the intertidal and nearshore subtidal regions. The size of these flow cells was variable in space and time and dependent on the elevation of the tide (spring or neap). At this location in the Eel River subterranean estuary, and offshore to at least 13.5 m, offshore flow of fresh groundwater apparently prevented a deeper saltwater wedge from discharging to the surface. Nitrate concentrations in shallow groundwater (30 to 122 centimeters [cm] depth) were variable in space and time, ranging from not detectable to 600 micromoles per liter (μmol/L) (8.4 milligrams per liter as N), and were highest in June 2016 at depths from 61 to 122 cm below the sediment/water interface and from 4 to 9 m offshore. Nitrate generally was not detectable in saline shallow groundwater at 30-cm depth or at any depth from 30 to 122 cm from 10 to 13.5 m offshore. Dissolved oxygen concentrations were suboxic (less than 16 μmol/L) in 60 percent of the sampled subterranean groundwater beneath the intertidal and subtidal zones. In the remaining sites, the range of dissolved oxygen concentrations was from 18 to 272 μmol/L and the median concentration was 43 μmol/L. Evidence for microbial nitrate reduction (denitrification and possibly anammox) was provided by the distribution of the reaction product nitrogen gas (excess N 2 , or N 2MIC ), as determined from analysis of the dissolved nitrogen gas and argon gas (Ar) concentrations in groundwater samples. Excess nitrogen gas provided evidence for nitrate reduction in shallow groundwater below the subtidal and, to a lesser extent, intertidal zones adjacent to the Seacoast Shores peninsula. These zones, where evidence for nitrate reduction was detected, were in fresh and brackish groundwater near subtidal or intertidal saltwater cells where discharging fresh groundwater mixed with infiltrating saline water. Infiltrating seawater may have supplied organic carbon, one of several potential electron donors that are required for denitrification. Other potential electron donors, such as organic carbon, iron, manganese, hydrogen, methane, ammonium, elemental sulfur, or sulfide phases, may have been supplied by the estuarine sediments. Drainage from surface runoff near the shore also may have supplied organic carbon to fresh groundwater near the intertidal saltwater cell. The highest amounts of nitrate converted to excess nitrogen gas were estimated to be in the range of 230 to 430 μmol/L in nearly fresh groundwater near the subtidal saltwater cell at depths of 61 to 122 cm below the sediment/water interface and from 10 to 13.5 m offshore. Evidence of denitrification within 10 m of the shore was sparse (generally limited to less than 50 μmol/L of N 2 -N) despite the presence of high nitrate concentrations. The spatial distribution of estimated nitrate reduction in the intertidal and nearshore subtidal fresh and brackish groundwater may be related to local variability in the distribution of reactive electron donors in those zones. Variations in the amount of nitrate reduction to nitrogen gas were not clearly related to potential aqueous electron donors such as dissolved organic carbon, nor to potential reaction products such as alkalinity, but may have been controlled by combinations of aqueous and solid-phase reactants. The distribution of relatively shallow fresh groundwater containing nitrate could indicate potential nitrate discharge areas in the lower intertidal zone and uncertain locations farther offshore; however, the data did not extend all the way to the sediment/water interface or to the offshore freshwater limit. This study confirmed substantial loss of nitrate from some of the fresh and brackish groundwater in shallow subestuarine sediments prior to discharge but did not quantify how much nitrate eventually discharged to the estuary.

Massachusetts

Preferential groundwater discharges along stream corridors are disregarded sources of greenhouse gases

Groundwater delivery of greenhouse gases (GHGs) to stream banks and riparian areas, before mixing with surface waters, has not been well quantified. We measured preferential groundwater delivery of GHGs to stream banks within three stream reaches, and found that stream banks with discharging groundwater emitted more CO 2 and were sources of N 2 O compared to stream banks without actively discharging groundwater, which emitted less CO 2 and were N 2 O sinks. At one of our stream reaches, groundwater CO 2 and N 2 O concentrations were 1.4–19.2 and 1.1–40.6 times higher than those in surface water, respectively, and groundwater delivery rates of CO 2 and N 2 O were 1.5 and 1.6 times higher than surface water emissions per unit area. On average, 21% (range 0%–100%) of CO 2 and N 2 O were emitted at the stream bank before mixing with surface waters. Preferential groundwater GHG emissions may contribute substantially to stream corridor emissions and may be underestimated when using a channel-centric approach to estimate riverine GHG budgets.

JGR Biogeosciences

Biodegradation of chlorinated ethenes at a karst site in middle Tennessee

This report presents results of field and laboratory investigations examining the biodegradation of chlorinated ethenes in a karst aquifer contaminated with trichloroethylene (TCE). The study site, located in Middle Tennessee, was selected because of the presence of TCE degradation byproducts in the karst aquifer and available site hydrologic and chlorinated-ethene information. Additional chemical, biological, and hydrologic data were gathered to evaluate whether the occurrence of TCE degradation byproducts in the karst aquifer was the result of biodegradation within the aquifer or simply transport into the aquifer. Geochemical analysis established that sulfate-reducing conditions, essential for reductive dechlorination of chlorinated solvents, existed in parts of the contaminated karst aquifer. Other areas of the aquifer fluctuated between anaerobic and aerobic conditions and contained compounds associated with cometabolism, such as ethane, methane, ammonia, and dissolved oxygen. A large, diverse bacteria population inhabits the contaminated aquifer. Bacteria known to biodegrade TCE and other chlorinated solvents, such as sulfate-reducers, methanotrophs, and ammonia-oxidizers, were identified from karst-aquifer water using the RNA-hybridization technique. Results from microcosms using raw karst-aquifer water found that aerobic cometabolism and anaerobic reductive-dechlorination degradation processes were possible when appropriate conditions were established in the microcosms. These chemical and biological results provide circumstantial evidence that several biodegradation processes are active in the aquifer. Additional site hydrologic information was developed to determine if appropriate conditions persist long enough in the karst aquifer for these biodegradation processes to be significant. Continuous monitoring devices placed in four wells during the spring of 1998 indicated that pH, specific conductance, dissolved oxygen, and oxidation-reduction potentials changed very little in areas isolated from active ground-water flow paths. These stable areas in the karst aquifer had geochemical conditions and bacteria conducive to reductive dechlorination of chlorinated ethenes. Other areas of the karst aquifer were associated with active ground-water flow paths and fluctuated between anaerobic and aerobic conditions in response to rain events. Associated with this dynamic environment were bacteria and geochemical conditions conducive to cometabolism. In summary, multiple lines of evidence developed from chemical, biological, and hydrologic data demonstrate that a variety of biodegradation processes are active in this karst aquifer.

Tennessee

Eruptive history of northern Harrat Rahat—Volume, timing, and composition of volcanism over the past 1.2 million years

Harrat Rahat, one of several large, basalt-dominated volcanic fields in the western part of the Kingdom of Saudi Arabia, is a prime example of continental, intraplate volcanism. Excellent exposure makes this an outstanding site to investigate changing volcanic flux and composition through time. We present 93 40 Ar/ 39 Ar ages and 6 36 Cl surface-exposure ages for volcanic deposits throughout northern Harrat Rahat that, integrated with a new geologic map, define 12 eruptive stages. Exposed volcanic deposits in the study area erupted less than 1.2 million years ago (Ma), and 214 of 234 identified eruptions occurred less than 570 thousand years ago (ka). Two eruptions were in the Holocene, including a historically described basaltic eruption in 1256 C.E. and a trachyte eruption newly recognized as Holocene (4.2±5.2 ka). An estimated approximately 82 cubic kilometers (km 3 ; dense rock equivalent) of volcanic products can be documented as having erupted since 1.2 Ma, though this is a lower limit because of concealment of deposits older than 570 ka. Over the last 570 thousand years (k.y.), the average eruption rate was 0.14 cubic kilometers per thousand years (km 3 /k.y.), but volcanism was episodic with periods alternating between low (0.04–0.06 km 3 /k.y.) and high (0.1–0.3 km 3 /k.y.) effusion rates. Before 180 ka, eruptions vented from the volcanic field’s dominant eastern vent axis and from a subsidiary, diffuse, western vent axis. After 180 ka, volcanism focused along the eastern vent axis, and the composition of volcanism varied systematically along its length from basalt dominated in the north to trachyte dominated in the south. We hypothesize that these compositional variations younger than 180 k.y. reflect the growth of a mafic intrusive complex beneath the southern part of the vent axis, which led to the development of evolved magmas. Lastly, these new age data allow for a reassessment of the volcanic recurrence interval at northern Harrat Rahat. Based on available data, volcanism in northern Harrat Rahat over the last 180 k.y. is poorly described using a Poisson distribution with a single recurrence interval. Instead, data for northern Harrat Rahat are better described using a mixed exponential distribution that is applicable for volcanic systems characterized by two different eruptive states, where one state with a longer recurrence interval corresponding to periods of low eruption frequency and one state with a shorter recurrence interval corresponding to periods of high eruption frequency. The preferred model for northern Harrat Rahat over the last 180 k.y. uses a long recurrence interval of 4.0 k.y. and a short recurrence interval of 0.22 k.y.

Professional Paper

Ground deformation and gravity for volcano monitoring

Introduction When magma accumulates or migrates, it can cause pressurization and related ground deformation. Characterization of surface deformation provides important constraints on the potential for future volcanic activity, especially in combination with seismic activity, gas emissions, and other indicators. A wide variety of techniques and instrument types have been applied to the study of ground deformation at volcanoes (sidebar, p. 2; Dzurisin, 2000, 2003, 2007). Geodetic instruments include continuously recording Global Navigation Satellite System (GNSS; of which the United States’ Global Positioning System is one example) stations (fig. D1), borehole tiltmeters, and interferometric synthetic aperture radar (InSAR) measurements (from satellites, occupied and unoccupied aircraft systems, and ground-based sensors). Additional geodetic measurements like continuous- and survey-mode gravity (fig. D2) can contribute substantially to interpreting these data. Borehole strainmeters (see chapter K , this volume, by Hurwitz and Lowenstern, 2024) also have outstanding utility for monitoring deformation, although because of cost and permitting challenges, we do not include them as part of standard volcano monitoring networks for U.S. volcanoes. Still other techniques like light detection and ranging (lidar), structure from motion, and optical satellite data can be used to derive gross topographic changes, which can be used to map volcanic deposits, infer eruption rates, and gain insights into the source processes associated with eruptive activity (see chapter G , this volume, on tracking surface changes caused by volcanic activity; Orr and others, 2024). Experience has shown that no single geodetic monitoring technique is adequate to detect and track the entire range of ground-motion patterns that occur at volcanoes, primarily because of the temporal and spatial diversity of volcano deformation (fig. D3). Similarly, the magnitude of surface deformation varies widely. Geodetic monitoring strategies should therefore include multiple techniques and instrument types to cover a wide range of spatial and temporal scales. In identifying recommendations for geodetic instrumentation for volcano monitoring networks, we attempted to maximize the diversity of instrument types to measure the full range of deformation signals and minimize their expense and number; thus, we do not include several well-known deformation-monitoring techniques in our recommendations. Extensometers, for example, measure strains over distances of a few meters and have an excellent record of success in detecting changes in preeruptive localized ground motion across existing cracks, including at Mount St. Helens, Washington (Iwatsubo and others, 1992), and Piton de la Fournaise, Réunion Island (Peltier and others, 2006). Despite being relatively inexpensive, extensometers are best used primarily when localized ground displacements (for example, ground cracks) need to be tracked, and are not necessary at all volcanoes. In considering volcano deformation monitoring strategies, two complicating factors are deserving of special attention. First, not all deformation is driven by subsurface magmatic activity—for example, at many large stratovolcanoes (for example, Mount Rainier), flank collapses and landslides are significant geologic hazards (Reid and others, 2001) that may occur even in the absence of magmatic activity. Monitoring the stability of volcanoes is thus another critical application of geodetic monitoring networks to inform hazard assessment. One of the most famous examples of edifice instability is the large flank collapse that initiated the May 18, 1980, eruption of Mount St. Helens. Deformation monitoring had detected a bulge on the north flank of the mountain in April 1980 that was expanding by several meters per day (Lipman and others, 1981). Given that flank collapses can happen at any time during a period of volcanic unrest (or even outside a period of unrest), the capability to assess edifice stability is critical. Second, although volcanoes are commonly treated as idealized structures that erupt from single points, like centralvent stratovolcanoes, many are characterized by long rift zones from which eruptions may originate, and distributed volcanic fields are characterized by broadly spaced vents. For example, linear dikes are common at Kīlauea, Mauna Loa, and between Mount Shasta and Medicine Lake in California. At Kīlauea, one of these linear dikes emerged more than 40 kilometers (km) away from the summit of the volcano during the lower East Rift Zone eruption in 2018. Other volcanic fields, like Lassen volcanic center, California, or the San Francisco Volcanic Field, Arizona, have many small vents spread over a wide area. Although the instrumentation guidelines presented in this chapter remain phrased for central-vent volcanoes, they should be modified as needed in the context of the eruptive characteristics of each individual volcanic system. Spatial analysis of geodetic network coverage could help to ensure adequate instrumentation in areas where volcanism can occur over a broad area as opposed to a central vent. As an example, consider the adjacent volcanoes Mount Shasta and Medicine Lake. If station locations are chosen based only on the distance from the centers of the volcanoes, then any geodetic anomalies between the two volcanoes—an area of potential volcanism as indicated by the presence of volcanic features—may remain undetected by ground-based instrumentation. The spatial analysis is accomplished via a grid of pressure point sources (Mogi, 1958) evenly distributed across the map area, at a depth of 5 km in this example (fig. D4). Each source is inflated until predicted deformations exceed the GNSS white noise uncertainty estimates at one site (Langbein, 2017; Murray and Svarc, 2017). This volume of detectable magma provides a measure of the quality of the coverage (fig. D4). The results indicate that, as of 2022, there is a large area between Mount Shasta and Medicine Lake volcano with existing mapped dikes in which a substantial amount of magma could intrude without being detected geodetically. Applying this style of analysis to individual volcanic systems can provide a guide for designing network geometry given the expected locations of future eruptions.

Scientific Investigations Report

High-precision earthquake catalog for Minto Flats fault zone, central Alaska, reveals complex and conjugate faulting

The Minto Flats fault zone (MFFZ) in central Alaska is a left‐lateral strike‐slip fault system situated between the continental‐scale right‐lateral Denali and Kaltag‐Tintina faults. The MFFZ has the potential to generate magnitude 7 earthquakes, and it hosted a magnitude 6 earthquake in 1995. It has also produced exotic events, such as very‐low‐frequency earthquakes and nucleation signals. We use network‐matched filtering and relative earthquake relocation techniques to derive a detailed catalog of earthquake locations for the MFFZ. The catalog spans from August 2014 to December 2019, a time period including 13 temporary seismic stations in the region. Our results provide the most complete catalog for the MFFZ and include deeper events, clusters of shallow seismicity, and a complex and segmented fault structure not observed in the original regional catalog. We document right‐lateral strike‐slip faulting, conjugate to the main northeast‐striking left‐lateral faults of the MFFZ. Below Nenana basin, the relocated seismicity reveals northwest‐dipping left‐lateral faults, supporting the inference that deep crustal active faulting is associated with recent basin deformation.

Alaska

Overview of the Cenozoic geology of the northern Harrat Rahat volcanic field, Kingdom of Saudi Arabia

The Harrat Rahat volcanic field, located in the west-central part of the Kingdom of Saudi Arabia, is one of the larger Cenozoic harrats among the more than 17 harrats situated upon the Arabia Plate. The map plate contained herein shows, at a scale of 1:100,000, the mapped volcanic geology of northern Harrat Rahat, which consists of the northernmost one-fifth of Harrat Rahat. Northern Harrat Rahat has an area of about 3,340 square kilometers (km 2 ), of which 2,567 km 2 is covered by Harrat Rahat volcanic rocks and deposits, and it encompasses more than 900 exposed vents (that is, craters, cryptodomes, maars, and scoria cones), 289 of which are isolated by younger volcanic deposits and have not been correlated with the 234 volcanic rock units distinguished by geologic mapping. Northern Harrat Rahat is of special interest owing to its proximity to the city of Al Madinah al Munawwarah, which sits within, and is continuing to expand southward over, the north end of the volcanic field. Al Madinah is home to an expanding population, currently at more than 2 million residents, together with the intermittent addition of approximately 3 million pilgrims during Hajj and Umrah (religious visitations). The center of Al Madinah is less than 8 km from lava flows of the only confirmed historically documented eruption, which occurred in 1256 C.E. (654 A.H.). Earlier prehistoric lava flows also encroached into the area of the present-day city limits, as demonstrated by volcanic rocks exposed widely throughout the city in roadcuts, parks, and excavations for new buildings, although no evidence has been found of any other than the 1256 C.E. lava having reached that area during times of human habitation. Geologic mapping was undertaken by the U.S. Geological Survey in collaboration with the Saudi Geological Survey. The features of primary interest within the map area are scoria cones, lava flows, lava domes, craters, and pyroclastic deposits from the Quaternary, which have compositions of basalt, hawaiite, mugearite, benmoreite, and trachyte. The geologic mapping was published by Downs and others (2019) on a single sheet at 1:75,000 scale and two detailed sheets at 1:25,000 scale, accompanied by detailed explanations of the geology. This report presents the geology on a single sheet at 1:100,000 scale and provides condensed geologic explanations for the convenience of readers of this volume. Some minor errors of nomenclature and lava source regions that do not change fundamental interpretations have been corrected herein.

Professional Paper

Hydrogeologic framework of the Mountain Home area, southern Idaho

In the arid western Snake River Plain around the City of Mountain Home, Idaho, declining groundwater levels concern agricultural, municipal, and other water users who rely on groundwater for sustenance because surface-water resources are limited. The U.S. Geological Survey developed this hydrogeologic framework to provide an updated characterization of groundwater resources in the western Snake River Plain around the City of Mountain Home. The hydrogeologic framework comprises: (1) a conceptual description of hydrogeologic units, (2) a three-dimensional hydrogeologic model and borehole database, (3) a map of groundwater levels and change, and (4) a discussion of groundwater occurrence and movement within the study area. Hydrogeologic units were defined based on existing literature and the borehole database compiled for this study; the five hydrogeologic units are granite, rhyolite, basalt, fine-grained sediments, and coarse-grained sediments. Each unit can bear water, but the main regional aquifer in the study area occurs in the basalt and fine-grained sediment units with depth to water ranging from 150 to 765 feet. A perched groundwater zone near the City of Mountain Home is primarily hosted in basalt and used domestically with most depths to water ranging from 30 to 100 feet. Interflow zones, scoria, and vertical fractures create heterogeneity within the basalt hydrogeologic unit that exerts strong control on groundwater movement, creating horizontal perching conditions and zones of enhanced vertical conductivity that facilitate downward groundwater percolation. In the fine- and coarse-grained sediments and rhyolite units, inferred faults both impede and enhance groundwater movement. The borehole database was constructed by digitizing 540 well-driller reports and was used to build a three-dimensional hydrogeologic framework model which reasonably represents the spatial distribution of hydrogeologic units in the study area. Generally, fine-grained sediments underlie much of the study area, with basalt concentrated in the central and western study area and rhyolite and granite in the uplands to the north. Groundwater levels were measured in 180 wells in March and November 2023; these data were used to develop water-table contour maps and describe groundwater-level change over an irrigation season. Groundwater generally flows south-southwest to the Snake River and groundwater levels declined across most of the study area (from 0.03 to 22.01 feet) between spring and autumn 2023, which is consistent with long-term declines in the Cinder Cone Butte Critical Groundwater Area and Mountain Home Groundwater Management Area. Groundwater levels rose (0.6 to 15.44 feet) over the irrigation season in most wells in the perched groundwater zone near the City of Mountain Home and near the Snake River, indicating the importance of surface-water recharge to groundwater in areas where surface water irrigation occurs. In aggregate, this hydrogeologic framework provides an updated characterization of and new insights into groundwater resources in the study area to help inform water resources management.

Idaho

An inset groundwater-flow model to evaluate the effects of layering configuration on model calibration and assess managed aquifer recharge near Shellmound, Mississippi

The U.S. Geological Survey has developed a high-resolution inset groundwater-flow model in the Mississippi Delta as part of an interdisciplinary collaboration coordinated by the Mississippi Alluvial Plain project to provide a tool that stakeholders can use to support water-resource management decisions. Groundwater withdrawals from the Mississippi River Valley alluvial (MRVA) aquifer have been vital to support agricultural production in the region, but substantial groundwater-level declines near Shellmound, Mississippi, have caused concerns for long-term sustainability of the aquifer. To better understand the subsurface and try to mitigate the long-term groundwater-level declines, stakeholders have undertaken actions including a Groundwater Transfer and Injection Pilot (GTIP) project using a riverbank filtration-based managed aquifer recharge approach. The pilot project consisted of extracting groundwater near the Tallahatchie River and reinjecting it into the aquifer 3 kilometers west where water levels have substantially declined. A high-resolution airborne electromagnetic (AEM) survey was also completed to collect electrical resistivity data to support the GTIP project and the development of the groundwater model. The inset groundwater-flow model was developed to (1) integrate the AEM data into the optimal layering configuration of the MRVA aquifer that the available observation data can support through calibration, and (2) assess the potential effect of the GTIP project on the groundwater levels. The AEM data were processed into three different layering configurations leading to the development of model A (18 layers), model B (16 layers), and model C (8 layers), all at a 100- x 100-meter cell spatial resolution using the U.S. Geological Survey modular finite-difference flow model 6 code with Newton-Raphson formulation. The model development process integrated recent advances in modeling, such as the incorporation of AEM data, the use of outputs from the soil-water-balance (SWB) model, and the Aquaculture and Irrigation Water-Use Model, and was facilitated by robust automation using the open-source python packages Modflow-setup and SFRmaker. Using Parameter Estimation ++ Iterative Ensemble Smoother, the three numerical groundwater-flow models (models A, B, and C) were calibrated against a set of observations, which included aquifer groundwater levels, streamflows, stream stage, and aquifer transmissivity. Results indicate that the detailed representation of MRVA aquifer layers in model A produced the best calibrated model by history matching, and the integration of data representing surficial connectivity played a key role in improving groundwater recharge and enhancing the ability of the model to match groundwater levels in the cone of depression. A forecast model simulated the managed aquifer recharge approach, and the results indicated that, given average irrigation and recharge conditions (2010–15), the GTIP project has the potential to induce groundwater-level increases of as much as 3 meters around the injection site, but a sustained increase would require repetition in subsequent years of water transfer at 2022 rates or above.

Mississippi

Complex landslide patterns explained by local intra-unit variability of stratigraphy and structure: Case study in the Tyee Formation, Oregon, USA

Lithology and geologic structure are important controls on landslide susceptibility and are incorporated into many regional landslide hazard models. Typically, metrics for mapped geologic units are used as model input variables and a single set of values for material strength are assumed, regardless of spatial heterogeneities that may exist within a map unit. Here we describe how differences in bedding thickness, grain size, inferred uniaxial compressive strength, and bedding dip control the inherent susceptibility of slopes to deep-seated failure within a single mapped geologic unit - the Tyee Formation of Oregon, USA. The Tyee, which covers over 15,000 km2 and underlies much of the Oregon Coast Range, comprises gently folded alternating beds of sandstone and siltstone deposited as turbidites, forming a 2-km thick Eocene submarine fan which has been uplifted and exhumed through the Cenozoic. Deep-seated landslides are widespread in the Tyee, but form a complex spatial pattern such that landslide density ranges from 0 to 24% of the total landscape area. These slides are often extensive and sufficiently deep to reduce local hillslope gradients, resulting in a strong negative correlation between landslide density and mean local slope. Mean annual precipitation and predicted strong ground motions from Cascadia earthquake scenarios also fail to explain the spatial distribution of deep-seated landslides. Consequently, landslide stability models, which are strongly influenced by landscape slope, pore-water pressure, and seismic acceleration, yield landslide susceptibility maps which are broadly anti-correlated with mapped deep-seated landslide density. Through a multivariable linear regression model, we show that much of the variance in deep-seated landslide density can be explained by variability of intra-unit stratigraphic and structural characteristics, which we measure at 128 sites across two study areas totaling ∼3000 km2. Our results suggest bedding dip is only weakly correlated to landslide density, but strongly influences landslide failure style. Subtle increases in bedding dip, even in the gently folded Tyee Formation, result in a substantially higher likelihood of a landslide being cataclinal, or parallel to bedding. Overall, we find a slight majority of landslides fail within these cataclinal slopes, and that these landslides tend to be larger than non-cataclinal landslides. We also show that the lithological and structural properties that influence landslide susceptibility are distinct for these two populations of landslides. Our results demonstrate how localized, intra-unit, geologic variability can exert strong control on landslide susceptibility and failure style. This suggests that in some locations, landslide hazard models could be significantly improved by incorporating detailed, spatially variable, geologic properties rather than relying solely on generalized geologic map units.

Oregon

Microclimate mediates the strength and direction of avian biotic interactions

Theory predicts that that the strength and direction of species interactions can shift from being competitive in benign environments toward being facilitative in stressful environments. However, the environmental context dependency of species interactions has rarely been tested in animal communities. We capitalized on a 15-year, landscape-scale dataset, collected annually in a relatively stable old-growth forest environment to test the long-held hypothesis that the strength and direction of species interactions might be mediated by climatic conditions. It is generally accepted that competitive and facilitative interactions drive the distributions of many species. Using multi-species dynamic occupancy models applied to long-term data, we tested whether annual settlement by bird species could affect either the persistence or settlement by other phylogenetically related species, and whether these interactions are mediated by microclimate. We found that species interactions were influenced by microclimate for some, but not all avian species pairs. Related species pairs more often showed settlement dynamics that were indicative of attraction rather than repulsion. In some cases, competitive interactions at warmer microclimates flipped to become facilitative in colder ones. However, the reverse was also true: facilitative interactions amplified for some species under warm conditions. Furthermore, species pairs that were closely related were more likely to exhibit competitive relationships along at least part of the microclimatic gradient. Our results highlight the importance of using long-term data to incorporate competitive and facilitative interactions into species distribution models and support the notion that the strength and direction of species interactions can be dependent on microclimatic environmental conditions.

Oregon

Simulation-based scenario ShakeMaps for large magnitude (MW6.5+) crustal earthquakes on the Seattle, Tacoma, and southern Whidbey Island faults, Washington, USA

Scenario ground‐motion maps based on empirical ground‐motion models (GMMs) provide a rapid and generally reliable means of estimating the amplitude and distribution of earthquake shaking. However, because GMMs are designed for broad applicability, they often rely on simplified representations of Earth structure, which can limit their accuracy in regions with complex source, path, and site effects. This can substantially impact the accuracy of predicted shaking in areas like western Washington State, where deep, interconnected basin structure exerts a strong influence on seismic‐wave propagation. In this study, we present a new suite of simulation‐based scenario ShakeMaps that characterize ground shaking from large‐magnitude ( ⁠ M W 6.5–7.5) earthquakes on the Seattle, Tacoma, and southern Whidbey Island faults. These maps are developed using results from recent 3D wave propagation simulations ( Stone et al. , 2022 , 2023 , 2025 ) that incorporate realistic rupture geometries, variable slip distributions, and a regional 3D seismic velocity model with shallow soils. Broadband ground motions are estimated by combining the low‐frequency (<1 Hz) deterministic seismograms from these studies with high‐frequency (1–10 Hz) stochastic seismograms. Simulated ground motions are corrected to account for the enforced minimum shear‐wave velocity and nonlinear site response. The resulting ShakeMaps represent median ground‐shaking estimates derived from multiple rupture scenarios with varying slip distributions and hypocenter locations for each fault. To extend ShakeMap coverage beyond the simulation domain (i.e., into eastern Washington, northern Oregon, and southwestern British Columbia), we scale GMM‐based ground‐motion estimates using amplification patterns observed in the simulations. These new ShakeMaps reveal the substantial influence of deep basin structure on shaking intensity, underscoring the importance of considering crustal structure complexity in regional hazard assessments for the Pacific Northwest.

Washington