Search NASA⌕ Search

SEARCH · Search NASA

Results for “climate model outputs”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

BSEC VPRM 10m Hourly Biogenic Fluxes in Baltimore (2021)

Model outputs from the Vegetation Photosynthesis and Respiration Model (VPRM: version from Horne et al. in prep). Model remote sensing inputs come from Sential 2-derived EVI and LSWI. Model meteorological inputs for two-meter air temperature and shortwave incoming come from the BSEC WRF 2021 Control Run (Foust, W. 2023). Plant functional Types (PFTs) are spatially classified using the Chesapeake Bay Program 2018 land use land cover product. The final biogenic flux (µmol CO2 m^-2 s^-1) outputs of NEE, RESP, and GEE are a weighted average based on the portion of PFTs within the cell. Individual PFT outputs are saved inside PFT directories (e.g., Crops, Grass, etc.) inside the specific month directory. Model outputs are denoted as a negative flux into the land system (i.e., photosynthesis) and a positive flux as a net release into the overlying atmosphere. Respiration (RESP) fluxes are positive and combine heterotrophic (only soil) and autotrophic sources. Gross ecosystem exchange (GEE) is a negative flux driven by only photosynthetic activity from vegetation, and the Net ecosystem exchange (NEE) is the sum of the two (i.e., NEE=RESP+GEE). Data Characteristics Spatial Resolution: 10m Temporal Resolution: Hourly File Format: VPRM_ _BSEC. .tif (Hour is in UTC) For more information on the model results, please email Jason Horne (jph6488@psu.edu). References: Foust, W. (2023). BSEC WRF 2021 Control Run Output (v0.1.0) [Data set]. MSD-LIVE Data Repository. https://data.msdlive.org/records/m0e6m-vvq17

Baltimore↗
DOE OSTI2023

Scoping review of Culex mosquito life history trait heterogeneity in response to temperature

Abstract Background Mosquitoes in the genus Culex are primary vectors in the US for West Nile virus (WNV) and other arboviruses. Climatic drivers such as temperature have differential effects on species-specific changes in mosquito range, distribution, and abundance, posing challenges for population modeling, disease forecasting, and subsequent public health decisions. Understanding these differences in underlying biological dynamics is crucial in the face of climate change. Methods We collected empirical data on thermal response for immature development rate, egg viability, oviposition, survival to adulthood, and adult lifespan for Culex pipiens, Cx. quinquefasciatus, Cx. tarsalis , and Cx. restuans from existing literature according to the PRISMA scoping review guidelines. Results We observed linear relationships with temperature for development rate and lifespan, and nonlinear relationships for survival and egg viability, with underlying variation between species. Optimal ranges and critical minima and maxima also appeared varied. To illustrate how model output can change with experimental input data from individual Culex species, we applied a modified equation for temperature-dependent mosquito type reproduction number for endemic spread of WNV among mosquitoes and observed different effects. Conclusions Current models often input theoretical parameters estimated from a single vector species; we show the need to implement the real-world heterogeneity in thermal response between species and present a useful data resource for researchers working toward that goal. Graphical Abstract

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Spatiotemporal evolution of melt ponds on Arctic sea ice

Melt ponds on sea ice play an important role in the Arctic climate system. Their presence alters the partitioning of solar radiation: decreasing reflection, increasing absorption and transmission to the ice and ocean, and enhancing melt. The spatiotemporal properties of melt ponds thus modify ice albedo feedbacks and the mass balance of Arctic sea ice. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition presented a valuable opportunity to investigate the seasonal evolution of melt ponds through a rich array of atmosphere-ice-ocean measurements across spatial and temporal scales. In this study, we characterize the seasonal behavior and variability in the snow, surface scattering layer, and melt ponds from spring melt to autumn freeze-up using in situ surveys and auxiliary observations. We compare the results to satellite retrievals and output from two models: the Community Earth System Model (CESM2) and the Marginal Ice Zone Modeling and Assimilation System (MIZMAS). During the melt season, the maximum pond coverage and depth were 21% and 22 ± 13 cm, respectively, with distribution and depth corresponding to surface roughness and ice thickness. Compared to observations, both models overestimate melt pond coverage in summer, with maximum values of approximately 41% (MIZMAS) and 51% (CESM2). This overestimation has important implications for accurately simulating albedo feedbacks. During the observed freeze-up, weather events, including rain on snow, caused high-frequency variability in snow depth, while pond coverage and depth remained relatively constant until continuous freezing ensued. Both models accurately simulate the abrupt cessation of melt ponds during freeze-up, but the dates of freeze-up differ. MIZMAS accurately simulates the observed date of freeze-up, while CESM2 simulates freeze-up one-to-two weeks earlier. This work demonstrates areas that warrant future observation-model synthesis for improving the representation of sea-ice processes and properties, which can aid accurate simulations of albedo feedbacks in a warming climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Data and Code for: Observation-constrained agroecosystem model inversion reveals continental-scale variation of winter wheat traits

