USGS ScienceSearch

SEARCH · USGS Science

Results for “PLoS ONE”

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

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes substantially to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides substantial volumes of agricultural products not just for California but for the United States. Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface-water and groundwater resources in the valley. This situation has resulted in water management challenges in the Salinas Valley that generally relate to the distribution of the water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. Historically, to deal with the distribution issue, water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface-water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface-water supplies are insufficient. The availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand for groundwater and surface-water resources throughout the Salinas Valley has resulted in undesirable effects from unsustainable water use, such as surface-water depletion, groundwater-level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface-water and groundwater resources. To meet water demands and reduce the undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and climate variability. To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resources Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley hydrogeologic system called the Salinas Valley System Model. The geologic framework is known as the Salinas Valley Geologic Framework and was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The surface-water model is called the Salinas Valley Watershed Model and simulates the Salinas River watershed. Monthly surface-water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

California

Low-frequency earthquakes track the motion of a captured slab fragment

Accurate tectonic models are essential for assessing seismic hazard and fault interactions. However, the plate configuration at the complex Mendocino triple junction, where the San Andreas Fault and the Cascadia subduction zone meet, remains uncertain. We analyzed fault slip associated with a recently identified zone of tectonic tremor and low-frequency earthquakes (LFEs) near the southern edge of the subducting Gorda slab. Based on tidal sensitivity and P-wave first motions, we show that the LFEs are generated by dipping, strike-slip motion. This suggests that a former Farallon slab fragment, now captured by the Pacific plate, is translating northward beneath westernmost North America. This geometry effectively extends the slab interface fault, challenging prevailing interpretations of slab window formation and creating a potential unaccounted earthquake hazard in this region.

Science

Preparing for today's and tomorrow's water-resources challenges in eastern Long Island, New York

Freshwater is a vital natural resource. Although New York is a water-rich State, the wise and economical use of water resources is needed to ensure that there is enough water of adequate quality for both human and ecological needs—both for today and for tomorrow. Nowhere in New York is this more evident than in Nassau and Suffolk Counties on Long Island, where the public water supply is obtained from the sole-source aquifers located directly beneath the nearly 3 million people who live there. In 2023, in eastern Long Island’s Suffolk County, groundwater was pumped from these aquifers by more than 1,100 public water-supply wells to meet the needs of about 1.5 million people.

New York

Age and colony variation in Adélie penguin metapopulation vital rates: Insights from a 25-year mark–recapture study

Understanding how vital rates vary with age, life-history stage, and among populations is fundamental for predicting the demographic consequences of environmental change, especially in longer-lived species with complex life histories. These species often exhibit delayed maturity and iteroparity, making nuanced demographic insights critical for assessing their long-term viability. This study investigated age- and colony-related variation in Adélie penguin vital rates including survival, recruitment, and breeding propensity. We used mark–recapture data collected over 25 years (1996–2020) from three Adélie penguin breeding colonies that differed in population sizes and trends but comprised a metapopulation located at capes Royds, Bird, and Crozier on Ross Island, Antarctica. We used multi-state models to estimate survival and detection rates relative to reproductive state and breeding colony and estimated transition probabilities reflecting movements between reproductive states and colonies. Apparent survival varied by reproductive state, colony, and age and averaged 0.80 (SD = 0.02) at Bird, and 0.72 (SD = 0.03) and 0.73 (SD = 0.03) at Crozier and Royds, respectively, for pre-breeders age 2–7 years with strong declines in pre-breeder survival after age 8. We observed less age-related variation in survival of breeders and non-breeders, but we observed differences between colonies with lower survival for breeders (0.72 to 0.80) compared to non-breeders (0.75 to 0.82). The average probability of surviving the first 2 years after fledging ranged from 0.43 (SD = 0.14) at Royds and Crozier (0.43, SD = 0.10) to 0.55 (SD = 0.15) at Bird. Movement between colonies was highest for pre-breeders (0.00%–12.00% depending on age and colony) and lowest for breeders (<0.20%). We observed the lowest age-related recruitment rates at Royds, with recruitment at Crozier almost twice as high, and intermediate at Bird. Breeding propensity was highest at Crozier and lowest at Bird. Colony-specific variation in vital rates likely contributed differently to population trajectories, suggesting that care must be taken to extrapolate vital rate estimates across colonies even within a metapopulation. These findings also highlight the importance of considering age, life-history stage, and geographic variation when assessing population vital rates in species with complex life histories.

