The signature of internal variability in the terrestrial carbon cycle

Uncertainty in model initial states produces uncertainty in climate simulations because of unforced variability internal to the climate system. Climate scientists use initial-condition ensembles to separate the forced signal of climate change from the unforced internal variability. Our analysis of an 11-member initial-condition ensemble from the Community Earth System Model Version 2 that spans the period 1850&ndash;2014 shows that a similar ensemble approach is needed to robustly assess trends in the terrestrial carbon cycle. Uncertainty in model initialization gives rise to internal variability that masks trends in carbon fluxes, and also creates spurious unforced trends, during the period 1960&ndash;2014 across North America, meaning that a single model realization can diverge from the observational record or from other models simply because of random behavior. The forced response is, however, evident in the ensemble mean and emerges from the noise of unforced variability at decadal timescales. Our results suggest that trends in the observational record must be interpreted with caution because of multiple possible histories that would have been observed if the sequence of internal variability had unfolded differently. Furthermore, internal variability produces irreducible uncertainty in the carbon cycle, leading to ambiguity in the magnitude and sign of carbon cycle trends, especially at small spatial scales and short timescales. The small spread in initial land carbon pools at 1850 suggests that internal climate variability arising from atmospheric and oceanic initialization, not the biogeochemical initialization, is the predominant cause of carbon cycle variability among ensemble members. Initial-condition ensembles with other Earth system models are needed to develop a multi-model understanding of internal variability in the terrestrial carbon cycle. &nbsp;</p

Implementation and Evaluation of Irrigation Techniques in the Community Land Model

Several previous studies have highlighted the irrigation-induced impacts on the global and regional water cycle, energy budget, and near-surface climate. While land models are widely used to address this question, the implementations of irrigation in these models vary in complexity. Here, we expand the representation of irrigation in Community Land Model to enable six different irrigation methods. We find that using a combination of irrigation methods, including default, sprinkler, flood and paddy techniques performs best as determined by evaluating the simulated irrigation water withdrawals against observations, and therefore select this combination as the new irrigation scheme. Then, the impact of the new irrigation scheme on surface fluxes is evaluated and detected using single-point simulations. Finally, the global and regional irrigation-induced impacts on surface energy and water fluxes are compared using both the original and the new irrigation scheme. The new irrigation scheme substantially reduces the bias and root-mean-square error of simulated irrigation water withdrawal in the USA and other countries, but considerably overestimates withdrawals in Central China. Results of single-point experiments show that different irrigation methods have different effects on surface fluxes, while the magnitudes are small. At the global scale, the new scheme enlarges the irrigation-induced impacts on water and energy variables relative to the original scheme, with varying magnitudes across regions. Overall, our results suggest that this newly developed scheme is a better tool for simulating irrigation-induced impacts on climate, and highlight the added value of incorporating human water management in Earth system models

Behavioral responses of the endemic shrimp Halocardina rubra (Malacostraca:Atyidae) to an introduced fish, Gambusia affinis (Actinopterygii: Poeciliidae) and implications for the trophic structure of Hawaiian anchialine ponds

In the Hawaiian Islands, intentionally introduced exotic fishes have been linked to changes in native biodiversity and community composition. In 1905, the mosquito fish Gambusia affinis was introduced to control mosquitoes. Subsequently, G. affinis spread throughout the Islands and into coastal anchialine ponds. Previous studies suggest that presence of invasive fishes in anchialine ponds may eliminate native species, including the endemic shrimp Halocaridina rubra. We examined effects of G. affinis on H. rubra populations in anchialine ponds on the Kona-Kohala coast of the island of Hawai/i. In the presence of G. affinis, H. rubra exhibited a diel activity pattern that was not seen in fishless ponds. Shrimp in ponds with fish were active only at night. This pattern was evident in anchialine ponds and in laboratory experiments. In laboratory predation experiments, G. affinis preferentially consumed smaller H. rubra, and in the field the H. rubra collected from invaded sites were larger than those from fishless ponds. Analysis of trophic position using stable isotope analyses showed that feeding of H. rubra was not significantly distinct from that of snails, assumed to feed at trophic level 2.0 on epilithic algae, but G. affinis was slightly omnivorous, feeding at tropic level 2.2. The mosquito fish diet was apparently composed primarily of algae when the defensive behavior of H. rubra made them substantially unavailable as prey. The effect of successful establishment of G. affinis on shrimp behavior has the potential to alter abundance of benthic algae and processing and recycling of nutrients in anchialine pond ecosystems

