Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Natural Resource Modeling”

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

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

At least 1,603 records · Page 89Linked to original sources
USGS2021

Occurrence and distribution of PFAS in sampled source water of public drinking-water supplies in the surficial aquifer in Delaware, 2018; PFAS and groundwater age-dating results

The U.S. Geological Survey, in cooperation with the Delaware Department of Natural Resources and Environmental Control and the Delaware Geological Survey, conducted a groundwater-quality investigation to (1) describe the occurrence and distribution of PFAS, and (2) document any changes in groundwater quality in the Columbia aquifer public water-supply wells in the Delaware Coastal Plain between 2000 and 2008 and between 2008 and 2018. Thirty public water-supply wells located throughout the Columbia aquifer of the Delaware Coastal Plain were sampled from August through November 2018. Groundwater collected from the wells was analyzed for the occurrence and distribution of 18 per- and polyfluorinated alkyl substances (PFAS) as well as groundwater age. Descriptive statistical analyses were performed to assess PFAS analytical results within the well network and the combined perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS) concentrations were compared to the U.S. Environmental Protection Agency’s (EPA) health advisory level (HAL) for informational purposes only and not for evidence of compliance or noncompliance with Federal regulations. The EPA’s HAL is a health-based reference level for public drinking water as supplied to customers and is not applied to source (raw) water. Groundwater-age data were compared for sites sampled in 2000, 2008, and 2018 to document any changes. All samples were analyzed for 18 PFAS using EPA Method 537 (modified). Forty-four percent of the analyzed PFAS were detected in the study well network. Sixteen of the sampled wells have one or more PFAS detections, and as many as eight different PFAS were found in a single sample. Wells with a higher number of PFAS detected (five or more) were in New Castle and Sussex Counties. The PFAS most frequently detected were PFOA, with 47 percent detection; perfluorohexanoic acid (PFHxA), with 33 percent detection; and PFOS and perfluorohexane sulfonate (PFHxS), with 27 percent detection each. PFAS concentrations were below 1,000 parts per trillion (ppt). Two wells exceeded the EPA’s lifetime-drinking water health advisory level of 70 ppt for combined concentrations of PFOA and PFOS. The average age of groundwater entering the screens of the supply wells sampled in 2018 ranged from 8.2 to 45.8 years, with a median groundwater age of 25.7 years. Groundwater age was positively correlated with well depth and negatively correlated with dissolved oxygen. Groundwater age and PFAS concentrations were negatively correlated in the Columbia aquifer. Data from the 23 resampled wells indicate a significant positive difference in the average modeled groundwater-sample-age results. The average groundwater age from samples collected in 2018 was generally 5 years older than the average groundwater age from samples collected in 2008. The same pattern was found during cycle two (2008) of this study, where the 2008 groundwater age was on average 7 years older than the samples collected in 2000. The distribution of groundwater sample ages among the 17 trend wells and during the three study cycles (2000, 2008, and 2018) indicates that sample-age medians were statistically different from zero; well-water sample-age data show a slight increase in groundwater sample age.

Delaware↗
USGS2003

Methodology for estimating times of remediation associated with monitored natural attenuation

Natural attenuation processes combine to disperse, immobilize, and biologically transform anthropogenic contaminants, such as petroleum hydrocarbons and chlorinated ethenes, in ground-water systems. The time required for these processes to lower contaminant concentrations to levels protective of human health and the environment, however, varies widely between different hydrologic systems, different chemical contaminants, and varying amounts of contaminants. This report outlines a method for estimating timeframes required for natural attenuation processes, such as dispersion, sorption, and biodegradation, to lower contaminant concentrations and mass to predetermined regulatory goals in groundwater systems. The time-of-remediation (TOR) problem described in this report is formulated as three interactive components: (1) estimating the length of a contaminant plume once it has achieved a steady-state configuration from a source area of constant contaminant concentration, (2) estimating the time required for a plume to shrink to a smaller, regulatoryacceptable configuration when source-area contaminant concentrations are lowered by engineered methods, and (3) estimating the time needed for nonaqueous phase liquid (NAPL) contaminants to dissolve, disperse, and biodegrade below predetermined levels in contaminant source areas. This conceptualization was used to develop Natural Attenuation Software (NAS), an interactive computer aquifers. NAS was designed as a screening tool and requires the input of detailed site information about hydrogeology, redox conditions, and the distribution of contaminants. Because NAS is based on numerous simplifications of hydrologic, microbial, and geochemical processes, the program may introduce unacceptable errors for highly heterogeneous hydrologic systems. In such cases, application of the TOR framework outlined in this report may require more detailed, site-specific digital modeling. The NAS software may be downloaded from the Web site http://www.cee.vt.edu/NAS/ Application of NAS illustrates several general characteristics shared by all TOR problems. First, the distance of stabilization of a contaminant plume is strongly dependent on the natural attenuation capacity of particular ground-water systems. The time that it takes a plume to reach a steady-state configuration, however, is independent of natural attenuation capacity. Rather, the time of stabilization is most strongly affected by the sorptive capacity of the aquifer, which is dependent on the organic matter content of the aquifer sediments, as well as the sorptive properties of individual contaminants. As a general rule, a high sorptive capacity retards a plume.s growth or shrinkage, and increases the time of stabilization. Finally, the time of NAPL dissolution depends largely on NAPL mass, composition, geometry, and hydrologic factors, such as ground-water flow rates. An example TOR analysis for petroleum hydrocarbon NAPL was performed for the Laurel Bay site in South Carolina. About 500 to 1,000 pounds of gasoline leaked into the aquifer at this site in 1991, and the NAS simulations suggested that TOR would be on the order of 10 years for soluble and poorly sorbed compounds, such as benzene and methyl tertiary-butyl ether (MTBE). Conversely, TOR would be on the order of 40 years for less soluble, more strongly sorbed compounds, such as toluene, ethylbenzene, and xylenes (TEX). These TOR estimates are roughly consistent with contaminant concentrations observed over 10 years of monitoring at this site where benzene and MTBE concentrations were observed to decrease rapidly and are approaching regulatory maximum concentration limits, whereas toluene, ethylbenzene, and xylene concentrations decreased at a slower rate and have remained relatively high. An example TOR analysis for petroleum hydrocarbon NAPL was performed for the Laurel Bay site in South Carolina.

South Carolina↗
USGS2023

Hydrogeology, karst, and groundwater availability of Monroe County, West Virginia

