Huimin Lei

Global patterns and controls of soil organic carbon dynamics as simulated by multiple terrestrial biosphere models: Current status and future directions

Abstract Soil is the largest organic carbon (C) pool of terrestrial ecosystems, and C loss from soil accounts for a large proportion of land‐atmosphere C exchange. Therefore, a small change in soil organic C (SOC) can affect atmospheric carbon dioxide (CO2) concentration and climate change. In the past decades, a wide variety of studies have been conducted to quantify global SOC stocks and soil C exchange with the atmosphere through site measurements, inventories, and empirical/process‐based modeling. However, these estimates are highly uncertain, and identifying major driving forces controlling soil C dynamics remains a key research challenge. This study has compiled century‐long (1901–2010) estimates of SOC storage and heterotrophic respiration (Rh) from 10 terrestrial biosphere models (TBMs) in the Multi‐scale Synthesis and Terrestrial Model Intercomparison Project and two observation‐based data sets. The 10 TBM ensemble shows that global SOC estimate ranges from 425 to 2111 Pg C (1 Pg = 1015 g) with a median value of 1158 Pg C in 2010. The models estimate a broad range of Rh from 35 to 69 Pg C yr−1 with a median value of 51 Pg C yr−1 during 2001–2010. The largest uncertainty in SOC stocks exists in the 40–65°N latitude whereas the largest cross‐model divergence in Rh are in the tropics. The modeled SOC change during 1901–2010 ranges from −70 Pg C to 86 Pg C, but in some models the SOC change has a different sign from the change of total C stock, implying very different contribution of vegetation and soil pools in determining the terrestrial C budget among models. The model ensemble‐estimated mean residence time of SOC shows a reduction of 3.4 years over the past century, which accelerate C cycling through the land biosphere. All the models agreed that climate and land use changes decreased SOC stocks, while elevated atmospheric CO2 and nitrogen deposition over intact ecosystems increased SOC stocks—even though the responses varied significantly among models. Model representations of temperature and moisture sensitivity, nutrient limitation, and land use partially explain the divergent estimates of global SOC stocks and soil C fluxes in this study. In addition, a major source of systematic error in model estimations relates to nonmodeled SOC storage in wetlands and peatlands, as well as to old C storage in deep soil layers.

Disentangling climatic and anthropogenic controls on global terrestrial evapotranspiration trends

We examined natural and anthropogenic controls on terrestrial evapotranspiration (ET) changes from 1982 to 2010 using multiple estimates from remote sensing-based datasets and process-oriented land surface models. A significant increasing trend of ET in each hemisphere was consistently revealed by observationally-constrained data and multi-model ensembles that considered historic natural and anthropogenic drivers. The climate impacts were simulated to determine the spatiotemporal variations in ET. Globally, rising CO2 ranked second in these models after the predominant climatic influences, and yielded decreasing trends in canopy transpiration and ET, especially for tropical forests and high-latitude shrub land. Increasing nitrogen deposition slightly amplified global ET via enhanced plant growth. Land-use-induced ET responses, albeit with substantial uncertainties across the factorial analysis, were minor globally, but pronounced locally, particularly over regions with intensive land-cover changes. Our study highlights the importance of employing multi-stream ET and ET-component estimates to quantify the strengthening anthropogenic fingerprint in the global hydrologic cycle.

Impact of large‐scale climate extremes on biospheric carbon fluxes: An intercomparison based on MsTMIP data

Understanding the role of climate extremes and their impact on the carbon (C) cycle is increasingly a focus of Earth system science. Climate extremes such as droughts, heat waves, or heavy precipitation events can cause substantial changes in terrestrial C fluxes. On the other hand, extreme changes in C fluxes are often, but not always, driven by extreme climate conditions. Here we present an analysis of how extremes in temperature and precipitation, and extreme changes in terrestrial C fluxes are related to each other in 10 state-of-the-art terrestrial carbon models, all driven by the same climate forcing. We use model outputs from the North American Carbon Program Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP). A global-scale analysis shows that both droughts and heat waves translate into anomalous net releases of CO2 from the land surface via different mechanisms: Droughts largely decrease gross primary production (GPP) and to a lower extent total respiration (TR), while heat waves slightly decrease GPP but increase TR. Cold and wet periods have a smaller opposite effect. Analyzing extremes in C fluxes reveals that extreme changes in GPP and TR are often caused by strong shifts in water availability, but for extremes in TR shifts in temperature are also important. Extremes in net CO2 exchange are equally strongly driven by deviations in temperature and precipitation. Models mostly agree on the sign of the C flux response to climate extremes, but model spread is large. In tropical forests, C cycle extremes are driven by water availability, whereas in boreal forests temperature plays a more important role. Models are particularly uncertain about the C flux response to extreme heat in boreal forests.

