USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Analytical and Applied Pyrolysis”

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

Global geologic map of Europa

Discovered by Galileo Galilei more than 400 years ago and imaged in detail by the Voyager 2 Galileo spacecraft, Jupiter’s icy moon Europa has been a source of intrigue. A range of science investigations indicate that it contains the key ingredients for habitability, notably energy, chemistry, and liquid water. Europa’s surface is geologically complex and, based on the dearth of impact craters, interpreted to be as young as ~60 Ma. The array of geologic features that characterize the surface include extensive ridged plains, regions of broad disruption termed chaotic terrain, long, quasi-linear ridges that span thousands of kilometers, and bands as much as 60 kilometers wide and that extend 100s of kilometers. These features, along with other geophysical measurements, indicate the presence of a global briny liquid water ocean beneath the ice shell. It was not until the arrival of the Galileo spacecraft in 1995 that the true nature and level of complexity of the surface was revealed. Although image data returned by Galileo provided insight into the structure of a variety of regions, the entire satellite has yet to be observed at a regional scale (less than 250 meters per pixel) and the detailed geologic nature of much of its surface remains a mystery. Establishing the global context of the distribution and timing of Europan geologic units forms a basis to understand regional and local scale processes, serves as a tool for the planning of future missions, and most of all is essential to gaining insight into the potential habitability of this icy world.

Scientific Investigations Map

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

U.S. Geological Survey science strategy to address white-nose syndrome and bat health in 2025–2029

Since its discovery in 2006, the fungal disease known as white-nose syndrome (WNS) has killed millions of bats. Of the 47 bat species native to the conterminous United States, Alaska, Hawaii, and Canada, 12 have been affected by WNS, including 3 endangered species and 1 proposed endangered species. WNS has also been detected in 40 States and 9 Canadian Provinces. U.S. Geological Survey (USGS) scientists have been critical in identifying the causal fungus for WNS ( Pseudogymnoascus destructans [Pd]), characterizing the effects of WNS, and tracking the spread of Pd in many bat populations in North America. The mission of the USGS WNS and Bat Health Science Team is to deliver integrated science in order to build resiliency into free-ranging bat populations through more effective WNS management, build capacity for bat health science, and enhance bat health information sharing across USGS science centers and cooperative research units as well as with stakeholders. The USGS can play an important role in supporting regional and national capacity building by providing resources and guidance to local, State, and Tribal management entities and by providing tools to enhance disease management. The USGS Ecosystems Mission Area’s Biological Threats and Invasive Species Research Program is the lead Federal program for free-ranging wildlife disease research and surveillance. As of 2024, guided by the science priorities set by the WNS Steering Committee, USGS scientists are engaged in a nationwide response to WNS. This work is done in close coordination with our partners at the U.S. Fish and Wildlife Service, National Park Service, Bureau of Land Management, U.S. Forest Service of the U.S. Department of Agriculture, U.S. Department of Defense, as well as State and Tribal agencies. In addition to conducting WNS research, the USGS is mapping the spread of WNS and coordinating the North American Bat Monitoring Program (NABat) to understand how WNS and other stressors affect the status and trends of native bats across their range. The USGS is supporting the national WNS response through four science goals: (1) provide situational awareness on the health of bat populations; (2) conduct ecological studies of bats along the gradient of disease vulnerability; (3) contribute actionable science to enhance the resiliency of bat populations; and (4) implement an adaptive, holistic approach to bat health.

Circular

Dissolved arsenic concentrations in surface waters within the upper portions of the Klamath River Basin, Oregon and California

