One plus one is greater than two: mixing litter types accelerates decomposition of low-quality alpine dwarf shrub litter

Litter decomposition patterns, non-additive effects, and spectral data of abundant alpine leaf litters were assessed in litterbag experiments containing single species and mixtures. We tested if low-quality shrub litter decomposes faster in mixtures with high-quality litter and if predictions on decomposed litter using spectral data are feasible. Chemical and physical traits and near-infrared reflectance (NIR) spectra of six alpine freshly fallen litter types were measured. A litterbag experiment (0.1 mm mesh size) with single and 2- and 3-species mixtures was conducted with three species from three functional groups (shrub, grass, forb). Decomposition rates, litter mass loss, non-additive effects, and NIR spectra were recorded after 6, 12, and 24 months. The six freshly fallen litter types differed significantly in leaf litter traits and NIR spectra. Decomposition rates steadily slowed during the 24 months, with shrub litter having the lowest on all sampling dates. In litter mixtures, shrub and grass litter showed higher decomposition rates after 12 and 24 months compared with the single-species treatments. The accelerated decomposition in litter mixtures indicates stimulating interactions between low- and high-quality litter types, most probably mediated by their associated microbiomes. By using NIR spectra, we successfully predicted all initial litter traits, but only total carbon content along the entire decomposition period.


Introduction
The decomposition of plant litter material on the soil surface is essential for terrestrial ecosystems by means of influencing vital ecosystem services such as nutrient cycling, carbon sequestration, and above-and below-ground biodiversity (Berg and McClaugherty 2014). The litter decomposition rates strongly depend on (i) temperature and (soil) moisture, (ii) the chemical and physical traits as well as the nutrient availability of the plant material (e.g. litter quality), (iii) site-specific local factors (e.g. microclimate), and (iv) the composition of the soil decomposer community (Coûteaux et al. 1995;Berg and McClaugherty 2014;Bradford et al. 2016).
In most natural habitats, leaves and stems from an array of plant species belonging to different functional groups (e.g. shrubs, grasses, and forbs) contribute to building up the litter layer. During decomposition, these litter typesand their associated microbiomesinteract, and it has increasingly been recognized that patterns are often not predictable from single-species dynamics (Gartner and Cardon 2004;Hättenschwiler et al. 2005). Litter types influence each other's decomposition rates by changing the chemical environment and by physically altering the litter surface, resulting in non-additive patterns (Gartner and Cardon 2004). Cuchietti et al. (2014), for example, found that the decomposition of slow-decomposing low-quality leaf litter can be accelerated by mixing with fastdecomposing high-quality litter.
In alpine and arctic ecosystems, slower decomposition rates compared with lower and warmer habitats were observed (Cornelissen et al. 2007). The lower temperatures along with a more shrubby and recalcitrant plant material, which contains high amounts of secondary plant products such as phenols and lignin, were found to be the main factors leading to slower decomposition rates (Aerts and Chapin 1999;Xu et al. 2010;Duan et al. 2013;Bradford et al. 2016). Possible consequences are litter accumulation and changes in soil properties (Seeber and Seeber 2005), affecting the above-and below-ground plant (Pruchniewicz 2017) and invertebrate species composition (Steinwandter et al. 2017). Increased temperatures due to global warming, which are predicted to be more severe in high-latitude and high-elevation regions (IPCC 2013), could on the one hand enhance decomposition in alpine ecosystems, but on the other hand limit it through reduced (soil) moisture (Aerts 2006) and through lowquality shrub litter replacing high-quality litter species (Cornelissen et al. 2007).
Some studies on litter decomposition in alpine grassland exist (e.g. Wang et al. 2010;Saccone et al. 2013;Liu et al. 2016), but little is known on interactions of alpine leaf litter mixtures (e.g. Xu et al. 2010;Duan et al. 2013;Chen et al. 2015), especially from the Alps. Beneficial interactions (i.e. positive non-additive effects) that accelerate litter decomposition and counteract litter accumulation might also apply to high-elevation ecosystems such as semi-natural grasslands in the Alps, but have not been investigated yet. Therefore, we conducted a litterbag experiment over the course of 24 months with alpine leaf litter from three functional groups (i.e. grass, forb, and dwarf shrub as singlespecies treatments and mixtures). Contrary to most alpine litterbag studies, we used fine-meshed litterbags that excluded soil meso-and macro-decomposers that might feed on the litter material, thus investigating the pure decomposition induced by microbial (Knapp et al. 2011) and physical processes. We assessed leaf litter traits and additionally measured the spectral composition by using near-infrared reflectance spectroscopy (NIRS). To date, only few litter decomposition studies included NIRS (e.g. Joffre et al. 1992;Parsons et al. 2011), although it has been proven an effective method to measure plant and litter chemical traits as well as soil material (Chodak 2008;Agelet and Hurburgh 2010).
In this study, we tested the hypothesis that along 24 months of decomposition, low-quality dwarf shrub litter degrades faster if mixed with high-quality litter. Further, we tested if NIRS can be a suitable tool for predicting litter traits of decomposed leaf litter material. Knowing the litter traits and spectral composition of key litter types from Alpine grassland will considerably enhance our understanding of alpine decomposition patterns. Further, it will help us to better anticipate the effects of ongoing Alpine land-use change (i.e. abandonment) and global warming (Brigham et al. 2018) on alpine litter decomposition as they might negatively affect soil processes, plant species compositions (i.e. due to shrub encroachment), the decomposer community (i.e. change in litter quality), and subsequently a variety of ecosystem services.

