Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Climate sensitivity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

The role of sea ice in 2 x CO2 climate model sensitivity. Part 1: The total influence of sea ice thickness and extent

As a first step in investigating the effects of sea ice changes on the climate sensitivity to doubled atmospheric CO2, the authors use a standard simple sea ice model while varying the sea ice distributions and thicknesses in the control run. Thinner ice amplifies the atmospheric temperature senstivity in these experiments by about 15% (to a warming of 4.8 C), because it is easier for the thinner ice to be removed as the climate warms. Thus, its impact on sensitivity is similar to that of greater sea ice extent in the control run, which provides more opportunity for sea ice reduction. An experiment with sea ice not allowed to change between the control and doubled CO2 simulations illustrates that the total effect of sea ice on surface air temperature changes, including cloud cover and water vapor feedbacks that arise in response to sea ice variations, amounts to 37% of the temperature sensitivity to the CO2 doubling, accounting for 1.56 C of the 4.17 C global warming. This is about four times larger than the sea ice impact when no feedbacks are allowed. The different experiments produce a range of results for southern high latitudes with the hydrologic budget over Antarctica implying sea level increases of varying magnitude or no change. These results highlight the importance of properly constraining the sea ice response to climate perturbations, necessitating the use of more realistic sea ice and ocean models.

Rind, D.↗
NASA NTRS1999

Sensitivity of Climate to Changes in NDVI

The sensitivity of global and regional climate to changes in vegetation density is investigated using a coupled biosphere-atmosphere model. The magnitude of the vegetation changes and their spatial distribution are based on natural decadal variability of the normalized difference vegetation index (ndvi). Different scenarios using maximum and minimum vegetation cover were derived from satellite records spanning the period 1982-1990. Albedo decreased in the northern latitudes and increased in the tropics with increased ndvi. The increase in vegetation density revealed that the vegetation's physiological response was constrained by the limits of the available water resources. The difference between the maximum and minimum vegetation scenarios resulted in a 46% increase in absorbed visible solar radiation and a similar increase in gross photosynthetic C02 uptake on a global annual basis. This caused the canopy transpiration and interception fluxes to increase, and reduced those from the soil. The redistribution of the surface energy fluxes substantially reduced the Bowen ratio during the growing season, resulting in cooler and moister near-surface climate, except when soil moisture was limiting. Important effects of increased vegetation on climate are : (1) A cooling of about 1.8 K in the northern latitudes during the growing season and a slight warming during the winter, which is primarily due to the masking of high albedo of snow by a denser canopy. and (2) A year round cooling of 0.8 K in the tropics. These results suggest that increases in vegetation density could partially compensate for parallel increases in greenhouse warming . Increasing vegetation density globally caused both evapotranspiration and precipitation to increase. Evapotranspiration, however increased more than precipitation resulting in a global soil-water deficit of about 15 %. A spectral analysis on the simulated results showed that changes in the state of vegetation could affect the low-frequency modes of the precipitation distribution and might reduce its low frequency variability in the tropics while increasing it in northern latitudes.

Bounoua, L.↗
NASA NTRS2010

Why Hasn't Earth Warmed as Much as Expected?

The observed increase in global mean surface temperature (GMST) over the industrial era is less than 40% of that expected from observed increases in long-lived greenhouse gases together with the best-estimate equilibrium climate sensitivity given by the 2007 Assessment Report of the Intergovernmental Panel on Climate Change (IPCC). Possible reasons for this warming discrepancy are systematically examined here. The warming discrepancy is found to be due mainly to some combination of two factors: the IPCC best estimate of climate sensitivity being too high and/or the greenhouse gas forcing being partially offset by forcing by increased concentrations of atmospheric aerosols; the increase in global heat content due to thermal disequilibrium accounts for less than 25% of the discrepancy, and cooling by natural temperature variation can account for only about 15 %. Current uncertainty in climate sensitivity is shown to preclude determining the amount of future fossil fuel CO2 emissions that would be compatible with any chosen maximum allowable increase in GMST; even the sign of such allowable future emissions is unconstrained. Resolving this situation, by empirical determination of the earth's climate sensitivity from the historical record over the industrial period or through use of climate models whose accuracy is evaluated by their performance over this period, is shown to require substantial reduction in the uncertainty of aerosol forcing over this period.

