Convergence of soil nitrogen isotopes across global climate gradients

Stable isotopes as a predictor for organic or conventional classification of berries and vegetables

Organic agriculture is expanding worldwide, driven by expectations of improving food quality and soil health. However, while organic certification by regulatory bodies such as the United States Department of Agriculture and the European Union confirms compliance with organic standards that prohibit synthetic chemical inputs, there is limited oversight to verify that organic practices, such as the use of authentic organic fertilizer sources, are consistently applied at the field level. This study investigated the elemental content of carbon (C) and nitrogen (N) and their stable isotopes (δ13C and δ15N) in seven different crops grown under organic or conventional practices to assess their applicability as a screening tool to verify the authenticity of organic labeled produce. Holm corrected Welch t-tests and a generalized linear mixed model (GLMM) were used to assess the potential of stable isotope or crop elemental content to differentiate organic vs. conventional production systems. Total C and N content or C/N ratio was not significantly different between production systems or among geographic origins for most crops. However, the average N stable isotope (δ15N) content differed, with conventional crops at 1.8 ±  2.2‰ and organic at 6.0 ±  3.4‰. A mixed model incorporating elemental contents and stable isotopes identified δ15N as the primary predictor in discriminating organic and conventional production systems. A δ15N threshold is suggested to differentiate conventional from organic grown raspberries (δ15N <  2.17‰) and strawberries (δ15N <  3.22‰), for an estimated false negative rate of 1%. Although further evaluation is needed, our extensive dataset (n = 791) captures key predictors of agricultural production systems and holds potential as a benchmark for future organic production verification.

Stable isotope analysis in soil prospection reveals the type of historic land-use under contemporary temperate forests in Europe

The determination of δ13C and δ15N values is a common method in archaeological isotope analysis-in studying botanical and human remains, dietary practices, and less typically soils (to understand methods of agricultural cultivation, including fertilization). Stable isotope measurements are also commonly used in ecological studies to distinguish different ecosystems and to trace diachronic processes and biogeochemical mechanisms, however, the application of this method in geochemical prospection, for determining historic land-use impact, remains unexplored. The study at hand focuses on a deserted site of a Cistercian manor, dating from the thirteenth to fifteenth centuries. Isotopic measurements of anthropogenically influenced soils have been compared to approximately 400 archaeobotanical, soil, and sediment samples collected globally. The results reveal the potential of isotope measurements in soil to study the impact of past land use as isotope measurements identify specific types of agricultural activities, distinguishing crop production or grazing. δ13C and δ15N ratios also likely reflect fertilization practices and-in this case-the results indicate the presence of cereal cultivation (C3 cycle plants) and fertilization and that the site of the medieval manor was primarily used for grain production rather than animal husbandry.

Soil δ15N spatial distribution is primarily shaped by climatic patterns in the semiarid Caatinga, Northeast Brazil

Soil nitrogen isotopic composition (δ15Nsoil) is an invaluable tool as it integrates nitrogen (N) transformations in soils. In addition to serving as a baseline to understand the N cycle, spatial representations of δ15Nsoil across landscapes (or isoscapes) is a multi-purpose tool useful to investigate, for example, plant-microbe interactions, animal migration and forensics. We investigate the climatic and edaphic controls of δ15Nsoil utilising data from 29 geographical locations sampled across the semiarid Brazilian Caatinga biome. The sampling covered a mean annual precipitation (PA) gradient ranging from 0.51 to 1.36 m a-1 and eight soil types originating from three different geological origins. Our data show that the combination of higher aridity and lower seasonality (ψ) leads to higher values of δ15Nsoil. Moreover, soil total carbon had a positive relationship with δ15Nsoil, appearing within the best-supported models according to the information-theoretic approach undertaken here. The contribution to the plant communities by the Fabaceae trees expressed as their basal area was not related to δ15Nsoil values, suggesting that the magnitude of biological N fixation in the Caatinga is not large enough to be reflected in the soil. In addition, considering PA in a categorical fashion, i.e., 'high' (> 0.8 m a-1) and 'low' PA (< 0.8 m a-1), we found that, within the wetter category, δ15Nsoil was positively related to several soil properties (i.e., clay content, effective cation exchange capacity, exchangeable calcium, silt content, pHH2O, total phosphorus and sum of bases) and negatively related to sand content. Our study provides new insights into the functioning of semiarid ecosystems from a pedo-isotopic perspective and contributes to the overall understanding of the N cycle in the Caatinga region, with the potential to support the development of new conceptualisation of biogeochemical process and testing of global models that simulate N and C cycles.

Patterns of δ15N in forest soils and tree foliage and rings between climate zones in relation to atmospheric nitrogen deposition: A review

The stable nitrogen (N) isotope ratio (δ15N) of forest samples (soils, tree foliage, and tree rings) has been used as a powerful indicator to explore the responses of forest N cycling to atmospheric N deposition. This review investigated the patterns of δ15N in forest samples between climate zones in relation to N deposition. Forest samples exhibited distinctive δ15N patterns between climate zones due to differences in site conditions (i.e., N availability and retention capacity) and the atmospheric N deposition characteristics (i.e., N deposition rate, N species, and δ15N of deposited N). For example, the δ15N of soil and foliage was higher for tropical forests than for other forests by >1.2 ‰ and 4 ‰, respectively due to the site conditions favoring N losses coupled with relatively low N deposition for tropical forests. This was further supported by the unchanged or increased δ15N of tree rings in tropical forests, which contrasts with other climate zones that exhibited a decreased wood δ15N since the 1920s. Subtropical forests under a high deposition of reduced N (NHy) had a lower δ15N by 2-5 ‰ in the organic layer compared with the other forests, reflecting high retention of 15N-depleted NHy deposition. At severely polluted sites in East Asia, the decreased δ15N in wood also reflected the consistent deposition of 15N-depleted NHy. Though our data analysis represents only a subset of global forest sites where atmospheric N deposition is of interest, the results suggest that the direction and magnitude of the changes in the δ15N of forest samples are related to both atmospheric N and site conditions particularly for tropical vs. subtropical forests. Site-specific information on the atmospheric N deposition characteristics would allow more accurate assessment of the variations in the δ15N of forest samples in relation to N deposition.

