USGS Science⌕ Search

SEARCH · USGS Science

Results for “Groundwater”

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 865 records · Page 48Linked to original sources
USGS2007

Role of aquifer heterogeneity in fresh groundwater discharge and seawater recycling: An example from the Carmel coast, Israel

A case study is shown in which the pattern of submarine groundwater discharge and of seawater recycling is controlled by local hydrogeological variability. The coastal aquifer in Dor Bay is composed of two units: a partly confined calcaranitic sandstone (Kurkar) and an overlying loose sand. Groundwater in the Kurkar has elevated activities of 222 Rn (∼390 dpm/L) and relatively low 224 Ra/ 223 Ra activity ratios (3–4), while the sand groundwater is significantly less radiogenic (6–90 dpm/L) and shows higher 224 Ra/ 223 Ra ratios. Groundwater discharging from sand-covered areas of the bay has salinities of 16–31 and an average 222 Rn activity of 168 dpm/L, which lies on a mixing line between Rn-rich Kurkar fresh water and Rn-poor seawater. Another key observation is that seawater infiltrates to some extent into onshore sand groundwater, while the fresh water within the submarine Kurkar can be traced up to 40 m offshore. This implies that while fresh water mainly discharges from the Kurkar unit, seawater recycling is limited to the loose sand, and that the discharge from sand-covered areas is a mixture of Kurkar water with recycled seawater. Advection rates from the bay floor were calculated from Rn time series and found to vary between 0 and 36 cm/d, correlating negatively with bay water depth. The average flux was 8.1 cm/d, and it did not seem to change much during March, May, and July 2006. The average amount of fresh water discharging to the bay was 5.0 m 3 /d per meter of shoreline. Radon activity in the sand groundwater also fluctuates due to influx of Kurkar-type groundwater.

Carmel Coast, Dor Bay↗
USGS2008

Influence of perched groundwater on base flow

Analysis with a three‐dimensional variably saturated groundwater flow model provides a basic understanding of the interplay between streams and perched groundwater. A simplified, layered model of heterogeneity was used to explore these relationships. Base flow contribution from perched groundwater was evaluated with regard to varying hydrogeologic conditions, including the size and location of the fine‐sediment unit and the hydraulic conductivity of the fine‐sediment unit and surrounding coarser sediment. Simulated base flow was sustained by perched groundwater with a maximum monthly discharge in excess of 15 L/s (0.6 feet 3 /s) over the length of the 2000‐m stream reach. Generally, the rate of perched‐groundwater discharge to the stream was proportional to the hydraulic conductivity of sediment surrounding the stream, whereas the duration of discharge was proportional to the hydraulic conductivity of the fine‐sediment unit. Other aspects of the perched aquifer affected base flow, such as the depth of stream penetration and the size of the fine‐sediment unit. Greater stream penetration decreased the maximum base flow contribution but increased the duration of contribution. Perched groundwater provided water for riparian vegetation at the demand rate but reduced the duration of perched‐groundwater discharge nearly 75%.

Water Resources Research↗
USGS2012

Influence of permafrost distribution on groundwater flow in the context of climate-driven permafrost thaw: Example from Yukon Flats Basin, Alaska, United States

