Search NASA⌕ Search

SEARCH · Search NASA

Results for “carbon accounting”

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 127 records · Page 7

Aerosol Absorption in Global Models From AeroCom Phase III

Aerosol-induced absorption of shortwave radiation can modify the climate through local atmospheric heating, which affects lapse rates, precipitation, and cloud formation. Presently, the total amount of aerosol absorption is poorly constrained, and the main absorbing aerosol species (black carbon (BC), organic aerosols (OA), and mineral dust) are diversely quantified in global climate models. As part of the third phase of the Aerosol Comparisons between Observations and Models (AeroCom) intercomparison initiative (AeroCom phase III), we here document the distribution and magnitude of aerosol absorption in current global aerosol models and quantify the sources of intermodel spread, highlighting the difficulties of attributing absorption to different species. In total, 15 models have provided total present-day absorption at 550 nm (using year 2010 emissions), 11 of which have provided absorption per absorbing species. The multi-model global annual mean total absorption aerosol optical depth (AAOD) is 0.0054 (0.0020 to 0.0098; 550 nm), with the range given as the minimum and maximum model values. This is 28 % higher compared to the 0.0042 (0.0021 to 0.0076) multi-model mean in AeroCom phase II (using year 2000 emissions), but the difference is within 1 standard deviation, which, in this study, is 0.0023 (0.0019 in Phase II). Of the summed component AAOD, 60 % (range 36 %–84 %) is estimated to be due to BC, 31 % (12 %–49 %) is due to dust, and 11 % (0 %–24 %) is due to OA; however, the components are not independent in terms of their absorbing efficiency. In models with internal mixtures of absorbing aerosols, a major challenge is the lack of a common and simple method to attribute absorption to the different absorbing species. Therefore, when possible, the models with internally mixed aerosols in the present study have performed simulations using the same method for estimating absorption due to BC, OA, and dust, namely by removing it and comparing runs with and without the absorbing species. We discuss the challenges of attributing absorption to different species; we compare burden, refractive indices, and density; and we contrast models with internal mixing to models with external mixing. The model mean BC mass absorption coefficient (MAC) value is 10.1 (3.1 to 17.7) m 2 g −1 (550 nm), and the model mean BC AAOD is 0.0030 (0.0007 to 0.0077). The difference in lifetime (and burden) in the models explains as much of the BC AAOD spread as the difference in BC MAC values. The difference in the spectral dependency between the models is striking. Several models have an absorption Ångstrøm exponent (AAE) close to 1, which likely is too low given current knowledge of spectral aerosol optical properties. Most models do not account for brown carbon and underestimate the spectral dependency for OA.

Aerosol absorption↗

Carbon isotope effects associated with autotrophic acetogenesis

The carbon kinetic isotope effects associated with synthesis of acetate from CO2 and H2 during autotrophic growth of Acetobacterium woodii at 30 degrees C have been measured by isotopic analyses of CO2, methyl-carbon, and total acetate. Closed systems allowing construction of complete mass balances at varying stages of growth were utilized, and the effects of the partitioning of carbon between CO2 and HCO3- were taken account. For the overall reaction, total carbonate --> total acetate, isotope effects measured in replicate experiments ranged from -59.0 +/- 0.9% to -57.2 +/- 2.3%. Taking into account all measurements, the weighted mean and standard deviation are -58.6 +/- 0.7%. There is no evidence for intramolecular ordering in the acetate. The carbon isotopic composition of sedimentary acetate, otherwise expected to be near that of sedimentary organic carbon, is likely to be depleted in environments in which autotrophic acetogenesis is occurring.

Non-NASA Center↗

Impact of a dynamic grid mix and climate on operational carbon emissions modeling for different building typologies and climate zones