Monroe County is in southeastern West Virginia, encompassing an area of 474 square miles. The area consists of karst and siliciclastic aquifers of Ordovician, Silurian, Devonian, and Mississippian age and is in parts of two physiographic provinces: the Valley and Ridge Province to the east of Peters Mountain, and the Appalachian Plateau Province to the west of Peters Mountain. This study was developed in response to inquiries from the Monroe County Commission requesting assessment of the water resources of the county to better understand the quantity of the county’s groundwater resources, for both current [2023] and future demand, and to provide information to support protection and management of the county’s valuable groundwater resources. Various products were developed for this study that provide knowledge with respect to water availability and contamination susceptibility of the karst aquifers within the county. U.S. Geological Survey (USGS) geologists conducted extensive geologic mapping in support of the project, producing (1) a countywide bedrock geologic map, (2) a countywide hydrogeologic map, and (3) a light detection and ranging (lidar)-derived countywide digital elevation model and associated sinkhole map. A significant part of this work was to map in detail the Greenbrier Group at the formation level, which prior to this study had only partially been completed. The report also includes (4) a description of the lithologic units identified as part of the geologic mapping process. U.S. Geological Survey hydrologists completed several additional products for the hydrology part of the effort, including development of (1) a countywide potentiometric surface (water-table) map, (2) a countywide base-flow stream assessment, (3) countywide water-budget estimates, (4) well log surveys for 15 wells to better understand subsurface controls on groundwater flow within the study area, (5) two groundwater tracer tests to better refine the groundwater divide from the northern and southern parts of the karst aquifer in Monroe County; and finally, based on all available data collected for the study including the potentiometric surface map, geologic map, current [2023] and legacy fluorometric groundwater tracer tests, and base-flow stream assessments, (6) groundwater-basin delineations were reassessed for principal groundwater basins within the Greenbrier aquifer. In Monroe County, four principal hydrogeologic settings produce large yields of water for residential, agricultural, and other uses. The most relied upon water-bearing zone with respect to current [2023] public water supply is from springs along Peters Mountain. These springs are derived from intervals of fractured sandstone and resultant alluvial deposits. Groundwater flows downslope through these permeable alluvial deposits and discharges at the contact with less permeable strata, such as the Reedsville Shale. The second most relied upon water-bearing zone in Monroe County is within the karstic Greenbrier Group aquifer, in which the basal Hillsdale Limestone overlies the less permeable Maccrady Shale. This geologic contact between the Hillsdale Limestone and Maccrady Shale is not only targeted as a source of water for agricultural supply but also is targeted as a source of water for residential supply. The third most relied upon water-bearing zone is composed of shallow perched aquifers within the Greenbrier Group. The discontinuous nature of these perched aquifers makes mapping their extent impossible, but they are related to permeable geologic strata, such as karstified limestones with solutionally enhanced permeability that overlies less permeable shale or chert bedrock. During geologic mapping of the county, several of these perched aquifers were documented in the Pickaway, Union, and Alderson Limestones. A fourth zone consists of springs from Ordovician carbonates at the base of Peters Mountain, which are influenced by sinking streams as well as upwelling along faults. In terms of water quantity, the most sustainable springs are those having deeper-sourced flows. Public supplies are a principal source of water used for residential and commercial supply in the region, accounting for 0.49 million gallons per day (Mgal/d) of fresh-water withdrawals (0.14 Mgal/d of groundwater and 0.35 Mgal/d of surface water) for residential and commercial use and serving 6,645 individuals (49.2 percent of the population). An estimated 6,861 people, (50.8 percent of the population) primarily rely on private wells or other unregulated sources, such as springs, and withdraw 0.55 Mgal/d of groundwater for their residential use. Public water supply in the region is primarily (71.4 percent) derived from springs and augmented by stream withdrawals (backup sources mainly during low-flow periods), with the remaining portion (28.6 percent) derived from groundwater withdrawals from wells. For rural residents, however, 100 percent of their withdrawals are derived from groundwater (wells or springs).

West Virginia↗
USGS2014

Executive summary

Fossil fuels from the Appalachian basin region have been major contributors to the Nation’s energy needs over much of the last three centuries. Early records indicate that Appalachian coal was first mined in the middle 1700s (Virginia and Pennsylvania) and was used sparingly to fuel colonial settlements and, later, a fledgling industrial-based economy along the eastern seaboard of the United States (de Witt and Milici, 1989). In 2011, central Appalachian basin coal production accounted for approximately 77 percent of all U.S. metallurgical (or coking) coal and 29 percent of total U.S. production (U.S. Energy Information Administration, 2013). Following initial discoveries and commercial use in western New York (1821) and Ohio and West Virginia (mid-1830s), the Appalachian petroleum (oil and gas) industry began in earnest in 1859 with the discovery of oil at the Drake well in northwestern Pennsylvania. Between 1860 and 1989, the Appalachian basin produced more than 2.5 billion barrels of oil (BBO) and more than 30 trillion cubic feet of gas (TCFG) from more than 500,000 wells (de Witt and Milici, 1989). Although both oil and gas continue to be produced in the Appalachian basin, most new wells in the region are drilled in shale reservoirs to produce natural gas. Appalachian coal and petroleum resources are still available in sufficient quantities to contribute significantly to the Nation’s energy needs. For example, the U.S. Energy Information Administration (2010) estimated that there are 6,484 million short tons of recoverable coal reserves in the Appalachian basin. Similarly, about 14.7 billion barrels of oil equivalent (BBOE) (1.2 BBO+81 TCFG [or 13.5 BBOE]) of recoverable Appalachian basin oil and gas remain available of an estimated ultimate endowment of approximately 25.5 billion BBOE (cumulative production + reserves + estimated recoverable undiscovered resources) (this volume, chap. C.1). U.S. Geological Survey (USGS) Professional Paper 1708 is a modern, indepth collection of reports, cross sections, and maps that describe the geology of the Appalachian basin and its fossil fuel resources. Several of the chapters have been published in outside journals or as other USGS publications. Although this volume is not a comprehensive regional treatment of all notable geologic and fossil fuel localities in the Appalachian basin, the selected study areas and topics presented in the chapters cover large segments of the basin and a wide range of stratigraphic intervals. As the title implies, this volume addresses topics that refer to the locations of coal and petroleum accumulations, the stratigraphic and structural framework, and the geochemical characteristics of the coal beds and petroleum in the basin, as well as the results and documentation of recent USGS assessments of coal, oil, and gas resources in the basin. Many of the maps and accompanying data supporting the reports in this volume are available as downloadable geographic information system (GIS) data files (such as selected coal beds, selected oil and gas fields, locations of oil and gas wells, coal production, coal chemistry, total petroleum system (TPS) boundaries, and bedrock geology). Log ASCII Standard (LAS) files for geophysical (gamma ray) wireline well logs also are included. This publication supplements and updates older USGS regional studies of Appalachian basin coal and petroleum resources such as those by Arndt and others (1968) and the numerous contributors to USGS Miscellaneous Map Series I−917 (for example, Harris and others, 1978), respectively. USGS Professional Paper 1708 is intended primarily for geoscientists in academia, industry, and government who are interested in Appalachian basin geology and its coal and petroleum resources. Other users, however, may find the wide variety of topics, papers, and digital images of value for landuse and policy planning issues. Among the anticipated benefits of the report are improvements in (1) resource assessment estimates and methodology, (2) exploration strategies, (3) basin models, and (4) energy use policies.

Appalachian Basin↗
USGS2009

Natural offshore oil seepage and related tarball accumulation on the California coastline — Santa Barbara Channel and the southern Santa Maria Basin; source identification and inventory