Schwartz, Stephen E.↗
NASA NTRS2024

Observational Constraint on A Feedback From Supercooled Clouds Reduces Projected Warming Uncertainty

The increase of carbon-dioxide-doubling-induced warming (climate sensitivity) in the latest climate models is primarily attributed to a larger extratropical cloud feedback. This is thought to be partly driven by a greater ratio of supercooled liquid-phase clouds to all clouds, termed liquid phase ratio. We use an instrument simulator approach to show that this ratio has increased in the latest climate models and is overestimated rather than underestimated as previously thought. In our analysis of multiple models, a greater ratio corresponds to stronger negative cloud feedback, in contradiction with single-model-based studies. We trace this unexpected result to a cloud feedback involving a shift from supercooled to warm clouds as climate warms, which corresponds to greater cloud amount and optical depth and weakens the extratropical cloud feedback. Better constraining this ratio in climate models – and thus this supercooled cloud feedback – impacts their climate sensitivities by up to 1 ˚C and reduces inter-model spread.

Climate and Earth system modelling↗
NASA NTRS2015

Modelling Bambara Groundnut Yield in Southern Africa: Towards a Climate-Resilient Future

Current agriculture depends on a few major species grown as monocultures that are supported by global research underpinning current productivity. However, many hundreds of alternative crops have the potential to meet real world challenges by sustaining humanity, diversifying agricultural systems for food and nutritional security, and especially responding to climate change through their resilience to certain climate conditions. Bambara groundnut (Vigna subterranea (L.) Verdc.), an underutilised African legume, is an exemplar crop for climate resilience. Predicted yield performances of Bambara groundnut by AquaCrop (a crop−water productivity model) were evaluated for baseline (1980−2009) and mid-century climates (2040−2069) under 20 downscaled Global Climate Models (CMIP5-RCP8.5), as well as for climate sensitivities (AgMIPC3MP) across 3 locations in Southern Africa (Botswana, South Africa, Namibia). Different land - races of Bambara groundnut originating from various semi-arid African locations showed diverse yield performances with diverse sensitivities to climate. S19 originating from hot-dry conditions in Namibia has greater future yield potential compared to the Swaziland landrace Uniswa Red-UN across study sites. South Africa has the lowest yield under the current climate, indicating positive future yield trends. Namibia reported the highest baseline yield at optimum current temperatures, indicating less yield potential in future climates. Bambara groundnut shows positive yield potential at temperatures of up to ~31degC, with further warming pushing yields down. Thus, many regions in Southern Africa can utilize Bambara groundnut successfully in the coming decades. This modelling exercise supports decisions on genotypic suitability for present and future climates at specific locations.

farm crops↗
NASA NTRS2020

Clouds and Convective Self-Aggregation in a Multi-Model Ensemble of Radiative-Convective Equilibrium Simulations

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective56equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climatecsensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique amongst intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than un-aggregated simulations. Plain Language Summary This study investigates tropical clouds and climate using results from more than 30 different numerical models set up in a simplified framework. The dataset of model simulations is unique in that it includes a wide range of model types configured in a consistent manner. We address some of the biggest open questions in climate science, including how cloud properties change with warming and the role that the tendency of clouds to form clusters plays in determining the average climate and how climate changes. While there are large differences in how the different models simulate average temperature, humidity, and cloudiness, in a majority of models, the amount of high clouds decreases as climate warms. Nearly all models simulate a tendency for clouds to cluster together. There is agreement that when the clouds are clustered, the atmosphere is drier with fewer clouds overall. We don’t find a conclusive result for how cloud clustering changes as the climate warms.

Allison A. Wing↗
NASA NTRS1997

Improved Upper Ocean/Sea Ice Modeling in the GISS GCM for Investigating Climate Change