Calculating operational carbon emissions through a building’s lifecycle is complex due to the dynamic nature of influencing factors such as climate and energy grid mix. This paper introduces a novel methodology for modeling 30-year operational carbon impacts of buildings and applies this method to mid-rise office and residential typologies across various US climate zones. The method accounts for these temporal variabilities using new and scarcely cited data sources. Key findings indicate that future changes in the climate, while impactful, play a relatively modest role in operational carbon emissions compared to significant reductions with modeling scenarios using the projected decarbonization of the electricity grid. Here, the study also finds that using annual, month-hourly, or hourly grid emission factors have a minimal impact on carbon accounting, except in certain climates and program types where emission patterns do not align with a building’s energy consumption. Warmer climates like Miami, Florida and Tucson, Arizona, which rely heavily on cooling, demonstrate larger variations in carbon emissions when using higher temporal resolution emission factors. Ultimately, this study underscores the critical role of grid decarbonization in reducing long-term emissions and the importance of incorporating this variable in life cycle assessment (LCA) modeling.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The cool DC white dwarf Stein 2051B

Visual scans from the Kitt-Peak 2.1-m telescope and UV spectra from the IUE satellite are combined in a spectroscopic study of the bright DC white dwarf, Stein 2051B. With a cooling age near 2 billion years, and a resulting atmospheric mass of hydrogen of the order of 6 x 10 to the 17th g, the accretion rate over the star's lifetime is less than 2 x 10 to the -25th solar masses/year. This is far below the predictions of standard accretion theories, and requires unreasonably strong suppression of hydrogen accretion in this star if the accretion-diffusion theory is to account for the carbon abundance found. It is suggested that the carbon may have been brought to the surface by convective mixing.

Wegner, G.↗

Carbon chain abundance in the diffuse interstellar medium

Thanks to the mid-IR sensitivities of the ISO and IRTS orbiting spectrometers it is now possible to search the diffuse interstellar medium for heretofore inaccessible molecular emission. In view of the recent strong case for the presence of C(7-) (Kirkwood et al. 1998, Tulej et al. 1998),and the fact that carbon chains possess prominent infrared active modes in a very clean portion of the interstellar spectrum, we have analyzed the IRTS spectrum of the diffuse interstellar medium for the infrared signatures of these species. Theoretical and experimental infrared band frequencies and absolute intensities of many different carbon chain species are presented. These include cyanopolyynes, neutral and anionic linear carbon molecules, and neutral and ionized, even-numbered, hydrogenated carbon chains. We show that--as a family--these species have abundances in the diffuse ISM on the order of 10(-10) with respect to hydrogen, values consistent with their abundances in dense molecular clouds. Assuming an average length of 10 C atoms per C-chain implies that roughly a millionth of the cosmically available carbon is in the form of carbon chains and that carbon chains can account for a few percent of the visible to near-IR diffuse interstellar band (DIB) total equivalent width (not DIB number).

Carbon/analysis/chemistry↗

The Impacts of Fossil Fuel Emission Uncertainties and Accounting for 3-D Chemical CO2 Production on Inverse Natural Carbon Flux Estimates from Satellite and In Situ Data

Atmospheric carbon dioxide (CO2) inversions for estimating natural carbon fluxes typically do not allow for adjustment of fossil fuel CO2 emissions, despite significant uncertainties in emission inventories and inadequacies in the specification of international bunker emissions in inversions. Also, most inversions place CO2 release from fossil fuel combustion and biospheric sources entirely at the surface. However, a non-negligible portion of the emissions actually occurs in the form of reduced carbon species, which are eventually oxidized to CO2 downwind. Omission of this ‘chemical pump’ can result in a significant redistribution of the inferred total carbon fluxes among regions. We assess the impacts of different prescriptions of fossil fuel emissions and accounting for the chemical pump on flux estimation, with a novel aspect of conducting both satellite CO2 observation-based and surface in situ-based inversions. We apply 3-D carbon monoxide (CO) loss rates archived from a state-of-the-art GEOS chemistry and climate model simulation in a forward transport model run to simulate the distribution of CO2 originating from oxidation of carbon species. We also subtract amounts from the prior surface CO2 fluxes that are actually emitted in the form of fossil and biospheric CO, methane, and nonmethane volatile organic compounds (VOCs). We find that the posterior large-scale fluxes are generally insensitive to the finer-scale spatial differences between the ODIAC and CDIAC fossil fuel CO2 gridded datasets and assumptions about international bunker emissions. However, accounting for 3-D chemical CO2 production and the surface correction shifts the global carbon sink, e.g., from land to ocean and from the tropics to the north, with a magnitude and even direction that depend on assumptions about the surface correction. A GOSAT satellite-based inversion is more sensitive to the chemical pump than one using in situ observations, exhibiting substantial flux impacts of 0.28, 0.53, and -0.47 Pg C yr(exp -1) over tropical land, global land, and oceans, due to differences in the horizontal and vertical sampling of the wo observation types. Overall, the biases from neglecting the chemical pump appear to be minor relative to the flux estimate uncertainties and the differences between the in situ and GOSAT inversions, but their relative importance will grow in the future as observational coverage further increases and satellite retrieval biases decrease.