A distributed scheme developed for eco-hydrological modeling in the upper Heihe River

Modeling the hydrological processes at catchment scale requires a flexible distributed scheme to represent the catchment topography, river network and vegetation pattern. This study has developed a distributed scheme for eco-hydrological simulation in the upper Heihe River. Based on a 1 km × 1 km grid system, the study catchment is divided into 461 sub-catchments, whose main streams form the streamflow pathway. Furthermore, a 1 km grid is represented by a number of topographically similar “hillslope-valley” systems, and the hillslope is the basic unit of the eco-hydrological simulation. This model is tested with a simplified hydrological simulation focusing on soil-water dynamics and streamflow routing. Based on a 12-year simulation from 2001 to 2012, it is found that variability in hydrological behavior is closely associated with climatic and landscape conditions especially vegetation types. The subsurface and groundwater flows dominate the total river runoff. This implies that the soil freezing and thawing process would significantly influence the runoff generation in the upper Heihe basin. Furthermore, the runoff components and water balance characteristics vary among different vegetation types, showing the importance of coupling the vegetation pattern into catchment hydrological simulation. This paper also discusses the model improvement to be done in future study.

Toward “optimal” integration of terrestrial biosphere models

Multimodel ensembles (MME) are commonplace in Earth system modeling. Here we perform MME integration using a 10-member ensemble of terrestrial biosphere models (TBMs) from the Multiscale synthesis and Terrestrial Model Intercomparison Project (MsTMIP). We contrast optimal (skill based for present-day carbon cycling) versus naive (“one model-one vote”) integration. MsTMIP optimal and naive mean land sink strength estimates (−1.16 versus −1.15 Pg C per annum respectively) are statistically indistinguishable. This holds also for grid cell values and extends to gross uptake, biomass, and net ecosystem productivity. TBM skill is similarly indistinguishable. The added complexity of skill-based integration does not materially change MME values. This suggests that carbon metabolism has predictability limits and/or that all models and references are misspecified. Resolving this issue requires addressing specific uncertainty types (initial conditions, structure, and references) and a change in model development paradigms currently dominant in the TBM community.

Estimation of evapotranspiration using a remote sensing model over agricultural land in the North China Plain

Evapotranspiration (ET) from agricultural land is the most common form of water consumption in the North China Plain. Based on in situ measurements by an eddy covariance (EC) system, the present study tested the surface energy balance system (SEBS) model for estimating surface energy fluxes from a wheat/maize rotation cropland. Applicability of the SEBS model using Moderate Resolution Imaging Spectroradiometer (MODIS) land products to estimate ET was further validated. The model proved to be appropriate to measure heat flux during the wheat-growing season, but underestimated the sensible heat flux during the maize-growing season. Although the SEBS model performed better during the wheat-growing season, the relative error of the latent heat flux was within 20% in both wheat- and maize-growing seasons. Results showed that the SEBS model using remote sensing data could provide a reasonable estimate of the surface energy fluxes from an irrigated cropland in the North China Plain.

Uncertainty in the response of terrestrial carbon sink to environmental drivers undermines carbon-climate feedback predictions

Terrestrial ecosystems play a vital role in regulating the accumulation of carbon (C) in the atmosphere. Understanding the factors controlling land C uptake is critical for reducing uncertainties in projections of future climate. The relative importance of changing climate, rising atmospheric CO2, and other factors, however, remains unclear despite decades of research. Here, we use an ensemble of land models to show that models disagree on the primary driver of cumulative C uptake for 85% of vegetated land area. Disagreement is largest in model sensitivity to rising atmospheric CO2 which shows almost twice the variability in cumulative land uptake since 1901 (1 s.d. of 212.8 PgC vs. 138.5 PgC, respectively). We find that variability in CO2 and temperature sensitivity is attributable, in part, to their compensatory effects on C uptake, whereby comparable estimates of C uptake can arise by invoking different sensitivities to key environmental conditions. Conversely, divergent estimates of C uptake can occur despite being based on the same environmental sensitivities. Together, these findings imply an important limitation to the predictability of C cycling and climate under unprecedented environmental conditions. We suggest that the carbon modeling community prioritize a probabilistic multi-model approach to generate more robust C cycle projections.

Responses of Natural Vegetation Dynamics to Climate Drivers in China from 1982 to 2011

