Linking Vegetation-Climate-Fire Relationships in Sub-Saharan Africa to Key Ecological Processes in Two Dynamic Global Vegetation Models

Africa is largely influenced by fires, which play an important ecological role influencing the distribution and structure of grassland, savanna and forest biomes. Here vegetation strongly interacts with climate and other environmental factors, such as herbivory and humans. Fire-enabled Dynamic Global Vegetation Models (DGVMs) display high uncertainty in predicting the distribution of current tropical biomes and the associated transitions, mainly due to the way they represent the main ecological processes and feedbacks related to water and fire. The aim of this study is to evaluate the outcomes of two state-of-the–art DGVMs, LPJ-GUESS and JSBACH, also currently used in two Earth System Models (ESMs), in order to assess which key ecological processes need to be included or improved to represent realistic interactions between vegetation cover, precipitation and fires in sub-Saharan Africa. To this end, we compare models and remote-sensing data, analyzing the relationships between tree and grass cover, mean annual rainfall, average rainfall seasonality and average fire intervals, using generalized linear models, and we compare the patterns of grasslands, savannas, and forests in sub-Saharan Africa. Our analysis suggests that LPJ-GUESS (with a simple fire-model and complex vegetation description) performs well in regions of low precipitation, while in humid and mesic areas the representation of the fire process should probably be improved to obtain more open savannas. JSBACH (with a complex fire-model and a simple vegetation description) can simulate a vegetation-fire feedback that can maintain open savannas at intermediate and high precipitation, although this feedback seems to have stronger effects than observed, while at low precipitation JSBACH needs improvements in the representation of tree-grass competition and drought effects. This comparative process-based analysis permits to highlight the main factors that determine the tropical vegetation distribution in models and observations in sub-Saharan Africa, suggesting possible improvements in DGVMs and, consequently, in ESM simulations for future projections. Given the need to use carbon storage in vegetation as a climate mitigation measure, these models represent a valuable tool to improve our understanding of the sustainability of vegetation carbon pools as a carbon sink and the vulnerability to disturbances such as fire

Coupling carbon allocation with leaf and root phenology predicts tree-grass partitioning along a savanna rainfall gradient

© Author(s) 2016. The relative complexity of the mechanisms underlying savanna ecosystem dynamics, in comparison to other biomes such as temperate and tropical forests, challenges the representation of such dynamics in ecosystem and Earth system models. A realistic representation of processes governing carbon allocation and phenology for the two defining elements of savanna vegetation (namely trees and grasses) may be a key to understanding variations in tree-grass partitioning in time and space across the savanna biome worldwide. Here we present a new approach for modelling coupled phenology and carbon allocation, applied to competing tree and grass plant functional types. The approach accounts for a temporal shift between assimilation and growth, mediated by a labile carbohydrate store. This is combined with a method to maximize long-term net primary production (NPP) by optimally partitioning plant growth between fine roots and (leaves + stem). The computational efficiency of the analytic method used here allows it to be uniquely and readily applied at regional scale, as required, for example, within the framework of a global biogeochemical model. We demonstrate the approach by encoding it in a new simple carbon-water cycle model that we call HAVANA (Hydrology and Vegetation-dynamics Algorithm for Northern Australia), coupled to the existing POP (Population Orders Physiology) model for tree demography and disturbance-mediated heterogeneity. HAVANA-POP is calibrated using monthly remotely sensed fraction of absorbed photosynthetically active radiation (fPAR) and eddy-covariance-based estimates of carbon and water fluxes at five tower sites along the North Australian Tropical Transect (NATT), which is characterized by large gradients in rainfall and wildfire disturbance. The calibrated model replicates observed gradients of fPAR, tree leaf area index, basal area, and foliage projective cover along the NATT. The model behaviour emerges from complex feedbacks between the plant physiology and vegetation dynamics, mediated by shifting above- versus below-ground resources, and not from imposed hypotheses about the controls on tree-grass co-existence. Results support the hypothesis that resource limitation is a stronger determinant of tree cover than disturbance in Australian savannas

Coupling carbon allocation with leaf and root phenology predicts tree-grass partitioning along a savanna rainfall gradient