James S Wang↗

A life cycle and product type based estimator for quantifying the carbon stored in wood products

Background: Timber harvesting and industrial wood processing laterally transfer the carbon stored in forest sectors to wood products creating a wood products carbon pool. The carbon stored in wood products is allocated to end-use wood products (e.g., paper, furniture), landfill, and charcoal. Wood products can store substantial amounts of carbon and contribute to the mitigation of greenhouse effects. Therefore, accurate accounts for the size of wood products carbon pools for different regions are essential to estimating the land-atmosphere carbon exchange by using the bottom-up approach of carbon stock change. Results: To quantify the carbon stored in wood products, we developed a state-of-the-art estimator (Wood Products Carbon Storage Estimator, WPsCS Estimator) that includes the wood products disposal, recycling, and waste wood decomposition processes. The wood products carbon pool in this estimator has three subpools: (1) end-use wood products, (2) landfill, and (3) charcoal carbon. In addition, it has a user-friendly interface, which can be used to easily parameterize and calibrate an estimation. To evaluate its performance, we applied this estimator to account for the carbon stored in wood products made from the timber harvested in Maine, USA, and the carbon storage of wood products consumed in the United States. Conclusion: The WPsCS Estimator can efficiently and easily quantify the carbon stored in harvested wood products for a given region over a specific period, which was demonstrated with two illustrative examples. In addition, WPsCS Estimator has a user-friendly interface, and all parameters can be easily modified.

54 ENVIRONMENTAL SCIENCES↗

How Much Blue Carbon do Hudson Estuary Marshes Sequester? A Full-Depth Carbon Stock Estimation of Iona Marsh, Lower Hudson, NY

Tidal marshes function as a key ecosystem, by serving as nurseries, improving water quality, protecting biodiversity, preventing erosion of coastlines, and preserving archives of paleoclimate and paleoecology. In the last few decades, researchers have come to appreciate coastal marshes’ carbon sequestration function. However, full peat depth is required for estimating the true carbon content of marshes, and very little research has been done to discover the full peat depth along coastlines. We selected Iona Marsh (73°58’44.972”W, 41°18’6.366”N), situated adjacent to Bear Mountain Park and part of the National Estuarine Research Reserve, to document the belowground storage of blue carbon. In order to calculate marsh volume, three full-depth transects were performed from 2020 to 2022, two from east to west in different regions of the marsh, and one from north to south at a minimum of 10 m intervals. We found that marsh depth ranged from 5.1 m to 13.2 m. We calculated marsh volume from the full peat depth and marsh area, obtained through the National Wetlands Inventory database, by applying the average marsh depth over its entire surface area. In order to calculate carbon content, we extracted a 9m sediment core and used loss-on-ignition (LOI) to determine organic matter, then estimate carbon content as half of the organic matter. We used 198 LOI samples to calculate the carbon, which formed at this site over 6000 years ago. We intend to apply this methodology framework to other marshes in the Hudson region. Because many previous carbon stock estimates only account for the top 1 m of peat, the role of these ecosystems in carbon sequestration is significantly underestimated. By determining the full-depth carbon stock of these marshes, we hope to emphasize the value of marshes to the community, beyond what is traditionally recognized, particularly in light of sea level rise and the vulnerability of marshes to erosion.

Tidal marshes↗

Carbon suboxide in Comet Halley?

The extremely dark nucleus of Comet Halley suggests that its volatile ices contain a few percent of carbonaceous material in the form of graphitic or amorphous carbon. The very high abundance of light elements in the coma dust and the emission feature near 3.4 microns further suggest the presence of a significant organic component, although the identified carbon-containing materials' parent species cannot account for all of such a component. It is presently proposed that an additional contribution from carbon suboxide can account for these observational data, assuming a production rate about 0.03-0.04 times that of water.