This repository contains the simulation outputs and processing scripts associated with the study of winter wheat traits across the United States, utilizing the Ecosys agroecosystem model. The dataset includes model results for both rainfed and irrigated winter wheat systems, supporting the findings presented in the manuscript titled "Observation-constrained agroecosystem model inversion reveals continental-scale variation of winter wheat traits." Data includes the original Ecosys simulation outputs (archived in .db format within the compressed .zip files) and extracted analysis data (stored in .pkl files for efficient processing). Python code for data processing and figure generation is provided in a Jupyter notebook. External Observational Datasets should refer to the following official repositories for the input and validation data used in this study. The eddy covariance data from the AmeriFlux network (https://ameriflux.lbl.gov/). Climate-forcing data of NLDAS-2 from NASA LDAS (https://ldas.gsfc.nasa.gov/nldas/nldas-2-forcing-data). Soil data from the Gridded Soil Survey Geographic Database (gSSURGO), available at (https://www.nrcs.usda.gov/resources/data-and-reports/gridded-soil-survey-geographic-gssurgo-database). Crop yields, planting and harvest dates from the USDA public databases (https://quickstats.nass.usda.gov/; https://webapp.rma.usda.gov/apps/actuarialinformationbrowser/CropCriteria.aspx). Satellite-derived SLOPE GPP data from ORNL DAAC (https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=1786). Land use and crop progress information from the USDA Crop Data Layer and Crop Progress and Condition Gridded Layers (https://www.nass.usda.gov/Research_and_Science/). The Ecosys model code is available online at https://github.com/jinyun1tang/ECOSYS.

Wheat↗
DOE OSTI2023

Little Change in Apparent Hydrological Sensitivity at Large CO 2 Forcing

Apparent hydrological sensitivity (η a ), the change in the global mean precipitation per degree K of global surface warming, is a key aspect of the climate system's response to increasing CO 2 forcing. To determine whether ηa depends on the forcing amplitude we analyze idealized experiments over a broad range of abrupt CO 2 forcing, from 2× to 8× preindustrial values, with two distinct climate models. We find little change in η a between 2× and 4×CO 2 , and almost no change beyond 5×CO 2 . We validate this finding under transient CO 2 forcing at 1%-per-year, up to 8×CO 2 . We further corroborate this result by analyzing the 1%-per-year output of more than 15 CMIP5/6 models. Lastly, we examine the 1,000-year long LongrunMIP model output, and again find little change in η a . This wealth of results demonstrates that η a is a very weak function of CO 2 forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Revealing the Statistics of Extreme Events Hidden in Short Weather Forecast Data

Extreme weather events have significant consequences, dominating the impact of climate on society. While high-resolution weather models can forecast many types of extreme events on synoptic timescales, long-term climatological risk assessment is an altogether different problem. A once-in-a-century event takes, on average, 100 years of simulation time to appear just once, far beyond the typical integration length of a weather forecast model. Therefore, this task is left to cheaper, but less accurate, low-resolution or statistical models. But there is untapped potential in weather model output: despite being short in duration, weather forecast ensembles are produced multiple times a week. Integrations are launched with independent perturbations, causing them to spread apart over time and broadly sample phase space. Collectively, these integrations add up to thousands of years of data. We establish methods to extract climatological information from these short weather simulations. Using ensemble hindcasts by the European Center for Medium-range Weather Forecasting archived in the subseasonal-to-seasonal (S2S) database, we characterize sudden stratospheric warming (SSW) events with multi-centennial return times. Consistent results are found between alternative methods, including basic counting strategies and Markov state modeling. By carefully combining trajectories together, we obtain estimates of SSW frequencies and their seasonal distributions that are consistent with reanalysis-derived estimates for moderately rare events, but with much tighter uncertainty bounds, and which can be extended to events of unprecedented severity that have not yet been observed historically. These methods hold potential for assessing extreme events throughout the climate system, beyond this example of stratospheric extremes.

58 GEOSCIENCES↗
DOE OSTI2022

Subseasonal Earth System Prediction with CESM2

Abstract Prediction systems to enable Earth system predictability research on the subseasonal time scale have been developed with the Community Earth System Model, version 2 (CESM2) using two configurations that differ in their atmospheric components. One system uses the Community Atmosphere Model, version 6 (CAM6) with its top near 40 km, referred to as CESM2(CAM6). The other employs the Whole Atmosphere Community Climate Model, version 6 (WACCM6) whose top extends to ∼140 km, and it includes fully interactive tropospheric and stratospheric chemistry [CESM2(WACCM6)]. Both systems are utilized to carry out subseasonal reforecasts for the 1999–2020 period following the Subseasonal Experiment’s (SubX) protocol. Subseasonal prediction skill from both systems is compared to those of the National Oceanic and Atmospheric Administration CFSv2 and European Centre for Medium-Range Weather Forecasts (ECMWF) operational models. CESM2(CAM6) and CESM2(WACCM6) show very similar subseasonal prediction skill of 2-m temperature, precipitation, the Madden–Julian oscillation, and North Atlantic Oscillation to its previous version and to the NOAA CFSv2 model. Overall, skill of CESM2(CAM6) and CESM2(WACCM6) is a little lower than that of the ECMWF system. In addition to typical output provided by subseasonal prediction systems, CESM2 reforecasts provide comprehensive datasets for predictability research of multiple Earth system components, including three-dimensional output for many variables, and output specific to the mesosphere and lower-thermosphere (MLT) region from CESM2(WACCM6). It is shown that sudden stratosphere warming events, and the associated variability in the MLT, can be predicted ∼10 days in advance. Weekly real-time forecasts and reforecasts with CESM2(CAM6) and CESM2(WACCM6) are freely available. Significance Statement We describe here the design and prediction skill of two subseasonal prediction systems based on two configurations of the Community Earth System Model, version 2 (CESM2): CESM2 with the Community Atmosphere Model, version 6 [CESM2(CAM6)] and CESM 2 with Whole Atmosphere Community Climate Model, version 6 [CESM2(WACCM6)] as its atmospheric component. These two systems provide a foundation for community-model based subseasonal prediction research. The CESM2(WACCM6) system provides a novel capability to explore the predictability of the stratosphere, mesosphere, and lower thermosphere. Both CESM2(CAM6) and CESM2(WACCM6) demonstrate subseasonal surface prediction skill comparable to that of the NOAA CFSv2 model, and a little lower than that of the ECMWF forecasting system. CESM2 reforecasts provide a comprehensive dataset for predictability research of multiple aspects of the Earth system, including the whole atmosphere up to 140 km, land, and sea ice. Weekly real-time forecasts, reforecasts, and models are publicly available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Biophysical Impact of Land-Use and Land-Cover Change on Subgrid Temperature in CMIP6 Models

In this study, we investigate the air temperature response to land-use and land-cover change (LULCC; cropland expansion and deforestation) using subgrid land model output generated by a set of CMIP6 model simulations. Our study is motivated by the fact that ongoing land-use activities are occurring at local scales, typically significantly smaller than the resolvable scale of a grid cell in Earth system models. It aims to explore the potential for a multimodel approach to better characterize LULCC local climatic effects. On an annual scale, the CMIP6 models are in general agreement that croplands are warmer than primary and secondary land (psl; mainly forests, grasslands, and bare ground) in the tropics and cooler in the mid–high latitudes, except for one model. The transition from warming to cooling occurs at approximately 40°N. Although the surface heating potential, which combines albedo and latent heat flux effects, can explain reasonably well the zonal mean latitudinal subgrid temperature variations between crop and psl tiles in the historical simulations, it does not provide a good prediction on subgrid temperature for other land tile configurations (crop vs forest; grass vs forest) under Shared Socioeconomic Pathway 5–8.5 (SSP5–8.5) forcing scenarios. A subset of simulations with the CESM2 model reveals that latitudinal subgrid temperature variation is positively related to variation in net surface shortwave radiation and negatively related to variation in the surface energy redistribution factor, with a dominant role from the latter south of 30°N. We suggest that this emergent relationship can be used to benchmark the performance of land surface parameterizations and for prediction of local temperature response to LULCC.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2016

Application of the CERES Flux-by-Cloud Type Simulator to GCM Output

The CERES Flux By CloudType data product produces CERES top-of-atmosphere (TOA) fluxes by region and cloud type. Here, the cloud types are defined by cloud optical depth (t) and cloud top pressure (pc), with bins similar to those used by ISCCP (International Satellite Cloud Climatology Project). This data product has the potential to be a powerful tool for the evaluation of the clouds produced by climate models by helping to identify which physical parameterizations have problems (e.g., boundary-layer parameterizations, convective clouds, processes that affect surface albedo). Also, when the flux-by-cloud type and frequency of cloud types are simultaneously used to evaluate a model, the results can determine whether an unrealistically large or small occurrence of a given cloud type has an important radiative impact for a given region. A simulator of the flux-by-cloud type product has been applied to three-hourly data from the year 2008 from the UK Met Office HadGEM2-A model using the Langley Fu-Lour radiative transfer model to obtain TOA SW and LW fluxes.

Eitzen, Zachary↗
DOE OSTI2020

Day-ahead photovoltaic power production forecasting methodology based on machine learning and statistical post-processing

A main challenge towards ensuring large-scale and seamless integration of photovoltaic systems is to improve the accuracy of energy yield forecasts, especially in grid areas of high photovoltaic shares. The scope of this paper is to address this issue by presenting a unified methodology for hourly-averaged day-ahead photovoltaic power forecasts with improved accuracy, based on data-driven machine learning techniques and statistical post-processing. More specifically, the proposed forecasting methodology framework comprised of a data quality stage, data-driven power output machine learning model development (artificial neural networks), weather clustering assessment (K-means clustering), post-processing output optimisation (linear regressive correction method) and the final performance accuracy evaluation. The results showed that the application of linear regression coefficients to the forecasted outputs of the developed day-ahead photovoltaic power production neural network improved the performance accuracy by further correcting solar irradiance forecasting biases. The resulting optimised model provided a mean absolute percentage error of 4.7% when applied to historical system datasets. Finally, the model was validated both, at a hot as well as a cold semi-arid climatic location, and the obtained results demonstrated close agreement by yielding forecasting accuracies of mean absolute percentage error of 4.7% and 6.3%, respectively. Finally, the validation analysis provides evidence that the proposed model exhibits high performance in both forecasting accuracy and stability.

42 ENGINEERING↗
DOE OSTI2019

Increasing impacts of extreme droughts on vegetation productivity under climate change

Terrestrial gross primary production (GPP) is the basis of food production and vegetation growth globally1 and plays a critical role in regulating atmospheric CO2 through its impact on ecosystem carbon balance. Even though higher CO2 concentrations in future decades can increase GPP2, low soil water availability, heat stress and disturbances associated with droughts could reduce the benefits of such CO2 fertilization. Here we analysed outputs of 13 Earth system models to show an increasingly stronger impact on GPP by extreme droughts than by mild and moderate droughts over the twenty-first century. Due to a dramatic increase in the frequency of extreme droughts, the magnitude of globally averaged reductions in GPP associated with extreme droughts was projected to be nearly tripled by the last quarter of this century (2075–2099) relative to that of the historical period (1850–1999) under both high and intermediate GHG emission scenarios. By contrast, the magnitude of GPP reduction associated with mild and moderate droughts was not projected to increase substantially. Our analysis indicates a high risk of extreme droughts to the global carbon cycle with atmospheric warming; however, this risk can be potentially mitigated by positive anomalies of GPP associated with favourable environmental conditions.

Xu, Chonggang↗
NASA NTRS2017

El Nino-Induced Tropical Ocean/Land Energy Exchange in MERRA-2 and M2AMIP

Studies have shown the correlation and connection of surface temperatures across the globe, ocean and land, related to Tropical SSTs especially El Nino. This climate variability greatly influences regional weather and hydroclimate extremes (e.g. drought and flood). In this paper, we evaluate the relationship of temperatures across the tropical oceans and continents in MERRA-2, and also in a newly developed MERRA-2 AMIP ensemble simulation (M2AMIP). M2AMIP uses the same model and spatial resolution as MERRA-2, producing the same output diagnostics over 10 ensemble members. Composite El Nino temperature data are compared with observations to evaluate the land/sea contrast, variations and phase relationship. The temperature variations are related to surface heat fluxes and the atmospheric temperatures and transport, to identify the processes that lead to the lagged redistribution of heat in the tropics and beyond. Discernable cloud, radiation and data assimilation changes accompany the onset of El Nino affecting continental regions through the progression to and following the peak values. While the model represents these variations in general, regional strengths and weaknesses can be identified.

Bosilovich, Michael G.↗
NASA NTRS2019

Increasing Impacts of Extreme Droughts on Vegetation Productivity Under Climate Change

Terrestrial gross primary production (GPP) is the basis of vegetation growth and food production globally and plays a critical role in regulating atmospheric CO2 through its impact on ecosystem carbon balance. Even though higher CO2 concentrations in future decades can increase GPP, low soil water availability, heat stress and disturbances associated with droughts could reduce the benefits of such CO2 fertilization. Here we analysed outputs of 13 Earth system models to show an increasingly stronger impact on GPP by extreme droughts than by mild and moderate droughts over the twenty-first century. Due to a dramatic increase in the frequency of extreme droughts, the magnitude of globally averaged reductions in GPP associated with extreme droughts was projected to be nearly tripled by the last quarter of this century (2075-2099) relative to that of the historical period (1850-1999) under both high and intermediate GHG emission scenarios. By contrast, the magnitude of GPP reductions associated with mild and moderate droughts was not projected to increase substantially. Our analysis indicates a high risk of extreme droughts to the global carbon cycle with atmospheric warming; however, this risk can be potentially mitigated by positive anomalies of GPP associated with favourable environmental conditions.

Xu, Chonggang↗
NASA NTRS2023

Changing Characteristics of Tropical Extreme Precipitation-Cloud Regimes in Warmer Climates

In this paper, we have investigated the changing characteristics of climatic scale (monthly) tropical extreme precipitation in warming climates using the Energy Exascale Earth System Model (E3SM). Results are from Atmospheric Model Intercomparison Project (AMIP)-type simulations driven by a) control experiment with present-day sea surface temperature (SST) and CO 2 concentration, b) P4K, same as in a) but with uniform increase by 4K in SST globally, and c) same as in a), but with imposed SST and CO 2 concentration from outputs of coupled E3SM forced by 4xCO 2 concentration. We find that as the surface warms under P4K and 4xCO 2 , both convective and stratiform rain increase. Importantly, there is increasing fractional contribution of stratiform rain as a function of precipitation intensity, with the most extreme but rare events occurring preferentially over land more than ocean, more so under 4xCO 2 than P4K. Extreme precipitation is facilitated by increased precipitation efficiency, reflecting accelerated rates of recycling of precipitation-cloud water (both liquid and ice phases) in regions with colder anvil cloud tops. Changes in vertical profiles of clouds, condensation heating and vertical motions indicate increasing precipitation-cloud-circulation organization from control, P4K to 4xCO 2 . Results suggest that large-scale ocean warming, i.e., P4K, is the primary cause contributing to organization structure resembling the well-known Mesoscale Convective System (MCS) for increased extreme precipitation on shorter (hourly to daily) time scales. Additional 4xCO 2 atmospheric radiative heating, and dynamically consistent anomalous SST further amplifies the MCS organization under P4K. Analyses of the surface moist static energy distribution show that increase in surface moisture (temperature) under P4K and 4xCO 2 is the key driver leading to enhanced convective instability over tropical ocean (land). However, the fast and large increase in land surface temperature and lack of available local moisture result in strong reduction in land surface relative humidity, reflecting severe drying and enhanced convective inhibition (CIN). It is argued that very extreme and rare “record-breaking” precipitation events found over land under P4K, more so under 4xCO 2 , are likely due to the delayed onset of deep convection, i.e., the longer the suppression of deep convection by CIN, the more severe the extreme precipitation when eventually occur, due to the release of large amount of stored surplus convective available potential energy in the lower troposphere during prolonged CIN.