Understanding the role of permafrost in controlling groundwater flow paths and fluxes is central in studies aimed at assessing potential climate change impacts on vegetation, species habitat, biogeochemical cycling, and biodiversity. Recent field studies in interior Alaska show evidence of hydrologic changes hypothesized to result from permafrost degradation. This study assesses the hydrologic control exerted by permafrost, elucidates modes of regional groundwater flow for various spatial permafrost patterns, and evaluates potential hydrologic consequences of permafrost degradation. The Yukon Flats Basin (YFB), a large (118,340 km 2 ) subbasin within the Yukon River Basin, provides the basis for this investigation. Model simulations that represent an assumed permafrost thaw sequence reveal the following trends with decreasing permafrost coverage: (1) increased groundwater discharge to rivers, consistent with historical trends in base flow observations in the Yukon River Basin, (2) potential for increased overall groundwater flux, (3) increased spatial extent of groundwater discharge in lowlands, and (4) decreased proportion of suprapermafrost (shallow) groundwater contribution to total base flow. These trends directly affect the chemical composition and residence time of riverine exports, the state of groundwater-influenced lakes and wetlands, seasonal river-ice thickness, and stream temperatures. Presently, the YFB is coarsely mapped as spanning the continuous-discontinuous permafrost transition that model analysis shows to be a critical threshold; thus, the YFB may be on the verge of major hydrologic change should the current permafrost extent decrease. This possibility underscores the need for improved characterization of permafrost and other hydrogeologic information in the region via geophysical techniques, remote sensing, and ground-based observations.

Alaska↗
USGS2012

Dissolved organic matter composition of winter flow in the Yukon River basin: Implications of permafrost thaw and increased groundwater discharge

Groundwater discharge to rivers has increased in recent decades across the circumpolar region and has been attributed to thawing permafrost in arctic and subarctic watersheds. Permafrost-driven changes in groundwater discharge will alter the flux of dissolved organic carbon (DOC) in rivers, yet little is known about the chemical composition and reactivity of dissolved organic matter (DOM) of groundwater in permafrost settings. Here, we characterize DOM composition of winter flow in 60 rivers and streams of the Yukon River basin to evaluate the biogeochemical consequences of enhanced groundwater discharge associated with permafrost thaw. DOC concentration of winter flow averaged 3.9 ± 0.5 mg C L −1 , yet was highly variable across basins (ranging from 20 mg C L −1 ). In comparison to the summer-autumn period, DOM composition of winter flow had lower aromaticity (as indicated by specific ultraviolet absorbance at 254 nm, or SUVA 254 ), lower hydrophobic acid content, and a higher proportion of hydrophilic compounds (HPI). Fluorescence spectroscopy and parallel factor analysis indicated enrichment of protein-like fluorophores in some, but not all, winter flow samples. The ratio of DOC to dissolved organic nitrogen, an indicator of DOM biodegradability, was positively correlated with SUVA 254 and negatively correlated with the percentage of protein-like compounds. Using a simple two-pool mixing model, we evaluate possible changes in DOM during the summer-autumn period across a range of conditions reflecting possible increases in groundwater discharge. Across three watersheds, we consistently observed decreases in DOC concentration and SUVA254 and increases in HPI with increasing groundwater discharge. Spatial patterns in DOM composition of winter flow appear to reflect differences in the relative contributions of groundwater from suprapermafrost and subpermafrost aquifers across watersheds. Our findings call for more explicit consideration of DOC loss and stabilization pathways associated with changing subsurface hydrology in watersheds underlain by thawing permafrost.

Global Biogeochemical Cycles↗
USGS2021

Post-drought groundwater storage recovery in California’s Central Valley

Groundwater depletion is a major threat to agricultural and municipal water supply in California's Central Valley. Recent droughts during 2007–2009 and 2012–2016 exacerbated chronic groundwater depletion. However, it is unclear how much groundwater storage recovered from drought-related overdrafts during post-drought years, and how climatic conditions and water management affected recovery times. We estimated groundwater storage change in the Central Valley for April 2002 through September 2019 using four methods: GRACE satellite data, a water balance approach, a hydrologic simulation model, and monitoring wells. We also evaluated the sensitivity of drought recovery to different climate scenarios (recent climate ± droughts and future climate change scenarios: 20 GCMs and 2 RCPs) using water balance method and statistical sampling of historical climate data. Estimated Central Valley groundwater loss during the two droughts ranged from 19 km 3 (2007–2009) to 28 km 3 (2012–2016) (median of four methods). Median aquifer storage recovery was 34% and 19% of the overdraft during the 2010–2011 and 2017–2019 post-drought years, respectively. Numerical experiments show that recovery times are sensitive to climate forcing, with longer recovery times for a future climate scenario that replicate historical climatology relative to historical forcing with no droughts. Overdraft recovery times decrease by ∼2× with implementation of pumping restrictions (30th to 50th percentiles of historical groundwater depletion) to constrain groundwater depletion relative to no restrictions with a no-drought future climatology. This study highlights the importance of considering water management implications for future drought recoveries within the context of climate change scenarios.