This study investigated the spatiotemporal variation of vegetation growth and the influence of climatic drivers from 1982 to 2011 across China using datasets from the Normalized Difference Vegetation Index (NDVI) and climatic drivers. Long term trends, significance and abrupt change points of interannual NDVI time series were analyzed. We applied both simple regression and multi-regression models to quantify the effects of climatic drivers on vegetation growth and compare their relative contributions. Results show that on average, the growing season NDVI significantly increased by 0.0007 year-1, with 76.5% of the research area showed increasing NDVI from 1982 to 2011. Seasonally, NDVI increased at high rates during the spring and autumn while changed slightly during the summer. At a national scale, the growing season NDVI was significantly and positively correlated to temperature and precipitation, with temperature being the dominant factor. At regional scales, the growing season NDVI was dominated by increasing temperature in most forest-covered areas in southern China and dominated by precipitation in most grassland in northern China. Within the past three decades, the increasing trend of national mean NDVI abruptly changed in 1994, slowing down from 0.0008 year-1 to 0.0003 year-1. To be regional specific, the growing season NDVI in forest covered southern China has accelerated together with temperature since mid 1990s, while parts of the grassland in northern China have undergone stalled NDVI trends corresponding to slowed temperature increment and dropped precipitation since around 2000. Typical region analysis suggested that apart from long term trends and abrupt change points of climatic drivers, the processes of NDVI variation were also affected by other external factors such as drought and afforestation. Further studies are needed to investigate the nonlinear responses of vegetation growth to climatic drivers and effects of non-climate factors on vegetation growth.

Combining the Crop Coefficient of Winter Wheat and Summer Maize with a Remotely Sensed Vegetation Index for Estimating Evapotranspiration in the North China Plain

AbstractEstimates of actual evapotranspiration (ETc) in the wheat and maize fields are essential in effective planning of irrigation water use in the North China Plain. A widely used method for ETc estimation in agriculture is crucially dependent on the determination of crop coefficient curves. Estimating the coefficient coefficients from vegetation index (VI) is useful for regional ETc simulation because the VI can represent the actual crop conditions and capture the spatial variability. In this study, the basal crop coefficient and soil evaporation coefficient were combined with the commonly used VI obtained from satellite sensor based on the observed data from a flux tower. The basal crop coefficient had a fairly good linear relationship with the VI, and the soil evaporation coefficient was well related to the vegetation fraction, which was calculated from the VI. Using the VI-derived crop coefficient curves, ETc can be well simulated by a widely used ETc estimation method. Moreover, simulation of ETc ...

Decadal trends in the seasonal-cycle amplitude of terrestrial CO2 exchange resulting from the ensemble of terrestrial biosphere models

The seasonal-cycle amplitude (SCA) of the atmosphere–ecosystem carbon dioxide (CO2) exchange rate is a useful metric of the responsiveness of the terrestrial biosphere to environmental variations. It is unclear, however, what underlying mechanisms are responsible for the observed increasing trend of SCA in atmospheric CO2 concentration. Using output data from the Multi-scale Terrestrial Model Intercomparison Project (MsTMIP), we investigated how well the SCA of atmosphere–ecosystem CO2 exchange was simulated with 15 contemporary terrestrial ecosystem models during the period 1901–2010. Also, we made attempt to evaluate the contributions of potential mechanisms such as atmospheric CO2, climate, land-use, and nitrogen deposition, through factorial experiments using different combinations of forcing data. Under contemporary conditions, the simulated global-scale SCA of the cumulative net ecosystem carbon flux of most models was comparable in magnitude with the SCA of atmospheric CO2 concentrations. Results from factorial simulation experiments showed that elevated atmospheric CO2 exerted a strong influence on the seasonality amplification. When the model considered not only climate change but also land-use and atmospheric CO2 changes, the majority of the models showed amplification trends of the SCAs of photosynthesis, respiration, and net ecosystem production (+0.19 % to +0.50 % yr−1). In the case of land-use change, it was difficult to separate the contribution of agricultural management to SCA because of inadequacies in both the data and models. The simulated amplification of SCA was approximately consistent with the observational evidence of the SCA in atmospheric CO2 concentrations. Large inter-model differences remained, however, in the simulated global tendencies and spatial patterns of CO2 exchanges. Further studies are required to identify a consistent explanation for the simulated and observed amplification trends, including their underlying mechanisms. Nevertheless, this study implied that monitoring of ecosystem seasonality would provide useful insights concerning ecosystem dynamics.