climate-scale extreme tropical precipitation↗
DOE OSTI2021

plutus: An R package to calculate electricity investments and stranded assets from the Global Change Analysis Model (GCAM)

plutus is an R package that post-processes outputs from the Global Change Analysis Model (GCAM) to calculate the power sector capital investment costs and stranded asset costs associated with GCAM projections of electricity generation by technology. GCAM is a market equilibrium model used to examine the dynamics of the coupled human-Earth system and the response of this system to global change, including socioeconomics, technology, climate, and policy. GCAM tracks electricity generation by technology and vintage over 32 geopolitical regions throughout the lifetime of each generating technology. plutus extends GCAM functionality by (1) estimating the foregone economic value of prematurely retired power plants as a result of economic-induced retirements compared to scheduled lifetimes; and (2) estimating the new installations and capital investments driven by future changes in economic, energy, agriculture, and land-use systems in GCAM.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2016

Estimating Model Prediction Error: Should You Treat Predictions as Fixed or Random?

Crop models are important tools for impact assessment of climate change, as well as for exploring management options under current climate. It is essential to evaluate the uncertainty associated with predictions of these models. We compare two criteria of prediction error; MSEP fixed, which evaluates mean squared error of prediction for a model with fixed structure, parameters and inputs, and MSEP uncertain( X), which evaluates mean squared error averaged over the distributions of model structure, inputs and parameters. Comparison of model outputs with data can be used to estimate the former. The latter has a squared bias term, which can be estimated using hindcasts, and a model variance term, which can be estimated from a simulation experiment. The separate contributions to MSEP uncertain (X) can be estimated using a random effects ANOVA. It is argued that MSEP uncertain (X) is the more informative uncertainty criterion, because it is specific to each prediction situation.

