Geology ReportsSearch

SEARCH · Geology Reports

Results for “JGR Earth Surface”

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.

641 records · Page 22Linked to original sources

Hydrologic response of groundwater and streamflow to natural and anthropogenic drivers of change in headwaters of the upper Colorado River basin during recent wet (1982–1999) and drought (2000–2022) conditions

Study region: Headwaters of the upper Colorado River basin (UCOL), USA Study focus: Surface-water and groundwater numerical models incorporating water-use information were used to investigate changes in climate, water use, and simulated hydrologic responses of snow processes, evapotranspiration, groundwater, and streamflow during recent wet (1982–1999) and drought (2000–2022) periods in the headwater subregions of the upper Colorado River basin. New hydrologic insights for the region: Decreases in average streamflow between wet and drought periods ranged from 20 % in the Colorado River headwaters subregion to 23 % in the Gunnison River headwaters subregion. Like streamflow, average surface runoff was statistically less during the drought than the wet period, with decreases from 24–31 % in the headwaters. On a volume basis, runoff decreases were greater than streamflow decreases in both the Colorado River and Gunnison River headwaters. Although the amount of water-year groundwater discharge to streams remained nearly the same between the wet and drought periods, groundwater as a percentage of streamflow increased between the wet and drought periods, highlighting the importance of groundwater in sustaining streamflow during drought conditions. Multiple linear regression analyses revealed that snowmelt-only models were better than the best precipitation and temperature models at explaining streamflow variability from all headwater subregions for both the wet and drought periods.

Arizona, Colorado, New Mexico, Utah, Wyoming

Continuous monitoring and temporal variability of fluorescence of dissolved organic matter in karst groundwater of the Edwards aquifer, south-central Texas, 2019–24

Widespread urbanization on the Edwards aquifer recharge zone in south-central Texas has prompted the development of new approaches and tools for evaluating the current and future status of water quality in the San Antonio and Barton Springs segments of the Edwards aquifer. The U.S. Geological Survey, in cooperation with the San Antonio Water System and the City of Austin, applied continuous monitoring of fluorescence of dissolved organic matter (fDOM) to characterize the sources and transport of dissolved organic matter and associated constituents in the Edwards aquifer. Continuous fDOM time-series data were adjusted using empirically derived temperature, turbidity, and inner-filter-effect correction factors prior to analysis and interpretation. Fully corrected fDOM time-series data from the study sites were evaluated for short- and long-term temporal variability in fDOM in the context of the hydrologic and climatic conditions that occurred during the 5-year study period to understand vulnerability of the aquifer to potential contaminants in recharge. Study results indicate that fDOM is an effective proxy for dissolved organic carbon that, in turn, is indicative of the influx of recent surface water and associated contaminants, and thus, aquifer vulnerability to those contaminants. Although other continuously monitored water-quality parameters provide insights into changes in water quality in response to varying hydrologic conditions, the quantitative relation between fDOM and dissolved organic carbon directly indicates the timing and magnitude of pulses of organic constituents derived from the land surface, whereas other water-quality parameters do not.

Texas

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri River near Kansas City, Missouri, August 8–9, 2023

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near 8 bridge crossings of the Missouri River near Kansas City, Missouri, on August 8–9, 2023. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches that extended about 1,550 to 1,640 feet longitudinally and generally extended laterally across the active channel from bank to bank during low floodflow to nonflood conditions. These surveys provided the channel geometry and hydraulic conditions of the river at the time of the surveys and provided characteristics of scour holes, which may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. The data collected from the surveys may also be useful to the Missouri Department of Transportation as a record of low floodflow conditions in regards to the stability and integrity of the bridges with respect to bridge scour. Bathymetric data were collected around every in-channel pier. Scour holes were at most piers where bathymetry could be obtained, except for those piers on banks or surrounded by riprap. All the bridge sites in this study were surveyed and documented in previous studies. The average difference between the bathymetric surfaces ranged from 0.07 to 4.16 feet higher in 2023 than 2019, which indicates overall deposition between the survey dates, as might be expected based purely on streamflow at the time of the survey. However, the average difference between the bathymetric surfaces ranged from 1.44 feet higher to 1.88 feet lower in 2023 than 2015, which indicates a dynamic equilibrium of scour and deposition overall between those surveys, despite the lower flow conditions in 2023. Similarly, the average difference between the bathymetric surfaces ranged from 3.18 feet higher to 5.19 feet lower in 2023 than 2011, which indicates a relative equilibrium between scour and deposition overall, albeit the trend was toward scour as might be expected because of the substantial flood event in 2011. Riprap blankets and alignment to flow had a substantial effect on the size of the scour hole for a given pier. Piers that were partially or fully surrounded by riprap blankets had scour holes that were substantially smaller (to nonexistent) compared to piers with no rock or riprap and effectively mitigated the scour holes historically observed at these piers. Several of the structures had piers that were skewed to primary approach flow. At most of the structures, the scour hole was deeper and longer on the side of the pier with impinging flow than the leeward side, with some amount of deposition on the leeward side, as typically observed at piers skewed to approach flow.

