- 1Soil Science, TUM School of Life Sciences Weihenstephan, Technical University of Munich, Munich, Germany
- 2Thünen Institute of Climate-Smart Agriculture, Braunschweig, Germany
- 3Georges Lemaître Centre for Earth and Climate Research-Earth and Life Institute, Université Catholique de Louvain, Ottignies-Louvain-la-Neuve, Belgium
- 4Institute for Advanced Study, Technical University of Munich, Munich, Germany
- 5Bavarian State Research Center for Agriculture, Institute for Organic Farming, Soil and Resource Management, Freising, Germany
Numerous approaches have been developed to isolate fast and slow cycling soil organic carbon (SOC) pools using physical and chemical fractionation. Most of these methods are complex, expensive, and time consuming and unsuited for high-throughput application, such as for regional scale assessments. For simpler and faster fractionation via particle size the key issue is the dispersion of soil. It is unclear how the initial dispersion of soil affects the turnover rates of isolated fractions. We investigated five commonly used dispersion methods using different intensities: shaking in water, shaking in water with glass beads, ultrasonication at 100 and 450 J ml−1 and sodium hexametaphosphate (Na-HMP). We used soils from long-term field experiments that included a change from C3 to C4 vegetation and adjacent control sites using δ13C isotope ratio mass spectrometry. We evaluated the degree of C3/C4 moieties of the fractions, mass and carbon recovery and reproducibility as well as the time expenditures of the dispersions, sieving and drying techniques to develop an efficient and cheap fractionation method. Our results indicate that ultrasonication as well as H2O treatment with and without glass beads resulted in fractions with different turnover. Moreover, isolation performances depended on soil texture. While the isolation of the fractions using water with and without glass beads was equivalent to ultrasonication in soils with low clay contents, these methods had limited potential for soils with high clay contents. Furthermore, treatment with water alone had less reproducible results than other tested methods. The SOC recovery was comparable and satisfactory amongst non-chemical dispersion methods and reached over 95% for each of these methods. The use of Na-HMP was unsuccessful due to high time expenditures and strong SOC leaching. We propose particle size fractionation combined with ultrasonic dispersion as a fast and highly reliable method to quantify slow and fast cycling SOC pools for a wide range of soil types and textures from agricultural sites in central Europe.
Introduction
Soil organic matter (SOM) plays an important role influencing soil functions and the global carbon cycle (1). It determines the nutrient cycle when it is decomposed and can also withdraw carbon (C) from the atmosphere and sequester it in the soil for long periods (2, 3). It is estimated that SOM stores two to three times as much C as the atmosphere (4). Different soil management practices strongly affect SOC stocks and appropriate soil management has the potential to improve nutrient availability and C sequestration (5, 6).
Physical and chemical fractionation methods are used to differentiate between functional SOC fractions containing SOC pools with different turnover times. Physical fractionation methods can be roughly categorized into aggregate and particle fractionations (7–9). Aggregate and particle fractionations isolate fractions by size, density or a mixture of both. Aggregate fractionations are mostly used to explore the influences of management practices on the ecosystem functions of soil aggregates, e.g., their respective pore spaces and microbial habitats (10–12). However, particle fractionations are based on the assumption that the associations between soil particles and their spatial arrangement are important for SOM dynamics because their bioaccessibility is a major driver of SOM decomposition (13). Therefore, particle fractionation methods are mostly used to isolate organic matter (OM) fractions with different turnover times (9, 14, 15). Particle density fractionation is applied to separate SOM that is not firmly associated with soil minerals. These light fractions (LFs) consist of plant residues that are either free in the soil matrix or occluded in particles (13, 15). Particle size fractionation is based on the concept, that decay induces fragmentation and OM formed during decomposition associates with minerals, preventing further decay (9, 13). Particle density fractionations are often combined with particle size fractionations in fractionation schemes to collect additional information on LFs (16–18). While density fractionations allow a more differentiated separation of heavy and light fractions, simple particle size fractionations can be conducted faster and cheaper (8, 9, 13, 14, 19–24).
In this regard, two fractions are most often considered: mineral-associated organic matter (MAOM), which can be assigned to different particle size classes of <20/<50/<63 μm, and particulate organic matter (POM), which is assigned to coarser particle size above these thresholds (25–27). The MAOM fraction is protected from biological decomposition by its binding to the mineral phase and has long mean residence times in the magnitude of decades to millennia (14, 28). However, the POM fraction is less protected from degradation and has short residence times (28–30). Therefore, MAOM fractions are suitable as indicators for SOC sequestration (14, 21), and POM fractions can be used to assess the of direct impacts of agricultural soil management on SOC content on a decadal timescale (31).
Methods for estimating the SOC storage potential of soils have emerged from the quantification of these two fractions, including the OC saturation concept (21, 32–34) and model-based concepts (35–37). For both concepts, measurements of POM and MAOM are needed. The substantial diversity of fractionation methods and their often time-consuming procedures calls for the development of simplified fractionation schemes (27). Poeplau et al. (24) used a large-scale approach to compare the isolation quality, time required and costs of different fractionation methods. Their results confirmed that the simple separation of silt and clay from sand-sized fractions provides strong contrasts in turnover rates. Therefore, simple separation of silt-, sand-, and clay-sized fractions is similarly successful at isolating SOC pools with distinct properties as much more complex methods using density fractionation, and it is faster.
All particle size fractionation protocols involve procedures that disperse soil aggregates, making them different from aggregate size fractionation protocols. However, little is known about the effect of different soil dispersion approaches on the isolation efficiency of distinct soil carbon fractions. Prevalent dispersion approaches are least destructive pre-treatments with water, sometimes with physically disruptive bodies such as glass beads (19, 38–40), the use of ultrasound (41–43), and treatment with chemical dispersants to peptize soil (44–46). In particular, ultrasonication at different intensities has received increasing attention over the past few decades. However, there were indications that the results were not comparable between laboratories. Such laboratory bias could be reduced when parameters are kept constant (47).
In this study, we evaluated a simplified fractionation scheme that allows the quantification of fast- and slow-cycling SOC pools for a wide range of soils considering time requirements and cost efficiency. Since soil dispersion is one of the key steps of particle size fractionation methods and it is still unclear which dispersion is best suited for large-scale data acquisition of POM and MAOM, the final aim of this study was to investigate and quantify the performance of various dispersion approaches for a particle size fractionation method (>20/ <20 μm). Dispersion efficiency was evaluated based on the separation of young and old SOC after C3 and C4 vegetation changes estimated via δ13C natural abundance measurements. The experiment was conducted on soils from five long-term field experiments in central Europe with a range in texture and SOC content. In addition to natural abundance measurements, further indicators with regard of recovery and mass balance, reproducibility and labor costs were used for the entire fractionation process, and subdivided into the fractionation modules dispersion, wet sieving and drying.
Materials and Methods
Experimental Sites and Soil Collection
We used five different soils from long-term field experiments (LTEs) in central Europe: Braunschweig (Bs), Hohenheim (Ho), Jyndevat (Jy), Rotthalmünster (Ro), and Tänikon (Ta) (Table 1). The experimental design included fields that had undergone a shift from C3 to C4 vegetation and adjacent control plots that were exclusively cultivated with C3 crops.
Table 1. Details of the five experimental sites: Mean annual temperature (MAT), mean annual precipitation (MAP), soil organic carbon content of C3 reference soil (SOCC3) and the C4 soil (SOCC4), proportions of sand, silt and clay, soil pH (H2O), years under C4 vegetation, C4 plant species, δ13C signatures of C3 and C4 soil, and the proportion of C4-derived C (fC4) in the bulk soil.
All LTEs were chosen with a focus on sampling a variety of soil types and a gradient in soil texture to create a dataset representative of central European soils. Soil clay contents ranged from 4% in Braunschweig to 29% in Tänikon, and silt contents ranged from 6% in Jyndevad to 75% in Rotthalmünster (Table 1). Mean annual temperatures were typical for temperate regions in central Europe ranging between 7.9°C in Jyndevad to 8.8°C in Braunschweig, whereas the mean annual precipitation was more diverse and ranged from 620 mm in Braunschweig to 1,185 mm in Tänikon. The SOC contents ranged from 14 mg g−1 in Rotthalmünster, 18 mg g−1 in Braunschweig, 20 mg g−1 in Jyndevad, and 19 mg g−1 (C3-site) to 34 mg g−1 (C4-site) in Tänikon. The SOC contents in Hohenheim were considerably higher than at the other sites due to its grassland history with approximately 80 mg g−1 (C3-site) and 58 mg g−1 (C4-site). Four out of five experimental sites were cultivated with the C4-plant Miscanthus x giganteus. The experimental site in Rotthalmünster was cultivated with the C4-plant Zea mays, and there was intermittent interruption of maize cultivation with millet (Sorghum spec.) in the years 2012 and 2014. The soil samples were collected after 17–36 years of C4/C3 plant cultivation (Table 1) as composite samples on ~ 300 m2 grids, which were described in (48), within a depth of 10 cm. Then, the samples were oven dried at 60°C and sieved to 2 mm.
Soil Fractionation
All experiments were conducted using~10 g of dried soil of each sample (Figure 1). We tested and compared the performance of five commonly used dispersion methods at isolating two fractions with distinct SOM turnover times. As we were looking for an approach that can easily be adapted by any laboratory, we tried to minimize complexity and the time and labor requirements. Therefore, we adapted and modified four physical dispersion methods with a gradient in dispersion intensity as follows:
1. Soils shaken in 100 ml of deionized water on a flatbed horizontal shaker for 25 min at 250 rpm [adapted and modified from van Wesemael et al. (31)]—This method is referred to as “H2O”
2. Soils shaken in 100 ml of deionized water on a flatbed horizontal shaker for 25 min at 250 rpm in presence of 10 glass beads with a diameter of 4 mm [adapted and modified from Rovira et al. (40) referred to as “Glassbeads”]
3. Ultrasonication (Sonopuls HD2200, Bandelin electronic GmbH, Germany) in 150 ml of deionized water with a low energy input of 100 J ml−1 (calorimetric calibration) [according to Hassink (21)] with a 13-mm-diameter sonotrode, an oscillation frequency of 20 kHz and an immersion depth of 15 mm (following Graf-Rosenfellner et al. (47)]–referred to as “US100”
4. Ultrasonication (Sonopuls HD2200, Bandelin electronic GmbH, Germany) in 150 ml deionized water with high energy input of 450 J ml−1 (calorimetric calibration) [according to Schmidt et al. (49)] with a 13–mm-diameter sonotrode, an oscillation frequency of 20 kHz and an immersion depth of 15 mm [according to Graf-Rosenfellner et al. (47)]–referred to as “US450”.
Figure 1. General scheme for a simplified fractionation to separate POM and MAOM with focus on the comparison of different dispersion procedures.
The initial experimental setup included a fifth dispersion method using sodium-hexametaphosphate (Na-HMP) as a chemical dispersion agent. However, this method was discontinued during the course of the experiment. The data acquired for this method was subsetted and separately evaluated from the other methods.
5. Chemical dispersion in 40 ml 0.5% Sodium-hexametaphosphate (NaPO3)6 shaken on a flatbed horizontal shaker [following (14, 45, 46)] referred to as “Na-HMP.”
In total, the main dataset (without the Na-HMP treatment) contained n = 80 samples [five sites × two vegetation types (C3/C4) × four dispersing methods × two replications], which were separated into fractions >20 and <20 μm. The Na-HMP subset consisted of eight samples dispersed with Na-HMP out of one site (Hohenheim C3) for additional analysis. The sieving process was modified from the other studies using this method. All samples were wet sieved using a 20 μm stainless steel sieve which was mounted on a horizontal shaking sieving machine (Analysette 03.502, Fritsch GmbH, Germany) with a frequency of ~50 Hz and an amplitude of 2–3 mm to standardize the sieving process.
Time Efficiency
Because time efficiency is of particular interest when working with large sample sets, we successively timed all actions that formed part of the fractionation schemes. We evaluated the time efficiency for all five dispersion methods. Additionally, we measured the time needed for manual wet sieving with a rubber spatula, automated wet sieving, and different drying procedures, including a drying cabinet, freeze dryers or rotary evaporators in a pre-experiment. We classified the required times into an active working time and a total time (comprising active working time and waiting time) needed for all fractionation steps. The measured time periods were averaged for the methods and rounded to 5 min intervals. A detailed scheme of the fractionation steps can be found in the Supplementary Figure 1.
Determination of C and δ13C
Total carbon (C) was measured in all fractions and bulk samples using dry combustion (HEKAtech EuroEA 3000). The measurement accuracy of the instrument was 2 mg C g−1 of soil. Total inorganic carbon (IC) was measured upon calcimetric reaction (4M HCl). We did not detect inorganic carbon in the samples and total C was taken as the organic C concentration. The δ13C (VPDB standard) was measured using an isotope ratio mass spectrometer (IRMS) (Delta V Advantage coupled to EA Isolink CN, Thermo Fischer Scientific, Waltham, USA), at the Thünen Institute isotope laboratory. For soils the accuracy of the system is δ 13C <0.5‰ based on routine analyses. All measurements were performed with at least two replicates, which is usually sufficient to reach the required accuracy as given above. Additional replicates were measured when the difference between the first two replicates reached thresholds of 2 mg C g−1 soil for total carbon and 0.5‰ δ13C for IRMS measurements.
Assessment of Performance Indicators
We assessed the incorporation of new C4-derived C into the fractions by using (19):
where fC4 is the factor of C4-derived C in the soil fractioni in the C4-vegetated soil and the δ13C values of the fractioni (C4 soil / C3 soil) and C4 plant refers to the δ13C values of the SOC fraction of interest in the C4 vegetated soil, the same SOC fraction in the C3 reference soil, and the δ13C value of the grown C4 plant. We fixed the δ13C reference value for the C4 plant at −12‰ (50). Thus, an fC4 of 1 corresponds to a sample or C fraction, where all C is new C4-plant derived C, whereas an fC4 of 0 corresponds to a C fraction with no input of C4-plant derived C. We used the range of the fC4 of the two fractions to quantify the degree of isolation for each dispersion treatment described as
where fC4max is the maximum fC4 value of both fractions of a soil while fC4min is the minimum value of both fractions of the same soil. As such, Range fC4 was calculated as the mean value for all replicates of a given method. Therefore, a high Range fC4 is an indicator for a well-performing fractionation method, while a low Range fC4 identifies a method that lacks performance in terms of isolating the C4/C3-C derived fractions within a soil (24). Additionally, we evaluated the mass loss/recovery of each method by adding up the mass contributions of the fractions divided by the initial mass of the sample to quantify the loss of soil during the application of our methods. Then. the mass recoveries were averaged for each method. Moreover, we evaluated the reproducibility of the methods by assessing the range of the mass contributions of a fraction between the replicates described as
where M is the mass of a corresponding fractioni and/or a given replicate or the total mass of the soil, i.e., the sum of the mass of both fractions of a given replicate, and |x| stands for the mathematical function of absolute values. The RangeMwas averaged for each method. Therefore, a smaller RangeMis an indicator of the reproducibility of the methods, where an increasing value indicates decreasing reproducibility. Additionally, we calculated C mass recovery and reproducibility by exchanging soil mass with soil C mass in the former equations.
The SOC loss caused by rinsing of the used Na-HMP was evaluated in the Na-HMP sample by comparing the SOC content against an unrinsed control group.
Additionally, we calculated the ∂13C-recovery as a cross check to verify the consistency and validity of our ∂13C measurements. The calculation and data can be seen in the Supplementary Material.
Statistical Analysis
All calculations were performed using the statistical software R. We used Kruskal-Wallis tests and analysis of variance models to test for average differences/influences in the Rangefc4 between the dispersion methods (H2O, Glassbeads, US100, US450) and the sites for significance at p < 0.05. Multiple comparisons were evaluated using the Tukey's Honestly Significance Difference test out of the package ”multcomp” (51). Moreover, we applied Fisher's t-test to test for mean value differences within the proportions of C4-derived carbon of the fractions <20 and >20 μm. The effects of Na-HMP induced SOC leaching were evaluated using Welsh's t-test for data with unequal variance. All data were checked for normality and homoscedasticity using a ShapiroWilk test on either the variables or the residues. We used a significance level of p < 0.05 for all tests.
Results
C4-Derived Organic Carbon Proportion in Fractions
The proportion of C4-derived organic carbon was significantly different (t = 5.9508, df = 33.696, p < 0.001) between the fractions >20 and <20 μm. Mean C4-derived OC was 26% in the fraction <20 μm and 54% in the fraction >20 μm (Figure 2). Mean C4-derived OC was dependent on the dispersion method ranging from 46 and 47%, respectively, for the Glassbeads and H2O treatments to 58 and 63%, respectively, for the US100 and US450 treatments in the fraction >20 μm. In the fraction <20 μm, it ranged from 22 and 23%, respectively, for H2O and Glassbeads treatments to 26 and 31% for the US100 and US450 treatments. Furthermore, despite there being a distinct increase in C4-derived C in the fractions > 20 μm, we also measured a slight shift of C4-derived C into the fractions <20 μm in almost all soils with increasing intensity of the used dispersion methods. The Rangefc4 was very constant for the dispersion methods in the soils of Braunschweig (0.20 ± 0.01), Jyndevad (0.22 ± 0.02), and Rotthalmünster (0.20 ± 0.01) as indicated by small standard deviations. In Hohenheim (0.27 ± 0.06), the US100 dispersion reached a higher Rangefc4 of 0.36 ± 0.06 compared to H2O 0.21 ± 0.06, Glassbeads 0.26 ± 0.01, and US450 0.25 ± 0.03 (Figure 3). The site with the highest clay content Tänikon (0.51 ± 0.17) had the highest differentiation between the dispersion methods. H2O had a Rangefc4 of 0.42 ± 0.05, Glassbeads 0.28 ± 0.06 and the fractions of US100 and US450 reached a higher grade of isolation with 0.65 ± 0.01 and 0.70 ± 0.02, respectively. The average Rangefc4 of the methods was 0.25 ± 0.09 for H2O, 0.23 ± 0.04 for Glassbeads, 0.33 ± 0.17 for US100 and 0.32 ± 0.19 for US450 (Figure 3). There were no significant differences in the Rangefc4 when comparing the mean Rangefc4 of the methods without differentiating by the sites (Kruskal-Wallis rank-sum test). However, there are highly significant influences (ANOVA) on the Rangefc4 of the methods (p < 0.001) when adding the sites as grouping variable (p < 0.001) as well as the interaction of sites and methods (p < 0.001). A subsequent Tukey-HSD test identified the Rangefc4 of the sites as significantly different for three groups (Figure 3).
Figure 2. Proportion of C4-derived C of the fractions >20 and <20 μm depicted for each dispersion method and site.
Figure 3. Rangefc4 of the fractions >20 and <20 μm depicted for each dispersion method and site. Different letters indicate significant differences at p < 0.05.
Dispersion Effects on Organic Carbon, C/N Ratio, and Mass Balance of Soil Fractions
The OC contents tended to slightly decrease in the >20 μm fraction for all methods with increasing intensity of dispersion, while the contents within the <20 μm fraction increased by roughly equal amounts (Figure 4). The soils with low clay content tended to have higher OC contents in the fraction <20 μm than in the >20 μm fraction. Stronger dispersion led to changes in the mass balances of the fraction. The mass of the >20 μm fraction decreased with increasing dispersion intensity for all tested soils whilst the <20 μm fraction increased (Figure 5). The differences were more pronounced in the silty soils of Hohenheim and Rotthalmünster and the clayey soil of Tänikon than the sandy soils of Braunschweig and Jyndevad. This shift from the >20 μm fraction to the <20 μm fraction with increasing dispersion intensity was also measurable for SOC when comparing the OC mass contribution of the fractions between the tested dispersion methods (Supplementary Table 1, Supplementary Figure 2).
Figure 4. Organic carbon (OC) content of the fractions >20 and <20 μm depicted for each dispersion method and site and their mean.
Figure 5. Mean mass proportion of the fractions >20 and <20 μm and loss of mass depicted for each experimental site and used dispersion method (H2O, Glassbeads, US100, US450).
The mean C/N ratio was higher in the >20 μm fraction (17 ± 0.3) than the <20 μm fraction (10 ± 0.1) (Figure 6). In the soil from Jyndevad, the C/N ratio did not differ across the dispersion methods of each fraction. However, there were differences between the dispersion methods in the >20 μm fractions of the soils from Braunschweig, Rotthalmünster, Hohenheim, and Tänikon. While the C/N ratio of the >20 μm fraction from these sites was 11–15 in samples treated with water or water and glass beads, it was higher using ultrasonication. The C/N ratio ranged between 14 and 21 for the US100 method and 16–29 using the US450 method.
Figure 6. C/N ratio of the fractions >20 and <20 μm depicted for each dispersion method and site and their mean.
The mean mass loss was <2% for all dispersion methods (Figure 7). With increasing intensity of the dispersion method, the mean mass loss increased slightly yet not significantly (p > 0.05). The reproducibility of the H2O method was lower than the other methods indicated by a remarkably higher RangeM of 18.25 mg g−1 than the “Glassbeads” (7.38 mg g−1), US100 (5.64 mg g−1) and US450 (5.98 mg g−1) methods (Figure 7). The mean OC loss ranged from 3.6 ± 2.9% to 4.3 ± 5.0 %, which is equivalent to 95.7 % to 96.8 % C recovery, and was not significantly different between the methods (p > 0.05). The reproducibility of the OC content within the fractions marginally differed between the methods. The RangeC between the replicates of each method was 0.94–1.35 mg g−1 C, and it was smaller than the mean RangeC of repeated bulk soil measurements (1.88 mg g−1 C).
Figure 7. Mass loss (left top) and C loss (left bottom) shown as indicators for mass and C recovery and Range mass/Range OC shown as indicators for mass/C reproducibility for the used dispersions and the mean of all dispersions (H2O, Glassbeads, US100, US450, mean). Low mass/C loss report high recovery whereas low Range C/mass report high reproducibility.
Workload and Procedural Time
For the dispersion methods, the working time ranged from 5 min for the H2O method to 15 min for the Glassbeads, US100 and US450 methods (Table 2). The total time until completing one sample was 15 min for US100, 30 min for US450, 30 min for H2O, and 35 min using the glass beads. The duration of the sieving process depended on soil texture, taking 30–60 min total time when using the sieving machine and 30–75 min when sieving by hand.
Table 2. Mean values of the quality indicators mass loss (recovery)/reproducibility, C loss (recovery)/reproducibility, RangefC4 and active, and total work time for each dispersion method.
Sodium Hexametaphosphate Subset Measurements
Chemical dispersion using Na-HMP without rinsing includes the risk of damaging the equipment or increased costs due to the higher wear and tear of the hose system within the elemental analyzer, especially with the high sample throughput as intended. For this reason, we evaluated the Na-HMP dispersion with subsequent rinsing. Rinsing the samples down to an electrical conductivity (EC) of <5 μS cm−1 took at least 2 days and was therefore not suitable for a fast application. Additionally, we measured a significant SOC difference (*, Welch's t-test, t = 3.7373, df = 3.1485, p = 0.03074) in SOC contents between the rinsed samples and the control soils. The mean SOC contents were 79.8 ± 1.5 mg g−1 in the control and 61.4 ± 9.8 mg g−1 in the rinsed soils.
Discussion
Isolation Efficiencies of Dispersion Methods Are Related to Soil Texture
All dispersion methods except the Na-HMP method revealed a high C recovery ranging from 95 to 97%. A recent study that compared different fractionation methods (24) quantified the recovery of bulk SOC to an average of 95 ± 4%. Although it had the highest RangeM among the methods, the H2O method was the less reproducible than the other methods (Figure 7), which indicates that this method lacks the mechanic power to disperse the soil adequately. Our data indicated that the ultrasonication dispersions performed better when comparing the efficiency to isolate fractions with distinct turnover rates. However, these findings are only significant when taking the site, including the specific texture, vegetation, and age of the experiment, and the interaction of the site and dispersion method into account. The results showed that the site with the highest clay content of 29% (Tänikon) required a stronger dispersion with ultrasonication than the other sites to achieve fractions with better isolation of C3- and C4-derived organic carbon. Clay particles tend to form very stable microaggregates (52–54) and increase the proportion of macro-aggregates in the soil (54, 55). These macroaggregates obviously resisted weaker dispersion methods with H2O and glass beads (19, 31). Our results showed that the parameter for the isolation quality (Rangefc4) is largely based on the increase of the moieties of C4-derived C and the C/N ratio of the > 20 μm fraction (Figure 6, Supplementary Figure 3). In the <20 μm fraction, the proportion of C4-derived C and the C/N ratio were hardly changed with increasing dispersion intensity (Figure 6, Supplementary Figure 3). As the <20 μm fraction remained almost unchanged with the increasing intensity of dispersion, we conclude that the most important criterion for a successful isolation of fractions with different turnover is the degree of purity of C4-derived OC in the >20 μm fraction which is determined by the clay content of the respective soil (Supplementary Figure 4). The higher C/N ratio of the > 20 μm fractions in the clayey soils with stronger dispersion is also explained by the destruction of microaggregates and the shift of MAOM into the <20 μm fraction. Moreover, differences in the OC content of the >20 μm fraction when using H2O and US450 might also be used to calculate the amount of C in the “occluded/aggregate protected” POM fraction (o-POM), obtaining data similar to that from density fractionation methods (14, 56), but avoiding their high cost. Such an approach would need to be tested to check the uncertainty of the results obtained from calculation compared to density fractionation.
Comparison of Isolation Qualities
We used the Rangefc4 indicator to compare the isolation quality of the investigated dispersion methods with methods investigated in a former study (24) that used soils of mostly the same long-term experiments. One of the tested fractionation methods within that study (40) combined the use of ultrasonication and glass beads but separated soils into four fractions (2,000–200, 200–50, 50–20, and <20 μm) and was thus methodically comparable to our approach. Despite being more elaborate, this method had no isolation advantage with a mean Rangefc4 indicator of 0.22 compared to 0.23–0.33 in our methods. Other methods evaluated in this study that had evidently higher Rangefc4 [(17, 19, 57–59); in Poeplau et al. (24)] all used chemical agents either for density fractionation (sodium polytungstate), alkaline extraction (sodium diphosphate) or OM oxidation (H2O2, NaOCl2). Besides higher work load through acquisition of additional fractions that have to be investigated and measured separately, chemical agents might change the chemical composition of POM (60, 61). However, water has no effect on the chemistry of POM (60). Additionally, to our knowledge, there is no evidence of ultrasonication altering the chemistry of POM.
Limited Performance of Chemical Dispersion With Sodium-Hexametaphosphate
Although Na-HMP is widely used for soil dispersion (16, 46, 62), our attempts to use it within a simplified fractionation scheme were not successful. Na-HMP needs to be rinsed out or maintained within the samples. It has a melting point of 628°C, and in the presence of moisture it is hydrolyzed to Na-trimetaphosphate with a melting point of 53°C. During combustion in an elemental analyzer or IRMS, it will dissolve and precipitate within devices, and it may damage equipment, especially when the number of measurements is high. Therefore, most laboratories filtrate the soil and rinse the salt with deionized water. Rinsing the Na-HMP increasing the work load, depending on soil texture and OC content. A sample required an average time of 2–8 days to be rinsed before reaching an EC <5 μS cm−1. Rinsing resulted in significant SOC extraction and thus SOC loss. This concurs with another study revealing critical changes in soils due to the use of sodium salts and rinsing (63). Therefore, we cannot recommend using Na-HMP as a fast and simple dispersion method. Furthermore, our results imply that the OC recovery rate is indispensable to evaluate the outcome of OC fractionation methods.
Evaluation of Performance Indicators and Time Expenditure
Our results showed that all dispersion methods have trade-offs in either isolation qualities, equipment and workload requirements or recovery and reproducibility criteria. In terms of the workload, Poeplau et al. (24) emphasized the correlation of increasing workloads with higher numbers of fractions. Additionally, the use of chemical agents significantly extends the fractionation process, as long as these chemicals are removed from the fractions afterwards. All dispersion methods tested, except the Na-HMP treatment, had a relatively short working time. The total time requirements of the methods was slightly shorter for ultrasonication than the H2O and glassbeads methods (Table 3). Costs depend on the workload and equipment used and can be divided into running costs and fixed costs. Ultrasonication devices have relatively high fixed costs and the sonotrodes need to be changed regularly, causing relatively higher running costs than other methods. Because costs are strongly related to market volatility and country-specific wage costs, we do not evaluate these costs in detail. Another vital factor regarding time expenditure is the amount of water used during wet sieving. Depending on soil texture, we used 3–10 l of deionized water, as we wanted to achieve the best possible separation indicated by rinsed water without any remaining colloids after sieving. However, other studies suggest using the lowest amount of deionized water possible to prevent DOM leaching (48) and recommend the use of 2 l of deionized water. We assume that both effects, either the leaching of DOM, or the insufficient translocation of < 20 μm particles after the use of only 2,000 ml of deionized water, will most likely not affect the OC balance significantly. Automated sieving is slightly faster than manual sieving. However, this effect diminishes with increasing difference in the mean particle size of the soil to the mesh size and recovery and reproducibility might differ between both methods. In terms of the time expenditures of the drying methods, drying cabinets are advantageous, but the drying process itself takes several days depending on the volume of water requiring evaporation and the power of the device. Big rotary evaporators working at low pressure are capable of evaporating large volumes of water, such as 10 l within 4 h in our case, but they also come with a higher workload due to the cumbersome cleaning of the vessels. Freeze dryers require both intermediate to high workloads and long drying times, but they have the advantage that they simplify the recovery from vessels. Therefore, they are beneficial when working with small amounts of soil. A complete scheme and timeframe of the experiment and the fractionation steps can be seen in Supplementary Figure 1.
Table 3. Evaluation of quality and efficiency criteria describing the performance of each work step of the modules dispersion, wet sieving, and drying.
Limitations
This study focused on agricultural soils in Central Europe. Although the methods might be usable on soils under different conditions, they were calibrated to isolate and quantify POM and MAOM content/stocks under criteria specific to Central Europe. Therefore, we cannot recommend them for use on soils with divergent land management or origin without former method calibration and optimization. Additionally, their efficiencies could be different for soils with a clay mineral composition different to that of central European soils, carbonaceous soils, organic soils, or OM-rich soils such as forest soils or soils with special conditions related to organic or inorganic carbon. Moreover, we did not compare different mesh sizes within the sieving process, such as 50 and 63 μm. Furthermore, we want to emphasize that the number of replicates (n = 2) used in this study is low, especially with regard to the reproducibility indicator, which is causing uncertainty. Nevertheless, we consider the number of replicates as adequate for the precision of the approach indicated by high reproducibility, low variance and high ∂13C-recovery as outlined in Figures 2, 7 and Supplementary Table 2.
Conclusions
Our study showed that a simple soil fractionation into two particle size fractions of >20 and <20 μm is a functional tool to assess information on SOC fractions with different turnover rates. However, the differentiation in SOC pools with distinct turnover rates strongly depends on soil texture and the dispersion method used. Clay-rich soils required higher dispersion intensities whereas soils with lower clay contents showed no differences in their dispersion methods. Ultrasonication with 450 J ml−1 has the overall best quality indicators in isolation of two fractions with different turnover time and in reproducibility compared to dispersion methods using water with or without glass beads and is therefore best suited for a simple fractionation. It enables a relatively quick and inexpensive, but nevertheless precise, assessment of functional SOM pools, which enables the estimation of management influences on soil functions such as nutrient turnover (POM) or C sequestration (MAOM). Due to the substantial potential of simple fractionation to estimate influences of cultivation or land use on associated soil functions, application in agriculture is promising, especially on a larger spatial (regional, national) scale with a large sample size. In particular, the method could be implemented within existing soil monitoring programs. Additionally, there is the option of using fractionation of archived soil samples to trace the development of the fast and slow cycling C pools of soils up to the present day.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
CJ conducted the experiments and prepared the first manuscript draft. CJ and MW analyzed the data. All authors conceived the experiment and contributed to its refinement.
Funding
This project work was funded by the German Federal Ministry of Education and Research (BMBF) in the framework of the funding measure Soil as a Sustainable Resource for the Bioeconomy—BonaRes, project BonaRes (Module B): BonaRes Center for Soil Research, subproject C (Grant No. 031B0511C).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are grateful to Anette Goeske and the staff of the Stable Isotope Laboratory of Thünen Institute for IRMS measurements and valuable advice.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fsoil.2021.692583/full#supplementary-material
References
1. Lal R. Soil health and carbon management. Food Energy Secur. (2016) 5:212–222. doi: 10.1002/fes3.96
2. Paustian K, Andrén O, Janzen HH, Lal R, Smith P, Tian G, et al. Agricultural soils as a sink. Soil Use Manag. (1997) 13:230–44. doi: 10.1111/j.1475-2743.1997.tb00594.x
3. Wiesmeier M, Urbanski L, Hobley E, Lang B, von Lützow M, Marin-Spiotta E, et al. Soil organic carbon storage as a key function of soils - a review of drivers and indicators at various scales. Geoderma. (2019) 333:149–62. doi: 10.1016/j.geoderma.2018.07.026
4. Ciais P, Sabine C, Bala G, Bopp L, Brovkin V, Canadell J, et al. Carbon and Other Biogeochemical Cycles. In: Stocker TF, Qin D, Plattner G-K, Tignor M, Allen SK, Boschung J, et al. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (New York, NY: Cambridge University Press).
5. Jarecki MK, Lal R. Crop management for soil carbon sequestration. Crit Rev Plant Sci. (2003) 22:471–502. doi: 10.1080/713608318
6. Vogel HJ, Eberhardt E, Franko U, Lang B, Lie ßM, Weller U, et al. Quantitative evaluation of soil functions: potential and state. Front Environ Sci. (2019) 7:164. doi: 10.3389/fenvs.2019.00164
7. Oades JM. The retention of organic matter in soils. Biogeochemistry. (1988) 5:35–70. doi: 10.1007/BF02180317
8. Moni C, Derrien D, Hatton PJ, Zeller B, Kleber M. Density fractions versus size separates: does physical fractionation isolate functional soil compartments? Biogeosciences Discuss. (2012) 9:8405–47. doi: 10.5194/bgd-9-8405-2012
9. Chenu C, Rumpel C, Lehmann J. Methods for Studying Soil Organic Matter. 4th ed. Waltham: Elsevier Inc (2015). doi: 10.1016/B978-0-12-415955-6.00013-X
10. Six J, Elliott ET, Paustian K. Aggregate and soil organic matter dynamics under conventional and No-tillage systems. Soil Sci Soc Am J. (1999) 63:1350–1358. doi: 10.2136/sssaj1999.6351350x
11. Helfrich M, Ludwig B, Buurman P, Flessa H. Effect of land use on the composition of soil organic matter in density and aggregate fractions as revealed by solid-state 13C NMR spectroscopy. Geoderma. (2006) 136:331–41. doi: 10.1016/j.geoderma.2006.03.048
12. Kravchenko AN, Wang ANW, Smucker AJM, Rivers ML. Long-term differences in tillage and land use affect intra-aggregate pore heterogeneity. Soil Sci Soc Am J. (2011) 75:1658–66. doi: 10.2136/sssaj2011.0096
13. von Lützow M, Kögel-Knabner I, Ekschmitt K, Flessa H, Guggenberger G, Matzner E, et al. SOM fractionation methods: Relevance to functional pools and to stabilization mechanisms. Soil Biol Biochem. (2007) 39:2183–207. doi: 10.1016/j.soilbio.2007.03.007
14. Cambardella CA, Elliot ET. Particulate soil organic-matter changes across a grassland cultivation sequence. Soil Sci Sock Am J. (1992) 56:777–83. doi: 10.2136/sssaj1992.03615995005600030017x
15. Christensen BT. Physical fractionation of soil and organic matter in primary particle size and density separates. Adv Soil Sci. (1992) 20:1–90. doi: 10.1007/978-1-4612-2930-8_1
16. Six J, Elliott ET, Paustian K, Doran JW. Aggregation and soil organic matter accumulation in cultivated and native grassland soils. Soil Sci Soc Am J. (1998) 62:1367–77. doi: 10.2136/sssaj1998.03615995006200050032x
17. Zimmermann M, Leifeld J, Schmidt MWI, Smith P, Fuhrer J. Measured soil organic matter fractions can be related to pools in the RothC model. Eur J Soil Sci. (2007) 58:658–67. doi: 10.1111/j.1365-2389.2006.00855.x
18. Steffens M, Kölbl A, Kögel-Knabner I. Alteration of soil organic matter pools and aggregation in semi-arid steppe topsoils as driven by organic matter input. Eur J Soil Sci. (2009) 60:198–212. doi: 10.1111/j.1365-2389.2008.01104.x
19. Balesdent J. The turnover of soil organic fractions estimated by radiocarbon dating. Sci Total Environ. (1987) 62:405–8. doi: 10.1016/0048-9697(87)90528-6
20. Balesdent J. The significance of organic separates to carbon dynamics and its modelling in some cultivated soils. Eur J Soil Sci. (1996) 47:485–93. doi: 10.1111/j.1365-2389.1996.tb01848.x
21. Hassink J. A model of the physical protection of organic matter in soils the capacity of soils to preserve organic C and N by their association with clay and silt particles. Plant Soil. (1997) 191:77–87. doi: 10.1023/A:1004213929699
22. Amelung W, Bol R, Friedrich C. Natural 13C abundance: A tool to trace the incorporation of dung-derived carbon into soil particle-size fractions. Rapid Commun Mass Spectrom. (1999) 13:1291–4. doi: 10.1002/(SICI)1097-0231(19990715)13:13<1291::AID-RCM637>3.0.CO;2-C
23. Rumpel C, Kögel-Knabner I, Knicker H, Hüttl RF. Composition and distribution of organic matter in physical fractions of a rehabilitated mine soil rich in lignite-derived carbon. Geoderma. (2000) 98:177–92. doi: 10.1016/S0016-7061(00)00060-4
24. Poeplau C, Don A, Six J, Kaiser M, Benbi D, Chenu C, et al. Isolating organic carbon fractions with varying turnover rates in temperate agricultural soils – a comprehensive method comparison. Soil Biol Biochem. (2018) 125:10–26. doi: 10.1016/j.soilbio.2018.06.025
25. Christensen BT. Matching measurable soil organic matter fractions with conceptual pools in simulation models of carbon turnover: revision of model structure. Eval Soil Org Matter Model. (1996) 38:143–159. doi: 10.1007/978-3-642-61094-3_11
26. Cotrufo MF, Ranalli MG, Haddix ML, Six J, Lugato E. Soil carbon storage informed by particulate and mineral-associated organic matter. Nat Geosci. (2019) 12:989–94. doi: 10.1038/s41561-019-0484-6
27. Lavallee JM, Soong JL, Cotrufo MF. Conceptualizing soil organic matter into particulate and mineral-associated forms to address global change in the 21st century. Glob Chang Biol. (2019) 26:261–273. doi: 10.1111/gcb.14859
28. von Lützow M, Kögel-Knabner I, Ludwig B, Matzner E, Flessa H, Ekschmitt K, et al. Stabilization mechanisms of organic matter in four temperate soils: development and application of a conceptual model. J Plant Nutr Soil Sci. (2008) 171:111–24. doi: 10.1002/jpln.200700047
29. Parton WJ, Stewart JWB, Cole CV. Dynamics of C. N, P and S in grassland soils: a model. Biogeochemistry. (1988) 5:109–31. doi: 10.1007/BF02180320
30. Bonde TA, Christensen BT, Cerri CC. Dynamics of soil organic matter as reflected by natural 13C abundance in particle size fractions of forested and cultivated oxisols. Soil Biol Biochem. (1992) 24:275–7. doi: 10.1016/0038-0717(92)90230-U
31. van Wesemael B, Chartin C, Wiesmeier M, von Lützow M, Hobley E, Carnol M, et al. An indicator for organic matter dynamics in temperate agricultural soils. Agric Ecosyst Environ. (2019) 274:62–75. doi: 10.1016/j.agee.2019.01.005
32. Six J, Conant RT, Paul EA, Paustian K. 2002 six stabilizatio nmechanisms of SOM implications for C saturation of soils.pdf. Plant Soil. (2002) 241:155–76. doi: 10.1023/A:1016125726789
33. Stewart CE, Paustian K, Conant RT, Plante AF, Six J. Soil carbon saturation: concept, evidence and evaluation. Biogeochemistry. (2007) 86:19–31. doi: 10.1007/s10533-007-9140-0
34. Castellano MJ, Mueller KE, Olk DC, Sawyer JE, Six J. Integrating plant litter quality, soil organic matter stabilization, and the carbon saturation concept. Glob Chang Biol. (2015) 21:3200–9. doi: 10.1111/gcb.12982
35. Coleman K, Jenkinson DS. RothC-26.3 - a model for the turnover of carbon in soil. Eval Soil Org Matter Model. (1996) 38:237–46. doi: 10.1007/978-3-642-61094-3_17
36. Chen S, Arrouays D, Angers DA, Chenu C, Barré P, Martin MP, et al. National estimation of soil organic carbon storage potential for arable soils: a data-driven approach coupled with carbon-landscape zones. Sci Total Environ. (2019) 666:355–67. doi: 10.1016/j.scitotenv.2019.02.249
37. Dechow R, Franko U, Kätterer T, Kolbe H. Evaluation of the RothC model as a prognostic tool for the prediction of SOC trends in response to management practices on arable land. Geoderma. (2019) 337:463–78. doi: 10.1016/j.geoderma.2018.10.001
38. Chenu C, Plante AT. Clay-sized organo-mineral complexes in a cultivation chronosequence: revisiting the concept of the “primary organo-mineral complex.” Eur J Soil Sci. (2006) 57:596–607. doi: 10.1111/j.1365-2389.2006.00834.x
39. Plante AF, Conant RT, Stewart CE, Paustian K, Six J. Impact of soil texture on the distribution of soil organic matter in physical and chemical fractions. Soil Sci Soc Am J. (2006) 70:287–96. doi: 10.2136/sssaj2004.0363
40. Rovira P, Jorba M, Romanyà J. Active and passive organic matter fractions in mediterranean forest soils. Biol Fertil Soils. (2010) 46:355–69. doi: 10.1007/s00374-009-0437-0
41. Tiessen H, Stewart JWB. Particle-size fractions and their use in studies of soil organic matter: II. Cultivation effects on organic matter composition in size fractions. Soil Sci Soc Am J. (1983) 47:509–14. doi: 10.2136/sssaj1983.03615995004700030023x
42. Kölbl A, Leifeld J, Kögel-Knabner I. A comparison of two methods for the isolation of free and occluded participate organic matter. J Plant Nutr Soil Sci. (2005) 168:660–7. doi: 10.1002/jpln.200521805
43. Kaiser M, Berhe AA, Sommer M, Kleber M. Application of ultrasound to disperse soil aggregates of high mechanical stability. J Plant Nutr Soil Sci. (2012) 175:521–6. doi: 10.1002/jpln.201200077
44. Tyner EH. The use of sodium metaphosphate for dispersion of soils for mechanical analysis. Soil Sci Soc Am J. (1940) 4:106–13. doi: 10.2136/sssaj1940.036159950004000C0020x
45. Six J, Merckx R, Kimpe K, Paustian K, Elliott ET. A re-evaluation of the enriched labile soil organic matter fraction. Eur J Soil Sci. (2000) 51:283–93. doi: 10.1046/j.1365-2389.2000.00304.x
46. Sanderman J, Fillery IRP, Jongepier R, Massalsky A, Roper MM, MacDonald LM, et al. Carbon sequestration under subtropical perennial pastures II: carbon dynamics. Soil Res. (2013) 51:771–80. doi: 10.1071/SR12351
47. Graf-Rosenfellner M, Kayser G, Guggenberger G, Kaiser K, Büks F, Kaiser M, et al. Replicability of aggregate disruption by sonication—an inter-laboratory test using three different soils from Germany. J Plant Nutr Soil Sci. (2018) 181:894–904. doi: 10.1002/jpln.201800152
48. Poeplau C, Don A. Sensitivity of soil organic carbon stocks and fractions to different land-use changes across Europe. Geoderma. (2013) 192:189–201. doi: 10.1016/j.geoderma.2012.08.003
49. Schmidt MWI, Rumpel C, Kögel-Knabner I. Evaluation of an ultrasonic dispersion procedure to isolate primary organomineral complexes from soils. Eur J Soil Sci. (1999) 50:87–94. doi: 10.1046/j.1365-2389.1999.00211.x
50. Menichetti L, Ekblad A, Kätterer T. Organic amendments affect δ13C signature of soil respiration and soil organic C accumulation in a long-term field experiment in Sweden. Eur J Soil Sci. (2013) 64:621–8. doi: 10.1111/ejss.12077
51. Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biometrical J. (2008) 50:346–63. doi: 10.1002/bimj.200810425
52. Wagner S, Cattle SR, Scholten T. Soil-aggregate formation as influenced by clay content and organic-matter amendment. J Plant Nutr Soil Sci. (2007) 170:173–80. doi: 10.1002/jpln.200521732
53. Kaiser K, Kalbitz K. Cycling downwards - dissolved organic matter in soils. Soil Biol Biochem. (2012) 52:29–32. doi: 10.1016/j.soilbio.2012.04.002
54. Krause L, Rodionov A, Schweizer SA, Siebers N, Lehndorff E, Klumpp E, et al. Microaggregate stability and storage of organic carbon is affected by clay content in arable Luvisols. Soil Tillage Res. (2018) 182:123–9. doi: 10.1016/j.still.2018.05.003
55. Schweizer SA, Bucka FB, Graf-Rosenfellner M, Kögel-Knabner I. Soil microaggregate size composition and organic matter distribution as affected by clay content. Geoderma. (2019) 355:113901. doi: 10.1016/j.geoderma.2019.113901
56. Christensen BT. Physical fractionation of soil and structural and functional complexity in organic matter turnover. Eur J Soil Sci. (2001) 52:345–53. doi: 10.1046/j.1365-2389.2001.00417.x
57. Leifeld J, Kögel-Knabner I. Organic carbon and nitrogen in fine soil fractions after treatment with hydrogen peroxide. Soil Biol Biochem. (2001) 33:2155–8. doi: 10.1016/S0038-0717(01)00127-4
58. Sollins P, Kramer MG, Swanston C, Lajtha K, Filley T, Aufdenkampe AK, et al. Sequential density fractionation across soils of contrasting mineralogy: evidence for both microbial- and mineral-controlled soil organic matter stabilization. Biogeochemistry. (2009) 96:209–31. doi: 10.1007/s10533-009-9359-z
59. Kaiser M, Zederer DP, Ellerbrock RH, Sommer M, Ludwig B. Effects of mineral characteristics on content, composition, and stability of organic matter fractions separated from seven forest topsoils of different pedogenesis. Geoderma. (2016) 263:1–7. doi: 10.1016/j.geoderma.2015.08.029
60. Chen X, Jin M, Xu Y, Chu W, Olk DC, Hu J, et al. Potential alterations in the chemical structure of soil organic matter components during sodium hydroxide extraction. J Environ. Qual. (2019) 48:1578–86. doi: 10.2134/jeq2019.02.0077
61. Ohno T, Hess NJ, Qafoku NP. Current understanding of the use of alkaline extractions of soils to investigate soil organic matter and environmental processes. J Environ Qual. (2019) 48:1561–4. doi: 10.2134/jeq2019.08.0292
62. Robertson AD, Paustian K, Ogle S, Wallenstein MD, Lugato E, Cotrufo MF. Unifying soil organic matter formation and persistence frameworks: the MEMS model. Biogeosci Discuss. (2018) 16:1–36. doi: 10.5194/bg-2018-430
Keywords: soil organic carbon, carbon sequestration, soil fractions, agricultural soils, particulate organic matter, mineral-associated organic matter
Citation: Just C, Poeplau C, Don A, van Wesemael B, Kögel-Knabner I and Wiesmeier M (2021) A Simple Approach to Isolate Slow and Fast Cycling Organic Carbon Fractions in Central European Soils—Importance of Dispersion Method. Front. Soil Sci. 1:692583. doi: 10.3389/fsoil.2021.692583
Received: 04 August 2021; Accepted: 03 September 2021;
Published: 08 October 2021.
Edited by:
Katell Quenea, Sorbonne Universités, FranceReviewed by:
Rota Wagai, National Agriculture and Food Research Organization (NARO), JapanWenting Feng, Institute of Agricultural Resources and Regional Planning (CAAS), China
Copyright © 2021 Just, Poeplau, Don, van Wesemael, Kögel-Knabner and Wiesmeier. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Christopher Just, Yy5oZW5rZSYjeDAwMDQwO3R1bS5kZQ==
†ORCID: Christopher Just 0000-0001-9727-5470
Christopher Poeplau 0000-0003-3108-8810
Axel Don 0000-0001-7046-3332
Bas van Wesemael 0000-0003-4007-0241
Ingrid Kögel-Knabner 0000-0002-7216-8326
Martin Wiesmeier 0000-0003-3981-5461