Reduced turnover rate of topsoil organic carbon in old-growth forests:a case study in subtropical China

Xin Xiong ,Juxiu Liu,Guoyi Zhou,Qi Deng,Huiling Zhang,Guowei Chu,Ze Meng and Deqiang Zhang*

Abstract

Keywords: Forest succession,Old-growth forests, Soil organic carbon, Turnover rate, 13C natural abundance

Background

Old-growth forests are considered the climax ecosystems(Clements 1916), and the age threshold varies from 100 to 500 years for different forest biomes (Jarvis 1989;Goulden et al. 2011; Liu et al. 2014). Increasing evidence shows that old-growth forests can continue to accumulate carbon (C) in soils (Zhou et al. 2006a; Luyssaert et al. 2008), contrary to the longstanding view that they are C neutral (Odum 1969; Jarvis 1989). However, the mechanism of the soil C sink in old-growth forests remains uncertain. Plant C input and the soil organic carbon (SOC) turnover rate are commonly considered two key factors in determining the C sequestration capacity of soil (Paul 2016; Nyawira et al. 2017). The relationship between litter inputs and soil C storage has been well examined in forest C cycle studies (Huang et al. 2011;Tefs & Gleixner 2012; Zhou et al. 2019), whereas the SOC turnover rate has become the dominant uncertainty in the C sequestration capacity of forest soils.Therefore, accurately quantifying the soil C turnover rate and its relationship with forest development is crucial to understanding the mechanisms of soil C sinks and predicting future C budget dynamics in the world’s forests.

A recent study revealed that old-growth forest soils have higher organic C inputs than young forest soils,which was caused by the higher proportion of products transferred to the soil during the decomposition of highquality litter in old-growth forests (Xiong et al. 2020a).Considering the emerging view that high-quality litter inputs more easily form stable soil organic matter through microbial assimilation and mineral protection(Cotrufo et al. 2013; Xu et al. 2014), whether SOC has a longer residence time in older forests becomes an intriguing question. Forest succession is the coevolution of above-and belowground systems.A rich body of literature demonstrates changes in the microclimate(Sánchez-Reyes et al. 2020), tree composition and stand structure (Gray et al. 2016), plant productivity (Goulden et al. 2011), soil microbial diversity and function(Ren et al. 2019),soil nutrient availability(Lu et al.2015),and soil physicochemical properties (Li et al. 2013) following forest succession.These variables inevitably exert direct and indirect impacts on the biogeochemical cycling of C in the soil, accompanied with the alteration of input and output processes (Soko?owska et al. 2020). It remains unclear,however, how the variables affect the turnover rate of SOC along forest succession.

As the current understanding of soil C cycle processes is mainly based on first-order kinetics(Manzoni &Porporato 2009), the soil C turnover rate is commonly defined as the ratio of the C flux (heterotrophic respiration, RH) to the C pool(SOC storage)(Chen et al.2013).Accurate estimation of RHin situ is challenging, which impedes related studies on SOC turnover. The alternative approach is to collect field soils for laboratory incubation(Vesterdal et al.2012;Li et al.2018),but it is difficult to reproduce the real scenarios in the field (Feng et al. 2016). The13C natural abundance(δ13C) is an elegant tracer to evaluate the SOC turnover rate in situ,as the changes in δ13C with soil depth could reveal some important information on soil C turnover not evident based on elemental concentration analysis alone(Garten 2006; Diochon & Kellman 2008). Soil δ13C values generally tend to increase with depth in well-drained and undisturbed soils, concomitant with decreasing C concentrations(Balesdent et al.1993;Wang et al.2017).The slope of the linear regression (β) between the log-transformed C concentration and δ13C value in the soil profile has been proposed as a proxy for the SOC turnover rate in a number of studies (Acton et al. 2013; Gautam et al. 2017; Wang et al.2018a).

In this study, we investigated the changes in C concentrations and δ13C values along soil profiles (0–30 cm) in an undisturbed old-growth forest and its two successional forests of natural restoration in subtropical China,which allowed us to examine the effects of forest succession on SOC turnover rates under the same regional climate. The three successional forests were a P.massoniana forest (PF; ～60 years old), coniferous and broadleaved mixed forest (MF; ～100 years old) and monsoon evergreen broadleaved forest (BF; ＞400 years old). The SOC accumulation rate increased in the order of PF ＜MF ＜BF (Huang et al. 2011), although soil respiration was shown to increase with forest succession at the same site (Yan et al. 2009). However, how the SOC turnover rate changes with forest succession remains unclear. We hypothesized that the SOC turnover rate declines with forest succession, as the specific RH(the amount of respired-C produced per unit microbial biomass C) decreases with the progressive maturity of the three forests (Huang et al. 2016).