Kansas, Missouri

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

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

Kansas, Missouri

Applying U-Th disequilbrium for dating siliceous sinters

Continental hydrothermal systems are critical avenues for the crustal transport of heat and mass captured for geothermal energy and mineral exploration. Thus, understanding their temporal evolution and longevity is important for resource characterization. Deposits of microlaminated siliceous sinter, common surface expressions of high temperature reservoirs (> 170 °C), have the potential to trace hydrothermal histories. Geothermal reservoirs are often located within uranium-bearing silicic volcanic rock where subsurface fluid-rock interactions extract U into hydrothermal fluids. U incorporated in the surface sinter deposit has the potential to provide a U— Th disequilibrium dating option. We focus on samples from El Tatio geyser field in the Altiplano of northern Chile, the largest geothermal system in the Andes. Our resulting 230 Th/U ages, along with the water and deposit elemental compositions, suggest concentrations of U and Th vary predictably along the sinter apron. While distal facies containing the highest U concentrations (> 50 μg/g) are least affected by detrital Th corrections, they can display suspected open-system behavior. In contrast, more medial facies, where bacterial mats and other porous textures are co mmonly concentrated, have only trace amounts of U (< 0.1 μg/g), which leads to unreliable or geologically improbable dates. Proximal facies tend to date most consistently. By comparing existing 14 C ages with 230 Th/U results, 230 Th/U ages tend to be younger than the 14 C ages, supporting the presence of a 14 C-dead carbon influence. New data confirm that the onset of geothermal activity at El Tatio goes back to the late Pleistocene.

Altiplano, El Tatio geyser field

Metal-rich lacustrine sediments from legacy mining perpetuate copper exposure to aquatic-riparian food webs

Historic copper mining left a legacy of metal-rich tailings resulting in ecological impacts along and within Torch Lake, an area of concern in the Keweenaw Peninsula, Michigan, USA. Given the toxicity of copper to invertebrates, this study assessed the influence of this legacy on present day nearshore aquatic and terrestrial ecosystems. We measured the metal (Co, Cu, Ni, Zn, Cd) and metalloid (As) concentrations in sediment, pore water, surface water, larval and adult insects, and two riparian spider taxa collected from Torch Lake and a nearby reference lake. Overall, elevated metal and metalloid concentrations, particularly Cu, were measured in all sediment samples and some surface and pore water samples collected from Torch Lake. For instance, Cu concentrations in the Torch Lake sediment were ∼200% higher than the reference lake and all measured concentrations exceeded predicted effects concentrations by at least ninefold. Within larval insect tissues, we observed 160% higher Cu concentrations than measured in the reference lake, and Cu was the only measured element above predicted effects concentrations in Torch Lake. Adult insects collected at both lakes had similar metal concentrations irrespective of exposure levels. Yet we found 100% higher copper concentrations in Torch Lake riparian spiders, demonstrating elevated exposure risk to insectivores across the aquatic-terrestrial boundary. Our results highlight that other metals in the mixture may not be as concerning to adjacent riparian ecosystems, but copper remains a contaminant of concern in Torch Lake 60 years after mining ceased.

Michigan

Assessing flood water infiltration and storage in a restored floodplain