This project built on our previous results in which we highlighted the importance of sea ice in overall climate sensitivity by determining that for both warming and cooling climates, when sea ice was not allowed to change, climate sensitivity was reduced by 35-40%. We also modified the Goddard Institute for Space Studies (GISS) 8 deg x lO deg atmospheric General Circulation Model (GCM) to include an upper-ocean/sea-ice model involving the Semtner three-layer ice/snow thermodynamic model, the Price et al. (1986) ocean mixed layer model and a general upper ocean vertical advection/diffusion scheme for maintaining and fluxing properties across the pycnocline. This effort, in addition to improving the sea ice representation in the AGCM, revealed a number of sensitive components of the sea ice/ocean system. For example, the ability to flux heat through the ice/snow properly is critical in order to resolve the surface temperature properly, since small errors in this lead to unrestrained climate drift. The present project, summarized in this report, had as its objectives: (1) introducing a series of sea ice and ocean improvements aimed at overcoming remaining weaknesses in the GCM sea ice/ocean representation, and (2) performing a series of sensitivity experiments designed to evaluate the climate sensitivity of the revised model to both Antarctic and Arctic sea ice, determine the sensitivity of the climate response to initial ice distribution, and investigate the transient response to doubling CO2.

Source record↗
NASA NTRS1998

Improved Upper Ocean/Sea Ice Modeling in the GISS GCM for Investigating Climate Change

This project built on our previous results in which we highlighted the importance of sea ice in overall climate sensitivity by determining that for both warming and cooling climates, when sea ice was not allowed to change, climate sensitivity was reduced by 35-40%. We also modified the GISS 8 deg x lO deg atmospheric GCM to include an upper-ocean/sea-ice model involving the Semtner three-layer ice/snow thermodynamic model, the Price et al. (1986) ocean mixed layer model and a general upper ocean vertical advection/diffusion scheme for maintaining and fluxing properties across the pycnocline. This effort, in addition to improving the sea ice representation in the AGCM, revealed a number of sensitive components of the sea ice/ocean system. For example, the ability to flux heat through the ice/snow properly is critical in order to resolve the surface temperature properly, since small errors in this lead to unrestrained climate drift. The present project, summarized in this report, had as its objectives: (1) introducing a series of sea ice and ocean improvements aimed at overcoming remaining weaknesses in the GCM sea ice/ocean representation, and (2) performing a series of sensitivity experiments designed to evaluate the climate sensitivity of the revised model to both Antarctic and Arctic sea ice, determine the sensitivity of the climate response to initial ice distribution, and investigate the transient response to doubling CO2.

Source record↗
NASA NTRS2019

Efficacy of Climate Forcings in PDRMIP Models

Quantifying the efficacy of different climate forcings is important for understanding the real‐world climate sensitivity. This study presents a systematic multimodel analysis of different climate driver efficacies using simulations from the Precipitation Driver and Response Model Intercomparison Project (PDRMIP). Efficacies calculated from instantaneous radiative forcing deviate considerably from unity across forcing agents and models. Effective radiative forcing (ERF) is a better predictor of global mean near‐surface air temperature (GSAT) change. Efficacies are closest to one when ERF is computed using fixed sea surface temperature experiments and adjusted for land surface temperature changes using radiative kernels. Multimodel mean efficacies based on ERF are close to one for global perturbations of methane, sulfate, black carbon, and insolation, but there is notable intermodel spread. We do not find robust evidence that the geographic location of sulfate aerosol affects its efficacy. GSAT is found to respond more slowly to aerosol forcing than CO2 in the early stages of simulations. Despite these differences, we find that there is no evidence for an efficacy effect on historical GSAT trend estimates based on simulations with an impulse response model, nor on the resulting estimates of climate sensitivity derived from the historical period. However, the considerable intermodel spread in the computed efficacies means that we cannot rule out an efficacy‐induced bias of ±0.4 K in equilibrium climate sensitivity to CO2 doubling when estimated using the historical GSAT trend.

T B Richardson↗
NASA NTRS2020

Effective radiative forcing and adjustments in CMIP6 models

