SEARCH · USGS Science

Results for “Geological Society of America Bulletin”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,666 records · Page 64Linked to original sources

USGS2024

Utility of an instantaneous salt dilution method for measuring streamflow in headwater streams

Streamflow records are biased toward large streams and rivers, yet small headwater streams are often the focus of ecological research in response to climate change. Conventional flow measurement instruments such as acoustic Doppler velocimeters (ADVs) do not perform well during low-flow conditions in small streams, truncating the development of rating curves during critical baseflow conditions dominated by groundwater inflow. We revisited an instantaneous solute tracer injection method as an alternative to ADVs based on paired measurements to compare their precision, efficiency, and feasibility within headwater streams across a range of flow conditions. We show that the precision of discharge measurements using salt dilution by slug injection and ADV methods were comparable overall, but salt dilution was more precise during the lowest flows and required less time to implement. Often, headwater streams were at or below the depth threshold where ADV measurements could even be attempted and transects were complicated by coarse bed material and cobbles. We discuss the methodological benefits and limitations of salt dilution by slug injection and conclude that the method could facilitate a proliferation of streamflow observation across headwater stream networks that are highly undersampled compared to larger streams.

Virginia↗

USGS1998

Water resources of Mellette and Todd counties, South Dakota

Mellette and Todd Counties are located in south-central South Dakota and have a combined area of 2,694 square miles. The White River and its tributaries, which include the Little White River, drain Mellette County and about one-half of Todd County. Tributaries to the Niobrara River, which include the Keya Paha River, drain the other one-half of Todd County. The average discharge of the Little White River is about 56 cubic feet per second as the river enters Todd County and is about 131 cubic feet per second as it discharges to the White River in northern Mellette County. The average discharge of the Keya Paha River just outside Todd County is about 39 cubic feet per second. The average annual runoff for Mellette and Todd Counties ranges from 0.94 to 2.36 inches based on records from nine streamflow-gaging stations in and near the counties. The average annual runoff is 1.62 inches, which compares with the average annual precipitation of about 19 inches. In Todd County, shallow wells completed in the alluvial, Ogallala, Arikaree, and White River aquifers generally can supply water that has low concentrations of dissolved solids, is fresh, and is soft to moderately hard. Ground water from shallow aquifers is limited in Mellette County; therefore, deep wells, often greater than 1,000 feet, are sometimes installed. The Pierre Shale often is used to supply rural domestic and stock wells in Mellette County even though well yields are low and the water has high dissolved solids, is moderately saline, and is very hard. Alluvial aquifers are present in both counties and store an estimated 1.6 million acre-feet of water. The water quality of the alluvial aquifers is dependent on the underlying deposits, and generally the water has low concentrations of dissolved solids, is fresh, and is soft to moderately hard where underlain by the Ogallala and Arikaree Formations; has moderate concentrations of dissolved solids, is slightly saline, and is hard where underlain by the White River Group; and has high concentrations of dissolved solids, is saline, and is very hard where underlain by the Pierre Shale. Also, yields often are lower where the alluvial aquifers are underlain by the Pierre Shale. The Ogallala aquifer is present in only Todd County, and the Arikaree aquifer is present throughout most of Todd County and southwestern and south-central Mellette County. The Ogallala aquifer contains an estimated 17 million acre-feet of water in storage, and the Arikaree aquifer contains an estimated 50 million acre-feet of water in storage. Both aquifers generally are suitable for irrigation, and yields from these aquifers are sometimes greater than 1,000 gallons per minute. Nitrate concentrations in 13 out of 92 water samples collected from the Ogallala aquifer exceeded the Primary Drinking Water Maximum Contaminant Level (MCL) of 10 milligrams per liter. In 11 out of 46 samples collected from the Arikaree aquifer, arsenic concentrations exceeded the MCL of 50 micrograms per liter. The White River aquifer, where present, is usually the shallowest source of ground water in Mellette County. The White River aquifer also is used in northern Todd County where the Ogallala and Arikaree aquifers are not present. The White River aquifer contains an estimated 50 million acre-feet of water in storage. Reported yields from the aquifer range from 1 to 30 gallons per minute, which generally is insufficient to support irrigation in most areas. However, yields are sufficient for livestock-watering and rural-domestic purposes. In both counties, the Pierre Shale is the shallowest bedrock aquifer and is exposed at the land surface throughout most of Mellette County. This aquifer is used primarily in Mellette County. Although the aquifer contains an estimated maximum of 1.5 million acre-feet of water in storage, it is not a viable source of ground water because the aquifer is relatively impermeable, yields are low, and water usually can be obtained from shallower sources, especially in Todd County. Reported yields from the Pierre Shale aquifer range from 1 to 8 gallons per minute. Because few test holes and wells penetrate below the Pierre Shale, little is known about the extent of the deeper bedrock aquifers. All wells completed in the Dakota Sandstone, Inyan Kara, and Minnelusa and Madison aquifers in the counties are used for stock-watering purposes. High concentrations of dissolved solids and hard water are characteristic of the water quality in the bedrock aquifers. Depths to the top of the deeper bedrock aquifers range from 1,270 feet to greater than 2,000 feet below land surface.