Frontiers in Ecology and Environment

Assessment and characterization of ephemeral stream channel stability and mechanisms affecting erosion in Grand Valley, western Colorado, 2018–21

The Grand Valley in western Colorado is in the semiarid Southwest United States. The north side of the Grand Valley has many ungaged ephemeral streams, which are of particular interest because (1) the underlying bedrock geology, Late Cretaceous Mancos Shale, is a sedimentary rock deposit identified as a major salinity contributor to the Colorado River and (2) despite infrequent streamflows of short duration, monsoon-derived floods in these ephemeral streams can carry substantial amounts of sediment downstream, affecting upstream and downstream banks and channel cross sections. The study area is of interest, because salinity, or the total dissolved solids concentration, in the Colorado River causes an estimated $300 million to $400 million per year in economic damages in the United States, and it is estimated 62 percent of the Upper Colorado River Basin’s total dissolved solid loads originate from geologic sources. In an effort to minimize salt contributions to the Colorado River from public lands administered by the Bureau of Land Management, a comprehensive salinity control approach is typically used to reduce nonpoint sources of salinity through land management techniques and practices. In 2018, the U.S. Geological Survey, in cooperation with the Bureau of Land Management, began an assessment of ephemeral streams located on the north side of the Grand Valley, western Colorado, to characterize stream channel stability and identify mechanisms affecting erosion. The U.S. Geological Survey developed a method for automatically extracting channel cross-section geometry from existing remotely sensed terrain models. Based on estimated flood stage and surrogate streamflows, hydraulic characteristics were calculated. Furthermore, the channel geometries and hydraulic characteristics were used to estimate channel stability using a statistical model. Cross-section stabilities were determined from a stream channel stability assessment for a subset of 1,406 visited (field observed) locations out of 13,415 cross sections, which were delineated from remotely sensed terrain models. The application of Manning’s resistance equation in combination with multiple logistic regression models demonstrated channel stability can be estimated with a 0.845 goodness of fit for a validation dataset when using a combination of drainage area, width-to-depth ratio, sinuosity, and shear stress as the explanatory variables. Stream channel stability was extrapolated for 13,415 unvisited (not field observed) cross sections using the multiple logistic regression model and defined explanatory variables. Mapping of the ephemeral streams and their associated stabilities may be used by the Bureau of Land Management to prioritize areas for remediation or changes in management strategies to reduce sediment and salinity loading to the Colorado River. The study found channel stability within the ephemeral streams to be spatially variable, longitudinally discontinuous, and dictated by changes in channel bed slope. The stable ephemeral streams were relatively wide and shallow and often had smaller drainage areas with less potential for producing shear stresses capable of overcoming channel adhesion. A change in channel bed slope can provide the means necessary to generate shear stresses appropriate to initiate erosion and a subsequent stability transition to incising channels. Channel widening happens when either or both banks of an incising channel reach a critical height for mass wasting, or when channel curvature causes higher sidewall stress. Regardless, widening channels can promote increases in sinuosity and subsequently reduce steep channel bed slopes. Consequently, stable and widening channels can have comparable bed slopes, making channel bed slope a poor explanatory variable to predict channel stability overall, despite its function to initiate channel instability. The results were based on a surrogate 0.10 annual exceedance probability (AEP; return period equal to the 10-year flood) interval streamflow, although it was recognized fluctuations in streamflow would also affect channel stability. Past and current changes within the study area affect streamflow; therefore, mechanisms affecting erosion include land use disturbances, soil compaction, loss of vegetation cover, drought, less frequent and more extreme precipitation, and fires—which all intensify the potential runoff and erosion within the study area.

Colorado

A joint Gaussian process model of geochemistry, geophysics, and temperature for groundwater TDS in the San Ardo Oil Field, California, USA