Arsenic toxicity is an environmental health problem. Levels of arsenic in surface waters at some locations in the Klamath River Basin in southern Oregon and northern California can exceed the U.S. Environmental Protection Agency (EPA) standard for drinking water. There are both anthropogenic and natural sources of arsenic. The Klamath River Basin consists primarily of volcanic deposits and contains an underground geothermal system with hot springs and warm water wells, all known natural sources of arsenic. Anthropogenic sources of arsenic are related to the agricultural use of herbicides, fungicides, and insecticides. Surface water arsenic levels can also be affected by fertilizer amendments, evaporative concentration, oxygen-level depletion, and various geochemical transformations that can increase arsenic mobilization. In this study by the U.S. Geological Survey and the Bureau of Reclamation, dissolved concentrations of arsenic, copper, and lead were measured in surface waters at 39 unique sites within the upper portions of the Klamath River Basin between 2018 and 2022. In every year, except 2022, sites were sampled four times between April and November. Surface-water arsenic concentrations varied up to four-orders of magnitude among sites. Median arsenic concentration was lowest at Cherry Creek (0.03 micrograms per liter [μg/L]) and highest at Wood Kimball Spring (36.7 μg/L), two sites located north of Upper Klamath Lake. The highest arsenic concentrations (17.4±4.9 μg/L, n =3) were found in drain sites (defined here as a waterbody returning used irrigation water) while the lowest arsenic concentrations were found in an artesian well (0.8 μg/L, n =1). The elevated arsenic concentrations of the drain sites suggest that arsenic might be concentrated or mobilized by agricultural activities, water re-use practices, and (or) by geochemical processes occurring around water stored in drains (that is, in the water column and across sediment water boundaries). A source of arsenic in drain water in the Klamath Strait Drain area includes water used for irrigation originating from Ady Canal. Other potential sources include groundwater, geothermal water, and local soils and sediments. Seasonal differences in surface-water arsenic concentrations were detected at 13 sites, 10 of which had higher arsenic concentrations in summer than in either spring or fall. The sites sampled around Upper Klamath Lake, the impounded rivers, one of the two canal sites, and 5 of the 14 river sites had higher surface-water arsenic concentrations in the summer than in either spring or fall. Surface-water arsenic concentrations from groundwater sources (that is, springs and in the artesian well) did not vary significantly among seasons (p-values greater than 0.1). Median surface-water concentrations of copper and lead ranged from 0.03 to 3.7 μg/L, and from 0.013 to 0.175 μg/L ( n =2–18), respectively. Dissolved concentrations of both metals were below acute toxicity endpoints reported by the EPA for freshwater animals. Surface-water arsenic concentrations varied independently from corresponding changes in surface-water lead or copper concentrations. However, arsenic concentrations measured in bed-sediment samples collected from a subset of sites located north of Upper Klamath Lake correlated strongly and significantly with the corresponding sedimentary lead concentrations ( p =0.015). Aqueous arsenic speciation measured in a subset of sites in 2019 and 2022 showed that all the arsenic existed as arsenic (V), the most oxidized arsenic species, and presumably, the least toxic. The highest proportions of arsenite (As(III)), the presumably most toxic arsenic species, relative to total arsenic concentrations were found at drain sites. Our assessment of dissolved arsenic concentrations in various surface-water bodies in the Upper Klamath River Basin reveals geographical areas of consistently low (below 2.1 μg/L), moderate (below 10 μg/L) and high (above 10 μg/L) surface-water arsenic concentrations. South of Upper Klamath Lake, surface-water arsenic concentrations were consistently higher than 20 μg/L at two drain sites located in an area of predominant agricultural land use with extensive water re-use practices. North of Upper Klamath Lake, surface-water arsenic concentrations greater than 20 μg/L were consistently measured at sites with limited nearby agricultural activities, suggesting a geogenic source. The consistently high arsenic levels from the Wood River at Jackson F. Kimball State Park, Fort Creek, and Crooked Creek, which are sites located at or near headwater spring sources, suggest a natural background source of arsenic. Water flowing downstream from this area could be a potential source of arsenic to Upper Klamath Lake and the Upper Klamath River.

California, Oregon

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

A scientific vision and roadmap for earthquake rupture forecast developments, a USGS perspective