South Dakota↗

USGS2025

Airborne radiometric data map alteration of porphyry copper systems in the Elkhorn district, MT

The Elkhorn district in southwestern Montana is known for past production of high-grade base and precious metals from skarn and carbonate replacement deposits. The district also hosts a porphyry Cu-Mo deposit. A radiometric survey flown over the district and surrounds provides maps of surficial concentrations of potassium, thorium, and uranium. Handheld gamma ray spectrometer measurements on altered and unaltered rocks in the district and surrounding region provide a means to ground truth the airborne data. We demonstrate that the airborne and ground-based measurements show similar map patterns. We calculate ratios of the radioelement concentrations and present a potassium enhancement map that combines potassium and ratios of potassium to thorium and uranium. The results highlight the rocks previously mapped as having K-feldspar-biotite and quartz-sericite alteration assemblages in the Elkhorn district and map additional hydrothermal systems in the region, including known alteration in the Radersburg district. The data enhancement techniques can be used as a screening tool for mapping additional porphyry copper systems.

Montana↗

USGS2026

Multi-objective optimization of a hydro-economic model in an over-allocated agricultural basin

Groundwater depletion for agricultural irrigation poses significant environmental and economic challenges. This study introduces a proof-of-concept that combines hydro-economic modeling, scenario-based modeling, and multi-objective optimization to manage pumping curtailment in an over-allocated basin in the western United States. Three optimization scenarios were evaluated, each offering different degrees of management flexibility. Results reveal that scenarios with finer spatial resolution achieved greater environmental benefits per unit profit loss. Additionally, strategies allowing fractional reductions in curtailed wells–rather than complete shutdowns based on water rights seniority–substantially improved efficiency, highlighting the value of increased decision-making flexibility. Although scenario testing can aid stakeholder engagement and strategy exploration, multi-objective optimization provides a systematic framework to quantify tradeoffs between competing objectives. This combined approach demonstrates promise for building consensus and supporting the design of sustainable water management strategies that balance agricultural livelihoods with ecosystem preservation.

Oregon↗

USGS2025

A framework for understanding the effects of subsurface agricultural drainage on downstream flows

Understanding controls on streamflow volume and magnitude is important to water resource management applications, such as critical water and transportation structure design and floodplain mapping. Changes in land use and agricultural practices, such as subsurface agricultural drainage, may be contributing to changes in streamflow characteristics. Subsurface agricultural drainage, also known as tile drainage, is the practice of installing drains in the subsurface of agricultural fields to improve productivity. Because of the complex interactions between subsurface drainage systems, precipitation, local soil conditions, and land management practices, it is difficult to determine how subsurface agricultural drainage affects downstream flow. Previously developed subsurface agricultural drainage conceptual models under dry, saturated, and winter conditions are summarized, and current literature on the effects of subsurface agricultural drainage on downstream flows, focusing on peak flow, non-event flow, and total flow to develop frameworks for discussing these systems is compiled. The effects that subsurface drainage has on hydrologic systems are expected to vary by site and are seasonally based on system design, soil type, moisture conditions, precipitation characteristics, and land conditions. Subsurface drainage can affect the magnitude of peak flow by converting surface runoff from a storm event to subsurface runoff. By increasing hydrologic connectivity of a catchment, subsurface drainage can increase non-event flow or the flow between two storm events, typically dependent on lateral flow through the subsurface and groundwater. Theoretically, by diverting water from groundwater recharge or by reducing water available for evapotranspiration, subsurface drainage may increase the total volume of flow. Precipitation changes may increase infiltration, excess overland flow, and flood risk regardless of the presence or absence of subsurface drainage.

Illinois, Iowa, Michigan, Minnesota, Missouri, Mon↗

USGS2024

Testing rate‐and‐state predictions of aftershock decay with distance

