Half a century of Scots pine forest ecosystem monitoring reveals long-term effects of atmospheric deposition and climate change.

At two forest sites in Germany (Pfaffenwinkel, Pustert) stocked with mature Scots pine (Pinus sylvestris L.), we investigated changes of topsoil chemistry during the recent 40 years by soil inventories conducted on replicated control plots of fertilization experiments, allowing a statistical analysis. Additionally, we monitored the nutritional status of both stands from 1964 until 2019 and quantified stand growth during the monitoring period by repeated stand inventories. Moreover, we monitored climate variables (air temperature and precipitation) and calculated annual climatic water balances from 1991 to 2019. Atmospheric nitrogen (N) and sulfur (S) deposition between 1964 and 2019 was estimated for the period 1969-2019 by combining annual deposition measurements conducted in 1985-1987 and 2004 with long-term deposition records from long-term forest monitoring stations. We investigated interrelations between topsoil chemistry, stand nutrition, stand growth, deposition, and climate trends. At both sites, the onset of the new millennium was a turning point of important biogeochemical processes. Topsoil acidification turned into re-alkalinization, soil organic matter (SOM) accumulation stopped, and likely turned into SOM depletion. In the new millennium, topsoil stocks of S and plant-available phosphorus (P) as well as S and P concentrations in Scots pine foliage decreased substantially; yet, age-referenced stand growth remained at levels far above those expected from yield table data. Tree P and S nutrition as well as climate change (increased temperature and drought stress) have replaced soil acidification as major future challenges for both forests. Understanding of P and S cycling and water fluxes in forest ecosystems, and consideration of these issues in forest management is important for successfully tackling the new challenges. Our study illustrates the importance of long-term forest monitoring to identify slow, but substantial changes of forest biogeochemistry driven by natural and anthropogenic global change.

A Novel Approach for the Quantification of Different Inorganic and Organic Phosphorus Compounds in Environmental Samples by P L2,3-Edge X-ray Absorption Near-Edge Structure (XANES) Spectroscopy.

Phosphorus (P) is an essential element for life on Earth, with an important and oftentimes unaccounted organic biogeochemical component. Current methods for the quantification of different organic P compounds in environmental samples (e.g., soils, sediments) are based on extraction techniques and often associated with incomplete P recovery or sample changes. In this study, we present a protocol for the quantification of different organic and inorganic P species in soils using synchrotron-based X-ray absorption near-edge structure (XANES) spectroscopy at the P L2,3-edge. Its accuracy and precision was evaluated by analyzing 40 standard mixtures composed of seven different inorganic and organic P compounds (with a mean of R2 = 0.85). In addition, we quantified the P species of two soils and two agro-industrial byproducts using P L2,3-edge XANES spectroscopy and the results were compared with those obtained by P K-edge XANES or 31P NMR spectroscopy. Using the P L2,3-edge, we identified different organic P species, including those not identified by the common P K-edge XANES. However, there is a consistent underestimation of organic polyphosphates. Overall, the application of P L2,3-edge XANES provides a higher level of information than by P K-edge XANES, although the ubiquitous use of this novel methodology is still limited to samples with a phosphorus content above 3 mg g-1.

Phosphorus K-edge XANES spectroscopy has probably often underestimated iron oxyhydroxide-bound P in soils.

Phosphorus (P) K-edge X-ray absorption near-edge structure (XANES) spectra of orthophosphate (oPO4) bound to soil FeIII minerals (e.g. ferrihydrite, goethite) show a pre-edge signal at 2148-2152 eV. It is unknown whether organic P bound to FeIII oxyhydroxides also show this feature. Otherwise, Fe-bound soil P may be underestimated by P K-edge XANES spectroscopy, because a large portion of Fe oxyhydroxide-bound P in soils is organic P. K-edge XANES spectra were obtained for different organic P compounds present in soils [inositol hexaphosphate (IHP), glucose-6-phosphate (G6P), adenosine triphosphate (ATP)] after sorption to ferrihydrite or goethite and compared with spectra of oPO4 adsorbed to these minerals. P sorption to ferrihydrite increased in the sequence IHP ≪ G6P < oPO4 < ATP. P sorption to goethite increased in the sequence G6P < oPO4 ≪ ATP = IHP. Pre-edge signals in P K-edge XANES spectra of organic P adsorbed to Fe oxyhydroxides were markedly smaller compared with those of oPO4 adsorbed to these minerals and absent for FeIII oxyhydroxide-bound ATP as well as goethite-bound IHP. Linear combination fitting (LCF) performed on spectra of IHP, G6P or ATP adsorbed to ferrihydrite or goethite, using only spectra of FeIII oxyhydroxide-bound oPO4 as reference compounds for Fe-bound P, erroneously assigned >93% (ferrihydrite) or >41% (goethite) of Fe-bound P to non-Fe-bound P species. Inclusion of FeIII oxyhydroxide-bound IHP as reference compounds markedly increased the recovery of oxyhydroxide-bound organic P. Thus, Fe-bound soil P has probably often been underestimated by LCF in soil XANES studies where IHP adsorbed to ferrihydrite and to goethite were not included as reference compounds.

