Stoichiometry of litter decomposition under the effects of climate change and nutrient enrichment: A meta‑analysis

Aims The cycling of nutrients from plant litter has key implications for the functioning of terrestrial eco-systems by controlling nutrient availability and net primary production. Despite extensive research on the effects of global change on ecosystem functioning, the direct implications of global change on stoichiometry and nutrient dynamics during litter decomposition remain poorly understood. To address this gap, we conducted a meta-analysis. Methods We analysed 178 experiments that simulated (i) warming, (ii) drought, (iii) increased water availability, (iv) N enrichment, (v) P enrichment, and (vi) combined N and P (N + P) enrichment. We compared earlier (approximately six months) and later (approximately one year) stages of decomposition and analysed the specific effects taking into account climate and plant type. Results The C:N and C:P ratios decreased in most warming and nutrient enrichment scenarios, leading to losses of litter C content, while the N:P ratio remained more resilient and affected by water availability. Furthermore, the abundance of resources (water and N + P) fosters the decomposition of litter. The nutrient mobilisation increases for both P and N under non-limited nutrient enrichment and it is faster for N than for P when water increases its availability. Nutrient enrichment was relevant in later stages of decomposition. Conclusions Our study provides insights into the fate of litter decomposition and its stoichiometric dynamics in response to drivers of global change. Concerning scenarios of C release and N and P immobilisation were identified. However, further experimentation and analysis are necessary to consider all interacting drivers.