We analyze aftershocks of the 2019 M 7.1 Ridgecrest mainshock and isolated M 5–6 mainshocks in southern California to test predictions made by the rate‐and‐state friction model of Dieterich (1994) . Rate‐and‐state friction predicts that the seismicity rate after a stress step follows Omori decay, where the Omori c ‐value, which is the saturation in aftershock rate observed at small times, is larger for smaller stress steps. Put in the context of an aftershock sequence, this predicts that the Omori c ‐value will be systematically larger at greater distances from the mainshock. To our knowledge, this predicted effect has not been observed. In part this may be because the Omori c ‐value is difficult to measure because it often reflects short‐term catalog incompleteness rather than a true saturation in aftershock rate. We explore the dependence of the Omori c ‐value on the distance to the mainshock by applying the “ a ‐positive” method ( van der Elst and Page, 2023 ). This method is insensitive to short‐term aftershock incompleteness and allows resolution of the true aftershock rate deep into the mainshock coda. For aftershocks of the Ridgecrest mainshock and stacked M 5–6 mainshocks, we observe systematic differences in early aftershock rates, relative to mainshock distance, consistent with the predictions of rate‐and‐state friction. Furthermore, for the larger Ridgecrest dataset, we observe that aftershocks nearer to the mainshock start earlier, and we resolve a flattening of the Omori curve consistent with a larger Omori c ‐value for the farthest aftershocks, as predicted by Dieterich (1994) .

Seismological Research Letters↗

USGS2025

Critical mineral inventory of select IOA-IOCG deposits, southwestern USA

Critical minerals are necessary for modern technology and strategic purposes. Their increasing importance requires finding new and nontraditional resources. Samples of ore, altered, and unaltered host rock were collected from 26 iron mines and prospects in California, Nevada, and Utah to assess the potential of these deposits to host economic quantities of different critical minerals. Geochemical analyses were conducted by 61 element ICP-OES-MS sodium peroxide fusion and major elements determined by WDXRF. These deposits concentrated many critical minerals beyond what is found in average upper crustal abundances, such as Sb, As, Bi, Co, Ga, Mg, Mn, Ni, Nb, Pd, REE, Sc, Te, Sn, Ti, W, V, and Zn. However, most of these are not concentrated enough in the ore to be considered as economic resources. Those critical minerals that are enriched enough in some of these deposits to possibly be considered as by-product commodities are Ni, REE, V, and potentially Co and Ga. These enrichments were not uniform, with REE more likely to be enriched in IOA deposits, whereas Co, Ga, Ni, and V could be found enriched in either IOA or IOCG deposits.

California, Nevada, Utah↗

USGS2026

A case study of maximum depositional ages from terrestrial sandstones near the Cretaceous—Paleogene transition, western Williston Basin, USA

We present a new Bayesian method for deriving maximum depositional ages from detrital K-feldspar using total-fusion 40 Ar/ 39 Ar data. Individual analyses weighted by K/Ca ratio, age uncertainty, and percent radiogenic argon prioritize a result to come from accurately measured volcanic sources. Results from this method applied to sandstones from the Hell Creek region of northeastern Montana show that detrital K-feldspar maximum depositional ages align closely with detrital U-Pb zircon ages as well as tephra based chronostratigraphic constraints from both 40 Ar/ 39 Ar and U-Pb data. An age model informed by maximum depositional ages and available tephra data provides an estimate of 1.74 ± 1.04 Ma for the duration of deposition for the Hell Creek Formation. Combined age model and maximum depositional age data suggest ≤ 1.84 Ma of missing time is represented by the Hell Creek and Fox Hills formational disconformity, coinciding with the final regression of the Western Interior seaway.

Montana↗

USGS2025

The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought

Quantitative understanding is lacking on how the depth of active groundwater circulation in bedrock affects mountain streamflow response to a multi‐year drought. We use an integrated hydrological model to explore the sensitivity of a variety of streamflow metrics to bedrock circulation depth and porosity under a plausible extreme drought scenario lasting up to 5 years. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. With drought, a deeper circulation system with higher drainable porosity more effectively buffers minimum flow and significantly limits perennial stream loss in comparison to a shallow circulation system. Streamflow buffering is accomplished through extensive groundwater storage loss. However, deeper circulation systems experience prolonged recovery from drought in comparison to storage‐limited shallow systems. Research highlights the importance of characterizing the deeper bedrock hydrogeology in mountainous watersheds to better understand and predict drought impacts on stream ecosystem health and water resource sustainability.