A distinctive latitudinal trend of nitrogen isotope signature across urban forests in eastern China

Rapid urbanization has greatly altered nitrogen (N) cycling from regional to global scales. Compared to natural forests, urban forests receive much more external N inputs with distinctive abundances of stable N isotope (δ15 N). However, the large-scale pattern of soil δ15 N and its imprint on plant δ15 N remain less well understood in urban forests. By collecting topsoil (0-20 cm) and leaf samples from urban forest patches in nine large cities across a north-south transect in eastern China, we analyzed the latitudinal trends of topsoil C:N ratio and δ15 N as well as the correlations between tree leaf δ15 N and topsoil δ15 N. We further explored the spatial variation of topsoil δ15 N explained by corresponding climatic, edaphic, vegetation-associated, and anthropogenic drivers. Our results showed a significant increase of topsoil C:N ratio towards higher latitudes, suggesting lower N availability at higher latitudes. Topsoil δ15 N also increased significantly at higher latitudes, being opposite to the latitudinal trend of soil N availability. The latitudinal trend of topsoil δ15 N was mainly explained by mean annual temperature, mean annual precipitation, and atmospheric deposition of both ammonium and nitrate. Consequently, tree leaf δ15 N showed significant positive correlations with topsoil δ15 N across all sampled plant species and functional types. Our findings reveal a distinctive latitudinal trend of δ15 N in urban forests and highlight an important role of anthropogenic N sources in shaping the large-scale pattern of urban forest 15 N signature.

Overestimated nitrogen loss from denitrification for natural terrestrial ecosystems in CMIP6 Earth System Models

Denitrification and leaching nitrogen (N) losses are poorly constrained in Earth System Models (ESMs). Here, we produce a global map of natural soil ¹⁵N abundance and quantify soil denitrification N loss for global natural ecosystems using an isotope-benchmarking method. We show an overestimation of denitrification by almost two times in the 13 ESMs of the Sixth Phase Coupled Model Intercomparison Project (CMIP6, 73 ± 31 Tg N yr^-1), compared with our estimate of 38 ± 11 Tg N yr^-1, which is rooted in isotope mass balance. Moreover, we find a negative correlation between the sensitivity of plant production to rising carbon dioxide (CO₂) concentration and denitrification in boreal regions, revealing that overestimated denitrification in ESMs would translate to an exaggeration of N limitation on the responses of plant growth to elevated CO₂. Our study highlights the need of improving the representation of the denitrification in ESMs and better assessing the effects of terrestrial ecosystems on CO₂ mitigation.

Nitrogen palaeo-isoscapes: Changing spatial gradients of faunal δ15N in late Pleistocene and early Holocene Europe

Nitrogen isotope ratio analysis (δ15N) of animal tissue is widely used in archaeology and palaeoecology to investigate diet and ecological niche. Data interpretations require an understanding of nitrogen isotope compositions at the base of the food web (baseline δ15N). Significant variation in animal δ15N has been recognised at various spatiotemporal scales and related to changes both in baseline δ15N, linked to environmental and climatic influence on the terrestrial nitrogen cycle, and animal ecology. Isoscapes (models of isotope spatial variation) have proved a useful tool for investigating spatial variability in biogeochemical cycles in present-day marine and terrestrial ecosystems, but so far, their application to palaeo-data has been more limited. Here, we present time-sliced nitrogen isoscapes for late Pleistocene and early Holocene Europe (c. 50,000 to 10,000 years BP) using herbivore collagen δ15N data. This period covers the Last Glacial-Interglacial Transition, during which significant variation in the terrestrial nitrogen cycle occurred. We use generalized linear mixed modelling approaches for interpolation and test models which both include and exclude climate covariate data. Our results show clear changes in spatial gradients of δ15N through time. Prediction of the lowest faunal δ15N values in northern latitudes after, rather than during, the Last Glacial Maximum is consistent with the Late Glacial Nitrogen Excursion (LGNE). We find that including climatic covariate data does not significantly improve model performance. These findings have implications for investigating the drivers of the LGNE, which has been linked to increased landscape moisture and permafrost thaw, and for understanding changing isotopic baselines, which are fundamental for studies investigating diets, niche partitioning, and migration of higher trophic level animals.

Climatic and edaphic controls over soil δ15N in temperate grassland of northern China: A PLS-PATH analysis

Identifying the impact path of climate and soil factors on soil δ15N is very crucial for better understanding the N turnover in soils and the integrated information about ecosystem N cycling. Many studies have showed that climate and soil variables influence the change of soil δ15N. However, most of the existing studies focused on the overall impact of factor on soil δ15N, without distinguishing between the direct and indirect effect. Although scholars have studied the relationships among temperature, precipitation, soil N, soil pH, and soil δ15N rather than estimating all the causal relationships simultaneously. To answer the above-mentioned questions, a regional-scale soil collection was conducted across a temperate grassland in northern China. Meanwhile, a PLS-PATH analysis was utilized to evaluate the direct and indirect effects of various factors on soil δ15N and to explore the causal relationships among variables. The results showed that along the transect, mean annual precipitation (MAP) and mean annual temperature (MAT) directly and significantly reduced soil δ15N, and indirectly affected soil δ15N through their effects on soil pH, soil clay, soil N and soil C/N. Soil C/N ratio has a significant direct impact on soil δ15N with a negative correlation. Soil clay, soil N content, and soil pH have a total positive effect on soil δ15N, but the total positive impact of soil pH is very weak because it has a negative indirect impact on soil δ15N by affecting soil clay, soil N and soil C/N ratio. The total influence is, in order, MAP > MAT > soil C/N > soil clay > soil N > soil pH (in absolute value). The above results will provide valuable information about ecosystem N cycle in temperate grassland of northern China.