Methods

Study sites

The present study was conducted at the Dinghushan Biosphere Reserve (23°09′21″–23°11′30″ N, 112°30′39″–112°33′41″ E) in Guangdong Province, China(Fig. 1). The reserve covers an area of 1155 ha with elevations ranging from 14 to 1000 m. The region is characterized by a typical subtropical monsoon climate. The mean annual temperature is 22.3°C, with a mean monthly temperature of 13.9°C in January and 28.9°C in August. The mean annual precipitation is approximately 1927 mm, of which nearly 80% falls in the hot-humid season (from April to September). The mean annual evaporation and relative humidity are 1115 mm and 78%, respectively. The bedrocks are sandstone and shale belonging to the Devonian Period, and clay minerals are mainly kaolinite and halloysite. Soils are classified in the Ultisol order and Udult suborder according to the USDA soil classification system (Buol et al. 2003).

Fig.1 The map of the study area and forest sites.PF,P.massoniana forest;MF, coniferous and broadleaved mixed forest; BF,monsoon evergreen broadleaved forest

There are three types of natural forest communities in this reserve: PF (～60 years old), MF (～100 years old)and BF (＞400 years old). The BF is the climax community in this climatic region, and the PF and MF represent its early- and mid-successional stages, respectively. The dominant species in the BF included Aporusa yunnanensis, Blastus cochinchinensis, Castanopsis chinensis, Cryptocarya chinensis, Gironniera subaequalis and Schima superba. The PF was dominated by P. massoniana. In the MF, the dominant species included Castanopsis chinensis, Schima superba and P. massoniana; the biomass of P. massoniana continuously declined naturally, and the stand structure gradually converged to the BF.

Field sample collection

In October 2015, soil samples of 0–30 cm were collected from eight random plots (10 m×10 m) in each forest. In each plot, ten soil cores were randomly collected at three depths (0–10, 10–20, and 20–30 cm) and then pooled by depth. Each soil sample was divided into two parts: one part was sieved through a 2-mm mesh while plant materials were removed, and stored at 4°C for the determination of soil microbial biomass carbon (MBC) and other chemical analyses; and the other part was air-dried for the fractionation of water-stable aggregates.

Water-stable aggregates

Soil aggregates were divided into four aggregate size classes:large macroaggregates(＞2000 μm),small macroaggregates (250–2000 μm), microaggregates (53–250 μm), and silt & clay particles (＜53 μm). Aggregate separation was performed by a wet sieving approach adapted from Six et al.(1998).

The stability of the soil aggregate structure was expressed as the mean weight diameter (MWD) of all aggregate size fractions, and the parameter was calculated with the following equation(Kemper&Rosenau 1986):

Laboratory analyses

The air-dried soil samples were measured for pH by a pH meter (PH510, Alalis, Shanghai, China) using a 1: 2.5 soil-water suspension. The soil water content (SWC)was measured by oven-drying 15 g of fresh soil sample at 105°C for 24 h. Soil MBC was determined by subjecting fresh soil samples to the chloroform fumigationextraction method (Vance et al. 1987). Soil readily oxidized organic carbon (ROC) was measured by 333 mM KMnO4oxidation-colorimetry (Blair et al. 1995). Soil available nitrogen (AN) was analyzed by alkali-diffusion method (Liu et al. 1996). The subsamples were analyzed for total organic C and total nitrogen (TN) contents and for C isotopic compositions by an isotope ratio mass spectrometer (IsoPrime 100, IsoPrime, Manchester, UK)connected to an elemental analyzer (Vario isotope cube,Elementar, Hanau, Germany). There was no treatment to remove inorganic C before elemental and isotopic analysis in the present study, since the soils were highly acidic (Table 2) and the inorganic C contents were negligible.

The C isotope ratio (13C/12C) was expressed in δ notation (in ‰ units):

where Rsamis the13C/12C ratio of the sample and Rstdis the13C/12C ratio of the Vienna Pee Dee Belemnite standard (Coplen et al. 2006). Analysis of internal laboratory standards ensured that the precision of the measurements was ±0.1‰ for δ13C.