Introduction
Nutrient stoichiometry is a crucial factor in controlling important ecological processes, such as the growth rate Responsible Editor: Ivika Ostonen.and life strategies of organisms, as well as ecosystem diversity, structure, and functionality (Sterner and Elser 2002;Sardans et al. 2012b;Peñuelas et al. 2019).In particular, litter decomposition is one of the key processes that determines ecosystem functioning (Zhang et al. 2018;Lanuza et al. 2019).Currently, the biosphere is facing various transformations, including climate-related changes like global warming and shifts in precipitation patterns (IPCC 2021), as well as changes in the availability of elements such as C and N due to anthropogenic activities (Gruber and Galloway 2008;Sundquist 1993).However, the availability of other essential nutrients like P and K are not increasing to the same extent (Sardans et al. 2012a(Sardans et al. , 2021;;Peñuelas et al. 2020).These changes in nutrient stoichiometry and imbalances are altering the "elementome" (i.e., the elemental composition of ecosystems) of soils (Peñuelas et al. 2019;Fernández-Martínez 2022;Peñuelas and Sardans 2022).Understanding how these shifts in elementome, together with climate change, affect nutrient cycles in plant-soil systems is crucial for the development of global change models, conservation strategies, and food security policies.
Litter decomposition is an important process for understanding plant-soil nutrient cycles, as litter serves as a crucial source of mineralizable nutrients for microbial and plant metabolism in terrestrial ecosystems (Silveira et al. 2011), that controls the availability of nutrients in soils, cation-exchange capacity, soil pH, and the formation of humus (Berg and McClaugherty 2008).The release of nutrients during litter decomposition allows microbes and plants to take them up.Litter decomposition also plays a significant role at the biospheric level, affecting the exchange of greenhouse gases and the carbon stocks in soil (Walkiewicz et al. 2021), which can, in turn, affect climate change (Parker et al. 2018).Drivers of global change such as increasing temperatures and precipitation can accelerate the rate of litter decomposition (Zhang et al. 2008), while droughts can slow it down and reduce the C content of litter (Deng et al. 2021).Additionally, nutrient enrichment can have diverse effects on litter decomposition rates, increasing them in grasslands and wetlands but decreasing them in forests (Su et al. 2022).
Litter decomposition and nutrient mineralisation also depends on the species that compose the vegetation cover, or the environmental factors, like climate, that characterise the ecosystem; and will determine the shifting of elementomes in soils.Temperature and precipitation can affect litter decomposition directly by regulating the activity of decomposers, and indirectly through changes in plant species composition and abundance, as a result, in variations in litter quality and quantity (Petraglia et al. 2019).Colder and drier climates present higher stabilisation of organic material (Mueller et al. 2018), while warmer and wetter climates are more likely to promote faster decomposition by the soil food webs (García-Palacios et al. 2013).On the other hand, litter quality is considered a key regulator of the decomposition dynamics in litter (Zhang et al. 2008).Decomposing plant litter from broadleaved forest beds of grass species is less constrained than conifer tree litter, whose lignin content and allelopathic substances are generally higher (Rahman et al. 2013;Liu et al. 2019).The effect of global change drivers then differ depending on the litter quality.N additions had been reported accelerating decomposition in low-lignin litters, while this effect was inhibiting in high-ligning litters (Hobbie et al. 2012;Knorr et al. 2005).Another aspect to consider is the stage of the litter decomposition, with earlier stages characterised by the decomposition of non-lignified cellulose and hemicellulose; and the leaching of soluble compounds; while later stages encompasses the degradation of the lignified tissues, requiring a more complex metabolism (Djukic et al. 2018).
Litter decomposition has been a widely studied field, even in relation to global change analyses (Zhang et al. 2008;Deng et al. 2021;Su et al. 2022).However, there is still a lack of an all-encompassing understanding of how litter stoichiometry during the decomposition process varies in relation to climate change and nutrient enrichment.Most review studies focus on one aspect of litter dynamics, such as litter decomposition rate (Zhang et al. 2008;Su et al. 2022) or the effect on a particular elemental content (Deng et al. 2021), and others, are barely focused on the dynamics of litter decomposition (Sardans et al. 2012a).Our objective with this analysis is to examine how nutrient addition and climate change, two important indicators of global change, affect the stoichiometry of litter decomposition, and to detect how it differs depending on and climate and the plant group producing the litter through a global assessment.To achieve this goal, we gathered data from experiments measuring the stoichiometric relationships of at least two elements (C, N, P; K was also contempled) in decomposing litter and where treatments simulated one of four drivers of global change: (i) warming, (ii) drought, (iii) irrigation (an increase in water availability), or (iv) N and P enrichment (separately or combined).We evaluated nutrient dynamics in the earlier (approximately six months) and later (approximately one year) stages of decomposition and also considered factors such as climate and plant type.By focusing only on studies that examined the simultaneous response of at least two elements, we reduced the number of studies available for the meta-analysis, but we were able to investigate elemental ratios.Additionally, we examined the effects on the decomposition rate and individual contents of C, N, and P.
Global change affects litter decomposition dynamics: an increment in the decomposition rate results in a preferential mineralisation of C to N and P; and N to P (Sardans et al. 2012a).When comparing experimental and control plots, we expect to find the following: (i) Increases in temperature will reduce the C:N, C:P and N:P ratios as decomposer activity is enhanced; (ii) Water scarcity will elevate C:N, C:P and N:P ratios, as decomposition is hampered, but the opposite effect is expected if water availability increases; (iii) Nutrient enrichments would also increase decomposition activity, leading to reductions of C:N and C:P ratios.However, the N:P ratio would differ depending on the enrichment -increasing if the enrichment in N dominates and decreasing if it is P dominated.This effect would be expected to be (iv) greater in earlier decomposition stages than in later, as the preferential mineralization of nutrients by decomposers may disappear with time when the carbon-rich -easy-to-decompose -litter component disappears.Additionally, we expect to find different responses to different climates and plant types.As an example, (v) we expect to find that increases of temperature and changes in precipitation are more likely to affect colder and dry climates compared to warm and wetter ones; and finally.(vi) we expect conifer, high-ligning, litter being less affected by nutrient enrichments compared to grasses of broadleaved species.

Data collection
To compile a dataset of the effects of global change drivers on litter stoichiometry, we conducted a bibliographic search on the 1st of February 2022 for peer-reviewed papers on Scopus, Web of Science, and Google Scholar.Three search strings were performed on title, keywords, and abstract: (i) "(litter AND decomposition) AND (warming OR (increase* AND temperature*)) AND experiment* AND ((N:P) OR (C:N) OR (C:P))" for global warming as an agent of change; (ii) "(litter AND decomposition) AND ((N OR P) AND (deposit* OR fertiliz* OR addition)) AND experiment* AND ((N:P) OR (C:N) OR (C:P))" for nutrient enrichment (e.g.fertilisation or enhanced deposition) as an agent of change, and (iii) "(litter AND decomposition) AND (precipitation OR rainfall OR drought) AND (reduc* OR addition OR increment* OR exclusion) AND experiment* AND ((N:P) OR (C:N) OR (C:P))" for changes in the availability of precipitation as an agent of global change.After screening the manuscripts, we selected studies that reported: (i) plots with experimental and control conditions, (ii) measurements of the variables in the experimental and control plots that were conducted under the same temporal and spatial scales, and where the magnitude of the treatment and experimental duration were clearly recorded and (iii) sufficientdata to calculate the C:N, C:P, and/or N:P ratios of the remaining litter.Our final database contained 43 studies (highlighted with an asterisk in the reference section; Table S1).All data for each study were obtainedfrom the associated published data sets provided by the author/s while extracting data from figures was done using WebPlotDigitizer (Rohatgi 2021).We considered different experiments when (i) the litter type differed, (ii) plots differed geographically within the same study, (iii) different treatments were applied independently, and (iv) plots varied in the intensity of treatment.A PRISMA formatted flow diagram (Moher et al. 2009) is provided in the supplementary material (Figure S1).
We recorded ancillary data for each study, including information for the study site such as geographic location of the sampling sites (latitude, longitude, and elevation), information for the litter species sampled, the litter sampling method (litterbag experiments vs litter collection), the methodological approach to measure each of the elements, and the initial litter biomass and elemental concentrations.The experimental conditions were characterised in three groups: (i) increase in temperature, the intensity of which was recorded in °C, (ii) enrichment of N, P, or both N and P, the intensity of which was recorded as g m −2 y −1 , and (iii) changes in the amount of precipitation, which was recorded as the percent decrease or increase in water availability relative to the control plots in a year.While we recorded this ancillary data from the paper, we also recorded the climate zone of the location (polar, cold, temperate, arid, or tropical) using the recorded longitude and latitude, and the kgc package of R (Bryant et al. 2017).If the plant type (angiosperm tree, gymnosperm tree, shrubs, and grasses; all herbaceous species were considered grasses) was not specified in the text, we filled the gaps using the Global Biodiversity Information Facility (GBIF 2021).We ultimately analysed six data sets: i) warming treatments (47 experiments), ii) enrichment with N only (55 experiments), iii) enrichment with P only (13 experiments), iv) enrichment with both N and P (38 experiments), vi) decreases in the amount of precipitation (drought treatment, 13 experiments), vii) increases in the amount of precipitation (irrigation treatment, 12 experiments).We rejected experiments that entailed more than one of the selected treatments (warming, irrigation, drought, N enrichment, P enrichment, and N + P enrichment) at the same time.After applying these requirements our database contained 178 experiments.The starting date and associated duration (days) for each experiment were compiled for 123 of the experiments.We then compiled both the percentage of remaning control and experimental litter masses (only for the litterbag experiments), the elemental concentrations, elemental contents (percentages of the initial quantity of the element remaining after the experiment) C, N and P; and the C:N and C:P litter ratios.The standard deviation and number of replicates for all these variables for both the control and treatment plots were also recorded.If provided, we also collected the K concentration.