Elevational shifts in foliar-soil δ15 N in the Hengduan Mountains and different potential mechanisms

The natural abundance of stable nitrogen isotopes (δ¹⁵ N) provides insights into the N dynamics of terrestrial ecosystems, the determination of which is considered an effective approach for gaining a better understanding ecosystem N cycling. However, there is currently little information available regarding the patterns and mechanisms underlying the variation in foliar-soil δ¹⁵ N among mountain ecosystems. In this study, we examined the determinants of foliar-soil δ¹⁵ N in association with N transportation rates along an elevational gradient in the Hengduan Mountains. Despite the relatively high levels of available N produced from high N fixation and mineralization, we detected the lowest levels of foliar δ¹⁵ N at 3500 m a.s.l., reflecting the stronger vegetation N limitation at medium high elevations. The enhanced vegetation N limitation was driven by the combined effects of higher microbial immobilization and inherent plant dynamic (the shifts of δ¹⁵ N in vegetation preference, including vegetation community) with changing climate along the elevational gradient. Unexpectedly, we established that soil δ¹⁵ N was characterized by an undulating rise and uncoupled correlation with foliar δ¹⁵ N with increasing elevation, thereby indicating that litter input might not be a prominent driver of soil δ¹⁵ N. Conversely, soil nitrification and denitrification were found to make a more pronounced contribution to the pattern of soil δ¹⁵ N along the elevational gradient. Collectively, our results serve to highlight the importance of microbial immobilization in soil N dynamics and provide novel insights that will contribute to enhancing our understanding of N cycling as indicated by foliar-soil δ¹⁵ N along elevational gradients.

Warming and redistribution of nitrogen inputs drive an increase in terrestrial nitrous oxide emission factor

Anthropogenic nitrogen inputs cause major negative environmental impacts, including emissions of the important greenhouse gas N2O. Despite their importance, shifts in terrestrial N loss pathways driven by global change are highly uncertain. Here we present a coupled soil-atmosphere isotope model (IsoTONE) to quantify terrestrial N losses and N2O emission factors from 1850-2020. We find that N inputs from atmospheric deposition caused 51% of anthropogenic N2O emissions from soils in 2020. The mean effective global emission factor for N2O was 4.3 ± 0.3% in 2020 (weighted by N inputs), much higher than the surface area-weighted mean (1.1 ± 0.1%). Climate change and spatial redistribution of fertilisation N inputs have driven an increase in global emission factor over the past century, which accounts for 18% of the anthropogenic soil flux in 2020. Predicted increases in fertilisation in emerging economies will accelerate N2O-driven climate warming in coming decades, unless targeted mitigation measures are introduced.

δ15N of Chironomidae: An index of nitrogen sources and processing within watersheds for national aquatic monitoring programs

Nitrogen (N) removal along flowpaths to aquatic ecosystems is an important regulating ecosystem service that can help reduce N pollution in the nation's waterways, but can be challenging to measure at large spatial scales. Measurements that integrate N processing within watersheds would be particularly useful for assessing the magnitude of this vital service. Because most N removal processes cause isotopic fractionation, δ¹⁵N from basal food-chain organisms in aquatic ecosystems can provide information on both N sources and the degree of watershed N processing. As part of EPA's National Aquatic Resource Surveys (NARS), we measured δ¹⁵N of Chironomidae collected from over 2000 lakes, rivers and streams across the continental USA. Using information on N inputs to watersheds and summer total N concentrations ([TN]) in the water column, we assessed where elevated chironomid δ¹⁵N would indicate N removal rather than possible enriched sources of N. Chironomid δ¹⁵N values ranged from -4 to +20‰, and were higher in rivers and streams than in lakes, indicating that N in rivers and streams underwent more processing and cycling that preferentially removes ¹⁴N than N in lakes. Chironomid δ¹⁵N increased with watershed size, N inputs, and water chemical components, and decreased as precipitation increased. In rivers and streams with high watershed N inputs, we found lower [TN] in streams with higher chironomid δ¹⁵N values, suggesting high rates of gaseous N loss such as denitrification. At low watershed N inputs, the pattern reversed; streams with elevated chironomid δ¹⁵N had higher [TN] than streams with lower chironomid δ¹⁵N, possibly indicating unknown sources elevated in δ¹⁵N such as legacy N, or waste from animals or humans. Chironomid δ¹⁵N values can be a valuable tool to assess integrated watershed-level N sources, input rates, and processing for water quality monitoring and assessment at large scales.

Foliar ẟ15N patterns in legumes and non-N fixers across a climate gradient, Hawai'i Island, USA

