SEARCH · USGS Science

Results for “Nature Based Solutions”

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.

At least 289 records · Page 16Linked to original sources

USGS2023

Integrated rangeland fire management strategy actionable science plan completion assessment: Invasives topic, 2015–20

Loss and degradation of sagebrush rangelands due to an accelerated invasive annual grass-wildfire cycle and other stressors are significant management, conservation, and economic issues in the western United States. These sagebrush rangelands comprise a unique biome spanning 11 states, support over 350 wildlife species, and provide important ecosystem services that include stabilizing the economies of western communities. Impacts to sagebrush ecosystem processes over large areas due to the annual grass-wildfire cycle necessitated the development of a coordinated, science-based strategy for improving efforts to achieve long-term protection, conservation, and restoration of sagebrush rangelands, which was framed in 2015 under the Integrated Rangeland Fire Management Strategy (IRFMS). Central to this effort was the development of an Actionable Science Plan (Plan) that identified 37 priority science needs (Needs) for informing the actions proposed under the 5 topics (Fire, Invasives, Restoration, Sagebrush and Sage-Grouse, Climate and Weather) that were part of the collective focus of the IRFMS. Notable keys to this effort were identification of the Needs co-produced by managers and researchers, and a focus on resulting science being “actionable.” Substantial investments aimed at fulfilling the Needs identified in the Plan have been made since its release in 2016. While the state of the science has advanced considerably, the extent to which knowledge gaps remain relative to identified Needs is relatively unknown. Moreover, new Needs have likely emerged since the original strategy as results from actionable science reveal new questions and possible (yet untested) solutions. A quantifiable assessment of the progress made on the original science Needs can identify unresolved gaps and new information that can help inform prioritization of future research efforts. This report details a systematic literature review that evaluated how well peer-reviewed journal articles and formal technical reports published between January 1, 2015, and December 31, 2020, addressed six needs (hereinafter “Needs”) identified under the Invasives topic in the Plan. The topic outlined research Needs related to the control of invasive plant species in sagebrush rangelands, with a special emphasis on invasive annual grasses. We established the level of progress towards addressing each Need following a standardized set of criteria, and developed summaries detailing how research objectives nested within Needs identified in the Plan (“Next Steps”) were either addressed well, partially addressed, or remain outstanding (that is, addressed poorly) in the literature through 2020. Our searches resulted in the inclusion of 198 science products that at least partially addressed a Need identified in the Invasives topic. The Needs that were well and partially addressed included: studies of natural and anthropogenic factors influencing the distribution and spread; methods of preventing, eradicating and controlling invasive plant species; development of mapping techniques that provide regularly updated annual grass and fine fuel projections; and assessment of the efficacy of potential cheatgrass biocontrol agents. Needs that were addressed poorly included (1) investigations of livestock grazing as a tool for managing invasive plants and (2) investigations of cheatgrass die-offs and identification and subsequent study of potential biocontrol agents associated with those die-offs. The information provided in this assessment will assist updating the Plan along with other science strategies.

western United States↗

USGS2000

Self-ordering and complexity in epizonal mineral deposits

Epizonal base and precious metal deposits makeup a range of familiar deposit styles including porphyry copper-gold, epithermal veins and stockworks, carbonate-replacement deposits, and polymetallic volcanic rock-hosted (VHMS) deposits. They occur along convergent plate margins and are invariably associated directly with active faults and volcanism. They are complex in form, variable in their characteristics at all scales, and highly localized in the earth’s crust. More than a century of detailed research has provided an extensive base of observational data characterizing these deposits, from their regional setting to the fluid and isotope chemistry of mineral deposition. This has led to a broad understanding of the large-scale hydrothermal systems within which they form. Low salinity vapor, released by magma crystallization and dispersed into vigorously convecting groundwater systems, is recognized as a principal source of metals and the gases that control redox conditions within systems. The temperature and pressure of the ambient fluid anywhere within these systems is close to its vapor-liquid phase boundary, and mineral deposition is a consequence of short timescale perturbations generated by localized release of crustal stress. However, a review of occurrence data raises questions about ore formation that are not addressed by traditional genetic models. For example, what are the origins of banding in epithermal veins, and what controls the frequency of oscillatory lamination? What controls where the phenomenon of mineralization occurs, and why are some porphyry deposits, for example, so much larger than others? The distinctive, self-organized characteristics of epizonal deposits are shown to be the result of repetitive coupling of fracture dilation consequent on brittle failure, phase separation (“boiling”), and heat transfer between fluid and host rock. Process coupling substantially increases solute concentrations and triggers fast, far-from-equilibrium depositional processes. Since these coupled processes lead to localized transient changes in fluid characteristics, paragenetic, isotope, and fluid inclusion data relate to conditions at the site of deposition and only indirectly to the characteristics of the larger-scale hydrothermal system and its longer-term behavior. The metal concentrations (i.e. grade) of deposits and their internal variation is directly related to the geometry of the fracture array at the deposit scale, whereas finer-scale oscillatory fabrics in ores may be a result of molecular scale processes. Giant deposits are relatively rare and develop where efficient metal deposition is spatially focused by repetitive brittle failure in active fault arrays. Some brief case histories are provided for epithermal, replacement, and porphyry mineralization. These highlight how rock competency contrasts and feedback between processes, rather than any single component of a hydrothermal system, govern the size of individual deposits. In turn, the recognition of the probabilistic nature of mineralization provides a firmer foundation through which exploration investment and risk management decisions can be made.

Annual Review of Earth and Planetary Sciences↗

USGS1996

Low-flow characteristics of streams in Maryland and Delaware

Hydrologic information on the variability of streamflow during low-flow periods is needed for the effective management of surface-water resources in Maryland and Delaware. Low-flow characteristics derived from streamflow under natural conditions are presented for 94 continuousrecord gaging stations in Maryland, Delaware, and surrounding States, and for 131 low-flow partial-record gaging stations in Maryland and Delaware. Methods are developed to estimate low-flow characteristics at ungaged stream sites for average 7-, 14-, and 30-consecutive-day lowflow discharges for recurrence intervals of 2, 10, and 20 years. The methods are applicable to unregulated streams in Maryland and Delaware for watersheds having drainage areas ranging from 2.0 to 875 square miles. Maryland and Delaware are divided into four regions on the basis of physiography and the results of regression analyses. The regions are the Eastern Shore region in Maryland and Delaware, the eastern-Piedmont region in Maryland, the western region in Maryland, and the southern region in Maryland. Equations are developed from regression analyses to estimate lowflow characteristics at ungaged sites on ungaged streams for the Eastern Shore region, the easternPiedmont region, and the western region. A generalized least-squares multiple-regression technique is used to develop the estimation equations that relate average 7-, 14-, and 30- consecutive-day low-flow discharges for recurrence intervals of 2, 10, and 20 years to physical and climatological features of drainage basins in the study area. Adjustments are necessary to low-flow discharges estimated from the equation for the western region because of the effects of carbonate rocks on low flows. The adjustment method is based on analysis of the residuals of regional estimation equations and the percentage of the basin underlain by carbonate rocks. An equation is developed on the basis of a distance-weighted average of low-flow discharges at gaging stations to estimate low flow characteristics at ungaged sites on ungaged streams for the southern region. Estimates of the accuracy of low-flow characteristics determined from all regional equations are provided. Two methods are presented for estimating low-flow characteristics at ungaged sites on gaged streams. One method combines low-flow characteristics determined from gaging station data with results from the regional regression equation at the gaging station and at the ungaged site to estimate the low-flow characteristics either upstream or downstream from the gaging station. The other method estimates low-flow characteristics at ungaged sites between gaging stations using a graphical solution. Estimates of accuracy of the low-flow characteristics determined by these methods also are provided.

Delaware, Maryland↗

USGS2016

The source of groundwater and solutes to Many Devils Wash at a former uranium mill site in Shiprock, New Mexico