Calculation of variables
Depending on the study, the elemental data were the concentration of the element, the content that remained relative to the initial value, and/or the ratios of the elements.We calculated the concentrations of C, N, P, and K, the contents of C, N, and P remaining relative to the initial values, and the C:N, C:P, and N:P ratios (on a mass basis, for the mass remaining) for all experiments when possible.We initially recorded all sampling dates for each experiment.First, we subsetted our database considering sampling dates resembling experiments that last for approximately half a year (180 d).This subset contained the samples representing the "earlier" stage of decomposition with a median of 233 ± 26 days (see below).To test changes in later decomposition stages (prediction iv), we subsetted the database filtering sampling dates closer to three years (1080 d), but the median resulted similar to one year: 412 ± 32 day.Thus, this database contained experiments where litter was at a "later" stage of decomposition.A detailed information about the days of experiments foe each study can be found in Table S2.We thus constructed a database for all experiments at an "earlier" stage of decomposition, and others at a "later" stage.Additionally, we calculated the decomposition rate based on the classical definition (Olson 1963): where k is the decomposition rate, t is the duration (in days) of the experiment, which was recorded for every experiment, and X 0 and X t are the initial litter mass and the litter mass remaining at the time of sampling, respectively (in grams).We also applied this formula to the initial and remaining litter elemental contents (for obtaining the rates of release of the elements), and to the initial and remaining litter elemental ratios (for obtaining the rates of change of the elemental ratios).
Before determining the stoichiometric effects of global change drivers, we assessed the effects of each treatment on the mass remaining at the end of the experiment and the decomposition rates of the litter biomasses.We determined the effects of the treatments on the litter concentrations of C, N, and P and calculated the elemental contents remaining (% from the initial) after the experiment (when the remaining contents were not provided).Due to the lower number of experiments providing initial litter masses compared to those providing elemental concentrations, the sample size used for the content analysis was smaller than that for the concentration analysis.Ratios were calculated as the quotient of the elemental contents, but we used the ratios provided when these data were not available.To calculate k (rate of change) associated with the shifting ratios, we needed the initial ratios at the start of the decomposition experiments.

The final data set
Our final data set contained 42 sampling locations, (the locations of the mesocosms were not considered) (Fig. 1); and data for litter from 71 species.
From the 178 experiments selected, and based on the Köppen-Geiger classification (Beck et al. 2018)