Study site
Figure 1 summarizes design and working steps of this study. About 500 g freshly fallen leaf litter of six common and abundant alpine plants as representatives of the functional groups grass, forb, and dwarf shrub were collected at the Alpine Long-Term Ecological Research (LTER) area BKaserstattalm^(Central Alps, Tyrol, Austria; site code LTER_EU_AT_015; 47. 1250°N , 11.2900°E) in autumn 2012. The sampling area, managed and abandoned pastureland, ranged from 1950 to 2050 m above sea level (a.s.l.) in elevation (Steinwandter et al. 2017). The mean annual temperature (MAT) and precipitation (MAP) at 1900 m a.s.l. are 3.0°C and 1097 mm, respectively (Ingrisch et al. 2017). During the experiment, climatic conditions were recorded (Fig. A1). The first winter (i.e. 2012/2013) was colder than the second winter (monthly averages of surface temperatures in December with −10.2°C and − 8.7°C, respectively). In both winters, relatively low snow covers (i.e. about 30 cm) were observed, with more snow occurring in late 2014 (ca. 60 cm, pers. observation). The summer 2013 was relatively dry, while the summer 2014 had more and evenly distributed precipitation (Fig. A1).
Ten chemical and physical litter traits of initial litter with four replicates each were analysed following the protocols of VDLUF method book volume III (2018). These were (i) crude ash by drying material at 105°C and ashing it at 550°C (see chapter 8.1 in the VDLUF method book), (ii) crude protein via the Dumas combustion method (ch. 4.1.2), (iii) crude fibre by processing the material with boiling hydrochloric acid and potash lye, washing, drying, and ashing it at 500°C (ch. 6.1.1), (iv) acid detergent fibre (ADF) by boiling the material for 1 h in acid detergents solution, filtering and washing it (ch. 6.5.2), (v) acid detergent lignin (ADL) by processing the filtered ADF material with 72% sulphuric acid for 3 h, washing and drying it (ch. 6.5.3), (vi) calcium (Ca), phosphorus (P), magnesium (Mg), potassium (K), and sodium (Na) by solving ashed material in hydrochloric acid, washing and drying it before determining it via inductively coupled plasma optical emission spectrometry (ICP-OES, ch. 10). Additionally, total carbon (C) and total nitrogen (N) content were analysed for all replicates and sampling dates using an organic elemental analyser (FlashEA Fig. 1 Scheme of the study design and working steps. Litter traits [1] and NIR spectra [2] data of freshly fallen alpine leaf litter types were combined in PLS models [3] to develop prediction models for the decomposed litter material. In the leaf litterbag experiment, LML, decomposition rates and non-additive effects [5] as well as NIR spectra [6] were measured. The NIR spectra of decomposed litter were used to test the prediction model based on the freshly fallen litter [7] 1112, Thermo Fisher S cientific, Waltham, Massachusetts). Based on these results (Fig. 2, Table 1) and due to their high abundance at the study site, the following three litter types were chosen for the litterbag experiment: V. vitis-idaea (further on referred to as SHRUB_Vv), D. glomerata (GRASS_Dg), and G. sylvaticum (FORB_Gs). Abbreviations are used to facilitate identification of this three litter types but do not represent the entire functional groups.