Decline in availability of fresh groundwater has expanded interest in brackish groundwater resources; however, the distribution of brackish groundwater is poorly understood. Water resources in sedimentary basins across the United States often overlie oil and gas development. Mapping of groundwater total dissolved solids (TDS) using data from oil well geophysical logs has become an important technique for identifying fresh and brackish groundwater. Existing geophysical log analysis methods use porosity and temperature to relate formation resistivity to TDS. Typically, natural geothermal gradients are used to estimate temperature at the location of collected resistivity. However, in thermally enhanced oil fields, steam is injected into the subsurface to mobilize high viscosity oil, creating variable temperature distributions. Furthermore, TDS derived from resistivity also depends on the fractions of dominant ions. Typically, chloride and bicarbonate fractions must be determined. It is also necessary to model TDS across many geologic units with heterogenous porosity distributions. Collectively, each quantity used to estimate TDS (resistivity, porosity, temperature, bicarbonate fraction) varies in space and time, and available data points are rarely collocated. Here, we present a new method of mapping groundwater TDS that continuously models each quantity together with a joint Gaussian process. This method enables mapping fresh and brackish water with practically available data. We apply this method to the San Ardo Oil Field in Monterey County, California, where steam injection occurs. In some areas of the aquifer system overlying the oil zone, the temperature is ∼75 °C, roughly twice the natural background value. Groundwater TDS is typically <1,500 mg/L in the aquifer and increases with depth to ∼9,000 mg/L in the oil-producing zone. A low-permeability clay layer delineates the fresh and brackish water, likely by inhibiting surface recharge from penetrating the deeper zones, allowing higher-TDS connate water to remain in place. Weaker lateral TDS trends may be controlled by recharge patterns associated with the Salinas River. Our model reveals with high certainty that groundwater has freshened in one localized part of the oil-producing zone and suggests with less certainty that more widespread freshening has also occurred. The lowering of TDS was possibly from decades of low-TDS steam injection and the associated fluid production and disposal operations.

California

Analysis of minimum 7-day discharges and estimation of minimum 7-day, 2-year discharges for streamflow-gaging stations in the Brazos River basin, Texas

The minimum 7-day, 2-year discharge (7Q2) was estimated for each of 97 streamflow-gaging stations in the Brazos River Basin, Texas, with at least 10 years of unregulated or 10 years of regulated discharge data using the Pearson Type III distribution fitted by L-moment statistics of the annual minimum 7-day discharges. Forty-six stations have at least 10 years of unregulated discharge data only, 26 stations have at least 10 years of regulated discharge data only, and 25 stations have at least 10 years of unregulated and 10 years of regulated discharge data. Statistical tests were applied to the annual minimum 7-day discharges to determine if the unregulated and regulated data have trends; and to determine if the data before regulation are significantly different from the data after regulation. Results of the Mann-Kendall test indicate a significant trend at the 95-percent confidence level for 8 of the 46 unregulated-only stations, for 6 of the 26 regulated-only stations, and for 5 stations before regulation and for 10 stations after regulation of the 25 stations with both unregulated and regulated data. Results of the Wilcoxon rank-sum test indicate a significant difference between the unregulated and regulated annual minimum 7-day discharges at the 95-percent confidence level for 16 of the 25 stations with both types of data. The 7Q2 is zero for 31 of the 46 unregulated-only stations and for 9 of the 26 regulated-only stations. The 7Q2 also is zero for 9 stations before regulation and for 3 stations after regulation of the 25 stations with both unregulated and regulated data. The 7Q2 ranges from 33 to 631 cubic feet per second for the eight regulated-only stations on the mainstem Brazos River. For the seven stations on the mainstem Brazos River with at least 10 years of unregulated and 10 years of regulated discharge data, the 7Q2 ranges from 0 to 568 cubic feet per second before regulation and from 0.30 to 670 cubic feet per second after regulation. The 7Q2 increased during regulation for 19 of the 25 stations with both unregulated and regulated data. The effect of regulation for most stations generally was an increase in the magnitude of the 7Q2.

Texas

Igneous differentiation in arcs