Recalcitrance index

In the present study, we defined the ‘recalcitrance index’(RI) of soil organic matter as follows:

Soil bulk density and SOC stock

Soil bulk density (BD) was determined in parallel with field sampling. Undisturbed soil cores were taken from three random plots at three depths (0–10, 10–20, and 20–30 cm) in each forest using a stainless steel corer(5.65 cm in diameter, 4 cm in depth). All soil cores were oven-dried at 105°C to a constant weight. Soil bulk density was calculated as the ratio of total dry weight to total soil volume.

SOC stock is the product of SOC concentration, layer thickness, and bulk density (Post et al. 2001). The SOC stock can be calculated as follows:

where H and D are the soil layer thickness (cm)and bulk density (g·cm-3),respectively.

β value

Enrichments in δ13C values with increasing depths in well-drained undisturbed soils concomitant with decreasing C concentrations have been broadly reported(Balesdent et al. 1993; Acton et al. 2013). The vertical changes in δ13C along soil profiles have been proven to be related to soil C turnover (Powers & Schlesinger 2002). The slope of a linear regression relating soil δ13C to the log-transformed C concentration was used to describe the SOC turnover rate at different spatial scales(Gautam et al. 2017; Wang et al. 2018a):

The slope of the regression (β) is a proxy of the SOC turnover rate; a lower β (steeper slope) is associated with a higher turnover rate, and b is a constant.

Data analyses

All statistical analyses were conducted with SPSS statistical software (SPSS Statistics 20, IBM Corp.). Data were tested for homogeneity of variances and normal distribution before applying analysis of variance. One-way ANOVA was conducted to examine the differences in δ13C, SOC, ROC, RI, MBC, TN, AN, C/N, SWC, pH,MWD, and BD of the upper 30-cm among different forests or different depths of same forest. Tukey’s multiple comparison test was further conducted to separate differences among means. Ordinary least squares regression was used to examine the relationships between the β value and soil properties. Statistical significance was determined at p ＜0.05 unless otherwise stated.

Results

SOC concentration and δ13C

Both the SOC concentration and soil δ13C value showed significant differences among forests and soil depths(p ＜0.05; Table 1). The SOC concentration at each depth in the BF was similar to that in the MF (p ＞0.05)but was significantly higher than that in the PF (p ＜0.05;Table 1). The SOC concentrations significantly decreased with depth in all three forests (p ＜0.05; Table 1).The soil δ13C value at each depth was lowest in the BF,and there was no significant difference between the PF and MF (p ＞0.05; Table 1). Contrary to the SOCconcentration, the soil δ13C values significantly increased with depth in each forest (p ＜0.05; Table 1).

Table 1 Concentrations and isotope abundance (δ13C)of soil organic carbon (SOC)in the upper 30 cm depth of different forests

SOC stock and β value

As shown in Fig. 2, the SOC stock of the upper 30-cm was lowest in the PF, and there was no significant difference between the MF and BF (p ＞0.05). The estimated β value increased in the order of PF (-7.04±0.48)＜MF(-6.66±0.42)＜BF (-5.22±0.41), and the value in the BF was significantly higher than that in the PF (p ＜0.05;Fig. 2).

Soil properties

Fig.2 Soil organic carbon(SOC)stocks and β values in the upper 30-cm of soils in different forests.PF,P.massoniana forest;MF,coniferous and broadleaved mixed forest; BF,monsoon evergreen broadleaved forest. Different lowercase letters indicate significant differences among different forests(p ＜0.05).The error bars indicate standard errors of the mean(n=8)

Table 2 Soil properties of the upper 30-cm in different forests

Significant differences were detected among different forests for most soil physicochemical properties (p ＜0.05; Table 2). The SOC concentration, ROC, RI, MBC,TN, SWC, and MWD were lowest in the PF, and there were no significant differences between the PF and MF(p ＞0.05; Table 2). Soil AN significantly increased in the order of PF ＜MF ＜BF (p ＜0.05; Table 2). The C/N ratio of soil was highest in the MF, and there was no significant difference between the PF and BF (p ＞0.05; Table 2). The soil pH and BD in the BF were similar to those in the MF (p ＞0.05) but were significantly lower than those in the PF (p ＜0.05; Table 2).

The relationships between the β value and soil properties