Reference spectra of important adsorbed organic and inorganic phosphate binding forms for soil P speciation using synchrotron-based K-edge XANES spectroscopy.

Direct speciation of soil phosphorus (P) by linear combination fitting (LCF) of P K-edge XANES spectra requires a standard set of spectra representing all major P species supposed to be present in the investigated soil. Here, available spectra of free- and cation-bound inositol hexakisphosphate (IHP), representing organic P, and of Fe, Al and Ca phosphate minerals are supplemented with spectra of adsorbed P binding forms. First, various soil constituents assumed to be potentially relevant for P sorption were compared with respect to their retention efficiency for orthophosphate and IHP at P levels typical for soils. Then, P K-edge XANES spectra for orthophosphate and IHP retained by the most relevant constituents were acquired. The spectra were compared with each other as well as with spectra of Ca, Al or Fe orthophosphate and IHP precipitates. Orthophosphate and IHP were retained particularly efficiently by ferrihydrite, boehmite, Al-saturated montmorillonite and Al-saturated soil organic matter (SOM), but far less efficiently by hematite, Ca-saturated montmorillonite and Ca-saturated SOM. P retention by dolomite was negligible. Calcite retained a large portion of the applied IHP, but no orthophosphate. The respective P K-edge XANES spectra of orthophosphate and IHP adsorbed to ferrihydrite, boehmite, Al-saturated montmorillonite and Al-saturated SOM differ from each other. They also are different from the spectra of amorphous FePO4, amorphous or crystalline AlPO4, Ca phosphates and free IHP. Inclusion of reference spectra of orthophosphate as well as IHP adsorbed to P-retaining soil minerals in addition to spectra of free or cation-bound IHP, AlPO4, FePO4 and Ca phosphate minerals in linear combination fitting exercises results in improved fit quality and a more realistic soil P speciation. A standard set of P K-edge XANES spectra of the most relevant adsorbed P binding forms in soils is presented.

Standard Protocol and Quality Assessment of Soil Phosphorus Speciation by P K-Edge XANES Spectroscopy.

Phosphorus (P) in soils is most often bound as phosphate to one or more of the following four elements or compounds: calcium, aluminum, iron, and soil organic matter. A promising method for direct P speciation in soils is synchrotron-based X-ray absorption near edge structure (XANES) spectroscopy at the K-edge of P. However, the quality of this method is debated controversially, partly because a standard protocol for reproducible spectrum deconvolution is lacking and minor modifications of the applied deconvolution procedure can lead to considerable changes in the P speciation results. On the basis of the observation that appropriate baseline correction and edge-step normalization are crucial for correct linear combination (LC) fitting results, we established a standard protocol for the deconvolution and LC fitting of P K-edge XANES spectra. We evaluated the quality of LC fits obtained according to this standard protocol with 16 defined dilute (2 mg P g(-1)) ternary mixtures of aluminum phosphate, iron phosphate, hydroxyapatite, and phytic acid in a quartz matrix. The LC fitting results were compared with the contribution of the different P compounds to total P in the various mixtures. Compared to using a traditional LC fitting procedure, our standard protocol reduced the fitting error by 6% (absolute). However, P portions smaller than 5% should be confirmed with other methods or excluded from the P speciation results. A publicly available database of P K-edge XANES reference spectra was initiated.