Oil spillage from natural sources is very common in the waters of southern California. Active oil extraction and shipping is occurring concurrently within the region and it is of great interest to resource managers to be able to distinguish between natural seepage and anthropogenic oil spillage. The major goal of this study was to establish the geologic setting, sources, and ultimate dispersal of natural oil seeps in the offshore southern Santa Maria Basin and Santa Barbara Basins. Our surveys focused on likely areas of hydrocarbon seepage that are known to occur between Point Arguello and Ventura, California. Our approach was to 1) document the locations and geochemically fingerprint natural seep oils or tar; 2) geochemically fingerprint coastal tar residues and potential tar sources in this region, both onshore and offshore; 3) establish chemical correlations between offshore active seeps and coastal residues thus linking seep sources to oil residues; 4) measure the rate of natural seepage of individual seeps and attempt to assess regional natural oil and gas seepage rates; and 5) interpret the petroleum system history for the natural seeps. To document the location of sub-sea oil seeps, we first looked into previous studies within and near our survey area. We measured the concentration of methane gas in the water column in areas of reported seepage and found numerous gas plumes and measured high concentrations of methane in the water column. The result of this work showed that the seeps were widely distributed between Point Conception east to the vicinity of Coal Oil Point, and that they by in large occur within the 3-mile limit of California State waters. Subsequent cruises used sidescan and high resolution seismic to map the seafloor, from just south of Point Arguello, east to near Gaviota, California. The results of the methane survey guided the exploration of the area west of Point Conception east to Gaviota using a combination of seismic instruments. The seafloor was mapped by sidescan sonar, and numerous lines of high -resolution seismic surveys were conducted over areas of interest. Biomarker and stable carbon isotope ratios were used to infer the age, lithology, organic matter input, and depositional environment of the source rocks for 388 samples of produced crude oil, seep oil, and tarballs mainly from coastal California. These samples were used to construct a chemometric fingerprint (multivariate statistics) decision tree to classify 288 additional samples, including tarballs of unknown origin collected from Monterey and San Mateo County beaches after a storm in early 2007. A subset of 9 of 23 active offshore platform oils and one inactive platform oil representing a few oil reservoirs from the western Santa Barbara Channel were used in this analysis, and thus this model is not comprehensive and the findings are not conclusive. The platform oils included in this study are from west to east: Irene, Hildago, Harvest, Hermosa, Heritage, Harmony, Hondo, Holly, Platform A, and Hilda (now removed). The results identify three “tribes” of 13 C-rich oil samples inferred to originate from thermally mature equivalents of the clayey-siliceous, carbonaceous marl, and lower calcareous-siliceous members of the Monterey Formation. Tribe 1 contains four oil families having geochemical traits of clay-rich marine shale source rock deposited under suboxic conditions with substantial higher-plant input. Tribe 2 contains four oil families with intermediate traits, except for abundant 28,30-bisnorhopane, indicating suboxic to anoxic marine marl source rock with hemipelagic input. Tribe 3 contains five oil families with traits of distal marine carbonate source rock deposited under anoxic conditions with pelagic but little or no higher-plant input. Tribes 1 and 2 occur mainly south of Point Conception in paleogeographic settings where deep burial of the Monterey Formation source rock favored generation from all three members or their equivalents. In this area, oil from the clayey-siliceous and carbonaceous marl members (Tribes 1 and 2) may overwhelm that from the lower calcareous-siliceous member (Tribe 3) because the latter is thinner and less oil-prone than the overlying members. Tribe 3 occurs mainly north of Point Conception, where shallow burial caused preferential generation from the underlying lower calcareous-siliceous member or another unit with similar characteristics. It is very desirable to be able to clearly distinguish the naturally occurring seep oils from the anthropogenically derived platform oils. Within the “training set” of oils and tars (388 samples), the biomarker parameters are sometimes sufficient to allow unique discrimination of individual platform oils. More often however, platform samples and seep samples with sources geographically close to each other are too similar to each other, with respect to the biomarker parameters, to definitively differentiate them on that basis alone. In some cases other parameters can be helpful. These other parameters are related to the degree of biogeochemical degradation or weathering that the oils or tars have experienced. These components include the typical oil distribution of n-alkane hydrocarbons and isoprenoids pristane and phytane. All of the platform oils in our sample set contain these components. On the other hand, the seep oils or tars have been exposed to significant biodegradation while in the near subsurface. The majority, but not all of seep oils or tars have been biodegraded up to or beyond the loss of n-alkanes and isoprenoids. Seep oils found in the vicinity of Coal Oil Point or Arroyo Burro are apparently the least weathered and are particularly likely to retain significant n-alkanes and isoprenoids. Therefore the combination of chemometric fingerprinting and the presence or absence of n-alkanes and isoprenoids help to differentiate anthropogenic production oils versus natural seeps oils and tars. The differentiation is not always definitive because of the close chemical similarity of some samples and the variability in the biodegradation progression. This is the case near Coal Oil Point, and near Platform A (Dos Cuadros Field) where seep oils and Platform Holly and Platform A oils are genetically very similar and cannot be definitively distinguished after a period of a few days of weathering. In contrast, oils from the Point Conception platforms can be distinguished on the basis of chemometric fingerprinting alone. In the middle of this spectrum are oils from Platforms Harmony, Heritage, and Hondo, where it is expected that oil weathering would take on the order of two weeks to a month to produce tarballs similar to those seen near Point Conception. In this case there is a much greater degree of weathering needed to proceed from produced oil to the biodegraded tar characteristic of tarball stranded on the beach. Tar deposition on beaches was monitored as part of cooperative with the County of Santa Barbara Energy Division and the U.S. Geological Survey during 2001-2003. We found tar deposition varies on a seasonal basis. In general, tarballs accumulate at a faster rate or remain longer on all beaches during the summer and fall months. The reasons for this are unclear based on our limited observations, however we speculate that factors such as prevailing winds and currents combined with more quiescent wave conditions favors the accumulation and preservation of tarballs on the beach during the summer and fall months. In contrast, winter storms, with much greater wave action remove beach sand and other materials, and stormy seas tend to break up oil that might weather into tarballs. Natural seepage is affected by the spring/neap tidal cycle; however, the link to tar deposition is unclear. Longer periods of monitoring are needed to address the variability in the data and provide a more robust statistical analysis.

California↗
USGS2026

Guidelines for producing integrated 210Pb and 14C age-models

Accurate reconstructions of past environmental changes are crucial in paleoecological research and require reliable chronologies of sedimentary archives. Establishing robust age-models and obtaining the most appropriate proxies for analysis is a complex scientific endeavor, requiring extensive resources and collaboration among specialists, including radiochronologists. Radiometric dating methods, such as 210 Pb and radiocarbon ( 14 C), are frequently employed to establish chronologies in aquatic sedimentary deposits and peat bogs. In this study, we review key aspects of sampling, analysis, and the principles underlying 210 Pb and 14 C age-models, focusing on methods for developing robust joint chronologies for paleoenvironmental research. Drawing largely from the authors' experiences and group discussions during and after a scientific workshop in 2022, we discuss important considerations for site selection, sampling strategies, and radiometric dating to construct integrated 210 Pb and 14 C age-models. Using expert consensus, this group – called Paleostats – aims to provide a set of best practices for other geochronologists with this methods paper. Among our conclusions, we emphasize the importance of accounting for site-specific factors such as prior information on sedimentation rates to establish appropriate sampling and analytical strategies. The use of appropriate coring devices can minimize disturbance to sediments and ensure the core surface remains intact and preserved until sectioning. Where excess 210 Pb is expected, sectioning at intervals of ≤1 cm provides an adequate sampling resolution for 210 Pb dating. Exceptions are possible, allowing for ~2–3 cm sections in areas with confirmed high sedimentation rates (e.g., > 1 cm yr −1 ). Recovering deeper core sections for 14 C dating with sufficient overlap allows for accounting errors in depth estimates made in the field. Special attention is advised during time intervals where validation proxies, such as the human-made radionuclides 137 Cs or post-bomb 14 C, are expected, and to determine the depth of secular equilibrium between 210 Pb and 226 Ra. Radiocarbon analyses are commonly performed by accelerator mass spectrometry, and age models are constructed mainly using Bayesian statistics with Markov Chain Monte Carlo techniques (e.g., Bacon ). A Bayesian approach ( Plum ) is now available for producing 210 Pb age-models, which infers the 210 Pb ex flux, eliminates the need for selecting an equilibrium depth, and allows dating cores with incomplete 210 Pb ex inventory. Plum offers improved chronologies by integrating raw 210 Pb and 14 C data, and these age-models can be enriched with other dating methodologies, such as identifying tephras and other well-recorded historical events. Harmonized reporting would contribute to making radiometric age-models reproducible, which would benefit from an international effort. Using 210 Pb and 14 C to produce integrated age-models may yield better insights into the interplay between natural and recent anthropogenic forcings on ecosystems. This can enhance our understanding of environmental processes and their impacts on climate change, ultimately supporting science-based assessments and decisions.

