SEARCH · Geology Reports

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.

250 records · Page 14Linked to original sources

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↗

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↗

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↗

USGS2012

Hydrological information products for the Off-Project Water Program of the Klamath Basin Restoration Agreement

The Klamath Basin Restoration Agreement (KBRA) was developed by a diverse group of stakeholders, Federal and State resource management agencies, Tribal representatives, and interest groups to provide a comprehensive solution to ecological and water-supply issues in the Klamath Basin. The Off-Project Water Program (OPWP), one component of the KBRA, has as one of its purposes to permanently provide an additional 30,000 acre-feet of water per year on an average annual basis to Upper Klamath Lake through "voluntary retirement of water rights or water uses or other means as agreed to by the Klamath Tribes, to improve fisheries habitat and also provide for stability of irrigation water deliveries." The geographic area where the water rights could be retired encompasses approximately 1,900 square miles. The OPWP area is defined as including the Sprague River drainage, the Sycan River drainage downstream of Sycan Marsh, the Wood River drainage, and the Williamson River drainage from Kirk Reef at the southern end of Klamath Marsh downstream to the confluence with the Sprague River. Extensive, broad, flat, poorly drained uplands, valleys, and wetlands characterize much of the study area. Irrigation is almost entirely used for pasture. To assist parties involved with decisionmaking and implementation of the OPWP, the U.S. Geological Survey (USGS), in cooperation with the Klamath Tribes and other stakeholders, created five hydrological information products. These products include GIS digital maps and datasets containing spatial information on evapotranspiration, subirrigation indicators, water rights, subbasin streamflow statistics, and return-flow indicators. The evapotranspiration (ET) datasets were created under contract for this study by Evapotranspiration, Plus, LLC, of Twin Falls, Idaho. A high-resolution remote sensing technique known as Mapping Evapotranspiration at High Resolution and Internalized Calibration (METRIC) was used to create estimates of the spatial distribution of ET. The METRIC technique uses thermal infrared Landsat imagery to quantify actual evapotranspiration at a 30-meter resolution that can be related to individual irrigated fields. Because evaporation uses heat energy, ground surfaces with large ET rates are left cooler as a result of ET than ground surfaces that have less ET. As a consequence, irrigated fields appear in the Landsat images as cooler than nonirrigated fields. Products produced from this study include total seasonal and total monthly (April-October) actual evapotranspiration maps for 2004 (a dry year) and 2006 (a wet year). Maps showing indicators of natural subirrigation were also provided by this study. "Subirrigation" as used here is the evapotranspiration of shallow groundwater by plants with roots that penetrate to or near the water table. Subirrigation often occurs at locations where the water table is at or above the plant rooting depth. Natural consumptive use by plants diminishes the benefit of retiring water rights in subirrigated areas. Some agricultural production may be possible, however, on subirrigated lands for which water rights are retired. Because of the difficulty in precisely mapping and quantifying subirrigation, this study presents several sources of spatially mapped data that can be used as indicators of higher subirrigation probability. These include the floodplain boundaries defined by stream geomorphology, water-table depth defined in Natural Resources Conservation Service (NRCS) soil surveys, and soil rooting depth defined in NRCS soil surveys. The two water-rights mapping products created in the study were "points of diversion" (POD) and "place of use" (POU) for surface-water irrigation rights. To create these maps, all surface-water rights data, decrees, certificates, permits, and unadjudicated claims within the entire 1,900 square mile study area were aggregated into a common GIS geodatabase. Surface-water irrigation rights within a 5-mile buffer of the study area were then selected and identified. The POU area was then totaled by water right for primary and supplemental water rights. The maximum annual volume (acre-feet) allowed under each water right also was calculated using the POU area and duty (allowable annual irrigation application in feet). In cases where a water right has more than one designated POD, the total volume for the water right was equally distributed to each POD listed for the water right. Because of this, mapped distribution of diversion rates for some rights may differ from actual practice. Water-right information in the map products was from digital datasets obtained from the Oregon Water Resources Department and was, at the time acquired, the best available compilation of water-right information available. Because the completeness and accuracy of the water-right data could not be verified, users are encouraged to check directly with the Oregon Water Resources Department where specific information on individual rights or locations is essential. A dataset containing streamflow statistics for 72 subbasins in the study area was created for the study area. The statistics include annual flow durations (5-, 10-, 25-, 50-, and 95-percent exceedances) and 7-day, 10-year (7Q10) and 7-day, 2-year (7Q2) low flows, and were computed using regional regression equations based on measured streamflow records in the region. Daily streamflow records used were adjusted as needed for crop consumptive use; therefore the statistics represent streamflow under more natural conditions as though irrigation diversions did not exist. Statistics are provided for flow rates resulting from streamflow originating from within the entire drainage area upstream of the subbasin pour point (referring to the outlet of the contributing drainage basin). The statistics were computed for the purpose of providing decision makers with the ability to estimate streamflow that would be expected after water conservation techniques have been implemented or a water right has been retired. A final product from the study are datasets of indicators of the potential for subsurface return flow of irrigation water from agricultural areas to nearby streams. The datasets contain information on factors such as proximity to surface-water features, geomorphic floodplain characteristics, and depth to water. The digital data, metadata, and example illustrations for the datasets described in this report are available on-line from the USGS Water Resources National Spatial Data Infrastructure (NSDI) Node Website http://water.usgs.gov/lookup/getgislist or from the U.S. Government website DATA.gov at http://www.data.gov with links provided in a Microsoft® Excel® workbook in appendix A.