This chapter assesses crystallization-differentiation of convergent-margin magmas both from the standpoints of natural rock suites and from laboratory phase-equilibrium experiments. Crystallization-differentiation of dry subalkaline and alkaline basalts are summarized at the beginning for reference and context. Experiments crystallizing arc basalts through andesites at middle to upper crustal pressures (100–400 MPa), saturated either with H 2 O or mixed H 2 O-CO 2 , produce liquid lines of descent similar to arc igneous suites, with higher f O 2 s giving results closer to calc-alkaline magma series. Early enrichment of melts in SiO 2 with little or no enrichment in FeO* (total Fe as FeO) comes about by the combined effects of H 2 O suppressing crystallization of plagioclase relative to olivine, shifting plagioclase to more calcic – therefore lower SiO 2 – compositions, and suppressing silicate minerals to greater extents than Fe-rich spinel. Also contributing to SiO 2 enrichment is crystallization near the liquidus of higher proportions of olivine relative to calcic clinopyroxene, compared with MORB, but it is unclear if this is brought about by H 2 O or elevated pressure. Exhumed arc terranes, however, generally fail to expose evidence for primitive liquids having reached and differentiated at middle to upper crustal depths, at least in any abundance. Instead, most differentiation in arcs appears to take place in the lower crust or near the crust-mantle transition, and although hydrous experiments at such pressures produce liquids broadly similar to arc suites in their SiO 2 enrichment, the experimental liquids differ in becoming peraluminous at early stages of differentiation while melts have basaltic, basaltic andesitic, or andesitic SiO 2 concentrations. Crystallization of abundant calcic clinopyroxene close to the liquidus causes the high-pressure experimental liquids to become peraluminous with only small increases in SiO 2 concentration. In contrast, natural arc volcanic and plutonic suites become gradually less metaluminous with increasing SiO 2 , crossing to peraluminous at dacitic to rhyolitic SiO 2 concentrations. Hybridism, encompassing magma-mixing, assimilation of old crust, and incorporation of near-solidus melts from only shortly older deep, antecedent arc intrusions may accompany natural crystallization-differentiation, preventing melts from attaining peraluminous compositions at low-SiO 2 concentrations. If so, then nearly all arc suites are to substantial extents composites, rather than preserving simple, progressive liquid lines of descent. Some comments and assessments are also made on the participation of garnet in arc differentiation and on the production of arc batholiths.

Book chapter

Water-quality assessment of south-central Texas — Descriptions and comparisons of nutrients, pesticides, and volatile organic compounds at three intensive fixed sites, 1996-98

Water-quality samples were collected during April 1996-April 1998 at three intensive fixed sites in the San Antonio region of the South-Central Texas study unit as part of the U.S. Geological Survey National Water-Quality Assessment Program. The sampling strategy for the intensive fixed-site assessment is centered on obtaining information about the occurrence and seasonal patterns of selected constituents including nutrients, pesticides, and volatile organic compounds. The three sites selected to determine the effects of agriculture and urbanization on surface-water quality in the study unit are Medina River at LaCoste (agriculture indicator site), Salado Creek (lower station) at San Antonio (urban indicator site), and San Antonio River near Elmendorf (integrator site). Concentrations of two nutrients, dissolved nitrite plus nitrate nitrogen and total phosphorus, were largest at the integrator site, which is downstream of municipal wastewater treatment plants. Nitrite plus nitrate nitrogen concentrations at this site often exceeded the U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL) for drinking water. All total phosphorus concentrations at the site exceeded the EPA recommended maximum concentration for streams not discharging directly into reservoirs. Nitrite plus nitrate nitrogen concentrations at the integrator site tended to be smaller, and total phosphorus concentrations at the urban site tended to be larger in samples collected during stormflow than during base flow. The most detections and largest concentrations of three pesticides (atrazine, diazinon, and prometon) were in samples collected at the urban site. Some pesticide concentrations at the agriculture site showed a seasonal pattern of increasing concentrations during spring, the peak application season. Four pesticides (atrazine, deethylatrazine, diazinon, and prometon) were detected in at least 38 percent of samples collected at all three sites. The concentrations of all detected pesticides that have an MCL were less than the MCL at the three sites. More volatile organic compounds (VOC) were detected at the urban indicator site than at the agriculture indicator site, mostly likely because more sources are located in urbanized areas. The most VOCs detected and the largest concentrations of two VOCs (chloroform and tetrahydrofuran) were in samples from the integrator site. More VOCs were detected in samples collected at the integrator site during stormflow than during base flow. The concentrations of all detected VOCs that have an MCL were less than the MCL at the three sites.

Texas

Estimation of dynamic geologic CO2 storage resources in the Illinois Basin, including effects of brine extraction, anisotropy, and hydrogeologic heterogeneity