The effective radiative forcing, which includes the instantaneous forcing plus adjustments from the atmosphere and surface, has emerged as the key metric of evaluating human and natural influence on the climate. We evaluate effective radiative forcing and adjustments in 17 contemporary climate models that are participating in the Coupled Model Intercomparison Project (CMIP6) and have contributed to the Radiative Forcing Model Intercomparison Project (RFMIP). Present-day (2014) global-mean anthropogenic forcing relative to pre-industrial (1850) levels from climate models stands at 2.00 (±0.23) W/sq. m, comprised of 1.81 (±0.09) W/sq. m from CO2, 1.08 (± 0.21) W/sq. m from other well-mixed greenhouse gases, −1.01 (± 0.23) W/sq. m from aerosols and −0.09 (±0.13) W/sq. m from land use change. Quoted uncertainties are 1 standard deviation across model best estimates, and 90 % confidence in the reported forcings, due to internal variability, is typically within 0.1 W/sq. m. The majority of the remaining 0.21 W/sq. m is likely to be from ozone. In most cases, the largest contributors to the spread in effective radiative forcing (ERF) is from the instantaneous radiative forcing (IRF) and from cloud responses, particularly aerosol–cloud interactions to aerosol forcing. As determined in previous studies, cancellation of tropospheric and surface adjustments means that the stratospherically adjusted radiative forcing is approximately equal to ERF for greenhouse gas forcing but not for aerosols, and consequentially, not for the anthropogenic total. The spread of aerosol forcing ranges from −0.63 to −1.37 W/sq. m, exhibiting a less negative mean and narrower range compared to 10 CMIP5 models. The spread in 4×CO2 forcing has also narrowed in CMIP6 compared to 13 CMIP5 models. Aerosol forcing is uncorrelated with climate sensitivity. Therefore, there is no evidence to suggest that the increasing spread in climate sensitivity in CMIP6 models, particularly related to high-sensitivity models, is a consequence of a stronger negative present-day aerosol forcing and little evidence that modelling groups are systematically tuning climate sensitivity or aerosol forcing to recreate observed historical warming.

effective radiative forcing↗
NASA NTRS2012