Earth-Science Reviews↗
USGS2001

Documentation and verification of VST2D; a model for simulating transient, Variably Saturated, coupled water-heat-solute Transport in heterogeneous, anisotropic 2-Dimensional, ground-water systems with variable fluid density

This report describes a model for simulating transient, Variably Saturated, coupled water-heatsolute Transport in heterogeneous, anisotropic, 2-Dimensional, ground-water systems with variable fluid density (VST2D). VST2D was developed to help understand the effects of natural and anthropogenic factors on quantity and quality of variably saturated ground-water systems. The model solves simultaneously for one or more dependent variables (pressure, temperature, and concentration) at nodes in a horizontal or vertical mesh using a quasi-linearized general minimum residual method. This approach enhances computational speed beyond the speed of a sequential approach. Heterogeneous and anisotropic conditions are implemented locally using individual element property descriptions. This implementation allows local principal directions to differ among elements and from the global solution domain coordinates. Boundary conditions can include time-varying pressure head (or moisture content), heat, and/or concentration; fluxes distributed along domain boundaries and/or at internal node points; and/or convective moisture, heat, and solute fluxes along the domain boundaries; and/or unit hydraulic gradient along domain boundaries. Other model features include temperature and concentration dependent density (liquid and vapor) and viscosity, sorption and/or decay of a solute, and capability to determine moisture content beyond residual to zero. These features are described in the documentation together with development of the governing equations, application of the finite-element formulation (using the Galerkin approach), solution procedure, mass and energy balance considerations, input requirements, and output options. The VST2D model was verified, and results included solutions for problems of water transport under isohaline and isothermal conditions, heat transport under isobaric and isohaline conditions, solute transport under isobaric and isothermal conditions, and coupled water-heat-solute transport. The first three problems considered in model verification were compared to either analytical or numerical solutions, whereas the coupled problem was compared to measured laboratory results for which no known analytic solutions or numerical models are available. The test results indicate the model is accurate and applicable for a wide range of conditions, including when water (liquid and vapor), heat (sensible and latent), and solute are coupled in ground-water systems. The cumulative residual errors for the coupled problem tested was less than 10–8 cubic centimeter per cubic centimeter, 10-5 moles per kilogram, and 102 calories per cubic meter for liquid water content, solute concentration and heat content, respectively. This model should be useful to hydrologists, engineers, and researchers interested in studying coupled processes associated with variably saturated transport in ground-water systems.

Water-Resources Investigations Report↗
USGS1992

Geohydrology and simulated ground-water flow, Plymouth-Carver aquifer, southeastern Massachusetts

The Plymouth-Carver aquifer underlies an area of 140 square miles and is the second largest aquifer in areal extent in Massachusetts. It is composed primarily of saturated glacial sand and gravel. The water-table and bedrock surface were mapped and used to determine saturated thickness of the aquifer, which ranged from less than 20 feet to greater than 200 feet. Ground water is present mainly under unconfined conditions, except in a few local areas such as beneath Plymouth Harbor. Recharge to the aquifer is derived almost entirely from precipitation and averages about 1.15 million gallons per day per square mile. Water discharges from the aquifer by pumping, evapotranspiration, direct evaporation from the water table, and seepage to streams, ponds, wetlands, bogs, and the ocean. In 1985, water use was about 59.6 million gallons per day, of which 82 percent was used for cranberry production. The Plymouth-Carver aquifer was simulated by a three-dimensional, finite difference ground-water-flow model. Most model boundaries represent the natural hydrologic boundaries of the aquifer. The model simulates aquifer recharge, withdrawals by pumped wells, leakage through streambeds, and discharge to the ocean. The model was calibrated for steady-state and transient conditions. Model results were compared with measured values of hydraulic head and ground-water discharge. Results of simulations indicate that the modeled ground-water system closely simulates actual aquifer conditions. Four hypothetical ground-water development alternatives were simulated to demonstrate the use of the model and to examine the effects on the ground-water system. Simulation of a 2-year period of no recharge and average pumping rates that occurred from 1980-85 resulted in water-level declines exceeding 5 feet throughout most of the aquifer and a decrease of 54 percent in average ground-water discharge to streams. In a second simulation, four wells in the northern part of the area were pumped at 10.4 million gallons per day in excess of rates simulated in the steady-state model for the four wells. This resulted in water-level declines of 2 feet or more in an area of 25 square miles and a decline in average ground-water discharge to streams of 6 percent. When this pumpage was simulated as recharge to the aquifer, water levels beneath the recharge area rose more than 40 feet, and ground-water discharge remained equal to average discharge in the calibrated steady-state model. In a third simulation, all 21 existing production wells were pumped at nearly the design capacity of 17.8 million gallons per day; this pumping rate produced water-level declines of less than 2 feet throughout most of the aquifer. When simulated pumpage was increased to 32.8 million gallons per day from existing wells and from 15 additional wells, the area where water-level declines exceeded 2 feet significantly increased. In another set of simulations, a well field close to a stream was pumped at rates of 2, 4, and 6 million gallons per day. At a pumping rate of 6 million gallons per day, ground-water discharge to the stream decreased 34 percent during periods of normal precipitation and 56 percent during drought conditions.

Massachusetts↗
USGS2011

Geohydrology, simulation of regional groundwater flow, and assessment of water-management strategies, Twentynine Palms area, California

The Marine Corps Air Ground Combat Center (MCAGCC) Twentynine Palms, California, overlies the Surprise Spring, Deadman, Mesquite, and Mainside subbasins of the Morongo groundwater basin in the southern Mojave Desert. Historically, the MCAGCC has relied on groundwater pumped from the Surprise Spring subbasin to provide all of its potable water supply. Groundwater pumpage in the Surprise Spring subbasin has caused groundwater levels in the subbasin to decline by as much as 190 feet (ft) from 1953 through 2007. Groundwater from the other subbasins contains relatively high concentrations of fluoride, arsenic, and (or) dissolved solids, making it unsuitable for potable uses without treatment. The potable groundwater supply in Surprise Spring subbasin is diminishing because of pumping-induced overdraft and because of more restrictive Federal drinking-water standards on arsenic concentrations. The U.S. Geological Survey, in cooperation with the MCAGCC, completed this study to better understand groundwater resources in the area and to help establish a long-term strategy for regional water-resource development. The Surprise Spring, Deadman, Mesquite, and Mainside subbasins are filled with sedimentary deposits of Tertiary age, alluvial fan deposits of Quaternary-Tertiary age, and younger alluvial and playa deposits of Quaternary age. Combined, this sedimentary sequence reaches a maximum thickness of more than 16,000 ft in the Deadman and Mesquite subbasins. The sedimentary deposits of Tertiary age yield a small amount of water to wells, and this water commonly contains high concentrations of fluoride, arsenic, and dissolved solids. The alluvial fan deposits form the principal water-bearing unit in the study area and have a combined thickness of 250 to more than 1,000 ft. The younger alluvial and playa deposits are unsaturated throughout most of the study area. Lithologic and downhole geophysical logs were used to divide the Quaternary/ Tertiary alluvial fan deposits into two aquifers (referred to as the upper and the middle aquifers) and the Tertiary sedimentary deposits into a single aquifer (referred to as the lower aquifer). In general, wells perforated in the upper aquifer yield more water than wells perforated in the middle and lower aquifers. The study area is dominated by extensive faulting and moderate to intense folding that has displaced or deformed the pre-Tertiary basement complex as well as the overlying Tertiary and Quaternary deposits. Many of these faults act as barriers to the lateral movement of groundwater flow and form many of the boundaries of the groundwater subbasins. The principal recharge to the study area is groundwater underflow across the western and southern boundaries that originates as runoff in the surrounding mountains. Groundwater discharges naturally from the study area as spring flow, as groundwater underflow to downstream basins, and as water vapor to the atmosphere by transpiration of phreatophytes and direct evaporation from moist soil. The annual volume of water that naturally recharged to or discharged from the groundwater flow system in the study area during predevelopment conditions was estimated to be 1,010 acre-feet per year (acre-ft/yr). About 90 percent of this recharge originated as runoff from the Little San Bernardino and the Pinto Mountains to the south, and the remainder originated as runoff from the San Bernardino Mountains to the west. Evapotranspiration by phreatophytes near Mesquite Lake (dry) was the primary form of predevelopment groundwater discharge. From 1953 through 2007, approximately 139,400 acre-feet (acre-ft) of groundwater was pumped by the MCAGCC from the Surprise Spring subbasin. A regional-scale numerical groundwater flow model was developed using MODFLOW-2000 for the Surprise Spring, Deadman, Mesquite, and Mainside subbasins. The aquifer system was simulated by using three model layers representing the upper, middle, and lower aquifers. Measured groundwater levels