In urban areas, floodplain restoration is gaining prominence as a strategy for restoring the natural functions of floodplain ecosystems and reducing flood risk. This has spurred research into potential interactions between floodwaters, the hyporheic zone, and the floodplain aquifer. An urban restored stream in Wisconsin, USA, was used as a case study to examine four methods to estimate floodplain infiltration and storage during overbank floods. We characterised flood-related infiltration over a 4-year period from 2018 through 2021 by simultaneously and continuously measuring groundwater levels and vertical temperature profiles with stream water levels linked to high-resolution flood inundation maps. High-resolution topographic data helped to quantify surface floodplain storage and the unsaturated soil volume relative to flood stage. Infiltration estimates from the simple methods align well with those from the more complex methods; however, the complex methods provide additional insights about the factors influencing infiltration. Results from all methods indicate that the volume of water that vertically infiltrates during floods is likely small relative to the total volume of the flood, with 0.08%–0.52% of flood water infiltrating into the floodplain, on average. Spatially variable vertical hydraulic gradients, driven by flood depth, groundwater level, and permeability, imply heterogeneous patterns of infiltration across the floodplain. Gradients favourable for infiltration typically occurred during the onset of flooding but, over the study period, were mostly (98% of the time) favourable for groundwater discharge to the channel (non-flood periods). These findings highlight the importance of considering surface-groundwater dynamics, floodplain soils, and unsaturated floodplain volume in defining the benefits of floodplain infiltration for flood attenuation.

Wisconsin

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

Site response in the Walnut Creek–Concord region of the San Francisco Bay, California: Ground motion amplification in a fault-bounded basin

Thirty‐seven portable accelerometers were deployed in the eastern San Francisco Bay communities of Walnut Creek and Concord to study site response in a fault‐bounded, urban, sedimentary basin. Local earthquakes were recorded for a period of two years from 2017 to 2019 resulting in 101 well‐recorded events. Site response is estimated by two methods: the reference site spectral ratio method and a source‐site spectral inversion method. The reference site spectral ratio method allows investigation of the variability of site amplification with source azimuth and frequency. The source‐site spectral inversion method yields the best least‐squares fit to site response for a database of ground‐motion records. Both methods show substantial amplification in the Walnut Creek–Concord basin below 2 Hz indicating strong surface‐wave development. Greater amplification is seen for sources aligned along the long axis of the basin. Inversion using close‐in sources at short distances yields lower amplification at longer periods than the entire data set due to reduced surface‐wave generation for steeper angles of incidence. Inversion of site response spectra for shallow shear‐wave velocity using a global search algorithm yields V S30 values consistent with generalized mapping results based on geology and topography but with greater variability due to local site variations. 3D finite‐element modeling shows greater amplification in the Walnut Creek–Concord basin with a basin‐edge effect likely contributing to higher ground motions. Topography is also seen to lead to increased scattering and shadowing effects.

California

Beyond the wedge: Impact of tidal streams on salinization of groundwater in a coastal aquifer stressed by pumping and sea-level rise

Saltwater intrusion (SWI) is a well-studied phenomenon that threatens the freshwater supplies of coastal communities around the world. The development and advancement of numerical models has led to improved assessment of the risk of salinization. However, these studies often fail to include the impact of surface waters as potential sources of aquifer salinity and how they may impact SWI. Based on field-collected data, we developed a regional, variable-density groundwater model using SEAWAT for east Dover, Delaware. In this location, major users of groundwater from the surficial aquifer are the City of Dover and irrigation for agriculture. Our model includes salinized marshland and tidal streams, along with irrigation and municipal pumping wells. Model scenarios were run for 100 years and included changes in pumping rates and sea-level rise (SLR). We examined how these drivers of SWI affect the extent and location of salinization in the surficial aquifer by evaluating differences in chloride concentration near surface waters and the subsurface freshwater-saltwater interface. We found the presence of the marsh inverts the typical freshwater-saltwater wedge interface and that the edge of the interface did not migrate farther inland. Additionally, we found that tidal streams are the dominant pathways of SWI at our site with salinization from streams being exacerbated by SLR. Our results also show that spatial distribution of pumping affects both the magnitude and extent of salinization, with an increase in concentrated pumping leading to more intensive salinization than a more widely distributed increase of the same total pumping volume.

Delaware

Land application of drill waste: A scope analysis