Recent studies from the Hawaiian Islands showed that pedogenic thresholds demarcate domains in which rock-derived nutrient dynamics remain similar across wide variations in rainfall. These thresholds appear related to certain aspects of N cycling, but the degree to which they correspond to patterns of biological N fixation (BNF)-the dominant input of N into less-managed ecosystems-remains unclear. We measured aboveground plant biomass, foliar nutrient concentrations, and foliar δ¹⁵N along a climate gradient on ~ 150,000-year-old basaltic substrate to characterize foliar N sources and spatially relate them to soil nutrients. Patterns in legume δ¹⁵N correspond to known pedogenic thresholds along the rainfall gradient, with low δ¹⁵N values (~ 0 to - 2‰) occurring in the dry, biologically inactive domain and the wet, highly weathered domain. Elevated δ¹⁵N in the middle, fertile domain suggests a greater reliance of legumes on soil N where it has accumulated over time. Non-legume face N deficiencies throughout most of the gradient while legumes maintain low C:N ratios via symbiotic BNF. However, legume abundance declines outside the fertile domain, limiting ecosystem N inputs. Breakpoints in legume δ¹⁵N data suggest that P (and potentially other nutrients) limits BNF and, by extension, legume abundance in wet region. Nutrients may also constrain legume abundance in the dry domain, but pedogenic effects could not be isolated from climatic constraints at the dry sites. We conclude that pedogenic thresholds defined by climate can be informative of foliar δ¹⁵N patterns in cases where legumes are not directly constrained by climate, land use, or other external factors.

Common mechanisms explain nitrogen-dependent growth of Arctic shrubs over three decades despite heterogeneous trends and declines in soil nitrogen availability

Heterogeneity has been observed in the responses of Arctic shrubs to climate variability over recent decades, which may reflect landscape-scale variability in belowground resources. At a northern fringe of tall shrub expansion (Yuribei, Yamal Peninsula, Russia), we sought to determine the mechanisms relating nitrogen (N) limitation to shrub growth over decadal time. We analysed the ratio of ¹⁵ N to ¹⁴ N isotopes in wood rings of 10 Salix lanata individuals (399 measurements) to reconstruct annual point-based bioavailable N between 1980 and 2013. We applied a model-fitting/model-selection approach with a suite of competing ecological models to assess the most-likely mechanisms that explain each shrub's individual time-series. Shrub δ¹⁵ N time-series indicated declining (seven shrubs), increasing (two shrubs) and no trend (one shrub) in N availability. The most appropriate model for all shrubs included N-dependent growth of linear rather than saturating form. Inclusion of plant-soil feedbacks better explained ring width and δ¹⁵ N for eight of 10 individuals. Although N trajectories were individualistic, common mechanisms of varying strength confirmed the N-dependency of shrub growth. The linear mechanism may reflect intense scavenging of scarce N; the importance of plant-soil feedbacks suggests that shrubs subvert the microbial bottleneck by actively controlling their environment.

Long-term nitrogen isotope dynamics in Encelia farinosa reflect plant demographics and climate

While plant δ¹⁵ N values have been applied to understand nitrogen (N) dynamics, uncertainties regarding intraspecific and temporal variability currently limit their application. We used a 28 yr record of δ¹⁵ N values from two Mojave Desert populations of Encelia farinosa to clarify sources of population-level variability. We leveraged > 3500 foliar δ¹⁵ N observations collected alongside structural, physiological, and climatic data to identify plant and environmental contributors to δ¹⁵ N values. Additional sampling of soils, roots, stems, and leaves enabled assessment of the distribution of soil N content and δ¹⁵ N, intra-plant fractionations, and relationships between soil and plant δ¹⁵ N values. We observed extensive within-population variability in foliar δ¹⁵ N values and found plant age and foliar %N to be the strongest predictors of individual δ¹⁵ N values. There were consistent differences between root, stem, and leaf δ¹⁵ N values (spanning c. 3‰), but plant and bulk soil δ¹⁵ N values were unrelated. Plant-level variables played a strong role in influencing foliar δ¹⁵ N values, and interannual relationships between climate and δ¹⁵ N values were counter to previously recognized spatial patterns. This long-term record provides insights regarding the interpretation of δ¹⁵ N values that were not available from previous large-scale syntheses, broadly enabling more effective application of foliar δ¹⁵ N values.

Magdalenian and Epimagdalenian chronology and palaeoenvironments at Kůlna Cave, Moravia, Czech Republic

Kůlna Cave is the only site in Moravia, Czech Republic, from which large assemblages of both Magdalenian and Epimagdalenian archaeological materials have been excavated from relatively secure stratified deposits. The site therefore offers the unrivalled opportunity to explore the relationship between these two archaeological phases. In this study, we undertake radiocarbon, stable isotope (carbon, nitrogen and sulphur), and ZooMS analysis of the archaeological faunal assemblage to explore the chronological and environmental context of the Magdalenian and Epimagdalenian deposits. Our results show that the Magdalenian and Epimagdalenian deposits can be understood as discrete units from one another, dating to the Late Glacial between c. 15,630 cal. BP and 14,610 cal. BP, and c. 14,140 cal. BP and 12,680 cal. BP, respectively. Stable isotope results (δ¹³C, δ¹⁵N, δ³⁴S) indicate that Magdalenian and Epimagdalenian activity at Kůlna Cave occurred in very different environmental settings. Magdalenian occupation took place within a nutrient-poor landscape that was experiencing rapid changes to environmental moisture, potentially linked to permafrost thaw. In contrast, Epimagdalenian occupation occurred in a relatively stable, temperate environment composed of a mosaic of woodland and grassland habitats. The potential chronological gap between the two phases, and their associations with very different environmental conditions, calls into question whether the Epimagdalenian should be seen as a local, gradual development of the Magdalenian. It also raises the question of whether the gap in occupation at Kůlna Cave could represent a change in settlement dynamics and/or behavioural adaptations to changing environmental conditions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12520-020-01254-4.

Plant response to N availability in permafrost-affected alpine wetlands in arid and semi-arid climate zones

