Stable carbon and nitrogen isotopic composition of leaves , litter , and soils of various ecosystems along an elevational and land-use gradient at Mount Kilimanjaro , Tanzania

Variations in the stable isotopic composition of carbon (δ13C) and nitrogen (δ15N) of fresh leaves, litter, and topsoils were used to characterize soil organic matter dynamics of 12 tropical ecosystems in the Mount Kilimanjaro region, Tanzania. We studied a total of 60 sites distributed along five individual elevational transects (860– 4550 m a.s.l.), which define a strong climatic and land-use gradient encompassing semi-natural and managed ecosystems. The combined effects of contrasting environmental conditions, vegetation, soil, and management practices had a strong impact on the δ13C and δ15N values observed in the different ecosystems. The relative abundance of C3 and C4 plants greatly determined the δ13C of a given ecosystem. In contrast, δ15N values were largely controlled by land-use intensification and climatic conditions. The large δ13C enrichment factors (δClitter− δCsoil) and low soil C/N ratios observed in managed and disturbed systems agree well with the notion of altered SOM dynamics. Besides the systematic removal of the plant biomass characteristic of agricultural systems, annual litterfall patterns may also explain the comparatively lower contents of C and N observed in the topsoils of these intensively managed sites. Both δ15N values and calculated δ15N-based enrichment factors (δNlitter− δNsoil) suggest the tightest nitrogen cycling at high-elevation (>3000 m a.s.l.) ecosystems and more open nitrogen cycling both in grass-dominated and intensively managed cropping systems. However, claims about the nature of the N cycle (i.e. open or closed) should not be made solely on the basis of soil δ15N as other processes that barely discriminate against 15N (i.e. soil nitrate leaching) have been shown to be quite significant in Mount Kilimanjaro’s forest ecosystems. The negative correlation of δ15N values with soil nitrogen content and the positive correlation with mean annual temperature suggest reduced mineralization rates and thus limited nitrogen availability, at least in high-elevation ecosystems. By contrast, intensively managed systems are characterized by lower soil nitrogen contents and warmer conditions, leading together with nitrogen fertilizer inputs to lower nitrogen retention and thus significantly higher soil δ15N values. A simple function driven by soil nitrogen content and mean annual temperature explained 68 % of the variability in soil δ15N values across all sites. Based on our results, we suggest that in addition to land-use intensification, increasing temperatures in a changing climate may promote soil carbon and nitrogen losses, thus altering the otherwise stable soil organic matter dynamics of Mount Kilimanjaro’s forest ecosystems. Published by Copernicus Publications on behalf of the European Geosciences Union. 410 F. Gerschlauer et al.: Stable isotopes of C and N at Mount Kilimanjaro


Introduction
Conversion of natural ecosystems to agriculture is a worldwide phenomenon, which is of particular significance in tropical regions where human population growth rates are currently the highest (FAO and JRC, 2012).Changes in climate and land use significantly alter vegetation composition and biogeochemical cycles, causing a strong impact on carbon (C) and nitrogen (N) turnover and stocks (Smith et al., 2014).Tropical forest biomes are particularly relevant in this context, as they are significant C storages and N turnover hotspots (Bai et al., 2012;Hedin et al., 2009;Lewis et al., 2009;Pan et al., 2011;Vitousek, 1984).Considering the increasing pressure on natural land, it is becoming even more crucial to understand how anthropogenic interventions affect ecosystem C and N cycling and gain better knowledge about the main drivers of nutrient cycling, as well as the associated exchange processes with the atmosphere and hydrosphere in tropical environments.
Research exploiting the natural abundance of stable isotopes has proved quite suitable for investigating potential impacts of land-use and/or climate change on C and N cycling in terrestrial systems (Michener and Lajtha, 2007;Panettieri et al., 2017;Saiz et al., 2015a).Variations in the stable isotopic composition of C (δ 13 C) and N (δ 15 N) in plants and soils are the result of fractionation processes occurring during ecosystem exchange of C and N. Thus, δ 13 C and δ 15 N can serve as valuable indicators about ecosystem state and provide useful insights on how these systems respond to biotic and abiotic factors (Dawson et al., 2002;Högberg, 1997;Ma et al., 2012;Pardo and Nadelhoffer, 2010;Peterson and Fry, 1987;Robinson, 2001).
Plants discriminate against 13 CO 2 (carbon dioxide) during photosynthetic CO 2 fixation depending on plant metabolism (i.e.C 3 and C 4 photosynthetic pathways).Most tropical grasses typically employ the C 4 photosynthetic pathway (δ 13 C values > − 15 ‰), while trees and shrubs use the C 3 photosynthetic pathway (δ 13 C values < − 24 ‰) (Bird et al., 1994;Bird and Pousai, 1997;Cernusak et al., 2013;Farquhar et al., 1980).The distribution of C 3 and C 4 vegetation shows clear patterns along elevational gradients, with increasing abundance of C 3 species towards high elevations (Bird et al., 1994;Körner et al., 1991;Tieszen et al., 1979).Environmental conditions such as water availability also exert a significant influence on isotopic discrimination during atmospheric CO 2 fixation.Accordingly, compared to optimal moisture conditions, water stress leads to the enrichment of 13 C in C 3 plants (Farquhar and Sharkey, 1982), while this isotopic fractionation is less obvious or even absent in C 4 plants (Ma et al., 2012;Swap et al., 2004).
The soil organic matter (SOM) pool integrates the isotopic signature of the precursor biomass over different spatiotemporal scales (Saiz et al., 2015a).Variation in soil δ 13 C values represents a valuable tool to better assess SOM dynamics, mineralization processes, or reconstruct past fire regimes (Saiz et al., 2015a;Wynn and Bird, 2007).The δ 13 C of SOM in a given ecosystem is greatly controlled by the relative abundance of C 3 and C 4 plants due to their contrasting C isotopic composition.Therefore, strong variations in soil δ 13 C can also be used to identify sources of particulate organic matter and vegetation shifts such as woody thickening.However, fractionation effects associated with differential stabilization of SOM compounds, microbial reprocessing of SOM, soil physico-chemical characteristics, and the terrestrial Seuss effect preclude a straightforward interpretation of soil δ 13 C values (Saiz et al., 2015a).
Plant and soil δ 15 N relate to environmental and management conditions controlling N turnover, availability, and losses.δ 15 N values in soils are generally more positive than those of vegetation due to the relatively large isotopic fractionation occurring during soil N transformations (Dawson et al., 2002).The N cycle of a given ecosystem may be characterized as closed if both efficient microbial N retention and the absence of external N inputs (e.g.atmospheric deposition and fertilizer additions) prevent substantial gaseous and/or leaching N losses.In contrast, open ecosystem N cycling is characterized by significant inputs and losses of N. On the one hand, gaseous N losses from soils are strongly depleted in 15 N due to the high fractionation factors associated with these processes (Denk et al., 2017).This results in high δ 15 N values of the residual substrate, which consequently leaves less importance to impacts of external N additions (Robinson, 2001;Zech et al., 2011).On the other hand, N leaching seems to only discriminate slightly against ecosystem 15 N.According to Houlton and Bai (2009) δ 15 N values of drained water agree well with those of soils across various natural ecosystems worldwide.Moreover, it is also important to consider the possibility that soil δ 15 N may also be influenced by other factors including rooting depth, uptake of different N compounds, and symbiotic N 2 fixation (Nardoto et al., 2014).Variations in δ 15 N values of plants and soils have been successfully applied to characterize N cycling across a large variety of ecosystems worldwide (Amundson et al., 2003;Booth et al., 2005;Craine et al., 2015a, b;Martinelli et al., 1999;Nardoto et al., 2014).This includes research work that has particularly focused on the study of N losses derived from land-use changes or intensification (Eshetu and Högberg, 2000;Piccolo et al., 1996;Zech et al., 2011).
Information on ecosystem C and N cycling is still scarce in many tropical ecosystems, particularly in remote regions of Africa (Abaker et al., 2016;2018;Saiz et al., 2012;Townsend et al., 2011).Furthermore, feedbacks between C and N cycles, such as limitations of N availability in ecosystem C sequestration and net primary productivity in tropical forest, require urgent investigations (Gruber and Galloway, 2008;Zaehle, 2013).In such a context, the Kilimanjaro region in Tanzania offers the rare possibility to study a broad range of tropical ecosystems across contrasting landuse management intensities and varying climatic conditions.This region hosts a large variety of semi-natural and managed ecosystems as a result of the strong elevational and land-use gradient.
We hypothesized that (i) vegetation composition (C 3 / C 4 ) is the main control for ecosystem δ 13 C values, whereas (ii) δ 15 N values are rather controlled by land-use management and climatic conditions.The main aim of this study is to evaluate the potential of δ 13 C and δ 15 N values in plant and soil material to assess C and N cycling across a broad variety of semi-natural and managed ecosystems under varying climatic conditions.