Statistical analysis
We used the metafor package (v2.4-0;Viechtbauer 2010) in R Studio (v4.1.2;R Core Team 2021) for the statistical analyses.We calculated the sizes of the effects of each treatment (warming, N enrichment, P enrichment, N + P enrichment, increase in water availability, and decrease in water availability) on the remaining litter biomass, the contents of C, N, P, and K, the release rate of these contents, the C:N, C:P, and N:P ratios, and the rates of change of these ratios.To calculate the effect size of a treatment, we calculated the standard mean difference (SMD; Hedges 1981) with the function escalc, where the numerator of SMD was set as the value of the variable under the treatment minus the value of the variable under the control (retaining negative signs).Thus, a positive treatment effect (SMD) indicates that the treatment group has, on average, a higher value than the control group.On the other hand, a negative treatment effect indicates that the treatment group has a lower value compared to the control group.We constructed meta-analytic models as random-effect models using the Restricted Maximum Likelihood method (REML; Cook 2014) with the rma.mv function.We rejected the meta-analytic models when the number of experiments was lower than 5, and those rejected models can be consulted in Table S3We calculated the general effects of the treatments over the experiments using the study as For the treatments that significantly affected the variables, we calculated the sensitivity of the effect size to the intensity of the treatment by including the intensity as a moderator.We also explored the effects of the treatments over the experiments taking into account factors such as climate (arid, cold, temperate, polar, or tropical; according to Köppen-Geiger classification) and plant type (broadleaved tree, conifer tree, grass, or shrub).We added the factor (climate orplant type) as a moderator of the random model, with the study as the random factor.In the case of plant type, the random factor was plant nested within the study.For those analyses, we only used data from the earlier decomposition database and the effect of decomposition rate, C, N, P, and K contents and ratios.
We also conducted a publication bias analysis using Egger's regression test (Egger et al 1997), adapted to the rma.mv function.We adjusted the metaregression models by including the standard error of the effect sizes as a moderator.If the intercept significantly differed from zero (P < 0.1, following Egger et al. 1997), we inferred that the relationship between precision and the size of the studies is asymmetrical, indicating bias.

General effects
We explored the general effects of the treatments over the experiments.The detailed data of the significant effects found in the general effect analysis are presented in Table S4.We display the data in C, N, and P ratios and contents in the litter remaining after decomposition.More analysis on C, N, and P concentration are available on Supplementary Material (Figure S2 and Table S5).
We identified publication bias in 26.1% of the analyses, which warrants some caution on the inferences presented here.'.The analysis affected by the bias are depicted in Figs. 2, 3 and 4 (see figure captions) and in Table S4.

Litter decomposition
The amount of mass remaining was lower in N + P enriched plots than in controls in both the earlier (6 months) and later stages (12 months).However, an increasing effect in litter decomposition rates was only found for later stages, while earlier stages did not show a significant effect.Plots under the irrigation treatment presented significantly higher decomposition rates than control plots in the later stages of  2).When analysing the sensitivity of the effects to the treatment intensity, we found that in the case of N + P enrichment experiments, the more N is added the higher the effect, but this was only found for later stages.Warming, drought, and N or P enrichment (alone) did not show any effects on litter mass decomposition.

Effects on the nutrient contents remaining
Compared to control plots, C content in litter decreased in the warmed treatments during earlier and later stages of decomposition.C content in litter also decreased in N + P enriched plots in the later stages of decomposition.N content in litter was incremented in later stages of decomposition under N + P enrichment conditions.No effect was found for N and P remaining content or release rate (Fig. 3).There were not enough (< 5) P enrichment experiments valid for assessing its effects on element contents remaining since the beggining of the experiments.

Effect on the nutrient ratios
Compared to control plots, litter C:N ratios decreased in P enriched plots in both earlier and later stages of  5).The sensitivity analysis revealed that the more P was accounted for in the P enrichment experiments, the higher the effect is over the C:N ratio.The same happened with N enrichment intensity in the N enrichment experiments, the more N was accounted for, the higher the effect was over the C:N ratio.Litter C:P ratio decreased, compared to control plots, under warming treatment in both earlier and later stages of decomposition.Warming also increased the C:N ratio change rate in earlier and later stages.N + P enrichment treatment decreased C:P ratios in the later stage of decomposition while increasing C:P change rate in the later stage of decomposition (Fig. 5).The sensitivity analysis revealed an indirect relationship between the intensity of the warming and the decreasing effect on the C:P ratio.The litter N:P ratio decreased, compared to control plots, with the irrigation treatment for the earlier stage of decomposition, however no other effects were found for N:P ratios of N:P change rates (Fig. 4).

Effects considering the climate of the sampling area
The warming experiments from cold climates presented a decreasing effect on C contents.Regarding N content, irrigation treatment had a decreasing effect on cold climates while drought had an increasing effect on arid climate.Although general changes in P content were not significant, different climate conditions led to significant changes in the P content significantly changed.In this regard, under warming experiments, we found an increasing effect on P content in cold climates, and decreasing in polar climates.In drought experiments, we found an increasing effect on P content in arid climates.In irrigation experiments, we found an increasing effect on P content in temperate climates.Last, in N + P enrichment experiments, we found an increasing effect on P content in temperate climates and a decreasing effect on P content in arid climates.
C:N ratio was also affected by climate in warming and N + P enrichment experiments despite not having any general significant effect.We found a decreasing effect of warming experiments on C:N ratio in cold climates.Also, in N + P enrichment experiments, we found decreasing effect on C:N for polar and temperate climates.In N enrichment experiments, which had a C:N ratio general decreasing effect, Climate-related effects were found for N enrichment experiments:C:N ratio, which showed a general decrease effect, decreasedin cold and tropical climates; C:P ratio was affected by climate in N + P and N enrichment experiments although no general effects were found.We found, for N + P enrichment experiments, a C:P ratio decreasedt in cold and tropical climates.We also found a N enrichment decreasing effect on C:P ratios in cold climates.For warming experiments, which showed a general decreasing of C:P ratios, we found that cold climates were driving the effect.All this information is summarised in Table 1.

