Geology Reports⌕ Search

USGS · 70273478

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Abstract

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system. Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer currently is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from areas outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2025, rates of withdrawal substantially lower than the historical maximum have resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or in adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up over 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Jason P. Pope, Alison D. Gordon, Ryan S. Frederiks. 2026-01-15. Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system. https://doi.org/10.31223/x5hb5d

Cite the original work for its findings. Save a collection to share your selection of sources.

Discover connections

Connections use source metadata and explicit phrase matches, not verified experimental comparisons.

KEEP EXPLORING

Related USGS reports

Evaluating bank restoration and stabilization in the Yellow River in northern Indiana: Mitigation of sediment and nutrient transport

The U.S. Geological Survey, in cooperation with the Kankakee River Basin and Yellow River Basin Development Commission, evaluated the effects of streambank remediation on the Yellow River in northwest Indiana. Two phases of remediation were completed on the Yellow River between 2021 and 2023, after assessments indicated it was a large contributor of sediment to the Kankakee River. Temporal effects were assessed using pre‑remediation (April 2019 –October 2021) and post‑remediation (February – October 2024) data from two long‑term monitoring sites, USGS 05516665 Yellow River near Oak Grove (Oak Grove) and USGS 05517000 Yellow River at Knox (Knox). Spatial patterns were evaluated among sites from upstream to downstream during the post-remediation period. Monitoring included the pre-existing Oak Grove and Knox sites, an additional monitoring site at USGS streamgage 05516500 Yellow River at Plymouth (Plymouth), and discrete sampling at two additional sites. Discrete sites include USGS 05516663 Yellow River at Upas Rd near Oak Grove (Upas Rd), located upstream of the remediation reach, and USGS 05516667 Yellow River at CR11100E near Oak Grove (CR1100E), located within the remediation reach. Generalized least squares and linear mixed‑effects regression models were used to evaluate temporal differences in turbidity-to-discharge relations for slope and turbidity units, discharge, and suspended‑sediment in before and after remediation data. Similar analysis using generalized least squares and linear mixed‑effects regression models were conducted for spatial analysis of discharge and sediment and nutrient concentrations and fluxes at each site post-remediation. Temporal comparisons showed that turbidity increased less with discharge increases after remediation, indicated by slope of relation, at both Oak Grove and Knox, indicating improved water clarity during higher flows. Mean turbidity and suspended‑sediment concentrations were not statistically different between remediation periods. Comparison of the post-remediation continuous water quality data from upstream to downstream sites through the remediated reaches indicated that downstream sites exhibited steeper turbidity‑to‑discharge slopes than the upstream control site at Plymouth. A statistically significant difference was detected between the slopes at Oak Grove and Plymouth, whereas Knox did not differ significantly from either Plymouth or Oak Grove. Downstream sites had significantly lower mean turbidity, suggesting improved water clarity within and downstream of the remediation reach. Nutrient concentrations generally decreased moving downstream and suspended‑sediment concentrations showed no consistent trend. Nutrient and sediment fluxes increased downstream, likely reflecting the influence of greater discharge volumes. Interpretation of remediation effects was limited by the absence of pre‑remediation data at upstream sites and by the short post‑remediation monitoring period. Continued monitoring could allow for assessment of longer‑term responses to streambank remediation on the Yellow River, as effects may become more pronounced over time.

EarthArXiv↗

Hydrogeology and water quality of the unconsolidated aquifer system in the Lower Fall Creek Valley, Tompkins County, New York