Huntress, Wesley T., Jr.↗

Model‐Based Interpretation of Solute Exports and Carbon Partitioning During Shale Weathering in a Mountainous Hillslope

The weathering of sedimentary rocks in high-elevation catchments influences freshwater quality and the global carbon cycle. While individual biogeochemical mechanisms involved in this process are relatively well understood, quantifying their contributions to solute export and carbon fluxes under natural, transient conditions remains challenging. Here, we implement a numerical multidimensional and multiphase model to simulate coupled hydrological and biogeochemical processes in a shale-underlain, snow-dominated hillslope in the Rocky Mountains, Colorado. The model captures the dynamic interplay between soil respiration, mineral weathering, and climate-driven hydrological forcing, reproducing observed soil CO 2 dynamics, groundwater chemistry, and subsurface flow. Our results reveal that seasonal snowmelt enhances carbonate weathering by promoting the infiltration of CO 2 -rich water to depth, while pyrite oxidation is primarily sensitive to low water saturation that facilitates O 2 diffusion through the regolith. Topography modulates the spatial distribution of shale weathering, as steeper slopes enhance lateral drainage, favoring the delivery of reactants to greater depths. While shale weathering at our site acts as a transient carbon sink, with silicates and carbonates buffering acidity and promoting atmospheric CO 2 consumption (1% of soil-derived CO 2 ), the exported dissolved inorganic carbon is predominantly geogenic (∼73%). Consequently, when accounting for long-term marine carbonate precipitation. The current weathering regime represents a net source of carbon to the atmosphere. The oxidation of pyrite and petrogenic organic carbon together release approximately 0.9 mol·m −2 ·yr −1 of CO 2 . Our findings highlight the role of topography, hydroclimate, and the coupling between acid-base reactions in shaping the carbon balance and the solute exports in mountainous critical zones.

carbon cycling↗

Planetary noble gas components in Orgueil

In the reported investigation, the approach employed by Jeffery and Anders (1970) and Eberhardt (1974), who used colloidal suspensions in separating minerals of the Orgueil carbonaceous chondrite (C1), was extended to an application involving the enrichment of a planetary noble gas carrier in several C1 and C2 meteorites. The thermal release pattern of Orgueil 'carbon' is considered, taking into account elemental effects, a comparison of noble gas release from 'carbon' separates and bulk meteorite samples, isotopic effects, and a comparison with Murray 'carbon'. An evaluation is conducted of the noble gas components in carbonaceous chondrites. The origin of primordial noble gases in carbonaceous chondrites is also discussed, giving attention to the origin of solar gases, the role of 'carbon' as the novel vehicle for planetary gases, and plausible astrophysical sites for the production of anomalous krypton and xenon.

Frick, U.↗

Energy Emissions Accounting Methods Can Determine Whether Direct Air Capture with Storage Achieves Net Removal

The voluntary carbon market within the United States has expanded rapidly in recent years and enabled private companies and other organizations to provide revenue streams to carbon dioxide removal (CDR) technologies. For a CDR technology to participate in the voluntary carbon market (VCM), the emissions associated with constructing and operating the technology must be less than the CO 2 captured from the atmosphere. Assessing the extent to which this is true for direct air capture with storage (DACS), a relatively energy-intensive CDR technology, strongly depends on the accounting method used to assess the emissions intensity of purchased energy. We simulate the hourly weather-dependent operation of sorbent- and solvent-based DACS in California, Louisiana, Texas, and Wyoming, representing a wide range of local weather and electric and natural gas grid compositions. In all cases, the single most important emissions accounting decision is the method used to estimate the emissions intensity of purchased grid electricity, which varies the calculated net removal by −1049% to +108%. All other factors influencing net removal introduce a variation of at most ±14%. No electricity emissions accounting method is universally conservative across all scenarios, and none is objectively more accurate. High-spatiotemporal-resolution, high-quality, publicly available data sets and models for electricity emissions accounting do not currently exist and are urgently needed to enable standardization of emissions accounting methods to more accurately determine the true emissions impacts of DACS and other energy-intensive facilities.