Effects considering the plant type
Although there was no general significant effect on litter decomposition rates under warming treatments, we found a increasing effect when looking at grass species Also, in N + P enrichment experiments under conifer trees, an increasing effect on decomposition rate.
We found a general decreasing effect of C contents in warming experiments (Table 1).Additionally, during warming experiments under grass species, a decreasing effect of the C content in the remaining litter was found.Although there was not a general effect for C content on drought experiments, we found a decreasing effect for C content when lookint at shrub species.Despite not finding any general effect for N content in irrigation experiments, we found that in grass species there is a decreasing effect on the N content.P content did not have any general effect for N + P enrichment Fig. 5 Graphical summary of the general effects of the metaanalysis on decomposition rate and C, N, and P contents and their ratios.General effects, while specific effects only from earlier stages of decomposition.Significant effects from the climate and plant type-specific analysis that were only found from one experiment are excluded from the figure and the Discussion section.Consult Table S4, Table 2, and Table 3 for p values and number of studies for each analysis experiments, but we found an increasing effect on P content driven by litter from broadleaved tree leaves.C:N ratio did not show any general effect for the warming, drought, and N + P enrichment treatments.However, C:N ratios decreased in grass species on warming experiments; in conifer species on N + P enrichment experiments and in grass species in drought experiments; while increased in shrub species in drought experiments.For warming experiments, which showed a general decreasing effect on C:P ratios, we found a decreasing effect for grass species and conifer trees.Irrigation experiments showed a decreasing effect of N:P ratio.This effect was also found for broadleaved trees.All this information is summarised in Table 2.

K concentration
With the recorded data of K concentration in litter (Table S5; Figure S2).We detected a general decreasing effect in K concentration for warmed plots in later stages of decomposition (Table S5).However, while there is not a general effect during early stages of decomposition, there are climate specific effects in the early stage, increasing its concentration in cold climates (p = 0.007, from 5 of the 6 experiments addressing K concentrations in warming experiments), whereas we found an opposite effect for temperate climates, where K concentrations increase (p = 0.044 from 1 of the 6 experiments addressing K concentrations in warming experiments).We also find K concentrations decreasing in irrigated plots from temperate climates (p < 0.001, from 1 of the 6 experiments addressing K concentrations in warming experiments).

Summary of the general effects
To facilitate the comprehension of the results of the meta-analysis, a summary figure was constructed with results groups by the distinct scenarios that the global change experiments portrayed (Fig. 5).

Discussion
Our meta-analysis provides an overview of the dynamics and stoichiometry of litter decomposition under drivers of global change such as warming, changes in water availability, and nutrient enrichment and limitation.The results portray a guide of the C:N:P stoichiometry under the fate of global change.