California↗
USGS2023

Paired Air and Stream Temperature Analysis (PASTA) to evaluate groundwater influence on streams

Groundwater is critical for maintaining stream baseflow and thermal stability; however, the influence of groundwater on streamflow has been difficult to evaluate at broad spatial scales. Techniques such as baseflow separation necessitate streamflow records and do not directly indicate whether groundwater inflow may be sourced from more dynamic shallow flowpaths. We present a web tool application PASTA (Paired Air and Stream Temperature Analysis; https://cuahsi.shinyapps.io/pasta/ ) that capitalizes on increased public stream temperature data availability and large-scale, gridded climate observations to provide new and efficient insights regarding relative groundwater influence on streams. PASTA analyzes paired air and stream water temperature signals to evaluate spatiotemporal patterns in stream thermal sensitivity and relative groundwater influence, including inference regarding the dominant source groundwater depth (shallow or deep (i.e., approximately >6 m depth)). The tool is linked to publicly available stream temperature datasets and accepts user-uploaded datasets. As local air temperature is not often monitored, PASTA pulls daily air temperature data from the comprehensive Daymet products when directly measured data are unavailable, allowing the repurposing of existing stream temperature data. After data are selected or uploaded, the tool (a) fits sinusoidal curves of daily stream and air temperatures by year (water or calendar) to indicate groundwater influence characteristics and (b) performs linear regressions for stream versus air temperatures to indicate stream thermal sensitivity. Results are exported in ASCII file format, creating an efficient and approachable analysis tool for the adoption of newly developed heat tracing analysis from stream reach to landscape scales.

Water Resources Research↗
USGS2023

Train, inform, borrow, or combine? Approaches to process-guided deep learning for groundwater-influenced stream temperature prediction

Although groundwater discharge is a critical stream temperature control process, it is not explicitly represented in many stream temperature models, an omission that may reduce predictive accuracy, hinder management of aquatic habitat, and decrease user confidence. We assessed the performance of a previously-described process-guided deep learning model of stream temperature in the Delaware River Basin (USA). We found lower accuracy (root mean square error [RMSE] of 1.71 versus 1.35°C) and stronger seasonal bias (absolute mean monthly bias of 1.06 vs. 0.68°C) for reaches primarily influenced by deep groundwater as compared to atmospheric conditions. We then tested four approaches for improving groundwater process representation: (a) a custom loss function leveraging the unique patterns of air and water temperature coupling characteristic of different temperature drivers, (b) inclusion of additional groundwater-relevant catchment attributes, (c) incorporation of additional process model outputs, and (d) a composite model. The custom loss function and the additional attributes significantly improved the predictive accuracy in groundwater-dominated reaches (RMSE of 1.37 and 1.26°C) and reduced the seasonal bias (absolute mean monthly bias of 0.44 and 0.48°C), but neither approach could identify holdout groundwater reaches. Variable importance analysis indicates the custom loss function nudges the model to use the existing inputs more efficiently, whereas with the added features the model relies on a broader suite of inputs. This analysis is a substantial step toward more accurately representing groundwater discharge processes in stream temperature models and will improve predictive accuracy and inform habitat management.

Delaware River Basin↗
USGS2025

Preferential groundwater discharges along stream corridors are disregarded sources of greenhouse gases

Groundwater delivery of greenhouse gases (GHGs) to stream banks and riparian areas, before mixing with surface waters, has not been well quantified. We measured preferential groundwater delivery of GHGs to stream banks within three stream reaches, and found that stream banks with discharging groundwater emitted more CO 2 and were sources of N 2 O compared to stream banks without actively discharging groundwater, which emitted less CO 2 and were N 2 O sinks. At one of our stream reaches, groundwater CO 2 and N 2 O concentrations were 1.4–19.2 and 1.1–40.6 times higher than those in surface water, respectively, and groundwater delivery rates of CO 2 and N 2 O were 1.5 and 1.6 times higher than surface water emissions per unit area. On average, 21% (range 0%–100%) of CO 2 and N 2 O were emitted at the stream bank before mixing with surface waters. Preferential groundwater GHG emissions may contribute substantially to stream corridor emissions and may be underestimated when using a channel-centric approach to estimate riverine GHG budgets.

JGR Biogeosciences↗
USGS2008

Groundwater flow and solute movement to drain laterals, western San Joaquin Valley, California: 1. Geochemical assessment

A study was undertaken to quantitatively evaluate the hydrologic processes affecting the chemical and isotopic composition of drain lateral water in a drained agricultural field in the western San Joaquin Valley, California. The results elucidate the process of mixing of deep and shallow groundwater (below and within 6 m from land surface) entering the drain laterals. The deep groundwater was subject to evapoconcentration prior to drainage system installation and has been displaced downward (to depths greater than 6 m) in the groundwater system. The proportions of deep and shallow groundwater entering the drain laterals was calculated from the end-member oxygen 18 compositions determined in groundwater samples. The percentage of total drain lateral flow which is deep groundwater flow is about 30% for the shallow drain lateral (1.8 m below land surface) (drain lateral 1)) and 60% for the deep drain lateral (2.7 m below land surface (drain lateral 2)). During irrigation, the percentages of deep groundwater flow decrease to 0 and 30% for the shallow and deep drain laterals, respectively. Selenium concentrations in drain lateral waters decrease during irrigation but selenium loads increase. Total estimated annual loads were 1.1 and 5.4 kg of selenium for drain laterals 1 and 2, respectively. Substantial percentages of the annual load occurred during 8 days of irrigation, 23 and 9% for drain laterals 1 and 2, respectively.