54 ENVIRONMENTAL SCIENCES↗

Uncertainty in land carbon budget simulated by terrestrial biosphere models: the role of atmospheric forcing

Global estimates of the land carbon sink are often based on simulations by terrestrial biosphere models (TBMs). The use of a large number of models that differ in their underlying hypotheses, structure and parameters is one way to assess the uncertainty in the historical land carbon sink. Here we show that the atmospheric forcing datasets used to drive these TBMs represent a significant source of uncertainty that is currently not systematically accounted for in land carbon cycle evaluations. We present results from three TBMs each forced with three different historical atmospheric forcing reconstructions over the period 1850–2015. We perform an analysis of variance to quantify the relative uncertainty in carbon fluxes arising from the models themselves, atmospheric forcing, and model-forcing interactions. We find that atmospheric forcing in this set of simulations plays a dominant role on uncertainties in global gross primary productivity (GPP) (75% of variability) and autotrophic respiration (90%), and a significant but reduced role on net primary productivity and heterotrophic respiration (30%). Atmospheric forcing is the dominant driver (52%) of variability for the net ecosystem exchange flux, defined as the difference between GPP and respiration (both autotrophic and heterotrophic respiration). In contrast, for wildfire-driven carbon emissions model uncertainties dominate and, as a result, model uncertainties dominate for net ecosystem productivity. At regional scales, the contribution of atmospheric forcing to uncertainty shows a very heterogeneous pattern and is smaller on average than at the global scale. We find that this difference in the relative importance of forcing uncertainty between global and regional scales is related to large differences in regional model flux estimates, which partially offset each other when integrated globally, while the flux differences driven by forcing are mainly consistent across the world and therefore add up to a larger fractional contribution to global uncertainty.

54 ENVIRONMENTAL SCIENCES↗

Intracellular pathways for lignin catabolism in white-rot fungi

White-rot fungi, the most efficient organisms at breaking down lignin from wood in Nature, utilize lignin degradation products as a carbon source. 15 Lignin is a biopolymer found in plant cell walls that accounts for 30% of the organic carbon in the biosphere. White-rot fungi (WRF) are considered the most efficient organisms at degrading lignin in Nature. While lignin depolymerization by WRF has been exhaustively studied, the possibility that WRF are able to utilize lignin as a carbon source is still a matter of controversy. Here we employ 13C-labeling and systems biology approaches to demonstrate that two WRF, Trametes versicolor and Gelatoporia subvermispora, funnel lignin-derived aromatic compounds into central carbon metabolism via intracellular catabolic pathways. These results provide insights into global carbon cycling in 20 soil ecosystems, and furthermore establishes a foundation for employing WRF in simultaneous lignin depolymerization and bioconversion to bioproducts – a key step towards enabling a sustainable bioeconomy.

Del Cerro, Carlos↗

A Global Modeling Study on Carbonaceous Aerosol Microphysical Characteristics and Radiative Effects

Recently, attention has been drawn towards black carbon aerosols as a short-term climate warming mitigation candidate. However the global and regional impacts of the direct, indirect and semi-direct aerosol effects are highly uncertain, due to the complex nature of aerosol evolution and the way that mixed, aged aerosols interact with clouds and radiation. A detailed aerosol microphysical scheme, MATRIX, embedded within the GISS climate model is used in this study to present a quantitative assessment of the impact of microphysical processes involving black carbon, such as emission size distributions and optical properties on aerosol cloud activation and radiative effects. Our best estimate for net direct and indirect aerosol radiative flux change between 1750 and 2000 is -0.56 W/m2. However, the direct and indirect aerosol effects are quite sensitive to the black and organic carbon size distribution and consequential mixing state. The net radiative flux change can vary between -0.32 to -0.75 W/m2 depending on these carbonaceous particle properties at emission. Taking into account internally mixed black carbon particles let us simulate correct aerosol absorption. Absorption of black carbon aerosols is amplified by sulfate and nitrate coatings and, even more strongly, by organic coatings. Black carbon mitigation scenarios generally showed reduced radiative fluxeswhen sources with a large proportion of black carbon, such as diesel, are reduced; however reducing sources with a larger organic carbon component as well, such as bio-fuels, does not necessarily lead to a reduction in positive radiative flux.