Effect of warming
Our initial predictions (i) posited that increases in temperature would lead to a decrease in C:N and C:P ratios.Examining the general effects, elevated temperatures resulted in C release and reduced C:P ratios, consistent with the expectations.However, in warming experiments there were no general effects on litter decomposition rate or C:N ratios.Water limitations may be hindering the effect or warming on litter decomposition.Similar results were obtained by Petraglia et al.(2019) in a national decomposition study reflecting a broad range of vegetation types and microclimatic conditions.The effects of temperature increase on litter decomposition were present in moist conditions.Alternatively, they pointed to vegetation type as a possible explanation since most of the moist soils were classified as grasslands.We found similar results: although a general effect of warming on decomposition rate was absent, there was an increasing decomposition rate effect in litter from grass species.Litter dynamics in cold climates are likely to be disturbed with increments in temperature and precipitation, potentially releasing the old and stabilised pools of organic matter (García-Palacios et al. 2013).When accounting for climate, we detected that C release under warming treatment was notably prominent in cold climates.The warming effect of litter decompostion in cold climates is enhanced since decomposition in these environments is limited by low temperatures which are raised in the warming treatments.Additionally, litter decomposition is promoted since cold climates are charecterised by precipitation exceeding evaporation (Beck et al. 2018), being non water-limited by definition.Consequently, our results confirmprevious predictions (Sardans et al. 2012a) that warming effects on C:N and C:P ratios are only Table 2 Statistical details of the meta-analysis of the significant effects considering the plant type of the sampling area as a moderator.The number of experiments using each plant type is expressed relative to the total number of experiments from which the general effect was calculated (i.e. from the 38 exper-iments addressing decomposition rates on warming treatments, 12 were done in grass litter).These analyses were only carried out in the earlier stage of decomposition for litter decomposition rate; C, N, and P contents; and C:N, C:P, and N:P ratios pronounced under water-abundant scenarios.In cold climates, P is immobilised in litter, explaining why the warming-induced decrease was only found for C:P.Concerning plant type, we found that warming in grasses significantly increased decomposition rates and reduced C contents, C:N, and C:P ratios, aligning with previous studies (Song et al. 2014).In the experiments with grass litter, which is generally characterised by lower N concentration than tree litter (Yuan and Chen 2009;Vallicrosa et al. 2022), decomposition was hindered.This could elucidate why warming effects on decomposition rate and C:N were exclusive to grasses.Warming in cold grasslands and conifer forests enhances C release and P immobilisation in litter.Released C may contribute to soil and atmospheric carbon, exacerbating greenhouse effects and global warming (Cornelissen et al. 2007;Hobbie et al 2002).P-limitation arises from P immobilisation in litter, rendering it unavailable as a resource for plants.Consequently, C:N and C:P ratios are reduced in these scenarios, potentially impacting soil trophic webs.Another concern, that was nos addressed in the meta-analysis, is the community shift towards a more recalcitrant leaf litter, shrub dominance, in response to increasing temperatures in cold climates (Cornelissen et al. 2007;Jackson et al. 2002).
Lastly, K release from the litter increased during the warming treatments.The effect was relevant in cold climates, that is non-water limited.However, we detected an opposite (increasing litter K concentration) effect in a case in a Mediterranean shrubland (Saura-Mas et al. 2012).Since K is an important element for plant water use efficiency (Xu et al 2021), it might be immobilised in litter as a result of biological strategies to cope with water limitation.

Effect of water availability
Our predictions (ii) postulated that water scarcity would elevate C:N and C:P ratios, hindering decomposition, with the opposite expected under increased water availability.C:N, C:P, and the contents of C, N, or P did not exhibit consistent increases or decreases in both drought and irrigation scenarios.General effects were observed for irrigation, increasing litter decomposition rate and decreasing the N:P ratio.Irrigation-induced reduction in the N:P ratio may be attributed to the low solubility and limited soil mobility of P (Manzoni et al. 2010), favoured by rising soil water content.
We also predicted (v) that increases in temperature and changes in precipitation will have greater effects in cold and dry climates (cold, arid and polar) compared to warm and wetter ones (temperate and tropical).In the case of warming, addressed in the section above, we did not find a greater effect of decomposition in cold climates, but we did for C release, P immobilisation and C:P ratio decrease.While there were no other specific effects related to climate.Regarding changes in the precipitation patterns, we did not find any effect regarding decomposition rates, but N and P were immobilised during drought treatment; and N released in irrigated plots from cold climates.Temperate and tropical climates were not affected specifically by any climate change driver.In cold climates (non water-limited), the accelerated N release may be caused by a facilitated hydrolysis of litter compounds rich in C and N, such as cellulose, lignin, and macroproteins.This could result in a higher mineralization rate for C and N compared to P. Conversely, under drought conditions in arid climates, N and P immobilisation suggesting an increased nutrient use efficiency under water-limiting conditions (Sardans et al. 2012a).
Examining the specific effects of plant type, the general decreasing effect on the N:P ratio was primarily driven by experiments involving litter from broadleaved trees.Additionally, a significant effect on N release was noted in grasses, potentially explained by the lower N and higher P concentrations in the senescent leaves of grasses compared to broadleaved trees (Yuan and Chen 2009;Vallicrosa et al. 2022).Broadleaved tree litter, with higher N:P ratios, may be more sensitive to changes in water availability than grass litter..