The Shiprock Disposal Site is the location of the former Navajo Mill (Mill), a uranium ore-processing facility, located on a terrace overlooking the San Juan River in the town of Shiprock, New Mexico. Following the closure of the Mill, all tailings and associated materials were encapsulated in a disposal cell built on top of the former Mill and tailings piles. The milling operations, conducted at the site from 1954 to 1968, created radioactive tailings and process-related wastes that are now found in the groundwater. Elevated concentrations of constituents of concern—ammonium, manganese, nitrate, selenium, strontium, sulfate, and uranium—have also been measured in groundwater seeps in the nearby Many Devils Wash arroyo, leading to the inference that these constituents originated from the Mill. These constituents have also been reported in groundwater that is associated with Mancos Shale, the bedrock that underlies the site. The objective of this report is to increase understanding of the source of water and solutes to the groundwater beneath Many Devils Wash and to establish the background concentrations for groundwater that is in contact with the Mancos Shale at the site. This report presents evidence on three working hypotheses: (1) the water and solutes in Many Devils Wash originated from the operations at the former Mill, (2) groundwater in deep aquifers is upwelling under artesian pressure to recharge the shallow groundwater beneath Many Devils Wash, and (3) the groundwater beneath Many Devils Wash originates as precipitation that infiltrates into the shallow aquifer system and discharges to Many Devils Wash in a series of springs on the east side of the wash. The solute concentrations in the shallow groundwater of Many Devils Wash would result from the interaction of the water and the Mancos Shale if the source of water was upwelling from deep aquifers or precipitation. In order to compare the groundwater from various wells to groundwater that has been affected by Mill activities, a classification system was developed to determine which wells were most likely to have been affected. Affects to groundwater by the Mill were determined by using the reported uranium alpha activity ratios measured in groundwater samples, along with the concentration of the uranium and the location of the wells relative to the Mill. Activity ratios of 1.2 or less were determined to be the most reliable indicator of Mill-affected groundwater. Wells with samples that had a reported activity ratio of 1.2 or less were classified as Mill affected. To compare groundwater with background water-quality, data from groundwater seeps and springs in the Upper Eagle Nest Arroyo and Salt Creek Wash, located north of the San Juan River, are also presented and analyzed. Based on groundwater elevations and tritium concentrations measured in wells located between the disposal cell and Many Devils Wash, Mill water is not likely to reach Many Devils Wash. The tritium concentrations also indicate that groundwater from the Mill has not substantially affected Many Devils Wash in the past. Upwelling from deep aquifers was also determined to be an unlikely source, primarily by comparing the composition of the stable isotopes of water in the shallow groundwater with those reported in groundwater samples from the deeper aquifers. The stable-isotope compositions of the shallow groundwater around the site are enriched relative to the San Juan River and local meteoric lines, which suggests that most of the shallow groundwater has been influenced by evaporation and therefore was recharged at the surface. Several observations indicate that focused recharge is the likely source of groundwater in the area of Many Devils Wash. The visible erosional features in Many Devils Wash provide evidence of piping and groundwater sapping, and the distribution and type of vegetation in Many Devils Wash suggest that the focused recharge of precipitation is occurring. The estimated recharge from precipitation was calculated to be 0.0008 inches per year (in/yr) by using the mass-balance approach from reported seep discharge and 0.0011 in/yr using the chloride mass-balance approach. A conceptual model of groundwater quality beneath Many Devils Wash is presented to explain the source of solutes in the groundwater beneath Many Devils Wash. The major-ion concentrations and geochemical evolution in the groundwater beneath Many Devils Wash and across the study area support the conceptual model that the underlying Mancos Shale is the source of solutes. Differences in the major-ion composition between groundwater samples collected around the site, result from the degree of weathering to the Mancos Shale. The cation distribution appears to be an indicator of effects from the Mill, with samples from the Mill-affected wells largely having a calcium/magnesium-sulfate composition that resembles the reported compositions of more weathered shale; however, that composition could change if the Mill-processed water flowed into areas where the Mancos Shale was less weathered. On the basis of the widespread presence of uranium in the Mancos Shale and the distribution of aqueous uranium in the analog sites and other sites in the region, it appears likely that uranium in the groundwater of Many Devils Wash is naturally sourced from the Mancos Shale.

New Mexico↗

USGS2004

Approaches to surface complexation modeling of Uranium(VI) adsorption on aquifer sediments

Uranium(VI) adsorption onto aquifer sediments was studied in batch experiments as a function of pH and U(VI) and dissolved carbonate concentrations in artificial groundwater solutions. The sediments were collected from an alluvial aquifer at a location upgradient of contamination from a former uranium mill operation at Naturita, Colorado (USA). The ranges of aqueous chemical conditions used in the U(VI) adsorption experiments (pH 6.9 to 7.9; U(VI) concentration 2.5 · 10 −8 to 1 · 10 −5 M; partial pressure of carbon dioxide gas 0.05 to 6.8%) were based on the spatial variation in chemical conditions observed in 1999–2000 in the Naturita alluvial aquifer. The major minerals in the sediments were quartz, feldspars, and calcite, with minor amounts of magnetite and clay minerals. Quartz grains commonly exhibited coatings that were greater than 10 nm in thickness and composed of an illite-smectite clay with occluded ferrihydrite and goethite nanoparticles. Chemical extractions of quartz grains removed from the sediments were used to estimate the masses of iron and aluminum present in the coatings. Various surface complexation modeling approaches were compared in terms of the ability to describe the U(VI) experimental data and the data requirements for model application to the sediments. Published models for U(VI) adsorption on reference minerals were applied to predict U(VI) adsorption based on assumptions about the sediment surface composition and physical properties (e.g., surface area and electrical double layer). Predictions from these models were highly variable, with results overpredicting or underpredicting the experimental data, depending on the assumptions used to apply the model. Although the models for reference minerals are supported by detailed experimental studies (and in ideal cases, surface spectroscopy), the results suggest that errors are caused in applying the models directly to the sediments by uncertain knowledge of: 1) the proportion and types of surface functional groups available for adsorption in the surface coatings; 2) the electric field at the mineral-water interface; and 3) surface reactions of major ions in the aqueous phase, such as Ca 2+ , Mg 2+ , HCO 3 − , SO 4 2− , H 4 SiO 4 , and organic acids. In contrast, a semi-empirical surface complexation modeling approach can be used to describe the U(VI) experimental data more precisely as a function of aqueous chemical conditions. This approach is useful as a tool to describe the variation in U(VI) retardation as a function of chemical conditions in field-scale reactive transport simulations, and the approach can be used at other field sites. However, the semi-empirical approach is limited by the site-specific nature of the model parameters.

Geochimica et Cosmochimica Acta↗

USGS2023

City-scale geothermal energy everywhere to support renewable resilience – A transcontinental cooperation

Cities have important and varying incentives to transform their energy sector to all-electric with low carbon emissions. However, they often encounter a number of impediments when attempting to implement such a change. For example, while urban areas have the highest energy demand-density, cities often lack the space for installing additional energy generation and/or long-duration energy storage systems. Cities also have existing environmental issues from energy sources (e.g., pollution from dust, waste heat or noise) that make residents sensitive to energy infrastructure development. Utilizing power from conventional sources, such as natural gas, biomass and hydropower, which usually are distanced from the urban areas, also make cities more vulnerable to supply disruptions. One promising de-carbonizing energy option for cities focuses on their heating and cooling needs, which constitutes around 1/3 of U.S. and 1/2 of European energy consumption (including industrial processes like drying, pasteurization, etc.; Jadun and others, 2017; EU Commission 2022). If heating and cooling loads can be met by geothermal direct-use technologies, then the need for new electric sources can be greatly lessened. Despite the proven efficacy of geothermal energy as a city/community-scale heating and cooling resource, it is currently only a niche resource in the heating and cooling sector, though has significant potential for future growth. Historically, emphasis has been placed on geothermal electricity generation potential that requires higher temperature (greater than 90 °C) resources at drillable depths, but potentially viable areas are geographically limited and typically well removed from urban centers. Key drivers for investments were represented by greater political interest in renewable electricity production, higher revenues and less effort in distributing the produced energy via grids. In contrast, low-temperature (less than 90 °C) geothermal resources can be used directly for heating and cooling almost everywhere and are cost-effective in urban/suburban settings. In addition, the increased prominence of renewable electricity sources, such as wind and solar onto city-scale electric grids, has led to new urgency around questions of energy storage. Underground thermal energy storage (UTES), wherein surplus or waste heat is stored underground for later use, could present a long-duration energy storage solution. From October 2022 through September 2024, a transcontinental consortium consisting of geological surveys, geoscience organizations, industry representatives and universities aims to develop an understanding of the global potential for city-scale geothermal, proposing guidelines to aid in promoting the economic utilization of low temperature geothermal resources. Efforts will focus on providing city managers and other decision makers with the information needed to evaluate and implement suitable city/community-scale geothermal technologies. Funded by the U.S. Geological Survey’s John Wesley Powell Center for Analysis and Synthesis, this interdisciplinary consortium will showcase tools, datasets, and scientific recommendations to accelerate the broader understanding and adoption of renewable energy systems that access geothermal resources. The collaborative research activities include standardization of nomenclatures, resource description and characterization strategies globally. The results from these activities will be combined with a preliminary climate-driven, city-based energy needs related analysis to perform energy supply/demand matching analysis. The identification of city-specific applications that would benefit from the geothermal technologies provides the basis to up-scale city-specific determinations to regional and national assessments of resource estimates. The city-scale geothermal energy research initiative will ultimately provide the synergies and management analysis that can address benefits, environmental impacts, regulatory frameworks, sustainability, and suitability in retrofitted buildings or new as well as existing heating networks.

Conference Paper↗

USGS2018

A comparison of synthetic flowpaths derived from light detection and ranging topobathymetric data and National Hydrography Dataset High Resolution Flowlines