California↗
USGS2010

Gravity-driven groundwater flow and slope failure potential: 1. Elastic effective-stress model

Hilly or mountainous topography influences gravity-driven groundwater flow and the consequent distribution of effective stress in shallow subsurface environments. Effective stress, in turn, influences the potential for slope failure. To evaluate these influences, we formulate a two-dimensional, steady state, poroelastic model. The governing equations incorporate groundwater effects as body forces, and they demonstrate that spatially uniform pore pressure changes do not influence effective stresses. We implement the model using two finite element codes. As an illustrative case, we calculate the groundwater flow field, total body force field, and effective stress field in a straight, homogeneous hillslope. The total body force and effective stress fields show that groundwater flow can influence shear stresses as well as effective normal stresses. In most parts of the hillslope, groundwater flow significantly increases the Coulomb failure potential Φ, which we define as the ratio of maximum shear stress to mean effective normal stress. Groundwater flow also shifts the locus of greatest failure potential toward the slope toe. However, the effects of groundwater flow on failure potential are less pronounced than might be anticipated on the basis of a simpler, one-dimensional, limit equilibrium analysis. This is a consequence of continuity, compatibility, and boundary constraints on the two-dimensional flow and stress fields, and it points to important differences between our elastic continuum model and limit equilibrium models commonly used to assess slope stability.

Water Resources Research↗
USGS2010

Chemical evolution of groundwater near a sinkhole lake, northern Florida: 2. Chemical patterns, mass-transfer modeling, and rates of chemical reactions

