Abstract
Biomass allocation in plants is fundamental for understanding and predicting terrestrial carbon storage. Yet, our knowledge regarding warming effects on root: shoot ratio (R/S) remains limited. Here, we present a meta-analysis encompassing more than 300 studies and including angiosperms and gymnosperms as well as different biomes (cropland, desert, forest, grassland, tundra, and wetland). The meta-analysis shows that average warming of 2.50 °C (median = 2 °C) significantly increases biomass allocation to roots with a mean increase of 8.1% in R/S. Two factors associate significantly with this response to warming: mean annual precipitation and the type of mycorrhizal fungi associated with plants. Warming-induced allocation to roots is greater in drier habitats when compared to shoots (+15.1% in R/S), while lower in wetter habitats (+4.9% in R/S). This R/S pattern is more frequent in plants associated with arbuscular mycorrhizal fungi, compared to ectomycorrhizal fungi. These results show that precipitation variability and mycorrhizal association can affect terrestrial carbon dynamics by influencing biomass allocation strategies in a warmer world, suggesting that climate change could influence belowground C sequestration.
Similar content being viewed by others
Introduction
Terrestrial carbon (C) dynamics primarily depend on plant carbon economics1,2,3. Plants allocate carbon to their above- and belowground organs, often studied as a ratio between root and shoot biomass (R/S)4. In general, the R/S reflects an optimal allocation of resources by plants and is a crucial parameter for estimating terrestrial C storage5. When global change factors induce alteration of resource supply for plants, their optimal allocation patterns may shift, resulting in deviations in R/S6. For example, drought and elevated CO2 generally enhance the biomass allocation of roots over shoots for acquiring water and nutrients7, whereas increased precipitation and nitrogen deposition usually elevate the biomass allocation to shoots for an increased light competition8,9. Climate warming has also been shown to affect net primary production in terrestrial ecosystems10,11,12, such as by increasing biomass production. If this triggers soil nutrient deficits, plant carbon economics are likely to change13,14,15. Accordingly, warming effects on R/S may depend on temperature effects on resource supply, which makes it challenging to predict warming effects on plants’ biomass allocations, particularly across different environmental conditions16,17.
The antecedent climatic conditions (e.g., mean annual temperature and precipitation) of a given area are often considered as the significant predictors of R/S variation across terrestrial plants4. However, whether the long-term correlation between R/S and climatic conditions conform to immediate R/S response to warmed environments remains unknown18,19. Apart from local climate, belowground properties of ecosystems (e.g., nutrients, root traits, and root-soil interface) can further influence the response of R/S to warming12,20,21. For instance, the mutualistic symbiosis between root and mycorrhizal fungi (MF), which is universally present in terrestrial ecosystems22, might affect patterns of plant carbon allocation depending on different nutrient foraging strategies of mycorrhizal root systems15,23,24. Plants with the same mycorrhizal association might respond in (ecologically) similar ways to global change factors, especially to alteration of resource availability in soils25. As two major types of mycorrhizas, arbuscular mycorrhizal fungi (AMF)’s hyphae penetrate the cell wall and invaginate the cell membrane within the host plant root, whereas the hyphae of ectomycorrhizal fungi (EMF) do not penetrate individual root cell. Moreover, AMF can increase the host-plant’s uptake of inorganic nutrients, while EMF provides their host plants with greater access to soil organic nutrient pools26,27. Studies have shown that the dominance of a given mycorrhizal fungi type (MFT) could vary across terrestrial biomes, which might explain the variable biomass allocations reported in different biomes28,29. Consequently, we suspect that warming effects on R/S across terrestrial biomes could also be regulated by biome-specific MFT dominance.
Here, we propose three potential hypotheses to better understand how biomass allocation in response to warming would alter the distribution of terrestrial carbon at a global scale30. First, biomass allocation between above- and belowground compartments of plants may differ due to warming. Using a meta-analysis, we test this hypothesis by examining how the intercept between R/S at warmed temperature and R/S at ambient temperature deviates from 0—we refer to it as a shift in vertical biomass allocation in terrestrial plants (Fig. 1a, b). For instance, when this intercept >0, it would imply that greater biomass is allocated downwards (i.e., vertically), meaning that there will be a greater R/S under warming than that at ambient temperature (Fig. 1a, b). Second, warming can homogenize R/S across biomes. To test this, we examine how the slope between R/S at warmed temperature and R/S at ambient temperature (as a proxy of biomes) deviates from 1—we refer to it as a shift in horizontal variability of R/S in terrestrial plants across biomes (Fig. 1c, d). For instance, when this slope < 1, it would imply that warming decreases the horizontal variability of R/S and thereby causes the homogenization of R/S among diverse biomes (Fig. 1c, d). Third, the mean annual precipitation (MAP) would be a key factor determining the response of R/S due to the direct dependency of belowground resources, such as nutrients and soil water, on MAP. To test these hypotheses, we collected results from 322 warming experiments (Supplementary Fig. 1) with 94 pairs of observations to explore the effects of warming on R/S (Fig. 1f). Using a meta-analytic technique, we show that, with a mean magnitude of 2.50 °C and the median of 2 °C (Supplementary Table 1), climate warming increases R/S (intercept > 0), whereas it decreases the variability in R/S among terrestrial biomes (slope < 1, Fig. 1e). These results suggest that plants with greater investment to shoot in ambient condition shift their allocation strategy by allocating more to root biomass in warmed conditions rather than to shoot biomass, and vice versa. This subsequently triggers warming-induced homogenization of R/S among diverse biomes (Fig. 1e). We further show that MAP and mycorrhizal association with plants correlate strongly with these patterns of R/S in terrestrial plants.
a–d Hypotheses of the warming effect on root:shoot ratio (R/S). e The actual relationship between log10-transformed R/S at experimental warming and ambient temperature. f The frequency distribution of response ratio of R/S. g The log10-transformed residuals of R/S under the warmed condition from 1:1 line along with log10-transformed R/S under ambient condition. The error bands in panel e represent the upper and lower 95% confidence intervals. Across studies in our meta-analysis, the intercept in panel e was larger than 0, implying greater biomass is allocated (vertically) downwards, i.e., greater R/S under warming than that at ambient temperature (a, e); the opposite scenario of downward allocation (i.e., upward allocation, intercept smaller than zero) is illustrated in panel (b). When the slope is smaller than 1, it would imply that warming decreases the horizontal variability of R/S and thus causes the homogenization of R/S among diverse biomes (c, e); the opposite scenario of homogenization (differentiation, slope >1) is illustrated in panel (d). The size of each dot in panel e indicates the relative weight of the individual response ratio of R/S. The p-value of the normal distribution test of RR(R/S) was 0.3818. The RR++ indicates the weighted response ratio of R/S in panel (f); the ‘n’ in panel (f) is the sample size of the response ratio of R/S.
Results
Based on 94 pairs of observations of root:shoot ratio (R/S), average warming of 2.5 °C increased the R/S by 6.1 to 8.8% (with a mean increase of 8.1%), which was significantly larger than 0 (p < 0.01, Fig. 1f). The intercept between R/S at warmed temperature and that at ambient temperature ranged from 0.044 to 0.107 (CI > 0, p < 0.05, Fig. 1e), suggesting higher biomass allocation to roots in vertical dimension in warmed environments (Fig. 1a). Meanwhile, the slope between R/S at warmed temperature and that at ambient temperature ranged from 0.908 to 0.997 (CI < 1, p < 0.05, Fig. 1e), indicating a horizontal homogenization of R/S among diverse biomes under warmed conditions (Fig. 1e). The biomes with higher R/S in ambient condition had greater variation in R/S in response to warming (Fig. 1g and Supplementary Fig. 2).
On the global scale, 164 species of vascular plants belonging to angiosperm (144) and gymnosperm taxa (20) were included in this study. Relative to annual species (23), perennial ones (141) were more common, whereas the numbers of woody (77) and herb species (87) were comparable in our meta-analysis (Fig. 2). Although the responses of R/S to warming across all the vascular plants displayed a significant phylogenetic signal (Blomberg’s K > 1, Supplementary Table 2), mean annual precipitation (MAP) was the most critical factor in determining the warming effects on R/S [i.e., response ratio (RR) of R/S, RR(R/S)] compared to all other predictors used in our analysis (Fig. 3). We found that MAP was negatively correlated with RR(R/S) (R2 = 0.104, Figs. 3 and 4, Supplementary Table 3). Moreover, the weighted RR(R/S) reversed from positive to negative in sites with MAP higher than 900 mm (Fig. 4).
This phylogenetic tree is based on 164 vascular plant species (consisting of both gymnosperms and angiosperms), which were included to estimate the weighted response ratio of root:shoot ratio, plant total biomass, above- and belowground biomass to warming. Plants in the phylogenetic tree are shown with their mycorrhizal association (AMF arbuscular mycorrhizal fungi, EMF ectomycorrhizal fungi), life form (woody or herb), and growth form (annual or perennial).
Model-averaged importance of various predictor variables for warming effects on root:shoot ratio [i.e., RR(R/S)]. b The proportion of total variance of RR(R/S) is explained by biomes, plant community composition, and mycorrhizal fungi type (MFT). The variables with an importance value >0.8 were considered as essential predictors. The importance value of a predictor is based on the sum of Akaike weights derived from the model selection using corrected Akaike’s information criteria. MFT mycorrhizal fungi type, PFT plant functional type, MAP mean annual precipitation, MAT mean annual temperature, LAT latitude, CLAY the proportion of clay in soils, BD bulk density, SOC soil organic carbon, WM warming magnitude, DUR warming duration (a). The CLAY, BD, and SOC were from the topsoil at 0–30 cm soil depth, the detailed information of these predictors is provided in Supplementary Table 3. Biomes include forest, grassland, cropland, tundra and wetland; Monocultures and mixed refer to plant communities with single plant species and multiple plant species, respectively. AMF, EMF, and AM-EMF indicated the biomes with dominant root symbiosis of arbuscular mycorrhizal (AMF), ectomycorrhizal fungi (EMF), mixed AMF, and EMF (AM-EMF), respectively (b).
a Correlations between mean annual precipitation (MAP, mm) and the response of root:shoot ratio [RR(R/S)]. b Correlations between MAP and RR(R/S) in biomes with dominant root symbiosis of arbuscular mycorrhizal and ectomycorrhizal fungi (AMF and EMF), respectively. The size of the data points indicates relative weights in each plot. The points above the line of RR(R/S) = 0 indicate that warming enhanced root biomass accumulation relative to shoot biomass (i.e., RR(R/S) > 0), and vice versa. The error bands in panels a, b represent the upper and lower 95% confidence intervals. The p-values were calculated from two-tailed tests.
Mycorrhizal fungi type (MFT: 121 for AMF and 43 for EMF, Fig. 2) was the second most important factor explaining warming effects on plant’s R/S (Fig. 3a). Besides MFT, no significant differences were found in any plant biomass variables (including R/S) in response to warming between evergreen vs. deciduous trees, broad vs. coniferous leaf trees, annual vs. perennial plants, and angiosperms vs. gymnosperms (Supplementary Table 4). Although there was a significant interaction between MFT and biomes (Supplementary Table 5), MFT still explained a greater proportion of variance for RR(R/S) relative to biome types (i.e., forest, grassland, cropland, tundra, and wetland, Fig. 3b). The sensitivity of RR(R/S) to MAP (i.e., the slopes in Fig. 4b) in biomes with ectomycorrhizal fungi (EMF) was smaller than in biomes with arbuscular mycorrhizal fungi (AMF). Moreover, both warming magnitude and warming duration didn’t affect the RR(R/S) across warming experiments (Fig. 3a and Supplementary Table 6).
Overall, warming enhanced plant total biomass (TB, +10.0%) mainly by increasing belowground biomass (BGB, +13.1%) even when it decreased aboveground biomass of plants (AGB, −8.9%, p > 0.05, Fig. 5). Warming did not change TB or BGB of plants associated with AMF but decreased their AGB by 13.8% and subsequently increased R/S by 9.9%. Warming stimulated TB, AGB, and BGB in biomes with EMF and with mixed AMF and EMF (AM-EMF), but only enhanced R/S in biomes with EMF-associated plants (Fig. 5). MFTs were significantly associated with the warming response of AGB. In contrast, the effects of biomes and plant functional types (PFTs) were true only for the response of BGB (Supplementary Table 5). In our path model, we further found that MFTS (i.e., AMF, EMF, and AM-EMF) explained 22.9% variation in AGB’s response to warming, whereas PFTS (i.e., herbs and woody plants) and MAP were the major significant pathways contributing to 25.4% of variation in BGB’s response (Fig. 6). Both responses of AGB and BGB to warming directly contributed (AGB negatively and BGB positively) to shifts in RR(R/S) in our path model, explaining a total of 46.3% of variation in RR(R/S) across studies (Fig. 6). Moreover, RR(R/S) was correlated with BGB responses to warming [RR(BGB)] (R2 = 0.31, p < 0.001), but not with AGB responses (R2 = 0.01, p = 0.392, Supplementary Fig. 2). The importance of PFTs on RR(R/S) was only true in studies with a warming-induced inhibition of TB (Supplementary Fig. 4). In addition, other experimental factors (e.g., methods of warming and belowground biomass measurement, Supplementary Table 7) and classifications of studies (e.g., field or laboratory experiments, Supplementary Table 4) did not explain any variation in RR(R/S).
a–d Weighted response ratio (RR++) of plant total biomass (TB). e–h RR++ of aboveground biomass (AGB). i–l RR++ of belowground biomass (BGB). m–p RR++ of root:shoot ratio (R/S). AMF (a, e, i, m), EMF (b, f, j, n), and AM-EMF (c, g, k, o) indicate biomes with dominant root symbiosis of arbuscular mycorrhizal fungi (AMF), ectomycorrhizal fungi (EMF), mixed AMF and EMF (AM-EMF), respectively (Supplementary Data 1). The numbers inside plots indicate sample sizes. Only the subgroups with a sample size >5 were included. The error bars indicated 95% confidence intervals (CI). When the sample size was larger than 20, the 95% CI was calculated as RR++ ± 1.96[S(RR++)], whereas when the sample size was lower than 20, we obtained 95% CI using the bootstrapping method. If 95% CI of effect sizes did not overlap with zero, the warming effect was considered to be statistically significant, and the symbol * with the blue and red color indicated statistical significance of increase and decrease, respectively, of variables in warmed conditions compared to control conditions (p-value < 0.05).
PFTS plant functional types (i.e., herbs or woody plants), MFTS dominant mycorrhizal fungi types of root symbiosis (i.e., AMF arbuscular mycorrhizal fungi, EMF ectomycorrhizal fungi, or AM-EMF mixed AMF and EMF), MAT mean annual precipitation, MAP mean annual precipitation, RR(R/S) response of root:shoot ratio. The solid and dashed lines indicate the significant (p-value < 0.05) and non-significant effects (p-value > 0.05), respectively, whereas the red and blue solid arrows indicate positive and negative effects (p-value < 0.05), respectively. We assigned ‘0’ and ‘1’ for herbs and woody plants, and assigned ‘0’, ‘0.5’, ‘1’ for AMF, AM-EMF, and EMF, respectively, in this analysis (Supplementary Data 2). Symbols *, **, *** indicate statistical significance with p-value < 0.05, <0.01, and <0.001, respectively, based on two-tailed tests. In Chi-square test for the path model, the P-value was calculated from one-tailed test.
Discussion
Our meta-analysis shows that climate warming enhances plants’ biomass allocation to belowground but mainly so in drier habitats, agreeing with the hypothesis of warming-induced vertical downward biomass allocation (Fig. 1)31. The contingency of warming effects on plant biomass via soil moisture (Supplementary Fig. 5) might be one of the primary reasons for the upregulation of R/S in many study sites, except for the tundra and alpine biomes where soils are often frozen17,32. Even in the studies with decreased total plant biomass in response to warming, we found that belowground biomass allocations were higher than aboveground ones (Supplementary Figs. 6 and 7)33,34. Such a plant strategy can have important implications for soil carbon dynamics, potentially due to a lower turnover rate of root-associated carbon (i.e., carbon in hyphae and root per se) compared to shoot-associated carbon33,34. Biomes with lower original R/S at ambient conditions showed greater allocation to belowground biomass when warming most likely had induced greater demand for resources35, agreeing with our hypothesis of horizontal homogenization of R/S across biomes in a warmer world (Fig. 1). We suspect that the upregulated R/S due to warming might alter the carbon transfer between atmosphere and soil21,36.
More specifically, our results suggest that in areas with mean annual precipitation (MAP) lower than ∼900 mm, warmed air and soil surface could induce greater belowground biomass allocation potentially for plants to exploit resources in deeper soil layers, where water is more readily available, and temperature is cooler (Fig. 7a)37. The long-term adaptation of plants to water deficiency, e.g., a higher water use efficiency, ensures more biomass accumulation under warmed conditions, thereby triggering an upregulation of R/S (Fig. 7a, and prolonging the C residence time in plants38. In areas with MAP greater than ∼900 mm, adequate soil moisture could enable positive effects of warming on nutrient turnover and availability, e.g., soil NH4+ and NO3− for AMF and EMF, respectively (Fig. 7b)15,39. Indeed, the allocation to belowground biomass in these areas declined consequently in our meta-analysis, but only for terrestrial plants associated with AMF (Fig. 7b)40.
a Biomes with low mean annual precipitation (MAP) irrespective of mycorrhizal association. b Biomes with high MAP and arbuscular mycorrhizal (AMF) or ectomycorrhizal fungi (EMF). TB plant total biomass, AGB aboveground biomass, BGB belowground biomass, MB microbial biomass, MB C/N microbial biomass C N ratio, SIN soil inorganic nitrogen, WUE water use efficiency (Supplementary Data 3, Supplementary Figs. 8 and 9). The up- and downward arrows indicate significant increment and decrement of a response variable (p-value < 0.05), whereas “ns” indicates non-significant changes under warmed conditions. “~” corresponds to no change, and “?” indicates not known.
The regulation of MAP on R/S under warmed conditions was mainly through the adjustment of belowground biomass accumulation (Supplementary Fig. 3), which responded to variable effects of warming on soil moisture among plant functional types (PFTs, Supplementary Fig. 6)12,41. Significant changes in belowground biomass (BGB) in response to warming indeed varied among PFTs compared to other categories used in our meta-analysis (Supplementary Table 5). Woody plants, which are characterized by deeper and evolutionarily mature roots42, can redistribute soil water by hydraulic lift, supporting an increment in BGB under warmed conditions (Fig. 5)43. Herbaceous plants, which are characterized by evolutionarily younger roots relative to woody plants, can potentially adjust root traits (e.g., specific root length) more readily than biomass12,41. These warming-induced changes in root traits may have contributed to their non-significant responses in terms of BGB under warmed conditions (Fig. 5). In addition to MAP, due to interaction between temperature and soil moisture, antecedent mean annual temperature (MAT) was another important background climatic factor influencing the response of R/S in our meta-analysis44. For instance, MAT displayed a negative relationship with R/S due to a warming-induced significant increment of aboveground biomass (AGB) in most of the biomes with MAT < 13 °C and plants associated with EMF (Supplementary Table 3)45.
Plant’s association with mycorrhizal fungi further regulated biomass allocation patterns of terrestrial plants in response to warming. The response of R/S to warming among 164 vascular host plants displayed a significant phylogenetic signal (Fig. 2 and Supplementary Table 2), whereas significant differences were not found between evergreen vs. deciduous trees, broad vs. coniferous trees, annual vs. perennial plants, and angiosperms vs. gymnosperms (Supplementary Table 4). Although both arbuscular mycorrhizal fungi (AMF-) and ectomycorrhizal fungi (EMF-) associated plants exhibited an overall increase in R/S, the underlying mechanisms are likely to be different (Figs. 5, 7). Benefiting from hydrolytic and oxidative extracellular enzymes produced by the hyphae of EMF46, EMF-associated plants are capable of growing better under nutrient-deficient conditions relative to AMF-associated plants13,45. EMF-associated plants showed an increasing trend of AGB, BGB, and R/S under warming (Fig. 5), likely due to hyphal-facilitated soil organic matter (SOM) degradation and root nutrient absorption47. By contrast, AMF-associated plants are often more abundant in warmer sites (e.g., in temperate grasslands and deciduous forests) with drier soils48,49,50,51. Although AMF association correlated with greater water use efficiency in plants under warmed conditions (WUE, Supplementary Fig. 7), AMF-associated plants seemed to slow down the shoot biomass accumulation due to warming-induced soil water deficiency (mainly in grasses, Supplementary Fig. 9)52 and the potential increment of plant-derived carbon to AMF53. Interestingly, the regulation of mycorrhizal fungi types on responses of R/S to warming were mainly related to shoot biomass but not to root biomass (Fig. 6), implying the important roles of mycorrhizal associations on aboveground biomass allocation in terrestrial plants in a warmer world24,54,55.
As for biomes with mixed AMF and EMF (AM-EMF, Fig. 5, Supplementary Fig. 8), the complementary resource use (e.g., water and nutrients) between plants with arbuscular and ectomycorrhizal fungi provides a plausible explanation for the observed stability of R/S and higher productivity in response to warming (Fig. 4)56,57. As we did not find any significant differences in plant biomass allocation (R/S) or biomass itself (AGB and BGB) to warming between evergreen vs. deciduous trees, broad vs. coniferous trees, annual vs. perennial plants, or angiosperms vs. gymnosperms (Supplementary Table 4) when we ignored the types of mycorrhizal fungi in our analysis, it further confirms the importance of plants’ association with mycorrhizal fungi in understanding their biomass responses to warming compared to many other plant biotic characteristics.
We suspect that a stronger association of warming effects on biomass allocation with MAP and mycorrhizal associations could potentially lead to a redistribution of R/S in the horizontal dimension with more homogenization and lower variability of biomass allocation patterns among diverse biomes (Fig. 1)58. Although we found some clear signals in how warming effects were observed through a greater allocation of root biomass in terrestrial plants relative to shoot biomass, we caution readers about the variability in effect sizes, which could be caused by the distinct effect of warming magnitude and warming duration in a given biome (Supplementary Table 3). For example, warming magnitude might negatively affect the response of R/S for biomes with EMF, while warming duration could induce negative and positive effects for biomes with woody and herbaceous plants, respectively, as revealed in our analysis (Supplementary Table 3).
In summary, climate warming with a mean magnitude of 2.5 °C can enhance belowground biomass allocation as a whole, but mainly so in drier habitats. The warming of 2.5 °C also resembles the predicted IPCC warming (2.4 °C rise of global mean temperature) under a high greenhouse gas emissions scenario (SSP5-8.5) in the following decades (2041–2060)59, making our results further relevant for the future of terrestrial carbon dynamics. Our results highlight that habitat dryness (or wetness), and types of mycorrhizal association are crucial factors to consider when we aim to predict biomass allocation strategies in terrestrial plants in response to warming across biomes.
Methods
Data collection
Using the ISI Web of Science, peer-reviewed journal articles (1950–2020) related to plant biomass, growth, carbon cycling, and/or soil nitrogen under experimental warming were searched with the following specific keywords: (warming OR temperature OR heat* OR greenhouse) AND (biomass OR plant OR allocation OR root OR leaf OR stem OR photosynthe* OR growth OR aboveground OR belowground OR respiration OR soil OR carbon OR nitrogen OR microb*). To avoid bias in publication selection, we used the following criteria to select the studies for our meta-analysis (Supplementary Data 1): (i) Warming treatments were conducted in terrestrial biomes and at least one of the following variables was measured: plant total biomass (TB), aboveground biomass (AGB), belowground biomass (or root biomass, BGB), root:shoot ratio (R/S), water use efficiency (WUE), soil inorganic nitrogen (SIN), soil NH4+, soil NO3−, microbial biomass (MB), and/or microbial biomass C/N (MB C/N). (ii) The experimental temperature, warming methods (e.g., open-top chamber, infrared heater, or soil heating cable and greenhouse), and dominant plant species were indicated clearly in both warming and control groups. (iii) Apart from the difference in experimental temperature, other initial environmental conditions and plant species compositions were the same between the control and warming groups. (iv) The duration of the warming experiment was at least longer than one growing season. (v) The mean, standard deviations/errors, and sample sizes of these variables could be extracted from the figures, tables, or the text directly.
Using the above-mentioned five selection criteria, 322 papers associated with the field (280) or laboratory (42) warming experiments were selected (Supplementary Table 4, Notes and References). The methods of air and/or soil warming in selected studies included open-top chamber, infrared heater, soil heating cable, and greenhouse (Supplementary Table 6). The observations (e.g., TB, AGB, BGB, or R/S, WUE, SIN, soil NH4+, soil NO3−, MB, and/or MB C/N) in both control and warming groups were extracted using software GetData (version 2.22) as well as treatment and environmental information [e.g., warming duration (DUR) and magnitude (WM), latitude (LAT), mean annual precipitation (MAP) and temperature (MAT)]. If the local climate condition was not reported in a paper, the WorldClim dataset (http://www.worldclim.org/) was used to obtain climate variables based on latitude and longitude for that given study60. Roots or belowground biomass measurements across studies included direct harvest, soil core, and ingrowth mesh bags (Supplementary Table 6). In order to analyze the effects of climate (only field studies) and treatment variable (e.g., warming magnitude) on plant biomass responses to warming, we extracted MAP, MAT, DUR, and WM as well as the soil properties (0–30 cm depth), i.e., bulk density (BD), clay content (CLAY), and soil organic carbon (SOC) at each site from Global Gridded Surfaces of Selected Soil Characteristics (BD from IGBP-DIS)61 and the Harmonized World Soil Database (CLAY and SOC, from version 1.2, https://daac.ornl.gov/SOILS/guides/HWSD.html, Supplementary Table 3). The sites of the selected studies were mainly located in East Asia, North America, and Europe (Supplementary Fig. 1). The MAP and MAT of sites in field experiments ranged from 27 mm to 2400 mm (median 636.2 mm), and −20 °C to 30 °C (median 7.4 °C), respectively. Experimental warming duration and magnitude ranged from one growing season to 25 years (median 2 years), and 0.26 °C to 12 °C (median 2 °C), respectively (Supplementary Fig. 1 and Table 1).
The biome types (including cropland, desert, forest, grassland, tundra, and wetland, Supplementary Table 5), plant functional types (PFTs, including woody plant vs. herb, tree vs. shrub, grass vs. forb, Supplementary Fig. 1), other taxonomy categories (evergreen vs. deciduous tree, broad vs. coniferous leaf tree, annual vs. perennial plant, angiosperm vs. gymnosperm, monoculture vs. mixed plant community, Figs. 2 and 3, Supplementary Table 4) and mycorrhizal fungi types (MFTs, Supplementary Fig. 1 and Supplementary Table 5) for dominant plants in each case study were confirmed based on the original publications and the latest FungalRoot database (https://nt.ars-grin.gov/fungaldatabases/)45. The MFTs of biomes were labeled as AMF, EMF, or AM-EMF if the root symbiosis of dominant plants were arbuscular mycorrhizal fungi (AMF), ectomycorrhizal fungi (EMF), or mixed AMF and EMF (AM-EMF). Most analyses in this study focused on AMF, EMF, and AM-EMF, as biomes with non-mycorrhizal fungi (NMF) were too few to conduct any meta-analysis (Supplementary Fig. 1).
Data analysis
In this study, we employed response ratio [RR, natural log (ln) of the ratio of the mean value in warming treatment (\(\bar{{X}_{t}}\)) to that in control (\(\bar{{X}_{{{{{{\rm{c}}}}}}}}\)), Eq. (1)] to reflect warming effects on various response variables62,63.
The weighted response ratio (RR++, Eq. (2)), and the standard error of RR++ [s(RR++), Eq. (3)] in each subgroup [e.g., different MFTs (AMF, EMF, or AM-EMF) or PFTs (woody plants or herbs)] were calculated using individual RR(RRij) and its weight (wij), which is the reciprocal of the variance (vij, Eq. (4), Fig. 4).
where m is the number of subgroups, and k is the number of RR in the ith subgroup (i = 1, 2, …, m; j = 1, 2, …, k).
The \({n}_{t}\), \({n}_{c}\), and \({S}_{t}\), \({S}_{c}\) are the numbers of replicates and standard deviations for a given variable in the warming and control group, respectively. The warming effect on a given response variable was considered statistically significant when the 95% CI of the effect size did not overlap with zero. The percentage change of a response variable under warmed conditions was estimated by [exp(RR++) − 1] × 100%.
The between-group heterogeneity (Qb) of effect sizes was estimated using the Q-statistic in MetaWin 2.1 (Supplementary Table 5)64. Stepwise linear regression was used to analyse the impacts of environmental conditions on the responses of root:shoot ratio [RR(R/S), Supplementary Table 4]. We used Blomberg’s K values to test the phylogenetic signals on the warming response of plant biomass variables using the “picante” package (version 1.8.2)65 in R (version 4.1.3, R Core Team, 2018, Supplementary Table 2). The effects of biomes, MFTs, PFTs, experimental methods (warming methods and methods to measure root or belowground biomass), warming magnitude (WM), and warming duration (DUR) were examined to explain the variation in response variables using the analysis of variance (ANOVA, Supplementary Tables 5 and 6). Both WM and DUR were sub-grouped into six classes according to the method explained in Lin et al. (Supplementary Table 6)12. The importance of each predictor for the response ratio of a given response variable was expressed as the sum of Akaike weights derived from the model selection using the “glmulti” package (version 1.0.8)66 in R (Fig. 3a, Supplementary Fig. 3). Nested analysis was performed using the “nlme” package (version 3.1-158)67 in R to estimate the variance explained by MFTs (Fig. 3b). We used the “metafor” package (version 3.4-0)68 in R for the linear regression analysis between MAP and RR(R/S) (Fig. 4). Structural equation model (SEM) was performed using the “lavaan” package (version 0.6-12)69 in R to examine the effects of PFTs, MFTs, background climate (MAT, MAP), and warming treatment (WM and DUR) on RR(R/S) via changes in AGB and BGB under warming (Fig. 6). The figures to depict correlations between specific response variables were drawn in SigmaPlot 12.5 (Systat Software, San Jose, CA, Supplementary Figs. 2 and 3).
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
Code availability
The code used in this study is available at https://figshare.com/articles/online_resource/Code1_2_txt/20390592/1.
References
Henneron, L., Cros, C., Picon-Cochard, C., Rahimian, V. & Fontaine, S. Plant economic strategies of grassland species control soil carbon dynamics through rhizodeposition. J. Ecol. 108, 528–545 (2020).
Arft, A. M. et al. Responses of tundra plants to experimental warming: meta-analysis of the international tundra experiment. Ecol. Monogr. 69, 491–511 (1999).
Bloom, A. A., Exbrayat, J.-F., van der Velde, I. R., Feng, L. & Williams, M. The decadal state of the terrestrial carbon cycle: global retrievals of terrestrial carbon allocation, pools, and residence times. Proc. Natl Acad. Sci. USA 113, 1285–1290 (2016).
Ma, H. et al. The global distribution and environmental drivers of aboveground versus belowground plant biomass. Nat. Ecol. Evol. 5, 1110-+ (2021).
Mokany, K., Raison, R. J. & Prokushkin, A. S. Critical analysis of root: shoot ratios in terrestrial biomes. Glob. Change Biol. 12, 84–96 (2006).
Shipley, B. & Meziane, D. The balanced-growth hypothesis and the allometry of leaf and root biomass allocation. Funct. Ecol. 16, 326–331 (2002).
Eziz, A. et al. Drought effect on plant biomass allocation: a meta‐analysis. Ecol. Evolution 7, 11002–11010 (2017).
Yan, Z. et al. Biomass allocation in response to nitrogen and phosphorus availability: Insight from experimental manipulations of Arabidopsis thaliana. Front. Plant Sci. 10, 598 (2019).
Li, C. et al. Precipitation and nitrogen addition enhance biomass allocation to aboveground in an alpine steppe. Ecol. Evol. 9, 12193–12201 (2019).
Piao, S. et al. Net carbon dioxide losses of northern ecosystems in response to autumn warming. Nature 451, 49–52 (2008).
Kim, J.-S. et al. Reduced North American terrestrial primary productivity linked to anomalous Arctic warming. Nat. Geosci. 10, 572–576 (2017).
Lin, D., Xia, J. & Wan, S. Climate warming and biomass accumulation of terrestrial plants: a meta-analysis. N. Phytol. 188, 187–198 (2010).
Fernandez, C. W. et al. Ectomycorrhizal fungal response to warming is linked to poor host performance at the boreal-temperate ecotone. Glob. Change Biol. 23, 1598–1609 (2017).
Keller, J. A. & Shea, K. Warming and shifting phenology accelerate an invasive plant life cycle. Ecology 102, e03219 (2020).
Cavagnaro, R. A., Oyarzabal, M., Oesterheld, M. & Grimoldi, A. A. Screening of biomass production of cultivated forage grasses in response to mycorrhizal symbiosis under nutritional deficit conditions. Grassl. Sci. 60, 178–184 (2014).
Rasheed, M. U. et al. The responses of shoot-root-rhizosphere continuum to simultaneous fertilizer addition, warming, ozone and herbivory in young Scots pine seedlings in a high latitude field experiment. Soil Biol. Biochem. 114, 279–294 (2017).
Xu, M., Liu, M., Xue, X. & Zhai, D. Warming effects on plant biomass allocation and correlations with the soil environment in an alpine meadow, China. J. Arid Land 8, 773–786 (2016).
Zhou, X., Talley, M. & Luo, Y. Biomass, litter, and soil respiration along a precipitation gradient in southern great plains, USA. Ecosystems 12, 1369–1380 (2009).
Hertel, D., Strecker, T., Mueller-Haubold, H. & Leuschner, C. Fine root biomass and dynamics in beech forests across a precipitation gradient - is optimal resource partitioning theory applicable to water-limited mature trees? J. Ecol. 101, 1183–1200 (2013).
Zhou, L. et al. Responses of biomass allocation to multi-factor global change: a global synthesis. Agriculture, Ecosyst. Environ. 304, 107115 (2020).
Poorter, H. et al. Biomass allocation to leaves, stems and roots: meta-analyses of interspecific variation and environmental control. N. Phytol. 193, 30–50 (2012).
Gorka, S., Dietrich, M., Mayerhofer, W., Gabriel, R. & Kaiser, C. Rapid transfer of plant photosynthates to soil bacteria via ectomycorrhizal hyphae and its interaction with nitrogen availability. Front. Microbiol. 10, 168 (2019).
Chen, W. et al. Root morphology and mycorrhizal symbioses together shape nutrient foraging strategies of temperate trees. Proc. Natl Acad. Sci. USA 113, 8741–8746 (2016).
Terrer, C. et al. A trade-off between plant and soil carbon storage under elevated CO2. Nature 591, 599–603 (2021).
Averill, C., Dietze, M. C. & Bhatnagar, J. M. Continental-scale nitrogen pollution is shifting forest mycorrhizal associations and soil carbon stocks. Glob. Change Biol. 24, 4544–4553 (2018).
Cheng, L. et al. Mycorrhizal fungi and roots are complementary in foraging within nutrient patches. Ecology 97, 2815–2823 (2016).
Averill, C. & Hawkes, C. V. Ectomycorrhizal fungi slow soil carbon cycling. Ecol. Lett. 19, 937–947 (2016).
Hollister, R. D. & Flaherty, K. J. Above- and below-ground plant biomass response to experimental warming in northern Alaska. Appl. Vegetation Sci. 13, 378–387 (2010).
Johnson, N. C., Rowland, D. L., Corkidi, L. & Allen, E. B. Plant winners and losers during grassland N-eutrophication differ in biomass allocation and mycorrhizas. Ecology 89, 2868–2878 (2008).
Xia, J., Yuan, W., Wang, Y. P. & Zhang, Q. Adaptive carbon allocation by plants enhances the terrestrial carbon sink. Sci. Rep. 7, 3341 (2017).
Litton, C. M. & Giardina, C. P. Below-ground carbon flux and partitioning: global patterns and response to temperature. Funct. Ecol. 22, 941–954 (2008).
Wang, P. et al. Belowground plant biomass allocation in tundra ecosystems and its relationship with temperature. Environ. Res. Lett. 11, 055003 (2016).
Hovenden, M. J. et al. Warming and elevated CO2 affect the relationship between seed mass, germinability and seedling growth in Austrodanthonia caespitosa, a dominant Australian grass. Glob. Change Biol. 14, 1633–1641 (2008).
Olszyk, D. M. et al. Whole-seedling biomass allocation, leaf area, and tissue chemistry for Douglas-fir exposed to elevated CO2 and temperature for 4 years. Can. J. For. Res. 33, 269–278 (2003).
Parmesan, C. & Hanley, M. E. Plants and climate change: complexities and surprises. Ann. Bot. 116, 849–864 (2015).
Hagedorn, F., Gavazov, K. & Alexander, J. M. Above- and belowground linkages shape responses of mountain vegetation to climate change. Science 365, 1119-+ (2019).
Pinto, R. S. & Reynolds, M. P. Common genetic basis for canopy temperature depression under heat and drought stress associated with optimized root distribution in bread wheat. Theor. Appl. Genet. 128, 575–585 (2015).
Rasse, D. P., Rumpel, C. & Dignac, M. F. Is soil carbon mostly root carbon? Mechanisms for a specific stabilisation. Plant Soil 269, 341–356 (2005).
Phillips, R. P., Brzostek, E. & Midgley, M. G. The mycorrhizal-associated nutrient economy: a new framework for predicting carbon-nutrient couplings in temperate forests. N. Phytol. 199, 41–51 (2013).
Wieder, W. R., Cleveland, C. C., Smith, W. K. & Todd-Brown, K. Future productivity and carbon storage limited by terrestrial nutrient availability. Nat. Geosci. 8, 441–444 (2015).
Wang, P., Huang, K. & Hu, S. Distinct fine‐root responses to precipitation changes in herbaceous and woody plants: a meta‐analysis. N. Phytol. 225, 1491–1499 (2020).
Ma, Z. et al. Evolutionary history resolves global organization of root functional traits. Nature 555, 94–97 (2018).
Prieto, I., Armas, C. & Pugnaire, F. I. Water release through plant roots: new insights into its consequences at the plant and ecosystem level. N. Phytol. 193, 830–841 (2012).
Bai, W. et al. Increased temperature and precipitation interact to affect root production, mortality, and turnover in a temperate steppe: implications for ecosystem C cycling. Glob. Change Biol. 16, 1306–1316 (2010).
Soudzilovskaia, N. A. et al. Global mycorrhizal plant distribution linked to terrestrial carbon stocks. Nat. Commun. 10, 5077 (2019).
Turner, B. L. Resource partitioning for soil phosphorus: a hypothesis. J. Ecol. 96, 698–702 (2008).
Phillips, L. A., Ward, V. & Jones, M. D. Ectomycorrhizal fungi contribute to soil organic matter cycling in sub-boreal forests. ISME J. 8, 699–713 (2014).
Gonzalez-Meler, M. A., Silva, L. B. C., Dias-De-Oliveira, E., Flower, C. E. & Martinez, C. A. Experimental air warming of a stylosanthes capitata, vogel dominated tropical pasture affects soil respiration and nitrogen dynamics. Front. Plant Sci. 8, 46 (2017).
Carrillo, Y., Pendall, E., Dijkstra, F. A., Morgan, J. A. & Newcomb, J. M. Response of soil organic matter pools to elevated CO2 and warming in a semi-arid grassland. Plant Soil 347, 339–350 (2011).
An, J. et al. Physiological and growth responses to experimental warming in first-year seedlings of deciduous tree species. Turkish J. Agriculture Forestry 41, 175–182 (2017).
Steidinger, B. S. et al. Climatic controls of decomposition drive the global biogeography of forest-tree symbioses. Nature 569, 404–408 (2019).
Liu, R., Li, Y., Wang, Y., Ma, J. & Cieraad, E. Variation of water use efficiency across seasons and years: Different role of herbaceous plants in desert ecosystem. Sci. Total Environ. 647, 827–835 (2018).
Duarte, A. G. & Maherali, H. A meta-analysis of the effects of climate change on the mutualism between plants and arbuscular mycorrhizal fungi. Ecol. Evol. 12, https://doi.org/10.1002/ece3.8518 (2022).
Bastos, A. & Fleischer, K. Fungi are key to CO2 response of soil. Nature 591, 532–534 (2021).
Wang, X., Peng, L. & Jin, Z. Effects of AMF inoculation on growth and photosynthetic physiological characteristics of Sinocalycanthus chinensis under conditions of simulated warming. Acta Ecologica Sin. 36, 5204–5214 (2016).
Tedersoo, L., Bahram, M. & Zobel, M. How mycorrhizal associations drive plant population and community biology. Science 367, 6480 (2020).
Jing, X. et al. The links between ecosystem multifunctionality and above- and belowground biodiversity are mediated by climate. Nat. Commun. 6, 8159–8159 (2015).
Pecl, G. T. et al. Biodiversity redistribution under climate change: Impacts on ecosystems and human well-being. Science 355, eaai9214 (2017).
IPCC. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change 1535 (Cambridge University Press, 2021).
Fick, S. E. & Hijmans, R. J. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315 (2017).
Task, G. Global Gridded Surfaces of Selected Soil Characteristics (IGBP-DIS) (2000).
Hedges, L. V., Gurevitch, J. & Curtis, P. S. The meta-analysis of response ratios in experimental ecology. Ecology 80, 1150–1156 (1999).
Luo, Y. Q., Hui, D. F. & Zhang, D. Q. Elevated CO2 stimulates net accumulations of carbon and nitrogen in land ecosystems: A meta-analysis. Ecology 87, 53–63 (2006).
Rosenberg, M. S., Adams, D. C. & Gurevitch, J. MetaWin: Statistical Software for Meta-analysis (Sinauer Associates, Incorporated, 2000).
Kembel, S. W. et al. Picante: integrating phylogenies and ecology. Bioinformatics 26, 1463–1464 (2018).
Calcagno, V. & De, C. M. Glmulti: an R package for easy automated model selection with (generalized) linear models. J. Stat. Softw. 34, https://doi.org/10.18637/jss.v034.i12 (2010).
Pinheiro, J. C., Bates, D. J., Debroy, S. D. & Sakar, D. nlme: Linear and nonlinear mixed effects models. R. package version 3, 1–117 (2009).
Viechtbauer, W. Metafor: meta-analysis package for R. J. Stat. Softw. 2010, 1–10 (2010).
Rosseel, Y. Lavaan: an R package for structural equation modeling. J. Stat. Softw. 48, https://doi.org/10.18637/jss.v048.i02 (2012).
Acknowledgements
We thank Pryianka Ambhavane (https://ambypriya.myportfolio.com/) for her help in drawing Fig. 7. X.Z. is supported by the National Natural Science Foundation of China (Grant No. 31930072). L.Z. is supported by the National Natural Science Foundation of China (Grant No. 32071593, 31600352), the National key Research and Development Program of China (Grant No. 2020YFA0608403), and the Fundamental Research Funds for the Central Universities. M.P.T. is supported by the Swiss State Secretariat for Education, Research, and lnnovation (SERI) under contract number M822.00029.
Author information
Authors and Affiliations
Contributions
L.Z. collected and analyzed the data and wrote the manuscript with substantial contributions from M.P.T. and X.Z. X.Z. conceived, designed, and oversaw the study. L.Z. and Y.H. did the statistical analysis with suggestions from M.P.T. Y.H., Y.F., Z.D., M.L., X.S., C.L., C.L., R.L., G.Z., and S.H.B. commented on the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zhou, L., Zhou, X., He, Y. et al. Global systematic review with meta-analysis shows that warming effects on terrestrial plant biomass allocation are influenced by precipitation and mycorrhizal association. Nat Commun 13, 4914 (2022). https://doi.org/10.1038/s41467-022-32671-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-022-32671-9