Bathymetric and topobathymetric light detection and ranging (lidar) digital elevation models created for the Delaware River were provided to the National Geospatial Program and used to evaluate synthetic flowpath extraction from bathymetric/topobathymetric lidar survey data as a data source for improving the density, distribution, and connectivity of the National Hydrography Dataset High Resolution Flowline Network. As the surface-water component of The National Map, the National Hydrography Dataset maintains the Nation’s drainage network flow information and geometries for surface-water features used in hydrologic, hydraulic, and other science and engineering disciplines. The regional lidar survey for the Delaware River between Hancock, New York, and Trenton, New Jersey, was collected for the U.S. Geological Survey using the Experimental Advanced Airborne Research Lidar sensor system and processed by the Coastal National Elevation Database Applications Program. Using 1 percent of the maximum flow accumulation value for the surveyed Delaware River corridor as the flow accumulation threshold for grid cells at 1-, 5-, and 10-meter resolution created 223 to 283 kilometers of synthetic flowpaths potentially representing the river channel thalweg, which is the deepest point in a riverbed cross-section. There was potential for improving the High Resolution National Hydrography Dataset (HR NHD) Flowline network in places where the Delaware River channel, depicted as an Artificial Path in the HR NHD, is offset from the extracted synthetic river flowpath which sometimes appeared better positioned than the Artificial Path to represent the river thalweg. For the same area, using 0.05 percent of the maximum flow accumulation at the 1-, 5-, and 10-meter resolutions extracted 744 to 1,317 kilometers of synthetic flowpaths, with extracted synthetic flowpaths representing the main river channel and additional synthetic flowpaths representing tributaries or streams adjacent to the main channel. Overlaying these results with the HR NHDFlowline Network indicates that some of the additional synthetic flowpaths are connected to or extend HR NHD stream/river feature types. Some disconnected or isolated synthetic flowpaths not included in stream/river feature types were validated in orthoimagery and U.S. Topo Maps and provide examples of how extracted synthetic flowpaths could be used to delineate new stream/river features. Other additional extracted synthetic flowpaths depict linear features such as canals, tree lines, roads, or linear topographic depressions. For some river reaches where obstructions to flow or where low-relief topographic or bathymetric surfaces alter the flow direction, the software tool used to develop the flow direction grid did not calculate a primary flowpath for the river channel. Based on the results of this analysis, site conditions for the Delaware River corridor did not affect the quality of lidar bathymetric survey data. However, depending on the resolution of the lidar bathymetric digital elevation models (BDEMs), site conditions do have different effects on results for extracted synthetic flowpaths. We found that synthetic flowpaths extracted from 1-meter resolution lidar DEMs had more varied flow directions around in-channel landforms that obstructed flow than synthetic flowpaths extracted from 5- or 10-meter resolution lidar DEMs. As a result the 1-meter resolution DEM created some isolated or discontinuous synthetic flowpath segments where the 5- and 10-meter DEMs developed more continuous flowpaths. In this case the river bed upstream from the in-channel obstruction is shallower than the river bed downstream. Under these conditions the 1-meter resolution DEM provided synthetic flowpaths delineating a potential river thalweg. In this same area, the software solution modified (virtually raised) the river bed in the 5- and 10-meter resolution DEMs and flattened the bathymetric surface to create a continuous downstream flow direction, which caused trellis-patterned synthetic flowpaths to form. Under different site conditions and converse to the above development of synthetic flowpaths at different resolutions, at an abandoned river flood plain (terrace) with low relief that is adjacent to the river channel, the flow direction grid for the 1-meter resolution DEM developed continuous synthetic flowpath corresponding to a HR NHD Flowline network stream/river feature that connected to the main river channel but the larger resolution DEMs created isolated or disconnected synthetic flowpaths. A project to continue an evaluation of benefits of or issues caused by extracting synthetic flowpaths to enhance the HR NHD could include a study to assess the potential for merging surface-water flowpaths extracted from lidar topobathymetry and 3D Elevation Program digital elevation models. The merged DEM approach to synthetic flowpath extraction could extend the HR NHDFlowline network and enhance flow accumulations that might develop better flow direction grids in low-relief areas. Because of the confined lateral extent of the Delaware River, the lidar DEMs were not used to create catchments or watersheds; however, the merged DEM approach could also be tested as a resource for enhancing HR NHD catchments and watersheds. This lidar DEM synthetic flowpath extraction project supports the National Geospatial Program efforts to collect and produce high-quality lidar data to provide 3-dimensional representations of natural feature and aligns with the National Spatial Data Infrastructure to improve utilization of geospatial data. The results also can be useful for understanding strategies that can help maintain quality data in the HR NHD programs. KEYWORDS: bathymetric, digital elevation model, extracted synthetic flowpath, lidar, High Resolution National Hydrography Dataset, topobathymetric

New Jersey↗

USGS1992

Chemical, crystallographic and stable isotopic properties of alunite and jarosite from acid-Hypersaline Australian lakes

Chemical, crystallographic and isotopic analyses were made on samples containing alunite and jarosite from the sediments of four acid, hypersaline lakes in southeastern and southwestern Australia. The alunite and jarosite are K-rich with relatively low Na contents based on chemical analysis and determination of unit cell dimensions by powder X-ray diffraction. Correcting the chemical analyses of fine-grained mineral concentrates from Lake Tyrrell, Victoria, for the presence of halite, silica and poorly crystalline aluminosilicates, the following formulas indicate best estimates for solid-solution compositions: for alunite, K 0.87 Na 0.04 (H 3 O) 0.09 (Al 0.92 Fe 0.08 ) 3 (SO 4 ) 2 (OH) 6 and for jarosite, K 0.89 Na 0.07 (H 3 O) 0.04 (Fe 0.80 Al 0.20 ) 3 (SO 4 ) 2 (OH) 6 . The δD-values of alunite are notably larger than those for jarosite from Lake Tyrrell and it appears that the minerals have closely approached hydrogen isotope equilibrium with the acidic regional groundwaters. The δD results are consistent with a fractionation ∼60–70‰ between alunite and jarosite observed in other areas. However, interpretation of δD results is complicated by large variability in fluid δD H2O from evaporation, mixing and possible ion hydration effects in the brine. δD-values of water derived from jarosite by step-wise heating tend to be smaller at 250°C, at which temperature hydronium and other non-hydroxyl water is liberated, than at 550°C, where water is derived from the hydroxyl site, but the differences are not sufficiently different to invalidate measurements of total δD obtained by conventional, single-step heating methods. δ 34 S-values for alunite and jarosite from the four lakes (+19.7 to +21.2‰ CDT) and for aqueous sulfate from Lake Tyrrell (+18.3 to +19.8‰) are close to the values for modern evaporites ( +21.5 ±0.3‰ ">‰ +21.5 ±0.3‰ ) and seawater ( +20±0.5‰ ">‰ +20±0.5‰ ) and are probably typical of seawater-derived aerosols in arid coastal environments. δ 34 -S-values slightly smaller than that for seawater may reflect a minor contribution of sulfate from pyrite oxidation in the Parilla Sand or a reservoir effect from removal of gypsum enriched in 34 S. δ 18 O SO4 -values for alunite from three Western Australia lakes (+17.8 to +18.3‰ V-SMOW), for alunite and jarosite from Lake Tyrrell (+22.6 to +24.9‰) and for aqueous sulfate from Lake Tyrrell (+17.3 to +19.0‰) are much larger than the average value for seawater (+9.6‰). The data suggest an approach to 18 O- 16 O equilibrium between aqueous sulfate and groundwater, which is known from experimental studies to be possible at low pH and low temperatures, but has not been previously documented in nature. A residence time of ∼0.1–1 kyr for sulfate in acidic water (pH 3–4) is needed to achieve the apparent partial oxygen exchange, using previously published data of R.M. Lloyd.

Chemical Geology↗

USGS1980

Handbook of acute toxicity of chemicals to fish and aquatic invertebrates : summaries of toxicity tests conducted at Columbia National Fisheries Research Laboratory, 1965-78

Acute toxicity is a major subject of research at Columbia National Fisheries Research Laboratory for evaluating the impact of toxic chemicals on fishery resources. The Laboratory has played a leading role in developing research technology for toxicity testing and data interpretation. In 1965-78, more than 400 chemicals were tested against a variety of invertebrates and fish species representative of both cold- and warm-water climates. The use of acute toxicity tests for assessing the potential hazard of chemical contaminants to aquatic organisms is well documented (Boyd 1957; Henderson et al. 1960; Sanders and Cope 1966; Macek and McAllister 1970). Static acute toxicity tests provide rapid and (within limits) reproducible concentration-response curves for estimating toxic effects of chemicals on aquatic organisms. These tests provide a database for determining relative toxicity of a large number of chemicals to a variety of species and for estimating acute effects of chemical spills on natural aquatic systems; they also assist in determining priority and design of additional toxicity studies. Acute toxicity tests usually provide estimates of the exposure concentration causing 50% mortality (LC50) to test organisms during a specified period of time. For certain invertebrates, the effective concentration is based on immobilization, or some other identifiable endpoint, rather than on lethality. The application of the LC50 has gained acceptance among toxicologists and is generally the most highly rated test for assessing potential adverse effects of chemical contaminants to aquatic life (Brungs and Mount 1978; American Institute for Biological Sciences 1978 a ). The literature contains numerous papers dealing with the acute toxicity of chemicals to freshwater organisms. However, there is a tremendous need for a concise compendium of toxicity data covering a large variety of chemicals and test species. This Handbook is a compilation of a large volume of acute toxicity data from the Columbia Laboratory and its field laboratories. It presents definitive acute toxicity data on 271 chemicals tested against a variety of freshwater invertebrates and fishes. The chemicals represent all major groups of pesticides, as well as numerous industrial chemicals. This compilation should serve as a useful database for the many agencies and organizations dealing with research and management programs concerned with the impact of chemicals on aquatic resources. The Columbia Laboratory has played a major role in developing currently used standard methodology for static acute toxicity testing. The use of standardized methodology greatly reduces variation in results. The data presented here have been carefully scrutinized to eliminate tests that failed to follow acceptable procedures. Handling of test organisms and procedures for static toxicity tests followed those described by Lennon and Walker (1964) and Macek and McAllister (1970), and conform well with those recommended by Brauhn and Schoettger (1975) and the Committee on Methods for Toxicity Tests with Aquatic Organisms (1975). The species of fish and invertebrates that were tested are listed in phylogenetic order in Tables 1 and 2. Fish were obtained from Federal and State hatcheries as either eggs or fry. Original stocks of invertebrates were collected and cultured from wild populations with no known source of contamination; these populations were replenished regularly. The invertebrates were cultured in the Laboratory by methods similar to those described by Sanders and Cope (1966). Test chemicals usually consisted of technical or analytical grade samples of known purity. Formulations of the chemicals were also tested when available. When purity of test chemicals was known, all calculated concentrations were based on percent active ingredients. Stock solutions were prepared immediately before each test, with commercial grade acetone as the carrier solvent. Occasionally, ethanol or dimethyl-formamide was substituted. Solvent concentrations did not exceed 0.5 mL/L in final dilution water. Test water (dilution water) was reconstituted from deionized water of at least 10 6 ohms resistivity by the addition of appropriate reagent grade chemicals (Marking 1969). Water was buffered to maintain a pH of 7.2 to 7.5, an alkalinity of 30 to 35 mg/L, and a hardness of 40 to 50 mg/L as CaCO 3 . Test water was mixed thoroughly and aerated before transfer into test chambers. Fish were acclimated to dilution water by gradually changing the water in acclimated tanks from 100% well water to 100% reconstituted water over a 1- to 3-day period at the desired testing temperature. Invertebrates were acclimated from well water to dilution water over a 4- to 6-h period. Toxicity tests were conducted under static conditions without aeration, and the organisms were not fed during acclimation or testing. Temperature of test solutions was maintained within ± 1°C of that required for a given test. Toxicity tests with fish were conducted in 18.9-liter (5-gal) wide-mouthed jars containing 15 liters of test solution. Fingerling fish weighing 0.2 to 1.5 g were tested at each concentration. Caution was taken not to exceed 0.8 g of test organisms per liter of solution. Duplicate test chambers were used to accommodate larger fish. Test chambers varied in size for invertebrates, depending on the species used; volume of test solution ranged from 0.25 to 4 liters. At least 10 organisms were exposed to each concentration for all definitive tests. At least six concentrations were used per toxicity test. The tests began upon initial exposure to the toxicant and continued for 96 h. Immobilization tests with invertebrates were conducted for only 48 h. The number of dead or affected organisms in each test chamber were recorded and the dead organisms were removed every 24 h; general observations on the condition of test organisms were also recorded at these times. Toxicity data were analyzed by a statistical method described by Litchfield and Wilcoxon (1949) to determine LC50 (theoretical estimate of the concentration lethal to 50% of the test animals) and 95% confidence intervals. This method is recommended by the American Public Health Association (1971) and by Sprague (1969) for determining median lethal concentrations. The procedure is easily modified for computing a single LC50 when replicate tests are performed.