The results of regression analyses showed that the β value of soil profile was positively correlated with RI(R2=0.38, p ＜0.01; Fig. 3a), TN (R2=0.27, p ＜0.05; Fig.3b), and AN (R2=0.26, p ＜0.05; Fig. 3c) but negatively correlated with pH (R2=0.34, p ＜0.01; Fig. 3d). However, the relationships between the β value and SOC,MBC, C/N, MWD and SWC were not significant (p ＞0.05; Figs.3e–i).

Discussion

SOC storage and turnover following forest succession

Forest age is an important factor affecting SOC storage(Marin-Spiotta & Sharma 2013). Our results revealed that the surface (0–30 cm) SOC stock in the 100-yearold MF reached the level of the primary forest, whereas the SOC stock in the 60-year-old PF was still at a low level. However, this result does not mean that forest soil C stocks will peak within 100 years. Increasing evidence has confirmed that old-growth forests can continue to accumulate C in soils for centuries or longer (Zhou et al.2006a; Luyssaert et al. 2008), even as net primary productivity declines (Xiong et al. 2020a). It is worth noting that the mechanisms underlying SOC accumulation in old-growth forests are distinct from those in young forests (Zhou et al. 2006b). In early-successional forests,the increase in the soil C pool mainly results from the continuous increase in biomass,and the gained C is concentrated in the topsoil layer (Xiong et al. 2020a). As the forest matures, SOC accumulation turns to the deeper soil layers, owing to the low C/N ratio of litter and rhizospheric material, and the moist soil environment. The former facilitates a high proportion of organic C yield input to the soil during decomposition, and the latter contributes to the migration of organic C along the whole soil profile (Tefs & Gleixner 2012; Xiong et al.2020a).

Our results showed the enrichment of soil δ13C with depth in all three forests, concomitant with decreasing C concentrations (Table 1). This result is consistent with earlier observations in well-drained forest soils (Powers& Schlesinger 2002; Acton et al. 2013; Wang et al.2018a). The increasing trend in δ13C with soil depth might be attributed to two main processes: mixing of different C sources and soil C decomposition (Diochon& Kellman 2008; Acton et al. 2013). The Suess effect describes the isotopic depletion of atmospheric CO2as a result of the combustion of13C-depleted fossil fuels and biomass burning (Friedli et al. 1986), which could have caused the increasing trend in δ13C with soil depth(more13C-enriched older SOC in the deep and more13C-depleted recent SOC at the surface) (Diochon &Kellman 2008). However, Torn et al. (2002) showed that the C isotopic signatures of a modern soil profile (collected in 1997)are similar to its 100-year-old archive soil(collected sometime between 1895 and 1903) from the same site, indicating that the change of soil δ13C with depth is not due to the isotopic depletion of atmospheric CO2. A few studies proposed that roots are generally more13C-enriched than leaves or branches from the same plant, and so the soil δ13C value at the surface may be lower than that in the deeper soil (Powers & Schlesinger 2002). However, we did not find significant differences in δ13C values between roots and above-ground litter at the present site by 10 years of observation(Xiong et al. 2020b). The existing evidence implies that the contribution of the mixing of different C sources to the vertical δ13C enrichment along soil profiles should be small in the present study.

The relationship between β value and SOC turnover rate has been tested in a number of studies (Gautam et al. 2017; Wang et al. 2018a; Xiong et al. 2020b), as a lower β value implies a higher turnover rate. The results showed that the β value ranged from -7.04 to -5.22 in the present study (Fig. 2), which is in accordance with the measurements in Brazilian tropical forests (Telles et al. 2003) but lower than other observations at similar latitudes (Chiang et al. 2004; Chen et al. 2005; Richards et al. 2007), the discrepancy probably was attributed to the deeper soils (0–100 cm) considered in the above studies.

Fig.3 Relationships between the β value and soil properties.Relationship between the β value and recalcitrance index(RI,a),total nitrogen(TN,b),available nitrogen(AN,c),soil pH (d),soil organic carbon(SOC,e), microbial biomass carbon(MBC,f),carbon to nitrogen ratio(C/N, g), mean weight diameter(MWD)of soil aggregates(h)and soil water content(SWC,i)