High predictability of terrestrial carbon fluxes from an initialized decadal prediction system

Interannual variations in the flux of carbon dioxide (CO2) between the land surface and the atmosphere are the dominant component of interannual variations in the atmospheric CO2 growth rate. Here, we investigate the potential to predict variations in these terrestrial carbon fluxes 1&ndash;10 years in advance using a novel set of retrospective decadal forecasts of an Earth system model. We demonstrate that globally-integrated net ecosystem production (NEP) exhibits high potential predictability for 2 years following forecast initialization. This predictability exceeds that from a persistence or uninitialized forecast conducted with the same Earth system model. The potential predictability in NEP derives mainly from high predictability in ecosystem respiration, which itself is driven by vegetation carbon and soil moisture initialization. Our findings unlock the potential to forecast the terrestrial ecosystem in a changing environment. </div

Impacts of shifting phenology on boundary layer dynamics in North America in the CESM

Published versio

Disentangling land model uncertainty via Matrix-based Ensemble Model Inter-comparison Platform (MEMIP)

Background Large uncertainty in modeling land carbon (C) uptake heavily impedes the accurate prediction of the global C budget. Identifying the uncertainty sources among models is crucial for model improvement yet has been difficult due to multiple feedbacks within Earth System Models (ESMs). Here we present a Matrix-based Ensemble Model Inter-comparison Platform (MEMIP) under a unified model traceability framework to evaluate multiple soil organic carbon (SOC) models. Using the MEMIP, we analyzed how the vertically resolved soil biogeochemistry structure influences SOC prediction in two soil organic matter (SOM) models. By comparing the model outputs from the C-only and CN modes, the SOC differences contributed by individual processes and N feedback between vegetation and soil were explicitly disentangled. Results Results showed that the multi-layer models with a vertically resolved structure predicted significantly higher SOC than the single layer models over the historical simulation (1900–2000). The SOC difference between the multi-layer models was remarkably higher than between the single-layer models. Traceability analysis indicated that over 80% of the SOC increase in the multi-layer models was contributed by the incorporation of depth-related processes, while SOC differences were similarly contributed by the processes and N feedback between models with the same soil depth representation. Conclusions The output suggested that feedback is a non-negligible contributor to the inter-model difference of SOC prediction, especially between models with similar process representation. Further analysis with TRENDY v7 and more extensive MEMIP outputs illustrated the potential important role of multi-layer structure to enlarge the current ensemble spread and the necessity of more detail model decomposition to fully disentangle inter-model differences. We stressed the importance of analyzing ensemble outputs from the fundamental model structures, and holding a holistic view in understanding the ensemble uncertainty

Evaluation of global terrestrial evapotranspiration using state-of-the-art approaches in remote sensing, machine learning and land surface modeling

Evapotranspiration (ET) is critical in linking global water, carbon and energy cycles. However, direct measurement of global terrestrial ET is not feasible. Here, we first reviewed the basic theory and state-of-the-art approaches for estimating global terrestrial ET, including remote-sensing-based physical models, machine-learning algorithms and land surface models (LSMs). We then utilized 4 remote-sensing-based physical models, 2 machine-learning algorithms and 14 LSMs to analyze the spatial and temporal variations in global terrestrial ET. The results showed that the ensemble means of annual global terrestrial ET estimated by these three categories of approaches agreed well, with values ranging from 589.6 mm yr−1 (6.56×104 km3 yr−1) to 617.1 mm yr−1 (6.87×104 km3 yr−1). For the period from 1982 to 2011, both the ensembles of remote-sensing-based physical models and machine-learning algorithms suggested increasing trends in global terrestrial ET (0.62 mm yr−2 with a significance level of p0.05), although many of the individual LSMs reproduced an increasing trend. Nevertheless, all 20 models used in this study showed that anthropogenic Earth greening had a positive role in increasing terrestrial ET. The concurrent small interannual variability, i.e., relative stability, found in all estimates of global terrestrial ET, suggests that a potential planetary boundary exists in regulating global terrestrial ET, with the value of this boundary being around 600 mm yr−1. Uncertainties among approaches were identified in specific regions, particularly in the Amazon Basin and arid/semiarid regions. Improvements in parameterizing water stress and canopy dynamics, the utilization of new available satellite retrievals and deep-learning methods, and model–data fusion will advance our predictive understanding of global terrestrial ET