Chemical patterns along evolutionary groundwater flow paths in silicate and carbonate aquifers were interpreted using solute tracers, carbon and sulfur isotopes, and mass balance reaction modeling for a complex hydrologic system involving groundwater inflow to and outflow from a sinkhole lake in northern Florida. Rates of dominant reactions along defined flow paths were estimated from modeled mass transfer and ages obtained from CFC-modeled recharge dates. Groundwater upgradient from Lake Barco remains oxic as it moves downward, reacting with silicate minerals in a system open to carbon dioxide (CO 2 ), producing only small increases in dissolved species. Beneath and downgradient of Lake Barco the oxic groundwater mixes with lake water leakage in a highly reducing, silicate-carbonate mineral environment. A mixing model, developed for anoxic groundwater downgradient from the lake, accounted for the observed chemical and isotopic composition by combining different proportions of lake water leakage and infiltrating meteoric water. The evolution of major ion chemistry and the 13 C isotopic composition of dissolved carbon species in groundwater downgradient from the lake can be explained by the aerobic oxidation of organic matter in the lake, anaerobic microbial oxidation of organic carbon, and incongruent dissolution of smectite minerals to kaolinite. The dominant process for the generation of methane was by the CO 2 reduction pathway based on the isotopic composition of hydrogen (δ 2 H(CH 4 ) = −186 to −234‰) and carbon (δ 13 C(CH 4 ) = −65.7 to −72.3‰). Rates of microbial metabolism of organic matter, estimated from the mass transfer reaction models, ranged from 0.0047 to 0.039 mmol L −1 yr −1 for groundwater downgradient from the lake.

Water Resources Research↗
USGS2010

Vulnerability of deep groundwater in the Bengal Aquifer System to contamination by arsenic

Shallow groundwater, the primary water source in the Bengal Basin, contains up to 100 times the World Health Organization (WHO) drinking-water guideline of 10g l 1 arsenic (As), threatening the health of 70 million people. Groundwater from a depth greater than 150m, which almost uniformly meets the WHO guideline, has become the preferred alternative source. The vulnerability of deep wells to contamination by As is governed by the geometry of induced groundwater flow paths and the geochemical conditions encountered between the shallow and deep regions of the aquifer. Stratification of flow separates deep groundwater from shallow sources of As in some areas. Oxidized sediments also protect deep groundwater through the ability of ferric oxyhydroxides to adsorb As. Basin-scale groundwater flow modelling suggests that, over large regions, deep hand-pumped wells for domestic supply may be secure against As invasion for hundreds of years. By contrast, widespread deep irrigation pumping might effectively eliminate deep groundwater as an As-free resource within decades. Finer-scale models, incorporating spatial heterogeneity, are needed to investigate the security of deep municipal abstraction at specific urban locations. ?? 2010 Macmillan Publishers Limited. All rights reserved.

Nature Geoscience↗
USGS2010

Optimal pump and recharge management model for nitrate removal in the Warren groundwater basin, California

The town of Yucca Valley located in the southwest part of the Mojave Desert in southern California relies on groundwater pumping from the Warren groundwater basin as its sole source of water supply. This significant dependency has resulted in a large imbalance between groundwater pumpage and natural recharge, causing groundwater levels in the basin to decline more than 90 m from the late 1940s to 1994. Consequently, an artificial recharge program proposed by the Hi-Desert Water District, which provides water service to the town of Yucca Valley, was implemented for the purpose of recovering the groundwater levels; however, the rise in groundwater levels has caused nitrate ( NO 3 ) "> ( NO 3 ) (NO3) concentration to increase simultaneously. The purpose of this study is to develop an optimal pump and recharge strategy for a planned conjunctive-use project to remove the high-nitrate concentration while maintaining groundwater levels at desired elevations at specified locations as well as meeting water demand. An optimization/management model is formulated with a linear objective function and nonlinear constraints. The response matrix approach is used to link the optimization model with the simulation model. Because of nonlinearity, the response matrix is updated and iteration is required for convergence. A systematic scheme is also developed for finding a feasible initial policy. Three different scenarios are considered in the management model. The results obtained from each scenario are analyzed and discussed.

