2.24. Fire

The fire parameterization in CLM contains four components: non-peat fires outside cropland and tropical closed forests, agricultural fires in cropland, deforestation fires in the tropical closed forests, and peat fires (see Li et al. 2012a, Li et al. 2012b, Li et al. 2013, Li and Lawrence 2017 for details). In this fire parameterization, burned area is affected by climate and weather conditions, vegetation composition and structure, and human activities. After burned area is calculated, we estimate the fire impact, including biomass and peat burning, fire-induced vegetation mortality, adjustment of the carbon and nitrogen (C/N) pools, and fire emissions.

2.24.1. Non-peat fires outside cropland and tropical closed forest

Burned area in a grid cell, \(A_{b}\) (km² s ^-1), is determined by

(2.24.1)\[A_{b} =N_{f} a\]

where \(N_{f}\) (count s^-1) is fire counts in the grid cell; \(a\) (km²) is average fire spread area of a fire.

2.24.1.1. Fire counts

Fire counts \(N_{f}\) is taken as

(2.24.2)\[N_{f} = N_{i} f_{b} f_{m} f_{se,o}\]

where \(N_{i}\) ( count s^-1) is the number of ignition sources due to natural causes and human activities; \(f_{b}\) and \(f_{m}\) (fractions) represent the availability and combustibility of fuel, respectively; \(f_{se,o}\) is the fraction of anthropogenic and natural fires unsuppressed by humans and related to the socioeconomic conditions.

\(N_{i}\) (count s^-1) is given as

(2.24.3)\[N_{i} = \left(I_{n} +I_{a} \right) A_{g}\]

where \(I_{n}\) (count km^-2 s^-1) and \(I_{a}\) (count km^-2 s^-1) are the number of natural and anthropogenic ignitions per km², respectively; \(A_{g}\) is the area of the grid cell (km²). \(I_{n}\) is estimated by

(2.24.4)\[I_{n} = \gamma \psi I_{l}\]

where \(\gamma\) =0.22 is ignition efficiency of cloud-to-ground lightning; \(\psi =\frac{1}{5.16+2.16\cos [3min(60,\lambda )]}\) is the cloud-to-ground lightning fraction and depends on the latitude \(\lambda\) (degrees); \(I_{l}\) (flash km^-2 s^-1) is the total lightning flashes. \(I_{a}\) is modeled as a monotonic increasing function of population density:

(2.24.5)\[I_{a} =\frac{\alpha D_{P} k(D_{P} )}{n}\]

where \(\alpha =0.01\) (count person^-1 mon^-1) is the number of potential ignition sources by a person per month; \(D_{P}\) (person km^-2) is the population density; \(k(D_{P} )=6.8D_{P} ^{-0.6}\) represents anthropogenic ignition potential as a function of human population density \(D_{P}\); n is the seconds in a month.

Fuel availability \(f_{b}\) is given as

(2.24.6)\[\begin{split}f_{b} =\left\{\begin{array}{c} {0} \\ {\frac{B_{ag} -B_{low} }{B_{up} -B_{low} } } \\ {1} \end{array} \begin{array}{cc} {} & {} \end{array}\begin{array}{c} {B_{ag} <B_{low} } \\ {\begin{array}{cc} {} & {} \end{array}B_{low} \le B_{ag} \le B_{up} } \\ {B_{ag} >B_{up} } \end{array}\right\} \ ,\end{split}\]

where \(B_{ag}\) (g C m^-2) is the biomass of combined leaf, stem, litter, and woody debris pools; \(B_{low}\) = 105 g C m ^-2 is the lower fuel threshold below which fire does not occur; \(B_{up}\) = 1050 g C m^-2 is the upper fuel threshold above which fire occurrence is not limited by fuel availability.

Fuel combustibility \(f_{m}\) is estimated by

(2.24.7)\[f_{m} = {f_{RH} f_{\beta}}, \qquad T_{17cm} > T_{f}\]