Drilling fluid waste land application, a process where drilling wastes are spread and tilled into the land surface, has become common in some petroleum-producing states, however, the potential benefits and risks of this practice are not well studied. Drilling fluids can be water- or oil-based and can have high concentrations of total soluble salts and total petroleum hydrocarbons. Comprehensive chemical characterization of these fluids is not well documented in the literature, and the extent of land application is largely unknown. We hypothesized that the land application of drill waste would fluctuate over time due to economic factors. To begin to understand the extent of historical and potential future land application, we analyzed data from over 5,800 drilling fluid land application permits collected by the Oklahoma Corporation Commission for years 2000, 2005, 2010, and 2015–2020. During the years studied, drilling fluid wastes were applied to more than 250,000 acres in Oklahoma, with over 54,000 thousand barrels (Mbbl) of liquids and nearly 21,000 Mbbl of solids applied. Land application is widespread (occurring in 59/77 counties), however recent drilling activity, land availability, and the economics of transportation have created conditions favorable for land application specifically in the Anadarko Basin. Land application can co-occur with sensitive areas, such as important groundwater and surface-water drinking sources and agricultural fields used for subsistence or feed crop production. Our approach for quantifying the extent of land application, along with further chemical characterization studies, can aid operators and land managers who are considering this practice in assessing the associated benefits and risks.

Oklahoma

Melt generation sources and conditions in the wake of a migrating slab window: Geochemistry and petrology of the million-year history of primitive volcanism at Clear Lake volcanic field, California

Clear Lake volcanic field (CLVF) is the northernmost and youngest (~2.2 Ma to 8 ka) of the volcanic centers distributed along the San Andreas transform fault in western California. The initial phase of CLVF volcanism (interval one) occurred between ~2.2 and 1.3 Ma and extends ~35 km southeast of Clear Lake, forming a semi-continuous upland plateau capped by lava flows, with isolated volcanic remnants on the periphery. This volcanism is broadly characterized by geochemically primitive compositions that reflect three source compositions and conditions of melt generation. (1) Partial melting of upwelling asthenospheric mantle lherzolite at moderate pressures (1.2–1.4 GPa) and temperatures (1297–1329 °C) produced high-CaO (9.8–11.3 wt %) basalts with high Al 2 O 3 (16.8–17.6 wt %), Mg#s (66–70), MgO (8–10 wt %), Ni (103–262 μg/g), and Cr (284–609 μg/g). These high-CaO basalts contain olivine (Fo 87–91 ) phenocrysts with Cr-spinel inclusions ± subordinate plagioclase and crop out only in the southern part of the CLVF. (2) Partial melting of depleted sub-continental lithospheric mantle harzburgite at variable pressures (0.7–1.5 GPa) and temperatures (1097–1299 °C) produced a compositional continuum of med-K 2 O, calc-alkaline, high-MgO basalts through high-MgO andesites with high Mg#s (67–77), MgO (8–14 wt %) and high Ni and Cr abundances (154–439 and 340–1124 μg/g, respectively). Mineral assemblages are olivine (Fo 88–93 ) with Cr-spinel inclusions ± subordinate clinopyroxene, orthopyroxene and plagioclase. Small (<2.5 cm) mantle harzburgite xenoliths and mantle olivine xenocrysts are also found in several of these samples. These high-MgO basalts through andesites represent the largest volume of primitive compositions and have erupted predominantly along the main, fault-controlled northwest-southeast trending axis of volcanism with peripheral outcrops to the north, west, and east. (3) Partial melting of the Gorda eclogite slab edge produced adakitic silicic slab melts with strong depletion in the heavy rare earth elements (Yb = 0.6 μg/g). Subsequent reaction of those melts with depleted ultramafic rocks during ascent imprinted the adakitic dacites with high Mg#s (65–78) and elevated Ni (117–210 μg/g) and Cr (191–283 μg/g). Phenocrysts of orthopyroxene (En 87–94 ) with spinel inclusions (Cr# = 80–88) and extremely Ni-rich (9483 μg/g) olivine cores (Fo 84–93 ) record those reactions. Small-volume outcrops of the adakites on the eastern periphery of the CLVF track the passing slab edge. The trio of melting sources recorded by early CLVF magmatism reflect the tectonically complex environment and the hot (1097–1329 °C), shallow (0.7–1.5 GPa) melting conditions for these primitive compositions and provide estimates of the heat delivered to the crust. Over time, this flux led to maturation of the CLVF magmatic system toward the more voluminous and silicic volcanism that characterizes the balance of its subsequent volcanic history and maintains the present-day anomalously high heat flow in the region. The current interval (interval four) of volcanic activity at CLVF is characterized by low-volume, fault-controlled eruptions of basaltic andesite and andesite suggestive of mantle magma and heat delivery to the crust, similar to interval one. This analogous activity provides motivation for the current study and begs the question of whether the system is undergoing thermal priming for renewed silicic volcanism.

California

Multiple-well monitoring site adjacent to the Midway- Sunset and Buena Vista Oil Fields, Kern County, California