California↗
USGS2014

Simulated effects of existing and proposed surface-water impoundments and gas-well pads on streamflow and suspended sediment in the Cypress Creek watershed, Arkansas

Cypress Creek is located in central Arkansas and is the main tributary to Brewer Lake, which serves as the primary water supply for Conway, Arkansas, and the surrounding areas. A model of the Cypress Creek watershed was developed and calibrated in cooperation with Southwestern Energy Company using detailed precipitation, streamflow, and discrete suspended-sediment data collected from 2009 through 2012. These data were used with a Hydrologic Simulation Program—FORTRAN model to address different potential gas-extraction activities within the watershed. The calibrated Hydrologic Simulation Program—FORTRAN model was used to simulate four land-use scenarios and examine the potential effects of these land-use changes on the streamflow and water quality within the Cypress Creek watershed. These simulated scenarios included (1) the conversion of all nonforested land to forest, representing a time period before extensive grazing activities and no gas-extraction activities; (2) a land-use change to that of 1949, representing a time period with some grazing activities and no gas-extraction activities; (3) a time period with current land-use conditions, but without any gas-extraction activities, that is, the exclusion of gas-well pads/pipelines, associated gravel roads, and surface-water impoundments; and (4) a time period with current land-use conditions, but with increased gas-extraction activities (for example, increased gas-well pad and surface-water impoundment activities) to represent a possible future natural gas full-development condition for the area. A current-conditions simulation also was built and calibrated and represents the current conditions (2013) within the watershed. This simulation was used as the comparison basis for the four land-use scenarios described above. The current-conditions simulation used the 2006 land-use conditions, which consisted primarily of forest and pasture, as well as the current (2013) 35 gas-well pads and pipelines and 6 surface-water impoundments, which account for approximately 1.6 percent of the land use. Simulating a time period before extensive-grazing activities and no gas-extraction activities for scenario 1 resulted in a decrease in suspended-sediment loads and volume of streamflow within the Cypress Creek watershed compared to the current-conditions simulation. Simulating a time period before any gas-extraction activities but with some grazing activities for scenario 2 also resulted in a decrease in suspended-sediment loads and volume of streamflow within the Cypress Creek watershed. Simulating current conditions, but without any natural gas-pad land use or related activities (including pipelines and associated gravel roads), for scenario 3 resulted in mostly unchanged suspended-sediment loads and volume of streamflow within the Cypress Creek watershed, as compared to the current-conditions simulation. Finally, simulating potential future conditions of increased gas-well pad and surface-water impoundment activities for scenario 4 resulted in a decrease (compared to the current-conditions simulation) in suspended-sediment loads and a slight increase of volume of streamflow within the Cypress Creek watershed. The Arkansas Natural Resources Commission and the Arkansas Department of Environmental Quality list suspended sediment from “poor pastures” as a primary source of nonpoint-source pollution in north-central Arkansas, but unpaved (gravel) roads are another important source of suspended sediment. Because of the high sediment-loading rates associated with gravel roads and the large amount of pasture within the watershed, the factors most responsible for suspended sediment within the Cypress Creek watershed are likely associated more with the pastureland and gravel roads, than factors associated with gas-well pads/pipelines.

Arkansas↗
USGS2017

Assessing the risk of non-native marine species in the Bering Sea

Invasive species are one of the leading global conservation concerns, which can have strong, negative impacts on ecosystems, vulnerable species, and valuable natural resources. Arctic regions have experienced a relatively low number of biological introductions to date. Their geographical remoteness, cold waters, and presence of sea ice present challenging conditions for both non-native organisms and the vessels that transport them, presumably leading to low rates of introduction and establishment. However, observed increases in water temperatures reductions in sea ice, and projected increases in shipping traffic are expected to render arctic marine regions more susceptible to the arrival and colonization of marine invasives. Risk assessments for these Arctic regions are important to inform management and monitoring priorities by determining which species pose the greatest risk. To this end, we developed a ranking system for non-native marine species and used this system to assess the risk of non-native species to the Bering Sea. Using species’ published physiological tolerances, we mapped habitat suitability under current and future climate scenarios to identify geographic areas of current and future concern. In addition, we described shipping traffic from commercial and fishing vessels to identify ports of entry for non-native species. Collectively, these analyses identify which marine species have the greatest risk for invasion, where in the Bering Sea invasion risk and species establishment is greatest, and which ports are most likely to serve as an entry point for marine invasives into Alaska’s Bering Sea. The ranking system we developed for non-native marine species consists of 33 questions grouped into five categories. The first four categories evaluate a species’ ability to arrive and establish in the Bering Sea, its reliance on humans for introductions, its biology, and its impacts on ecological and human systems. The fifth category is not included in the total ranking score, but provides information on management considerations. The ranking system has methods to account for data deficiencies and calculates these deficiencies to allow readers to weigh the lack of knowledge with the ranking score. We prioritized non-native species for ranking based on their geographic proximity to the Bering Sea. We evaluated 46 species and ranking scores ranged from 29.1 to 74.3 (out of a possible 100), with highest scores indicating greatest risk. Taxonomy at the level of phylum did not explain variation in ranking values, likely due to the substantial biological variation relative to our ranking criteria among members of the same phylum. To investigate where non-native species may survive and persist in the Bering Sea, we compared species’ temperature and salinity thresholds to environmental conditions of the Bering Sea. Environmental conditions were obtained from three Regional Ocean Modeling Systems (ROMS) and investigated under two time periods: current (2003-2012) and mid-century (2030-2039). We identified potential habitat for survival for 42 species, and potential habitat for reproduction for 29 species. Under current conditions, all species had temperature and salinity thresholds that would allow survival in the Bering Sea for at least part of the year, and most species (79% to 83%) had thresholds that would allow for survival year-round. For species with temperature and salinity thresholds unsuitable for survival in the Bering Sea, winter temperatures appear to be the limiting factor. Most species had six to nine weeks of suitable conditions for reproduction. Future increases in water temperatures are expected to open more habitat for marine invasives. Two of the three ROMs project an increase in the number of non-native species that would be able to survive year-round by mid-century. Moreover, models project between 37% and 60% of the Bering Sea shelf habitat to become more suitable under mid-century climate condition.

Alaska↗
USGS2007

Applying radar technology to migratory bird conservation and management: Strengthening and expanding a collaborative