Nutrient cycling in alpine permafrost-affected wetlands remains insufficiently studied, as it is influenced by a complex network of interrelated climatic and environmental factors, at both regional and local scale. Therefore, we applied mathematical models to examine relationship between environmental factors and plant functional traits reflecting N availability in wetland communities developed under locally variable conditions in a geographic and climatic gradient of high-altitude habitats. Moreover, we assessed impact of local differences in soil chemistry on plant fractionation of N isotopes as a response to N availability. Based on environmental data and chemistry of biomass from 192 study sites from the Pamir Mountains (Tajikistan) and Khangai and Khentei Mountains (Mongolia), a matrix of rank correlations was prepared for regional and local factors and community level plant functional traits. For the traits that were highly correlated either with regional or with local drivers (that is plant N:P ratio and plant δ¹⁵N), linear models were built, with a limited set of predictors selected according to the Risk Inflation Criterion and the SOS algorithm. The models were fitted for each of the studied regions. Presented regional models indicated significant influence of soil NH₄⁺ and/or PO₄^3- content on plant N:P ratio, which showed increase with altitude and lowering precipitation. Thus, its values clearly distinguished between the Pamir Mountains (high N:P) and the Mongolian ranges (low N:P). Models for plant δ¹⁵N showed its strong positive correlations with soil δ¹⁵N and soil salinity. Average values of plant δ¹⁵N were comparable for both study areas. The studied plant functional traits showed different response to regional and local drivers. Plant N:P ratio was controlled by regional drivers via their influence on soil NH₄⁺ content. Contrastingly, plant δ¹⁵N was significantly affected by local factors, namely soil δ¹⁵N and soil salinity expressed as Na:EC.

Climatic, Edaphic and Biotic Controls over Soil δ13C and δ15N in Temperate Grasslands

Soils δ13C and δ15N are now regarded as useful indicators of nitrogen (N) status and dynamics of soil organic carbon (SOC). Numerous studies have explored the effects of various factors on soils δ13C and δ15N in terrestrial ecosystems on different scales, but it remains unclear how co-varying climatic, edaphic and biotic factors independently contribute to the variation in soil δ13C and δ15N in temperate grasslands on a large scale. To answer the above question, a large-scale soil collection was carried out along a vegetation transect across the temperate grasslands of Inner Mongolia. We found that mean annual precipitation (MAP) and mean annual temperature (MAT) do not correlate with soil δ15N along the transect, while soil δ13C linearly decreased with MAP and MAT. Soil δ15N logarithmically increased with concentrations of SOC, total N and total P. By comparison, soil δ13C linearly decreased with SOC, total N and total P. Soil δ15N logarithmically increased with microbial biomass C and microbial biomass N, while soil δ13C linearly decreased with microbial biomass C and microbial biomass N. Plant belowground biomass linearly increased with soil δ15N but decreased with soil δ13C. Soil δ15N decreased with soil δ13C along the transect. Multiple linear regressions showed that biotic and edaphic factors such as microbial biomass C and total N exert more effect on soil δ15N, whereas climatic and edaphic factors such as MAT and total P have more impact on soil δ13C. These findings show that soil C and N cycles in temperate grasslands are, to some extent, decoupled and dominantly controlled by different factors. Further investigations should focus on those ecological processes leading to decoupling of C and N cycles in temperate grassland soils.

Radiocarbon chronology and environmental context of Last Glacial Maximum human occupation in Switzerland

Central Europe during the Last Glacial Maximum (LGM) was dominated by polar desert and steppe-tundra biomes. Despite this, a human presence during this time period is evident at several locations across the region, including in Switzerland, less than 50 km from the Alpine ice sheet margin. It has been hypothesised that such human activity may have been restricted to brief periods of climatic warming within the LGM, but chronological information from many of these sites are currently too poorly resolved to corroborate this. Here we present a revised chronology of LGM human occupation in Switzerland. AMS radiocarbon dating of cut-marked reindeer (Rangifer tarandus) bones from the sites of Kastelhöhle-Nord and Y-Höhle indicates human occupation of Switzerland was most likely restricted to between 23,400 and 22,800 cal. BP. This timeframe corresponds to Greenland Interstadial 2, a brief warming phase, supporting the hypothesis that human presence was facilitated by favourable climatic episodes. Carbon, nitrogen and sulphur stable isotope analysis of the fauna provides palaeoenvironmental information for this time period. These findings contribute to our understanding of human activity in ice-marginal environments and have implications for understanding cultural connections across central Europe during the LGM.

Exploring cycad foliage as an archive of the isotopic composition of atmospheric nitrogen

Molecular nitrogen (N₂ ) constitutes the majority of Earth's modern atmosphere, contributing ~0.79 bar of partial pressure (pN₂ ). However, fluctuations in pN₂ may have occurred on 10⁷ -10⁹  year timescales in Earth's past, perhaps altering the isotopic composition of atmospheric nitrogen. Here, we explore an archive that may record the isotopic composition of atmospheric N₂ in deep time: the foliage of cycads. Cycads are ancient gymnosperms that host symbiotic N₂ -fixing cyanobacteria in modified root structures known as coralloid roots. All extant species of cycads are known to host symbionts, suggesting that this N₂ -fixing capacity is perhaps ancestral, reaching back to the early history of cycads in the late Paleozoic. Therefore, if the process of microbial N₂ fixation records the δ¹⁵ N value of atmospheric N₂ in cycad foliage, the fossil record of cycads may provide an archive of atmospheric δ¹⁵ N values. To explore this potential proxy, we conducted a survey of wild cycads growing in a range of modern environments to determine whether cycad foliage reliably records the isotopic composition of atmospheric N₂ . We find that neither biological nor environmental factors significantly influence the δ¹⁵ N values of cycad foliage, suggesting that they provide a reasonably robust record of the δ¹⁵ N of atmospheric N₂ . Application of this proxy to the record of carbonaceous cycad fossils may not only help to constrain changes in atmospheric nitrogen isotope ratios since the late Paleozoic, but also could shed light on the antiquity of the N₂ -fixing symbiosis between cycads and cyanobacteria.