We articulate a scientific vision and roadmap for the development of improved Earthquake Rupture Forecast models, which are one of the two main modeling components used in modern seismic hazard and risk analysis. One primary future objective is to provide fully time-dependent models that include both elastic rebound and spatiotemporal clustering nationwide, which is particularly important for shorter-term hazard and risk considerations (e.g., earthquake insurance products). We also discuss the importance and perennial challenges associated with quantifying epistemic uncertainties, including those associated with deformation-model slip rates, un-quantified sampling errors with respect to off-fault seismicity, and any spatial covariances. The need for more physics-based approaches is also emphasized, as is the benefit of adding model valuation (quantifying usefulness) to our verification and validation protocols. Given the multidisciplinary and system-level nature of this activity, modular design is critical. Future updates will also draw from best-available science by both the United States Geological Survey and the external community. The primary goal of this paper is to highlight plans that guide research and facilitate community engagement with model development, especially with respect to lowering the entry barrier for early career scientists and engineers. The paper is written so readers can focus on the sections that interest them most (see table of contents), with the Introduction and Discussion providing a stand-alone overview and summary.

Bulletin of the Seismological Society of America

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

Seismicity zoning at Coso geothermal field and stress changes from fluid production and migration

The Coso geothermal field is a major geothermal power production site in the western United States. It has been observed that low-magnitude seismic events ( M < 3.71) are unevenly distributed in three distinct zones, namely, nearfield (<3 km), midfield (3–6 km), and farfield (> 6 km) from the Coso geothermal plant. These zones exhibit distinct changes in earthquake location before and during geothermal production episodes that began in 1986. After 1986, the midfield region of the main flank experiences a significantly lower seismicity rate than the surrounding areas before production episodes. During 2014–2019, the farfield earthquakes cluster in the eastern and western parts of the greater Coso area, which is discernably different from how those pre-production earthquake events were distributed along the conjugate NW-SE and SW-NW trending structures across the main flank. Here, we analyze the stage of stress with finite-element-based poroelastic simulations to illustrate how the spatiotemporal evolution of the seismicity is associated with the pattern of stress perturbations caused by fluid migration amid the operations of geothermal power plants. Generally, ∼70% of co-production seismicity is found in zones of increased Coulomb stress between 2014 and 2019 at >99% confidence. Meanwhile, the midfield zone of seismic paucity overlaps with the zone of decreasing pore-fluid pressure. Overall, the results provide a physical explanation of how decadal geothermal operations at Coso have perturbed stress-field changes and contributed to the evolving characteristic seismic pattern, shedding insights into assessing the seismic hazard in other geothermal settings.

California

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

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

Colored shaded relief bathymetry and acoustic backscatter of Ozette Lake, Washington

Offshore of the Pacific Northwest of the United States is the Cascadia Subduction Zone, a 1,000-kilometer-long tectonic boundary defined by a large fault, called a megathrust, that extends from the Mendocino Junction off northern California to the Nootka Fracture Zone off Vancouver Island, Canada (U.S. Geological Survey, 2023). The Juan de Fuca and Gorda oceanic plates to the west of this boundary subduct under the North America continental plate to the east. Several other smaller faults that cut through the North America plate crust also affect the region. Although their effects upon Ozette Lake are uncertain, geological evidence for past earthquakes, such as underwater landslides, may be found in Pacific Northwest lakes. Underwater landslides caused by past earthquakes should be well preserved in these relatively undisturbed lake environments. The floor of Ozette Lake, Washington, located along the Pacific coast of the United States, west of the Puget Sound region and about 140 kilometers east of the megathrust was mapped by the U.S. Geological Survey in July of 2019 to search for evidence of past earthquakes. Mapping was completed using a SWATHplus-M 234-kHz interferometric side scan sonar system pole-mounted on the U.S. Geological Survey research vessel San Lorenzo. The system collected full-coverage bathymetric and acoustic backscatter data that were processed to 2-meter spatial resolution (Dartnell and others, 2024). This two-map series displays the results of this mapping. A colored shaded relief bathymetry map (sheet 1) and an acoustic backscatter map (sheet 2) show the lake floor morphology and backscatter intensities, respectively, that can be analyzed for evidence of past earthquakes.

Washington