The U.S. Geological Survey, in cooperation with Tompkins County and the Town of Dryden, conducted a hydrogeologic investigation of the lower Fall Creek unconsolidated aquifer system in Tompkins County, New York. This study is part of an ongoing effort to characterize major valley-fill aquifers in the region to support sustainable groundwater management. The lower Fall Creek valley, encompassing the Village of Freeville, the unincorporated communities of Etna, and the hamlets of McLean and Varna, has experienced moderate development and relies heavily on groundwater for residential, agricultural, and commercial use. The study area spans approximately 51 square miles along the boundary of the Allegheny Plateau and North Rim of the Allegheny Plateau physiographic provinces. The aquifer system, composed primarily of glacially derived sand and gravel deposits, includes both unconfined and confined aquifer units. For the approximately 2,300 people who live over the aquifers, groundwater is the sole source of drinking water. In some areas, such as the Village of Freeville, artesian conditions exist, with potentiometric heads exceeding 25 feet above land surface. This report describes the hydrogeologic framework, aquifer distribution, groundwater-flow dynamics, and water quality of the lower Fall Creek aquifer system. Groundwater is generally hard and contains elevated concentrations of redox-sensitive constituents, such as iron and manganese relative to U.S. Environmental Protection Agency secondary maximum contaminant levels. Arsenic was not detected above the U.S. Environmental Protection Agency maximum contaminant level during this study but remains a regional concern due to past detections in the study area and nearby aquifers. The findings provide data for local planners and resource managers to use in understanding and protecting this water resource.

EarthArXiv↗

Investigation of drivers and sources of potential harmful algal blooms in Lake Anna, Virginia, 2023-24

Algal blooms in freshwater reservoirs are an increasing concern because of the risk to aquatic ecosystem health, recreational use, and water quality. Among the taxa that comprise an algal bloom, cyanobacteria are of particular concern due to the potential to produce toxic compounds which can cause acute and chronic illness in humans if ingested through contaminated water or shellfish. When toxins are above a toxin-specific threshold in the water column, the Commonwealth of Virginia categorizes the bloom as a harmful algal bloom (HAB), distinguishing it from an algal bloom without toxins present. Toxins are of particular concern in bodies of water which serve as major recreational areas, such as Lake Anna in central Virginia. Lake Anna is a 27 km-long and 47 km2 reservoir in central Virginia. Since the Commonwealth of Virginia revised the HAB monitoring framework in 2018, algal bloom advisories for potentially toxigenic blooms have been issued in the western, riverine zone of the reservoir. There were detections of microcystin and anatoxin production but no detections of cylindrospermopsin and saxitoxin production in Lake Anna; however, toxin concentrations were low, at values less than 0.65 µg/L. In 2023, the U.S. Geological Survey, in cooperation with Virginia Department of Environmental Quality, initiated a 19-month intensive and extensive study of western Lake Anna and its tributaries (North Anna River and Pamunkey Creek) to (1) characterize the algal community, biomass, and toxin production; (2) identify potential drivers leading to algal bloom initiation, persistence, and decline; and (3) assess the primary sources of the drivers contributing to algal blooms. Continuous and discrete monitoring of water quality within the lake, watershed inputs, meteorological conditions, and algal communities were monitored. Multiple modeling scenarios identified that the formation and persistence of algal blooms in Lake Anna are driven by interactions among chemical, macro- and micronutrient, and physical factors. Key drivers include total nitrogen, total phosphorus, water temperature, suspended sediment concentration, wind speed, copper, iron, molybdenum, nickel, sodium, and alkalinity. During this study, winter stormflow events delivered the highest loading of macro- and micronutrients from the tributaries to the lake. The lake is a sink for these nutrients in the water column or in the lake-bed sediments. During high-occupancy and recreation periods on the lake in the summer months, boating activity can promote internal loading by the resuspension of the lake-bed sediments, thereby increasing the nutrient availability for algal uptake. Thermal stratification and hypoxic conditions were observed, which can contribute to the passive release of nutrients into the water column, further facilitating algal bloom proliferation. These processes may create a positive feedback loop that supports the formation, persistence, and annual recurrence of algal blooms in the western region of Lake Anna. Continued water-quality monitoring of early warning indicators may improve the detection and forecasting of algal growth in the lake and help managers proactively manage water quality and algal growth.

EarthArXiv↗