Soil phosphorus status and P nutrition strategies of European beech forests on carbonate compared to silicate parent material.

Sustainable forest management requires understanding of ecosystem phosphorus (P) cycling. Lang et al. (2017) [Biogeochemistry, https://doi.org/10.1007/s10533-017-0375-0] introduced the concept of P-acquiring vs. P-recycling nutrition strategies for European beech (Fagus sylvatica L.) forests on silicate parent material, and demonstrated a change from P-acquiring to P-recycling nutrition from P-rich to P-poor sites. The present study extends this silicate rock-based assessment to forest sites with soils formed from carbonate bedrock. For all sites, it presents a large set of general soil and bedrock chemistry data. It thoroughly describes the soil P status and generates a comprehensive concept on forest ecosystem P nutrition covering the majority of Central European forest soils. For this purpose, an Ecosystem P Nutrition Index (ENI P ) was developed, which enabled the comparison of forest P nutrition strategies at the carbonate sites in our study among each other and also with those of the silicate sites investigated by Lang et al. (2017). The P status of forest soils on carbonate substrates was characterized by low soil P stocks and a large fraction of organic Ca-bound P (probably largely Ca phytate) during early stages of pedogenesis. Soil P stocks, particularly those in the mineral soil and of inorganic P forms, including Al- and Fe-bound P, became more abundant with progressing pedogenesis and accumulation of carbonate rock dissolution residue. Phosphorus-rich impure, silicate-enriched carbonate bedrock promoted the accumulation of dissolution residue and supported larger soil P stocks, mainly bound to Fe and Al minerals. In carbonate-derived soils, only low P amounts were bioavailable during early stages of pedogenesis, and, similar to P-poor silicate sites, P nutrition of beech forests depended on tight (re)cycling of P bound in forest floor soil organic matter (SOM). In contrast to P-poor silicate sites, where the ecosystem P nutrition strategy is direct biotic recycling of SOM-bound organic P, recycling during early stages of pedogenesis on carbonate substrates also involves the dissolution of stable Ca-Porg precipitates formed from phosphate released during SOM decomposition. In contrast to silicate sites, progressing pedogenesis and accumulation of P-enriched carbonate bedrock dissolution residue at the carbonate sites promote again P-acquiring mechanisms for ecosystem P nutrition. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s10533-021-00884-7.

Sulfur speciation in soil by S K-Edge XANES spectroscopy: comparison of spectral deconvolution and linear combination fitting.

Defined, quartz-diluted mixtures of sulfur (S) compounds with different oxidation state (OS) were analyzed by K-edge XANES spectroscopy using linear combination fitting (LCF) and spectrum deconvolution by fitting several Gaussian and arctangent functions (GCF). Additionally, for different soils the S speciation as calculated by both methods was compared with results of a wet-chemical S speciation. For mixtures of FeS, L-cysteine, and Na2SO4, the S speciation was recovered with satisfactory accuracy and precision by both methods at the 2 and 0.2 mg S g(-1) level. For GCF, white-line peaks must be normalized with respect to their OS-specific absorption cross-section. LCF must be conducted with dilute reference compounds to avoid self-absorption effects. For mixtures of FeS, FeS2, S°, and L-cysteine, both procedures showed poor accuracy. For the soils, similar percentages of reduced inorganic S, organic S, and sulfate were calculated by LCF, GCF, and wet chemical S speciation. GCF allows a fair estimation of S species groups with different OS (inorganic reduced S, organic reduced S, organic intermediate S, oxidized S) in soils without standards. If dilute standards of all S compounds assumed to be present in a sample are available, LCF is more objective and allows a more detailed S speciation.

Influence of Physical-Chemical Soil Parameters on Microbiota Composition and Diversity in a Deep Hyperarid Core of the Atacama Desert.

The extreme environmental conditions and lack of water on the soil surface in hyperarid deserts hamper microbial life, allowing only highly specialized microbial communities to the establish colonies and survive. Until now, the microbial communities that inhabit or have inhabited soils of hyperarid environments at greater depths have been poorly studied. We analyzed for the first time the variation in microbial communities down to a depth of 3.4 m in one of the driest places of the world, the hyperarid Yungay region in the Atacama Desert, and we related it to changes in soil physico-chemical characteristics. We found that the moisture content changed from 2 to 11% with depth and enabled the differentiation of three depth intervals: (i) surface zone A (0-60 cm), (ii) intermediate zone B (60-220 cm), and (iii) deep zone C (220-340 cm). Each zone showed further specific physicochemical and mineralogical features. Likewise, some bacterial phyla were unique in each zone, i.e., members of the taxa Deinococcota, Halobacterota, and Latescibacterota in zone A; Crenarchaeota, Fusobacteriota, and Deltaproteobacterium Sva0485 in zone B; and Fervidibacteria and Campilobacterota in zone C, which indicates taxon-specific preferences in deep soil habitats. Differences in the microbiota between the zones were rather abrupt, which is concomitant with abrupt changes in the physical-chemical parameters. Overall, moisture content, total carbon (TC), pH, and electric conductivity (EC) were most predictive of microbial richness and diversity, while total sulfur (TS) and total phosphorous (TP) contents were additionally predictive of community composition. We also found statistically significant associations between taxa and soil properties, most of which involved moisture and TC contents. Our findings show that under-explored habitats for microbial survival and existence may prevail at greater soil depths near water or within water-bearing layers, a valuable substantiation also for the ongoing search for biosignatures on other planets, such as Mars.

Assessment of sulfur and iron speciation in a soil aggregate by combined S and Fe micro-XANES: microspatial patterns and relationships.

To test whether synchrotron-based spectromicroscopy can be used to identify spatial patterns of sulfur (S) and iron (Fe) speciation as well as relationships between the speciation of S and Fe in soil colloids or aggregates at the micrometre and sub-micrometre level, an anoxically prepared dissected soil aggregate (size approximately 1 mm(3)) was analyzed by micro-XANES at the K-edges of S (2472 eV) and Fe (7112 eV). The experiment included (i) elemental mapping at the S K-edge (S, Si, Al) and the Fe K-edge (Fe, Si), (ii) acquisition of 300 microm x 300 microm images of the region of interest with X-ray energies of 2474 eV (addressing reduced organic and inorganic S), 2483 eV (total S), 7121 eV (divalent Fe) and 7200 eV (total Fe), as well as (iii) acquisition of S and Fe micro-XANES spectra at two different positions, where image analysis suggested the dominance of reduced and oxidized S and Fe, respectively. Image analysis revealed a heterogeneous distribution of total Si, S and Fe as well as of different S and Fe species in the aggregate. Microregions which were either enriched in reduced or in oxidized S and Fe could be identified. A microregion with a large contribution of oxidized S (sulfate, sulfonate) to total S contained exclusively Fe(III) oxyhydroxides (probably ferrihydrite) as S-bearing phase, whereas another microregion with a large contribution of reduced organic S (thiol, organic disulfide) to total S contained a small amount of Fe(II)-bearing silicate in addition to the dominating Fe(III) oxyhydroxides. Our results show that combined S and Fe micro-XANES is a powerful tool for studying microscale spatial patterns of S and Fe speciation as well as microscale relationships between the speciation of S and Fe in soil aggregates.

Comparison of soil organic carbon speciation using C NEXAFS and CPMAS 13C NMR spectroscopy.

We compared synchrotron-based C near-edge X-ray absorption fine structure (NEXAFS) and CPMAS 13C nuclear magnetic resonance (NMR) spectroscopy with respect to their precision and accuracy to quantify different organic carbon (OC) species in defined mixtures of soil organic matter source compounds. We also used both methods to quantify different OC species in organic surface horizons of a Histic Leptosol as well as in mineral topsoil and subsoil horizons of two soils with different parent material, stage of pedogenesis, and OC content (Cambisol: 15-30 OC mgg-1, Podzol: 0.9-7 OC mgg-1). CPMAS 13C NMR spectroscopy was more accurate and precise (mean recovery of different C functional groups 96-103%) than C NEXAFS spectroscopy (mean recovery 92-113%). For organic surface and topsoil samples, NMR spectroscopy consistently yielded larger O-alkyl C percentages and smaller alkyl C percentages than C NEXAFS spectroscopy. For the Cambisol subsoil samples both methods performed well and showed similar C speciation results. NEXAFS spectroscopy yielded excellent spectra with a high signal-to-noise ratio also for OC-poor Podzol subsoil samples, whereas this was not the case for CPMAS 13C NMR spectroscopy even after sample treatment with HF. Our results confirm the analytical power of CPMAS 13C NMR spectroscopy for a reliable quantitative OC speciation in soils with >10mgOCg-1. Moreover, they highlight the potential of synchrotron-based C NEXAFS spectroscopy as fast, non-invasive method to semi-quantify different C functional groups in soils with low C content (0.9-10mgg-1).

Micro-scale heterogeneity of soil phosphorus depends on soil substrate and depth.

Soils comprise various heterogeneously distributed pools of lithogenic, free organic, occluded, adsorbed, and precipitated phosphorus (P) forms, which differ depending on soil forming factors. Small-scale heterogeneity of element distributions recently has received increased attention in soil science due to its influence on soil functions and soil fertility. We investigated the micro-scale distribution of total P and different specific P binding forms in aggregates taken from a high-P clay-rich soil and a low-P sandy soil by combining advanced spectrometric and spectroscopic techniques to introduce new insights on P accessibility and availability in soils. Here we show that soil substrate and soil depth determine micro-scale P heterogeneity in soil aggregates. In P-rich areas of all investigated soil aggregates, P was predominantly co-located with aluminium and iron oxides and hydroxides, which are known to strongly adsorb P. Clay minerals were co-located with P only to a lesser extent. In the low-P topsoil aggregate, the majority of the P was bound organically. Aluminium and iron phosphate predominated in the quartz-rich low-P subsoil aggregate. Sorbed and mineral P phases determined P speciation in the high-P top- and subsoil, and apatite was only detected in the high-P subsoil aggregate. Our results indicate that micro-scale spatial and chemical heterogeneity of P influences P accessibility and bioavailability.

Cumulative impact of 40 years of industrial sulfur emissions on a forest soil in west-central Alberta (Canada).

The impact of 40 years of sulfur (S) emissions from a sour gas processing plant in Alberta (Canada) on soil development, soil S pools, soil acidification, and stand nutrition at a pine (Pinus contorta x Pinus banksiana) ecosystem was assessed by comparing ecologically analogous areas subjected to different S deposition levels. Sulfur isotope ratios showed that most deposited S was derived from the sour gas processing plant. The soil subjected to the highest S deposition contained 25.9 kmol S ha(-1) (uppermost 60 cm) compared to 12.5 kmol S ha(-1) or less at the analogues receiving low S deposition. The increase in soil S pools was caused by accumulation of organic S in the forest floor and accumulation of inorganic sulfate in the mineral soil. High S inputs resulted in topsoil acidification, depletion of exchangeable soil Ca2+ and Mg2+ pools by 50%, podzolization, and deterioration of N nutrition of the pine trees.

Parameters influencing sulfur speciation in environmental samples using sulfur k-edge x-ray absorption near-edge structure.

This paper aims to enhance the credibility of applying the sulfur K-edge XANES spectroscopy as an innovative "fingerprint" for characterizing environmental samples. The sensitivities of sulfur K-edge XANES spectra of ten sulfur compound standards detected by two different detectors, namely, Lytle detector (LyD) and Germanium detector (GeD), were studied and compared. Further investigation on "self-absorption" effect revealed that the maximum sensitivities of sulfur K-edge XANES spectra were achieved when diluting sulfur compound standards with boron nitride (BN) at the mixing ratio of 0.1%. The "particle-size" effect on sulfur K-edge XANES spectrum sensitivities was examined by comparing signal-to-noise ratios of total suspended particles (TSP) and particulate matter of less than 10 millionths of a meter (PM(10)) collected at three major cities of Thailand. The analytical results have demonstrated that the signal-to-noise ratios of sulfur K-edge XANES spectra were positively correlated with sulfate content in aerosols and negatively connected with particle sizes. The combination of hierarchical cluster analysis (HCA) and principal component analysis (PCA) has proved that sulfur K-edge XANES spectrum can be used to characterize German terrestrial soils and Andaman coastal sediments. In addition, this study highlighted the capability of sulfur K-edge XANES spectra as an innovative "fingerprint" to distinguish tsunami backwash deposits (TBD) from typical marine sediments (TMS).