Microbial Community, Newly Sequestered Soil Organic Carbon, and δ15N Variations Driven by Tree Roots

Rhizosphere microbes in forests are key elements of the carbon sequestration of terrestrial ecosystems. To date, little is known about how the diversity and species interactions of the active rhizomicrobial community change during soil carbon sequestration and what interactions drive these changes. In this study, we used a combination of DNA and stable isotope method to explore correlations between the composition of microbial communities, N transformation, and the sequestration de novo of carbon in soils around Pinus tabuliformis and Quercus variabilis roots in North China. Rhizosphere soils from degraded lands, primary stage land (tree roots had colonized in degraded soil for 1 year), and nature forest were sampled for analyses. The results showed that microbial communities and newly sequestered soil organic carbon (SOC) contents changed with different tree species, environments, and successive stages. The fungal unweighted and weighted UniFrac distances could better show the different microbial species structures and differences in successive stages. Newly sequestered SOC was positively correlated with the bacterial order Rhizobiales (in P. tabuliformis forests), the fungal order Russulales (in Q. variabilis forests), and δ15N. Consequently, the bacterial order Rhizobiales acted as an important taxa for P. tabuliformis root-driven carbon sequestration, and the fungal order Russulales acted as an important taxa for Q. variabilis root-driven carbon sequestration. The two plant species allocated root exudates to different portion of their root systems, which in turn altered microbial community composition and function. The δ15N of soil organic matter could be an important indicator to estimate root-driven carbon sequestration.

Leaf and Soil δ15N Patterns Along Elevational Gradients at Both Treelines and Shrublines in Three Different Climate Zones

The natural abundance of stable nitrogen (N) isotope (δ15N) in plants and soils can reflect N cycling processes in ecosystems. However, we still do not fully understand patterns of plant and soil δ15N at alpine treelines and shrublines in different climate zones. We measured δ15N and N concentration in leaves of trees and shrubs and also in soils along elevational gradients from lower altitudes to the upper limits of treelines and shrublines in subtropical, dry- and wet-temperate regions in China. The patterns of leaf δ15N in trees and shrubs in response to altitude changes were consistent, with lower values occurring at higher altitude in all three climate zones, but such patterns did not exist for leaf Δδ15N and soil δ15N. Average δ15N values of leaves (−1.2‰) and soils (5.6‰) in the subtropical region were significantly higher than those in the two temperate regions (−3.4‰ and 3.2‰, respectively). Significant higher δ15N values in subtro4pical forest compared with temperate forests prove that N cycles are more open in warm regions. The different responses of leaf and soil δ15N to altitude indicate complex mechanisms of soil biogeochemical process and N sources uptake with environmental variations.

In Situ Quantification of Biological N2 Production Using Naturally Occurring 15N15N

We describe an approach for determining biological N₂ production in soils based on the proportions of naturally occurring ¹⁵N¹⁵N in N₂. Laboratory incubation experiments reveal that biological N₂ production, whether by denitrification or anaerobic ammonia oxidation, yields proportions of ¹⁵N¹⁵N in N₂ that are within 1‰ of that predicted for a random distribution of ¹⁵N and ¹⁴N atoms. This relatively invariant isotopic signature contrasts with that of the atmosphere, which has ¹⁵N¹⁵N proportions in excess of the random distribution by 19.1 ± 0.1‰. Depth profiles of gases in agricultural soils from the Kellogg Biological Station Long-Term Ecological Research site show biological N₂ accumulation that accounts for up to 1.6% of the soil N₂. One-dimensional reaction-diffusion modeling of these soil profiles suggests that subsurface N₂ pulses leading to surface emission rates as low as 0.3 mmol N₂ m^-2 d^-1 can be detected with current analytical precision, decoupled from N₂O production.

Soil organic carbon stability in forests: Distinct effects of tree species identity and traits

Rising atmospheric CO₂ concentrations have increased interest in the potential for forest ecosystems and soils to act as carbon (C) sinks. While soil organic C contents often vary with tree species identity, little is known about if, and how, tree species influence the stability of C in soil. Using a 40 year old common garden experiment with replicated plots of eleven temperate tree species, we investigated relationships between soil organic matter (SOM) stability in mineral soils and 17 ecological factors (including tree tissue chemistry, magnitude of organic matter inputs to the soil and their turnover, microbial community descriptors, and soil physicochemical properties). We measured five SOM stability indices, including heterotrophic respiration, C in aggregate occluded particulate organic matter (POM) and mineral associated SOM, and bulk SOM δ¹⁵ N and ∆¹⁴ C. The stability of SOM varied substantially among tree species, and this variability was independent of the amount of organic C in soils. Thus, when considering forest soils as C sinks, the stability of C stocks must be considered in addition to their size. Further, our results suggest tree species regulate soil C stability via the composition of their tissues, especially roots. Stability of SOM appeared to be greater (as indicated by higher δ¹⁵ N and reduced respiration) beneath species with higher concentrations of nitrogen and lower amounts of acid insoluble compounds in their roots, while SOM stability appeared to be lower (as indicated by higher respiration and lower proportions of C in aggregate occluded POM) beneath species with higher tissue calcium contents. The proportion of C in mineral associated SOM and bulk soil ∆¹⁴ C, though, were negligibly dependent on tree species traits, likely reflecting an insensitivity of some SOM pools to decadal scale shifts in ecological factors. Strategies aiming to increase soil C stocks may thus focus on particulate C pools, which can more easily be manipulated and are most sensitive to climate change.