Study sites
This study was conducted on the southern slopes of Mount Kilimanjaro (3.07 • S, 37.35 • E; 5895 m a.s.l.) in north-east Tanzania.The climate is characterized by a bimodal precipitation pattern with a major rainy season between March and May and the other peak between October and November.Recently, Appelhans et al. (2016) used a network of 52 meteorological stations strategically deployed in the Kilimanjaro region to measure air temperature and precipitation.They then used geo-statistical and machine-learning techniques for the gap filling of the recorded meteorological time series and their regionalization, which provides the means to calculate the meteorological data used for the complete set of sites (60) used in our work.Please refer to Appelhans et al. (2016) for more details.Maximum mean annual precipitation (MAP) of 2552 mm occurs at an elevation of around 2260 m a.s.l., decreasing towards lower and higher elevations, reaching 657 and 1208 mm yr −1 at 871 and 4550 m, respectively (Table 1).Variations in air temperature are dominated by diurnal rather than seasonal patterns (Duane et al., 2008).Mean annual temperature (MAT) decreases with increasing elevation, ranging from 24.8 • C at 860 m to 3.5 • C at 4550 m (Table 1).
Five altitudinal transects ranging from 860 to 4550 m a.s.l. were established along the mountain slopes.At each transect, 12 ecosystems occurring over a strong land-use gradient encompassing intensively managed cropping systems and semi-natural stands were investigated.Hence, the total number of plots studied was 60 (5 transects × 12 ecosystems; Table 1 and Fig. 1).The cropping systems comprised multilayer and multi-crop agroforestry home gardens (Hom), monoculture coffee plantations (Cof) with dispersed shading trees, and maize fields (Mai) subject to regular albeit moderate fertilizer and pesticide applications.Plant litter is regularly removed from Cof and Mai sites.Home gardens are manually ploughed, while combustion engine machinery is used for ploughing coffee plantations and maize fields.Coffee plantations are irrigated with drip irrigation systems.Both Hom and Cof sites still host indigenous forest trees that include Albizia schimperi, a species that may potentially fix atmospheric N.This is one of the five most abundant species in 1517 (54) 6.2 (0.6) Andosol 3.9 (0.2) 29.5 (5.1) 24.1 (6.2) 28.1 (2.4) 1.5 (0.1) 18.9 (0.7) Helichrysum vegetation (Hel) (S-N) alpine scrub vegetation 4250 (100) 1293 (31) 4.2 (0.4) Andosol 5.7 (0.3) 7.9 (1.4) 69.9 (9.5) 6.1 (3.3) 0.3 (0.2) 12.0 (1.1) Land uses are generically classified as managed (M) and semi-natural ecosystems (S-N).MAP and MAT stand for mean annual precipitation and temperature, respectively.Climatic values are according to Appelhans et al. (2016).Data represent mean values (n = 5 ± SE) for different ecosystems.The most representative soil type is shown for each ecosystem.Soil properties are given for topsoil (0-10 cm for pH and soil texture, 0-5 cm for soil organic carbon and total nitrogen).
www.biogeosciences.net/16/409/2019/Biogeosciences, 16, 409-424, 2019 two and four of the Hom and Cof sites, respectively, making up less than 25 % of the vegetation cover in all cases.Grasslands (Gra) and savannas (Sav) are extensively managed by means of domestic grazing and occasional grass cutting, thus having significantly lower anthropogenic disturbances than cropping systems.Semi-natural ecosystems include several montane forest stands.These include lower montane (Flm), Ocotea (Foc), Podocarpus (Fpo), Erica (Fer), and alpine shrub vegetation Helichrysum (Hel).Even though lower montane forests are currently under protection they are still subject to sporadic illegal logging.In addition to sampling undisturbed forest ecosystems of Ocotea and Podocarpus, we purposely studied sites that had been affected by logging activities and fire events prior to the establishment of the Kilimanjaro National Park (Soini, 2005): Ocotea (Fod) and Podocarpus (Fpd) (Table 1).Erica forests represent Africa's highest forests in the subalpine zone.Higher above is the alpine zone, the realm of Helichrysum vegetation that is dominated by cushion plants and tussock grasses (Ensslin et al., 2015;Hemp, 2005).Potential ecosystem productivity and decomposition rates show a hump-shaped pattern resembling that of precipitation (Supplement Fig. S1).It is interesting to see the close match between the two variables along the elevation range, although this trend weakens slightly towards higher-elevation sites.Optimum growth and decomposition conditions are shown between 1800 and 2500 m a.s.l.These locations correspond to low-altitude forest ecosystems (Flm and Foc) that do not experience severe seasonal limitations in moisture or temperature as is otherwise the case in lowerand higher-elevation systems that are moisture and temperature limited, respectively (Becker and Kuzyakov, 2018).Detailed physico-chemical characteristics of the dominant soils are listed in Table 1.Soils in the Mount Kilimanjaro region are mainly derived from volcanic rocks and ashes.The wide array of climatic conditions present along the elevational gradient influence soil genesis, which results in the occurrence of andosols at high elevations and soils of more advanced genesis at lower elevations (e.g.nitosols) (Majule, 2003).
It is extremely difficult to provide reliable estimates of both fertilizer and pesticide rates used in small household farms in sub-Saharan Africa.This is because the actual use of these products is strongly dependent on both its availability in the local and/or regional market, the economic circumstances of each individual farmer, and individual perceptions about their use (Saiz and Albrecht, 2016).The only sites receiving fertilizer are the two monocultures, maize (Mai) fields and coffee (Cof) plantations, and to a lesser extent the home garden (Hom) sites.In the latter sites Gütlein et al. (2018) report that weed control is mainly done by hand, and the use of mineral or organic N fertilizers is low or nonexistent.Extensively managed sites (i.e.Sav and Gra) receive varying amounts of organic inputs as a result of grazing activities, but again, their actual rates are unknown.A more detailed explanation on fertilizer and pesticide inputs used in the region is provided in the Supplement (Table S1).