where \(f_{RH}\) and \(f_{\beta }\) represent the dependence of fuel combustibility on relative humidity \(RH\) (%) and root-zone soil moisture limitation \(\beta\) (fraction); \(T_{17cm}\) is the temperature of the top 17 cm of soil (K) and \(T_{f}\) is the freezing temperature. \(f_{RH}\) is a weighted average of real time \(RH\) (\(RH_{0}\)) and 30-day running mean \(RH\) (\(RH_{30d}\)):

(2.24.8)\[f_{RH} = (1-w) l_{RH_{0}} + wl_{RH_{30d}}\]

where weight \(w=\max [0,\min (1,\frac{B_{ag}-2500}{2500})]\), \(l_{{RH}_{0}}=1-\max [0,\min (1,\frac{RH_{0}-30}{80-30})]\), and \(l_{{RH}_{30d}}=1-\max [0.75,\min (1,\frac{RH_{30d}}{90})]\). \(f_{\beta}\) is given by

(2.24.9)\[\begin{split}f_{\beta } =\left\{\begin{array}{cccc} {1} & {} & {} & {\beta\le \beta_{low} } \\ {\frac{\beta_{up} -\beta}{\beta_{up} -\beta_{low} } } & {} & {} & {\beta_{low} <\beta<\beta_{up} } \\ {0} & {} & {} & {\beta\ge \beta_{up} } \end{array}\right\} \ ,\end{split}\]

where \(\beta _{low}\) =0.85 and \(\beta _{up}\) =0.98 are the lower and upper thresholds, respectively.

For scarcely populated regions (\(D_{p} \le 0.1\) person km ^-2), we assume that anthropogenic suppression on fire occurrence is negligible, i.e., \(f_{se,o} =1.0\). In regions of \(D_{p} >0.1\) person km^-2, we parameterize the fraction of anthropogenic and natural fires unsuppressed by human activities as

(2.24.10)\[f_{se,o} =f_{d} f_{e}\]

where \({f}_{d}\) and \({f}_{e}\) are the effects of the demographic and economic conditions on fire occurrence. The demographic influence on fire occurrence is

(2.24.11)\[f_{d} =0.01 + 0.98 \exp (-0.025D_{P} ).\]

For shrub and grass PFTs, the economic influence on fire occurrence is parameterized as a function of Gross Domestic Product GDP (k 1995US$ capita^-1):

(2.24.12)\[f_{e} =0.1+0.9\times \exp [-\pi (\frac{GDP}{8} )^{0.5} ]\]

which captures 73% of the observed MODIS fire counts with variable GDP in regions where shrub and grass PFTs are dominant (fractional coverage of shrub and grass PFTs \(>\) 50%). In regions outside tropical closed forests and dominated by trees (fractional coverage of tree PFTs \(>\) 50%), we use

(2.24.13)\[\begin{split}f_{e} =\left\{\begin{array}{c} {0.39} \\ {0.79} \\ {1} \end{array} \begin{array}{cc} {} & {} \end{array}\begin{array}{c} {GDP > 20 } \\ { 8 < GDP \le 20 } \\ { GDP \le 8 } \end{array}\right\} \ ,\end{split}\]

to reproduce the relationship between MODIS fire counts and GDP.

2.24.1.2. Average spread area of a fire

Fire fighting capacity depends on socioeconomic conditions and affects fire spread area. Due to a lack of observations, we consider the socioeconomic impact on the average burned area rather than separately on fire spread rate and fire duration:

(2.24.14)\[a=a^{*} F_{se}\]

where \(a^{*}\) is the average burned area of a fire without anthropogenic suppression and \(F_{se}\) is the socioeconomic effect on fire spread area.

Average burned area of a fire without anthropogenic suppression is assumed elliptical in shape with the wind direction along the major axis and the point of ignition at one of the foci. According to the area formula for an ellipse, average burned area of a fire can be represented as:

