Regional Carbon Fluxes from Land Use and Land Cover Change in Asia, 1980-2009

We present a synthesis of the land-atmosphere carbon flux from land use and land cover change (LULCC) in Asia using multiple data sources and paying particular attention to deforestation and forest regrowth fluxes. The data sources are quasi-independent and include the U.N. Food and Agriculture Organization-Forest Resource Assessment (FAO-FRA 2015; country-level inventory estimates), the Emission Database for Global Atmospheric Research (EDGARv4.3), the 'Houghton' bookkeeping model that incorporates FAO-FRA data, an ensemble of 8 state-of-the-art Dynamic Global Vegetation Models (DGVM), and 2 recently published independent studies using primarily remote sensing techniques. The estimates are aggregated spatially to Southeast, East, and South Asia and temporally for three decades, 1980–1989, 1990-1999 and 2000-2009. Since 1980, net carbon emissions from LULCC in Asia were responsible for 20%-40% of global LULCC emissions, with emissions from Southeast Asia alone accounting for 15%-25% of global LULCC emissions during the same period. In the 2000s and for all Asia, three estimates (FAO-FRA, DGVM, Houghton) were in agreement of a net source of carbon to the atmosphere, with mean estimates ranging between 0.24 to 0.41 Pg C yr?1, whereas EDGARv4.3 suggested a net carbon sink of ?0.17 Pg C yr?1. Three of 4 estimates suggest that LULCC carbon emissions declined by at least 34% in the preceding decade (1990-2000). Spread in the estimates is due to the inclusion of different flux components and their treatments, showing the importance to include emissions from carbon rich peatlands and land management, such as shifting cultivation and wood harvesting, which appear to be consistently underreported

C4MIP - The Coupled Climate-Carbon Cycle Model Intercomparison Project: experimental protocol for CMIP6

Coordinated experimental design and implementation has become a cornerstone of global climate modelling. Model Intercomparison Projects (MIPs) enable systematic and robust analysis of results across many models, by reducing the influence of ad hoc differences in model set-up or experimental boundary conditions. As it enters its 6th phase, the Coupled Model Intercomparison Project (CMIP6) has grown significantly in scope with the design and documentation of individual simulations delegated to individual climate science communities. The Coupled Climate-Carbon Cycle Model Intercomparison Project (C4MIP) takes responsibility for design, documentation, and analysis of carbon cycle feedbacks and interactions in climate simulations. These feedbacks are potentially large and play a leading-order contribution in determining the atmospheric composition in response to human emissions of CO2 and in the setting of emissions targets to stabilize climate or avoid dangerous climate change. For over a decade, C4MIP has coordinated coupled climate-carbon cycle simulations, and in this paper we describe the C4MIP simulations that will be formally part of CMIP6. While the climate-carbon cycle community has created this experimental design, the simulations also fit within the wider CMIP activity, conform to some common standards including documentation and diagnostic requests, and are designed to complement the CMIP core experiments known as the Diagnostic, Evaluation and Characterization of Klima (DECK). C4MIP has three key strands of scientific motivation and the requested simulations are designed to satisfy their needs: (1) pre-industrial and historical simulations (formally part of the common set of CMIP6 experiments) to enable model evaluation, (2) idealized coupled and partially coupled simulations with 1% per year increases in CO2 to enable diagnosis of feedback strength and its components, (3) future scenario simulations to project how the Earth system will respond to anthropogenic activity over the 21st century and beyond. This paper documents in detail these simulations, explains their rationale and planned analysis, and describes how to set up and run the simulations. Particular attention is paid to boundary conditions, input data, and requested output diagnostics. It is important that modelling groups participating in C4MIP adhere as closely as possible to this experimental design

Impact of the numerical solution approach of a plant hydrodynamic model (v0.1) on vegetation dynamics

Numerous plant hydrodynamic models have started to be implemented in vegetation dynamics models, reflecting the central role of plant hydraulic traits in driving water, energy, and carbon cycles, as well as plant adaptation to climate change. Different numerical approximations of the governing equations of the hydrodynamic models have been documented, but the numerical accuracy of these models and its subsequent effects on the simulated vegetation function and dynamics have rarely been evaluated. Using different numerical solution methods (including implicit and explicit approaches) and vertical discrete grid resolutions, we evaluated the numerical performance of a plant hydrodynamic module in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES-HYDRO version 0.1) based on single-point and global simulations. Our simulation results showed that when near-surface vertical grid spacing is coarsened (grid size >10 cm), the model significantly overestimates aboveground biomass (AGB) in most of the temperate forest locations and underestimates AGB in the boreal forest locations, as compared to a simulation with finer vertical grid spacing. Grid coarsening has a small effect on AGB in the tropical zones of Asia and South America. In particular, coarse surface grid resolution should not be used when there are large and prolonged water content differences among soil layers at depths due to long dry-season duration and/or well-drained soil or when soil evaporation is a dominant fraction of evapotranspiration. Similarly, coarse surface grid resolution should not be used when there is lithologic discontinuity along the soil depth. This information is useful for uncertainty quantification, sensitivity analysis, or the training of surrogate models to design the simulations when computational cost limits the use of ensemble simulations

Tree crown damage and its effects on forest carbon cycling in a tropical forest.