Since the vast majority of carbon dioxide (CO 2 ) storage resources in the United States are in deep saline aquifers, optimizing the use of these saline storage resources could be crucial for efficient development of geologic CO 2 storage (GCS) resources and basin- or larger-scale deployment of GCS in the country. Maximum CO 2 injection rates can be enhanced by extracting brine from the CO 2 storage unit. However, disposal of the extracted brine is both a technological and economic challenge. The lowest-cost option would likely be reinjection of the extracted brine into another formation above or below the CO 2 storage unit. Therefore, it is important to estimate brine injectivity as it will constrain the potential to increase CO 2 injectivity at an injection site that has access to multiple geologic storage units where either CO 2 or brine can be injected. Using a simulation-optimization framework, coupled with a non-isothermal, multiphase CO 2 -water-salt equation-of-state module, we developed a computationally efficient method for evaluating optimization of simultaneous CO 2 injection, brine extraction, and brine (re)injection at hypothetical injection sites deployed across a geologic basin. The Illinois basin is ideal for testing our methodology because it contains multiple geologic storage units with seals in between them to isolate injection of CO 2 in one unit from interfering with the injection of either brine or CO 2 in another unit above or below it. In addition, we investigated the relative effects of variation in key geologic parameters as well as two reservoir structures (hydrogeologic heterogeneity/anisotropy and homogeneity/isotropy) on CO 2 injectivities and enhancement of CO 2 injectivity through extracting brine. Results suggest that permeability, depth, and especially thickness of the storage unit could be the most influential parameters determining CO 2 injectivity. They also suggest that only injecting CO 2 into the storage unit with the greatest injectivity, enhancing that unit’s injectivity by extracting brine, and disposing of the produced brine in other suitable units could maximize total CO 2 injectivity in limited regions of the basin. At the majority of simulated injection sites, however, we found that injecting CO 2 into all of the accessible and suitable storage units was more likely to maximize the CO 2 storage resource.

Illinois, Kentucky, Indiana

Footprints of past mining in Alaska (USA) derived from high-resolution satellite imagery

Mapping the land area used for mining in the past is essential for guiding the remediation of affected landscapes and assessing the resource potential of related waste products. Despite significant recent progress delineating footprints of active and inactive mining globally, the known inventory of such mine lands remains incomplete. Here, I describe a new map dataset of footprints of land surface disturbance and waste at sites of past mining in Alaska (USA) based on visual interpretation of satellite imagery. This dataset maps 6–14 times the area of previous regional and global mine footprint maps in Alaska and is the first in the region to explicitly delineate mine waste landforms (e.g., tailings piles). The data are publicly available from the U.S. Geological Survey under a “no rights reserved” Creative Commons (CC0) license agreement.

Alaska

Water resources related to breccia pipe uranium mining in the Grand Canyon region

Introduction In the arid Grand Canyon region, water resources are limited to primarily the Colorado River and associated tributaries and to groundwater in the form of seeps and springs. Groundwater resources in the region supply water for human use and support diverse and rich ecosystems in the locations immediately surrounding the seeps and springs. Throughout the region, uranium resources occur and may interact with water resources in both mined and unmined uranium deposits. There is a need to better understand groundwater in the region and the effects from uranium mining in order to better manage the limited water resources in the area. This Fact Sheet summarizes results from U.S. Geological Survey studies that were conducted on this topic from 2012 to 2023.

Arizona

A history of cryohydrogeology modeling and recent advancements through the integration of solute transport

Groundwater flow systems and permafrost are interrelated because permafrost thaw enhances permeability, while groundwater flow can advect heat and accelerate permafrost thaw (McKenzie et al. 2021). Given amplified climate change in cold regions, there is renewed interest in ‘cryohydrogeology’, the study of groundwater in cold regions. Many data-driven studies have shown that permafrost thaw is leading to activated aquifers and increased baseflow across the pan-Arctic region (e.g. Walvoord and Striegl 2007, Evans et al. 2020). Empirical evidence of a subsurface ‘replumbing’ (Walvoord and Kurylyk 2016) in permafrost regions raises questions about the fate of sequestered contaminants in the North (Langer et al. 2023). We will discuss the history of and emerging opportunities in cryohydrogeological modeling, with a focus on recent contaminant transport modeling.

Conference Paper

Creation of aerial camera system boresight installations in Bangor, Maine; Monroe, Louisiana; and Bassville Park, Florida