Sampling and analyses
Fieldwork took place in February and March in 2011 and 2012.Sampling was conducted on 50 × 50 m plots established at each of the 60 studied sites (12 ecosystems × 5 transects).Surface litter and mineral topsoil (0-5 cm) were sampled at five locations (four corners and the central point) at each plot.Additionally, fresh mature leaves of the five most abundant plant species covering 80 % of total plant biomass per site were collected (Schellenberg Costa et al., 2017).All sampled materials (leaves, litter, and soil) were air-dried until constant weight, and leaf material was subsequently ovendried at 70 • C for 60 h prior to grinding.Soil was sieved to 2 mm with visible root fragments being further removed prior to grinding with a mixer mill (MM200, Retsch, Haan Germany).Soil pH was determined with a pH meter (MultiCal SenTix61, WTW, Weilheim, Germany) in a 0.01 M CaCl 2 solution, with a CaCl 2 -to-soil ratio of 2 : 1. Particle size distribution was determined gravimetrically using the pipette method (van Reeuwijk, 2002).
All soil, litter, and leaf samples were analysed with a dry combustion elemental analyser (Costech International S.p.A., Milano, Italy) fitted with a zero-blank autosampler coupled to a ThermoFinnigan DeltaPlus XL using continuous-flow isotope ratio mass spectrometry (CF-IRMS) for the determination of the abundance of elemental C and N and their stable isotopic composition (δ 13 C, δ 15 N).Precisions (standard deviations) on internal standards for elemental C and N concentrations and stable isotopic compositions were better than 0.08 % and 0.2 ‰, respectively.
Natural 13 C or 15 N abundances are expressed in δ units according to Eq. ( 1): where R sample denotes the ratio 13 C/ 12 C or 15 N/ 14 N in the sample, and R standard denotes the ratios in Pee Dee Belemnite or atmospheric N 2 (international standards for C and N, respectively).The average values for the plant samples were weighted considering their relative abundance at each site.Individual values for soil, litter, and leaves were averaged for each plot.
In addition, both δ 13 C-and δ 15 N-based enrichment factors (ε) were calculated following Eqs.( 2) and (3): ε N = δ 15 N litter − δ 15 N soil . (3) These were used as indicators for SOM decomposition dynamics and ecosystem N status (Garten et al., 2008;Mariotti et al., 1981).Note that we use the stable isotopic values of litter material rather than fresh leaves from various species to calculate enrichment factors, since litter provides a more unbiased representation of the quality, quantity, and spatiotemporal dynamics of organic inputs entering the SOM pool (Saiz et al., 2015a).

Statistical analysis
The normal distribution of the data was confirmed with the Shapiro-Wilk test.One-way ANOVA was performed to test for significant differences between ecosystems, while Tukey's HSD was used as a post hoc procedure to test for significant differences across sites (P ≤ 0.05).Correlation analyses were performed to identify soil, foliar, and climatic variables influencing soil δ 15 N values.Subsequently, a principal component analysis (PCA) was conducted to reveal relationships between the main variables affecting soil δ 15 N values.The PCA was based on a correlation matrix including soil (C and N concentrations, C/N ratio, δ 13 C, pH values, sand and clay contents) and climatic parameters (MAT and MAP).A stepwise multiple regression was used to identify the main driving parameters determining soil δ 15 N across the elevational transect.All statistical analyses were conducted with R (version 3.2.2;R Core Team, 2015).