California↗
USGS2012

Groundwater depletion and sustainability of irrigation in the US High Plains and Central Valley

Aquifer overexploitation could significantly impact crop production in the United States because 60% of irrigation relies on groundwater. Groundwater depletion in the irrigated High Plains and California Central Valley accounts for ∼50% of groundwater depletion in the United States since 1900. A newly developed High Plains recharge map shows that high recharge in the northern High Plains results in sustainable pumpage, whereas lower recharge in the central and southern High Plains has resulted in focused depletion of 330 km 3 of fossil groundwater, mostly recharged during the past 13,000 y. Depletion is highly localized with about a third of depletion occurring in 4% of the High Plains land area. Extrapolation of the current depletion rate suggests that 35% of the southern High Plains will be unable to support irrigation within the next 30 y. Reducing irrigation withdrawals could extend the lifespan of the aquifer but would not result in sustainable management of this fossil groundwater. The Central Valley is a more dynamic, engineered system, with north/south diversions of surface water since the 1950s contributing to ∼7× higher recharge. However, these diversions are regulated because of impacts on endangered species. A newly developed Central Valley Hydrologic Model shows that groundwater depletion since the 1960s, totaling 80 km 3 , occurs mostly in the south (Tulare Basin) and primarily during droughts. Increasing water storage through artificial recharge of excess surface water in aquifers by up to 3 km 3 shows promise for coping with droughts and improving sustainability of groundwater resources in the Central Valley.

California, Colorago, Kansas, Nebraska, New Mexico↗
USGS2017

Groundwater declines are linked to changes in Great Plains stream fish assemblages

Groundwater pumping for agriculture is a major driver causing declines of global freshwater ecosystems, yet the ecological consequences for stream fish assemblages are rarely quantified. We combined retrospective (1950–2010) and prospective (2011–2060) modeling approaches within a multiscale framework to predict change in Great Plains stream fish assemblages associated with groundwater pumping from the United States High Plains Aquifer. We modeled the relationship between the length of stream receiving water from the High Plains Aquifer and the occurrence of fishes characteristic of small and large streams in the western Great Plains at a regional scale and for six subwatersheds nested within the region. Water development at the regional scale was associated with construction of 154 barriers that fragment stream habitats, increased depth to groundwater and loss of 558 km of stream, and transformation of fish assemblage structure from dominance by large-stream to small-stream fishes. Scaling down to subwatersheds revealed consistent transformations in fish assemblage structure among western subwatersheds with increasing depths to groundwater. Although transformations occurred in the absence of barriers, barriers along mainstem rivers isolate depauperate western fish assemblages from relatively intact eastern fish assemblages. Projections to 2060 indicate loss of an additional 286 km of stream across the region, as well as continued replacement of large-stream fishes by small-stream fishes where groundwater pumping has increased depth to groundwater. Our work illustrates the shrinking of streams and homogenization of Great Plains stream fish assemblages related to groundwater pumping, and we predict similar transformations worldwide where local and regional aquifer depletions occur.

Colorado, Kansas, Nebraska↗
USGS2022

Sea-level rise and warming mediate coastal groundwater discharge in the Arctic