Aerial camera systems rely on boresights for defining internal and external metrics. The U.S. Fish and Wildlife Service’s Migratory Bird Surveys Branch required boresights to calibrate their camera systems for creating georeferenced imagery that can then be used in geographic information systems. The U.S. Geological Survey created high-accuracy boresight installations in Bangor, Maine, in July 2024; Monroe, Louisiana, in September 2024; and Bassville Park, Florida, in January 2025 to provide this metric-defining capacity to the U.S. Fish and Wildlife Service. This U.S. Geological Survey Data Report details the methods and results of the boresight installation process and provides camera operators with the data needed to use each of the installations.

Florida, Louisiana, Maine

Hydrologic investigation of water level fluctuations at Moreau Lake, Moreau Lake State Park, town of Moreau, New York

The causes of water level fluctuations at Moreau Lake, within Moreau Lake State Park in the town of Moreau, New York, were investigated from 2016 to 2021 after lake water levels dropped between 2015 and 2016, raising concerns about the loss of a shallow swimming area at the park beach. Annual variation in precipitation records from the area did not account for the lake water level decline. Two possible causes for the low lake water levels were investigated: the increase in groundwater withdrawals from new residential development since about 2000 and seasonal changes (nongrowing and growing seasons) in precipitation. Investigation of the potential effects of nearby groundwater withdrawals required the compilation and collection of well-log data, seismic surveys, and measurements of lake and groundwater levels, field chemical parameters, and water isotopes to define the hydrogeologic system and to estimate water use. The net result of this work was the determination that Moreau Lake is a “flow though” lake with no surface water outlet; groundwater enters the lake on the upgradient side and exits through the downgradient side, however, groundwater does not flow southward from the lake toward nearby groundwater withdrawals from the semiconfined aquifer, and thus groundwater withdrawals were unlikely to have an effect on lake water levels. Investigation of the historic precipitation records during nongrowing (November through April) and growing (May through October) indicated that (1) nongrowing season precipitation from 2011–12 to 2015–16 was more deficient than any similar period during the past 78 years and (2) since about 2000, nongrowing seasons have been drier overall and growing seasons have been considerably wetter. Initiation of lake water level monitoring in 2016 provided an opportunity to compare seasonal precipitation with seasonal lake water level changes. Nongrowing season lake water levels are very sensitive to precipitation, such that high precipitation (40 percent above the seasonal median) resulted in a 5-foot rise in lake water level. In contrast, the growing season lake water levels are sensitive to dry conditions; for example, deficient rainfall (about 6 percent below the seasonal median) resulted in a decline of lake water levels of about 3.5 feet. However, lake water levels are insensitive to high growing season rainfall inputs (about 10 to 47 percent above the seasonal median); lake water levels consistently declined about by 0.8 feet above this range of seasonal excessive precipitation.

New York

Inset groundwater-flow models for the Cache and Grand Prairie Critical Groundwater Areas, northeastern Arkansas

The water resources in the Mississippi alluvial plain, located in parts of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas, supports a multibillion-dollar agricultural industry that relies heavily on pumping of groundwater for irrigation of crops and aquaculture. The primary source of groundwater for agricultural-related pumping is the Mississippi River Valley alluvial aquifer, which has declined in storage for decades; secondary groundwater sources include the middle Claiborne aquifer and Wilcox aquifer system. Two areas in northeastern Arkansas that lie within the Mississippi alluvial plain, part of the Cache and Grand Prairie regions, have been designated as Critical Groundwater Areas owing to decades of groundwater declines that resulted from past and current water use. The multidisciplinary Mississippi Alluvial Plain project, led by the U.S. Geological Survey, and funded by their Water Availability and Use Science Program, included objectives to develop numerical groundwater models in focus regions, including the part of the Cache and Grand Prairie regions of northeastern Arkansas. Two inset models were developed using the child model capabilities of MODFLOW 6, the U.S. Geological Survey’s Modular Hydrologic Model simulation software. Both models, called the Cache model and Grand Prairie model, simulated the groundwater system and surface-water/groundwater interactions for the Mississippi River Valley alluvial aquifer and underlying Tertiary-age aquifers and confining units to the Midway confining unit. Each model was spatially discretized into 500-meter x 500-meter orthogonal cells on a grid with 5-meter constant-thickness vertical layers that represented the Mississippi River Valley alluvial aquifer and increasing thickness layers for the aquifers and confining units below the alluvial aquifer. The Cache and Grand Prairie models were calibrated with the PEST++ iterative ensemble smoother Version 5 and employed high dimensional parameterization schemes of 13,740 and 30,436 parameters, respectively. The Cache mean absolute residual for groundwater-level observations within each model domain for the priority well was 1.58 meters. Grand Prairie mean absolute residuals for the alluvial aquifer and middle Claiborne aquifer groundwater-level observations were 2.71 and 10.78 meters, respectively. The groundwater budgets for the Cache and Grand Prairie models were characterized by substantial outflows to irrigation wells, which constituted about 52 and 54 percent of all outflows, with the primary source of water to those wells being releases from unconfined aquifer storage.