General soil characteristics
Soil C and N contents were the highest in forest ecosystems and showed a decreasing trend towards managed sites (i.e. home gardens, grasslands, coffee, and maize fields) (Table 1).Also, natural savannas and Helichrysum ecosystems had lower soil C and N values compared to forest ecosystems.The low temperatures and sandy nature of the Helichrysum sites play a strong role in their characteristically low productivity and moderate decomposition potentials (Table 1; Fig. S1), which unquestionably affects the comparatively low soil C and N contents of these alpine systems.
An opposite trend to that of soil C and N abundance was observed for soil C/N ratios, whereby managed sites showed significantly lower values compared to those of semi-natural ecosystems.Soil pH values revealed acidic conditions at all sites, with the lowest values observed in forest sites having comparatively higher MAP (Table 1).

Variation of δ 13 C values along the elevational and land-use gradient
There were large variations in δ 13 C values along the elevational and land-use gradient, with distinct differences between managed and semi-natural ecosystems (Fig. 2).Compared to soils and litter, leaves invariably showed the lowest δ 13 C values in all the studied ecosystems, with the exception of grasslands and savannas that exhibited lower soil δ 13 C values than plant material.The δ 13 C values of semi-natural ecosystems ranged between −32.8 and −24.1 ‰ (mean ± SE: soil −26.0±0.2 ‰; litter −27.2 ± 0.2 ‰; leaves −29.3 ± 0.3 ‰), showing a progressive reduction with decreasing elevation (i.e. from 4500 to 1750 m a.s.l.; Fig. S2).The variation in δ 13 C values was much higher (−29.7 ‰ to −13.3 ‰) in managed ecosystems located at lower elevations (i.e. between 860 and 1750 m a.s.l.; Fig. S2).The highest δ 13 C values were observed in C 4 -dominated ecosystems (i.e.savannas, maize fields, and grasslands; soil −16.8 ± 0.6 ‰, litter −19.3 ± 0.8 ‰, leaves −18.8 ± 1.1 ‰), while lower δ 13 C values were obtained for coffee plantations and home gardens (soil −24.8±0.5 ‰, litter −27.2±0.4 ‰, leaves −27.3±0.4 ‰).Coffee plantations showed a slight influence of C 4 vegetation in the soil data as a result of grasses growing between the rows of coffee plants.No   Lower-case letters show significant differences between sampled materials within each ecosystem (one-way ANOVA followed by Tukey's HSD test as a post hoc procedure, P ≤ 0.05).The ecosystem acronyms used are as per Table 1.Mai, Cof, and Hom are managed cropping sites, Gra and Sav are extensively managed grasslands and savannas, and the rest represent semi-natural ecosystems.
Sites are ordered by increasing altitude.
litter and leaves in the ecosystems with a predominance of C 4 vegetation (savannas, maize fields, and grasslands).Exploratory data analyses revealed that in most cases soil, litter, leaf, and climatic variables cross-correlated with each other (Table S2 in the Supplement).
Figure 3 shows relatively small variations in δ 13 C enrichment factors (> − 1.25 ‰) both in undisturbed semi-natural and extensively managed sites along the elevational gradient, while managed and disturbed sites show higher and more variable δ 13 C enrichment factors.

Variation of δ 15 N values along the elevational and land-use gradient
Significantly higher δ 15 N values were observed for all sampled materials in the intensively managed (cropping) systems compared to semi-natural and grass-dominated ecosystems (Fig. 4a).The δ 15 N values for managed systems ranged between −2.6 and 7.8 ‰ (mean ± SE: soil 5.6 ± 0.3 ‰, litter 1.7 ± 0.5 ‰, leaves 2.0 ± 0.5 ‰).By contrast, seminatural ecosystems had considerably lower δ 15 N values, which ranged from −5.0 ‰ to 3.6 ‰ (soil 1.5 ± 0.2 ‰, litter  Capital letters indicate significant differences between ecosystems (one-way ANOVA followed by Tukey's HSD test as a post hoc procedure, P ≤ 0.05).The ecosystem acronyms used are the same as those in Table 1.Sites are ordered by increasing altitude.
−2.1 ± 0.2 ‰, leaves −1.3 ± 0.3 ‰).Soil δ 15 N values were significantly higher than those of leaves and litter across all the ecosystems studied, with the only exception being agroforestry home gardens (Fig. 4a).δ 15 N values of leaves and litter did not show significant differences within any given ecosystem.Calculated δ 15 N-based enrichment factors showed high variability across all ecosystems with values ranging from −7.5 ‰ to −1.6 ‰ (Fig. 4b).A differentiation between managed and natural ecosystems was less clear than for δ 15 N values.The most negative enrichment factors (< − 4.0 ‰) were observed for Helichrysum, Erica, Podocarpus disturbed, and grass-dominated ecosystems (savannas and grasslands).These enrichment factors were significantly less negative for montane forests at lower elevations (Podocarpus, Ocotea, and lower montane) and intensively managed (cropping) systems (i.e. home garden, coffee, and maize; Fig. 4b).

Impacts of soil and climatic variables on soil δ 15 N values
Two principal components (PCs) explained 78.3 % of the total soil δ 15 N variation (Fig. 5).The first component explained 55.8 % of the variability, and included soil chemistry and climatic variables (soil C and N concentrations, soil C/N ratio, soil pH, soil δ 13 C, MAP, and MAT).Highly significant correlations (P < 0.001) were obtained between PC 1 and the above factors (r = 0.93, 0.93, 0.61, −0.87, −0.76, 0.87, and −0.63, respectively; Table S3).The sec-ond component explained an additional 22.5 % of soil δ 15 N variability and included soil texture (clay and sand contents) and MAT.These variables were highly correlated with PC 2 (r = −0.84,0.82, and −0.65;Table S3).The principal component bi-plot showed a strong grouping between managed and semi-natural ecosystems (Fig. 5).Managed sites clustered around MAT, soil δ 13 C, and soil pH, while C 4 -dominated ecosystems (grassland, savannas, and maize fields) were preferentially influenced by the latter two variables.In contrast, semi-natural montane forest ecosystems rather grouped around soil chemical properties such as C and N contents, C/N ratio, and MAP, while alpine Helichrysum ecosystems clustered around soil sand content.
In addition to PCA, multiple regression analyses were performed using a stepwise procedure that identified soil N content and MAT as the main driving variables explaining the variation in soil δ 15 N. A paraboloid model explained 68 % of this variability (P < 0.05; Fig. 6).The combination of relatively high soil N contents (1 % to 3 %) and low MAT (up to 14 • C) invariably corresponded to low soil δ 15 N values (< 2 ‰) characteristic of semi-natural ecosystems.Conversely, the relatively high soil δ 15 N values (> 2 ‰) observed in managed ecosystems corresponded to low soil N contents (< 1 %) and comparatively high MAT (17 to 25 • C).
The relationship between soil δ 15 N values and climatic and edaphic variables provided valuable information about potentially different SOM dynamics in the various ecosystems studied, with data showing a clear differentiation between semi-natural and managed ecosystems (Fig. S4).The www.biogeosciences.net/16/409/2019/Biogeosciences, 16, 409-424, 2019 former is characterized by comparatively higher C/N ratios and lower δ 15 N values (averaging 15.5 ‰ and 1.5 ‰, respectively), while the latter showed lower C/N ratios and higher soil δ 15 N values (averaging 11.9 ‰ and 3.5 ‰, respectively).Managed ecosystems further grouped into intensively cropped (home gardens, maize fields, and coffee plantations) and extensively managed grass-dominated ecosystems (savannas and grasslands).

Factors influencing the variation of δ 13 C values along the elevational and land-use gradient
The δ 13 C values of leaves in C 3 -dominated (semi-natural) ecosystems at Mount Kilimanjaro increased with elevation (Figs. 1 and S2), which is in agreement with findings from other mountainous ecosystems in the tropics, Europe, and North America (Bird et al., 1994;Körner et al., 1991;Ortiz et al., 2016;Zhou et al., 2011;Zhu et al., 2009).The wider scatter of δ 13 C values observed in leaves relative to soils is most certainly due to the inherently large (interspecific and intraspecific) variability of δ 13 C in plants (Bird et al., 1994).
Different tissues within the plant can present widely divergent δ 13 C values as a result of fractionation processes associated with the C compounds involved in their construc- tion (Dawson et al., 2002).Moreover, other factors including light intensity, humidity, and the reutilization of previously respired low 13 C-CO 2 within the canopy may further contribute to the variability of δ 13 C in leaf tissues (Ometto et al., 2006;van der Merwe and Medina, 1989).
While fractionation effects preclude a straightforward interpretation of δ 13 C of SOM, this variable provides an integrated measure of the isotopic composition of the precursor biomass at the ecosystem level (Bird et al., 2004;Saiz et al., 2015a).Mass balance calculations that assume (i) 5 % (w/w) average root mass (< 2 mm) in soil samples and (ii) leaves having similar isotopic signals as roots show that the removal of visible sieved roots might cause a very small effect on soil isotopic values.This would amount to values ∼ 0.15 ‰ higher than the original soil isotopic values, with such discrepancy being even smaller if root samples were considered to have values 0.5 ‰-1 ‰ higher than leaves, as is commonly reported in the literature (calculations not shown).Besides the natural variability of soil δ 13 C values observed in C 3 -dominated semi-natural ecosystems, there were distinct patterns in the δ 13 C values of soil samples collected in extensively managed, low-elevation ecosystems where woody and grass vegetation coexist (i.e.grasslands and savannas), which indicates the strong influence exerted by C 4 vegetation on the C isotopic composition of all sampled materials (Fig. 2).The results obtained in semi-natural ecosystems at Mount Kilimanjaro fit well within the interpretative framework for ele-vational soil δ 13 C data proposed by Bird et al. (1994).These authors suggest that besides temperature and atmospheric pressure, other primary factors influencing soil δ 13 C values are the age and degree of decomposition of SOM and variables related to the characteristics of the canopy, including the proportion of respired CO 2 that is recycled during photosynthesis, the relative contribution of leaf and woody litter to SOM, and soil moisture.
Besides the factors explained above, soil δ 13 C values are strongly influenced by the balance between ecosystem C inputs and outputs.It seems reasonable to assume that in the case of natural ecosystems there may be a steady state between SOM inputs and decomposition rates.This should be in contrast with the typically altered nutrient dynamics of disturbed systems, particularly those under agricultural management (Wang et al., 2018).Low fractionation factors in δ 13 C are commonly reported between plant material and topsoils in natural systems, mainly because of the relatively limited humification of recent organic matter prevalent in topsoils (Acton et al., 2013;Wang et al., 2018).Thus, we hypothesized that if C inputs and outputs were roughly in balance, then the difference in δ 13 C values between plant material and topsoil would be smaller in undisturbed sites compared to managed or disturbed sites.The results shown in Fig. 3 agree well with this notion.
Soil δ 13 C values decreased with increasing MAP and decreasing MAT, which also corresponded with higher SOC contents (Fig. S3).This suggests that the relatively cooler and wetter conditions of high-elevation semi-natural forest ecosystems (i.e.Foc, Fpo) promote the accumulation of SOM, which is similar to previous findings of work conducted along elevational gradients (Bird et al., 1994;Kohn, 2010).Compared to high-elevation locations, the climatic conditions of mid-elevation ecosystems are more favourable for the activities of SOM decomposers, as these sites are consistently warmer and drier than the characteristically cool and occasionally waterlogged high-altitude ecosystems (Fig. S1; Becker and Kuzyakov, 2018;Borken and Matzner, 2009;Garten et al., 2009;Kirschbaum, 1995;Leirós et al., 1999).The comparatively high soil δ 13 C values observed in the disturbed Podocarpus (Fpd) and Erica forest (Fer) plots may have been partly caused by recurrent fire events (Hemp, 2005) leading to reduced SOC contents and higher C/N ratios (Saiz et al., 2015a).Further variations in soil δ 13 C values could also be related to the biochemical composition of the precursor biomass.For instance, herbaceous vegetation is pervasive at high elevations and contains relatively low amounts of lignin -an organic compound characteristically depleted in 13 C (Benner et al., 1987).This may contribute to explaining the higher δ 13 C values observed in plant and soil materials in alpine ecosystems dominated by Helichrysum vegetation compared to forest ecosystems at lower elevations (Fig. 2).
Elevation also has a strong influence on the seasonal litterfall dynamics observed in Mount Kilimanjaro and thus may have significant implications in the SOM cycling across the various ecosystems.Becker et al. (2015) suggest that the large accumulation of particulate organic matter observed at the end of the dry season in low-and mid-altitude ecosystems may result in the increased mineralization of easily available substrates (Mganga and Kuzyakov, 2014) and nutrient leaching (Gütlein et al., 2018) during the following wet season.Agricultural practices such as the removal of biomass or ploughing deplete SOM, particularly in intensively managed systems (i.e.maize, home gardens, and coffee plantations), thus leading to lower SOC contents and C/N ratios and slightly higher soil δ 13 C values than those observed in semi-natural ecosystems at comparable elevations (e.g.lower montane forests; Fig. S3).Indeed, the relationship between δ 13 C enrichment factors and soil C/N ratios shown in Fig. 3 is quite informative regarding SOM dynamics.As previously mentioned, soil C/N ratios provide a good indication of SOM decomposition processes, typically showing comparatively low values in managed and disturbed systems.These correspond well with sites having large enrichment factors (< − 1.25 ‰; i.e. intensively managed and disturbed sites), which agrees with the notion of altered SOM dynamics.Therefore, besides the systematic removal of plant biomass characteristic of agricultural systems, annual litterfall patterns may also explain the comparatively lower contents of C and N observed in the topsoils of intensively managed sites (Table 1; Figs.S3, S4).Moreover, low-elevation ecosystems contain a variable mixture of C 3 and C 4 vegetation, which have been shown to have differential mineralization dynamics as demonstrated by incubation experiments (Wynn and Bird, 2007) and field-based research (Saiz et al., 2015a).
Our data show strong relationships between temperature and variables directly related to SOM dynamics such as soil δ 13 C, C, N, and C/N ratios (Table S2).These results agree well with recent findings by Becker and Kuzyakov (2018), who studied SOM decomposition dynamics at these very sites.An important finding revealed by that study is seasonal variation in temperature being a major factor controlling litter decomposition.Their study shows that small seasonal variations in temperature observed at high-elevation sites exert a strong effect on litter decomposition rates.Therefore, the authors argue that the projected increase in surface temperature may result in potentially large soil C losses at these sites due to the comparatively strong temperature sensitivity to decomposition that is commonly observed at low temperatures and at high-elevation sites (Blagodatskaya et al., 2016).
Savannas and grasslands are subject to recurrent fire events, and thus the soils of these ecosystems may potentially contain significant amounts of fire-derived (pyrogenic) C (Saiz et al., 2015b).This can be partly demonstrated by the higher soil C/N ratios observed in these ecosystems compared to C 4 -dominated agricultural systems protected from fire (e.g.maize plantations; Fig. S3d).Moreover, the δ 13 C values of soils in grasslands and savannas were lower than www.biogeosciences.net/16/409/2019/Biogeosciences, 16, 409-424, 2019 those of leaves, which may be due to the savanna isotope disequilibrium effect (SIDE) (Bird and Pousai, 1997;Saiz et al., 2015b).The latter concept explains the difference in C isotopic composition between the precursor vegetation and pyrogenic C compounds produced during the combustion of biomass.Saiz et al. (2015b) have demonstrated that savanna fires produce pyrogenic C that is relatively 13 C depleted with respect to the precursor biomass.Furthermore, the combustion of C 4 vegetation produces finer pyrogenic C particles than woody biomass, resulting in the preferential export of grass-derived pyrogenic particles from the site of burning, which further enhances the depletion of 13 C in these soils (Saiz et al., 2018).

Variation of δ 15 N values along the elevational and land-use gradient
The δ 15 N values of leaves, litter, and topsoil presented here (Fig. 4a) agree well with the range of data reported from earlier investigations in the same study region (Amundson et al., 2003;Zech et al., 2011), but with our study involving more ecosystems, replicate sites, and a far larger spatial sampling domain.Overall, the δ 15 N values for montane tropical forest ecosystems in Mount Kilimanjaro are considerably lower than the mean values reported for a broad variety of tropical lowland forests worldwide (soil values ranging from 3 ‰ to 14 ‰; de Freitas et al., 2015;Martinelli et al., 1999;Nardoto et al., 2014;Piccolo et al., 1996;Sotta et al., 2008).Rather, the δ 15 N values observed in the montane forests investigated are in the same range of temperate forest ecosystems reported in a comprehensive literature review by Martinelli et al. (1999).These authors argue that, compared to tropical lowland forests, the lower δ 15 N values of temperate and montane tropical forests result from their lower N availability and thus lower ecosystem N losses.However, this hypothesis may not completely hold for the montane forest ecosystems of our study, since Gütlein et al. (2018) reported elevated soil NO − 3 and DON concentrations at deep soil solution (80 cm) and significant nitrogen leaching rates of 10-15 kg N ha −1 yr −1 .The relatively low δ 15 N-based enrichment factors observed in the lower montane, Ocotea, and undisturbed Podocarpus forest (Fig. 4b) were probably due to the prevalence of biological dinitrogen fixation (BNF) at these ecosystems.The assumption of significant BNF is supported by leaf δ 15 N values close to 0 ‰ (Fig. 4a) and is in line with previous works (Craine et al., 2015a;Nardoto et al., 2014;Robinson, 2001).Furthermore, sporadic measurements of N compounds in rainfall and throughfall conducted at our forest sites showed substantial input of N via atmospheric deposition, which may be of the order of N leaching losses (unpublished results).This agrees well with findings from Bauters et al. (2018) reporting 18 kg N ha −1 yr −1 N inputs via wet deposition into tropical forests of the Congo Basin, which are predominantly derived from biomass burning and long-range atmospheric transport.High N inputs into these forest ecosystems are likely to be in a similar range as N outputs (prevailed by leaching losses, particularly where MAP is highest; Gütlein et al., 2018), and therefore they would not translate to strong effects on ecosystem δ 15 N values.The significantly more negative enrichment factors observed in the disturbed Podocarpus and Erica forests (Fig. 4b) may be related to past fire events (Hemp, 2005;Zech et al., 2011).Burning of vegetation may cause losses of 15 N-depleted NO χ gas and N leachate, resulting in higher soil δ 15 N values, thus producing variations in δ 15 N-based enrichment factors (Zech et al., 2011).
Previous studies have shown that δ 15 N values generally increase with land-use intensification (Martinelli et al., 1999;Stevenson et al., 2010), which corresponds well with the more positive δ 15 N values observed in the intensively managed agricultural systems occurring at the mountain's foot slope (Fig. 4a).Indeed, agronomic practices such as fertilization, the removal of plant material after harvest, and ploughing are factors known to impact N turnover processes that strongly affect δ 15 N values (Bedard-Haughn et al., 2003;Saiz et al., 2016).However, our values are in the lower range of published data for other land-use gradients (Aranibar et al., 2008;Eshetu and Högberg, 2000;Traoré et al., 2015) and may partly be the result of comparably low to moderate organic and inorganic N fertilization rates currently applied in the region (anecdotal evidence gathered by the authors and found in the Supplement).Additionally, the nitrogen isotopic signal of mineral fertilizers commonly used in the region is ∼ 0 ‰ (Bateman and Kelly, 2007), and thus it may not exert a significant additional bias on the interpretation of soil δ 15 N values.However, the addition of manure (δ 15 N ∼ 8 ‰) in Hom systems, although used in low quantities (Gütlein et al., 2018), may well have contributed to the high δ 15 N values observed in this ecosystem (Fig. 4).Also, we suggest that the use of pesticides may not pose a strong bias in our isotopic results since their use is limited to intensively managed sites, and the actual isotopic values of pesticides work in the opposite direction to the observed data (Fig. 4; Supplement).
Compared to other low-elevation managed stands such as home gardens and coffee plantations, the higher δ 15 N-based enrichment factors observed in maize fields and in grassdominated ecosystems (grasslands and savannas) (Fig. 4b) may be related to both the organic inputs resulting from grazing activities and the influence of C 4 vegetation.Both Aranibar et al. (2008) and Wang et al. (2010) have suggested that variations in δ 15 N values within a given ecosystem could be due to C 3 and C 4 plants preferentially absorbing chemical forms of N with differing 15 N abundances.Moreover, recurrent fires characteristic of tropical grasslands and savannas may have also influenced their comparatively high soil δ 15 N, causing the relatively high δ 15 N-based enrichment factors.

Factors controlling soil δ 15 N along the elevational and land-use gradient
The strong controlling effects exerted by climatic and edaphic factors on soil δ 15 N values agree well with numerous previous works (Amundson et al., 2003;Conen et al., 2013;Eshetu and Högberg, 2000;Martinelli et al., 1999;Stevenson et al., 2010).The principal component analysis of factors controlling soil δ 15 N revealed a strong clustering between managed and semi-natural ecosystems (Fig. 5), which was also reflected in the multiple regression analysis and graphical representation depicting soil δ 15 N as a function of soil N concentration and MAT (Fig. 6).Semi-natural ecosystems were characterized by relatively low soil δ 15 N values and occurred across a broad range of soil N contents in locations with low to medium MAT.By contrast, intensively managed ecosystems had higher soil δ 15 N values and corresponded to locations with low soil N contents and high MAT.The negative correlation of δ 15 N values with soil nitrogen content and the positive correlation with mean annual temperature suggest reduced mineralization rates and thus limited nitrogen availability, at least in high-elevation ecosystems.The sharp contrast observed both in soil C/N ratios and δ 15 N values between managed and semi-natural ecosystems offers additional useful information about their potentially contrasting SOM dynamics (Fig. S4d).Intensively managed sites consistently showed low soil C/N ratios and high soil δ 15 N values, which may initially suggest a more open N cycle and potentially greater N losses as reported by Gerschlauer et al. (2016) for some of these ecosystems.This may be due to the C limitation of heterotrophic microbial N retention under low C/N ratios (Butterbach-Bahl and Dannenmann, 2012).However, nitrate leaching is quite a relevant process that discriminates only slightly against 15 N (Denk et al., 2017), which may confound the interpretation of soil δ 15 N values.Indeed, Gütlein et al. (2018) have recently shown that nitrate leaching may be quite significant in Mount Kilimanjaro's semi-natural forests.Therefore, at least in these ecosystems, claims about the nature of the N cycle (i.e.open or closed) should not be made solely on the basis of soil δ 15 N.
Grass-dominated ecosystems (grasslands and savannas) were noticeably different to the intensively managed croplands, as demonstrated by the higher soil C/N ratios and lower soil δ 15 N of the former, which suggests a lower degree of decomposition of organic matter and potentially lower N turnover rates (Saiz et al., 2016).Within the intensively managed sites, the stands under maize cultivation show an interesting case of enhanced SOM dynamics.These sites are under an intensive management regime that involves the removal of aboveground vegetation after harvest.This fact combined with the faster decomposition rates reported for C 4 -derived SOM (Saiz et al., 2015a(Saiz et al., , 2016;;Wynn and Bird, 2007) may invariably lead to their characteristically low SOC and N contents (Table 1; Figs.S3, S4).Furthermore, low soil C/N ratios have been reported to enhance gaseous losses in semi-arid systems, which leads to increased soil δ 15 N values (Aranibar et al., 2004) and may explain why maize stands showed the highest soil δ 15 N values of all the land uses studied.
Semi-natural ecosystems showed rather high soil C/N ratios and low soil δ 15 N values compared to managed sites (Fig. S4d).The more humid and cooler conditions prevalent in forest ecosystems may limit decomposition processes, thereby contributing significantly to their higher SOM abundance (Table 1).A small variation range in soil δ 15 N values was also reported by Zech et al. (2011) for semi-natural ecosystems (Foc and Fpo) when working along the same land-use and elevation gradient.Like us, these authors also observed a strong significant correlation of soil δ 15 N with MAT, but not with MAP (Table S2).Additionally, sitespecific soil characteristics and the structural composition of vegetation have a strong influence on ecosystem nutrient dynamics (Saiz et al., 2012(Saiz et al., , 2015a)).Ecosystem disturbances (e.g.fire, selective logging, etc.) cause changes in vegetation cover that affect SOM cycling and may translate into variations in soil C/N ratios (Saiz et al., 2016).Both Ocotea and Podocarpus forests contain disturbed (Fod, Fpd) and undisturbed stands (Foc, Fpo), though only the Podocarpus ecosystems allow for a general overview of disturbance impacts on SOM-related properties.While changes in the isotopic composition of C and N were not significant, soil C/N ratios were heavily influenced by disturbance (Fig. S4).Compared to non-disturbed sites, the lower C and N contents observed in the soil of disturbed ecosystems indicate reduced OM inputs to the soil and/or enhanced decomposition of SOM (Table 1).The higher soil C/N ratios observed in the Podocarpus disturbed and Erica forests may well be the result of fire, which may preferentially promote N losses while accruing relatively recalcitrant C forms (i.e.pyrogenic C).Woody biomass combustion produces pyrogenic C that accumulates preferentially close to the site of production (Saiz et al., 2018), thus likely contributing to the higher soil C/N ratios observed at these disturbed ecosystems.The lowest soil C/N ratios among all semi-natural ecosystems were observed at the alpine Helichrysum sites, which may relate to their characteristically sparse vegetation and extremely low MAT.Under such circumstances, soil development, biomass inputs, decomposition processes, and thus soil N turnover may be strongly limited, as was confirmed by a recent study conducted at one of these sites (Gütlein et al., 2017).

Conclusions
The variations in δ 13 C and δ 15 N values combined with the interpretation of other indices, such as δ 13 C-and δ 15 N-based enrichment factors and soil C/N ratios, enabled a qualitative characterization of regional differences in C and N dynam- Our data show that SOM contents are higher in cold and wet high-elevation ecosystems than at low-elevation managed sites.Management practices such as tillage, harvest, and vegetation burning promote the loss of OM, with SOM decomposition being further enhanced by the warm and moderately wet conditions of the mountain's foot slope.Based on our results, we suggest that besides management, increasing temperatures in a changing climate may promote C and N losses, thus altering the otherwise stable SOM dynamics of Mount Kilimanjaro's forest ecosystems.Moreover, the current situation of low N inputs in managed systems of sub-Saharan Africa is likely to change, since national efforts that aim to increase fertilizer use are currently at < 10 % of recommended rates (Hickman et al., 2014).Therefore, our data may also be valuable as a generic reference for low-elevation tropical agrosystems managed under low N inputs, while it may also allow for the monitoring of expected changes in agricultural management and associated impacts on the ecosystem N cycle through the study of variation in δ 15 N values.
In addition to climatic and edaphic factors, δ 15 N values of plant and soil material can largely depend on both the amount δ 15 N signal of atmospheric deposition and BNF, which highlights the importance of conducting additional measurements of site-specific N cycling when comparing ecosystem δ 15 N values across different biomes and regions.The combination of qualitative isotope natural abundance studies at a large number of sites (this study) with more elaborated quantitative process studies using enriched isotope labelling and N losses on a lower number of selected sites represents an ideal approach to characterize the ecosystem C and N cycling of the larger Mount Kilimanjaro region with its diverse ecosystems, climate, and management.
Author contributions.FG contributed to design, performed the study, and co-wrote the paper; GS contributed to analyses and cowrote the paper; DSC and MK provided plant samples and contributed to writing; MD contributed to writing; and RK designed the study and contributed to analyses and writing.
Competing interests.The authors declare that they have no conflict of interest.

Figure 1 .
Figure 1.Geographical distribution of investigated ecosystems: (a) along the elevational and land-use gradient.MAP denotes mean annual precipitation and MAT mean annual temperature.The colours of the boxes framing the ecosystem names match the colours of symbols in the GeoTIFF panel below; (b) along the southern slope of Mount Kilimanjaro.Symbols represent individual ecosystems (12) replicated five times (60 study sites in total).

Figure 2 .
Figure 2. Variation in δ 13 C values for leaves, litter, and soil along the Kilimanjaro elevational and land-use gradient.Ecosystem data represent the average values of five sites (one per each transect), with each site being composed of five samples (n = 5).Box plots show median values per ecosystem with whiskers representing the first and third quartiles.Dots represent outliers.The shaded region represents managed ecosystems (both intensively and extensively), while those unshaded indicate semi-natural ecosystems.Lower-case letters show significant differences between sampled materials within each ecosystem (one-way ANOVA followed by Tukey's HSD test as a post hoc procedure, P ≤ 0.05).The ecosystem acronyms used are as per Table1.Mai, Cof, and Hom are managed cropping sites, Gra and Sav are extensively managed grasslands and savannas, and the rest represent semi-natural ecosystems.Sites are ordered by increasing altitude.

Figure 3 .
Figure 3. (a) Variation in δ 13 C-based enrichment factors (δ 13 C itter-soil ) with elevation; (b) relationship between δ 13 C-based enrichment factors (δ 13 C itter-soil ) and SOC concentration (log SOC); and (c) relationship between δ 13 C-based enrichment factors (δ 13 C itter-soil ) and soil C/N ratios.Note: a savanna site with large C 3 influence was removed from the figure for clarity.

Figure 4 .
Figure 4. Variation in δ 15 N values and δ 15 N-based enrichment factors along the Kilimanjaro elevational and land-use gradient.(a) Variation in δ 15 N values for leaves, litter, and soil material sampled along the Kilimanjaro elevational and land-use gradient.Box plots show median values per ecosystem with whiskers representing the first and third quartiles.Dots represent outliers.Ecosystem data represent the average values of five sites (one per each transect), with each site being composed of five samples.Lower-case letters show significant differences between sampled materials within each ecosystem (one-way ANOVA followed by Tukey's HSD test as a post hoc procedure, P ≤ 0.05).(b) Variation in δ 15 N-based enrichment factors (δ 15 N litter-soil ) calculated for the different ecosystems along the elevational and land-use gradient.Dotted box plots indicate ecosystems dominated by C 4 vegetation.Capital letters indicate significant differences between ecosystems (one-way ANOVA followed by Tukey's HSD test as a post hoc procedure, P ≤ 0.05).The ecosystem acronyms used are the same as those in Table1.Sites are ordered by increasing altitude.

Figure 5 .
Figure 5. Principal component analysis bi-plot for soil and climate variables potentially controlling soil δ 15 N. Symbols are as per all previous figures.Acronyms are as per Table 1.C/N: soil C/N ratio, C: soil carbon content, N: soil nitrogen content, MAP: mean annual precipitation, clay: soil clay content, MAT: mean annual temperature, δ 13 C: soil δ 13 C, and pH: soil pH.

Table 1 .
General characteristics of ecosystems investigated at Mount Kilimanjaro, Tanzania.