Nutrient resorption strategies of three oak tree species in response to interannual climate variability

Nutrient resorption is critical for plants toward balancing their nutritional requirements and adapting to environmental variabilities, which further impacts litter quality and nutrient cycling. However, the interannual variability of nutrient resorption under climate change remains unclear. We investigated the five-year nutrient resorption efficiencies (NuRE, %) of 14 elements in three deciduous oak tree species (Quercus aliena var. acuteserrata, Q. glandulifera, and Q. variabilis) in a warm-temperate forest of Central China and assessed their relationships with interannual climate and soil factors. Nutrient resorption did not differ between species but varied significantly between different years. For each year, N, P, S, K, C, Mg, and Zn were preferentially resorbed in all of the oak species in contrast to Ca, Na, Mn, Ba, Al, Fe, Cu, which were to some extent discriminated. Among the 14 elements, the NuRE of C, N, P, S, Ca, and Mg was more sensitive to interannual climate variations in the three oak species. The carbon resorption efficiency was significantly increased during the driest year of the study (2014); N resorption efficiency was reduced with temperature; whereas N and P resorption efficiency initially decreased and then increased with precipitation. Moreover, the elements with higher NuREs typically had lower coefficient of variation (CV) in all three oak species. Different oak species exhibited analogous nutrient conservation strategies in response to annual climate variabilities, and interannual climate variations strongly impacted plant nutrient resorption. Deciduous plants may establish a tradeoff mechanism to rebalance somatic nutrients for regrowth at the end of the growing season.