model structure uncertainty↗
NASA NTRS2021

Tonlé Sap Food Security and Agriculture II: Evaluating Changes in Ecosystem Vitality and Freshwater Healthin the Tonlé Sap Basin using Remotely Sensed Data

The Tonlé Sap Lake and river basin in central Cambodia provide critical ecosystem services to the region, including fisheries, agricultural irrigation, hydropower, and biodiverse habitats. Deforestation, increased pumping for farming, and effects of climate change such as droughts and forest fires threaten the health of the lake and food security in the region. This project built upon the previous term through a partnership with Conservation International (CI), the Cambodian Ministry of Water Resources and Meteorology, and the Tonlé Sap Authority to assess ecosystem vitality and implement CI’s Freshwater Health Index (FHI) tool, in an effort to prioritize resource expenditure and highlight areas of concern. Due to the COVID-19 pandemic and related travel restrictions, partners had not been able to readily collect in situ data for the past year, which make up the majority of FHI inputs. To help fill this data gap, we developed a methodology for using Gravity Recovery and Climate Experiment (GRACE) satellite data to calculate groundwater storage depletion, and a Python Application Programming Interface for processing and formatting remotely-sensed data for the Soil and Water Assessment Tool (SWAT) model. We then used SWAT to model nutrient flows and of phosphorous, nitrogen and suspended sediments amounts in the basin from October 2000 to December 2020. These outputs served as inputs for the FHI and provided policy makers with robust monitoring information to aid decision-making in the area and safeguard the lake’s vital fisheries and biodiversity.