California;Oregon↗

USGS2017

Geology and mining history of the Southeast Missouri Barite District and the Valles Mines, Washington, Jefferson, and St. Francois Counties, Missouri

The Southeast Missouri Barite District and the Valles Mines are located in Washington, Jefferson, and St. Francois Counties, Missouri, where barite and lead ore are present together in surficial and near-surface deposits. Lead mining in the area began in the early 1700’s and extended into the early 1900’s. Hand mining of lead in the residuum resulted in widespread pits (also called shafts or diggings), and there was some underground mining of lead in bedrock. By the 1860’s barite was recovered from the residuum by hand mining, also resulting in widespread diggings, but generally not underground mines in bedrock. Mechanized open-pit mining of the residuum for barite began in the 1920’s. Barite production slowed by the 1980’s, and there has not been any barite mining since 1998. Mechanized barite mining resulted in large mined areas and tailings ponds containing waste from barite mills. The U.S. Environmental Protection Agency (EPA) has determined that lead is present in surface soils in Washington and Jefferson Counties at concentrations exceeding health-based screening levels. Also, elevated concentrations of barium, arsenic, and cadmium have been identified in surface soils, and lead concentrations exceeding the Federal drinking-water standard of 15 micrograms per liter have been identified in private drinking-water wells. Potential sources of these contaminants are wastes associated with barite mining, wastes associated with lead mining, or unmined natural deposits of barium, lead, and other metals. As a first step in helping EPA determine the source of soil and groundwater contamination, the U.S. Geological Survey (USGS), in cooperation with the EPA, investigated the geology and mining history of the Southeast Missouri Barite District and the Valles Mines. Ore minerals are barite (barium sulfate), galena (lead sulfide), cerussite (lead carbonate), anglesite (lead sulfate), sphalerite (zinc sulfide), smithsonite (zinc carbonate), and chalcopyrite (copper-iron sulfide). The Cambrian Potosi Dolomite is the most important formation for the ore deposits, followed by the Eminence Dolomite. Because galena, sphalerite, and barite are less soluble than dolomite, chemical weathering of the ore-bearing dolomite bedrock resulted in the concentration of ore minerals in the residuum. Most of the barite and lead mining was in the residuum, which averages 10 to 15 feet thick. Lead mining by French explorers may have begun in 1719 along Old Mines Creek at Cabanage de Renaudiere, which was followed shortly by the discovery of lead and the development of lead mines at Mine Renault (also called Forche a Renault Mine), Old Mines, and at other places along the Big River, Mineral Fork, and Forche a Renault Creek. Lead mining began sometime between 1775 and 1780 at Mine a Breton, the name of which was later changed to Potosi. Other mining areas were developed in the early part of the 19th century, including Fourche a Courtois (Palmer Mines), the French Diggings, and the Richwoods Mines. Zinc became a valuable resource after the Civil War, and the Valles Mines was an important supplier of zinc as well as lead, with at least some production up until the 1920’s. Lead mining declined in the early part of the 20th century as mining in the Old Lead Belt, Mine La Motte, and the Tri-State District expanded. The earliest lead mines were diggings in the residuum and were round holes (shafts) about 4 feet in diameter dug with pick and shovel about 15–20 feet deep, with drifts dug a short distance laterally from the bottom of the shafts. This mining process was repeated a short distance away until a large area was covered with pits. Some mining in bedrock began by about 1800, with shafts as deep as 170 feet and as much as several hundred feet of lateral drifts. Smelting of the lead ore to elemental lead was first done using a log furnace, which was inefficient; estimates have been made that only about 50 percent of the lead was recovered, and the remainder was lost to the ashes (slags) and to volatilization. Starting in 1798, ash furnaces were used to smelt the ashes from the log furnaces. These two furnaces were worked in tandem for many years but were gradually replaced by other furnaces, including the Scotch hearth. Estimates of lead recovery as high as 80–90 percent have been made for the Scotch hearth. By the mid-1870’s the air furnace was being used, also with estimated lead recovery as high as 80–90 percent. Zinc furnaces were built when zinc became a valuable commodity, but much of the zinc ore was shipped out of the area, either to a smelter in St. Louis, Missouri, or to other smelters. The total lead and zinc production from the Southeast Missouri Barite District and the Valles Mines is estimated at 180,000 tons of lead and 60,000 tons of zinc. An estimated 97,000 tons of lead and an estimated 120,000 tons of zinc were lost during smelting. The estimated losses do not include losses at the mine site during mining and preparation for smelting, such as the loss of fine-grained galena during hand cleaning or the discarding of zinc ore before its value was known, for which no estimates are available. Hand mining for barite in the residuum was active by at least the 1860’s and peaked from 1905 to the 1930’s when several thousand people were engaged in barite mining. Hand mining (diggings) and cleaning of the ore was done in much the same way as earlier lead mining, with the additional use of a rattle box to further clean the barite. Mechanized open-pit mining of old barite diggings began in 1924 to recover barite left behind by hand mining, and washing plants were used to clean the clay from the barite. Hand mining, however, continued to thrive, and washer plants began to close temporarily in 1931; nearly all of the barite produced before 1937 was by hand mining. By the 1940’s, however, all barite mining was mechanized. Mechanized mining used shovels powered by steam, gasoline, or electricity (and by the 1950’s draglines and front-end loaders) to mine the residuum. The ore was loaded onto rail cars (and by the 1940’s, trucks) for shipment to washer plants. Clay was removed from the barite using a log washer, and a jig was used to concentrate the barite. Overflow from the log washers was waste and went to a mud (tailings) pond. The coarse jig tailings went to tailings piles or were used as railroad ballast and, later, to create roads within the mine pit. Some barite was ground, depending on its final use, and some ground barite was bleached using a hot solution of sulfuric acid to remove impurities such as iron minerals and lead sulfide (galena). An earlier bleaching process used lead-lined tanks. Large quantities of water were required for milling the barite; some was recirculated water and the remainder came from dammed streams or was pumped from wells. Tailings and wastewater were impounded behind dikes that were built across small valleys and were increased in height as necessary using washer waste and any overburden that had been stripped. In some cases, dikes were built across valleys that had already been mined for barite. The total production of barite from the Southeast Missouri Barite District and the Valles Mines is estimated to have been about 13.1 million tons. Most of the barite production was from Washington County. Hand mining and processing of barite was inefficient. Estimates of barite recovery range from less than one-fourth to about one-half because pillars between the shafts in the residuum needed to be left unmined for stability. With mechanized mining, large amounts of barite were lost during the milling process. It has been estimated that about 30 percent of the barite was lost and that about two-thirds of the lost barite was fine-grained and was discharged to the tailings ponds. Some galena was lost to the tailings ponds. A 1972 inventory of tailings ponds by the Missouri Geological Survey identified 67 ponds in the Southeast Missouri Barite District (there are more than this currently documented). Results from samples from four ponds that were drilled were used to estimate that the 67 ponds contained almost 39 million tons (or cubic yards) of tailings averaging about 5 percent barite, for a potential reserve of 1.935 million tons of barite. It is not known how much lead was removed during barite mining, either by hand or mechanized mining and processing, how much lead was recovered, or how much lead went as fines to the tailing ponds or as coarse material to mine roads or was otherwise lost.

