Modeling long-term fire impact on ecosystem characteristics and surface energy using a process-based vegetation–fire model SSiB4/TRIFFID-Fire v1.0
Summary (4 min read)
1 Introduction
- Many empirically determined parameters are included in the complicated process-based models, which leads to large uncertainties.
- Moreover, current fire models have deficiencies in simulating the peak fire month, fire season length, and interannual variability, as reported by the Fire Model Intercomparison Project (FireMIP; Hantson et al., 2020; Li et al., 2019).
- The SSiB4/TRIFFID-Fire model updates fire-induced carbon loss every 10 d, which has been rarely employed in current process-based fire models, and is used to provide a quantitative assessment of fire impact on ecosystem characteristics and surface energy at subseasonal, seasonal, interannual, and long-term scales.
2.1 Land and vegetation model
- The Simplified Simple Biosphere Model (SSiB, Xue et al., 1991; Zhan et al., 2003) is a biophysical model which simulates fluxes of surface radiation, momentum, sensible and Geosci.
- Latent heat, runoff, soil moisture, surface temperature, and vegetation gross and net primary productivity (GPP and NPP) based on energy and water balance.
- The SSiB was coupled with a dynamic vegetation model, the Top-down Representation of Interactive Foliage and Flora Including Dynamics Model , to calculate leaf area index (LAI), canopy height, and PFT fractional coverage according to the carbon balance (Cox, 2001; Zhang et al., 2015; Harper et al., 2016; Liu et al., 2019).
- The authors have improved the PFT competition strategy and plant physiology processes to make the SSiB4/TRIFFID suitable for seasonal, interannual, and decadal studies (Zhang et al., 2015; Liu et al., 2019).
2.2 Fire model and modifications
- A process-based fire model of intermediate complexity has been implemented in the SSiB4/TRIFFID, called SSiB4/TRIFFID-Fire.
- The fire model is comprised of three parts: fire occurrence, fire spread, and fire impact.
2.2.1 Fire occurrence
- Nf is the product of the number of potential ignition counts due to both natural causes (In, count s−1 km−2) and human ignitions (Ia, count s−1 km−2), fuel availability (fb), fuel combustibility (fm), and human suppression factors (feo).
- The authors only consider non-crop fire by excluding the cropland fraction from burning: Nf = (In+ Ia)fbfmfeo(1− )Ag, (2) where Ag is the land area of the grid cell (km2).
- Fires in the croplands are excluded here due to their small extent within the major fire regions and their relatively low intensity (Bistinas et al., 2014).
2.2.2 Average spread area after fire ignition
- The average spread area of a fire is assumed to be elliptical in shape, with the ignition point located at one of the foci and the fastest spread occurring along the major axis.
- The human suppression factor, Fse, reflects the human impact on fire spread through firefighting activities and is parameterized following Li et al. (2013).
2.2.3 Carbon emissions, post-fire mortality, and emissions of aerosols and trace gases
- As the carbon cycle in current SSiB4/TRIFFID does not explicitly represent the litter carbon storage and decomposition, the authors assume the litter carbon and woody debris account for about 25 % of aboveground biomass for global forest and about 30 % for savanna and grassland based on previous studies (Pan et al., 2011; de Oliveira et al., 2019).
- CCj is the corresponding combustion completeness for leaf, stem, root, and litter of the j th PFT (Table S2).
- Meanwhile, fire-induced mortality transfers carbon from uncombusted leaf, stem, and root to litter: ψj = AbCj · (1−CCj ) ·Mj , (14) where Mj = (Mleaf,Mstem,Mroot)j is the corresponding mortality factor for the j th PFT(Table S2).
- The emissions of trace gases and aerosols can be applied in the atmospheric chemistry model to calculate the production of secondary aerosols, transport of pollutants, and the resultant aerosol direct and indirect effects on climate.
2.2.4 Including the fire effect on the carbon pool
- When the Li et al. (2012) fire model is coupled with CLM, the vegetation distribution is prescribed using satellite-based Geosci.
- In SSiB4/TRIFFID-Fire, the fire-related carbon loss due to combustion and post-fire mortality is transferred to changes of PFT fraction based on carbon balance.
- Cvj is the carbon density (g C km−2); cij is the competition coefficient between the ith and j th PFTs; and γj (s−1) is the constant disturbance representing the loss of PFT fraction due to fires, pests, windthrow, and many other processes.
- The fire-caused PFT fraction loss results in bare soil for vegetation spreading decided by the competition strategy in TRIFFID.
2.3 Implementing the fire model in SSiB4/TRIFFID
- In SSiB4/TRIFFID, SSiB4 provides GPP, autotrophic respiration, and other physical variables such as canopy temperature and soil moisture every 3 h for TRIFFID (Fig. 1).
- TRIFFID accumulates the 3-hourly GPP and respiration and provides biotic carbon, PFT fractional coverage, vegetation height, and LAI every 10 d, which are used to update surface properties (albedo, canopy height, roughness length, and aerodynamic/canopy resistances) in SSiB4.
- When the fire model is included, it uses the meteorological forcings and physical variables provided by SSiB4 every 3 h and the biophysical properties (PFT fraction and biotic carbon) provided by TRIFFID every 10 d. TRIFFID updates the vegetation dynamics based on carbon balance on day 10, using the net primary production, fireinduced carbon loss, and PFT competition strategy.
- The updated vegetation dynamics are transferred to SSiB4 to reflect fire effects on surface properties.
3.1 Experimental design
- Two sets of offline experiments have been conducted using SSiB4/TRIFFID-Fire, which consist of FIRE-ON (SSiB4/TRIFFID-Fire with fire model switched on) and FIRE-OFF (SSiB4/TRIFFID-Fire with fire model switched off).
- To obtain the initial conditions for these two experiments, similar to their previous SSiB4/TRIFFID experiments (Zhang et al., 2015; Liu et al., 2019), the authors conducted spinup simulations (SPFIRE-ON and SPFIRE-OFF) for 100 years to reach a quasi-equilibrium PFT distribution with and without fire disturbance.
- These spin-up simulations were initialized using the quasi-equilibrium state from Liu et al. (2019) (SPINIT in Fig. 2) and were driven by climatology forcing averaged over 1948–1972 and atmospheric CO2 concentration, population density, and GDP in 1948 (Fig. 2).
- FIRE-ON produced the fire regime, ecosystem, and surface conditions during 1948–2014.
- A FIRE-OFF run, based on SPFIRE-OFF, was carried out with the fire model switched off during 1948–2014.
3.2 Model input and validation data
- The meteorological forcings used to drive SSiB4/TRIFFIDFire for the period of 1948–2014 are from the Princeton global meteorological dataset for land surface modeling (Sheffield et al., 2006), including surface air temperature, surface pressure, specific humidity, wind speed, downward shortwave radiation flux, downward longwave radiation flux, and precipitation (Table 3).
- The dataset is constructed by combining global observation-based datasets with the NCEP/NCAR reanalysis.
- The 2-hourly climatology lightning flashes data from NASA LIS/OTD v2.2 at 2.5◦× 2.5◦ resolution are used to calculate the number of natural ignitions.
- The GFED fire product provides the burned area and fire emissions on the global scale and has been widely used for fire model validation and calibration (van Marle et al., 2017; Li et al., 2019).
4 Results
- This section evaluates the model simulation of burned area, carbon emissions, PFT fraction, and GPP by comparing FIRE-ON results with GFED4s, GLC2000, and FLUXNET-MTE data.
- Specifically, the authors will focus on the model description of fire seasonality.
- -Fire is applied to assess the long-term fire effect on the ecosystem and surface energy budget using the differences between the FIRE-ON and FIRE-OFF.
4.1 Burned area
- The simulations of burned area are evaluated using satellitebased product GFED4s for the period of 2000–2014.
- In 8 out of the 14 sub-regions, SSiB4/TRIFFID-Fire reproduces the observed interannual variability (IAV) of burned area well, with the correlation between simulations and observations significant at p < 0.05 (Fig. 4c).
- The simulated IAV of SHAF burned area is not as good as other savanna fire regions (e.g., NHAF, SHSA, and AUST), although the IAV is small there (Fig. 4c).
- As lightning flash is the predominant ignition source in the Northern Hemisphere high latitudes, the application of climatology lightning has a greater impact in boreal regions than in other parts of the globe.
- -Fire generates two fire seasons in June–July–August and December–January–February, capturing the peak month in August but underestimating the burned area in December.
4.2 Fire emissions
- Biomass burning emissions are determined by burned area, fuel combustion rate per unit area, and emission factors per unit mass of fuel burned (van der Werf et al., 2017).
- -Fire captures the high carbon emissions in tropical savannas, the intermediate emissions in Northern Hemisphere boreal forests, and the low emissions in hu- https://doi.org/10.5194/gmd-13-6029-2020.
- The simulated global CO emission is 433.7.
- In NHAF, SSiB4/TRIFFID-Fire reproduces the large CO emission in DJF, although the model slightly underestimates CO emis- sion in December and overestimates it in February.
4.3 PFT distribution and GPP
- The simulation of vegetation coverage, which represents model description of biomass allocation and influences the fuel availability and flammability in fire modeling, is evaluated against GLC2000 (Bartholome and Belward, 2005).
- Compared with the observations (Fig. 10a), SSiB4/TRIFFID-Fire captures the fractional coverage of trees in the Amazon rainforest, tropical Africa, equatorial Asia, Southeast Asia, southeastern North America, and the Northern Hemisphere Boreal regions (Fig. 10c).
- Shrubs are primarily located in the semiarid regions and the pan-Arctic area and tundra is limited to the pan-Arctic area and Tibetan Plateau (Fig. S3f–g).
- The highest GPP occurs in the tropical evergreen forest and decreases with latitude in both observations and model.
- In addition, the correlation of IAV of global GPP is 0.68 (p < 0.05) between SSiB4/TRIFFID-Fire and FLUXNET-MTE, indicating that the model has reasonably captured the terrestrial ecosystem variability during the historical period.
4.4 Fire effects on ecosystem characteristics and surface properties
- The authors investigate long-term fire effects on ecosystem characteristics, surface properties, and surface energy budget using the differences between FIREON and FIRE-OFF (FIRE-ON minus FIRE-OFF).
- In Africa, the simulated fire effects on vegetation structure (tree and grass cover, LAI, and vegetation height) peak in the tropical savanna surrounding the forests and gradually decrease towards the deserts (Fig. 11).
- Therefore, the latent heat can be increased due to enhanced soil evaporation.
5 Conclusions and discussions
- The authors have implemented a process-based fire model of intermediate complexity into a DGVM, SSiB4/TRIFFID, which is based on the surface water, carbon, and energy balances, as well as the PFT competition.
- The SSiB4/TRIFFID-Fire is then applied to study the long-term fire effects on ecosystem characteristics and surface energy.
- As Li et al. (2017) is the only modeling study investigating the long-term fire effects on the land energy budget, their simulation provides another approach that quantifies fire effects using a different land surface model with different approaches in parameterizing some land surface processes and vegetation dynamics.
- All authors (HH, YX, FL, and YL) have contributed to the analysis methods and to the text.
- The authors acknowledge the use of the Cheyenne supercomputer (https://doi.org/10.5065/D6RX99HX, Computational and Information Systems Laboratory, 2019), provided by NCAR CISL, for providing HPC resources.
Did you find this useful? Give us your feedback
Citations
948 citations
21 citations
17 citations
12 citations
6 citations
References
12 citations
12 citations
6 citations
4 citations
Related Papers (5)
Frequently Asked Questions (17)
Q2. What have the authors stated for future works in "Modeling long-term fire impact on ecosystem characteristics and surface energy using a process-based vegetation–fire model ssib4/triffid-fire v1.0" ?
It reasonably reproduces the global GPP and PFT distribution, which is important to study fire effects on the ecosystem. The SSiB4/TRIFFID-Fire is then applied to study the long-term fire effects on ecosystem characteristics and surface energy. HH conducted the simulation with suggestions from FL and YL.
Q3. What are the sources of uncertainties in the simulation of fire effects?
Other sources of uncertainties include the differences in the partitioning between latent heat and sensible heat fluxes in land surface models, the differences in the parameterization of the evaporation processes, and the changes due to atmospheric feedbacks, such as cloud cover and precipitation changes.
Q4. How much vegetation has been reduced by fire?
Over the African and South American savanna, the authors find fire has reduced the area-averaged LAI and vegetation height by 0.52 m2 m−2 (12.5 %) and 5.76 m (49.1 %), respectively.
Q5. How have fire models been used to assess the long-term impact on the terrestrial carbon cycle?
Some fire models have been used to assess long-term fire impact on the terrestrial carbon cycle by comparing a reference simulation with fire and a sensitivity simulation representing “a world without fire”.
Q6. What is the carbon density vector for the j th PFT?
In TRIFFID (Cox, 2001), the fractional change of the j th PFT ( dfjdt ) is governed by the Lotka–Volterra equation:dfj dt = λjNPPj fj Cvj[ 1− ∑ j cijfj ] − γjfj , (16)where fj is the fractional coverage of the j th PFT; λjNPPj is the carbon available for spreading;
Q7. What is the reason for the lack of a realistic fire season in some regions?
The inaccurate simulation of fire season in several fire regions could come from deficiency of the forcing data, the inaccuracy in dynamic vegetation processes, or some processes that control the fire but are not represented in the model.
Q8. What is the corresponding mortality factor for the j th PFT?
When the fire model is coupled to SSiB4/TRIFFID, the loss of PFT fraction due to fires (βj ) can be explicitly derived from the fire-induced carbon loss:βj = (ϕj +ψj ) · fjCvj , (17)where ϕj and ψj are PFT-dependent carbon loss due to combustion and post-fire mortality, respectively.
Q9. What is the average spread area of a fire?
The average spread area of a fire is assumed to be elliptical in shape, with the ignition point located at one of the foci and the fastest spread occurring along the major axis.
Q10. What is the feo factor used in the Li et al. (2012) fire?
In the Li et al. (2012) fire scheme, this factor (fθ ) is parameterized using root zone soil moisture potential factor β (0–1.0), a model-dependent variable used to calculate transpiration in CLM (Li and Lawrence, 2017).
Q11. What is the carbon density of the ith and j th PFT?
Cvj is the carbon density (g C km−2); cij is the competition coefficient between the ith and j th PFTs; and γj (s−1) is the constant disturbance representing the loss of PFT fraction due to fires, pests, windthrow, and many other processes.
Q12. What is the carbon density vector for leaf, stem, root, and litter of the j?
(13)Cj = (Cleaf,Cstem,Croot,Clitter)j is carbon density vector (g C km−2) for leaf, stem, root, and litter of the j th PFT calculated in TRIFFID.
Q13. What is the average burned area in a savanna?
The burned area in African savanna accounts for more than 60 % of the global burned area in both GFED4s and SSiB4/TRIFFID-Fire (Fig. 4b).
Q14. What are the main reasons why fire models have been used to reconstruct fire history?
Fire models have been used to reconstruct fire history before the satellite era (Yang et al., 2015; van Marle et al., 2017; Li et al., 2019).
Q15. What is the variable to describe the dependence of fuel combustibility on the preceding climate?
In SSiB4/TRIFFID-Fire, the root zone soil moisture θ is found to be the best variable to describe the dependence of fuel combustibility on the preceding climate.
Q16. What is the effect of fire on vegetation structure in the savanna?
Their results are consistent with the long-term fire experiments that reported that fire strongly affected vegetation structure, lowering the proportions of trees to fire-resistant grasses and reducing the vegetation height and aboveground biomass (Shackleton and Scholes, 2000; Higgins et al., 2007; van Wilgen et al., 2007; Furley et al., 2008; Smit et al., 2010; Devine et al., 2015), and that fire impact is more significant in wetter savanna than in drier savanna (Moreira, 2000; Sankaran et al., 2005).
Q17. What is the spatial correlation of the GFED4s and TRIFFID-Fi?
In general, the spatial distribution of carbon emissions coincides with that of the burned area: SHAF, NHAF, and SHSA are the major fire emission regions and they contribute to 65.4 % of the total emission in both GFED4s and SSiB4/TRIFFID-Fire (Fig. 8a).