Executive Summary Understanding the factors affecting migratory bird and bat populations during all three phases of their life cycle—breeding, nonbreeding, and migration—is critical to species conservation planning. This includes the need for information about these species’ responses to natural challenges, as well as information about the effects of human activities and structures. Habitats and other resources critical to migrants during passage and stopover are being destroyed, degraded, and threatened by human activities. Birds and bats are also uniquely susceptible to human use of the airspace. Wind turbines, communication and power transmission towers, and other tall structures, known to cause bird and bat mortality, are being erected or proposed in increasing numbers across the country. In addition, the potential for bird/aircraft collisions poses human safety threats. Management and regulatory agencies, conservation organizations, and industry currently lack the information they need to meet their missions and statutory responsibilities. The biological data available from various radar technologies offer a unique opportunity to learn more about the spatiotemporal distribution patterns, flight characteristics, and habitat use of “aero-fauna.” Recognizing the opportunities presented by radar technologies, the U.S. Geological Survey (USGS), the U.S. Fish and Wildlife Service (USFWS), and university partners collaborated first on individual projects and then in a broader, informal “collaborative” to coordinate their radar-related research and work together to develop the suite of products needed for conservation of birds and bats. Having produced two summary documents (Sojda and others, 2005; Ruth and others, 2005), the next objective was to convene a workshop for researchers, management and regulatory agencies, and other interested parties. The focus of this initial workshop was on strengthening the existing USGS-USFWS-university partnership and expanding the “collaborative” to include new Federal agency partners. The subject matter was centered on discussing available technologies, appropriate applications, management-related needs, and ways to strengthen collaborative research and conservation efforts. The workshop opened with presentations about the history of the “radar collaborative,” a description of the workshop objectives and focuses, and a summary of resource management and regulatory needs. Scientific presentations describing current research projects or subjects followed, given by USGS scientists, as well as scientists from other Federal agencies, academia, conservation and ornithological organizations, and a private contracting firm. Presenters addressed a wide variety of management issues including siting of wind-power facilities, bird/aircraft collisions, effects of hurricanes Katrina and Rita on bird migration, bird use of Conservation Reserve Program land, defining bird migration patterns at a broad regional scale, and associating migrant birds with their stopover habitats. Presentations described a variety of radar technologies including NEXRAD weather surveillance radar, modified mobile marine radar, military tracking radar, pencil beam radar, and dual polarization radar, as well as complementary techniques and analysis methods such as acoustic monitoring, thermal imaging, artificial intelligence, and individual-based modeling. Key issues, themes, and questions identified during the open discussions that followed fell into five main categories: (1) agency needs and challenges; (2) radar technology and applications—technical questions and issues; (3) tools and resources for managers and researchers; (4) standardization of protocols; and (5) collaborative opportunities. Participants identified the following management, regulatory, or business issues facing them which may be addressed with radar technologies: tall structures; wind turbines; identification and protection of key habitats; assessment of management activities; and bird/aircraft strikes. Participants frequently expressed the need for specific information about which radar technologies are best used for answering particular questions. User groups emphasized the importance of clear, defensible scientific information on which they can base their activities. In turn, researchers emphasized their need for clearly defined, specific questions from managers so that they can design and conduct the required research. Discussions about technical issues requiring further research and collaboration included target identity, ground-truthing, linking migrants to habitat, and standardized protocols for applied research. Workshop participants identified and endorsed a series of seven action items that would promote collaboration and begin to address key issues identified at the workshop: Action Item #1 : Establish a working subgroup to address large-scale surveillance radar standardization issues. Action Item #2 : Establish a working subgroup to address small-scale radar standardization issues. Action Item #3 : Bring management and regulatory agencies together to identify the three most important information needs for each key management issue relating to radar technologies. Action Item #4 : Develop Fact Sheet(s) to provide information about radar technology applications to migratory bird and bat conservation issues. Action Item #5 : Create a “radar collaborative” Website to provide information about radar biology applications, contacts, publications, and so forth. Action Item #6 : Formalize and expand the USGS-USFWS “radar collaborative.” Action Item #7 : Advance basic research, such as target identity and validation, which will support and improve our abilities to apply radar technologies to conservation objectives. There was considerable interest in expanding the “radar collaborative” to include those agencies, organizations, and industries represented at the workshop. It was felt that the publication of the workshop proceedings, implementation of action items, and additional future meetings or workshops will be crucial in strengthening the “radar collaborative” effort and promoting the use of these valuable technologies for conserving migratory species.

Open-File Report↗
USGS2003

Predicting changes in hydrologic retention in an evolving semi-arid alluvial stream

Hydrologic retention of solutes in hyporheic zones or other slowly moving waters of natural channels is thought to be a significant control on biogeochemical cycling and ecology of streams. To learn more about factors affecting hydrologic retention, we repeated stream-tracer injections for 5 years in a semi-arid alluvial stream (Pinal Creek, Ariz.) during a period when streamflow was decreasing, channel width increasing, and coverage of aquatic macrophytes expanding. Average stream velocity at Pinal Creek decreased from 0.8 to 0.2 m/s, average stream depth decreased from 0.09 to 0.04 m, and average channel width expanded from 3 to 13 m. Modeling of tracer experiments indicated that the hydrologic retention factor ( R h ), a measure of the average time that solute spends in storage per unit length of downstream transport, increased from 0.02 to 8 s/m. At the same time the ratio of cross-sectional area of storage zones to main channel cross-sectional area ( A s / A ) increased from 0.2 to 0.8 m 2 /m 2 , and average water residence time in storage zones ( t s ) increased from 5 to 24 min. Compared with published data from four other streams in the US, Pinal Creek experienced the greatest change in hydrologic retention for a given change in streamflow. The other streams differed from Pinal Creek in that they experienced a change in streamflow between tracer experiments without substantial geomorphic or vegetative adjustments. As a result, a regression of hydrologic retention on streamflow developed for the other streams underpredicted the measured increases in hydrologic retention at Pinal Creek. The increase in hydrologic retention at Pinal Creek was more accurately predicted when measurements of the Darcy–Weisbach friction factor were used (either alone or in addition to streamflow) as a predictor variable. We conclude that relatively simple measurements of channel friction are useful for predicting the response of hydrologic retention in streams to major adjustments in channel morphology as well as changes in streamflow.

Advances in Water Resources↗
USGS1989

Analysis of water surface and flow distribution for the design flood at a proposed highway crossing of the Sabine River near Tatum, Texas

The hydraulic effects of the proposed Texas Highway 43 crossing of the Sabine River near Tatum, Texas, were determined on the basis of results from a two-dimensional finite-element surface-water-flow model. In planning the replacement crossing by the Texas State Department of Highways and Public Transportation, approximations of apportionment of flow among the openings and velocities within the openings were of concern. The model was used to simulate flow in the river floodplain system for the proposed design, an alternate design, and for the natural condition. The proposed bridge design by the State Department of Highways and Public Transportation consisted of a 320-foot main channel opening, four left overflow bridges with widths of 120, 320, 320, and 280 feet, and one right overflow opening with a width of 440 feet. The alternate design consisted of a 950-foot main channel opening, three left overflow bridges with widths of 320, 320, and 280 feet, and one right overflow opening with a width of 200 feet. Preliminary one-dimensional computations were used as an aid in establishing the boundary conditions for the two-dimensional analysis. The analysis evaluates the interaction of complex hydraulic characteristics of varied features found at many multiple-opening highway crossings. Valid applications of accepted standard one-dimensional techniques would not reveal some aspects of the overall hydraulics of the crossing. The simulation of the proposed design shows, among other things, which openings are the most productive per foot of bridge length. The results of the two-dimensional simulation of the proposed design indicate some differences in the apportionment of flow among the openings when compared to the one-dimensional proposed crossing computations. The alternate design results in computed water-surface altitudes which were slightly lower than those of the proposed design. The alternate design would require less modification to the existing embankment. Velocities computed within the openings and at the bridge abutments, using the two-dimensional model, were within the design specifications of the State Department of Highways and Public Transportation. The simulations of the proposed and alternate designs indicate a lateral component of the water-surface slope at the embankment. Redistribution of flow across the floodplain also is indicated in both simulations. Some of the differences in the response between the two designs are affected by geometric features of the floodplain other than the embankment-opening geometry.