Controls of nitrogen cycling evaluated along a well-characterized climate gradient

The supply of nitrogen (N) constrains primary productivity in many ecosystems, raising the question "what controls the availability and cycling of N"? As a step toward answering this question, we evaluated N cycling processes and aspects of their regulation on a climate gradient on Kohala Volcano, Hawaii, USA. The gradient extends from sites receiving <300 mm/yr of rain to those receiving >3,000 mm/yr, and the pedology and dynamics of rock-derived nutrients in soils on the gradient are well understood. In particular, there is a soil process domain at intermediate rainfall within which ongoing weathering and biological uplift have enriched total and available pools of rock-derived nutrients substantially; sites at higher rainfall than this domain are acid and infertile as a consequence of depletion of rock-derived nutrients, while sites at lower rainfall are unproductive and subject to wind erosion. We found elevated rates of potential net N mineralization in the domain where rock-derived nutrients are enriched. Higher-rainfall sites have low rates of potential net N mineralization and high rates of microbial N immobilization, despite relatively high rates of gross N mineralization. Lower-rainfall sites have moderately low potential net N mineralization, relatively low rates of gross N mineralization, and rates of microbial N immobilization sufficient to sequester almost all the mineral N produced. Bulk soil δ¹⁵ N also varied along the gradient, from +4‰ at high rainfall sites to +14‰ at low rainfall sites, indicating differences in the sources and dynamics of soil N. Our analysis shows that there is a strong association between N cycling and soil process domains that are defined using soil characteristics independent of N along this gradient, and that short-term controls of N cycling can be understood in terms of the supply of and demand for N.

Leaf stable isotopes suggest shared ancestry is an important driver of functional diversity

Plant physiological strategies of carbon (C) and nitrogen (N) uptake and metabolism are often regarded as outcomes of environmental selection. This is likely true, but the role of evolutionary history may also be important in shaping patterns of functional diversity. Here, we used leaf C and N stable isotope ratios (δ¹³C, δ¹⁵N) as integrators of physiological processes to assess the relative roles of phylogenetic history and environment in a diverse group of Ericaceae species native to North America. We found strong phylogenetic signal in both leaf δ¹³C and δ¹⁵N, suggesting that close relatives have similar physiological strategies. The signal of phylogeny was generally stronger than that of the local environment. However, within some specialized environments (e.g., wetlands, sandy soils), we found environmental effects and/or niche conservatism. Phylogenetic signal in δ¹³C appears to be most closely related to the constraints on metabolic demand and supply of C, and δ¹⁵N appears to be most strongly related to mycorrhizal associations within the family.

Namib Desert primary productivity is driven by cryptic microbial community N-fixation

Carbon exchange in drylands is typically low, but during significant rainfall events (wet anomalies) drylands act as a C sink. During these anomalies the limitation on C uptake switches from water to nitrogen. In the Namib Desert of southern Africa, the N inventory in soil organic matter available for mineralisation is insufficient to support the observed increase in primary productivity. The C4 grasses that flourish after rainfall events are not capable of N fixation, and so there is no clear mechanism for adequate N fixation in dryland ecosystems to support rapid C uptake. Here we demonstrate that N fixation by photoautotrophic hypolithic communities forms the basis for the N budget for plant productivity events in the Namib Desert. Stable N isotope (δ¹⁵N) values of Namib Desert hypolithic biomass, and surface and subsurface soils were measured over 3 years across dune and gravel plain biotopes. Hypoliths showed significantly higher biomass and lower δ¹⁵N values than soil organic matter. The δ¹⁵N values of hypoliths approach the theoretical values for nitrogen fixation. Our results are strongly indicative that hypolithic communities are the foundation of productivity after rain events in the Namib Desert and are likely to play similar roles in other arid environments.

Plant controls on Late Quaternary whole ecosystem structure and function

Plants and animals influence biomass production and nutrient cycling in terrestrial ecosystems; however, their relative importance remains unclear. We assessed the extent to which mega-herbivore species controlled plant community composition and nutrient cycling, relative to other factors during and after the Late Quaternary extinction event in Britain and Ireland, when two-thirds of the region's mega-herbivore species went extinct. Warmer temperatures, plant-soil and plant-plant interactions, and reduced burning contributed to the expansion of woody plants and declining nitrogen availability in our five study ecosystems. Shrub biomass was consistently one of the strongest predictors of ecosystem change, equalling or exceeding the effects of other biotic and abiotic factors. In contrast, there was relatively little evidence for mega-herbivore control on plant community composition and nitrogen availability. The ability of plants to determine the fate of terrestrial ecosystems during periods of global environmental change may therefore be greater than previously thought.

Soil parent material-A major driver of plant nutrient limitations in terrestrial ecosystems