Reply to ''Comments on 'Why Hasn't Earth Warmed as much as Expected?'''

In response to our article, Why Hasnt Earth Warmed as Much as Expected? (2010), Knutti and Plattner (2012) wrote a rebuttal. The term climate sensitivity is usually defined as the change in global mean surface temperature that is produced by a specified change in forcing, such as a change in solar heating or greenhouse gas concentrations. We had argued in the 2010 paper that although climate models can reproduce the global mean surface temperature history over the past century, the uncertainties in these models, due primarily to the uncertainty in climate forcing by airborne particles, mean that the models lack the confidence to actually constrain the climate sensitivity within useful limits for climate prediction. Knutti and Plattner are climate modelers, and they argued essentially that because the models could reproduce the surface temperature history, the issue we raised was moot. Our response amounts to straightening out this confusion; for the models to be constraining, they must be able to reproduce the surface temperature history with sufficient confidence, not just to match the measurements, but to exclude alternative histories. As before, we concluded that if we can actually make the aerosol measurements using currently available, state-of-the-art techniques, we can determine the aerosol climate forcing to the degree required to constrain that aspect of model climate sensitivity. A technical issue relating to the timescale over which a change in CO2 emissions would be equilibrated in the environmental energy balance was also discussed, again, a matter of differences in terminology.

Schwartz, Stephen E.↗
NASA NTRS2012

Making Sense of Palaeoclimate Sensitivity

Many palaeoclimate studies have quantified pre-anthropogenic climate change to calculate climate sensitivity (equilibrium temperature change in response to radiative forcing change), but a lack of consistent methodologies produces a wide range of estimates and hinders comparability of results. Here we present a stricter approach, to improve intercomparison of palaeoclimate sensitivity estimates in a manner compatible with equilibrium projections for future climate change. Over the past 65 million years, this reveals a climate sensitivity (in K W−1 m2) of 0.3-1.9 or 0.6-1.3 at 95% or 68% probability, respectively. The latter implies a warming of 2.2-4.8 K per doubling of atmospheric CO2, which agrees with IPCC estimates.

Rohling, E. J.↗
NASA NTRS2013

Fast Atmosphere-Ocean Model Runs with Large Changes in CO2

How does climate sensitivity vary with the magnitude of climate forcing? This question was investigated with the use of a modified coupled atmosphere-ocean model, whose stability was improved so that the model would accommodate large radiative forcings yet be fast enough to reach rapid equilibrium. Experiments were performed in which atmospheric CO2 was multiplied by powers of 2, from 1/64 to 256 times the 1950 value. From 8 to 32 times, the 1950 CO2, climate sensitivity for doubling CO2 reaches 8 C due to increases in water vapor absorption and cloud top height and to reductions in low level cloud cover. As CO2 amount increases further, sensitivity drops as cloud cover and planetary albedo stabilize. No water vapor-induced runaway greenhouse caused by increased CO2 was found for the range of CO2 examined. With CO2 at or below 1/8 of the 1950 value, runaway sea ice does occur as the planet cascades to a snowball Earth climate with fully ice covered oceans and global mean surface temperatures near 30 C.

ocean models↗
NASA NTRS2010

Contributions to Future Stratospheric Climate Change: An Idealized Chemistry-Climate Model Sensitivity Study

Within the framework of an idealized model sensitivity study, three of the main contributors to future stratospheric climate change are evaluated: increases in greenhouse gas concentrations, ozone recovery, and changing sea surface temperatures (SSTs). These three contributors are explored in combination and separately, to test the interactions between ozone and climate; the linearity of their contributions to stratospheric climate change is also assessed. In a simplified chemistry-climate model, stratospheric global mean temperature is most sensitive to CO2 doubling, followed by ozone depletion, then by increased SSTs. At polar latitudes, the Northern Hemisphere (NH) stratosphere is more sensitive to changes in CO2, SSTs and O3 than is the Southern Hemisphere (SH); the opposing responses to ozone depletion under low or high background CO2 concentrations, as seen with present-day SSTs, are much weaker and are not statistically significant under enhanced SSTs. Consistent with previous studies, the strength of the Brewer-Dobson circulation is found to increase in an idealized future climate; SSTs contribute most to this increase in the upper troposphere/lower stratosphere (UT/LS) region, while CO2 and ozone changes contribute most in the stratosphere and mesosphere.

Hurwitz, M. M.↗
NASA NTRS2022

Future Climate Change Under SSP Emission Scenarios with GISS-E2.1

This paper presents the response to anthropogenic forcing in the GISS-E2.1 climate models for the 21st century Shared Socioeconomic Pathways emission scenarios within the Coupled Model Intercomparison Project Phase 6 (CMIP6). The experiments were performed using an updated and improved version of the NASA Goddard Institute for Space Studies (GISS) coupled general circulation model that includes two different versions for atmospheric composition: A non-interactive version (NINT) with prescribed composition and a tuned aerosol indirect effect and the One-Moment Aerosol model (OMA) version with fully interactive aerosols which includes a parameterized first indirect aerosol effect on clouds. The effective climate sensitivities are 3.0°C and 2.9°C for the NINT and OMA models, respectively. Each atmospheric version is coupled to two different ocean general circulation models: The GISS ocean model (E2.1-G) and HYCOM (E2.1-H). We describe the global mean responses for all future scenarios and spatial patterns of change for surface air temperature and precipitation for four of the marker scenarios: SSP1-2.6, SSP2-4.5, SSP4-6.0, and SSP5-8.5. By 2100, global mean warming ranges from 1.5°C to 5.2°C relative to 1,850–1,880 mean temperature. Two high-mitigation scenarios SSP1-1.9 and SSP1-2.6 limit the surface warming to below 2°C by the end of the 21st century, except for the NINT E2.1-H model that simulates 2.2°C of surface warming. For the high emission scenario SSP5-8.5, the range is 4.6–5.2°C at 2100. Due to about 15% larger effective climate sensitivity and stronger transient climate response in both NINT and OMA CMIP6 models compared to CMIP5 versions, there is a stronger warming by 2100 in the SSP emission scenarios than in the comparable Representative Concentration Pathway (RCP) scenarios in CMIP5. Changes in sea ice area are highly correlated to global mean surface air temperature anomalies and show steep declines in both hemispheres, with the largest sea ice area decreases occurring during September in the Northern Hemisphere in both E2.1-G (−1.21 × 106 km2/°C) and E2.1-H models (−0.94 × 106 km2/°C). Both coupled models project decreases in the Atlantic overturning stream function by 2100. The largest decrease of 56%–65% in the 21st century overturning stream function is produced in the warmest scenario SSP5-8.5 in the E2.1-G model, comparable to the reduction in the corresponding CMIP5 GISS-E2 RCP8.5 simulation. Both low-end scenarios SSP1-1.9 and SSP1-2.6 also simulate substantial reductions of the overturning (9%–37%) with slow recovery of about 10% by the end of the 21st century (relative to the maximum decrease at the middle of the 21st century).

climate models↗
NASA NTRS2018

PACE Technical Report Series, Volume 4: Cloud Retrievals in the PACE Mission: PACE Science Team Consensus Document

Earth is a complex dynamical system exhibiting continuous change in its atmosphere, ocean,and surface elements. Nearly all (99.97%) of the energy driving these systems is linked to the Sun. Measurements of reflected sunlight contain a unique signature of wavelength-specific scattering and absorption interactions occurring between incoming solar energy and atmospheric (molecules, aerosols,clouds) and surface features Clouds can affect significantly both shortwave and long wave radiation, depending on altitude/vertical structure, thermodynamic phase, and optical properties. Low, warm, and optically thick clouds predominantly have a cooling effect, while high, cold, optically thin clouds can cause warming by absorbing warmer radiation emitted from the surface and lower atmosphere.When the net difference between outgoing and incoming solar radiation is matched by the net infrared radiation emitted to space, the Earth's climate is in radiative balance. While radiative forcing components (GHGs, aerosols - direct and indirect) contribute to a net radiative imbalance, climate sensitivity is ultimately determined by the contribution of various system feed backs. The role of cloud feedback in a warming climate is currently the largest inter-model uncertainty in climate sensitivity and therefore in climate prediction [Bony and Dufresne 2005]. A comprehensive understanding of current cloud propertiesand dynamic/microphysical processes requires a global perspective from satellites.

Plankton, Aerosol, Cloud, ocean Ecosystem (PACE)↗
NASA NTRS2021

Snow Reconciles Observed and Simulated Phase Partitioning and Increases Cloud Feedback

The surprising increase of Earth's climate sensitivity in the most recent Coupled Model Intercomparison Project (CMIP) models has been largely attributed to extratropical cloud feedback, which is thought to be driven by greater supercooled water in present-day cloud phase partitioning (CPP). Here we report that accounting for precipitation in the Goddard Institute for Space Studies ModelE3 radiation scheme, neglected in more than 60% of CMIP6 and 90% of CMIP5 models, systematically changes its apparent CPP and substantially increases its cloud feedback, consistent with results using CMIP models. Including precipitation in the comparison with Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) measurements and in model radiation schemes is essential to faithfully constrain cloud amount and phase partitioning, and simulate cloud feedbacks. Our findings suggest that making radiation schemes precipitation-aware (missing in most CMIP6 models) should strengthen their positive cloud feedback and further increase their already high mean climate sensitivity.

Earth's climate sensitivity↗
NASA NTRS2017

Enhanced Groundwater Recharge Rates and Altered Recharge Sensitivity to Climate Variability Through Subsurface Heterogeneity

Our environment is heterogeneous. In hydrological sciences, the heterogeneity of subsurface properties, such as hydraulic conductivities or porosities, exerts an important control on water balance. This notably includes groundwater recharge, which is an important variable for efficient and sustainable groundwater resources management. Current large-scale hydrological models do not adequately consider this subsurface heterogeneity. Here we show that regions with strong subsurface heterogeneity have enhanced present and future recharge rates due to a different sensitivity of recharge to climate variability compared with regions with homogeneous subsurface properties. Our study domain comprises the carbonate rock regions of Europe, Northern Africa, and the Middle East, which cover 25 of the total land area. We compare the simulations of two large-scale hydrological models, one of them accounting for subsurface heterogeneity. Carbonate rock regions strongly exhibit karstification, which is known to produce particularly strong subsurface heterogeneity. Aquifers from these regions contribute up to half of the drinking water supply for some European countries. Our results suggest that water management for these regions cannot rely on most of the presently available projections of groundwater recharge because spatially variable storages and spatial concentration of recharge result in actual recharge rates that are up to four times larger for present conditions and changes up to five times larger for potential future conditions than previously estimated. These differences in recharge rates for strongly heterogeneous regions suggest a need for groundwater management strategies that are adapted to the fast transit of water from the surface to the aquifers.

water resources↗