Arkansas

Ungulate migrations of the Western United States, volume 6

This report, volume 6 in the “Ungulate Migrations of the Western United States” report series, showcases the migrations of 23 ungulate herds in the Western United States. The report series is produced by the Corridor Mapping Team (CMT). Led by the U.S. Geological Survey, the CMT is a collaboration among 11 State agencies, as well as regional and Federal partners, and an expanding number of Tribal wildlife agencies. The CMT was initiated in response to the U.S. Department of the Interior Secretarial Order 3362, which was signed in 2018 and provided Federal support to expand existing research efforts to study ungulate populations and conserve their migrations throughout the Western United States. Including this volume, the report series has detailed the migrations of 237 unique ungulate herds throughout the Western United States and continues to serve as a valuable resource to guide local and regional management, policy, and on-the-ground work necessary to maintain intact and functional ungulate migrations. This report highlights several guiding principles of the CMT that facilitate collaboration among the diverse set of partners and contribute to the program’s continued successes. Notably, raw global positioning system data are not shared among participating agencies and the U.S. Geological Survey, delineating migration corridors and seasonal ranges relies on empirical data, the CMT provides flexible approaches to participating State and Tribal partners, and regular CMT meetings create a framework for open communication among agency partners that supports transboundary mapping of migrations. The 237 ungulate migrations that have been included in the report series are an expanding inventory, which can help maintain ungulate migrations in perpetuity.

Arizona, California, Colorado, Idaho, Montana, Nev

Land-use interactions, Oil-Field infrastructure, and natural processes control hydrocarbon and arsenic concentrations in groundwater, Poso Creek Oil Field, California, USA

Like many hydrocarbon production areas in the U.S., the Poso Creek Oil Field in California includes and is adjacent to other land uses (agricultural and other developed lands) that affect the hydrology and geochemistry of the aquifer overlying and adjacent to oil development. We hypothesize that the distributions of hydrocarbons and arsenic in groundwater in such areas will be controlled by complex interactions between mixed land uses, oil-field infrastructure, and natural processes. In 2020–2021, samples of groundwater and surface water were collected and analyzed for a large suite of inorganic and organic chemicals and isotope and gas tracers to test this hypothesis. Those data are supplemented with ancillary data on historical geochemistry, hydrology, geology, and oil-field infrastructure. Hydrocarbons in groundwater (e.g., methane through pentane gases and benzene) are associated with natural processes (e.g., fault offsets or transition in sediment depositional environment) and oil-field infrastructure (e.g., fluid-migration pathways associated with uncemented annulus in oil wells or unlined pits). Arsenic concentrations >10 μg per liter (μg/L; maximum concentration 12.9 μg/L) are associated with natural processes in old, high-pH groundwater, and more recent recharge of water from natural and/or engineered recharge processes. Along the southwest margin of the oil field, pumping for drinking-water and irrigation supplies in combination with engineered groundwater recharge produce a depression in groundwater elevations where groundwater with elevated sulfate concentrations from agricultural areas and groundwater with hydrocarbons from the oil field mix to produce a zone of sulfate reduction that removes hydrocarbons and arsenic from groundwater but produces elevated sulfide (S 2- ) concentrations (maximum concentration 29 mg per liter, mg/L). In this study, multiple approaches were required to resolve the overlapping effects of land uses, oil-field infrastructure, and natural processes on the distributions of hydrocarbons and arsenic in groundwater. The combined use of geographic, historical, physical, chemical, isotopic, and other information to constrain processes could be a useful approach for studies in other hydrocarbon-production areas. This is particularly important where land uses affect aquifer hydrology to an extent that causes mixing of groundwaters with different chemical compositions.

California