Abstract. The relative complexity of the mechanisms underlying savanna ecosystem dynamics, in comparison to other biomes such as temperate and tropical forests, challenges the representation of such dynamics in ecosystem and Earth system models. A realistic representation of processes governing carbon allocation and phenology for the two defining elements of savanna vegetation (namely trees and grasses) may be a key to understanding variations in tree/grass partitioning in time and space across the savanna biome worldwide. Here we present a new approach for modelling coupled phenology and carbon allocation, applied to competing tree and grass plant functional types. The approach accounts for a temporal shift between assimilation and growth, mediated by a labile carbohydrate store. This is combined with a method to maximise long-term net primary production (NPP) by optimally partitioning plant growth between fine roots and (leaves + stem). The computational efficiency of the analytic method used here allows it to be uniquely and readily applied at regional scale, as required, for example, within the framework of a global biogeochemical model. We demonstrate the approach by encoding it in a new simple carbon/water cycle model that we call HAVANA (Hydrology and Vegetation-dynamics Algorithm for Northern Australia), coupled to the existing POP (Population Orders Physiology) model for tree demography and disturbance-mediated heterogeneity. HAVANA-POP is calibrated using monthly remotely-sensed fraction of absorbed photosynthetically active radiation (fPAR) and eddy-covariance-based estimates of carbon and water fluxes at 5 tower sites along the Northern Australian Tropical Transect (NATT), which is characterized by large gradients in rainfall and wildfire disturbance. The calibrated model replicates observed gradients of fPAR, tree leaf area index, basal area and foliage projective cover along the NATT. The model behaviour emerges from complex feed-backs between the plant physiology and vegetation dynamics, mediated by shifting above- vs. below-ground resources, and not from imposed hypotheses about the controls on tree/grass co-existence. Results support the hypothesis that resource limitation is a stronger determinant of tree cover than disturbance in Australian savannas. </jats:p

Coupling carbon allocation with leaf and root phenology predicts tree-grass partitioning along a savanna rainfall gradient

The relative complexity of the mechanisms underlying savanna ecosystem dynamics, in comparison to other biomes such as temperate and tropical forests, challenges the representation of such dynamics in ecosystem and Earth system models. A realistic representation of processes governing carbon allocation and phenology for the two defining elements of savanna vegetation (namely trees and grasses) may be a key to understanding variations in tree–grass partitioning in time and space across the savanna biome worldwide. Here we present a new approach for modelling coupled phenology and carbon allocation, applied to competing tree and grass plant functional types. The approach accounts for a temporal shift between assimilation and growth, mediated by a labile carbohydrate store. This is combined with a method to maximize long-term net primary production (NPP) by optimally partitioning plant growth between fine roots and (leaves + stem). The computational efficiency of the analytic method used here allows it to be uniquely and readily applied at regional scale, as required, for example, within the framework of a global biogeochemical model. We demonstrate the approach by encoding it in a new simple carbon–water cycle model that we call HAVANA (Hydrology and Vegetation-dynamics Algorithm for Northern Australia), coupled to the existing POP (Population Orders Physiology) model for tree demography and disturbance-mediated heterogeneity. HAVANA-POP is calibrated using monthly remotely sensed fraction of absorbed photosynthetically active radiation (fPAR) and eddy-covariance-based estimates of carbon and water fluxes at five tower sites along the North Australian Tropical Transect (NATT), which is characterized by large gradients in rainfall and wildfire disturbance. The calibrated model replicates observed gradients of fPAR, tree leaf area index, basal area, and foliage projective cover along the NATT. The model behaviour emerges from complex feedbacks between the plant physiology and vegetation dynamics, mediated by shifting above- versus below-ground resources, and not from imposed hypotheses about the controls on tree–grass co-existence. Results support the hypothesis that resource limitation is a stronger determinant of tree cover than disturbance in Australian savannas.The contributions of V. Haverd and P. Briggs were made possible by the support of the Australian Climate Change Science Program. B. Smith acknowledges funding as an OCE Distinguished Visiting Scientist to the CSIRO Oceans & Atmosphere Flagship, Canberr

Evaluation of global fire simulations in CMIP6 Earth system models

Fire is the primary form of terrestrial ecosystem disturbance on a global scale and an important Earth system process. Most Earth system models (ESMs) have incorporated fire modeling, with 19 out of them submitting model outputs of fire-related variables to the Coupled Model Intercomparison Project Phase 6 (CMIP6). This study provides the first comprehensive evaluation of CMIP6 historical fire simulations by comparing them with multiple satellite-based products and charcoal-based historical reconstructions. Our results show that most CMIP6 models simulate the present-day global burned area and fire carbon emissions within the range of satellite-based products. They also capture the major features of observed spatial patterns and seasonal cycles, the relationship of fires with precipitation and population density, and the influence of El Ni&ntilde;o-Southern Oscillation (ENSO) on the interannual variability of tropical fires. Regional fire carbon emissions simulated by the CMIP6 models from 1850 to 2010 generally align with the charcoal-based reconstructions, although there are regional mismatches, such as in southern South America and eastern temperate North America prior to the 1910s and in temperate North America, eastern boreal North America, Europe, and boreal Asia since the 1980s. The CMIP6 simulations have addressed three critical issues identified in the CMIP5: (1) the simulated global burned area less than half of the observations, (2) the failure to reproduce the high burned area fraction observed in Africa, and (3) the weak fire seasonal variability. Furthermore, the CMIP6 models exhibit improved accuracy in capturing the observed relationship between fires and both climatic and socioeconomic drivers, and better align with the historical long-term trends indicated by charcoal-based reconstructions in most regions worldwide. However, the CMIP6 models still fail to reproduce the decline in global burned area and fire carbon emissions observed over the past two decades, mainly attributed to an underestimation of anthropogenic fire suppression, and the spring peak in fires in the Northern Hemisphere mid-latitudes, mainly due to an underestimation of crop fires. In addition, the model underestimates the fire sensitivity to wet-dry conditions, indicating the need to improve fuel wetness estimation. Based on these findings, we present specific guidance for fire scheme development and suggest the post-processing methodology for using CMIP6 multi-model outputs to generate reliable fire projection products