Water-Resources Investigations Report↗
USGS2014

Hydrostratigraphic interpretation of test-hole and borehole geophysical data, Kimball, Cheyenne, and Deuel Counties, Nebraska, 2011-12

Recently (2004) adopted legislation in Nebraska requires a sustainable balance between long-term supplies and uses of surface-water and groundwater and requires Natural Resources Districts to understand the effect of groundwater use on surface-water systems when developing a groundwater-management plan. The South Platte Natural Resources District (SPNRD) is located in the southern Nebraska Panhandle and overlies the nationally important High Plains aquifer. Declines in water levels have been documented, and more stringent regulations have been enacted to ensure the supply of ground-water will be sufficient to meet the needs of future generations. Because an improved understanding of the hydrogeologic characteristics of this aquifer system is needed to ensure sustainability of groundwater withdrawals, the U.S. Geological Survey, in cooperation with the SPNRD, Conservation and Survey Division of the University of Nebraska-Lincoln, and the Nebraska Environmental Trust, began a hydrogeologic study of the SPNRD to describe the lithology and thickness of the High Plains aquifer. This report documents these characteristics at 29 new test holes, 28 of which were drilled to the base of the High Plains aquifer. Herein the High Plains aquifer is considered to include all hydrologically connected units of Tertiary and Quaternary age. The depth to the base of aquifer was interpreted to range from 37 to 610 feet in 28 of the 29 test holes. At some locations, particularly northern Kimball County, the base-of-aquifer surface was difficult to interpret from drill cutting samples and borehole geophysical logs. The depth to the base of aquifer determined for test holes drilled for this report was compared with the base-of-aquifer surface interpreted by previous researchers. In general, there were greater differences between the base-of-aquifer elevation reported herein and those in previous studies for areas north of Lodgepole Creek compared to areas south of Lodgepole Creek. The largest difference was at test hole 5-SP-11, where an Ogallala-filled paleovalley prevously had been interpreted based on relatively sparse test-hole data west of 5-SP-11. The base of aquifer near test hole 5-SP-11 reported herein is approximately 230 ft higher in elevation than previously interpreted. Among other test holes that are likely to have been drilled in Ogallala-filled paleovalleys, the greatest difference in the interpreted base of aquifer was for test hole 7-CC-11, northeast of Potter, Nebraska, where the base of aquifer is 180 feet deeper than previously interpreted. Interpretation of test-hole and borehole geophysical data for 29 additional test holes will improve resource managers’ understanding of the hydrogeologic characteristics, including aquifer thickness. Aquifer thickness, which is related to total water in storage, is not well quantified in the north and south tablelands. The additional hydrostratigraphic interpretations provided in this report will improve the hydrogeologic framework used in current (2014) and future groundwater models, which are the basis for many water-management decisions.

Nebraska↗
USGS2004

Hydrogeology and simulation of regional ground-water-level declines in Monroe County, Michigan

Observed ground-water-level declines from 1991 to 2003 in northern Monroe County, Michigan, are consistent with increased ground-water demands in the region. In 1991, the estimated ground-water use in the county was 20 million gallons per day, and 80 percent of this total was from quarry dewatering. In 2001, the estimated ground-water use in the county was 30 million gallons per day, and 75 percent of this total was from quarry dewatering. Prior to approximately 1990, the ground-water demands were met by capturing natural discharge from the area and by inducing leakage through glacial deposits that cover the bedrock aquifer. Increased ground-water demand after 1990 led to declines in ground-water level as the system moves toward a new steady-state. Much of the available natural discharge from the bedrock aquifer had been captured by the 1991 conditions, and the response to additional withdrawals resulted in the observed widespread decline in water levels. The causes of the observed declines were explored through the use of a regional ground-water-flow model. The model area includes portions of Lenawee, Monroe, Washtenaw, and Wayne Counties in Michigan, and portions of Fulton, Henry, and Lucas Counties in Ohio. Factors, including lowered water-table elevations because of below average precipitation during the time period (1991 - 2001) and reduction in water supply to the bedrock aquifer because of land-use changes, were found to affect the regional system, but these factors did not explain the regional decline. Potential ground-water capture for the bedrock aquifer in Monroe County is limited by the low hydraulic conductivity of the overlying glacial deposits and shales and the presence of dense saline water within the bedrock as it dips into the Michigan Basin to the west and north of the county. Hydrogeologic features of the bedrock and the overlying glacial deposits were included in the model design. An important step of characterizing the bedrock aquifer was the determination of inputs and outputs of water—leakage from glacial deposits and flows across model boundaries. The imposed demands on the groundwater system create additional discharge from the bedrock aquifer, and this discharge is documented by records and estimates of water use including: residential and industrial use, irrigation, and quarry dewatering. Hydrologic characterization of Monroe County and surrounding areas was used to determine the model boundaries and inputs within the ground-water model. MODFLOW-2000 was the computer model used to simulate ground-water flow. Predevelopment, 1991, and 2001 conditions were simulated with the model. The predevelopment model did not include modern water use and was compared to information from early settlement of the county. The 1991 steady-state model included modern demands on the ground-water system and was based on a significant amount of data collected for this and previous studies. The predevelopment and 1991 simulations were used to calibrate the numerical model. The simulation of 2001 conditions was based on recent data and explored the potential ground-water levels if the current conditions persist. Model results indicate that the ground-water level will stabilize in the county near current levels if the demands imposed during 2001 are held constant.

Michigan↗
USGS2009

Sand resources, regional geology, and coastal processes of the Chandeleur Islands Coastal System: An evaluation of the Breton National Wildlife Refuge