Effect of nutrient enrichment
Our prediction (iii) suggested that nutrient enrichments would increase decomposition activity, leading to reductions in C:N and C:P ratios.However, the N:P ratio would differ depending on the enrichment, increasing in N-enriched treatments and decreasing in P-enriched treatments.For C:N ratios, the prediction was confirmed, generally decreasing under N, P, and N + P enrichment treatments.Under N + P enrichment, the decomposition rate increased, C:P ratio decreased (in contrast to N or P enrichment scenarios), and C was released while N was immobilised.A meta-analysis by Ferreira et al. (2015) reported increased litter decomposition only when both N and P were applied, suggesting a co-limitation of these elements for faster decomposition.Imbalanced nutrient enrichment, favouring N over P, could lead to N-limitation, explaining the lack of general effects of N enrichment on decomposition (Colin and Waring 2017).The observed release of C and immobilisation of N in N + P enrichment scenarios suggest a preferential metabolism of C-rich compounds by decomposing microorganisms (Manzoni et al. 2021;Hobbie et al. 2012)).
In climate-specific analysis, cold climates exhibited a prominent decrease in C:N and C:P ratios under N enrichment and C:P decrease under N + P enrichment.The combination of nutrient enrichment and water-abundant environments intensified the activity of fast-decomposing communities, leading to a preferential mobilisation of C. A 17 year old experiment on N addition in cold climates (Maaroufi et al. 2017) suggests that C mobilisation from litter in response to N enrichment was driven by changes in soil microbiota.C:P ratio also prominently decreased in tropical climates, where water is not scarce (Beck et al. 2018).Under N + P enrichment, P immobilisation occurred in temperate climates but was released in arid climates.Further analysis addressing the interaction of nutrient enrichment and water availability is necessary to fully understand this effect.
We predicted (vi) that conifer, high-lignin, litter would be less affected by nutrient enrichments compared to broadleaved trees or grass litter.However, we detected the opposite: conifer tree litter had a prominent decomposition rate increment and C:N decrease.Broadleaved tree litter exhibited a significant P immobilisation effect in N + P treatment.The differing N:P ratios in conifer and broadleaved tree litter, with lower N:P ratios and higher lignin content in conifer leaves (higher P contents, Rahman et al. 2013), may explain why conifer tree litter is less P-limited, fostering higher decomposition rates.The correlation between climate and plant type as moderators could contribute to these effects, especially in temperate climates where broadleaved tree experiments predominated.Opposite effects between grass and tree species were not found, as could be expected by previous meta-analysis comparing the effects of N addition in forest and grasslands (Su et al. 2022).
Earlier vs later stages Contrary to our predictions (iv) expecting a decrease in C:N and C:P ratios only in earlier stages of decomposition, we did not find a consistent pattern.Instead, distinct effects between earlier and later stages were observed, particularly under N + P enrichment.In cases where there was a difference between stages, the later stage exhibited significant effects.Earlier stages of decomposition are less specific, as they consist on the mineralisation of labile C.However, later stages of decomposition require complex metabolites to degrade the recalcitrant tissues (Djukic et al. 2018).Impacts on the microbiota due to nutrient enrichments (Maaroufi and De Long 2020) may be more critical during a less generalistc later stage of decompositon.Moreover, the systemic impact of nutrient enrichment on litter decomposition stoichiometry appears relevant in the long term, supporting the notion that N + P enrichment effects are cumulative (Hanley et al. 2021;Seabloom et al. 2021).This aligns with explanations for N immobilisation under N + P enrichment, where N-limitation results from a prolonged metabolic process.

Limitations of the study
Our initial question was about how drivers of global change interacted with the elemental stoichiometry in leaf litter, so we only selected articles that sampled at least two elements simultaneously in litter samples for control and treated plots.The number of studies/ experiments to test our hypothesis was thus reduced.We, therefore, used data from experiments of decomposition in litterbags but also from litter sampling without decomposition experiments, from which we could only extract data for the differences in elemental concentrations or ratios between control and treated plots but not for the rates of mineralisation or decomposition.Caution is needed when interpreting these results, some aspects of the meta-analysis are low in number of experiments.Regarding the sampled sites (Fig. 1), there is a clear example of the geographic bias, frequently found in meta-analyses, where most of the field studies are within the latitudinal ranges of China, Europe, and North America.Also, all arid places are found in Asia.More experiments in other areas are needed in order to infer more precise global patterns.
Another issue to consider is that the processes of decomposition vary across ecosystems, chemical composition of litter, and stages of decomposition.The days of experiments were not consistent across studies; and even if they were, the experiments from each study began in different seasons, which made the appropriate comparison of experiments difficult.Thus, our results can serve as general information for determining how nutrient dynamics may behave under current global change and may differ across ecosystems and decomposing time.The methods used for analysing elemental compositions were not homogeneous and could introduce bias to our results.We provide in Table S2 information about the measurement methodologies applied on each study, and the experiment days that were considered for the earlier and later stages.For further information about specific climate and litter species, we detail the characteristics of each study used in Table S1. is important to highlight the presence of publication bias in 26.1% of the analysis.However, the models took into account the possible non-independence of experiments from the same study variability between studies.The use of plant type and climate, which are correlated variables, should also be carefully interpreted for specific analysis independently, avoiding crossreferencing.The use of the rates of elemental release and change of elemental ratios was useful for identifying effects that were not apparent when focusing only on the elemental contents in the litter.
Another aspect that was not contemplated in this study is the soil biota, crucial for the litter decomposition process (Njoroge et al. 2021), which is affected by nutrient enrichment (Nessel et al 2021) and climate change (García-Palacios et al. 2013;Maaroufi and De Long 2020).Finally, we focused on experiments and observations of one driving factor, but all drivers of global change interact in nature.Therefore, multifactorial experiments and observations are strongly recommended.