Colorado↗

USGS2025

Dietary bioavailability of uranium to a model freshwater invertebrate

Uranium (U) mining increases environmental exposures. Understanding how U is taken up by organisms can aid in evaluating the potential for bioaccumulation and toxicity. Although the importance of aqueous geochemical speciation is well recognized for U bioavailability after dissolved exposures, far less is known about the processes controlling U bioavailability after dietary exposures. This study characterizes the biogeochemical drivers of dietary U uptake in the freshwater snail Lymnaea stagnalis in laboratory experiments. Solids tested included benthic diatoms pre-exposed to dissolved U(VI), soils from contaminated U mine sites, and colloidal hydrous ferric oxide (HFO) synthesized in the presence of dissolved U(VI) or with U complexed by natural organic matter (NOM). Results showed that U was bioavailable from all solids. Uranium assimilation efficiency (AE), a proxy for dietary U bioavailability, varied among solids. AE was lowest for the U-contaminated soils (25 ± 17%) and highest for the U-laden diatoms (71 ± 13%). AE varied slightly among HFO preparations, suggesting modest influences of NOM and iron on U bioavailability. Increases in dietary U exposures reduced feeding rates, and the extent of feeding inhibition appeared inversely related to U bioavailability. The high U assimilation and range of bioavailability have implications for toxicity risks inferred without considering dietary uptake.

Environmental Science and Technology↗

USGS2025

Reaction kinetics and accelerant effects of sulfides in early mature hydrocarbon generation using hydrous pyrolysis

Hydrocarbon generation in organic-rich sediments is influenced by the molecular organic composition and relative abundance of associated minerals. Certain mineral-derived elements act as catalysts and reaction intermediaries, facilitating early-stage hydrocarbon formation in potential source rocks. This study investigated the role of sulfur contributed from pyrite as an accelerant in thermal reaction, focusing on its effects on early maturation and consequent hydrocarbon generation from gilsonite (low-sulfur solid petroleum). Hydrous pyrolysis (HP) experiments were conducted on mixtures of gilsonite and pyrite in varying ratios (1:0.1, 1:0.5, 1:1, 1:2, and 1:10 w/w gilsonite:pyrite) at 320, 350, and 370 °C for 72 h. Untreated and thermally altered residues were analyzed using solid bitumen reflectance (BR o , %), total organic carbon (TOC) content, programmed temperature pyrolysis, scanning electron microscopy with energy-dispersive spectroscopy (SEM-EDS), and X-ray diffraction (XRD) to evaluate the potential accelerant role of pyritic sulfur in hydrocarbon formation. The results show HP residues at 320 and 350 °C with greater pyrite concentrations had higher BR o , while reflectance values were similar in the 370 °C residues, regardless of pyrite concentration, suggesting enhanced reaction at lower thermal conditions. Increasing pyrite content systematically decreased hydrogen index (HI) values while increasing the transformation ratio (TR) and production index (PI), indicating enhanced conversion of organic matter to hydrocarbons with increasing pyrite concentrations. Gas yields increased with pyrite addition, particularly at 350 °C, confirming secondary cracking effects. However, gas production stabilized or declined at higher pyrite loadings (1:10), suggesting alternative reaction pathways such as coke formation. Our data indicate the presence of pyrite lowers the activation energy for thermal cracking, shifting peak experimental hydrocarbon generation temperatures downward by 20–30 °C, with the most pronounced accelerant effects observed at moderate pyrite concentrations (1:0.5 and 1:1). The thermodynamic framework reveals that pyrite stability is influenced by experimental conditions, with pyrrhotite formation favored in the presence of gilsonite due to reduced oxygen fugacity. Pyrite transformation to pyrrhotite, as observed through XRD, SEM-EDS, and predicted by thermodynamic data, further supports the accelerant role of S, as pyrrhotite exhibits a higher hydrogen transfer potential, promoting early oil generation. These findings highlight the importance of pyrite in modulating hydrocarbon generation pathways in organic-rich systems.

Journal of Analytical and Applied Pyrolysis↗

USGS2025

Time-varying rates of organic and inorganic mass accumulation in southeast Louisiana marshes: Relationships to sea-level anomalies and tropical storms