Crown damage can account for over 23% of canopy biomass turnover in tropical forests and is a strong predictor of tree mortality; yet, it is not typically represented in vegetation models. We incorporate crown damage into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), to evaluate how lags between damage and tree recovery or death alter demographic rates and patterns of carbon turnover. We represent crown damage as a reduction in a tree's crown area and leaf and branch biomass, and allow associated variation in the ratio of aboveground to belowground plant tissue. We compare simulations with crown damage to simulations with equivalent instant increases in mortality and benchmark results against data from Barro Colorado Island (BCI), Panama. In FATES, crown damage causes decreases in growth rates that match observations from BCI. Crown damage leads to increases in carbon starvation mortality in FATES, but only in configurations with high root respiration and decreases in carbon storage following damage. Crown damage also alters competitive dynamics, as plant functional types that can recover from crown damage outcompete those that cannot. This is a first exploration of the trade-off between the additional complexity of the novel crown damage module and improved predictive capabilities. At BCI, a tropical forest that does not experience high levels of disturbance, both the crown damage simulations and simulations with equivalent increases in mortality does a reasonable job of capturing observations. The crown damage module provides functionality for exploring dynamics in forests with more extreme disturbances such as cyclones and for capturing the synergistic effects of disturbances that overlap in space and time

Functionally Assembled Terrestrial Ecosystem Simulator (FATES) for Hurricane Disturbance and Recovery

Abstract Tropical cyclones are an important cause of forest disturbance, and major storms caused severe structural damage and elevated tree mortality in coastal tropical forests. Model capabilities that can be used to understand post‐hurricane forest recovery are still limited. We use a vegetation demography model, the Functionally Assembled Terrestrial Ecosystem Simulator, coupled with the Energy Exascale Earth System Model Land Model (ELM‐FATES) to study the processes and the key factors regulating post‐hurricane forest recovery. We implemented hurricane‐induced forest damage, including defoliation, structural biomass reduction, and tree mortality, performed ensemble model simulations, and used random forest feature importance. For the simulation in the Luquillo Experimental Forest, Puerto Rico, we identified factors controlling the post‐hurricane forest recovery, and quantified the sensitivity of key model parameters to the post‐hurricane forest recovery. The results indicate a tendency for the Bisley forests to shift toward the light demanding plant functional type (PFT) when the pre‐hurricane biomass between the light demanding and shade tolerant PFTs is nearly equal and forests experience hurricane disturbance with mortality >60% for both the two PFTs. Under more realistic conditions where the shade tolerant PFT is initially dominant, mortality >80% is required for a shift toward dominance of the light demanding PFT at Bisley. Hurricane mortality and background mortality are the two major factors regulating post‐hurricane forest recovery in simulations. This research improves understanding of the ELM‐FATES model behavior associated with hurricane disturbance and provides guidance for dynamic vegetation model development in representing hurricane induced forest damage with varied intensities

A review of global wetland carbon stocks and management challenges

Wetlands have unique soil, vegetation, and biogeochemistry that arises from their landscape position and wetland hydrology, which creates low oxygen levels in the soil. With reduced oxygen availability, plants develop adaptations to survive, such as aerenchyma, that allow transport of atmospheric oxygen to their roots, and soil microbial communities become dominated by anaerobic respiration processes that are less efficient in oxidizing carbon. Combined, the above- and belowground carbon stocks of wetlands play a key role in the global carbon cycle at varying time scales. This chapter provides a comprehensive assessment of wetland carbon stocks, research methodologies, and their historical and future trajectories. We estimate wetland carbon stocks range between 520-710 PgC (and 1792 to 1882 PgC with permafrost carbon) globally.</p

Drivers and mechanisms of tree mortality in moist tropical forests.

Tree mortality rates appear to be increasing in moist tropical forests (MTFs) with significant carbon cycle consequences. Here, we review the state of knowledge regarding MTF tree mortality, create a conceptual framework with testable hypotheses regarding the drivers, mechanisms and interactions that may underlie increasing MTF mortality rates, and identify the next steps for improved understanding and reduced prediction. Increasing mortality rates are associated with rising temperature and vapor pressure deficit, liana abundance, drought, wind events, fire and, possibly, CO2 fertilization-induced increases in stand thinning or acceleration of trees reaching larger, more vulnerable heights. The majority of these mortality drivers may kill trees in part through carbon starvation and hydraulic failure. The relative importance of each driver is unknown. High species diversity may buffer MTFs against large-scale mortality events, but recent and expected trends in mortality drivers give reason for concern regarding increasing mortality within MTFs. Models of tropical tree mortality are advancing the representation of hydraulics, carbon and demography, but require more empirical knowledge regarding the most common drivers and their subsequent mechanisms. We outline critical datasets and model developments required to test hypotheses regarding the underlying causes of increasing MTF mortality rates, and improve prediction of future mortality under climate change

Toward more realistic projections of soil carbon dynamics by Earth system models

International audienceSoil carbon (C) is a critical component of Earth system models (ESMs), and its diverse representations are a major source of the large spread across models in the terrestrial C sink from the third to fifth assessment reports of the Intergovernmental Panel on Climate Change (IPCC). Improving soil C projections is of a high priority for Earth system modeling in the future IPCC and other assessments. To achieve this goal, we suggest that (1) model structures should reflect real-world processes, (2) parameters should be calibrated to match model outputs with observations, and (3) external forcing variables should accurately prescribe the environmental conditions that soils experience. First, most soil C cycle models simulate C input from litter production and C release through decomposition. The latter process has traditionally been represented by first-order decay functions, regulated primarily by temperature, moisture, litter quality, and soil texture. While this formulation well captures macroscopic soil organic C (SOC) dynamics, better understanding is needed of their underlying mechanisms as related to microbial processes, depth-dependent environmental controls, and other processes that strongly affect soil C dynamics. Second, incomplete use of observations in model parameterization is a major cause of bias in soil C projections from ESMs. Optimal parameter calibration with both pool-and flux-based data sets through data assimilation is LUO ET AL. SOIL CARBON MODELING 40 PUBLICATION