Leaf litterbags
Each of the following treatments was replicated 15 times: (i) single-species treatments of SHRUB_Vv, GRASS_Dg, and FORB_Gs; (ii) 2-species mixtures of SHRUB_Vv + GRASS_Dg, SHRUB_Vv + FORB_Gs, and GRASS_Dg + FORB_Gs; and (iii) a 3-species mixture SHRUB_Vv + GRASS_Dg + FORB_Gs. Each polyamide litterbag (size of 15 × 15 cm, mesh size of 0.1 mm) was filled with 10.00 g (±0.01) air-dried mass of the corresponding plant material; 10.00 g, 5.00 g and 3.33 g in the single-species treatments, 2-species mixtures, and 3-species mixtures, respectively. In the mixtures, a 1 mm mesh size fly net separated the litter types while at the same time ensuring direct contact. The fine mesh size excluded soil meso-and macroinvertebrates and hindered roots growing into the litter material.
At the incubation site (area of approx. 3 × 10 m of an abandoned pasture at 2000 m), the vegetation was cut to 1 cm height to ensure maximum contact of the litterbags with the surface. The pasture was abandoned since 1983, the vegetation was described as Seslerio-Caricetum with scattered dwarf shrubs (Ingrisch et al. 2017;Steinwandter et al. 2017). The litterbags were placed randomly on the ground and secured with a coarse metal grid. The cut vegetation was scattered onto the litterbags to simulate vegetation presence and to prevent desiccation. In total, 105 litterbags were incubated in October 2012.
After 6, 12, and 24 months, five replicates of each treatment were sampled, cleaned, sorted, and separated into the single-litter types, air-dried, and finally weighed. The litter mass loss (LML) was depicted as the percentage of loss from the initial air-dried mass. The litter decomposition rate k was calculated using the The shaded areas represent 2D-normal-based ellipses covering 95% of data distribution estimated from the variance-covariance matrix of their X and Y coordinates. Legend: triangles represent grass, circles forb, and rhombi the dwarf shrub leaf litter. The arrows represent the litter traits. Eigenvalues: axis 1 = 0.4825, axis 2 = 0.2309 Included are results of a one-way ANOVA and Tukey's HSD post-hoc test of leaf litter traits. Superscript letters indicate significant differences at p < .05. The PLS model explained 97.75% of the variance: Factor-1 (77%, 29%), Factor-2 (20%, 65%). Root-mean-square error and regression coefficient of calibration (RMSEC and R 2 Cal.), root-mean-square error of crossvalidation (RMSECV), and root-mean-square error of prediction (RMSEP) are given SHRUB_Vv Vaccinium vitis-idaea; GRASS_Dg Dactylis glomerata; FORB_Gs Geranium sylvaticum # For carbon (C), nitrogen (N), and the C:N ratio a separate PLS model was used with 4 factors explaining 98.19% of the variance: Factor-1 (77%, 50%), Factor-2 (21%, 40%) § Potassium (K) and sodium (Na) were excluded from the final PLS model due to low regression coefficients of determination (R 2 ) equation ln(M t /M 0 ) = −kt, in which M 0 is the initial litter biomass and M t the biomass at a certain time point t (Olson 1963). In calculating decomposition rates between sampling dates, M 0 was the litter mass at the earlier date, t the time between the dates (i.e. 0.5 and 1.0 years; for this calculation, months were converted to years), and M t the litter mass remaining after t.

Near-infrared reflectance spectra
Near-infrared reflectance spectroscopy (NIRS) represents a fast, low-cost, and non-destructive technology of increasing popularity in ecology (Foley et al. 1998), agriculture (Martinez et al. 2010), biodiversity studies (Kinzner et al. 2015), pharmacology (Reich 2005), and more recently also in soil studies (Chodak 2008;Zangerlé et al. 2014). It analyses elements and compounds with certain chemical bonds that reflect the radiation in different ways. Litter material was stored in an oven for 48 h at 60°C before being milled to powder in a ball mill (Retsch, Haan, Germany) and sieved through a 1 mm screen. Each sample consisted of one to two grams of ovendried litter material which was transferred to a small round plastic container and flattened gently with a spatula to obtain a smooth measuring surface (approx. 1 cm 2 ). Constant room conditions (i.e. temperature, humidity, and light) and the distance from samples to the NIRS light emitting probe (diameter of 3 mm) were identical across all NIRS measurements. Each sample was scanned five times according to a crosswise scheme (i.e. centre and cardinal points) with a portable VIS/NIR spectrophotometer (LabSpec 5000, ASD, Boulder, Colorado); each spectrum was achieved by averaging 50 measurements. Each replicate was measured five times, resulting in 25 spectra per replicate, 100 spectra per initial, and 125 spectra per decomposed litter type. We measured the reflectance ranging from 300 to 2500 nm and 1 nm resolution; regular calibrations were made using a white Spectralon® plate with 2.2 mm distance to the probe (Kinzner et al. 2015). The reflectance spectra (R) were converted to absorbance spectra (A) using the eq. A = log (1/R).