Louisiana's coastal wetlands are complex systems that require a continuous input of organic and inorganic material to keep pace with relative sea-level rise. Coastal restoration projects such as sediment diversions are being implemented to mitigate land loss and increase availability of inorganic sediment to coastal wetlands, and marshes specifically rely on organic material to build soil volume and maintain surface elevation. Interannual-to-decadal sea-level anomalies such as hurricanes can affect marsh accretion, mineral deposition, and plant productivity. In this light, complex ecogeomorphic feedback controls whether a marsh surface is sustainable or eroded/drowns. This study performs some of the first differential vertical accretion rates (VARs) and organic and inorganic mass accumulation rates (MARs) over time in SE Louisiana marshes determined from the 210 Pb Constant Rate of Supply model, coupled with standard 137 Cs VARs. These accumulation rates over the past ∼100 years were measured from a total of six brackish and salt marsh locations in Barataria Basin near the proposed Mid-Barataria Sediment Diversion. They were then related to interannual sea-surface elevations at Grand Isle, Louisiana, over the last ∼60 years and recorded hurricane activity in the delta. Results show VARs range from 0.63 cm/y to 1.69 cm/y and total MARs range from 0.11 to 0.43 g/cm 2 /y. Temporally, VARs and MARs (total, inorganic, and organic) are characterized by gradual increases in rates with decreasing age along with episodic peaks in VARs and MARs. The findings of this study indicate that no relationship occurs between sea-level anomalies and VARs or organic and inorganic MARs; however, a strong relationship appears to occur between major hurricanes to VAR and MAR contributions. Furthermore, high water content (81 ± 8%) and organic-rich soils in the sediment cores highlight the significance of belowground biomass and associated pore volume in maintaining marsh elevation in the study area.

Louisiana↗

USGS2026

Landslide volume estimation from seismic waveform features using a global catalog of seismogenic events

Seismic waves generated by large, rapid landslides encode information about the source and can be analyzed rapidly following an event. Even remote landslides can pose hazards to downstream communities, so rapid detection and characterization using existing seismic monitoring networks could be beneficial. In this study, we expand on past regionally limited work by presenting a globally applicable method for estimating landslide volume from seismic features that could be integrated into future landslide seismic monitoring frameworks. We train the model using multivariable linear regression and five seismic features derived from recordings of 129 landslide events of a range of styles and locations with independently estimated volumes. We present two preferred models, one that combines long‐period (LP) and high‐frequency (HF) features and one for use on smaller landslides without observed long‐period signals. We find that our best‐performing model, applicable to landslides larger than 100,000 m 3 with signals containing observable long‐period energy, requires only two features: LP (20–100 s) absolute maximum amplitude and HF (1–5 Hz) rise time (time between the signal start and the maximum envelope amplitude) and has an R 2 score of 0.79. This model predicted volumes within one order of magnitude for 55 out of 58 events. We find that this combination mitigates the trade‐off between mass and acceleration and other variations of landslide style that limit methods based on amplitude alone.

Seismological Research Letters↗

USGS2025

From critical minerals to food security, the benefits of data collaboration

The volume of data in the public geoscience sphere is rapidly and continually expanding. At Geoscience Australia (GA) we saw an over 500% increase in data points within our relational databases between 2018 and 2024, over the life of the Exploring for the Future (EFTF) program. With the Resourcing Australia’s Prosperity initiative, a continued increase in data quantity will be seen for the next 10 to 35 years. At the same time, a broadening audience for geoscience data is increasing the desire to enhance the diversity of delivery streams. This ranges from data-dense highly technical outputs for geoscience specialists to curated interpretive products for people who are non-geoscientists. Development of these curated outputs has contributed to our awareness of the need for data to be collected and compiled in a way that ensures its reuse, with a focus on quality metadata and data provenance.

Conference Paper↗

USGS2025

Chemostratigraphy of the Cretaceous Hue Shale in Arctic Alaska: Exploring paleoceanographic controls on trace element enrichment, organic matter accumulation, and source-rock evolution