Groundwater discharge is an important mechanism through which fresh water and associated solutes are delivered to the ocean. Permafrost environments have traditionally been considered hydrogeologically inactive, yet with accelerated climate change and permafrost thaw, groundwater flow paths are activating and opening subsurface connections to the coastal zone. While warming has the potential to increase land-sea connectivity, sea-level change has the potential to alter land-sea hydraulic gradients and enhance coastal permafrost thaw, resulting in a complex interplay that will govern future groundwater discharge dynamics along Arctic coastlines. Here, we use a recently developed permafrost hydrological model that simulates variable-density groundwater flow and salinity-dependent freeze-thaw to investigate the impacts of sea-level change and land and ocean warming on the magnitude, spatial distribution, and salinity of coastal groundwater discharge. Results project both an increase and decrease in discharge with climate change depending on the rate of warming and sea-level change. Under high warming and low sea-level rise scenarios, results show up to a 58% increase in coastal groundwater discharge by 2100 due to the formation of a supra-permafrost aquifer that enhances freshwater delivery to the coastal zone. With higher rates of sea-level rise, the increase in discharge due to warming is reduced to 21% as sea-level rise decreased land-sea hydraulic gradients. Under lower warming scenarios for which supra-permafrost groundwater flow was not established, discharge decreased by up to 26% between 1980 and 2100 for high sea-level rise scenarios and increased only 8% under low sea-level rise scenarios. Thus, regions with higher warming rates and lower rates of sea-level change (e.g. northern Nunavut, Canada) will experience a greater increase in discharge than regions with lower warming rates and higher rates of sea-level change. The magnitude, location and salinity of discharge have important implications for ecosystem function, water quality, and carbon dynamics in coastal zones.

Environmental Research Letters↗
USGS2017

Forecasting the probability of future groundwater levels declining below specified low thresholds in the conterminous U.S.

We present a logistic regression approach for forecasting the probability of future groundwater levels declining or maintaining below specific groundwater-level thresholds. We tested our approach on 102 groundwater wells in different climatic regions and aquifers of the United States that are part of the U.S. Geological Survey Groundwater Climate Response Network. We evaluated the importance of current groundwater levels, precipitation, streamflow, seasonal variability, Palmer Drought Severity Index, and atmosphere/ocean indices for developing the logistic regression equations. Several diagnostics of model fit were used to evaluate the regression equations, including testing of autocorrelation of residuals, goodness-of-fit metrics, and bootstrap validation testing. The probabilistic predictions were most successful at wells with high persistence (low month-to-month variability) in their groundwater records and at wells where the groundwater level remained below the defined low threshold for sustained periods (generally three months or longer). The model fit was weakest at wells with strong seasonal variability in levels and with shorter duration low-threshold events. We identified challenges in deriving probabilistic-forecasting models and possible approaches for addressing those challenges.

Journal of the American Water Resources Associatio↗
USGS2012

Glaciation and regional groundwater flow in the Fennoscandian shield

Regional-scale groundwater flow modeling of the Fennoscandian shield suggests that groundwater flow can be strongly affected by future climate change and glaciation. We considered variable-density groundwater flow in a 1500-km-long and approximately 10-km-deep cross-section through southern Sweden. Groundwater flow and shield brine transport in the cross-sectional model were analyzed under projected surface conditions for the next 140 ka. Simulations suggest that blockage of recharge and discharge by low-permeability permafrost or cold-based ice causes sinking of brine and consequent freshening of near-surface water in areas of natural discharge. Although recharge of basal meltwater is limited by the requirement that water pressure at the base of the ice sheet not exceed the pressure exerted by the weight of the ice, warm-based ice with basal melting creates a potential for groundwater recharge rates much larger than those of present, ice-free conditions. In the simulations, regional-scale redistribution of recharged water by subsurface flow is minor over the duration of a glacial advance (approximately 10 ka). During glacial retreat, significant upward flow of groundwater may occur below the ice sheet owing to pressure release. If the mechanical loading efficiency of the rocks is high, both subsurface penetration of meltwater during glacial advance and up-flow during glacial retreat are reduced because of loading-induced pressure changes. The maximum rate of groundwater discharge in the simulations occurs at the receding ice margin, and some discharge occurs below incursive postglacial seas. Recharge of basal meltwater could decrease the concentration of dissolved solids significantly below present-day levels at depths of up to several kilometers and may bring oxygenated conditions to an otherwise reducing chemical environment for periods exceeding 10 ka.

Fennoscandian Shield↗