(2.24.15)\[a^{*} =\pi \frac{l}{2} \frac{w}{2} \times 10^{-6} =\frac{\pi u_{p}^{2} \tau ^{2} }{4L_{B} } (1+\frac{1}{H_{B} } )^{2} \times 10^{-6}\]

where \(u_{p}\) (m s^-1) is the fire spread rate in the downwind direction; \(\tau\) (s) is average fire duration; \(L_{B}\) and \(H_{B}\) are length-to-breadth ratio and head-to-back ratio of the ellipse; 10 ^-6 converts m ² to km ².

According to Arora and Boer (2005),

(2.24.16)\[L_{B} =1.0+10.0[1-\exp (-0.06W)]\]

where \(W\)(m s^-1) is the wind speed. According to the mathematical properties of the ellipse, the head-to-back ratio \(H_{B}\) is

(2.24.17)\[H_{B} =\frac{u_{p} }{u_{b} } =\frac{L_{B} +(L_{B} ^{2} -1)^{0.5} }{L_{B} -(L_{B} ^{2} -1)^{0.5} } .\]

The fire spread rate in the downwind direction is represented as

(2.24.18)\[u_{p} =u_{\max } C_{m} g(W)\]

(Arora and Boer, 2005), where \(u_{\max }\) (m s^-1) is the PFT-dependent average maximum fire spread rate in natural vegetation regions; \(C_{m} =\sqrt{f_{m}}\) and \(g(W)\) represent the dependence of \(u_{p}\) on fuel wetness and wind speed \(W\), respectively. \(u_{\max }\) is set to 0.33 m s ^-1for grass PFTs, 0.28 m s ^-1 for shrub PFTs, 0.26 m s^-1 for needleleaf tree PFTs, and 0.25 m s^-1 for other tree PFTs. \(g(W)\) is derived from the mathematical properties of the ellipse and equation (2.24.16) and (2.24.17).

(2.24.19)\[g(W)=\frac{2L_{B} }{1+\frac{1}{H_{B} } } g(0).\]

Since g(W)=1.0, and \(L_{B}\) and \(H_{B}\) are at their maxima \(L_{B} ^{\max } =11.0\) and \(H_{B} ^{\max } =482.0\) when \(W\to \infty\), g(0) can be derived as

(2.24.20)\[g(0)=\frac{1+\frac{1}{H_{B} ^{\max } } }{2L_{B} ^{\max } } =0.05.\]

In the absence of globally gridded data on barriers to fire (e.g. rivers, lakes, roads, firebreaks) and human fire-fighting efforts, average fire duration is simply assumed equal to 1 which is the observed 2001–2004 mean persistence of most fires in the world (Giglio et al. 2006).

As with the socioeconomic influence on fire occurrence, we assume that the socioeconomic influence on fire spreading is negligible in regions of \(D_{p} \le 0.1\) person km^-2, i.e., \(F_{se} = 1.0\). In regions of \(D_{p} >0.1\) person km^-2, we parameterize such socioeconomic influence as:

(2.24.21)\[F_{se} =F_{d} F_{e}\]

where \({F}_{d}\) and \({F}_{e}\) are effects of the demographic and economic conditions on the average spread area of a fire, and are identified by maximizing the explained variability of the GFED3 burned area fraction with both socioeconomic indices in grid cells with various dominant vegetation types. For shrub and grass PFTs, the demographic impact factor is

(2.24.22)\[F_{d} =0.2+0.8\times \exp [-\pi (\frac{D_{p} }{450} )^{0.5} ]\]

and the economic impact factor is

(2.24.23)\[F_{e} =0.2+0.8\times \exp (-\pi \frac{GDP}{7} ).\]

For tree PFTs outside tropical closed forests, the demographic and economic impact factors are given as

(2.24.24)\[F_{d} =0.4+0.6\times \exp (-\pi \frac{D_{p} }{125} )\]

and