Resource Publication↗

USGS2008

Development, Testing, and Sensitivity and Uncertainty Analyses of a Transport and Reaction Simulation Engine (TaRSE) for Spatially Distributed Modeling of Phosphorus in South Florida Peat Marsh Wetlands

Alterations to the predevelopment delivery of water and nutrients into the Everglades of southern Florida have been occurring for nearly a century. Major regional drainage projects, large-scale agricultural development, and changes to the hydrology of the Kissimmee River-Lake Okeechobee watershed have resulted in substantial phosphorus transport increases by surface waters. Excess phosphorus has accumulated in the soils of northern Everglades marshes to levels that have impaired the natural resources of the region. Regulations now limit the amount of phosphorous that enters the Everglades through an extensive network of water-control structures. This study involved the development and application of water-quality modeling components that may be applied to existing hydrologic models of southern Florida to evaluate the effects of different management scenarios. The result of this work is a spatially distributed water-quality model for phosphorus transport and cycling in wetlands. The model solves the advection-dispersion equation on an unstructured triangular mesh and incorporates a wide range of user-selectable mechanisms for phosphorus uptake and release parameters. In general, the phosphorus model contains transfers between stores; examples of stores that can be included are soil, water column (solutes), pore water, macrophytes, suspended solids (plankton), and biofilm. Examples of transfers are growth, senescence, settling, diffusion, and so forth, described with first order, second order, and Monod types of transformations. Local water depths and velocities are determined from an existing two-dimensional, overland-flow hydrologic model. The South Florida Water Management District Regional Simulation Model was used in this study. The model is applied to three case studies: intact cores of wetland soils with water, outdoor mesocosoms, and a large constructed wetland; namely, Cell 4 of Stormwater Treatment Area 1 West (STA-1W Cell 4). Different levels of complexity in the phosphorus cycling mechanisms were simulated in these case studies using different combinations of phosphorus reaction equations. Changes in water column phosphorus concentrations observed under the controlled conditions of laboratory incubations, and mesocosm studies were reproduced with model simulations. Short-term phosphorus flux rates and changes in phosphorus storages were within the range of values reported in the literature, whereas unknown rate constants were used to calibrate the model output. In STA-1W Cell 4, the dominant mechanism for phosphorus flow and transport is overland flow. Over many life cycles of the biological components, however, soils accrue and become enriched in phosphorus. Inflow total phosphorus concentrations and flow rates for the period between 1995 and 2000 were used to simulate Cell 4 phosphorus removal, outflow concentrations, and soil phosphorus enrichment over time. This full-scale application of the model successfully incorporated parameter values derived from the literature and short-term experiments, and reproduced the observed long-term outflow phosphorus concentrations and increased soil phosphorus storage within the system. A global sensitivity and uncertainty analysis of the model was performed using modern techniques such as a qualitative screening tool (Morris method) and the quantitative, variance-based, Fourier Amplitude Sensitivity Test (FAST) method. These techniques allowed an in-depth exploration of the effect of model complexity and flow velocity on model outputs. Three increasingly complex levels of possible application to southern Florida were studied corresponding to a simple soil pore-water and surface-water system (level 1), the addition of plankton (level 2), and of macrophytes (level 3). In the analysis for each complexity level, three surface-water velocities were considered that each correspond to residence times for the selected area (1-kilometer long) of 2, 10, and 20

Scientific Investigations Report↗

USGS2025

Refining the earthquake history of south-central Alaska through lake records

The Alaska–Aleutian subduction zone (AASZ) is one of the world's most seismically active plate boundaries and the source of the 1964 M w 9.2 Great Alaska earthquake–the second largest instrumentally recorded earthquake in the world. Understanding the nature and frequency of such earthquakes is necessary for seismic and tsunami hazard assessment, but instrumental and historical records that span less than 150 years are too short to allow a statistically reliable analysis of earthquake recurrence times. This calls for studies of evidence of past earthquakes, extending the earthquake catalog further back in time. Subduction-zone paleoseismology in south-central Alaska is predominantly based on coastal evidence of land-level changes and tsunamis generated by megathrust earthquakes and preserved in the geological record. A complementary approach is lacustrine paleoseismology, which is still a relatively young discipline in Alaska. However, globally, lake basins are well-established high-resolution and continuous recorders of paleoseismic activity along subduction zones, relying on the identification of underwater landslide deposits and turbidites generated by seismic shaking. As a result, lake basins not only register ground shaking from megathrust earthquakes, but also from intraslab and crustal earthquakes, which are typically not accompanied by significant land-level changes. In this review paper, we combine coastal and lacustrine paleoseismology approaches to refine the south-central Alaskan earthquake history by comparing the paleoseismic records from two lakes (i.e., Eklutna Lake, located in the Chugach Mountain Range, and Skilak Lake, situated on the Kenai Peninsula) with the coastal and crustal earthquake catalog in Alaska. The resulting age ranges of all known megathrust earthquakes involving the Alaskan megathrust between the Kodiak and Prince William Sound (PWS) sections are more precise and accurate for the last 1.3 kyrs BP than the previously published age ranges from coastal records. As a result, this study supports the following key conclusions: (1) The 1964 CE earthquake was an exceptionally strong and unique event in the last 2000 years, rupturing the PWS, Kenai, Barren Islands, and Kodiak sections simultaneously. (2) The high-resolution and seasonal markings of the varved lake records now disentangle for the first time closely timed earthquakes, which was not possible based on the coastal evidence alone. (3) No persistent megathrust rupture boundaries exist. So, the possibility of a full rupture of the entire eastern AASZ, from PWS to Semidi cannot be excluded. (4) The rupture pattern in the eastern AASZ reveals superimposed cycles of multi-asperity ruptures (1964 earthquake) and clustered complementary partial ruptures, or rupture cascades. (5) The PWS section hosts the largest asperity in the eastern AASZ. (6) The shaking record of megathrust earthquakes indicates a time-dependent (quasiperiodic) behavior for the study area, but the observation of complementary clusters means that the hazards will not drop to zero but instead may even increase for a neighboring section. (7) The time-independent behavior of intraplate earthquakes implies that the intraslab hazard did not decrease following the 2016 and 2018 earthquakes. This study utilizes an integrated approach for subduction zone paleoseismology as a solution for unraveling recurrence and rupture patterns in Alaska, which can be applied worldwide.

Alaska↗

USGS2003

Anaerobic degradation of 1,1,2,2-tetrachloroethane and association with microbial communities in a freshwater tidal wetland, Aberdeen Proving Ground, Maryland: Laboratory experiments and comparisons to field data