Missouri↗

USGS2022

Fishway Entrance Palisade

This technical report summarizes the work that was conducted by the University of Massachusetts Amherst and the United States Geological Survey (USGS), along with other project partners, on the Fishway Entrance Palisade (EP), a projected funded through the Department of Energy’s (DOE) funding opportunity titled ‘Innovative Solutions for Fish Passage at Hydropower Dams’ (DE‐FOA‐0001662). The period of performance ranged from September 1, 2018 through September 30, 2021. The EP is a novel fish passage engineering technology designed to provide more favorable entry conditions for fish and to reduce costs relative to conventional fishway auxiliary water systems (AWS). The EP project has four primary components. First, the Northeast United States Auxiliary Water Systems Database was created (Northeast Fishway Auxiliary Water Systems Database Section). The database, developed with material provided by the U.S. Fish and Wildlife Service, contains information on fishway type (e.g., lift, Denil, pool and weir) and Auxiliary Water System (AWS) details (e.g., water conveyance method, diffuser type) for 60 hydroelectric sites in the region. Findings indicate that nearly 4 out of every 10 fishway in the region is a fish lift and approximately 1 out of every 4 is a Denil ladder. The remainder are a mix of vertical slot fishways, pool and weirs, and Ice Harbor fishways. Furthermore, over half of all AWS systems use floor diffusers to discharge the auxiliary (or attraction) water into the entrance of a fishway, whereas only 14% use wall diffusers. Second, limited experiments on a conventional AWS with live, actively migrating fish were conducted at the USGS Easter Ecological Science Center (EESC) S.O. Conte Research Laboratory (Conventional Auxiliary Water System Experiments Section). This study determined how water velocity through a wall diffuser, without turning vanes or timber baffles to distribute the flow, affects the behavior and passage of adult American shad, a conservative surrogate species for migratory fish on the East Coast. Two gross diffuser velocity treatments were examined, 0.5 ft/s and 1.0 ft/s. These wall diffuser velocities represented current (0.5 ft/s) and past (1.0 ft/s) design criteria guidelines set forth by the USFWS North Atlantic-Appalachian Region (Rojas 2020; USFWS 2019). Six trials with a total of 151 American Shad were conducted in June of 2019 for the two treatments. No differences in American shad passage efficiency were discovered between the two treatments, while approximately 3 in every 4 attempts were successful at passing the diffuser. While these results may appear to indicate that the generally accepted gross wall diffuser velocity criteria for American shad of 0.5 ft/s could be safely increased to 1.0 ft/s, further analysis is warranted. Furthermore, it is unknown how other migratory and resident fish species that traverse these structures would be impacted by such a change. Studying the wall diffuser hydraulics led to an important AWS observation. Without turning vanes or timber baffles in this study, doubling the diffuser area was insufficient at producing the type of flow field change one may expect by halving the gross diffuser velocity. Instead, the flow fields throughout each treatments study area were similar, which led to similar results in shad performance. This not only highlights the importance of installing flow guidance devices like turning vanes, but also to the importance of properly maintaining them, which can be costly. Third, more expansive experiments on the novel EP were conducted in the spring of 2019 and 2021 (Fishway Entrance Palisade Experiments). The goal of this study was to determine how adult American shad responded to a variety of conditions at a full-scale EP. A total of six treatments were examined by changing the average auxiliary channel velocity between 1.0 and 5.0 ft/s in intervals of 1.0 ft/s and by inserting/removing an entrance gate at the opening of the fishway. Thirty trials with a total of 1,273 shad were conducted over the two years. In all treatments, at least ~7 out of every 10 fish successfully passed the EP diffuser and swam into the entrance channel within the 3.5-hour long trial, highlighting the general effectiveness of the novel AWS technology. In both study years, lower velocities through the EP diffuser led to increased shad performance, though performance peaked for the 2 ft/s velocity treatment. This treatment condition represents an approximate six-fold increase in gross diffuser velocity relative to conventional auxiliary water systems, which in turn presents opportunities for cost savings (e.g., reduction in diffuser size). Shad performance, in general, was worse in 2019 than in 2021, potentially due to the different run timing when our trials were conducted (2019 trials occurred near the end of the migration season, unlike in 2021). Treatments in 2019 had approximately a 20% reduction in entrance efficiency by the trial end, including a 16.7% drop for the 3 ft/s velocity treatment in 2019 relative to 2021 (the only carryover treatment between years). Lastly, adding an entrance gate caused a significant delay to entry. The time to 25% entry raised ~20 minutes from the near instantaneous 25% entry that was reported for the other treatments conducted in the same year (2021). Though by the end of the 3.5-hour trial, the overall entrance efficiency nearly matched those of the other 2021 treatments. The fourth and final component of the EP project was an economic analysis that focused on the cost of attraction and environmental flows (Modeling Power Generation Losses Due to Environmental and Fish Passage Attraction Flows at a Run-Of-River Hydroelectric Operation in the Northeast). The study assessed the economic impact of meeting environmental flow requirements at a representative hydroelectric facility and fish lift in the Northeast. An initial finding of the study was that there is a paucity of published data on the costs of meeting attraction and environmental flows. This is due, in part, to the proprietary nature of this data. To explore the costs associated with these flows, three types of environmental flows were assessed: upstream fishway attraction flows, downstream fishway attraction flows, and habitat maintenance minimum flows. A physics-based model was developed and calibrated with three years of hourly generation and flow data as inputs. Gage flow inputs were adjusted and used to calculate power generated. To address hydrologic variability, the model was executed to simulate 30 years of historical flows. Results indicate that both interannual and seasonal climatic factors impact the costs of meeting environmental flow requirements. Generation potential is most strongly curtailed during dry years in terms of maximizing the capacity factor (the percent of time a plant generates at capacity). Dry years, and especially dry summers, have the most significant costs associated with mitigation flows. Of the three types of flows, habitat flows are most costly in terms of power production, followed by upstream attraction flows. Downstream attraction flows are least costly. This finding is the likely result of differences in both flow rates and duration of the seasonal requirement for each flow. Overall, environmental flows represented a 2-12% loss in annual generation, but losses during a dry summer can reach over 20%.

