Geology topics
Kurt J. McCoy
Publications and source records attributed to Kurt J. McCoy.
Modeling groundwater-level responses to multiple stresses using transfer-function models and wavelet analysis in a coastal aquifer system
In coastal aquifers, dynamic stresses such as climate forcings, groundwater withdrawals, and ocean tidal fluctuations cause nonlinear responses to groundwater levels. Such responses to the stresses impact groundwater resources and related flooding and infrastructure risks at multiple scales. We used time-series models such as transfer-function models and wavelet analysis to quantify the relative contribution of these stresses to groundwater-level fluctuation in wells from the unconfined and confined aquifers in an Atlantic coastal aquifer. Climate forcings, such as precipitation and temperature, explained most of the groundwater-level variation for wells in the unconfined aquifer, whereas groundwater withdrawals were the dominant driver of groundwater levels for wells in the confined aquifer. The impact of groundwater withdrawals also was detected in several wells in the unconfined aquifer. Although the influence of ocean tides on groundwater levels commonly is observed in coastal aquifers, we found that daily groundwater withdrawals can obscure the semi-diurnal coherence signal of the two series. The magnitude of groundwater-level fluctuation that could be explained solely by tides was minor compared to that explained by climate or withdrawal stresses. Transfer-function modeling showed seasonal withdrawals from wells in confined aquifers had a significant, yet heterogeneous influence on groundwater levels in coastal aquifers, which highlights climate and withdrawals as key compounding stresses in coastal hydrology. This study demonstrates the value of time-series approaches to advance characterization of groundwater systems in areas with limited hydrogeologic parameter information.
GPS data from 2019 and 2020 campaigns in the Chesapeake Bay region towards quantifying vertical land motions
The Chesapeake Bay is a region along the eastern coast of the United States where sea-level rise is confounded with poorly resolved rates of land subsidence, thus new constraints on vertical land motions (VLM) in the region are warranted. In this paper, we provide a description of two campaign-style Global Positioning System (GPS) datasets, explain the methods used in data collection and validation, and present the experiment designed to quantify a new baseline of VLM in the Chesapeake Bay region of eastern North America. Data from GPS campaigns in 2019 and 2020 are presented as ASCII RINEX2.11 files and logsheets for each observation from the campaigns. Data were quality checked using the open-source program TEQC, resulting in average multipath 1 and 2 values of 0.68 and 0.57, respectively. All data are archived and publicly available for open access at the geodesy facility UNAVCO to abide by Findable, Accessible, Interoperable, Reusable (FAIR) data principles.
Groundwater-level contour map of Fauquier County, Virginia, October-November 2018
Groundwater withdrawals provide most public-water supplies and all private-domestic users in Fauquier County, Virginia, a fast-growing rural area southwest of Washington, D.C. Groundwater levels were measured in 129 wells during a county-wide synoptic survey from October 29 through November 2, 2018. Field measurements, combined with datapoints from the National Hydrography Dataset, were used to develop a county-wide groundwater-level contour map. Groundwater levels and withdrawals during the synoptic survey were near or slightly above long-term medians. Error analysis indicated that the estimated groundwater-level contours generally were lower than observed measurements, with a root-mean-squared error of 33.52 feet. Groundwater levels in Fauquier County are controlled largely by topography: low levels in the crystalline Blue Ridge aquifers in the northwestern part of the county contrast markedly with higher levels in the sedimentary Mesozoic Basin aquifers in the southeast. At current levels of groundwater withdrawal, and at the scale and scope of the synoptic survey, no cones of depression in the groundwater surface were detected. The Fauquier County groundwater-level contour map is available as a U.S. Geological Survey data release.
Groundwater hydrology and chemistry of Jamestown Island, Virginia—Potential effects of tides, storm surges, and sea-level rise on archaeological, cultural, and ecological resources
As the site of the first permanent English settlement in North America in 1607, Jamestown Island, Colonial National Historical Park (COLO), Virginia, contains a rich archaeological record that extends from the Paleoindian period (15,000 to 8,000 years ago) through the 20th century. The island is located on the lower James River near the mouth of Chesapeake Bay. Jamestown Island vegetation is dominated by upland forests surrounded by tidal, freshwater-to-oligohaline marshes. Along the Virginia coast, relative sea-level rise was more than 2.5 times the global average during the 20th century. Consequently, the National Park Service (NPS) has identified COLO as one of the 25 national parks most threatened by climate change. Surface waters across the island are hydraulically connected to the laterally continuous Surficial aquifer. The land-surface altitude of the island is low, with two-thirds of the island less than 5 feet (ft) above the North American Vertical Datum of 1988 (NAVD 88). Consequently, sea-level rise, combined with tides and storm surges, threatens the island and its resources as surface-water and groundwater levels rise, saltwater enters the Surficial aquifer, and groundwater chemistry changes. The impact of sea-level rise on the island’s surface-water resources has been well studied, but groundwater effects have been largely ignored. Quantifying the effects of tides, storm surges, and sea-level rise on groundwater levels and chemistry is essential to developing an effective strategy for managing climate-induced changes. The first step in developing a response strategy includes a parkwide general risk assessment for archaeological sites on the island, so that sites can be prioritized for management actions. The U.S. Geological Survey and the NPS began a study in 2015 to develop a long-term groundwater-monitoring program to evaluate this risk and to develop an updated management strategy. The groundwater-monitoring program consists of 45 wells and piezometers in two individual clusters and three transects across the island in different hydrologic and chemical settings. Samples for water quality were collected from the wells and piezometers from October 2015 through September 2018 at variable time intervals. Results of the monitoring identified disparate hydrologic and chemical responses to saltwater intrusion across the island. Specific conductance (an indicator of salinity) of groundwater beneath several marshes responded differently to changes in James River salinity. Groundwater response to changes in James River specific conductance appeared to be controlled by land-surface altitude and slope, differences in lateral and vertical sediment characteristics, distance from surface waters, and the degree of surface water/groundwater connectivity between channels and the aquifer. Groundwater chemistry data from monitoring wells at Black Point, a low-altitude, upland setting, are in contrast with conditions observed in Island House observation wells, a high-altitude, upland setting. Specific conductance (less than 200 microsiemens per centimeter [μS/cm]) and pH (greater than 5.0) of groundwater beneath much of the uplands that characterize the Island House observation wells are typical of groundwater in noncarbonate sedimentary aquifers recharged by precipitation. At Black Point, specific conductance ranged from 2,490 to 15,200 μS/cm, and pH ranged from 3.1 to 6.6 standard units. At the Black Point observation wells, the most saline and dense water was at the water table rather than deeper in the aquifer, causing a density inversion that persisted throughout the study. The density inversion likely resulted from differences in permeability between the shallow clay and fine-grained sands and the deeper coarse-grained sand and gravel. Groundwater with the lowest pH was at the water table. As saline groundwater flows through organic sediment beneath the marshes, bacterial biodegradation of organic matter creates anoxic conditions. Continued biodegradation concomitantly reduces iron-oxide minerals in the sediment and sulfate in saline water. When oxygen is reintroduced into groundwater, iron and sulfur can reoxidize to form sulfuric acid, locally lowering the pH of the water. This report describes the groundwater monitoring network design, rationale for site selection, monitoring approach, and results of monitoring from October 2015 through September 2018. Maps of inundation at selected water-level altitudes are included to identify the risk to archaeological, cultural, and ecological resources. The monitoring results of the hydrology and chemistry data are interpreted, and the different hydrologic and chemical settings are described. The implications of the study results for management decisions are presented, and suggestions for improving the monitoring network are included.
Documentation of a Soil-Water-Balance Model to estimate recharge to Blue Ridge, Piedmont, and Mesozoic Basin fractured-rock aquifers, Fauquier County, Virginia, 1996 through 2015
This report documents a Soil-Water-Balance (SWB) model that was developed for an area covering the Blue Ridge, Piedmont, and Mesozoic basin fractured-rock aquifers in Fauquier County, Virginia, for the calendar years 1996–2015. The SWB model includes an area of 1,498 square miles, divided into 1,076-square-foot (100-square-meter) grid cells on which daily groundwater recharge was estimated using existing elevation, meteorological, land-use, and soil property datasets. Daily groundwater recharge estimates obtained from the model were summarized annually, and annual model output was compared to the results of the hydrograph separation method, PART, on streamflow data from two streamgages in Fauquier County with periods of continuous record overlapping those of the SWB model period (01643700 Goose Creek near Middleburg, Virginia, and 01656000 Cedar Run near Catlett, Virginia). Spatially distributed groundwater recharge results from the SWB model represent annual conditions and the 20-year average values for the years 1996–2015, including estimated recharge during a previously defined drought in 2001. The 20-year average recharge in Fauquier County from the SWB model ranged from 8.1 inches per year (in/yr) in Blue Ridge aquifers to 5.3 in/yr in Mesozoic basin aquifers. Although mean annual precipitation volumes vary slightly across the County, the contrast in recharge among the Blue Ridge and western Piedmont aquifers with that of the Mesozoic basin aquifers is largely a result of differences in soil infiltration capacity. Precipitation totals 20 percent below mean annual precipitation from 1996–2015 produced drought recharge rates that were less than 50 percent of mean annual recharge. The SWB model and model output, including spatially distributed annual estimates of groundwater recharge, evapotranspiration, and gross precipitation for the 1996 through 2015 model period, are publicly available as a U.S. Geological Survey data release.
GIS-based method for estimating surficial groundwater levels in coastal Virginia using limited information
In many coastal areas, high water tables are present, complicating installation of some stormwater best management practices (BMPs) that rely on infiltration. Regional estimates of the seasonal high water table (SHWT) often rely on sources such as soil surveys taken over a decade ago; these data are static and do not account for groundwater withdrawals or other anthropogenic impacts. To improve estimates of the SHWT, we developed a GIS-based methodology relying on surface water elevations. Data sources included a 1.5-m (5.0 ft) resolution Lidar-derived digital elevation model (DEM), aerial imagery, and publicly available shapefiles of water boundaries. Twenty-six groundwater monitoring wells were screened to eliminate well locations influenced by pumping, yielding 22 wells. In coastal Virginia, tidal water bodies and ditches form terminal boundaries for discharge from the water-table aquifers and permit water table elevations to be fixed at the landward boundaries of surface water bodies. Water table elevations interpolated from well data and boundary elevations were used to create a triangulated irregular network representing the water table elevations for November 2012, which was the date of the DEM. An adjustment factor, calculated from the highest recorded April water table depth from long-term groundwater monitoring data, was added to estimate the SHWT elevation. SHWT elevations were subtracted from the DEM to yield SHWT depth, which was compared with long-term monitoring well data, yielding an R2 value of 0.91. Residual errors were random, although the method underpredicted the highest expected SHWT and overpredicted the median SHWT. The SHWT depth map was validated by using water table depths from 57 soil borings at 10 different sites, and consistently matched observations better than available soil survey estimates. The SHWT depth map could be useful for BMP siting and feasibility studies in similar hydrogeological settings.
The role of uplift and erosion in the persistence of saline groundwater in the shallow subsurface
In many regions of the world, the shallow (<300 m) subsurface is replenished with meteoric recharge within a few centuries or millennia, but in some regions saline groundwater persists despite abundant rainfall. Analyses of the flushing rate of shallow groundwater usually consider the permeability and recharge rate and a static landscape. The influence of landscape evolution can become important over millions of years, however. Here we present numerical simulations of fluid flow and transport in the top 1 km of a sedimentary foreland basin dominated by aquitards, where the rate of uplift and erosion (20 m Ma −1 ) balances that of meteoric flushing. Paleozoic age saline groundwater and brine persist at shallow depths that might otherwise have contained potable water. Similar hydrogeologic conditions, and uplift and erosion rates, likely exist in many other regions of the world, where a moving landscape has probably never been considered as an important contributor to groundwater quality.
Hydrogeology, groundwater flow, and groundwater quality of an abandoned underground coal-mine aquifer, Elkhorn Area, West Virginia
The Pocahontas No. 3 coal seam in southern West Virginia has been extensively mined by underground methods since the 1880’s. An extensive network of abandoned mine entries in the Pocahontas No. 3 has since filled with good-quality water, which is pumped from wells or springs discharging from mine portals (adits), and used as a source of water for public supplies. This report presents results of a three-year investigation of the geology, hydrology, geochemistry, and groundwater flow processes within abandoned underground coal mines used as a source of water for public supply in the Elkhorn area, McDowell County, West Virginia. This study focused on large (> 500 gallon per minute) discharges from the abandoned mines used as public supplies near Elkhorn, West Virginia. Median recharge calculated from base-flow recession of streamflow at Johns Knob Branch and 12 other streamflow gaging stations in McDowell County was 9.1 inches per year. Using drainage area versus mean streamflow relationships from mined and unmined watersheds in McDowell County, the subsurface area along dip of the Pocahontas No. 3 coal-mine aquifer contributing flow to the Turkey Gap mine discharge was determined to be 7.62 square miles (mi2), almost 10 times larger than the 0.81 mi2 surface watershed. Results of this investigation indicate that groundwater flows down dip beneath surface drainage divides from areas up to six miles east in the adjacent Bluestone River watershed. A conceptual model was developed that consisted of a stacked sequence of perched aquifers, controlled by stress-relief and subsidence fractures, overlying a highly permeable abandoned underground coal-mine aquifer, capable of substantial interbasin transfer of water. Groundwater-flow directions are controlled by the dip of the Pocahontas No. 3 coal seam, the geometry of abandoned mine workings, and location of unmined barriers within that seam, rather than surface topography. Seven boreholes were drilled to intersect abandoned mine workings in the Pocahontas No. 3 coal seam and underlying strata in various structural settings of the Turkey Gap and adjacent down-dip mines. Geophysical logging and aquifer testing were conducted on the boreholes to locate the coal- mine aquifers, characterize fracture geometry, and define permeable zones within strata overlying and underlying the Pocahontas No. 3 coal-mine aquifer. Water levels were measured monthly in the wells and showed a relatively static phreatic zone within subsided strata a few feet above the top of or within the Pocahontas No. 3 coal-mine aquifer (PC3MA). A groundwater-flow model was developed to verify and refine the conceptual understanding of groundwater flow and to develop groundwater budgets for the study area. The model consisted of four layers to represent overburden strata, the Pocahontas No. 3 coal-mine aquifer, underlying fractured rock, and fractured rock below regional drainage. Simulation of flow in the flooded abandoned mine entries using highly conductive layers or zones within the model, was unable to realistically simulate interbasin transfer of water. Therefore it was necessary to represent the coal-mine aquifer as an internal boundary condition rather than a contrast in aquifer properties. By representing the coal-mine aquifer with a series of drain nodes and optimizing input parameters with parameter estimation software, model errors were reduced dramatically and discharges for Elkhorn Creek, Johns Knob Branch, and other tributaries were more accurately simulated. Flow in the Elkhorn Creek and Johns Knob Branch watersheds is dependent on interbasin transfer of water, primarily from up dip areas of abandoned mine workings in the Pocahontas No. 3 coal-mine aquifer within the Bluestone River watershed to the east. For the 38th, 70th, and 87th percentile flow duration of streams in the region, mean measured groundwater discharge was estimated to be 1.30, 0.47, and 0.39 cubic feet per square mile (ft3/s/mi2
Groundwater hydrology and estimation of horizontal groundwater flux from the Rio Grande at selected locations in Albuquerque, New Mexico, 2009–10
The Albuquerque area of New Mexico has two principal sources of water: (1) groundwater from the Santa Fe Group aquifer system, and (2) surface water from the Rio Grande. From 1960 to 2002, pumping from the Santa Fe Group aquifer system caused groundwater levels to decline more than 120 feet while water-level declines along the Rio Grande in Albuquerque were generally less than 40 feet. These differences in water-level declines in the Albuquerque area have resulted in a great deal of interest in quantifying the river-aquifer interaction associated with the Rio Grande. In 2003, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, acting as fiscal agent for the Middle Rio Grande Endangered Species Collaborative Program, and the U.S. Army Corps of Engineers, began a study to characterize the hydrogeology of the Rio Grande inner valley alluvial aquifer in the Albuquerque area of New Mexico. The study provides hydrologic data in order to enhance the understanding of rates of water leakage from the Rio Grande to the alluvial aquifer, groundwater flow through the aquifer, and discharge of water from the aquifer to riverside drains. The study area extends about 20 miles along the Rio Grande in the Albuquerque area. Piezometers and surface-water gages were installed in paired transects at eight locations. Nested piezometers, completed at various depths in the alluvial aquifer, and surface-water gages, installed in the Rio Grande and riverside drains, were instrumented with pressure transducers. Water-level and water-temperature data were collected from 2009 to 2010. Water levels from the piezometers indicated that groundwater movement was usually away from the river towards the riverside drains. Annual mean horizontal groundwater gradients in the inner valley alluvial aquifer ranged from 0.0024 (I-25 East) to 0.0144 (Pajarito East). The median hydraulic conductivity values of the inner valley alluvial aquifer, determined from slug tests, ranged from 30 feet per day (ft/d) (Montaño) to 120 ft/d (Central) for paired transects, with a median hydraulic conductivity for all transects of 50 ft/d. Daily mean groundwater fluxes from the river through the inner valley alluvial aquifer computed using Darcy’s Law and the slug test results ranged from about 0.01 ft/d (Montaño West) to between 1.0 and 2.0 ft/d (Central East). Median annual groundwater fluxes from the river through the inner valley alluvial aquifer determined using the Suzuki-Stallman method was greatest at Alameda East (0.50 ft/d) and lowest at Alameda West (0.25 ft/d). The results from both methods agreed reasonably well. Seepage investigations conducted by measuring discharge in the east and west riverside drains provided information for computing changes in flow within the drains and for evaluating results from Darcy’s Law and Suzuki-Stallman method flux calculations. Discharge measured in the east riverside drain between the Barelas Bridge and the I-25 bridge indicated that the flow in the east riverside drain increased by an average of 56.5 cubic feet per day per linear foot (ft 3 /d/ft) of drain. Discharge measured in the west riverside drain between the Central bridge and the I-25 bridge indicated that flow increased between west drain miles 0 and 4, an average of 53.8 ft 3 /d/ft of drain, and that flow increased between west drain miles 7 and 10, an average of 44.9 ft 3 /d/ft of drain. In comparison to the seepage measurement results, the groundwater fluxes from the river through the inner valley alluvial aquifer calculated from Darcy’s Law (q slug ) and by the Suzuki-Stallman method (q heat ) would account for 20–36 percent or 53–95 percent, respectively, of the total flow in the east riverside drain and 22–31 percent or 19–26 percent, respectively, of the total flow in the west drain. These results indicate that the drains likely also receive water from outside the inner valley. The spatial variability of horizontal hydraulic gradients and groundwater fluxes can be primarily attributed to variability in the distances between the river and riverside drains throughout the study area and geologic heterogeneities in the alluvial aquifer. Temporal variability in the water levels, which control the horizontal hydraulic gradients and fluxes between the Rio Grande and the riverside drains, can be primarily attributed to seasonal fluctuations in river stage and irrigation practices.
Hydrogeology and simulation of groundwater flow in fractured-rock aquifers of the Piedmont and Blue Ridge Physiographic Provinces, Bedford County, Virginia
An annual groundwater budget was computed as part of a hydrogeologic characterization and monitoring effort of fractured-rock aquifers in Bedford County, Virginia, a growing 764-square-mile (mi 2 ) rural area between the cities of Roanoke and Lynchburg, Virginia. Data collection in Bedford County began in the 1930s when continuous stream gages were installed on Goose Creek and Big Otter River, the two major tributaries of the Roanoke River within the county. Between 2006 and 2014, an additional 2 stream gages, 3 groundwater monitoring wells, and 12 partial-record stream gages were operated. Hydrograph separation methods were used to compute base-flow recharge rates from the continuous data collected from the continuous stream gages. Mean annual base-flow recharge ranged from 8.3 inches per year (in/yr) for the period 1931–2012 at Goose Creek near Huddleston (drainage area 188 mi 2 ) to 9.3 in/yr for the period 1938–2012 at Big Otter River near Evington (drainage area 315 mi 2 ). Mean annual base-flow recharge was estimated to be 6.5 in/yr for the period 2007–2012 at Goose Creek at Route 747 near Bunker Hill (drainage area 125 mi 2 ) and 8.9 in/yr for the period 2007–2012 at Big Otter River at Route 221 near Bedford (drainage area 114 mi 2 ). Base-flow recharge computed from the partial-record data ranged from 5.0 in/yr in the headwaters of Goose Creek to 10.5 in/yr in the headwaters of Big Otter River. A steady-state groundwater-flow simulation for Bedford County was developed to test the conceptual understanding of flow in the fractured-rock aquifers and to compute a groundwater budget for the four major drainages: James River, Smith Mountain and Leesville Lakes, Goose Creek, and Big Otter River. Model results indicate that groundwater levels mimic topography and that minimal differences in aquifer properties exist between the Proterozoic basement crystalline rocks and Late Proterozoic-Cambrian cover crystalline rocks. The Big Otter River receives 40.8 percent of the total daily groundwater outflow from fractured-rock aquifers in Bedford County; Goose Creek receives 25.8 percent, the James River receives 18.2 percent, and Smith Mountain and Leesville Lakes receive 15.2 percent. The remaining percentage of outflow is attributed to pumping from the aquifer (consumptive use).
Hydrologic budget and conditions of Permian, Pennsylvanian, and Mississippian aquifers in the Appalachian Plateaus physiographic province
In response to challenges to groundwater availability posed by historic land-use practices, expanding development of hydrocarbon resources, and drought, the U.S. Geological Survey Groundwater Resources Program began a regional assessment of the Appalachian Plateaus aquifers in 2013 that incorporated a hydrologic landscape approach to estimate all components of the hydrologic system: surface runoff, base flow from groundwater, and interaction with atmospheric water (precipitation and evapotranspiration). This assessment was intended to complement other Federal and State investigations and provide foundational groundwater-related datasets in the Appalachian Plateaus. A regional Soil-Water-Balance model was constructed for a 160,000-square-mile study area that extended to the topographic divide of all streams originating outside but flowing into areas underlain by Appalachian Plateaus aquifers. The model incorporated soil, landscape, and climate variables to estimate an annual water budget for the 32-year period from 1980 to 2011 and was calibrated using base-flow data estimated by hydrograph separation techniques from 20 streamflow gaging stations across the study area. Over this period, an average of 47 inches per year (in/yr) of precipitation fell on Appalachian Plateaus aquifers. Simulations from the regional Soil-Water-Balance model indicate that only 19 percent of the precipitation or an average 9 in/yr recharged aquifers, and 19 percent resulted in surface runoff to streams. The remaining 62 percent, an average of 27 in/yr of water, was returned to the atmosphere via evapotranspiration. Because withdrawals from aquifers due to pumping equated to less than 1 percent of the water budget, differences in predevelopment and postdevelopment regional water budgets of the Appalachian Plateaus were minimal. Storage changes caused by filling of abandoned coal-mine aquifers and long-term differences in aquifer storage resulting from climate fluctuations constitute a small portion of the overall water budget. The percentage of precipitation that results in recharge, runoff, or evapotranspiration from the landscape varies annually by up to a factor of two depending on temporal changes in prevailing climate conditions and spatial changes in basin characteristics, precipitation patterns, and sources of atmospheric moisture over a large study area. A comparison of water-budget estimates from the regional Soil-Water-Balance model for a dry year (1988) and wet year (2004) showed that evapotranspiration accounts for most of the annual differences in precipitation. As a portion of annual precipitation, evapotranspiration ranged from 69 percent (dry year) to 52 percent (wet year), a range four times greater than the 15 percent (dry year) to 18 percent (wet year) range estimated for recharge. Evapotranspiration as a percentage of precipitation peaks during dry periods, whereas base flow and runoff tend to reach minimum values. During wet periods, this relationship is reversed and base flow and runoff as a percentage of precipitation generally peak while evapotranspiration percentages reach minimum values. Annual recharge in the Appalachian Plateaus reaches a maximum at near 20 percent of annual precipitation, regardless of the severity of wet conditions. Hydrograph separation data from 849 streamflow gaging stations in the study area were used to assess trends in streamflow, base flow, surface runoff, and base-flow index, or ratio of base flow to streamflow, in the Appalachian Plateaus for the period from 1930 to 2011. Annual data anomalies for each of the four variables were individually defined as the annual standard deviation from the mean at all 849 streamflow gaging stations. Annual data anomalies confirm the close relation of annual precipitation to both base flow and runoff components of streamflow, and both components increased during the period of analysis. Around 1970, conditions shifted streamflow from values generally below to above long-term means. At a regional scale, increases in base flow account for most of these observed increases in mean annual streamflow. The independence of the base-flow index to annual climate trends indicate that changes in the components of streamflow of the Appalachian Plateaus are probably in response to shifts in seasonal precipitation or widespread land-use practices. A subset of 77 index streamgages, defined as having 60 or more years of complete record between the years 1930 and 2011 with no more than 20 percent missing data, was selected to show spatial patterns of change in the water budget. Data from the index streamgages showed that the overall trends in base flow are dependent upon the period of evaluation. Long-term (1930–2011) increases in base flow were observed throughout the study area. For two shorter periods (1930–1969 and 1970–2011) trends in base flow were largely negative. In general, spatial patterns of change in streamflow, base flow, and runoff were mixed but generally consistent with prevailing climate patterns and land-use changes.
Annual and average estimates of water-budget components based on hydrograph separation and PRISM precipitation for gaged basins in the Appalachian Plateaus Region, 1900-2011
As part of the U.S. Geological Survey’s Groundwater Resources Program study of the Appalachian Plateaus aquifers, annual and average estimates of water-budget components based on hydrograph separation and precipitation data from parameter-elevation regressions on independent slopes model (PRISM) were determined at 849 continuous-record streamflow-gaging stations from Mississippi to New York and covered the period of 1900 to 2011. Only complete calendar years (January to December) of streamflow record at each gage were used to determine estimates of base flow, which is that part of streamflow attributed to groundwater discharge; such estimates can serve as a proxy for annual recharge. For each year, estimates of annual base flow, runoff, and base-flow index were determined using computer programs—PART, HYSEP, and BFI—that have automated the separation procedures. These streamflow-hydrograph analysis methods are provided with version 1.0 of the U.S. Geological Survey Groundwater Toolbox, which is a new program that provides graphing, mapping, and analysis capabilities in a Windows environment. Annual values of precipitation were estimated by calculating the average of cell values intercepted by basin boundaries where previously defined in the GAGES–II dataset. Estimates of annual evapotranspiration were then calculated from the difference between precipitation and streamflow.
Assessing the magnitude and timing of anthropogenic warming of a shallow aquifer: example from Virginia Beach, USA
Groundwater temperature measurements in a shallow coastal aquifer in Virginia Beach, Virginia, USA, suggest groundwater warming of +4.1 °C relative to deeper geothermal gradients. Observed warming is related to timing and depth of influence of two potential thermal drivers—atmospheric temperature increases and urbanization. Results indicate that up to 30 % of groundwater warming at the water table can be attributed to atmospheric warming while up to 70 % of warming can be attributed to urbanization. Groundwater temperature readings to 30-m depth correlate positively with percentage of impervious cover and negatively with percentage of tree canopy cover; thus, these two land-use metrics explain up to 70 % of warming at the water table. Analytical and numerical modeling results indicate that an average vertical groundwater temperature profile for the study area, constructed from repeat measurement at 11 locations over 15 months, is consistent with the timing of land-use change over the past century in Virginia Beach. The magnitude of human-induced warming at the water table (+4.1 °C) is twice the current seasonal temperature variation, indicating the potential for ecological impacts on wetlands and estuaries receiving groundwater discharge from shallow aquifers.
Hydrogeology, water resources, and water budget of the upper Rio Hondo Basin, Lincoln County, New Mexico, 2010
The upper Rio Hondo Basin occupies a drainage area of 585 square miles in south-central New Mexico and comprises three general hydrogeologic terranes: the higher elevation “Mountain Block,” the “Central Basin” piedmont area, and the lower elevation “Hondo Slope.” As many as 12 hydrostratigraphic units serve as aquifers locally and form a continuous aquifer on the regional scale. Streams and aquifers in the basin are closely interconnected, with numerous gaining and losing stream reaches across the study area. In general, the aquifers are characterized by low storage capacity and respond to short-term and long-term variations in recharge with marked water-level fluctuations on short (days to months) and long (decadal) time scales. Droughts and local groundwater withdrawals have caused marked water-table declines in some areas, whereas periodically heavy monsoons and snowmelt events have rapidly recharged aquifers in some areas. A regional-scale conceptual water budget was developed for the study area in order to gain a basic understanding of the magnitude of the various components of input, output, and change in storage. The primary input is watershed yield from the Mountain Block terrane, supplying about 38,200 to 42,300 acre-feet per year (acre-ft/yr) to the basin, as estimated by comparing the residual of precipitation and evapotranspiration with local streamgage data. Streamflow from the basin averaged about 21,200 acre-ft/yr, and groundwater output left the basin at an estimated 2,300 to 5,700 acre-ft/yr. The other major output (about 13,500 acre-ft/yr) was by public water supply, private water supply, livestock, commercial and industrial uses, and the Bonito Pipeline. The residual in the water budget, the difference between the totals of the input and output terms or the potential change in storage, ranged from -2,200 acre-ft/yr to +5,300 acre-ft/yr. There is a high degree of variability in precipitation and consequently in the water supply; small variations in annual precipitation can result in major changes in overall watershed yield. Changing water-use patterns, concentrated areas of groundwater withdrawal, and variations in precipitation have created localized areas where water-table declines and diminished surface flow are of concern.
Hydrogeology and hydrologic conditions of the Northern Atlantic Coastal Plain aquifer System from Long Island, New York, to North Carolina
The seaward-dipping sedimentary wedge that underlies the Northern Atlantic Coastal Plain forms a complex groundwater system. This major source of water provides for public and domestic supply and serves as a vital source of freshwater for industrial and agricultural uses throughout the region. Population increases and land-use and climate changes, however, have led to competing demands for water. The regional response of the aquifer system to these stresses poses regional challenges for water-resources management at the State level because hydrologic effects often extend beyond State boundaries. In response to these challenges, the U.S. Geological Survey Groundwater Resources Program began a regional assessment of the groundwater availability of the Northern Atlantic Coastal Plain aquifer system in 2010. The initial phase of this investigation included a refinement of the hydrogeologic framework and an updated hydrologic budget of this aquifer system from the last regional aquifer system assessment completed by the U.S. Geological Survey in the 1980s. Refinements to the hydrogeologic framework include revision of the regional aquifer names to be more consistent with local names in New York, New Jersey, Delaware, Maryland, and Virginia, the primary States included in the study area. Other revisions to the framework include characterization of the aquifers of the regional Potomac aquifer system. The regional Potomac aquifer system is subdivided for this report into two regional aquifers. These aquifers include the single Potomac aquifer in Virginia and two aquifers in Maryland, Delaware, and New Jersey, where the Potomac aquifer system thickens within the Salisbury Embayment. The two regional aquifers making up the Potomac aquifer system include the Potomac-Patapsco aquifer and the underlying Potomac-Patuxent aquifer. The Potomac-Patuxent aquifer includes the Lower Potomac-Raritan-Magothy aquifer in southern New Jersey and the Patuxent aquifers in Delaware and Maryland. In northern New Jersey and on Long Island, New York, the PotomacPatuxent aquifer is absent, but the Late Cretaceous fluvialdeltaic aquifer that is laterally equivalent with the upper part of the Potomac Formation now is considered part of the regional Potomac-Patapsco aquifer. This aquifer includes the Middle Potomac-Raritan-Magothy aquifer in New Jersey and the Lloyd aquifer on Long Island. The name “Upper Potomac aquifer” has been removed as part of this regional framework revision. The local aquifer previously considered part of the Upper Potomac aquifer now are part of the regional Magothy aquifer. These units include the Upper Potomac-Raritan-Magothy aquifer in New Jersey, the Magothy aquifers on Long Island, Delaware, and Maryland, and the Virginia Beach aquifer in Virginia. Updates to the regional hydrologic budget include revised estimates of aquifer recharge, water use and streamflow data. Inflow to the aquifer system of about 20,000 million gallons per day (Mgal/d) includes 19,600 Mgal/d from recharge from precipitation, 200 Mgal/d of recharge from wastewater via onsite domestic septic systems, and 200 Mgal/d from the release of water from aquifer storage. Outflow from the aquifer system includes groundwater discharge to streams (11,900 Mgal/d), groundwater withdrawals (1,500 Mgal/d), and groundwater discharge to coastal waters (6,600 Mgal/d). A numerical modeling analysis is required to improve this hydrologic budget calculation and to forecast future changes in water levels and aquifer storage caused by groundwater withdrawals, land-use changes, and the effects of climate variability and change.
Groundwater hydrology and estimation of horizontal groundwater flux from the Rio Grande at selected locations in Albuquerque, New Mexico, 2003-9
The Albuquerque, New Mexico, area has two principal sources of water: groundwater from the Santa Fe Group aquifer system and surface water from the San Juan-Chama Diversion Project. From 1960 to 2002, groundwater withdrawals from the Santa Fe Group aquifer system have caused water levels to decline more than 120 feet in some places within the Albuquerque area, resulting in a great deal of interest in quantifying the river-aquifer interaction associated with the Rio Grande. In 2003, the U.S. Geological Survey in cooperation with the Bureau of Reclamation, the Middle Rio Grande Endangered Species Collaborative Program, and the U.S. Army Corps of Engineers began a detailed characterization of the hydrogeology of the Rio Grande riparian corridor in the Albuquerque, New Mexico, area to provide hydrologic data and enhance the understanding of rates of water leakage from the Rio Grande to the alluvial aquifer, groundwater flow through the aquifer, and discharge of water from the aquifer to the riverside drains. A simple conceptual model of flow indicates that the groundwater table gently slopes from the Rio Grande towards riverside drains and the outer boundaries of the inner valley. Water infiltrating from the Rio Grande initially moves vertically below the river, but, as flow spreads farther into the Rio Grande inner valley alluvial aquifer, flow becomes primarily horizontal. The slope of the water-table surface may be strongly controlled by the riverside drains and influenced by other more distal hydrologic boundary conditions, such as groundwater withdrawals by wells. Results from 35 slug tests performed in the Rio Grande inner valley alluvial aquifer during January and February 2009 indicate that hydraulic-conductivity values ranged from 5 feet per day to 160 feet per day with a median hydraulic-conductivity for all transects of 40 feet per day. Median annual horizontal hydraulic gradients in the Rio Grande inner valley alluvial aquifer ranged from 0.011 to 0.002. Groundwater fluxes through the alluvial aquifer calculated by using median slug-test results (qm slug ) and Darcy's law ranged from about 0.1 feet per day to about 0.7 feet per day. Groundwater fluxes calculated by using the Suzuki-Stallman method (qm heat ) ranged from 0.52 feet per day to 0.23 feet per day. Results from the Darcy's law and Suzuki-Stallman flux calculations were compared to discharge measured in riverside drains on both sides of the river north of the Montaño Bridge on February 26, 2009. Flow in the Corrales Riverside Drain increased by 1.4 cubic feet per second from mile 2 to mile 4, about 12 cubic feet per day per linear foot of drain. Flow in the Albuquerque Riverside Drain increased by 15 cubic feet per second between drain miles 0 and 3, about 82 cubic feet per day per linear foot of drain. The flux of water from the river to the aquifer was calculated to be 2.2 cubic feet per day per linear foot of river by using the median qm slug of 0.09 feet per day at Montaño transects west of the river. The total flux was calculated to be 6.0 cubic feet per day per linear foot of river by using the mean(qm heat of 0.24 feet per day for the Montaño transects west of the river. Assuming the Corrales Riverside Drain intercepted all of this flow, the qm slug or qm heat fluxes account for 18 to 50 percent, respectively, of the increase of flow in the drain. The flux of water from the river to the aquifer was calculated to be 15 cubic feet per day per linear foot of river by using the median qm slug of 0.30 feet per day at the Montaño transects east of the river. The flux of water from the river to the aquifer was calculated to be 17 cubic feet per day per linear foot of river by using the mean flux calculated from the Suzuki-Stallman method for the Montaño East transects of 0.34 feet per day. Assuming the Albuquerque Riverside Drain intercepted all this flow, the qm slug or (qm heat fluxes would only account for 18 to 21 percent, respectively, of the increase in flow in the drain. The comparison of these results with those of previous investigations suggests that calculated flux through the Rio Grande inner valley alluvial aquifer is strongly scale dependent and that the thickness of aquifer through which river water flows may be greater than indicated by the vertical temperature profiles.
Use of modflow drain package for simulating inter-basin transfer in abandoned coal mines
Simulation of groundwater flow in abandoned mines is difficult, especially where flux to and from mines is unknown or poorly quantified, and inter-basin transfer of groundwater occurs. A 3-year study was conducted in the Elkhorn area, West Virginia to better understand groundwater-flow processes and inter-basin transfer in above drainage abandoned coal mines. The study area was specifically selected, as all mines are located above the elevation of tributary receiving streams, to allow accurate measurements of discharge from mine portals and tributaries for groundwater model calibration. Abandoned mine workings were simulated in several ways, initially as a layer of high hydraulic conductivity bounded by lower permeability rock in adjacent strata, and secondly as rows of higher hydraulic conductivity embedded within a lower hydraulic conductivity coal aquifer matrix. Regardless of the hydraulic conductivity assigned to mine workings, neither approach to simulate mine workings could accurately reproduce the inter-basin transfer of groundwater from adjacent watersheds. To resolve the problem, a third approach was developed. The MODFLOW DRAIN package was used to simulate seepage into and through mine workings discharging water under unconfined conditions to Elkhorn Creek, North Fork, and tributaries of the Bluestone River. Drain nodes were embedded in a matrix of uniform hydraulic conductivity cells that represented the coal mine aquifer. Drain heads were empirically defined from well observations, and elevations were based on structure contours for the Pocahontas No. 3 mine workings. Use of the DRAIN package to simulate mine workings as an internal boundary condition resolved the inter-basin transfer problem, and effectively simulated a shift from a topographic- dominated to a dip-dominated flow system, by dewatering overlying unmined strata and shifting the groundwater drainage divide up dip within the Pocahontas No. 3 coal seam several kilometers into the adjacent Bluestone River Watershed. Model simulations prior to use of the DRAIN package for simulating mine workings produced estimated flows of 0.32 to 0.34 m3/s in each of the similar sized Elkhorn Creek and North Fork Watersheds, but failed to estimate inter-basin transfer of groundwater from the adjacent Bluestone River Watershed. The simulation of mine entries and discharge using the MODFLOW DRAIN package produced estimated flows of 0.46 and 0.26 m3/s for the Elkhorn Creek and North Fork watersheds respectively, which matched well measured flows for the respective watersheds of 0.47 and 0.26 m3/s.