(2.24.25)\[\begin{split}F_{e} =\left\{\begin{array}{cc} {0.62,} & {GDP>20} \\ {0.83,} & {8<GDP\le 20} \\ {1,} & {GDP\le 8} \end{array}\right. .\end{split}\]

Equations (2.24.22) - (2.24.25) reflect that more developed and more densely populated regions have a higher fire fighting capability.

2.24.1.3. Fire impact

In post-fire regions, we calculate PFT-level fire carbon emissions from biomass burning of the \(j\)th PFT, \({\phi}_{j}\) (g C s^-1), as

(2.24.26)\[\phi _{j} =A_{b,j} \mathbf{C}_{j} \bullet \mathbf{CC}_{j}\]

where \(A_{b,j}\) (km² s^-1) is burned area for the \(j\)th PFT; C_j =(\(C_{leaf}\), \(C_{stem}\), \(C_{root}\), \(C_{ts}\)) is a vector with carbon density (g C km ^-2) for leaf, stem (live and dead stem), root (fine, live coarse and dead coarse root), and transfer and storage carbon pools as elements; \(\mathbf{CC}_{j}\) = (\(\mathbf{CC}_{leaf}\), \(\mathbf{CC}_{stem}\), \(\mathbf{CC}_{root}\), \(\mathbf{CC}_{ts}\)) is the corresponding combustion completeness factor vector (Table 2.24.1). Moreover, we assume that 50% and 28% of column-level litter and coarse woody debris are burned and the corresponding carbon is transferred to atmosphere.

Tissue mortality due to fire leads to carbon transfers in two ways. First, carbon from uncombusted leaf, live stem, dead stem, root, and transfer and storage pools \(\mathbf{C^{'} _{j1}} ={(C_{{leaf}} (1-CC_{{leaf}} ),C_{{livestem}} (1-CC_{{stem}} ),C_{{deadstem}} (1-CC_{{stem}} ),C_{{root}} (1-CC_{{root}} ),C_{{ts}} (1-CC_{{ts}} ))}_{j}\) (g C km^-2) is transferred to litter as

(2.24.27)\[\Psi _{j1} =\frac{A_{b,j} }{f_{j} A_{g} } \mathbf{C^{'} _{j1}} \bullet M_{j1}\]

where \(M_{j1} =(M_{{leaf}},M_{{livestem,1}},M_{{deadstem}},M_{{root}},M_{{ts}} )_{j}\) is the corresponding mortality factor vector (Table 2.24.1). Second, carbon from uncombusted live stems is transferred to dead stems as:

(2.24.28)\[\Psi _{j2} =\frac{A_{b,j} }{f_{j} A_{g} } C_{livestem} (1-CC_{stem} )M_{livestem,2}\]

where \(M_{livestem,2}\) is the corresponding mortality factor (Table 2.24.1).

Fire nitrogen emissions and nitrogen transfers due to fire-induced mortality are calculated the same way as for carbon, using the same values for combustion completeness and mortality factors. With CLM’s dynamic vegetation option enabled, the number of tree PFT individuals killed by fire per km² (individual km^-2 s^-1) is given by

(2.24.29)\[P_{disturb,j} =\frac{A_{b,j} }{f_{j} A_{g} } P_{j} \xi _{j}\]

where \(P_{j}\) (individual km^-2) is the population density for the \(j\) th tree PFT and \(\xi _{j}\) is the whole-plant mortality factor (Table 2.24.1).

2.24.2. Agricultural fires

The burned area of cropland (km² s^-1) is taken as \({A}_{b}\):

(2.24.30)\[A_{b} =a_{1} f_{se} f_{t} f_{crop} A_{g}\]

where \(a_{1}\) (s^-1) is a constant; \(f_{se}\) represents the socioeconomic effect on fires; \(f_{t}\) determines the seasonality of agricultural fires; \(f_{crop}\) is the fractional coverage of cropland. \(a_{1}\) = 1.6x10^-4 hr^-1is estimated using an inverse method, by matching 1997-2004 simulations to the analysis of van der Werf et al. (2010) that shows the 2001-2009 average contribution of cropland fires is 4.7% of the total global burned area.

The socioeconomic factor \(f_{se}\) is given as follows:

(2.24.31)\[f_{se} =f_{d} f_{e} .\]

Here

(2.24.32)\[f_{d} =0.04+0.96\times \exp [-\pi (\frac{D_{p} }{350} )^{0.5} ]\]

and

(2.24.33)\[f_{e} =0.01+0.99\times \exp (-\pi \frac{GDP}{10} )\]

are the effects of population density and GDP on burned area, derived in a similar way to equation (2.24.32) and (2.24.33). \(f_{t}\) is set to 1 at the first time step during the climatological peak month for agricultural fires (van der Werf et al. 2010); \({f}_{t}\) is set to 0 otherwise. Peak month in this dataset correlates with the month after harvesting or the month before planting. In CLM we use this dataset the same way whether the CROP option is active or not, without regard to the CROP option’s simulated planting and harvesting dates.

In the post-fire region, fire impact is parameterized similar to section 2.24.1.3 but with combustion completeness factors and tissue mortality factors for crop PFTs (Table 2.24.1).

2.24.3. Deforestation fires

CLM focuses on deforestation fires in tropical closed forests. Tropical closed forests are defined as grid cells with tropical tree (BET and BDT tropical) coverage \(>\) 60% according to the FAO classification. Deforestation fires are defined as fires caused by deforestation, including escaped deforestation fires, termed degradation fires. Deforestation and degradation fires are assumed to occur outside of cropland areas in these grid cells. Burned area is controlled by the deforestation rate and climate:

(2.24.34)\[A_{b} = b \ f_{lu} f_{cli,d} f_{b} A_{g}\]

where \(b\) (s^-1) is a global constant; \(f_{lu}\) (fraction) represents the effect of decreasing fractional coverage of tree PFTs derived from land use data; \(f_{cli,d}\) (fraction) represents the effect of climate conditions on the burned area.

Constants \(b\) and \({f}_{lu}\) are calibrated based on observations and reanalysis datasets in the Amazon rainforest (tropical closed forests within 15.5 °S \(\text{-}\) 10.5 °N, 30.5 ° W \(\text{-}\) 91 ° W). \(b\) = 0.033 d^-1 and \(f_{lu}\) is defined as

(2.24.35)\[f_{lu} = \max (0.0005,0.19D-0.001)\]

where \(D\) (yr^-1) is the annual loss of tree cover based on CLM land use and land cover change data.

The effect of climate on deforestation fires is parameterized as:

(2.24.36)\[\begin{split}\begin{array}{ll} f_{cli,d} \quad = & \quad \max \left[0,\min (1,\frac{b_{2} -P_{60d} }{b_{2} } )\right]^{0.5} \times \\ & \quad \max \left[0,\min (1,\frac{b_{3} -P_{10d} }{b_{3} } )\right]^{0.5} \times \\ & \quad \max \left[0,\min (1,\frac{0.25-P}{0.25} )\right] \end{array}\end{split}\]

where \(P\) (mm d ^-1) is instantaneous precipitation, while \(P_{60d}\) (mm d^-1) and \(P_{10d}\) (mm d ^-1) are 60-day and 10-day running means of precipitation, respectively; \(b_{2}\) (mm d ^-1) and \(b_{3}\) (mm d ^-1) are the grid-cell dependent thresholds of \(P_{60d}\) and \(P_{10d}\); 0.25 mm d ^-1 is the maximum precipitation rate for drizzle. Le Page et al. (2010) analyzed the relationship between large-scale deforestation fire counts and precipitation during 2003 \(\text{-}\)2006 in southern Amazonia where tropical evergreen trees (BET Tropical) are dominant. Figure 2 in Le Page et al. (2010) showed that fires generally occurred if both \(P_{60d}\) and \(P_{10d}\) were less than about 4.0 mm d ^-1, and fires occurred more frequently in a drier environment. Based on the 30-yr (1985 to 2004) precipitation data in Qian et al. (2006). The climatological precipitation of dry months (P < 4.0 mm d ^-1) in a year over tropical deciduous tree (BDT Tropical) dominated regions is 46% of that over BET Tropical dominated regions, so we set the PFT-dependent thresholds of \(P_{60d}\) and \(P_{10d}\) as 4.0 mm d ^-1 for BET Tropical and 1.8 mm d ^-1 (= 4.0 mm d ^-1 \(\times\) 46%) for BDT Tropical, and \(b\)₂ and \(b\)₃ are the average of thresholds of BET Tropical and BDT Tropical weighted bytheir coverage.

The post-fire area due to deforestation is not limited to land-type conversion regions. In the tree-reduced region, the maximum fire carbon emissions are assumed to be 80% of the total conversion flux. According to the fraction of conversion flux for tropical trees in the tree-reduced region (60%) assigned by CLM4-CN, to reach the maximum fire carbon emissions in a conversion region requires burning this region about twice when we set PFT-dependent combustion completeness factors to about 0.3 for stem [the mean of 0.2\({-}\)0.4 used in van der Werf et al. (2010). Therefore, when the burned area calculated from equation (2.24.36) is no more than twice the tree-reduced area, we assume no escaped fires outside the land-type conversion region, and the fire-related fraction of the total conversion flux is estimated as \(\frac{A_{b} /A_{g} }{2D}\). Otherwise, 80% of the total conversion flux is assumed to be fire carbon emissions, and the biomass combustion and vegetation mortality outside the tree-reduced regions with an area fraction of \(\frac{A_{b} }{A_{g} } -2D\) are set as in section 2.24.1.3.

2.24.4. Peat fires

The burned area due to peat fires is given as \({A}_{b}\):

(2.24.37)\[A_{b} = c \ f_{cli,p} f_{peat} (1 - f_{sat} ) A_{g}\]

where \(c\) (s^-1) is a constant; \(f_{cli,p}\) represents the effect of climate on the burned area; \(f_{peat}\) is the fractional coverage of peatland in the grid cell; and \(f_{sat}\) is the fraction of the grid cell with a water table at the surface or higher. \(c\) = 0.17 \(\times\) 10 ^-3 hr^-1 for tropical peat fires and \(c\) = 0.9 \(\times\) 10 ^-5 hr ^-1 for boreal peat fires are derived using an inverse method, by matching simulations to earlier studies: about 2.4 Mha peatland was burned over Indonesia in 1997 (Page et al. 2002) and the average burned area of peat fires in Western Canada was 0.2 Mha yr ^-1 for 1980-1999 (Turetsky et al. 2004).

For tropical peat fires, \(f_{cli,p}\) is set as a function of long-term precipitation \(P_{60d}\) :

(2.24.38)\[f_{cli,p} = \ max \left[0,\min \left(1,\frac{4-P_{60d} }{4} \right)\right]^{2} .\]

For boreal peat fires, \(f_{cli,p}\) is set to

(2.24.39)\[f_{cli,p} = \exp (-\pi \frac{\theta _{17cm} }{0.3} )\cdot \max [0,\min (1,\frac{T_{17cm} -T_{f} }{10} )]\]

where \(\theta _{17cm}\) is the wetness of the top 17 cm of soil.

Peat fires lead to peat burning and the combustion and mortality of vegetation over peatlands. For tropical peat fires, based on Page et al. (2002), about 6% of the peat carbon loss from stored carbon is caused by 33.9% of the peatland burned. Carbon emissions due to peat burning (g C m^-2 s^-1) are therefore set as the product of 6%/33.9%, burned area fraction of peat fire (s^-1), and soil organic carbon (g C m^-2). For boreal peat fires, the carbon emissions due to peat burning are set as 2.2 kg C m^-2 peat fire area (Turetsky et al. 2002). Biomass combustion and vegetation mortality in post-fire peatlands are set the same as section 2.24.1.3 for non-crop PFTs and as section 2.24.2 for crops PFTs.

2.24.5. Fire trace gas and aerosol emissions

CESM2 is the first Earth system model that can model the full coupling among fire, fire emissions, land, and atmosphere. CLM5, as the land component of CESM2, calculates the surface trace gas and aerosol emissions due to fire and fire emission heights, as the inputs of atmospheric chemistry model and aerosol model.

Emissions for trace gas and aerosol species x and the j-th PFT, \(E_{x,j}\) (g species s^-1), are given by

(2.24.40)\[E_{x,j} = EF_{x,j}\frac{\phi _{j} }{[C]}.\]

Here, \(EF_{x,j}\) (g species (g dm)^-1) is PFT-dependent emission factor scaled from biome-level values (Li et al., in prep, also used for FireMIP fire emissions data) by Dr. Val Martin and Dr. Li. \([C]\) = 0.5 (g C (g dm)^-1) is a conversion factor from dry matter to carbon.

Emission height is PFT-dependent: 4.3 km for needleleaf tree PFTs, 3 km for other boreal and temperate tree PFTs, 2.5 km for tropical tree PFTs, 2 km for shrub PFTs, and 1 km for grass and crop PFTs. These values are compiled from earlier studies by Dr. Val Martin.

Table 2.24.1 PFT-specific combustion completeness and fire mortality factors.

PFT	CCleaf	CCstem	CCroot	CCts	Mleaf	Mlivestem,1	Mdeadstem	Mroot	Mts	Mlivestem,2	\(\xi\)j
NET Temperate	0.80	0.30	0.00	0.50	0.80	0.15	0.15	0.15	0.50	0.35	0.15
NET Boreal	0.80	0.30	0.00	0.50	0.80	0.15	0.15	0.15	0.50	0.35	0.15
NDT Boreal	0.80	0.30	0.00	0.50	0.80	0.15	0.15	0.15	0.50	0.35	0.15
BET Tropical	0.80	0.27	0.00	0.45	0.80	0.13	0.13	0.13	0.45	0.32	0.13
BET Temperate	0.80	0.27	0.00	0.45	0.80	0.13	0.13	0.13	0.45	0.32	0.13
BDT Tropical	0.80	0.27	0.00	0.45	0.80	0.10	0.10	0.10	0.35	0.25	0.10
BDT Temperate	0.80	0.27	0.00	0.45	0.80	0.10	0.10	0.10	0.35	0.25	0.10
BDT Boreal	0.80	0.27	0.00	0.45	0.80	0.13	0.13	0.13	0.45	0.32	0.13
BES Temperate	0.80	0.35	0.00	0.55	0.80	0.17	0.17	0.17	0.55	0.38	0.17
BDS Temperate	0.80	0.35	0.00	0.55	0.80	0.17	0.17	0.17	0.55	0.38	0.17
BDS Boreal	0.80	0.35	0.00	0.55	0.80	0.17	0.17	0.17	0.55	0.38	0.17
C3 Grass Arctic	0.80	0.80	0.00	0.80	0.80	0.20	0.20	0.20	0.80	0.60	0.20
C3 Grass	0.80	0.80	0.00	0.80	0.80	0.20	0.20	0.20	0.80	0.60	0.20
C4 Grass	0.80	0.80	0.00	0.80	0.80	0.20	0.20	0.20	0.80	0.60	0.20
Crop	0.80	0.80	0.00	0.80	0.80	0.20	0.20	0.20	0.80	0.60	0.20

Leaves (\(CC_{leaf}\) ), stems (\(CC_{stem}\) ), roots (\(CC_{root}\) ), and transfer and storage carbon (\(CC_{ts}\) ); mortality factors for leaves (\(M_{leaf}\) ), live stems (\(M_{livestem,1}\) ), dead stems (\(M_{deadstem}\) ), roots (\(M_{root}\) ), and transfer and storage carbon (\(M_{ts}\) ) related to the carbon transfers from these pools to litter pool; mortality factors for live stems (\(M_{livestem,2}\) ) related to the carbon transfer from live stems to dead stems; whole-plant mortality factor (\(\xi _{j}\) ).