Background
Plants obtain the nutrients they require, on one hand, through soil nutrient cycling, which includes inputs from litter decomposition, biological fixation, atmospheric deposition, weathering, removal by root uptake, and leaching losses (Pallardy 2010;Cronan 2017), and on the other hand through the reallocation of nutrients to different structures (Killingbeck 2004;Marschner and Marschner 2012). Nutrient resorption is a process through which plants reallocate nutrients from senescing structures to storage and/or developing tissues, which contributes to a reduction in their dependence on nutrient uptake (Aerts and Chapin 1999;Killingbeck 2004;Brant and Chen 2015). Nutrient resorption can assist plants with adapting to environmental variabilities (May and Killingbeck 1992;Killingbeck 1993;Prieto and Querejeta 2020), which can further affect plant reproduction, litter decomposition, and ecosystem nutrient cycling through altering nutrient storage and litter quality (Attiwill and Adams 1993;Aerts 1996;. Therefore, the study of nutrient resorption under climate change has gained increasing attention in recent years (Fujita et al. 2014;Yan et al. 2016;Drenovsky et al. 2019).
Nutrient resorption varies substantially with temperature and precipitation along large geographical gradients (Kobe et al. 2005;Yuan and Chen 2009;Reed et al. 2012;Vergutz et al. 2012;Du et al. 2017;Tong et al. 2021). Climate change issues including global warming, extreme drought, and N deposition also influence nutrient resorption (Suseela et al. 2015;Zheng et al. 2018;Huang et al. 2018;Prieto and Querejeta 2020;Lü et al. 2021). For instance, climatic warming may reduce nutrient resorption (Suseela et al. 2015), while drought can advance leaf senescence, which increases (Lobo-do-Vale et al. 2019) or reduces (Zhao et al. 2020) resorption. Aside from climatic factors, nutrient resorption is also determined through the availability of soil nutrients, plant nutrient status, and functional groups (Chapin and Moilanen 1991;Wright and Westoby 2003;Liu et al. 2014;Bahamonde et al. 2019;Gerdol et al. 2019;, as well as source-sink relations (Killingbeck 2004), and the mobility of nutrients in the phloem (Bukovac and Wittwer 1957;Hill 1980). Since most species are unlikely to migrate quickly enough to keep pace with environmental changes (IPCC 2014), the capacity to adapt to variations in localized environments might assist them with surviving under the altered climatic environments that are anticipated in the near future (Jump and Peñuelas 2005; Prieto and Querejeta 2020). However, it is largely unknown how local plants might respond to disparities in environmental factors that are driven by interannual climate variations, which is critical to elucidate the physiological and evolutionary adaptation of plants under climate change.
The nutrient resorption levels of plants are associated with their constituent elements (Liu et al. 2014). Generally, N, P, S, and K are preferentially withdrawn from senescing leaves and exhibit high resorption, while Ca, Mg, Al, and Fe are accreted in senesced leaves, thus, showing low or negative resorption (Killingbeck 1993;Liu et al. 2014;Estiarte and Peñuelas 2015;Du et al. 2017;Prieto and Querejeta 2020). Various elements are often distributed in a coordinated pattern during the nutrient resorption process in a tradeoff manner (Killingbeck 1993(Killingbeck , 2004, i.e. some elements are taken up preferentially (NuRE ≥ 0) and others are discriminated (NuRE < 0) (Du et al. 2017). Further, to adapt to environmental changes, plants have also evolved a homeostasis mechanism that is based on the stability of the limiting elements hypothesis, i.e. with higher elemental concentrations and lower variability (low coefficient of variation) (Han et al. 2011). However, less is known about how different nutrients interact during the resorption process. Moreover, plants growing at the same sites may exhibit parallel or convergent adaptations in terms of functional traits, physiology, and genetics (Aspinwall et al. 2016;Yeaman et al. 2016;Bauters et al. 2017;Bahamonde et al. 2019). Whether different plants also exhibit analogous responses in nutrient resorption remains unclear. These findings will assist ecologists with optimizing their biogeochemical models, while improving their understanding of nutrient cycling.
The genus Quercus is widely distributed in temperate forests of the Northern Hemisphere and provides important ecological and economic roles (Aldrich and Cavender-Bares 2011). In our study area, three Quercus species emerged to dominate an area. Previous studies have focused primarily on the impacts of warming and drought on soil carbon dynamics , spatiotemporal variations in plant stoichiometry, and water use efficiencies in different oak forests (Du et al. 2017(Du et al. , 2021), yet, less is known about the effects of annual climatic variations on nutrient resorption. For this study, we investigated the interannual variations in multiple element stoichiometry and resorption for three deciduous oak species over five years. Our objectives were to examine the nutritional responses of different plants to climatic differences driven by interannual variability, and to assess the regulation patterns and driving factors of nutrient resorption based on climate change and multiple additional elements. To achieve our objectives, we measured the NuRE variations of seven macroelements (C, N, P, S, K, Ca, and Mg), four microelements (Fe, Mn, Zn, and Cu), and three trace elements (Na, Al, and Ba) (CRE referenced Vergutz et al. 2012). We evaluated the effects of interannual temperature, precipitation, as well as green leaf concentration and soil nutrient status, on the nutrient resorption of multiple elements.

Study area
This study was conducted from 2013 to 2016, and 2018 at the Forest Ecological Research Station of the Baotianman Natural Reserve (111°47′-112°04′ E, 33°20′-33°36′ N, 500 to 1,845 m a.s.l) in Henan Province, Central China. The missing study in 2017 was caused by accidental damage to the sampling equipment. This area belongs to a temperate-subtropical ecotone zone with mean annual precipitation of 936 mm, and air temperature of 12.4°C over the past thirty years (from 1981 to 2010) (http://climateap.net/). With the main soil type comprised of haplic luvisol , the forest primarily consists of deciduous oak, including Q. aliena var. acuteserrata, Q. glandulifera, and Q. variabilis, in addition to several minor components of other deciduous woodland species and shrubs.

Soil and leaf sampling
Three study sites were selected, where each site represented a stand of a subject tree species. For each site, three 20 m × 20 m plots were demarcated and the latitude, longitude, and altitude were recorded with a global positioning system (Thales Navigation, Santa Barbara, CA, USA) ( Table 1). Within each plot, five dominant trees were selected for the collection of leaf samples in the upper central south-facing part of crown. Two hundred leaves were collected in August (green leaves) and again in December (senesced leaves). The green and senesced leaves were obtained by cutting the upper middle branches with a retrofitted fishing rod on the ground. A total of five soil cores (0-10 cm) were randomly extracted using a 2-cm diameter stainless-steel earth borer in each plot, which were combined to produce a composite sample. The soil samples were airdried, while the green leaf and senesced leaf samples both were oven dried (65°C) for 72 h. The soil and leaf samples were ground and sieved through a 60-mesh sieve (0.25 mm diameter) for chemical analysis.
Nutrient resorption was expressed as the difference in the elemental concentrations between the green (X gr ) and senescent (X sen ) leaves by dry mass.

Nutrient resorption efficiency
where MLCF is the mass loss correction factor and~0.784 for deciduous broadleaved species (Vergutz et al. 2012).

Chemical analyses
The total C and N concentrations of the leaves and soil were determined using elemental analysis (Vario EL cube, Elementar, Germany) (Sun et al. 2012), whereas the total P, S, K, Na, Ca, Mg, Al, Fe, Mn, Zn, Cu, and Ba were assessed using an inductively coupled plasma optical emission spectrometer (ICP-OES) (Iris Advantage 1000, Thermo Jarrell Ash, Franklin, MA, USA) following digestion with hydrofluoric and nitrohydrochloric acid for the soil samples, and nitric and perchloric acid for the leaf samples.

Data analyses
Soil element concentrations of the three sampled oak species sites are summarized in Table S1. Meteorological data was obtained from the Neixiang Meteorological Bureau and Baotianman Forest Ecological Research Station ( Fig. 1). One-way ANOVA was conducted to test the differences in nutrient resorption between the different years.
One-way repeated measures ANOVA was employed to determine whether resorption of nutrients and their interacting effects differed between species and years. Cluster analysis was run to group elements by NuRE. Subsequently, linear regressions were employed to examine the nutrient resorption correlations with MAT, MAP, and the concentrations of elements in the green leaves.   Hierarchical partitioning (HP) analysis was performed to test the combined effects of MAT, MAP, as well as leaf and soil element concentrations on nutrient resorption efficiencies. One-way ANOVA, one-way repeated measures ANOVA, and cluster analyses were performed using SPSS 16.0 (IBM, Chicago, IL, USA). The HP analysis were carried using R version 3.6.1 (R Development Core Team, accessed 2019). The linear regressions and all figures were conducted using SigmaPlot 10.0 (Systat Software, Inc., Richmond, CA, USA).

Interannual and interspecies variations in nutrient resorption efficiencies
For all elements examined, NuRE generally followed the order of P, K, N, S, Mg, C, Zn, Mn, Na, Ca, Ba, Fe, Cu, and Al (from high to low) for the three oak species. Generally, the NuRE of P, K, N, S, Mg, C, and Zn was greater than zero, while that of Mn, Na, Ca, Ba, Fe, Cu, and Al was less than zero (Table 2). In the three oak species, with the exception of Zn, the NuREs of other elements did not exhibit significant differences, and over the five years of the study the NuREs of all elements (except for K, Mn, and Ba) showed significant differences (Table 3). Cluster analysis indicated five similar clusters in the three oak species (Fig. 2). Canonical discriminant analysis further indicated that nutrient resorption did not differ between species, but varied significantly between years (P < 0.001) (Fig. S1). For all three species, the coefficient of variation (CV, %) of elements between the different years generally decreased with higher NuREs (Fig. 3, Table S2).

Effects of climatic and nutritional conditions on plant nutrient resorption
During the driest year of the study (2014), CRE was significantly enhanced in the oak trees. It was observed that NRE and SRE were generally reduced with higher MAT, while CaRE and MgRE initially increased and then decreased with MAT (Fig. 4). The CRE and CuRE were typically decreased with higher MAP. Additionally, NRE, PRE, and SRE initially decreased and then increased with MAP, whereas CaRE and MgRE first increased and then decreased with MAP (Fig. 5). The NuREs of most elements were positively correlated with the elemental concentrations of green leaves; however, PRE was negatively correlated with the P concentrations of green leaves (Fig. S2). Hierarchical partitioning analysis indicated temporal variations in NuRE over the years under study, from 18.13 % KRE to 96.04 % PRE in the three oak species (Table S3). The NRE was strongly associated with MAT and leaf N, CuRE with MAP in the three species, PRE with MAT and leaf P, MgRE with MAT in Q. aliena var. acuteserrata and Q. glandulifera, and CaRE, MnRE, and BaRE with MAT and leaf elements  in Q. glandulifera and Q. variabilis. The NuRE of the remaining elements was mostly related to leaf element concentrations, except for K, Na, and Cu. Soil elements did not influence NuRE, except for CRE in Q. glandulifera and ZnRE in Q. variabilis (Table S3).

Analogous response of interspecies NuREs to interannual climate variability
In contrast to deciduous plants on a global scale (Vergutz et al. 2012), the three oak species under study demonstrated higher N, P, and K resorption efficiencies, but lower C and Mg resorption efficiencies. The NuRE of N and P in this study was higher than that of deciduous plants in a subtropical karst region (Liu et al. 2014), and higher than that of Q. variabilis across China (Sun et al. 2015). These results signified that plants may have more N and P demand for growth in this region. This may have been because the plants grew at high altitudes where due to low temperatures, they required higher levels of nutrient resorption to store nutrients for the following growing season. This was demonstrated by Du et al. (2017), who found that the resorption of N, P, and S (at least in Q. variabilis) was the greatest at high altitudes. Moreover, PRE in the three oaks was greater than NRE (except in 2018). This was in alignment with the research of Chai et al. (2015) and Du et al. (2017), who revealed that in the same or neighboring regions, the growth of oaks was commonly restrained by P in contrast to N.
As a conservation mechanism, nutrient resorption plays a key role in balancing plant nutrient demands and acquisition (Killingbeck 2004;Brant and Chen 2015). Our results revealed that NuRE varied significantly among the different years but did not differ between species. Interestingly, the three oak species displayed analogous NuRE responses to interannual climate variability, particularly for C, N, P, S, Mg, Fe, Mn, and Cu (Table 2), and cluster analysis also indicated five similar clusters in the three oak species (Fig. 2). The similarity in NuRE between different plants indicated that the same plant types (only deciduous plants were investigated in the study) growing in the same region might have parallel or convergent mechanisms in terms of nutrient conservation (Yeaman et al. 2016). This aligned well with a study by Bahamonde et al. (2019), who found that the NuRE did not vary between Nothofagus pumilio and N. antarctica at the sites they investigated. Similar results were reported following a five year manipulative field study by Prieto and Querejeta (2020), and demonstrated by Lü et al. (2013), who found that with N inputs, plants exhibited a convergent response in the resportion of N and P. These results indicated that plants growing under the same conditions, or occupying the same niche, may experience analogous limitations for nutrient absorption and conservation processes. Thus, our findings may have important implications for the study of localized plant adaptations and evolution.  Table S2 In addition, our results showed that the NuRE of P, K, N, S, Mg, C, and Zn was greater than zero, while that of Mn, Na, Ca, Ba, Fe, Cu, and Al was less than zero in the three oak species (Table 2), which was demonstrated by the study of Du et al. (2017) for Q. variabilis. Macroelements (P, K, N, S, Mg, C, and Zn) that are constituents of amino acids, proteins, nucleic acids, and enzymes (Marschner and Marschner 2012) are essential for new growth, and therefore had a high resorption efficiency (except for Ca) (Liu et al. 2014;Du et al. 2017). Microelements and trace elements (Mn, Na, Ca, Ba, Fe, Cu, and Al) are closely related to enzymes, cell structure and detoxification. The negative resorption efficiencies of these elements have also been reported previously (Killingbeck 1993;Liu et al. 2014;Du et al. 2017). This is likely due to these nutrients being freely loaded into xylem but highly restricted in phloem (Bukovac and Wittwer 1957;Hill 1980), which limits their inputs into perennial tissues, and therefore had a low even negative resorption efficiency. As expected, seasonal variations in the elements of Q. variabilis of the same region also demonstrated low nutrients for microelements and trace elements and high nutrients for macroelements at the onset of the growing season (Du et al. 2021). These findings indicated that the oak species engaged in resorption trade-offs as the end of the growing season approached, which may affect plant regrowth. Furthermore, we found that the elemental CVs between different years generally decreased with higher NuREs in all three oak species (Fig. 3) as supported by the study of Karimi and Folt (2006). This was because high elemental concentrations with low variability benefitted plants by maintaining their homeostasis (Han et al. 2011).

Effects of interannual climate, and leaf and soil nutrient status on NuRE
Significant fluctuations in nutrient resorption due to interannual variabilities in temperature and precipitation were reported in two prior studies (May and Killingbeck 1992;Killingbeck 1993). Of the fourteen elements examined the NuRE of C, N, P, S, Ca, and Mg was closely associated with MAT and MAP (Figs. 4 and 5). The results suggested that these elements were more sensitive to interannual climate variability compared to others, which due to their importance translated to limitations in photosynthesis and growth in plants (Elser et al. 2007;Marschner and Fig. 4 Relationships between element resorption efficiencies and mean annual temperature Marschner 2012). This finding signified that the responses of resorption to climatic variability were contingent on element type (Vergutz et al. 2012;Sun et al. 2015;Du et al. 2017), where some elements that are directly related to growth (such as nucleic acid-proteins and photosynthetic elements) may be more sensitive to climate change.
There was negative relationship between NRE and temperature across all species (Fig. 4), which was demonstrated in numerous previous studies (Yuan and Chen 2009;Sun et al. 2015;Prieto and Querejeta 2020). This was primarily because colder temperatures not only inhibited the metabolic activity of plants, but also restricted nutrient release and root N uptake. Consequently, plants need to extract additional N to enhance metabolic activity and growth rates, while reducing their dependence on available soil nutrients. This result also uncovered the interesting finding that temperature had a consistent influence on NRE regardless of spatial or temporal scales. However, precipitation had more complex influence on nutrient resorption. For instance, NRE and PRE initially decreased and then increased with MAP (Fig. 5). Comparatively, Yuan and Chen (2009) found that P resorption increased with higher precipitation, while Vergutz et al. (2012) reported that N and P resorption decreased with precipitation on a global scale. Moreover, Du et al. (2017) found that N and P resorption was enhanced with precipitation along an altitude gradient in the same area. Our study indicated that excessive or low precipitation may have restricted the growth of plants, which resulted in high N and P resorption. These results revealed that precipitation may play a predominant role in the determination of nutrient demands and plant growth in different oak forests.
During the driest year of the study (2014) CRE was significantly enhanced in the oak trees (Table 2), which suggested that plants may store additional energy for reproduction under drought conditions (Vergutz et al. 2012). Correspondingly, the concentrations of C in the litter of the forest floor were reduced (Fig. S3), as was demonstrated by Sardans et al. (2008). This might influence future litter quality and decomposition, and even ecosystem C cycling. Based on the predicted changes in local temperature and precipitation (0.9°C higher annual mean temperature and 6.4 mm lower total precipitation over the next twenty years) (Zhao et al. 2012), the resorption efficiency of C, N, and P will likely (to some extent) decrease over the next twenty years (Table S4), as demonstrated by Prieto and Querejeta (2020). This decrease likely reduced nutrient demands and plant growth and may have some influence on litter quality and decomposition (Prescott 2010; Prieto et al. 2019).
NuRE was also associated with the elemental concentrations of green leaves. Surprisingly, there were positive relationships between the concentrations of leaf element and NuREs (except for P) (Fig. S2). Our results were different from previous studies, where the high nutrient contents of green leaves typically led to low nutrient resorption (Kobe et al. 2005;Ratnam et al. 2008;Vergutz et al. 2012), or there was no significant relationship between them (Chapin and Moilanen 1991;Aerts 1996). For our study, the relationships between the elemental concentrations of green leaves and climatic factors were also analyzed, which were consistent with the relationships between NuRE and climatic factors. This revealed that both the elemental concentrations of green leaves and the resorption of nutrients were impacted by climatic factors, which may form coupling relationships. When plants possessed high nutrient concentrations during the growing season, increased nutrient resorption was likely to occur at the end of the growing season to maintain a source-sink balance (Killingbeck 2004). Furthermore, we found that soil had no influence on NuRE during the years under study (Table S3), which agreed with the findings of Aerts (1996) and Aerts and Chapin (1999). This was likely caused by the low fluctuation of available soil nutrients between years and low sensitivity to changes in the soil due to the long-term adaptation to the local environment.

Conclusions
Our study described a comprehensive representation relating to the interannual patterns of nutrient resorption of multiple elements in three oak species. Our conclusions were that (1) different plants growing in the same region may have analogous nutrient conservation strategies in response to interannual climate change; (2) at the conclusion of the growing season deciduous plants may engage a tradeoff mechanism to rebalance nutrients for regrowth; (3) interannual climate variations strongly impacted plant nutrient resorption, particularly, C, N, P, S, Ca, and Mg, which play important roles in adapting to interannual climate variabilities and reproduction. Our findings have clear implications for understanding localized plant adaptation strategies and biogeochemical cycling processes under global climate change.