Bauer, S. E.↗

A Comprehensive Greenhouse Gas Assessment of Biomass-Based Carbon Dioxide Removal and Storage

Biomass-based Carbon Removal and Storage (BiCRS) is a suite of technological pathways that provide a relatively affordable and high-potential way to remove carbon dioxide from the atmosphere and simultaneously re-purpose residual biomass. BiCRS technologies are individually optimized to different biomass feedstocks, temperatures and oxygen, with the primary goal of capturing a high percentage of biomass-derived carbon dioxide and storing it in permanent below-ground storage (i.e. geologic storage), and a secondary goal of producing renewable energy. The current standard for BiCRS life cycle assessments does not account for greenhouse gas (carbon dioxide, methane and nitrous oxide) fluxes from surface soil amendment of BiCRS byproducts (e.g. char, ash) and their alternate fates. This is partly due to a paucity of empirical data due to the novelty of BiCRS conversion technologies, the range in feedstocks, and the heterogeneity of soils to which they might be amended. To fill this gap in in the life-cycle assessment and provide a realistic range of parameters for a soil amendment component of the BiCRS life cycle assessment, we conducted an incubation experiment to measure soil carbon changes, microbial respiration, methane fluxes, and nitrous oxide fluxes from two different soils amended with char- and ash- byproducts from biomass after gasification to hydrogen, fast pyrolysis to bio-oil, and torrefaction. Our results suggest that nitrous oxide and methane emission or consumption from BiCRS char amendments to soil are context dependent. Nitrous oxide emissions from amendments are higher in soils with higher pH, whereas some char amendments to low pH soil can reduce N 2 O emissions relative to control. In contrast, methane is emitted from BiCRS char amended to low pH soil, but consumed or neutral relative to control when amended to high pH soil. We present carbon and nitrogen mass balance throughout the experiment to help understand short-term durability of BiCRS chars and their counterfactuals, and preliminary suggestions for revisions to the broader BiCRS Measurement, Reporting, and Verification framework.

54 ENVIRONMENTAL SCIENCES↗

Consistency Between Sun-Induced Chlorophyll Fluorescence and Gross Primary Production of Vegetation in North America

Accurate estimation of the gross primary production (GPP) of terrestrial ecosystems is vital for a better understanding of the spatial-temporal patterns of the global carbon cycle. In this study,we estimate GPP in North America (NA) using the satellite-based Vegetation Photosynthesis Model (VPM), MODIS (Moderate Resolution Imaging Spectrometer) images at 8-day temporal and 500 meter spatial resolutions, and NCEP-NARR (National Center for Environmental Prediction-North America Regional Reanalysis) climate data. The simulated GPP (GPP (sub VPM)) agrees well with the flux tower derived GPP (GPPEC) at 39 AmeriFlux sites (155 site-years). The GPP (sub VPM) in 2010 is spatially aggregated to 0.5 by 0.5-degree grid cells and then compared with sun-induced chlorophyll fluorescence (SIF) data from Global Ozone Monitoring Instrument 2 (GOME-2), which is directly related to vegetation photosynthesis. Spatial distribution and seasonal dynamics of GPP (sub VPM) and GOME-2 SIF show good consistency. At the biome scale, GPP (sub VPM) and SIF shows strong linear relationships (R (sup 2) is greater than 0.95) and small variations in regression slopes ((4.60-5.55 grams Carbon per square meter per day) divided by (milliwatts per square meter per nanometer per square radian)). The total annual GPP (sub VPM) in NA in 2010 is approximately 13.53 petagrams Carbon per year, which accounts for approximately 11.0 percent of the global terrestrial GPP and is within the range of annual GPP estimates from six other process-based and data-driven models (11.35-22.23 petagrams Carbon per year). Among the seven models, some models did not capture the spatial pattern of GOME-2 SIF data at annual scale, especially in Midwest cropland region. The results from this study demonstrate the reliable performance of VPM at the continental scale, and the potential of SIF data being used as a benchmark to compare with GPP models.

photosynthesis model↗

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↗