Groundwater quality in and around oil fields in the Southern San Joaquin Valley is of interest to many California residents that rely heavily on groundwater for domestic, commercial, and agricultural use. To help assess the effects of historical oil-field activities and natural geologic sources on groundwater near the southwest margins of the Kern County Groundwater Subbasin, a multiple-well monitoring site was installed near the administrative boundary between the Midway-Sunset and Buena Vista Oil Fields in Kern County, California. The installation of the Midway-Sunset Buena Vista multiple-well monitoring site (MSBV) supports regional analysis of the relations of oil and gas sources to groundwater quality by providing information about the geology, hydrology, geophysical properties, and water quality of the alluvial and upper Tulare aquifers in areas where groundwater data were limited. Data collected from the site included drill cuttings, whole core samples, sidewall core samples, mud-gas analysis, borehole geophysical logs, depth to water measurements, and water quality samples. Whole cores were scanned using dual energy computed tomography. Subsamples of selected cores were analyzed for density, porosity, specific retention, and bulk minerology. Thin sections of the subsamples were prepared, photographed, and examined. Two samples were analyzed using scanning electron microscope technology to examine the microporosity of diatomite laden sediment. Instrumentation installed in the wells collect hourly depth to water measurements. Analysis of the data show there is 355 feet of alluvium overlying the Tulare Formation at the well site. The contact between the two formations is an aquitard resulting in a perched aquifer in the alluvium and unconfined aquifer in the Tulare Formation. The alluvium is more heterogenous and finer grained than the Tulare Formation resulting in markedly higher porosity in the alluvium compared to the Tulare Formation. Higher specific retention observed in the alluvium is attributed to the finer grained sediment and greater abundance of reworked diatomite (as represented by opal-CT [cristobalite-tridymite]) compared to the Tulare Formation. Total dissolved solids (TDS) approached or exceeded 10,000 milligrams per liter (mg/L) in the alluvium from approximately 176 to 242 feet below land surface and at the top of the Amnicola clay at approximately 670 feet below land surface within the Tulare Formation. Elevated TDS, chloride, and boron concentrations in the alluvium and on top of the Amnicola clay likely reflect groundwater that is mixed with oil-field water. Water chemistry and modern-aged groundwater in the alluvial monitoring well (MSBV #3) are consistent with the oil-field water in the alluvium being derived from documented historical surface disposal of oil-field water upslope (northwest) of the site. Water chemistry and pre-modern groundwater age in the deeper Tulare monitoring well (MSBV #1) on top of the Amnicola clay are consistent with oil-field fluids derived from upslope natural geologic sources or old oil wells that leak in the subsurface. Shallow groundwater in the Tulare (MSBV #2) is not affected by mixing with oil-field sources.

EarthArXiv

Simulated mean monthly groundwater-transported nitrogen loads in watersheds on the north shore of Long Island Sound, 1993–2022

Elevated nitrogen loads are pervasive in the Long Island Sound, an estuary that receives freshwater and nutrients from both surface-water and groundwater discharge. Surface-water nitrogen loads to the Long Island Sound are relatively well characterized, but less is known about groundwater-transported nitrogen loads. Prior work on the northern shore of Long Island Sound (Connecticut and areas of New York and Rhode Island) suggested that groundwater travel times are relatively short (median less than 2 years) and that decade-long nutrient legacies are not widespread. Because the travel times are short, groundwater flow and nutrient loads likely vary substantially between months. In the current study, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency’s Long Island Sound Study and the Connecticut Department of Energy and Environmental Protection, developed a set of models to better characterize spatial and temporal patterns of groundwater-transported nitrogen loading from atmospheric deposition, septic systems, and fertilizers within the study area. The models provide an estimate, with uncertainty, of groundwater-transported nitrogen loads in the study area, filling a key gap in the nitrogen budget for Long Island Sound. The models also highlight the spatial and temporal variation in nitrogen loading throughout the study area. The modeling workflow involved four models. (1) A soil-water-balance model was developed by using the Soil-Water-Balance software to simulate groundwater recharge across the study area for water years 2005 through 2022. The simulated mean monthly recharge from the soil-water-balance model was used as input into a groundwater-flow model. (2) The groundwater-flow model was developed by using the MODFLOW 6 software and data for water years 1993 through 2022 and simulates average monthly hydrologic conditions. The groundwater-flow model was calibrated by using the Iterative Ensemble Smoother method within the PEST++ software. The Iterative Ensemble Smoother method generates an ensemble of sets of parameter values, with each set producing reasonable simulated hydrologic parameter values. (3) An ensemble of MODPATH particle-tracking simulations were run to generate particle flow paths and travel times, with each simulation using a different set of the flow model parameters. (4) A nitrogen load model uses the MODPATH simulation outputs to track nitrogen from the land surface through multiple attenuation zones until it discharges into fresh or saline surface water. As with the groundwater-flow model, the nitrogen model simulated average monthly groundwater-transported nitrogen loads for water years 1993 through 2022. One novel aspect of the nitrogen load model is that the nitrogen attenuation parameters were calibrated to observed nitrogen loads. Across the ensemble of simulated nitrogen loads, the median study-area-wide monthly simulated nitrogen loads from the aquifer to Long Island Sound throughout the year ranged from 900 to 18,600 kilograms of nitrogen per day, with a median load of 5,100 kilograms of nitrogen per day. The simulated loads were based on average monthly conditions for water years 1993 through 2022. Loads were highest during the winter and early spring and lowest during the late summer. However, simulated travel times for groundwater and nitrogen loads discharged to Long Island Sound during summer were longer than travel times for groundwater and loads discharged during the winter, indicating that, on average, groundwater discharged during summer traveled along different, and longer, flow paths, than groundwater discharged during winter. This indicates that summer loads would respond more slowly to changes in nitrogen inputs at the water table than winter loads. Over the entire study area, approximately 15 percent of the simulated load is from atmospheric deposition sources, 30 to 40 percent is from fertilizer, and 50 to 60 percent is from septic systems. The final analysis of the study involved simulating the change in groundwater-transported nitrogen load in response to upgrading septic systems or reducing fertilizing inputs to areas of turf grass. Both management interventions reduced the groundwater-transported nitrogen load, and reductions were greater in areas with greater loads from septic systems or turf-grass fertilizers. The delay between management actions and substantial reductions in groundwater-transported nitrogen loads varied seasonally; loads during the late summer months remained elevated longer than the winter loads.

Connecticut, Rhode Island

Baseflow and snowmelt sustained streamflow in the Upper Colorado River Basin, 1986-2020

The Upper Colorado River Basin (UCRB) faces substantial water availability limitations. Although most streamflow originates as snowmelt, the partitioning of snowmelt between surface runoff and groundwater recharge and subsequent groundwater discharge to streams is highly uncertain. On average, over half of the streamflow in the UCRB is estimated to originate from groundwater discharge to streams, highlighting the importance of baseflow in sustaining surface water. However, the historical patterns of baseflow and streamflow, along with their variability over space and time and their specific sources, remain unknown at the basin scale. This study addresses those gaps by characterizing the sources and transport pathways of both baseflow and streamflow in the UCRB at a seasonal timestep from 1986 to 2020, including the lagged delivery of subsurface water to streams beyond the current season, using coupled models of baseflow and streamflow. Between 1986 and 2020, on average 63% of UCRB streamflow originated from baseflow. About half of this baseflow took longer than one season to reach streams, and outside the snowmelt season, baseflow was the dominant source of streamflow. Snowmelt was a key source of both baseflow and streamflow. Current season snowmelt contributed 33% of streamflow via runoff, and 22% of the 29% of streamflow that originated as current season baseflow via subsurface flow to streams. Over the study period, baseflow index (BFI) declined in headwaters and increased at mid-elevations. Springtime increases in BFI demonstrate the increasingly important role baseflow plays in water supply. Identifying the sources, locations, and timing of water that contributed to the UCRB outlet can inform management of water resources in the basin.

Arizona, Colorado, New Mexico, Utah, Wyoming

Base flow and ground water in upper Sweetwater Valley, Tennessee

The upper Sweetwater Valley area has a flow system with complex interaction between surface and ground water. A water budget study indicated that during dry years approximately three-fourths of the annual flow to Sweetwater Creek may be derived from ground-water sources. Hydrograph analysis showed seasonal variation of recharge to the ground-water flow system. Streamflow records were analyzed to estimate the frequency of low flow of Sweetwater Creek at river mile 16.7, and indicated the lowest average flow for 1 day in 20 years to be about 5.1 cubic feet per second. Two periods of base-flow measurements of Sweetwater Creek identified channel reaches with significant gains and losses of streamflow. Base flow measurements also showed interbasin transfer of water among sub-basins of the valley. Major flow surpluses were associated with areas in which the majority of flow originated at a spring. Topographically low areas adjacent to the main stem of Sweetwater Creek generally have surplus flow. Topographically higher areas generally have deficient surface outflow unless significant spring flow occurs in the basin. Ground-water recharge occurs by water draining into sinkholes, faults, and fractures. Ground-water flow is regionally diffused across formation strikes from the topographically low areas unless the water is exposed to highly permeable formations or impervious formations. Ground water infiltrates the highly permeable formations and flows along strike. Ground water encountering impervious formations may discharge at small springs at the contact, or may reroute along the contact if the upgradient rock is sufficiently permeable or has well developed secondary porosity. Groundwater discharges to streams at innumerable seeps and at a few large springs. Areas of ground-water flow up-gradient of large springs are hypothesized as likely areas of significant ground-water reservoirs.

Tennessee

Dynamic baseflow storage estimates and the role of topography, geology and evapotranspiration on streamflow recession characteristics in the Neversink Reservoir Watershed, New York

Estimates of dynamic groundwater volumes supplying baseflow to streams are important for water availability projections during extended periods of drought. The primary goals of this study were to provide dynamic storage volume estimates, inferred from streamflow recession analysis, for baseflow regimes within seven gaged catchments within the Neversink Reservoir Watershed (NRW), a critical municipal water source for New York City. Additionally, geomorphological properties, surficial geology and hydro-meteorological processes were quantified and described in relation to time and spatially variable recession behaviour and storage estimates across the NRW. To explore these relationships, we (1) evaluated seasonal trends in streamflow recession behaviour in relation to modelled potential evapotranspiration (PET) and catchment runoff rates, (2) derived empirical streamflow models for cool-season runoff using both linear and nonlinear reservoir assumptions for baseflow and (3) calculated metrics related to the geology and geomorphology of each catchment and compared these metrics to area normalised baseflow dynamic storage estimates. Results show that baseflow recession behaves as a nonlinear reservoir, and applying linear groundwater reservoir assumptions may underestimate the total dynamic storage volumes compared to what would be predicted for a nonlinear reservoir. Increases in PET caused decreases in storage conditions that resulted in increased recession rates and nonlinearity in streamflow recession during the growing season. Additionally, we found that while no single physical catchment characteristic solely predicted catchment storage dynamics, sediment volume and stream gradients were stronger predictors of normalised storage volumes than catchment surface area or surface topography alone. Within the NRW, catchments with the highest sediment volume exhibited the lowest recession rates and higher dynamic storage volumes, while the smallest catchment, mostly devoid of sediment, had the fastest recession rate and lowest dynamic storage volume.

New York

Water table rise sustains carbon release from soils in wetland-dominated landscapes: An intact soil core study

Wetland-dominated landscapes influence carbon cycling through their potential to act as both carbon sinks and sources. Wetlands in low-relief landscapes have dynamic terrestrial-aquatic interfaces that change seasonally with variable surface water and groundwater levels. However, few studies have directly quantified dissolved organic matter (DOM) release and greenhouse gas (CO 2 , CH 4 ) fluxes from wetland soils along terrestrial-aquatic interfaces as they are seasonally re-saturated by groundwater. To estimate groundwater-mediated soil DOM and gas fluxes, we performed laboratory simulations of vertical groundwater rise on intact soil cores collected from four Delmarva bay wetlands located in the Mid-Atlantic United States. At each wetland, one core was collected from within the wetland basin and the other from the transitional zone at the basin edge. Cores were re-saturated with groundwater over 15 days and then kept fully saturated for an additional 25 days. Source groundwater, soil porewater, and exfiltrated surface water samples were collected and analyzed for pH, ORP, DOM concentration, and DOM optical indices. In both the wetland and transition zone cores, porewater DOM concentrations increased over the wet up and were sustained during prolonged saturation. Optical indices shifted from recently produced, microbial-like signatures towards aromatic, terrestrial-like signatures. Fluxes of CO 2 decreased as the duration of soil saturation increased and soil cores switched from CH 4 sinks to sources upon full soil core saturation. Results indicate that groundwater rise sustains carbon mobilization from soils in wetland-dominated landscapes, emphasizing the need to understand how climate-driven changes to groundwater dynamics may affect carbon fluxes along terrestrial-aquatic interfaces.

Delaware, Maryland