Statistical analyses
Differences in plant components and LML were analysed using Analyses of Variance (ANOVA, at significance level p < .05 and Tukey's HSD post-hoc tests) with the factors sampling date * litter type. Linear mixed-effects models with random factor time were calculated for LML across the entire decomposition period using the package LME4 (version 1.1-19, Bates et al. 2014) in the open-source statistical programming language R (version 3.5.2, R Core Team 2018) in RStudio (version 1.1.463, RStudio Team 2016). Principal Component Analyses (PCA) of plant components were calculated using Canoco 5.11 (Ter Braak and Šmilauer 2018).
Non-additive effects on LML are defined as the deviance from the sum of the purely LML observed for each individual plant species in litter mixtures (Hättenschwiler et al. 2005). They were calculated as the difference of the individual observed LML (of the single litter types in the mixtures) minus the expected mean LML (of the corresponding single-species treatments). If the difference is a positive value, the effects are called synergistic; when negative, antagonistic.
All NIR spectra data transformations and calculations were performed with the multivariate data analysis software The Unscrambler X 10.5 (CAMO Software, Oslo, Norway). All spectra were pre-processed: First, the wavelength range was reduced from 500 to 2300 nm (i.e. 1800 data points) due to high distortion at the endpoints. Then, the spectra were smoothed using Savitzky-Golay derivates (polynomial order 2, 11 smoothing points) and Extended Multiplicative Scatter Correction (EMSC) (Agelet and Hurburgh 2010; CAMO Software 2017). Finally, the spectra were reduced by factor 5 corresponding to the average spectra of the five crosswise measurements per plastic container. Thus, we obtained 20 spectra for each initial (i.e. 0 month) and 25 spectra for each decomposed litter type. The covariance of litter traits and NIR spectra data was maximised in Partial Least Squares (PLS) regression models with the Wide Kernel PLS algorithm; the litter traits data were centred and standardised. One outlier was removed according the Hoteling T 2 outlier detection (Agelet and Hurburgh 2010; CAMO Software 2017). We finally proceeded with the best PLS model using just the three main litter types data and 12 out of 14 plant components (i.e. excluding K and Na), resulting in 97.75% average explanation of variance using four factors (for details see Table 1). Additionally, we validated the PLS model with blind samples gaining satisfying results (i.e. R 2 greater than 0.90).

Plant components and near-infrared reflectance spectra
The six freshly fallen alpine litter types significantly differed in their litter traits ( Fig. 2; Table 1 and A1) and NIR spectra ( Fig. A3 and A4). In FORB_Gs Ca concentration was very high with almost 45 mg g −1 , while all grasses and dwarf shrubs remained below 14 mg g −1 ; for Mg we found a similar but less pronounced result (9.48 mg g −1 vs. 2.08 mg g −1 in FORB_Gs and SHRUB_Vv, respectively). However, FORB_Gs had the lowest P, total C and N content, the latter resulting in the widest and most variable mean C:N ratio of 47.67 ± 6.81. The two dwarf shrub species showed the highest total C content exceeding 500 mg g −1 . The three litter types used in the litterbag experiment differed significantly in all measured plant components except the Na content (Table 1).
Along the decomposition process, total C and N content changed differently for each litter and mixture type (Table A3). The total C content of the SHRUB_Vv single-species treatments increased only slightly, whereas the total N content significantly increased stepwise at each sampling date. Therefore, also the C:N ratio significantly decreased from 39.82 ± 2.09 to 25.21 ± 1.17. The total C and N content from the SHRUB_Vv material within litter mixtures, in most cases, were similar to the single-species treatments. However, the SHRUB_Vv + F (i.e. SHRUB_Vv litter in the SHRUB_Vv + FORB_Gs mixture) and the 3-species treatments (SHRUB_Vv + GF) had the significantly lowest total C content after 6 and 12 months, and therefore also a significantly narrower C:N ratio compared with the single-species treatments. In contrast, the C and N content of GRASS_Dg and FORB_Gs litter increased significantly after six months and levelled off at this concentration for the rest of the experiments. As a consequence, the C:N ratio of GRASS_Dg and FORB_Gs dropped from 30.92 ± 4.12 and 47.67 ± 6.81, respectively, to values around 20.
A clear separation of four out of six alpine litter types can be seen; the grasses N. stricta and L. lutea showed low correspondence to the plant components in the PCA plot (Fig. 2). The intra-specific variation was highest in D. glomerata (i.e. GRASS_Dg) showing the largest grouping area (green ellipse). The result of the PCA plot supports the selection of the abundant species V. vitis-idaea, D. glomerata, and G. sylvaticum (i.e. SHRUB_Vv, GRASS_Dg, and FORB_Gs, respectively) for the litterbag experiment, as they are located the farthest away in the plot (Fig. 2).
In addition, SHRUB_Vv, GRASS_Dg, and FORB_Gs differed significantly regarding their NIR spectra during the 24 months of decomposition (Fig. A5). The four sampling dates followed a similar pattern for all three litter types: considering their NIR spectra, the initial freshly fallen litter was always the most different one located on the opposite side from the decomposed litter in the PCA plots (Fig. A5). The NIR spectra of the latter were found to cluster together but were still clearly separated.

Mean litter mass loss and decomposition rates
The decomposition rates differed significantly between treatments within each sampling date and along the entire decomposition period, and slowed down with time (Tables 2 and 3; Fig. A6). SHRUB_Vv litter decomposed significantly slower compared with the other two litter types with a mean LML of only 25.4% ±1.7 after 24 months. Further, it had the significantly lowest mean decomposition rates at each sampling point; the initial rate after six months (k = 0.25 ± 0.04) was even less than half of that of FORB_Gs after 12 months. FORB_Gs showed the highest (i.e. fastest) mean litter decomposition rate with 64.6% ±2.9 LML after 24 months, followed by GRASS_Dg (LML of 47.0% ±1.6, Table 2 and Fig. A6).
In the litterbags with 2-and 3-species mixtures, we observed a deviation of the collective LML from the average LML of the single-species decomposition rates. Accelerated decomposition rates were found especially for the 12 and 24 months treatments as observed by a positive deviation from the mean calculated from the single-species treatments (Table 3). For example, after 12 months a mean of 19.2% of the SHRUB_Vv and 35.4% of the GRASS_Dg single-species litter were decomposed. In the SHRUB_Vv + GRASS_Dg litter mixture, 33.2% were decomposed after 12 months, i.e. 5.9% more than the average of both single-species treatments (i.e. 27.3%). Again, decomposition rates gradually slowed down with time, except for the mixture GRASS_Dg + FORB_Gs, which had the same mean k-values for the 6 and 12 months treatments (Table 3).

Non-additive effects
The individual decomposition rates of litter in the 2-and 3-species mixtures generally showed highly variable patterns within and among treatments (Table 2). For example, after six months of decomposition no clear pattern could be found in the 2-and 3-species mixtures, with both lower (e.g. significant for SHRUB_Vv + GF) as well as higher decomposition rates compared with the single-species treatments. At this point, the litter types seemed not to profit from one another, especially GRASS_Dg and FORB_Gs even seemed to hamper each other. However, after 12 months of decomposition 12 and 15 out of the 18 mixtures showed accelerated (but not always significantly different) mean decomposition rates and LML, respectively, compared with their corresponding single-species treatments ( Table 2). The SHRUB_Vv litter decomposed slightly faster in the 2and 3-species mixtures, resulting in mean synergistic effects for LML of up to +3.5% and +1.9%, respectively, which was highly significant for SHRUB_Vv + G considering the entire decomposition period (p < .001). Additionally, the mean decomposition rate k indicated slightly higher values for each SHRUB_Vv mixture material after 12 months; the effect was significant in the 3-species mixture. After 24 months, a similarly positive but less pronounced effect of the decomposition rate was found for both the SHRUB_Vv and the GRASS_Dg mixture but not for FORB_Gs ( Table 2).
The non-additive effects of SHRUB_Vv litter in the 3species mixture results in the steepest positive slope of the linear trend line (Fig. 3), implying the highest acceleration of decomposition when more species were present. We found a similar picture for GRASS_Dg with a highly significant positive effects for GRASS_Dg + SF but not for FORB_Gs, in that decomposition was only slightly accelerated when combined with GRASS_Dg (Fig. A7).

Partial least squares predictions
Including all six alpine litter types resulted in lower calibration, cross-validation, and prediction scores compared with including only the three main species which were used for the litterbag experiment (i.e. SHRUB_Vv, FORB_Gs, and GRASS_Dg).
Using the PLS models to predict the values for decomposed litter worked well only for some litter types and plant components (Table A2). For example, Table 2 Estimated deviants from intercepts (i.e. single-species treatments) for mean litter decomposition rates k and mean litter mass loss (LML) of single-species treatments and 2-and 3-species mixtures of three alpine leaf litter types after three time points along the decomposition process Standard deviations are given in parentheses. For litter from mixtures, letters after the plus sign indicate the influencing litter type. The results of a two-way ANOVA (with Tukey's HSD post-hoc test) within sampling dates are given with k * influence and LML * influence, respectively. Superscript letters indicate differences between the groups at p < .05 level. For litter from the mixtures, positive (green) and negative (red) deviations from the expected values (i.e. the mean of the corresponding single-species treatments) are presented. The results of linear mixed-effect models over the entire decomposition period are given (lme). n = 5. SHRUB_Vv Vaccinium vitis-idaea; GRASS_Dg Dactylis glomerata; FORB_Gs Geranium sylvaticum the total C content values of the 12 months old SHRUB_Vv and GRASS_Dg single-species litter were predicted very well (average deviation +0.3% and + 2.7%, respectively), while we obtained high deviations for the FORB_Gs single-species treatment (+12.7%) that exceeded the 95% confidence limit. On the other hand, total N content and thus also C:N ratio predictions were both highly underestimated for the former (up to −25.8% for SHRUB_Vv + G) and highly overestimated for the latter (up to +91.1%, FORB_Gs). We obtained a similar outcome for the six and 24 month old litter predictions, the total C content showing proper values, while those of total N content and C:N ratio were again predicted incorrectly (Table A2).

Discussion
To our knowledge, we are the first to combine both litter traits and spectral data of six common alpine plant species; three of them even over the course of 24 months of decomposition. Generally, more litter decomposition studiesideally multiple-method studies like this oneare needed to better understand decomposition patterns of multi-species litter, not only in alpine ecosystems. The investigated Alpine decomposition patterns were in line with our expectations based on previous studies (e.g. Rief et al. 2012) and the literature as in this study recalcitrant SHRUB_Vv litter showed the slowest and both high-quality GRASS_Dg and FORB_Gs litter a significantly faster decomposition (Montané et al. 2010;Duan et al. 2013;Berg 2014). Further, we found both synergistic and antagonistic nonadditive effects on LML in litter mixtures, similar to the few-species mixtures in the study of Duan et al. (2013). Analyses of the NIR spectra of litter material successfully separated litter types, functional groups as well as the decomposition stages. However, reliable predictions were possible only for the total C content due to methodological constraints.
Near-infrared reflectance spectra and partial least squares predictions The successful prediction of soil properties by using NIRS such as total C and N content, soil moisture, and   Standard deviations are given in parentheses. The results of a two-way ANOVA (with Tukey's HSD post-hoc test) within sampling dates are given with k * influence and LML * influence, respectively. For the mixtures, positive (green) and negative (red) deviations from the expected values (i.e. the mean of the corresponding single-species treatments) are presented at the bottom. Superscript letters indicate differences between the groups at p < .05 level. n = 5. SHRUB_Vv Vaccinium vitis-idaea; GRASS_Dg Dactylis glomerata; FORB_Gs Geranium sylvaticum others, has been shown by Chang et al. (2001), resulting in correlations of measured and predicted values that exceeded R 2 > 0.80. In this study, we successfully validated our PLS model with the initial litter material. Further, we were able to validate PLS predictions for total C and total N content, and thus also for the C:N ratio of the decomposed SHRUB_Vv, GRASS_Dg, and FORB_Gs litters (Table A2). These parameters were measured for all the litter material along the entire decomposition experiment. However, when attempting to predict the other litter traits of the decomposed from the freshly fallen litter material, we encountered difficulties in terms of extremely high deviations of the predicted values. Generally, litter traits considerably change during the decomposition process (Cornelissen and Thompson 1997;Berg 2014). Therefore, the PLS model which was built on the initial litter material (i.e. 0 months) was not able to predict the litter traits of the decomposed litter (6 to 24 months). In other words, PLS predictions work properly only for the same or very similar kind of Bproduct^(CAMO Software 2017). We assume that, as the decomposed litterbag material differed significantly in their spectral composition after 6, 12, and 24 months in the field (Fig. A3 and A4), also their litter traits differed too much to be classified as similar Bproduct^. However, we believe that measuring material at each decomposition stage and use that data for building corresponding PLS models would result in highly accurate predictions, as we could already prove for the total C and N contents.
Thus, NIRS offers a powerful tool in soil ecology and in decomposition studies (Albrecht et al. 2008). According to our experiences with NIRS and the PLS predictions of alpine litter traits along the decomposition process of 24 months, we advise everyone to collect and estimate generously the needed initial freshly fallen litter material and to be aware that both the litterbag approach and the NIRS measurements need a sufficient amount of litter material left at the end of the decomposition experiments (at least 20 g air-dried mass).

Leaf litter traits
For freshly fallen litter material, we found four out of the six alpine plant species to be significantly different in their litter traits and therefore well separable (Fig. 2). Our results are in line with the study of Rief et al. (2012), which was conducted at the same study site and compared the palatability of eight alpine plant species to invertebrate decomposers; five of them were also used in our study. Regarding decomposability, the plant components can be divided into three groups: rapidly degrading carbohydrates and proteins (the latter expressed here as N), slower degrading structural components such as cellulose and hemicellulose (ADF), and recalcitrant lignin and polyphenols (ADL) (Berg and Laskowski 2005).
The three litterbag species, each belonging to a different functional groups, differed significantly in their initial chemical traits (Fig. 2). The evergreen dwarf shrub V. vitis-idaea (i.e. SHRUB_Vv) can be found mainly in abandoned pastures at the study site (Steinwandter et al. 2017) along with other Vaccinium species; it is wide-spread in natural alpine heathland of Eurasia (Meusel et al. 1978) and North America (Flora of North America Editorial Committee 2009). The freshly fallen SHRUB_Vv litter was characterised by the highest lignin content (i.e. ADL), which was also reflected by a high total C content that exceeded 50% of their total litter biomass; on average, lignin is composed of 65% of carbon (Montané et al. 2010). Additionally, freshly fallen litter of Vaccinium spp. and other dwarf shrubs was reported to contain higher amounts of secondary metabolites such as phenols and tannins compared with grass and forb species (Montané et al. 2010;Berg and McClaugherty 2014), thus being responsible for low palatability to large decomposers such as earthworms (Rief et al. 2012). These recalcitrant compounds impair the establishment and growth of microorganisms by forming resistant complexes with proteins and nitrogen (Hättenschwiler and Vitousek 2000), further inhibiting and regulating decomposition processes. However, Rief et al. (2012) reported that after 12 months in the field, the total phenolic concentration in aged dwarf shrub litter decreased considerably (e.g. approximately 90%), most likely making aged SHRUB_Vv litter again more attractive to the decomposer community.
We found a relatively low total N content in freshly fallen litter types, with SHRUB_Vv exceeding that of FORB_Gs but not GRASS_Dg. Generally, a high N concentration was found to boost the decomposition of hemicellulose (i.e. ADF) in the first phase of decomposition, but to hamper that of lignin (i.e. ADL) in the long term, thus leading to lower decomposition rates (Berg 2014). This might be true for the SHRUB_Vv litter material where we found high initial lignin:N, lignin, and N values (Table 1). Further, as the total C content of SHRUB_Vv treatments was still high after 12 and 24 months in the field, it might indicate that these litter still contained large concentrations of lignin (Berg 2014). Personal observations confirm this, as tough and recalcitrant Vaccinium spp. leaves of different decomposition stages and ages could be found in the litter layer at our study site (see also Seeber and Seeber 2005). The lignin:N ratio was found to be a good predictor of decomposability in tropical and Mediterranean climates, but not in the temperate region (Aerts 1997). Even though we did not test it, our data may indicate that lignin:N could be a useful and applicable predictor for the decomposability of alpine litter.
Regarding the C:N ratios that can be used as index for palatability and litter quality (Amelung et al. 2018), we obtained relatively high values for all freshly fallen litter types ranging from 30 to almost 50. This would suggest the initial litter material to be generally of lowquality mainly due a relatively low total N versus a high total C content (Amelung et al. 2018). However, the SHRUB_Vv C:N ratios decreased only slowly along decomposition, therefore becoming only slowly more favourable to the decomposer community. On the other hand, the C:N ratios of perennial GRASS_Dg and FORB_Gs significantly dropped after six month, therefore becoming more favourable to decomposers already in the early stages.

Decomposition rates
The speed and rates of decomposition for SHRUB_Vv, GRASS_Dg, and FORB_Gs were in line with our expectations and the literature (Montané et al. 2010;Duan et al. 2013;Berg 2014). All the litter material degraded rapidly in the beginning followed by gradually decelerating decomposition. SHRUB_Vv as an evergreen plant showed the lowest decomposition rate, while the perennials GRASS_Dg and FORB_Gs degraded significantly faster. Compared with studies from the lowland, we observed lower decomposition rates most probably due to the harsher environmental conditions at our Alpine study site (Gavazov 2010). Generally, decomposition rates highly depend on climate (i.e. temperature and (soil) moisture), the composition of the litter material (i.e. nutritional quality and physical toughness), sitespecific local factors, and the decomposer community (Cornelissen and Thompson 1997;Berg and McClaugherty 2014). A fast decomposition can be observed when both temperature and moisture are high, as found in the tropics (Aerts 1997;Bradford et al. 2016). However, cold climatesfor example high-elevation and high-latitude environmentsonly partly fulfil these criteria, as they are characterised by high amounts of precipitation but low temperatures. Additionally, snow can persist for a long period and can return in summer months (Körner 2003). In these cold biomes, temperature, when not limited by soil moisture (Aerts 2006), was found to be the main driver of decomposition rates (Gavazov 2010;Bradford et al. 2016) along with other factors such as the decomposer community composition (Gessner et al. 2010;Kitz et al. 2015). As larger decomposers were excluded in this study, the observed decomposition rates can be accounted solely to the leaf-litter associated microbiome (Knapp et al. 2011).
Interestingly, litter decomposition was found to be quite high and active under the snow cover (Sjögersten and Wookey 2004;Baptist et al. 2010), due to the insulation of the snow pack that keep temperatures close to the freezing point. We recorded decomposition rates up to 42% LML after the first six months (Table 2) when the litterbags were covered by a 30 cm thick snow cover for most of the time, which is in line with the literature (Saccone et al. 2013).

Non-additive effects
Non-additive effects are defined as effects that are greater or less than the sum of two or more net (i.e. additive) effects (Hättenschwiler et al. 2005), the former being called synergistic and the latter antagonistic. We observed both non-additive effects, depending on the functional groups present, with an inconsistent pattern at the beginning of the decomposition processes and a synergistic trend after six months. This is consistent with Duan et al. (2013) who studied litter mixtures of up to 25 species in alpine meadows on the Tibetan Plateau and found antagonistic as well as synergistic effects in few-species mixtures, especially when shrubs were included. However, with an increasing number of litter types mixed with shrubs, mostly synergistic effects were recorded. In our study, GRASS_Dg showed the highest synergistic effects of all treatments followed by the SHRUB_Vv, hinting that adding at least one other functional litter type already increases the overall decomposition rate considerably (Fig. 2). The effect was even more pronounced for SHRUB_Vv and GRASS_Dg in the 3-species mixtures with a highly significant effect for GRASS_Dg + SF. Cuchietti et al. (2014) reported that high-quality litter can stimulate the decomposition of low-quality litter and vice-versa, hence benefit the decomposability of all interacting litter types. Our results for SHRUB_Vv and GRASS_Dg in mixtures fit into this scheme as both litter types significantly accelerated each other's decomposition (i.e. SHRUB_Vv + G, GRASS_Dg + S, and GRASS_Dg + SF). The reason might be a nutrient transfer from the high to the low-quality litter, leading to an elevated attraction of the species-specific microbiome (Hättenschwiler et al. 2005). Rief et al. (2012) found well separable communities of microorganisms inhabiting litter of eight alpine plant species from the same study site, which also included five species investigated from our study. A different picture was found for FORB_Gs, for which only adding SHRUB_Vv slightly increased decomposition (Table 2 and Fig. S6). In fact, we assume that in the first phase of decomposition GRASS_Dg and FORB_Gs negatively influence each other, the latter hampering decomposition in mixtures more severely. However, FORB_Gs had the highest mean decomposition rate of all tested litter types, therefore adding one or two litter types might slow down the decomposition processes due to a deficit of N and/or P (Zheng et al. 2017).
The finding that alpine low-quality dwarf shrub litter decomposition is accelerated when present in mixtures with high-quality litter types is of actual interest regarding ongoing land-use changes in Alpine pasturelands (Schirpke et al. 2013). These were shaped by extensive low-input agriculture for centuries, but are increasingly taken out of management due to socio-economic reasons (MacDonald et al. 2000), although being of high conservation status in the Habitats Directive (European Commission 1992). In sub-alpine pastureland, abandonment is followed by dwarf shrub encroachment, aggravated by global warming (Brigham et al. 2018), leading to high amounts of litter accumulation (Seeber and Seeber 2005;Steinwandter et al. 2017) and slow rates of decomposition (Montané et al. 2010), with mostly negative effects on various ecosystem services. However, these effects could be partly mitigated by the observed synergistic effects when low-quality litter occurs in mixtures with high-quality litter and by generally higher net-decomposition rates induced by global warming (Cornelissen et al. 2007;Luo et al. 2010;Duan et al. 2013).