We document chemostratigraphy in an outcrop of late Albian to early Campanian (∼103–82 Ma) marine source rocks to better understand paleoenvironmental controls on trace element (TE) enrichment and organic matter accumulation in the distal Colville foreland basin of Arctic Alaska and how those drivers are linked to arc volcanism and successions of Cretaceous oceanographic and climatic biogeochemical events. This unique, 113-m-thick section of Cretaceous Hue Shale deposited during a series of previously undocumented Arctic Cretaceous oceanic anoxic events (Lease et al., 2024) is the only known exposure of thermally immature (0.48–0.52% R o , random vitrinite reflectance) Hue Shale in Arctic Alaska. Strata comprise mainly clay-rich mudstone with elevated total organic carbon (TOC) and hydrogen index values reaching 26.3 wt% (mean = 7.5 wt%) and 689 mg hydrocarbon (HC)/g TOC (mean = 385 mg HC/g TOC), respectively. Maceral composition consists predominantly of fluorescent amorphous organic matter, with abundant brightly fluorescent alginite, including Tasmanites , acritarchs, and Leiosphaeridia . Discrete layers of volcanic ash (preserved as bentonite) are present throughout the section and provide quantitative age control based on U–Pb dates. Chemostratigraphic trends are investigated to advance our understanding of local oceanographic conditions and controls on stratigraphic and temporal heterogeneity of Brookian source rocks. Concurrent sedimentary enrichment in Mo, U, V, Pb, and Cu across the Albian–Cenomanian boundary of the exposed basal gamma-ray zone, may reflect anoxic to euxinic benthic redox conditions favoring organic matter accumulation and preservation. Fluctuating degrees of anoxia-euxinia are inferred throughout the overlying Hue Shale succession, reflected by varying patterns of TE enrichment and TE–TOC covariation. Elevated C org /P molar ratios (>250) across most of the section, with several values exceeding 690, signify that enhanced biological productivity is sustained throughout deposition. Enhanced productivity, recorded by both C org /P and excess Ba, also parallels increases in source rock richness (elevated TOC and S 2 values) during the late Albian–early Cenomanian and late Cenomanian – Turonian. Enhanced productivity and variations in oceanic circulation/stratification likely both drove changes in benthic redox conditions that favored organic carbon accumulation and preservation. Increased continental arc volcanism (e.g., Okhotsk–Chukotka volcanic belt) and High Arctic Large Igneous Province magmatic eruptions throughout the Cretaceous, inferred to have influenced nutrient cycling and local aqueous nutrient availability, also have been invoked as potential drivers of organic carbon burial and source-rock development across the sedimentary sequence. Results presented here document the organic-rich and oil-prone source-rock quality of the Hue Shale in the distal part of the Colville foreland basin and bolster the potential for a Cretaceous petroleum system beneath the eastern North Slope.

Alaska↗

USGS2025

Impact of land subsidence on housing sale values: Evidence from the San Joaquin Valley, California

This study assesses the impact of land subsidence on housing sale values in the San Joaquin Valley, California. The study utilizes home sale transactions and vertical land-surface displacement data from Interferometric Synthetic Aperture Radar techniques. Using fine-scale fixed effects, matching, as well as a repeat-sales approach, our results indicate that land subsidence resulted in a 2.4% to 5.8% reduction in housing sale values, with the largest reductions occurring in areas where substantial subsidence occurred. Such findings may have implications for groundwater management and can potentially help inform policy design to help mitigate the causes and impacts of land subsidence.

California↗

USGS2025

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

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

Colorado↗

USGS2024

Agricultural tile drains increase the susceptibility of streams to longer and more intense streamflow droughts

Streamflow droughts are receiving increased attention worldwide due to their impact on the environment and economy. One region of concern is the Midwestern United States, whose agricultural productivity depends on subsurface pipes known as tile drains to improve trafficability and soil conditions for crop growth. Tile drains accomplish this by rapidly transporting surplus soil moisture and shallow groundwater from fields, resulting in reduced watershed storage. However, no work has previously examined the connection between tile drainage and streamflow drought. Here, we pose the question: does the extent of watershed-level tile drainage lead to an increased susceptibly and magnitude of streamflow droughts? To answer this, we use daily streamflow data for 122 watersheds throughout the Midwestern United States to quantify streamflow drought duration, frequency, and intensity. Using spatial multiple regression models, we find that agricultural tile drainage generates statistically significant ( p < 0.05) increases in streamflow drought duration and intensity while significantly reducing drought frequency. The magnitude of the effect of tile drainage on streamflow drought characteristics is similar to that of water table depth and precipitation seasonality, both of which are known to influence streamflow droughts. Furthermore, projected changes in regional precipitation characteristics will likely drive the installation of additional tile drainage. We find that for each 10% increase in tile-drained watershed area, streamflow drought duration and intensity increase by 0.03 d and 12%, respectively, while frequency decreases by 0.10 events/year. Such increases in tile drainage may lead to more severe streamflow droughts and have a detrimental effect on the socio-environmental usage of streams throughout the Midwest.

Illinois, Indiana, Iowa, Michigan, Minnesota, Ohio↗