Our results also showed that the β value increased with the progressive maturity of the three forests, indicating that the SOC turnover rate decreased during forest succession. This is consistent with a previous study that measured soil respiration components by a trenching method in the three forests (Huang et al. 2016).Although higher soil RHwas measured in the latesuccessional forest, this could be mainly caused by higher C input and SOC stock, whereas the specific RHdeclined with forest succession, and soil microbes became more efficient in conserving C resources with forest succession (Huang et al. 2016; Wang et al. 2019).Our finding is supported by another related study(Wang et al. 2018b). Based on the data consolidation of 1087 forest plots in eastern China, the researchers reported that the SOC turnover time was positively correlated with forest age. Reduced SOC turnover, combined with continuous high C input (Xiong et al. 2020a),sustains the soil C sinks in old-growth forests.

Influences of soil properties on the SOC turnover rate

Our results showed that soil physical and chemical properties differed in the different forests, but the differences in most soil properties were slight between the mid- and late-successional stages (Table 2). Soil properties are considered to play an important role in regulating the dynamics of SOC turnover across different scales (Torn et al. 1997; Xu et al. 2016). In the present study, the correlations between the β value and soil physicochemical variables were analyzed (Fig. 3), and the results showed that soil RI, nitrogen status and pH were the main controls that regulated the SOC turnover rates.

Our results showed that the β value was positively correlated with soil RI (Fig. 3a), which is in accordance with the conventional opinion that reduced SOC turnover is accompanied by enhanced soil chemical recalcitrance(Sollins et al. 1996). For a long time, the formation of recalcitrant components (e.g., humus, lignin, and long chain n-alkanes) has been considered the major pathway for SOC stabilization (Derenne & Largeau 2001; Piccolo 2001). According to the concept of ‘selective preservation’, soil organic matter is composed of both labile and relatively recalcitrant compounds, the latter being used by microorganisms only when the former are exhausted(von Lützow et al. 2006). Despite recent analytical and experimental advances demonstrating that intrinsic molecular characteristics alone do not control SOC stability:environmental and biological controls act jointly (Kleber 2010; Schmidt et al. 2011), we cannot ignore the chemical nature of organic matter in specific studies. Lavallee et al. (2019) reported that pyrogenic C decomposition is controlled by a different mechanism than biogenic C,and the latter has a shorter residence time. Nisar &Bendi (2020) conducted a field management experiment and found that the addition of compost, in addition to leading to C build-up, improved C stabilization by imparting recalcitrance to SOC.In the three successional forests, the proportion of labile constituents in litter increased with forest maturity (Huang et al.2011),and this change contributed to the formation of recalcitrant components in the soil, which was shown in other studies(Cotrufo et al. 2013;Xu et al. 2014).

In the present study, we observed that the rates of SOC turnover decreased with increasing soil nitrogen status (Figs. 3b–c). This finding is supported by previous studies in the three successional forests. Based on a seven-year nitrogen manipulation experiment, Lu et al.(2013) reported that nitrogen addition increased soil C accumulation in the upper soils by decreasing dissolved organic C leaching. Mo et al. (2008) found that high nitrogen addition (15 g N·m-2·year-1) could reduce soil respiration accompanied with a decrease in soil microbial activities. As the forest matures, the soil nitrogen status increases, and the combination of increasing soil C storage and decreasing soil respiration leads to a reduction in SOC turnover rates. Nitrogen addition has been found to impede organic matter decomposition and thus stimulate C sequestration across different ecosystems (Reay et al. 2008). A meta-analysis by Janssens et al.(2010)revealed that,averaged over 36 nitrogen manipulation studies in forest ecosystems, soil RHdeclined by 15% when nitrogen was added. Another metaanalysis showed that nitrogen fertilization decreased the soil C turnover rate in agroecosystems (Lu et al. 2011).The positive coupling relationship between soil C storage and nitrogen storage has also been reported in numerous studies (Batjes 1996; Carter et al. 1997). In the light of our results and previous studies, we speculate that the increase in global nitrogen deposition may mitigate the rise of atmospheric CO2to some extent by increasing the retention time of C in soils.