The Fire Modeling Intercomparison Project (FireMIP), phase 1: experimental and analytical protocols

The important role of fire in regulating vegetation community composition and contributions to emissions of greenhouse gases and aerosols make it a critical component of dynamic global vegetation models and Earth system models. Over two decades of development, a wide variety of model structures and mechanisms have been designed and incorporated into global fire models, which have been linked to different vegetation models. However, there has not yet been a systematic examination of how these different strategies contribute to model performance. Here we describe the structure of the first phase of the Fire Model Intercomparison Project (FireMIP), which for the first time seeks to systematically compare a number of models. By combining a standardized set of input data and model experiments with a rigorous comparison of model outputs to each other and to observations, we will improve the understanding of what drives vegetation fire, how it can best be simulated, and what new or improved observational data could allow better constraints on model behavior. Here we introduce the fire models used in the first phase of FireMIP, the simulation protocols applied, and the benchmarking system used to evaluate the models

Quantitative assessment of fire and vegetation properties in simulations with fire-enabled vegetation models from the Fire Model Intercomparison Project

Global fire-vegetation models are widely used to assess impacts of environmental change on fire regimes and the carbon cycle and to infer relationships between climate, land use and fire. However, differences in model structure and parameterizations, in both the vegetation and fire components of these models, could influence overall model performance, and to date there has been limited evaluation of how well different models represent various aspects of fire regimes. The Fire Model Intercomparison Project (FireMIP) is coordinating the evaluation of state-of-the-art global fire models, in order to improve projections of fire characteristics and fire impacts on ecosystems and human societies in the context of global environmental change. Here we perform a systematic evaluation of historical simulations made by nine FireMIP models to quantify their ability to reproduce a range of fire and vegetation benchmarks. The FireMIP models simulate a wide range in global annual total burnt area (39–536 Mha) and global annual fire carbon emission (0.91–4.75 Pg C yr−1) for modern conditions (2002–2012), but most of the range in burnt area is within observational uncertainty (345–468 Mha). Benchmarking scores indicate that seven out of nine FireMIP models are able to represent the spatial pattern in burnt area. The models also reproduce the seasonality in burnt area reasonably well but struggle to simulate fire season length and are largely unable to represent interannual variations in burnt area. However, models that represent cropland fires see improved simulation of fire seasonality in the Northern Hemisphere. The three FireMIP models which explicitly simulate individual fires are able to reproduce the spatial pattern in number of fires, but fire sizes are too small in key regions, and this results in an underestimation of burnt area. The correct representation of spatial and seasonal patterns in vegetation appears to correlate with a better representation of burnt area. The two older fire models included in the FireMIP ensemble (LPJ–GUESS–GlobFIRM, MC2) clearly perform less well globally than other models, but it is difficult to distinguish between the remaining ensemble members; some of these models are better at representing certain aspects of the fire regime; none clearly outperforms all other models across the full range of variables assessed

The EC-Earth3 Earth system model for the Coupled Model Intercomparison Project 6

The Earth system model EC-Earth3 for contributions to CMIP6 is documented here, with its flexible coupling framework, major model configurations, a methodology for ensuring the simulations are comparable across different high-performance computing (HPC) systems, and with the physical performance of base configurations over the historical period. The variety of possible configurations and sub-models reflects the broad interests in the EC-Earth community. EC-Earth3 key performance metrics demonstrate physical behavior and biases well within the frame known from recent CMIP models. With improved physical and dynamic features, new Earth system model (ESM) components, community tools, and largely improved physical performance compared to the CMIP5 version, EC-Earth3 represents a clear step forward for the only European community ESM. We demonstrate here that EC-Earth3 is suited for a range of tasks in CMIP6 and beyond