Breton National Wildlife Refuge, the Chandeleur Islands chain in Louisiana, provides habitat and nesting areas for wildlife and is an initial barrier protecting New Orleans from storms. The U.S. Geological Survey (USGS) in partnership with the University of New Orleans Pontchartrain Institute for Environmental Sciences undertook an intensive study that included (1) an analysis of island change based on historical maps and remotely sensed shoreline and topographic data; (2) a series of lidar surveys at 3- to 4-month intervals after Hurricane Katrina to determine barrier island recovery potential; (3) a discussion of sea level rise and effects on the islands; (4) an analysis of sea floor evolution and sediment dynamics in the refuge over the past 150 years; (5) an assessment of the local sediment transport and sediment resource availability based on the bathymetric and subbottom data; (6) a carefully selected core collection effort to groundtruth the geophysical data and more fully characterize the sediments composing the islands and surrounds; (7) an additional survey of the St. Bernard Shoals to assess their potential as a sand resource; and (8) a modeling study to numerically simulate the potential response of the islands to the low-intensity, intermediate, and extreme events likely to affect the refuge over the next 50 years. Results indicate that the islands have become fragmented and greatly diminished in subaerial extent over time: the southern islands retreating landward as they reorganize into subaerial features, the northern islands remaining in place. Breton Island, because maintenance of the Mississippi River-Gulf Outlet (MRGO) outer bar channel requires dredging, is deprived of sand sufficient to sustain itself. Regional sediment transport trends indicate that large storms are extremely effective in transporting sand and controlling the shoreline development and barrier island geometry. Sand is transported north and south from a divergent zone near Monkey Bayou at the southern end of the Chandeleur Islands. Numerical simulation of waves and sediment transport supports the geophysical results and indicates that vast areas of the lower shoreface are affected and are undergoing erosion during storm events, that there is little or no fair weather mechanism to rework material into the littoral system, and that as a result, there is a net loss of sediment from the system. Lidar surveys revealed that the island chain immediately after Hurricane Katrina lost about 84 percent of its area and about 92 percent of its prestorm volume. Marsh platforms that supported the islands’ sand prior to the storm were reduced in width by more than one-half. Repeated lidar surveys document that in places the shoreline has retreated about 100 m under the relatively low-energy waves since Hurricanes Katrina and Rita; however, this retreat is nonuniform. Recent high-resolution geophysical surveys of the sea floor and subsurface within 5–6 km of the Chandeleur Islands during 2006 and 2007 show that, in addition to the sand that is rebuilding portions of the island chain, a large volume of sand is contained in Hewes Point, in an extensive subtidal spit platform that has formed at the northern end of the Chandeleur Islands. Hewes Point appears to be the depositional terminus of the alongshore transport system. In the southern Chandeleurs, sand is being deposited in a broad tabular deposit near Breton Island called the southern offshore sand sheet. These two depocenters account for approximately 70 percent of the estimated sediment volume located in potential borrow sites. An additional large potential source of sand for restoration lies in the St. Bernard Shoals, which are estimated to contain approximately 200 × 10 6 m 3 of sand. Successful restoration planning for the Breton National Wildlife Refuge should mimic the natural processes of early stages of barrier island evolution including lateral transport to the flanks of the island chain from a centralized sand source that will ultimately enhance the ability of the islands to naturally build backbarrier marsh, dunes, and a continuous sandy shoreline. Barrier island sediment nourishment should be executed with the understanding that gulf shoreline erosion is inevitable but that island area can be maintained and enhanced during retreat (thus significantly prolonging the life of the island chain) with strategic sand placement.

Louisiana↗
USGS2018

Evaluation of the Source and Transport of High Nitrate Concentrations in Ground Water, Warren Subbasin, California

Ground water historically has been the sole source of water supply for the Town of Yucca Valley in the Warren subbasin of the Morongo ground-water basin, California. An imbalance between ground-water recharge and pumpage caused ground-water levels in the subbasin to decline by as much as 300 feet from the late 1940s through 1994. In response, the local water district, Hi-Desert Water District, instituted an artificial recharge program in February 1995 using imported surface water to replenish the ground water. The artificial recharge program resulted in water-level recoveries of as much as 250 feet in the vicinity of the recharge ponds between February 1995 and December 2001; however, nitrate concentrations in some wells also increased from a background concentration of 10 milligrams per liter to more than the U.S. Environmental Protection Agency (USEPA) maximum contaminant level (MCL) of 44 milligrams per liter (10 milligrams per liter as nitrogen). The objectives of this study were to: (1) evaluate the sources of the high-nitrate concentrations that occurred after the start of the artificial-recharge program, (2) develop a ground-water flow and solute-transport model to better understand the source and transport of nitrates in the aquifer system, and (3) utilize the calibrated models to evaluate the possible effect of a proposed conjunctive-use project. These objectives were accomplished by collecting water-level and water-quality data for the subbasin and assessing changes that have occurred since artificial recharge began. Collected data were used to calibrate the ground-water flow and solute-transport models. Data collected for this study indicate that the areal extent of the water-bearing deposits is much smaller (about 5.5 square miles versus 19 square miles) than that of the subbasin. These water-bearing deposits are referred to in this report as the Warren ground-water basin. Faults separate the ground-water basin into five hydrogeologic units: the west, the midwest, the mideast, the east and the northeast hydrogeologic units. Water-quality analyses indicate that septage from septic tanks is the primary source of the high-nitrate concentrations measured in the Warren ground-water basin. Water-quality and stable-isotope data, collected after the start of the artificial recharge program, indicate that mixing occurs between imported water and native ground water, with the highest recorded nitrate concentrations in the midwest and the mideast hydrogeologic units. In general, the timing of the increase in measured nitrate concentrations in the midwest hydrogeologic unit is directly related to the distance of the monitoring well from a recharge site, indicating that the increase in nitrate concentrations is related to the artificial recharge program. Nitrate-to-chloride and nitrogen-isotope data indicate that septage is the source of the measured increase in nitrate concentrations in the midwest and the mideast hydrogeologic units. Samples from four wells in the Warren ground-water basin were analyzed for caffeine and selected human pharmaceutical products; these analyses suggest that septage is reaching the water table. There are two possible conceptual models that explain how high-nitrate septage reaches the water table: (1) the continued downward migration of septage through the unsaturated zone to the water table and (2) rising water levels, a result of the artificial recharge program, entraining septage in the unsaturated zone. The observations that nitrate concentrations increase in ground-water samples from wells soon after the start of the artificial recharge program in 1995 and that the largest increase in nitrate concentrations occur in the midwest and mideast hydrogeologic units where the largest increase in water levels occur indicate the validity of the second conceptual model (rising water levels). The potential nitrate concentration resulting from a water-level rise in the midwest and mideast hydrogeologic units was estimated using a simple mixing-cell model. The estimated value is within the range of concentrations measured in samples from wells, further indicating the validity of the second conceptual model. A ground-water flow model and a solute-transport model were developed for the Warren ground-water basin for the period 1956-2001. MODFLOW-96 was used for the ground-water flow model and MOC3D was used for the solute-transport model. The model cell size is about 500 feet by 500 feet and the models were discretized vertically into three layers. The models were calibrated using a trial-and-error approach using water-level and nitrate-concentration data collected between 1956-2001. In order to better match the measured data, low fault hydraulic characteristic values were required, thereby compartmentalizing the ground-water basin. In addition, it was necessary to parameterize the specific yield distribution for the top model layer where unconfined ground-water conditions occur into three homogeneous zones. Separate sets of specific- yield values were needed to simulate the drawdown and subsequent water-level recovery. In addition, the calibrated natural recharge was about 83 acre-feet per year. The entrainment of unsaturated-zone septage was simulated as recharge having an associated nitrate concentration. The volume of recharge was a function of the measured water-level rise between 1994-98 and the moisture content of the unsaturated zone. The nitrate concentration of the recharge water was a weighted function of the assumed nitrate concentration in the infiltrating water associated with the overlying land use. The simulated hydraulic head and nitrate concentration results were in good agreement with the measured data indicating that the mechanism for the increase in nitrate concentrations was rising water levels entraining high-nitrate septage in the unsaturated-zone. The calibrated models were used to simulate the possible effects of a planned conjunctive-use project in the western part of the ground-water basin. The simulated project included the addition of a new recharge pond and a new extraction well. In addition, recharge at two existing recharge ponds was increased and three existing production wells were pumped, treated in a nitrate-removal facility, and used for water supply. The simulated hydraulic heads increased in the west, the mideast, and parts of the east hydrogeologic units; however, the simulated hydraulic heads decreased in the midwest and northeast hydrogeologic units. The simulated nitrate concentrations increased to above the MCL of 44 milligrams per liter (10 milligrams per liter as nitrogen) in parts of the west as a result of the increase in simulated hydraulic head. The simulated nitrate concentrations decreased in part of the midwest hydrogeologic unit as a result of the artificial recharge and pumping from the nitrate-removal wells. The simulated nitrate concentrations increased to above the MCL of 44 milligrams per liter in part of the mideast and parts of the east hydrogeologic units beneath commercial land-use areas.

California↗