Our results showed that the β value was inversely associated with soil pH (Fig. 3d), indicating that the SOC turnover rate decreased with increasing soil acidity. This finding is in line with previous studies in cropland and grassland (Kemmitt et al. 2006; Leifeld et al. 2013). The regulation of soil pH on SOC turnover rate may be explained by the changes in substrate availability and microbial activity. A low pH could reduce substrate availability through the combination operation of physical and chemical protection, such as, decreasing surface charge density and degree of dissociation, and increasing co-precipitation with aluminum (Leifeld et al. 2013). Soil microbial growth is retarded in an acidic environment,followed by the reduction of soil respiration (Malik et al.2018). Increased soil acidity could reduce the intrinsic activity of the microbial community by a number of direct and indirect routes including: (1) reducing the amount of substrate entering the soil, (2) inhibiting certain members of the community (e.g., nitrifiers), and (3)increasing aluminum in the soil, which induces physiological toxicity (Kemmitt et al. 2006). Our previous studies are in favor of the above inferences. Wu et al. (2016)showed that prolonged acid rain facilitated organic C accumulation in the topsoil, following the inhibition of litter decomposition, soil microbial biomass, enzyme activities, and soil respiration. Furthermore, the depressed effect on soil respiration showed an increasing trend with the progressive acidification in soils (Liang et al. 2013). On the other hand, the proportion of recalcitrant components was higher in soils with lower pH,which also favors the retention and accumulation of organic C in mineral soils (Wu et al. 2020).

The results of this study showed that soil MBC had no direct effect on the β value (Fig. 3f), probably because SOC turnover is not affected by soil microbial biomass but is governed by microbial community composition and C use efficiency (Geyer et al. 2016). Unexpectedly,we did not detect direct evidence for the impact of soil moisture on the SOC turnover rate (Fig. 3i), whereas it is widely accepted as a critical regulator of SOC dynamics(Falloon et al. 2011; Green et al. 2019).Soil water content has great temporal variability (Hu et al. 2013),and the results of a single sampling may not be able to present the true soil moisture status. In the present study, soil MWD did not show a linear trend with forest development(Table 2).The formation of soil aggregate structure is controlled by the complex interplay of multiple factors, such as organic input, microbial activity, and soil chemical properties (Bronick & Lal 2005; Regelink et al. 2015), and the relationship between soil aggregates and the SOC turnover rate needs further exploration. We also recognized the importance of climatic factors in the process of SOC turnover; unfortunately, there is a lack of long-term meteorological data in individual forests. Generally, climatic factors are considered the primary controls in regulating β values (Garten et al. 2000; Acton et al. 2013).However, in a particular climate area, soil properties may play a more important role than climatic factors in SOC turnover dynamics (Doetterl et al. 2015; Wang et al.2017).Yang et al.(2015) also reported that edaphic rather than climatic variables were better determinants of δ13C enrichment from vegetation to soil on the Tibetan Plateau.

Conclusions

Our results showed that the β value in the upper 30-cm of soils increased with forest maturity, which was mainly regulated by the soil recalcitrance index, nitrogen status,and soil pH. These findings demonstrate that SOC has lower turnover rates in old-growth forests, which helps to elucidate the mechanism underlying C sequestration in old-growth forest soils. The present study emphasizes the important value of old-growth forests as a global C sink, and provides guidance in predicting the responses of forest soil C pools to global change.

Abbreviations

AN: Available nitrogen; BD: Bulk density;BF: Monsoon evergreen broadleaved forest; C:Carbon; C/N: Carbon to nitrogen ratio; MBC: Microbial biomass carbon; MF: Coniferous and broadleaved mixed forest; MWD: Mean weight diameter; PF: P. massoniana forest; RH: Heterotrophic respiration;RI: Recalcitrance index; ROC: Readily oxidized organic carbon; SOC: Soil organic carbon; SWC: Soil water content; TN: Total nitrogen; δ13C:13C natural abundance

Acknowledgements

Not applicable.

Authors’contributions

DZ conceived and designed the study;XX analyzed and interpreted the data, and wrote the primary draft; JL,GZ, and QD revised the draft; HZ and XX did the measurements;GC and ZM collected the samples. All authors read and approved the final manuscript.

Funding

The research was jointly supported by the China Postdoctoral Science Foundation (No. 2020 M682951), the National Natural Science Foundation of China (No.NSFC41773088), and the Key Research Program of the Chinese Academy of Sciences (No.QYZDJ-SSW-DQC003).

Availability of data and materials

The datasets used during the current study are available from the corresponding author on reasonable request.

Declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Author details

1Key Laboratory of Vegetation Restoration and Management of Degraded Ecosystem, South China Botanical Garden, Chinese Academy of Sciences,Guangzhou 510650, China.2Center for Plant Ecology, Core Botanical Gardens,Chinese Academy of Sciences, Guangzhou 510650, China.3School of Applied Meteorology, Nanjing University of Information Science &Technology,Nanjing 210044, China.

Received: 21 February 2021 Accepted: 29 April 2021