Defining biodegradation rates and processes is a critical part of assessing the feasibility of monitored natural attenuation as a remediation method for ground water containing organic contaminants. During 1998–2001, the U.S. Geological Survey conducted a microbial study at a freshwater tidal wetland along the West Branch Canal Creek, Aberdeen Proving Ground, Maryland, as part of an investigation of natural attenuation of chlorinated volatile organic compounds (VOCs) in the wetland sediments. Geochemical analyses and molecular biology techniques were used to investigate factors controlling anaerobic degradation of 1,1,2,2-tetrachloroethane (TeCA), and to characterize the microbial communities that potentially are important in its degradation. Rapid TeCA and daughter product degradation observed in laboratory experiments and estimated with field data confirm that natural attenuation is a feasible remediation method at this site. The diverse microbial community that seems to be involved in TeCA degradation in the wetland sediments varies with changing spatial and seasonal conditions, allowing continued effective natural attenuation throughout the year. Rates of TeCA degradation in anaerobic microcosm experiments conducted with wetland sediment collected from two different sites (WB23 and WB30) and during three different seasons (March–April 1999, July–August 1999, and October–November 2000) showed little spatial variability but high seasonal variability. Initial first-order degradation rate constants for TeCA ranged from 0.10±0.01 to 0.16±0.05 per day (halflives of 4.3 to 6.9 days) for March–April 1999 and October–November 2000 microcosms incubated at 19 degrees Celsius, whereas lower rate constants of 0±0.03 and 0.06±0.03 per day were obtained in July–August 1999 microcosms incubated at 19 degrees Celsius. Microbial community profiles showed that low microbial biomass and microbial diversity in the summer, possibly due to competition for nutrients by the wetland vegetation, could account for these unexpectedly low degradation rates. In microcosms incubated at 5 degrees Celsius, about 50 percent of the initial TeCA in solution was converted to daughter products within a 35-day incubation period, indicating that biodegradation in the wetland sediments can continue during cold winter temperatures. Initial pathways of TeCA degradation were the same in the wetland sediment microcosms regardless of the season or sediment collection site, the reduction-oxidation conditions, and the previous exposure of the sediment to contamination. Immediate and simultaneous dichloroelimination and hydrogenolysis, producing 1,2-dichloroethene (12DCE) and 1,1,2-trichloroethane (112TCA), respectively, were the initial TeCA degradation pathways in all live microcosm experiments. The production and degradation of vinyl chloride (VC), which is the most toxic of the TeCA daughter compounds, was affected by spatial and seasonal variability, reduction-oxidation condition, and pre-exposure of the wetland sediment. TeCA-amended microcosms constructed with WB30 sediment showed approximately twice as much VC production as those constructed with WB23 sediment. Results of 112TCA-amended microcosms indicated that the greater production of VC in the WB30 sediment resulted from a greater predominance of the 112TCA dichloro elimination pathway in these sediments. VC degradation also was substantially higher in microcosms constructed with WB30 sediment than those constructed with WB23 sediment, resulting in lower VC concentrations at the end of WB30 microcosms. Enrichment experiments in which microcosm slurry was amended with high initial VC concentrations showed that the spatial difference in VC degradation was negligible after prolonged incubation under methanogenic conditions. Inhibition of methanogenic activity in microcosms by addition of sulfate or of 2-bromoethanesulfonic acid inhibited production and degradation of VC. Inhibition of methanogenesis by addition of ferric iron or of 2-bromoethanesulfonic acid also completely inhibited VC degradation in VC-amended enrichment experiments. Pre-exposure to VC substantially increased degradation in VC-amended enrichment experiments. A microbial consortium, rather than one microbial species or group, likely is involved in the degradation of TeCA, as indicated by the occurrence of multiple degradation pathways and the variability in VC production and degradation. A bacterial peak at 90 base pair (bp) fragment length in terminal- restriction fragment length polymorphism (TRFLP) profiles was associated with TeCA hydrogenolysis to 112TCA, and bacterial species represented by 198 and 170 bp fragment lengths were associated with TeCA dichloroelimination to 12DCE. Dichloroelimination of 112TCA to VC was associated with increasing dominance of the 198 bp bacterial peak in March–April 1999 and October–November 2000 microcosms, whereas an 86 bp or the 170 bp bacterial peak was associated with 112TCA dichloroelimination in the summer experiment. Hydrogenolysis of 12DCE to VC was associated with a carbon dioxide-utilizing methanogen at 307 bp in the March–April 1999 and October–November 2000 microcosm experiments, whereas production of VC occurred despite low methanogen biomass and methane production in the July–August 1999 experiments. Production of VC in the absence of methane production also occurred in 12DCE-amended enrichment cultures. The exponential production of VC in the 12DCE-amended enrichment cultures after an initial lag indicated growth of a microbial species or group, possibly one of the known dehalorespiring bacteria. Molecular analyses using specific primers targeting dehalorespiring bacteria of the Dehalococcoides group (Dehalococcoides ethenogenes and Dehalococcoides sp. strain FL2) and of the acetate-oxidizing Desulfuromonas group (Desulfuromonas sp. strain BB1 and Desulfuromonas chloroethenica) showed the presence of these bacteria in microcosm slurry from site WB30 but not from site WB23. Addition of hydrogen, which is the favored substrate of Dehalococcoides, tripled VC production in 12DCE-amended enrichment cultures. VC degradation showed a marked association with an increase in the relative proportion of Methanosarcinaceae, a family of methanogens that includes all those capable of utilizing acetate as a substrate, in the total methanogen community. Half-lives for TeCA and TCE estimated from field data were in the range of 60 to 100 days, which agrees well with laboratory estimates of degradation rates considering the inherent differences in the laboratory and field systems. Both laboratory microcosm experiments and field data showed that 12DCE and VC are the predominant, persistent daughter compounds from TeCA degradation. In addition, porewater chemistry showed higher accumulation of VC in the wetland sediment at site WB30 than at site WB23, as was observed in the microcosm experiments. Molecular analyses of grab samples of surficial wetland sediment showed that all the microbial species or groups linked to TeCA degradation in the microcosm experiments were present in all sediment samples. Microbial biomass and diversity were lowest in an area of the wetland (transect C-C') where porewater VOC concentrations are highest, indicating that the higher VOC concentrations could result from lower degradation rates. The lower microbial biomass and diversity in this area could be caused by toxic effects of the contaminants, or possibly from differences in frequency and duration of tidal inundation.

Maryland↗

USGS2016

Application of hydrogeology and groundwater-age estimates to assess the travel time of groundwater at the site of a landfill to the Mahomet Aquifer, near Clinton, Illinois

The U.S. Geological Survey used interpretations of hydrogeologic conditions and tritium-based groundwater age estimates to assess the travel time of groundwater at a landfill site near Clinton, Illinois (the “Clinton site”) where a chemical waste unit (CWU) was proposed to be within the Clinton landfill unit #3 (CLU#3). Glacial deposits beneath the CWU consist predominantly of low-permeability silt- and clay-rich till interspersed with thin (typically less than 2 feet in thickness) layers of more permeable deposits, including the Upper and Lower Radnor Till Sands and the Organic Soil unit. These glacial deposits are about 170 feet thick and overlie the Mahomet Sand Member of the Banner Formation. The Mahomet aquifer is composed of the Mahomet Sand Member and is used for water supply in much of east-central Illinois. Eight tritium analyses of water from seven wells were used to evaluate the overall age of recharge to aquifers beneath the Clinton site. Groundwater samples were collected from six monitoring wells on or adjacent to the CLU#3 that were open to glacial deposits above the Mahomet aquifer (the upper and lower parts of the Radnor Till Member and the Organic Soil unit) and one proximal production well (approximately 0.5 miles from the CLU#3) that is screened in the Mahomet aquifer. The tritium-based age estimates were computed with a simplifying, piston-flow assumption: that groundwater moves in discrete packets to the sampled interval by advection, without hydrodynamic dispersion or mixing. Tritium concentrations indicate a recharge age of at least 59 years (pre-1953 recharge) for water sampled from deposits below the upper part of the Radnor Till Member at the CLU#3, with older water expected at progressively greater depth in the tills. The largest tritium concentration from a well sampled by this study (well G53S; 0.32 ± 0.10 tritium units) was in groundwater from a sand deposit in the upper part of the Radnor Till Member; the shallowest permeable unit sampled by this study. That result indicated that nearly all groundwater sampled from well G53S entered the aquifer as recharge before 1953. Tritium was detected in a trace concentration in one sample from a second monitoring well open to the upper part of the Radnor Till Member (well G07S; 0.11 ± 0.09 tritium units), and not detected in samples collected from two monitoring wells open to a sand deposit in the lower part of the Radnor Till Member, from two samples collected from two monitoring wells open to the Organic Soil unit, and in two samples collected from a production well screened in the middle of the Mahomet aquifer (a groundwater sample and a sequential replicate sample). The lack of tritium in five of the six groundwater samples collected from the shallow permeable units beneath CLU#3 site and the two samples from the one Mahomet aquifer well indicates an absence of post-1952 recharge. Groundwater-flow paths that could contribute post-1952 recharge to the lower part of the Radnor Till Member, the Organic Soil unit, or the Mahomet aquifer at the CLU#3 are not indicated by these data. Hypothetical two-part mixtures of tritium-dead, pre-1953 recharge water and decay-corrected tritium concentrations in post-1952 recharge were computed and compared with tritium analyses in groundwater sampled from monitoring wells at the CLU#3 site to evaluate whether tritium concentrations in groundwater could be represented by mixtures involving some post-1952 recharge. Results from the hypothetical two-part mixtures indicate that groundwater from monitoring well (G53S) was predominantly composed of pre-1953 recharge and that if present, younger, post-1955 recharge, contributed less than 2.5 percent to that sample. The hypothetical two-part mixing results also indicated that very small amounts of post-1952 recharge composing less than about 2.5 percent of the sample volume could not be distinguished in groundwater samples with tritium concentrations less than about 0.15 TU. The piston-flow based age of recharge determined from the tritium concentration in the groundwater sample from monitoring well G53S yielded an estimated maximum vertical velocity from the land surface to the upper part of the Radnor Till Member of 0.85 feet per year or less. This velocity, ifassumed to apply to the remaining glacial till deposits above the Mahomet aquifer, indicates that recharge flows through the 170 feet of glacial deposits between the base of the proposed chemical waste unit and the top of the Mahomet aquifer in a minimum of 200 years or longer. Analysis of hydraulic data from the site, constrained by a tritium-age based maximum groundwater velocity estimate, computed minimum estimates of effective porosity that range from about 0.021 to 0.024 for the predominantly till deposits above the Mahomet aquifer. Estimated rates of transport of recharge from land surface to the Mahomet aquifer for the CLU#3 site computed using the Darcy velocity equation with site-specific data were about 260 years or longer. The Darcy velocity-based estimates were computed using values that were based on tritium data, estimates of vertical velocity and effective porosity and available site-specific data. Solution of the Darcy velocity equation indicated that maximum vertical groundwater velocities through the deposits above the aquifer were 0.41 or 0.61 feet per year, depending on the site-specific values of vertical hydraulic conductivity (laboratory triaxial test values) and effective porosity used for the computation. The resulting calculated minimum travel times for groundwater to flow from the top of the Berry Clay Member (at the base of the proposed chemical waste unit) to the top of the Mahomet aquifer ranged from about 260 to 370 years, depending on the velocity value used in the calculation. In comparison, plausible travel times calculated using vertical hydraulic conductivity values from a previously published regional groundwater flow model were either slightly less than or longer than those calculated using site data and ranged from 230 to 580 years. Tritium data from 1996 to 2011 USGS regional sampling of groundwater from domestic wells in the confined part of the Mahomet aquifer—which are 2.5 to about 40 miles from the Clinton site—were compared with site-specific data from a production well at the Clinton site. Tritium-based groundwater-age estimates indicated predominantly pre- 1953 recharge dates for USGS and other prior regional samples of groundwater from domestic wells in the Mahomet aquifer. These results agreed with the tritium-based, pre-1953 recharge age estimated for a groundwater sample and a sequential replicate sample from a production well in the confined part of the Mahomet aquifer beneath the Clinton site. The regional tritium-based groundwater age estimates also were compared with pesticide detections in samples from distal domestic wells in the USGS regional network that are about 2.5 to 40 miles from the Clinton site to identify whether very small amounts of post-1952 recharge have in places reached confined parts of the Mahomet aquifer at locations other than the Clinton site in an approximately 2,000 square mile area of the Mahomet aquifer. Very small amounts of post-1952 recharge were defined in this analysis as less than about 2.5 percent of the total recharge contributing to a groundwater sample, based on results from the two-part mixing analysis of tritium data from the Clinton site. Pesticide-based groundwater-age estimates based on 22 detections of pesticides (13 of these detections were estimated concentrations), including atrazine, deethylatrazine (2-Chloro-4-isopropylamino-6-amino- s-triazine), cyanazine, diazinon, metolachlor, molinate, prometon, and trifluralin in groundwater samples from 10 domestic wells 2.5 to about 40 miles distant from the Clinton site indicate that very small amounts of post-1956 to post-1992 recharge can in places reach the confined part of the Mahomet aquifer in other parts of central Illinois. The relative lack of tritium in these samples indicate that the amounts of post-1956 to post-1992 recharge contributing to the 10 domestic wells were a very small part of the overall older groundwater sampled from those wells. The flow process by which very small amounts of pesticide-bearing groundwater reached the screened intervals of the 10 domestic wells could not be distinguished between well-integrity related infiltration and natural hydrogeologic features. Potential explanations include: (1) infiltration through man-made avenues in or along the well, (2) flow of very small amounts of post-1956 to post-1992 recharge through sparsely distributed natural permeable aspects of the glacial till and diluted by mixing with older groundwater, or (3) a combination of both processes. Presuming the domestic wells sampled by the USGS in 1996–2011 in the regional study of the confined part of the Mahomet aquifer are adequately sealed and produce groundwater that is representative of aquifer conditions, the regional tritium and pesticide-based groundwater-age results indicate substantial heterogeneity in the glacial stratigraphy above the Mahomet aquifer. The pesticide-based groundwater-age estimates from the domestic wells distant from the Clinton site also indicate that parts of the Mahomet aquifer with the pesticide detections can be susceptible to contaminant sources at the land surface. The regional pesticide and tritium results from the domestic wells further indicate that a potential exists for possible contaminants from land surface to be transported through the glacial drift deposits that confine the Mahomet aquifer in other parts of central Illinois at faster rates than those computed for recharge at the Clinton site, including CLU#3. This analysis indicates the potential value of sub-microgram-per-liter level concentrations of land-use derived indicators of modern recharge to indicate the presence of very small amounts of modern, post-1952 age recharge in overall older, pre-1953 age groundwater.

Illinois↗

USGS2002

Geology and hydrogeology of the Caribbean Islands aquifer system of the Commonwealth of Puerto Rico and the U.S. Virgin Islands

Poorly lithified to unconsolidated carbonate and clastic sedimentary rocks of Tertiary (Oligocene to Pliocene) and Quaternary (Pleistocene to Holocene) age compose the South Coast aquifer and the North Coast limestone aquifer system of Puerto Rico; poorly lithified to unlithified carbonate rocks of late Tertiary (early Miocene to Pliocene) age make up the Kingshill aquifer of St. Croix, U.S. Virgin Islands. The South Coast aquifer, North Coast limestone aquifer system, and Kingshill aquifer are the most areally extensive and function as the major sources of ground water in the U.S. Caribbean Islands Regional Aquifer-System Analysis (CI-RASA) study area. In Puerto Rico's South Coast ground-water province, more than 1,000 meters of clastic and carbonate rocks of Oligocene to Pliocene age infill the South Coast Tertiary Basin. The pattern of lithofacies within this basin appears to have been controlled by changes in base level that were, at times, dominated by tectonic movement (uplift and subsidence), but were also influenced by eustasy. Deposition of the 70-kilometer long and 3- to 8-kilometer wide fan-delta plain that covers much of the South Coast ground-water province occurred largely in response to glacially-induced changes in sea level and climate during the Quaternary period. Tectonic movement played a much less important role during the Quaternary. The North Coast ground-water province of Puerto Rico is underlain by homoclinal coastal plain wedge of carbonate and siliciclastic rocks that infill the North Coast Tertiary Basin and thicken to more than 1,700 meters. A thin basal siliciclastic sequence of late Oligocene age is overlain by a thick section of mostly carbonate rocks of Oligocene to middle Miocene age. Globigerinid limestone of late Miocene to Pliocene age crops out and lies in the shallow subsurface areas of northwestern Puerto Rico. Oligocene to middle Miocene age rocks tentatively can be divided into five depositional sequences and associated systems tracts; these rocks record carbonate and minor siliciclastic deposition that occurred in response to changes in relative sea level. The Cibao Formation represents the most complex of these sequences and contains a varied facies of carbonate, mixed carbonate-siliciclastic, and siliciclastic rocks that reflect differential uplift, subsidence, and transgression of the sea. Uplift, graben formation, and gradual shallowing of the sea are reflected within the bathyal-dominated sedimentary facies of the Kingshill Limestone in St. Croix, U.S. Virgin Islands. Reef-tract limestone beds of Pliocene age were subject to exposure, resubmergence, and meteoric leaching of aragonitic skeletal debris; these beds contain patchy lenses of dolomite that are restricted to a small, structurally-controlled embayment. The South Coast aquifer, the principal water-bearing unit of Puerto Rico's South Coast ground-water province, consists of boulder- to silt-size detritus formed by large and small coalescing fan deltas of Pleistocene to Holocene age. Deep well data indicates that it is possible to vertically separate and group a highly complex and irregular-bedded detrital sequence that underlies distal parts of the fan-delta plain into discrete water-bearing units if correlated with 30- to 40-meter thick, eustatically-controlled depositional cycles. Lithofacies maps show that greatest hydraulic conductivity within the fan-delta plain is generally associated with proximal fan and midfan areas. Distal and interfan areas are least permeable. Alluvial valley aquifers located in the western part of the South Coast ground-water province are important local sources of water supply and appear to contain some of the same physical and hydraulic characteristics as the South Coast aquifer. Older sedimentary rocks within the basin are poor aquifers; conglomeratic beds are well-cemented, and carbonate beds do not contain well-developed solution features, except locally where the beds are overlain by alluvium. Ground-water occurs under unconfined conditions in proximal and midfan areas. Confined conditions within deeper parts of the system and in interfan and some midfan areas are created largely by the intercalated nature of discontinuous fine-grained beds that retard vertical ground-water movement. The development of water resources in southern Puerto Rico has modified the hydrologic system of the South Coast aquifer considerably. Under predevelopment conditions, the South Coast aquifer was recharged in the unconfined, proximal fan and some midfan areas by infrequent rainfall and seepage from streams near the fan apex. Discharge occurred as seabed seepage, baseflow discharge along the lower coastal reach of streams, seepage to coastal wetlands, or evapotranspiration in areas underlain by a shallow water table. Under development conditions, seepage from irrigation canals and areal recharge from furrow irrigation represented a principal mechanism for recharge to the aquifer. Increased ground-water withdrawals in the 1960's and 1970's resulted in declines in the water table to below sea level in some places and intrusion of salt water into the aquifer. By the middle 1980's, a reduction in ground-water withdrawals and a shift from furrow irrigation to drip-irrigation techniques resulted in the recovery of water levels. Under present-day (1986) conditions, regional ground-water flow is coastward but with local movement to some well fields. In addition to the discharge mechanisms described above, ground-water discharges also to coastal canals. The North Coast limestone aquifer system consists of limestone, lesser amounts of dolomite, and minor clastic detritus of Oligocene to Pliocene age that form an unconfined upper aquifer and a confined lower aquifer; these aquifers are separated by a clay, mudstone, and marl confining unit. Topographic relief and incision of carbonate coastal plain rocks by streams are the principal factors controlling the direction of ground-water flow. The North Coast limestone aquifer system is recharged principally by precipitation that enters the upper and lower aquifers where they crop out. Regional groundwater movement from the upper aquifer is to the major rivers, wells, coastal wetlands, coastal, nearshore, and offshore springs, or as seabed seepage. Regional discharge from the lower aquifer is to the major rivers along its unconfined parts or where the confining unit has been breached by streams. Discharge from the lower aquifer also occurs in the San Juan area where the Mucarabones Sand provides an avenue for diffuse upward ground-water flow. Transmissivity within the upper limestone aquifer appears to be largely regulated by the thickness of the freshwater lens. The lens is thickest and transmissivity is greatest in interstream areas that lie in a zone that closely corresponds to the landwardmost extent of the underlying saltwater wedge. Hydraulic conductivity of the upper aquifer generally increases in a coastward direction and reflects lithologic control, karstification in the upper 30 to 100 meters of the section, and enhanced permeability in a zone of freshwater and saltwater mixing. Transmissivity of the lower aquifer is an order of magnitude smaller than that of the upper aquifer; highest transmissivities in the lower aquifer largely correspond to a coarse grainstone-packstone and coral-patch-reef depositional facies contained within the outcropping parts of the Montebello Limestone Member and its subsurface equivalents. Porosity within the North Coast limestone aquifer system is high in grainstone-packstones and low in wackestone and marl. Dolomitized zones and moldic grainstone-packstone strata are the most porous carbonate rocks, but occur in thin beds that usually are only a few meters thick. Processes of karstification that include the development of caverous zones and large vugs, and dissolution along possible regional fracture sets has enhanced permeability within the upper part of the aquifer system. Stratigraphic and lithologic control play an important role controlling permeability within the lower part of the system. The Kingshill aquifer of St. Croix, in large part, is composed of deepwater limestone that contains only microscopic pores and is poorly permeable; however, the upper part of the aquifer, a shallow-water skeletal and reef limestone, is fairly permeable, but restricted in areal extent. Permeability within these uppermost beds of the aquifer has been enhanced by meteoric leaching, dissolution within a mixing zone of saltwater and fresh water, and dolomitization. However, most large-yield wells completed in the Kingshill aquifer are also screened in alluvium that overlies or infills incised channels. The alluvial deposits serve as a temporary storage zone for rainfall, runoff, and ground water slowly entering the Kingshill aquifer.

Puerto Rico↗

USGS2023

A field test of R package GPSeqClus: For establishing animal location clusters

The ability to track animals with Global Positioning System (GPS) collars opened an enormous potential for studying animal movements and behaviour in their natural environment. One such endeavour is to identify clusters of GPS locations as a way to estimate predator kill rate. Clapp et al. (2021) developed an R package ( GPSeqClus ) to assess a location dataset based on user-defined parameters to identify clusters and their characteristics. These characteristics can then help to distinguish resting-site clusters from kill sites of their large (>50 kg) prey. We identified location clusters of an adult male wolf Canis lupus on Ellesmere Island, Nunavut, Canada in July 2009 and tracked him until he died in April 2010. Identifying location clusters was challenging because the collar only obtained two GPS locations per day (12 h apart). In July 2010, we searched 30 of 52 location-clusters we identified as kill/scavenge sites and found 17 of them as such, given they had muskox Ovibos moschatus or caribou Rangifer tarandus pearyi remains nearby. We also documented five wolf rendezvous sites, two den sites, and the wolf's death site to total 60 location-clusters in all. We used a two-step process in testing the R Package GPSeqClus (hereafter GPSeqClus ): (1) compare the number of clusters our method discerned with the number identified by the new algorithm, and (2) compare the number of biologically significant clusters (e.g. den sites, kill/feeding sites) we found with the number the new algorithm located. We made these tests with GPSeqClus by varying the search radius, number of days at a site, and minimum number of locations required for a cluster. GPSeqClus compared well to our technique, with the best sub-algorithm among the 25 we tested only missing three of our identified clusters and yielding six additional clusters. GPSeqClus identified 16 of the 17 confirmed sites of remains, all wolf home sites, and the wolf's carcass site. Identifying clusters using a 500-m search radius, a 1.5-day window, and a minimum of two GPS locations per cluster was suitable for a coarse GPS acquisition rate of two locations per day when prey are large, such as muskox or caribou. Given that GPSeqClus performed well with our coarse location dataset, we expect it will also perform even better with a collar acquiring more than two locations per day. Having a field-tested utility such as GPSeqClus will enhance carnivore predation studies elsewhere.

Ecological Solutions and Evidence↗

USGS2019

Demonstrating the value of Earth observations—methods, practical applications, and solutions—group on Earth observations side event proceedings

Executive Summary The U.S. Geological Survey, the National Oceanic and Atmospheric Administration, the European Association for Remote Sensing Companies, and the European Space Agency in coordination with the GEOValue Community hosted a side event to the Group on Earth Observations Plenary on October 23–24, 2017, in Washington, D.C. The workshop, entitled “Demonstrating the Value of Earth Observations: Methods, Practical Applications and Solutions,” brought together more than 60 international experts including economists, scientists, and engineers to consider the state of the science and applications of valuing Earth observations (EO). This 2-day workshop built upon previous activities developed under the GEOValue initiative. This workshop brought together expert analysts from multiple disciplines and backgrounds who are developing methods to identify and measure the value of information generated from the use of satellite and in-situ data. The mix of government agencies, international financial institutions, and independent consultants who participated in the workshop blended to develop a rich mix of views, approaches, and outcomes. During the first part of the workshop, the focus was on the latest science in valuing EO. A number of methodologies were described. Approaches generally assess the societal benefits of specific actions (for example, investments in EO). Some methods focus on broad measures of economic activity (for example, gross domestic product) or methods to assess total economic value such as contingent valuation surveys. Alternatively, use-case approaches (a use case is defined as an evaluation in which one or more decisions, applications, or other uses of data, information, and information products are specifically considered) start with the specific actions and how information is used to support decision making and affect outcomes. The second part of the meeting was focused on the use and development of value chains and decision trees. A value chain can be defined as the set of value-adding activities that one or more organizations perform in creating and distributing goods and services. In terms of EO, the value chain approach can be applied to consider societal benefits of the data and assess the value of data and data features. The EO value chain considers the geospatial data sources and the processing of the data into value added information to be incorporated into decision-support systems, leading to decision makers’ actions. To understand the value of EO, one would also need to recognize the demand side of the equation or how EO benefits users. Extending the value chain concept and incorporating tenets of Bayesian decision making, a decision tree would include one or more use cases. The value provided by the marginal increase in information could flow from one or several parts of the supply side of the value chain. The decision tree is based on the premise that information has no value if it is not used in at least one decision. By connecting the value chain and the decision tree, a framework is created that allows for conceptualizing the value of EO in its many uses. One can then apply economic techniques to monetize the marginal benefit of an outcome with information versus one without. A third part of the meeting applied the value chain and decision-tree frameworks to five specific thematic areas, each with the focus of using information for a decision point: Effect of increasing temperatures on human health; Flooding—Mitigating, managing, and avoiding impacts to safety and property damage; Harmful algal blooms—Effects on human health, recreation, and tourism; Energy and mineral supply—Mitigating, managing, and avoiding impacts of shortfalls on the economy; and Effects of natural hazards on transportation systems—Effects on mobility, safety, and the economy. During the working session, five separate groups worked to define and delineate the value chains and decision trees associated with each topic, discussing the related challenges and data needs. The outcomes were reported back to the full group. Because of the complexity of the topics, most groups first identified a network of value chains and then narrowed the scope to develop a single value chain to address their group’s topic. Although they worked separately and on different topics, the groups came to similar conclusions, concurring that the value chain and decision-tree frameworks are very effective for informing quantitative impact assessments and developing a relatable narrative to assist the public in understanding the link between EO and citizens.

Open-File Report↗

USGS2014

Hydrology and numerical simulation of groundwater movement and heat transport in Snake Valley and surrounding areas, Juab, Miller, and Beaver Counties, Utah, and White Pine and Lincoln Counties, Nevada

Snake Valley and surrounding areas, along the Utah-Nevada state border, are part of the Great Basin carbonate and alluvial aquifer system. The groundwater system in the study area consists of water in unconsolidated deposits in basins and water in consolidated rock underlying the basins and in the adjacent mountain blocks. Most recharge occurs from precipitation on the mountain blocks and most discharge occurs from the lower altitude basin-fill deposits mainly as evapotranspiration, springflow, and well withdrawals. The Snake Valley area regional groundwater system was simulated using a three-dimensional model incorporating both groundwater flow and heat transport. The model was constructed with MODFLOW-2000, a version of the U.S. Geological Survey’s groundwater flow model, and MT3DMS, a transport model that simulates advection, dispersion, and chemical reactions of solutes or heat in groundwater systems. Observations of groundwater discharge by evapotranspiration, springflow, mountain stream base flow, and well withdrawals; groundwater-level altitudes; and groundwater temperatures were used to calibrate the model. Parameter values estimated by regression analyses were reasonable and within the range of expected values. This study represents one of the first regional modeling efforts to include calibration to groundwater temperature data. The inclusion of temperature observations reduced parameter uncertainty, in some cases quite significantly, over using just water-level altitude and discharge observations. Of the 39 parameters used to simulate horizontal hydraulic conductivity, uncertainty on 11 of these parameters was reduced to one order of magnitude or less. Other significant reductions in parameter uncertainty occurred in parameters representing the vertical anisotropy ratio, drain and river conductance, recharge rates, and well withdrawal rates. The model provides a good representation of the groundwater system. Simulated water-level altitudes range over almost 2,000 meters (m); 98 percent of the simulated values of water-level altitudes in wells are within 30 m of observed water-level altitudes, and 58 percent of them are within 12 m. Nineteen of 20 simulated discharges are within 30 percent of observed discharge. Eighty-one percent of the simulated values of groundwater temperatures in wells are within 2 degrees Celsius (°C) of the observed values, and 55 percent of them are within 0.75 °C. The numerical model represents a more robust quantification of groundwater budget components than previous studies because the model integrates all components of the groundwater budget. The model also incorporates new data including (1) a detailed hydrogeologic framework, and (2) more observations, including several new water-level altitudes throughout the study area, several new measurements of spring discharge within Snake Valley which had not previously been monitored, and groundwater temperature data. Uncertainty in the estimates of subsurface flow are less than those of previous studies because the model balanced recharge and discharge across the entire simulated area, not just in each hydrographic area, and because of the large dataset of observations (water-level altitudes, discharge, and temperatures) used to calibrate the model and the resulting transmissivity distribution. Groundwater recharge from precipitation and unconsumed irrigation in Snake Valley is 160,000 acre-feet per year (acre-ft/yr), which is within the range of previous estimates. Subsurface inflow from southern Spring Valley to southern Snake Valley is 13,000 acre-ft/yr and is within the range of previous estimates; subsurface inflow from Spring Valley to Snake Valley north of the Snake Range, however, is only 2,200 acre-ft/yr, which is much less than has been previously estimated. Groundwater discharge from groundwater evapotranspiration and springs is 100,000 acre-ft/yr, and discharge to mountain streams is 3,300 acre-ft/yr; these are within the range of previous estimates. Current well withdrawals are 28,000 acre-ft/yr. Subsurface outflow from Snake Valley moves into Pine Valley (2,000 acre-ft/yr), Wah Wah Valley (23 acre-ft/yr), Tule Valley (33,000 acre-ft/yr), Fish Springs Flat (790 acre-ft/yr), and outside of the study area towards Great Salt Lake Desert (8,400 acre-ft/yr); these outflows, totaling about 44,000 acre-ft/yr, are within the range of previous estimates. The subsurface flow amounts indicate the degree of connectivity between hydrographic areas within the study area. The simulated transmissivity and locations of natural discharge, however, provide a better estimate of the effect of groundwater withdrawals on groundwater resources than does the amount and direction of subsurface flow between hydrographic areas. The distribution of simulated transmissivity throughout the study area includes many areas of high transmissivity within and between hydrographic areas. Increased well withdrawals within these high transmissivity areas will likely affect a large part of the study area, resulting in declining groundwater levels, as well as leading to a decrease in natural discharge to springs and evapotranspiration.

Nevada, Utah↗

USGS2012

Hydrogeology, water chemistry, and transport processes in the zone of contribution of a public-supply well in Albuquerque, New Mexico, 2007-9

The National Water-Quality Assessment Program (NAWQA) of the U.S. Geological Survey began a series of groundwater studies in 2001 in representative aquifers across the Nation in order to increase understanding of the factors that affect transport of anthropogenic and natural contaminants (TANC) to public-supply wells. One of 10 regional-scale TANC studies was conducted in the Middle Rio Grande Basin (MRGB) in New Mexico, where a more detailed local-scale study subsequently investigated the hydrogeology, water chemistry, and factors affecting the transport of contaminants in the zone of contribution of one 363-meter (m) deep public-supply well in Albuquerque. During 2007 through 2009, samples were collected for the local-scale study from 22 monitoring wells and 3 public-supply (supply) wells for analysis of major and trace elements, arsenic speciation, nutrients, dissolved organic carbon, volatile organic compounds (VOCs), dissolved gases, stable isotopes, and tracers of young and old water. To study groundwater chemistry and ages at various depths within the aquifer, the monitoring wells were divided into three categories: (1) each shallow well was screened across the water table or had a screen midpoint within 18.3 m of the water level in the well; (2) each intermediate well had a screen midpoint between about 27.1 and 79.6 m below the water level in the well; and (3) each deep well had a screen midpoint about 185 m or more below the water level in the well. The 24-square-kilometer study area surrounding the "studied supply well" (SSW), one of the three supply wells, consists of primarily urban land within the MRGB, a deep alluvial basin with an aquifer composed of unconsolidated to moderately consolidated deposits of sand, gravel, silt, and clay. Conditions generally are unconfined, but are semiconfined at depth. Groundwater withdrawals for public supply have substantially changed the primary direction of flow from northeast to southwest under predevelopment conditions, to west to east under modern conditions. Analysis of age tracers indicates that groundwater from most sampled wells is dominated by old (pre-1950) water, ranging in mean age from about 4,000 years to more than 22,000 years, but includes a fraction of young (post-1950) recharge. Patterns in chemical and isotopic data are consistent with the conclusions that shallow groundwater in the area typically includes a fraction that evaporated prior to recharge and (or) flushed accumulated solutes out of the unsaturated zone during recharge, and that shallow groundwater has mixed to deeper parts of the aquifer, which receives recharge mainly by seepage from the Rio Grande. Among shallow and intermediate wells that produced water with a fraction of young recharge, that fraction ranged between 1.5 and 46 percent. Samples from the two deep wells had groundwater ages exceeding 18,000 years, with no fraction of young recharge. Two supply wells (including the SSW) had a fraction of young recharge, which ranged between about 3 and 11 percent, despite mean groundwater ages exceeding 10,000 years. The fraction of young recharge to the SSW varied seasonally, probably because seasonal pumping patterns affected local hydraulic gradients and (or) because of flow through the well bore when the SSW is not pumping. Well-bore flow data collected during winter (low-pumping season) indicated that about 61 percent of the water pumped from the SSW entered the well from the intermediate part of the aquifer, and that the remaining 39 percent entered from the deep part of the aquifer. Volatile organic compounds (VOCs) were detected in samples from most shallow and intermediate monitoring wells and from two of three supply wells, including the SSW. Detected VOCs were primarily chlorinated solvents or their degradation products. Many of the wells in which most of these VOCs were detected are located near known sites of solvent contamination that were targeted for sampling because trichloroethylene (TCE) and cis-1,2-dichloroethylene had been detected in the SSW, and several of these wells may have become contaminated at least partly because of enhanced vertical migration associated with the pumping of and (or) direct migration down deep well bores. Except for TCE in the sample from a shallow monitoring well, all detections of VOCs were at concentrations below Maximum Contaminant Levels (MCLs) set by the U.S. Environmental Protection Agency. Concentrations of all VOCs detected in the supply wells were less than one-tenth of the corresponding MCLs. However, the presence of VOCs in all but deep groundwater, including the detection of chloroform (a chlorination byproduct) in several shallow wells, indicates that groundwater in the study area commonly is affected by human activities, even to substantial depths. The only natural contaminant detected at concentrations near or above its MCL was arsenic, which has been detected at elevated concentrations across broad areas of the MRGB. Concentrations of arsenic, present primarily as arsenate, exceeded the MCL of 10 micrograms per liter (μg/L) in water from the two deep wells (one of which had the highest concentration, 35 μg/L), from one intermediate well, and from two supply wells, including the SSW. Water-quality and solid-phase data from this study are consistent with elevated arsenic concentrations in groundwater being related to pH-dependent desorption of arsenic from ferric oxyhydroxides in sediments in deep parts of the aquifer. Concentrations of nitrate ranged between 1.3 and 5.4 milligrams per liter (mg/L) in water from shallow wells screened across the water table, but were less than 0.9 mg/L in water from all but one deeper well. Nitrogen isotopes and chloride/bromide ratios for shallow wells were consistent with natural soil nitrogen. Nitrate concentrations and nitrogen isotopes indicated that denitrification is occurring at intermediate aquifer depths, and that the progress of the denitrification reaction typically is greatest for wells that include a fraction of groundwater associated with particular recharge sources or with known sites of contamination contributing organic compounds that can provide a carbon source for microbial respiration. Overall, hydrologic and chemical data from the study area indicate that young recharge is reaching the aquifer across broad areas and is migrating from shallow to intermediate depths of the aquifer as a result of mixing that is associated with human development of groundwater. Consequently, groundwater that human activities in the urban study area have affected is present at depths that are within the screened intervals of public-supply wells, resulting in detections of VOCs and implying greater vulnerability to anthropogenic contamination than might be assumed based on the dominantly old age of the regional groundwater. However, the fractions of old groundwater that public-supply wells produce substantially dilute the anthropogenic contaminants, while contributing natural contaminants (primarily arsenic) to the wells. Based on data from the SSW, vulnerability of public-supply wells to natural and anthropogenic contaminants in the area changes through time, including with seasonal changes in pumping stresses that alter the fractions of young and old water being contributed to wells.

New Mexico↗