Conclusions
Our meta-analysis explored the effects of warming, changes in water availability, and nutrient enrichment, which are known drivers of global change, on the processes of elemental mineralisation/immobilisation and on the stoichiometric relationships during the decomposition of litter.The effects can vary across scales, but they can help us determine the generalised effects of global change on the nutrient dynamics of terrestrial ecosystems.
As the main general effects, we found that the abundance of resources (water and N + P) fosters the decomposition of litter.Warming and nutrient enrichment led to losses in litter C content and reduction of the C:N and C:P ratios.C can move from the litter to the soil but also to the atmosphere as CO 2 , contributing to the greenhouse effect and generating a feedback loop, because more warming derived from the greenhouse effect may intensify the release of C from litter.Subsequently, the C: N and C: P ratios were reduced in these scenarios, with potential impacts on the soil trophic webs.This effect is prominent in grassland and coniferous forests in cold climates.
The N:P ratios were only affected by the increase in the availability of water in broadleaved tree leaf litter, which may be caused by the higher solubility of N with respect to P. This effect was only apparent at earlier stages of decomposition.Our analysis found that N:P ratios, in contrast with C:N and C:P ratios, were more resilient to changes in temperature, water scarcity, and nutrient enrichment.
We identified immobilisation of P in the litter, which may hinder P-limitation scenarios, in warming and cold climates; arid environments experiencing droughts; temperate environments with increasing water availability; and in temperate or on broadleaved tree leaf litters experiencing N and P enrichment.N immobilisation was detected under the combined enrichment of N and P at later stages of decomposition.Arid environments experiencing droughts are subject to both N and P immobilisation.This exploratory analysis is useful for distinguishing between various impacts of climate change and nutrient enrichment on ecosystem nutrient dynamics.We aim to inspire new and warranted research on the unknown mechanisms responsible for the results and interactions amongst multiple drivers.Future research with experimental designs is encouraged research to simulate scenarios of global change for exploring the stoichiometric implications of the simulated scenarios, with the maximum possible number of drivers interacting and in different stages of decomposition.

Fig. 1
Fig. 1 Global map of the sampling locations (coloured dots) used in the meta-analysis

Fig. 2
Fig. 2 Effect sizes (± 95% CI) of the treatments on the amount of the initial litter mass remaining after approximately half a year (earlier stage, green bars, above) or approximately a year (later stage, pink bars, below) of decomposition (a) and the decomposition rates (b).Coloured bars indicate significant effect sizes.The boxes depict the number of experiments (simi-

Fig. 3
Fig. 3 Effect sizes (± 95% CI) of the treatments on the content and k for the contents and release rates of C (a, b), N (c, d), and P (e, f) approximately half a year (earlier stage, green bars, above) or approximately a year (later stage, pink bars, below) of decomposition (a) and the decomposition rates (b).Coloured bars indicate significant effect sizes.The boxes depict the number of experiments (similar to earlier and later stages).Treatments are only represented in the plots if the number of experiments is ≥ 5. Publication bias is highlighted with warning symbols (triangle with exclamation mark inside) in the title of the treatment.Green warning symbols depict publication analysis for the earlier stage; while pink warning symbols are for the later stages

Fig. 4
Fig. 4 Effect sizes (± 95% CI) of the treatments on the stoichiometric ratios and associated change rates for C:N (a, b), C:P (c, d), and N:P (e, f) approximately half a year (earlier stage, green bars, above) or approximately a year (later stage, pink bars, below) of decomposition (a) and the decomposition rates (b).Coloured bars indicate significant effect sizes.The boxes depict the number of experiments (similar to earlier and later stages).Treatments are only represented in the plots if the number of experiments is ≥ 5. Publication bias is highlighted with warning symbols (triangle with exclamation mark inside) in the title of the treatment.Green warning symbols depict publication analysis for the earlier stage; while pink warning symbols are for the later stages , we found that 7 of the experiments were carried out in tropical climates, 48 in arid climates, 55 in temperate climates, 33 in cold climates, and 32 in polar climates.The starting months for the litterbag experiments differed across studies: 10 experiments started from January to March, 63 from April to June, 12 from July to September, and 21 from October to December.Plant species and climate were correlated variables.

Table 1
Statistical details of the meta-analysis of the significant effects considering the climate of the sampling area (according to Köppen-Geiger classification; in the brackets) as a moderator.The number of experiments under each climate is expressed relative to the total of experiments from which the general effect was calculated (i.e. from the nine experiments that addressed C content in warming experiments, eight were from cold climates).These analyses were only carried out in the earlier stage of decomposition for litter decomposition rate; C, N, and P contents; and C:N, C:P, and N:P ratios Vol.: (0123456789)