Because the capability of terrestrial ecosystems to fix carbon is constrained by nutrient availability, understanding how nutrients limit plant growth is a key contemporary question. However, what drives nutrient limitations at global scale remains to be clarified. Using global data on plant growth, plant nutritive status, and soil fertility, we investigated to which extent soil parent materials explain nutrient limitations. We found that N limitation was not linked to soil parent materials, but was best explained by climate: ecosystems under harsh (i.e., cold and or dry) climates were more N-limited than ecosystems under more favourable climates. Contrary to N limitation, P limitation was not driven by climate, but by soil parent materials. The influence of soil parent materials was the result of the tight link between actual P pools of soils and physical-chemical properties (acidity, P richness) of soil parent materials. Some other ground-related factors (i.e., soil weathering stage, landform) had a noticeable influence on P limitation, but their role appeared to be relatively smaller than that of geology. The relative importance of N limitation versus P limitation was explained by a combination of climate and soil parent material: at global scale, N limitation became prominent with increasing climatic constraints, but this global trend was modulated at lower scales by the effect of parent materials on P limitation, particularly under climates favourable to biological activity. As compared with soil parent materials, atmospheric deposition had only a weak influence on the global distribution of actual nutrient limitation. Our work advances our understanding of the distribution of nutrient limitation at global scale. In particular, it stresses the need to take soil parent materials into account when investigating plant growth response to environment changes.

Spatial variations in larch needle and soil δ15N at a forest-grassland boundary in northern Mongolia

The spatial patterns of plant and soil δ¹⁵N and associated processes in the N cycle were investigated at a forest-grassland boundary in northern Mongolia. Needles of Larix sibirica Ledeb. and soils collected from two study areas were analysed to calculate the differences in δ¹⁵N between needle and soil (Δδ¹⁵N). Δδ¹⁵N showed a clear variation, ranging from -8 ‰ in the forest to -2 ‰ in the grassland boundary, and corresponded to the accumulation of organic layer. In the forest, the separation of available N produced in the soil with ¹⁵N-depleted N uptake by larch and ¹⁵N-enriched N immobilization by microorganisms was proposed to cause large Δδ¹⁵N, whereas in the grassland boundary, small Δδ¹⁵N was explained by the transport of the most available N into larch. The divergence of available N between larch and microorganisms in the soil, and the accumulation of diverged N in the organic layer control the variation in Δδ¹⁵N.

Land-use and hydroperiod affect kettle hole sediment carbon and nitrogen biogeochemistry

Kettle holes are glaciofluvially created depressional wetlands that collect organic matter (OM) and nutrients from their surrounding catchment. Kettle holes mostly undergo pronounced wet-dry cycles. Fluctuations in water table, land-use, and management can affect sediment biogeochemical transformations and perhaps threaten the carbon stocks of these unique ecosystems. We investigated sediment and water of 51 kettle holes in NE Germany that differ in hydroperiod (i.e. the duration of the wet period of a kettle hole) and land-use. Our objectives were 1) to test if hydroperiod and land management were imprinted on the isotopic values (δ¹³C, δ¹⁵N) and C:N ratios of the sediment OM, and 2) to characterize water loss dynamics and kettle hole-groundwater connectivity by measuring the stable δ¹⁸O and δD isotope values of kettle hole water over several years. We found the uppermost sediment layer reflected recent OM inputs and short-term processes in the catchment, including land-use and management effects. Deeper sediments recorded the degree to which OM is processed within the kettle hole related to the hydroperiod. We see clear indications for the effects of wet-dry cycles for all kettle holes, which can lead to the encroachment of terrestrial plants. We found that the magnitude of evaporation depended on the year, season, and land-use type, that kettle holes are temporarily coupled to shallow ground water, and, as such, kettle holes are described best as partially-closed to open systems.

Plant δ15 N reflects the high landscape-scale heterogeneity of soil fertility and vegetation productivity in a Mediterranean semiarid ecosystem

We investigated the magnitude and drivers of spatial variability in soil and plant δ¹⁵ N across the landscape in a topographically complex semiarid ecosystem. We hypothesized that large spatial heterogeneity in water availability, soil fertility and vegetation cover would be positively linked to high local-scale variability in δ¹⁵ N. We measured foliar δ¹⁵ N in three dominant plant species representing contrasting plant functional types (tree, shrub, grass) and mycorrhizal association types (ectomycorrhizal or arbuscular mycorrhizal). This allowed us to investigate whether δ¹⁵ N responds to landscape-scale environmental heterogeneity in a consistent way across species. Leaf δ¹⁵ N varied greatly within species across the landscape and was strongly spatially correlated among co-occurring individuals of the three species. Plant δ¹⁵ N correlated tightly with soil δ¹⁵ N and key measures of soil fertility, water availability and vegetation productivity, including soil nitrogen (N), organic carbon (C), plant-available phosphorus (P), water-holding capacity, topographic moisture indices and normalized difference vegetation index. Multiple regression models accounted for 62-83% of within-species variation in δ¹⁵ N across the landscape. The tight spatial coupling and interdependence of the water, N and C cycles in drylands may allow the use of leaf δ¹⁵ N as an integrative measure of variations in moisture availability, biogeochemical activity, soil fertility and vegetation productivity (or 'site quality') across the landscape.

Nitrogen cycling in an extreme hyperarid environment inferred from δ(15)N analyses of plants, soils and herbivore diet

Climate controls on the nitrogen cycle are suggested by the negative correlation between precipitation and δ¹⁵N values across different ecosystems. For arid ecosystems this is unclear, as water limitation among other factors can confound this relationship. We measured herbivore feces, foliar and soil δ¹⁵N and δ¹³C values and chemically characterized soils (pH and elemental composition) along an elevational/climatic gradient in the Atacama Desert, northern Chile. Although very positive δ¹⁵N values span the entire gradient, soil δ¹⁵N values show a positive correlation with aridity as expected. In contrast, foliar δ¹⁵N values and herbivore feces show a hump-shaped relationship with elevation, suggesting that plants are using a different N source, possibly of biotic origin. Thus at the extreme limits of plant life, biotic interactions may be just as important as abiotic processes, such as climate in explaining ecosystem δ¹⁵N values.