skip to main content

Title: Uncertainty of simulated groundwater recharge at different global warming levels: a global-scale multi-model ensemble study
Abstract. Billions of people rely on groundwater as being an accessible source of drinking water and for irrigation, especially in times of drought. Its importance will likely increase with a changing climate. It is still unclear, however, how climate change will impact groundwater systems globally and, thus, the availability of this vital resource. Groundwater recharge is an important indicator for groundwater availability, but it is a water flux that is difficult to estimate as uncertainties in the water balance accumulate, leading to possibly large errors in particular in dry regions. This study investigates uncertainties in groundwater recharge projections using a multi-model ensemble of eight global hydrological models (GHMs) that are driven by the bias-adjusted output of four global circulation models (GCMs). Pre-industrial and current groundwater recharge values are compared with recharge for different global warming (GW) levels as a result of three representative concentration pathways (RCPs). Results suggest that projected changes strongly vary among the different GHM–GCM combinations, and statistically significant changes are only computed for a few regions of the world. Statistically significant GWR increases are projected for northern Europe and some parts of the Arctic, East Africa, and India. Statistically significant decreases are simulated in southern Chile, parts of Brazil, central USA, more » the Mediterranean, and southeastern China. In some regions, reversals of groundwater recharge trends can be observed with global warming. Because most GHMs do not simulate the impact of changing atmospheric CO2 and climate on vegetation and, thus, evapotranspiration, we investigate how estimated changes in GWR are affected by the inclusion of these processes. In some regions, inclusion leads to differences in groundwater recharge changes of up to 100 mm per year. Most GHMs with active vegetation simulate less severe decreases in groundwater recharge than GHMs without active vegetation and, in some regions, even increases instead of decreases are simulated. However, in regions where GCMs predict decreases in precipitation and where groundwater availability is the most important, model agreement among GHMs with active vegetation is the lowest. Overall, large uncertainties in the model outcomes suggest that additional research on simulating groundwater processes in GHMs is necessary. « less
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Hydrology and Earth System Sciences
Page Range or eLocation-ID:
787 to 810
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Assessing impacts of climate change on hydrologic systemsis critical for developing adaptation and mitigation strategies for waterresource management, risk control, and ecosystem conservation practices. Suchassessments are commonly accomplished using outputs from a hydrologic modelforced with future precipitation and temperature projections. The algorithmsused for the hydrologic model components (e.g., runoff generation) canintroduce significant uncertainties into the simulated hydrologic variables.Here, a modeling framework was developed that integrates multiple runoffgeneration algorithms with a routing model and associated parameteroptimizations. This framework is able to identify uncertainties from bothhydrologic model components and climate forcings as well as associatedparameterization. Three fundamentally different runoff generationapproaches, runoff coefficient method (RCM, conceptual), variableinfiltration capacity (VIC, physically based, infiltration excess), andsimple-TOPMODEL (STP, physically based, saturation excess), were coupledwith the Hillslope River Routing model to simulate surface/subsurface runoffand streamflow. A case study conducted in Santa Barbara County, California,reveals increased surface runoff in February and March but decreasedrunoff in other months, a delayed (3 d, median) and shortened (6 d,median) wet season, and increased daily discharge especially for theextremes (e.g., 100-year flood discharge, Q100). The Bayesian modelaveraging analysis indicates that the probability of such an increase can be up to85 %. For projected changes in runoff and discharge, general circulationmodels (GCMs) and emission scenariosmore »are two major uncertainty sources,accounting for about half of the total uncertainty. For the changes inseasonality, GCMs and hydrologic models are two major uncertaintycontributors (∼35 %). In contrast, the contribution ofhydrologic model parameters to the total uncertainty of changes in thesehydrologic variables is relatively small (<6 %), limiting theimpacts of hydrologic model parameter equifinality in climate change impactanalysis. This study provides useful information for practices associatedwith water resources, risk control, and ecosystem conservation and forstudies related to hydrologic model evaluation and climate change impactanalysis for the study region as well as other Mediterranean regions.« less
  2. Abstract

    Global groundwater resources are stressed and the effects of climate change are projected to further disrupt recharge processes. Therefore, we must identify the buffers to climate change in hydrogeologic systems in order to understand which groundwater resources will be disproportionally affected by these changes. Here, we utilize a novel combination of remote sensing (e.g. Landsat) and groundwater residence time data (3H,36Cl) to identify the factors controlling the hydrogeologic stability of aridland mountain-front springs in response to a major climate event, the 2011–2017 California drought. Desert springs within Owens Valley (CA) support unique ecosystems that are surrounded by lush, green vegetation sustained only by discharging groundwater and are not reliant on localized precipitation. Therefore, the health or ecological response of this vegetation is a direct reflection of the hydrogeologic stability of the mountain-block groundwater system since water is the limiting resource for riparian plant growth in arid regions. We compared spring water residence times to vegetation health metrics computed from Landsat imagery leading up to and during the drought interval. We observe that the vegetation surrounding springs discharging a high fraction of modern and bomb-pulse groundwater (<100 years) showed evidence of increased drying and desiccation as the drought progressed. Inmore »comparison, springs discharging a higher fraction of old groundwater (>100 years) showed little response thereby supporting the conceptual model where old groundwater, i.e. a distribution of deep and stable groundwater flowpaths, buffers short- to long-term climate perturbations and may provide hydrogeologic resistance to future effects from climate change.

    « less
  3. This study presents near future (2020–2044) temperature and precipitation changes over the Antarctic Peninsula under the high-emission scenario (RCP8.5). We make use of historical and projected simulations from 19 global climate models (GCMs) participating in Coupled Model Intercomparison Project phase 5 (CMIP5). We compare and contrast GCMs projections with two groups of regional climate model simulations (RCMs): (1) high resolution (15-km) simulations performed with Polar-WRF model forced with bias-corrected NCAR-CESM1 (NC-CORR) over the Antarctic Peninsula, (2) medium resolution (50-km) simulations of KNMI-RACMO21P forced with EC-EARTH (EC) obtained from the CORDEX-Antarctica. A further comparison of historical simulations (1981–2005) with respect to ERA5 reanalysis is also included for circulation patterns and near-surface temperature climatology. In general, both RCM boundary conditions represent well the main circulation patterns of the historical period. Nonetheless, there are important differences in projections such as a notable deepening and weakening of the Amundsen Sea Low in EC and NC-CORR, respectively. Mean annual near-surface temperatures are projected to increase by about 0.5–1.5 ◦ C across the entire peninsula. Temperature increase is more substantial in autumn and winter ( ∼ 2 ◦C). Following opposite circulation pattern changes, both EC and NC-CORR exhibit different warming rates, indicating a possible continuation ofmore »natural decadal variability. Although generally showing similar temperature changes, RCM projections show less warming and a smaller increase in melt days in the Larsen Ice Shelf compared to their respective driving fields. Regarding precipitation, there is a broad agreement among the simulations, indicating an increase in mean annual precipitation ( ∼ 5 to 10%). However, RCMs show some notable differences over the Larsen Ice Shelf where total precipitation decreases (for RACMO) and shows a small increase in rain frequency. We conclude that it seems still difficult to get consistent projections from GCMs for the Antarctic Peninsula as depicted in both RCM boundary conditions. In addition, dominant and common changes from the boundary conditions are largely evident in the RCM simulations. We argue that added value of RCM projections is driven by processes shaped by finer local details and different physics schemes that are introduced by RCMs, particularly over the Larsen Ice Shelf.« less
  4. Reliable projections of future changes in tropical cyclone (TC) characteristics are highly dependent on the ability of global climate models (GCMs) to simulate the observed characteristics of TCs (i.e., their frequency, genesis locations, movement, and intensity). Here, we investigate the performance of a suite of GCMs from the U.S. CLIVAR Working Group on Hurricanes in simulating observed climatological features of TCs in the Southern Hemisphere. A subset of these GCMs is also explored under three idealized warming scenarios. Two types of simulated TC tracks are evaluated on the basis of a commonly applied cluster analysis: 1) explicitly simulated tracks, and 2) downscaled tracks, derived from a statistical–dynamical technique that depends on the models’ large-scale environmental fields. Climatological TC properties such as genesis locations, annual frequency, lifetime maximum intensity (LMI), and seasonality are evaluated for both track types. Future changes to annual frequency, LMI, and the latitude of LMI are evaluated using the downscaled tracks where large sample sizes allow for statistically robust results. An ensemble approach is used to assess future changes of explicit tracks owing to their small number of realizations. We show that the downscaled tracks generally outperform the explicit tracks in relation to many of the climatologicalmore »features of Southern Hemisphere TCs, despite a few notable biases. Future changes to the frequency and intensity of TCs in the downscaled simulations are found to be highly dependent on the warming scenario and model, with the most robust result being an increase in the LMI under a uniform 2°C surface warming.

    « less
  5. Abstract

    Simultaneous heatwaves affecting multiple regions (referred to as concurrent heatwaves) pose compounding threats to various natural and societal systems, including global food chains, emergency response systems, and reinsurance industries. While anthropogenic climate change is increasing heatwave risks across most regions, the interactions between warming and circulation changes that yield concurrent heatwaves remain understudied. Here, we quantify historical (1979–2019) trends in concurrent heatwaves during the warm season [May–September (MJJAS)] across the Northern Hemisphere mid- to high latitudes. We find a significant increase of ∼46% in the mean spatial extent of concurrent heatwaves and ∼17% increase in their maximum intensity, and an approximately sixfold increase in their frequency. Using self-organizing maps, we identify large-scale circulation patterns (300 hPa) associated with specific concurrent heatwave configurations across Northern Hemisphere regions. We show that observed changes in the frequency of specific circulation patterns preferentially increase the risk of concurrent heatwaves across particular regions. Patterns linking concurrent heatwaves across eastern North America, eastern and northern Europe, parts of Asia, and the Barents and Kara Seas show the largest increases in frequency (∼5.9 additional days per decade). We also quantify the relative contributions of circulation pattern changes and warming to overall observed concurrent heatwave day frequencymore »trends. While warming has a predominant and positive influence on increasing concurrent heatwave frequency, circulation pattern changes have a varying influence and account for up to 0.8 additional concurrent heatwave days per decade. Identifying regions with an elevated risk of concurrent heatwaves and understanding their drivers is indispensable for evaluating projected climate risks on interconnected societal systems and fostering regional preparedness in a changing climate.

    Significance Statement

    Heatwaves pose a major threat to human health, ecosystems, and human systems. Simultaneous heatwaves affecting multiple regions can exacerbate such threats. For example, multiple food-producing regions simultaneously undergoing heat-related crop damage could drive global food shortages. We assess recent changes in the occurrence of simultaneous large heatwaves. Such simultaneous heatwaves are 7 times more likely now than 40 years ago. They are also hotter and affect a larger area. Their increasing occurrence is mainly driven by warming baseline temperatures due to global heating, but changes in weather patterns contribute to disproportionate increases over parts of Europe, the eastern United States, and Asia. Better understanding the drivers of weather pattern changes is therefore important for understanding future concurrent heatwave characteristics and their impacts.

    « less