Justine Spore↗
DOE OSTI2021

A data-driven global soil heterotrophic respiration dataset and the drivers of its inter-annual variability

Soil heterotrophic respiration (SHR), one of the primary carbon fluxes from terrestrial ecosystems to the atmosphere, is important for carbon-climate feedbacks because of its sensitivity to available litter and soil carbon, climatic conditions, and nutrient availability. However, until recently limited SHR data were available, and most published global SHR estimates have either a short time span, coarse spatial resolution, or reply on overly-simple model formulations. To better understand and quantify the global distribution of SHR and its sensitivity to climate variability, we produced a new global SHR dataset using Random Forest algorithms, up-scaling 455 point data from the Global Soil Respiration Database (SRDB 4.0) with gridded fields of climatic, edaphic and productivity as explanatory variables. We estimated a global total SHR of 46.8 Pg C yr-1 over 1985-2013 (95% confidence interval: 38.6-56.3 Pg C yr-1), with a significant increasing trend of 0.03 Pg C yr-2 during this period. We found that the choice of soil moisture datasets contributes more to the difference among these data-driven SHR members rather than that of productivity, temperature and precipitation data sources. We also analyzed the influence of climatic variables on the inter-annual variability (IAV) of our SHR product. Water availability was the dominant driver of IAV at global scales, although the inferred sensitivity depends on the choice of the soil moisture gridded dataset. At the ecosystem scale, temperature strongly controls the IAV of SHR in tropical forests, while water availability dominates in extra-tropical forest and semi-arid regions. Our machine-learning gridded SHR dataset and outputs from process-based land surface models (TRENDYv6) show agreement for a strong association between water variability and SHR IAV at the global scale, but the two approaches lead to different temporal trend globally and different controlling variables for IAV at the ecosystem scale. Our study provides evidence for the pervasive and important role of water availability in driving SHR, indicating both a direct effect limiting decomposition rates and an indirect effect through the amount of fresh organic matter made available to SHR from productivity. In consideration of potential limitations and uncertainties remaining in our data-driven SHR datasets, we call for a more scientifically designed observation network for SHR, more observation data compilation, and increased use of deep learning methods making maximum use of observation data in hand. This will benefit process-based models, and improve our understanding of SHR response to future anomalous environmental conditions.

Yao, Yitong↗