Final Technical Report↗

USGS1967

Swatara Creek basin of southeastern Pennsylvania: An evaluation of its hydrologic system

Local concentrations of population in the Swatara Creek basin of Pennsylvania find it necessary to store, transport, and treat water because local supplies are either deficient or have been contaminated by disposal of wastes in upstream areas. Water in the basin is available for the deficient areas and for dilution of the coal-mine drainage in the northern parts and the sewage wastes in the southern parts. Swatara Creek drains 576 square miles just east of Harrisburg, Pa., and is the largest tributary to the Susquehanna River from the north side below Harrisburg. It rises in the southern Pocono Mountains and flows southwestward across the Lebanon Plateau. On an average day Swatara Creek discharges more than 630 million gallons into the Susquehanna River at Middletown, Pa. In a year this amounts to about 23 inches of water over the entire basin and is the residual from an average annual precipitation of 45.5 inches. During an average year the flow in Swatara Creek from the upper third of the basin above Harper Tavern is always greater than 1,300 mgd (million gallons per day) for at least 15 days and is always greater than 25 mgd for at least 350 days. The daily streamflow from the basin averages 1.1 mgd per sq mi, but yields from different areas range from 0.97 to 1.22 mgd per sq mi. These variations are caused chiefly by differences in precipitation and land cover. The area of lowest yield is in the valleys west of Tremont, and the highest yields are in the Upper and Lower Little Swatara Creek subbasins. At high and medium stages the chemical character of the water in the streams is suitable for public and private supplies. At lower stages, defending on the areas and the amounts of contamination by coal-mine drainage and sewage pollution, the natural flow may require some treatment. At low stages the chemical characteristics of the natural flow not affected by man is almost identical with that of the ground water in the area drained by the stream. In general, the total dissolved solids range from about 25 to 400 parts per million and the hardness is as much as about 300 parts per million. The ground-water increment to the base flow of Swatara Creek averages about 240 mgd, or about 8.8 inches annually, for the basin. Generally, ground-water supplies in amounts of less than 0.5 mgd can be developed south of Blue Mountain. Supplies of several million gallons per day have been developed for industrial use from the permeable limestones in the south-central part of the basin. More intensive investigation in other parts of the basin would indicate areas where supplies of more than 0.5 mgd could be developed from properly spaced wells. The chemical character of water from wells depends largely on the host rock. In highly soluble rocks water contains large amount of dissolved solids; in more resistant rocks concentrations are lower. The chemical character of unpolluted ground water generally reflects the composition of the more readily soluble minerals in the local geologic environment. Areas contaminated by septic- tank effluent may have above normal amounts of nitrate and detergent products. Except where polluted, most ground water is suitable for public and industrial uses without extensive treatment. Sites for storage of surface water exist in the part of the basin lying in the valley and ridge area. As much as 30 to 40 percent of the annual flow could be impounded for release as low-flow augmentation for dilution of mine drainage and other wastes in the basin. Low sediment yields of supplying drainage areas would ensure a long life expectancy of reservoirs at these sites. Overbank flooding of the main stem of the Swatara Creek and its tributaries has occurred many times in the past. However, it has not been a hazard because urban development has not encroached on the flood plain. An inundation map of the August 1933 flood provides a basis that urban planners may use to avoid future damage. As water in the Swatara Creek moves downstream to the Susquehanna River, the flow is influenced consecutively by a large annual rainfall on the northern valley and ridge area, the wastes of surface and subsurface coal-mining activities, and less annual rainfall on the part of the basin lying in the Lebanon Plateau area; the flow is supplemented and further influenced by many tributaries and by the industrial and domestic wastes that are carried by these secondary streams. The annual precipitation ranges from 52 inches at the east edge and 49 inches at the west edge of the mountainous part of the basin to about 41 inches at the southwestern part at Middletown. The rainfall generally is adequate during the growing season to mature the crops. The mean annual temperature at Lebanon is about 52°F, and the growing season is about 180 days. In this report the basin has been divided into eight hydrologic zones, leased on runoff, natural use of water, and chemical character of water. Four zones lie in the valley and ridge area, three lie in the Lebanon Plateau area, and one lies in the highland along the southeastern basin boundary. In each of the zones the hydrologic characteristics are virtually the same, but they may be completely different from those in adjacent zones. The boundaries of the zones generally coincide with boundaries between geologic formations, and the areas in each zone include rocks of similar influence on water. Streams in zone 4 at the northeast edge of the plateau have the highest average surface runoff from 1.2 to 1.1 mgd per sq mi whereas those in zone 2 at the northwest edge of the valley and ridge area have the lowest, about 1.0 mgd. Streams in zone 8, along the southeast edge of the basin, have the largest sustained low-flow yield, about 0.26 to 0.19 mgd per sq mi; those in zone 5 overlying the Martinsburg Shale east of Harrisburg have the smallest sustained low-flow yields, 0.03 to 0.01 mgd. Streams in the limestone area of zone 7 have the greatest range in low-flow yields in any one zone from 0.60 to 0 mgd per sq mi. Low-flow yields in zones 1 through 4 range from 0.13 to 0.03 mgd per sq mi. Surface flows from zones 1 and 2 are generally acidic and contain high concentrations of sulfate, iron, and total dissolved solids especially where contaminated with mine wastes. Surface flows from zones 3 and 4 are dilute, slightly alkaline, and suitable for public water supplies. Surface flows from zones 5, 6, and 7 are alkaline and contain moderate concentrations of dissolved solids with waters of highest hardness occurring in zone 7. Surface flows from zone 8 are dilute to moderately mineralized and are relatively high in silica concentration. Nitrate concentrations are high in surf Fee flows below sewage outfalls and in ground water contaminated by septic tank effluent and industrial wastes. Average annual sediment yields of 550 to 650 tons per square mile are characteristic of zones 1 and 2 where strip mining has destroyed the forest cover and coal culm is carried into the streams. From agricultural lands on the Martinsburg Shale in zones 5 and 6, annual sediment yields range from 300 to 350 tons per square mile; but from agricultural lands on the siliceous rocks in zone 8 and zones 3 and 4 in the valley and ridge area, the sediment yield ranges from 200 to 250 tons annually per square mile. Lowest annual sediment yields in the basin are in the forested areas of siliceous rocks in zones 2, 3, 4, and 5, and in the sinkhole topography of the limestones in zone 7 where the yield ranges from 30 to 35 tons and 50 to 60 tons per square mile, respectively. The amount of ground water that can be developed in the basin is dependent on the ability of the underlying rocks to yield water to wells. More than 300 gpm (gallons per minute) can be obtained from wells in alluvial materials in the valley bottoms and in some of the limestones where large solution channels and fractures are penetrated by the wells. From 50 to 300 gpm can be obtained from wells in loosely cemented sandstones and in fractured limestones. From 10 to 50 gpm can be developed from wells in the shales and harder sandstones. The most dense rocks will yield from 1 to 10 gpm from fractures and crevices. Most wells yield water from the upper 350 feet of the formation, for this part contains the most fractures or solution channels. Studies show that the velocity at which a contaminant will move downstream in the basin is related to the discharge of the stream at the time. At a stream discharge of about 400 mgd at Pine Grove, a contaminant in Swatara Creek would require about 40 hours to move from Pine Grove to Middletown. As a result of dispersion and dilution, the maximum concentration of the contaminant at Middletown would be less than 20 percent the concentration at Pine Grove under these conditions. An evaluation of the availability of water in the basin indicates that about I,239 mgd enters as precipitation, 630 mgd leaves as streamflow, 580 mgd is evaporated and transpired, and 56 mgd is diverted for use by man. Not all the diversions for man's use are lost to the basin, as about 27 mgd is returned as sewage for reuse. About one-fourth of the waste water is returned to the ground and the remainder to stream drainageways. Of that diverted by man, 11.6 mgd is used for public supply and 44.4 mgd for industrial and private supplies. Diversions of streamflow furnish 86 percent of the public supply and 27 percent of the industrial supply, and ground-water sources yield the remainder. Municipal and private sewage treatment plants are upgrading the waste water in many places, but no provisions are being made for treatment other than natural dilution and assimilation for the 15 mgd of coal-nine drainage in the northern part of the basin. Technology for economic treatment of mine water is not available at this time, although research in this field is being done. Urbanization eastward from Harrisburg and around Lebanon has increased the population density of the basin. Densities of 500 people per square mile and water use exceeding 2.0 mgd per sq mi can be expected in the future. By the year 2000 the population of the basin may increase 60 percent; and if the per capita rate of use increases 0.5 percent per year the domestic requirements for water will be about two times the present use, or 23 mgd. Similarly, if the present 1:4 ratio of domestic use to industrial use of water continues, at least 89 mgd will be needed for industry in the future. Although an increase to twice the present use of water can be foreseen, or 112 mgd, water for the dilution and assimilation of wastes from treatment systems are not included. Providing water for dilution of wastes from treatment plants has not been a problem, but in the future the amounts needed for this purpose will be greater as the population increases. As water becomes more valuable, treatment of sewage wastes to reduce the biochemical-oxygen-demand load by at least 80 to 90 percent will be necessary to conserve water for more productive uses. As much as 100 mgd may be needed for waste dilution in the basin by year 2000. The present trends in suburban and light industrial development will probably persist in the basin. Problems arising through changes in economic value of water, conflicts in use, and alternatives in development are typical of those confronting the manager of a water-resource system.

Pennsylvania↗