Mean nutrient uptake depths of cereal crops change with compost incorporation into subsoil – evidence from 87Sr/86Sr ratios

Root restricting layers often hinder crops from accessing the large reservoir of bioavailable mineral nutrients situated in subsoil. This study aims to explore changes in the mean nutrient uptake depth of cereal crops when removing root restricting layers through subsoil management. Subsoil management was performed by deep loosening, cultivation of lucerne as deep-rooting pre-crop, and their combination with compost incorporation. Management effects were evaluated by means of shoot biomass and element concentrations in shoots and soil compartments. The mean nutrient uptake depth was fingerprinted by graphically matching the 87Sr/86Sr ratios in shoots with the 87Sr/86Sr ratios in the exchangeable fraction in soil. Nutrient uplift from depth to topsoil was inferred from element concentrations in the exchangeable fraction in soil. Shoot biomass remained constant in management and control plots. The mean nutrient uptake depth changed with subsoil management in the order: deep loosening < control < deep loosening with compost incorporation. The latter coincided with a reallocation of compost-derived Na and hence resulted in increased levels of bioavailable Na below the depth of compost incorporation, which may have led to an improved water use efficiency of the crops. Thus, Na relocation triggered the deepening of the mean uptake depth of water and nutrients. Moreover, nutrient uplift from depth to topsoil was evident 21 months after subsoiling. Subsoil management by deep loosening with compost incorporation provides a sustainable use of soil resources because otherwise unused deep geogenic-derived nutrient reservoirs were additionally involved in crop nutrition.


Introduction
The global population is rising to about 10 billion people by 2050 (Adam 2021). This population growth dramatically increases human demand for agricultural products. From 1961 to 2017, for example, global cereal production was boosted by about 240% (IPCC 2019). Negative effects of this development are omnipresent and include land use changes like deforestation, groundwater pollution, and freshwater depletion (Campbell et al. 2017). To avoid further conversions of natural ecosystems into intensely managed agricultural land and to cope with prolonged summer droughts resulting from climate warming, an efficient use of soil resources is urgently required.
Traditionally, the productivity of agricultural land was almost exclusively maintained in the plough horizon, the agricultural topsoil. However, the subsoil beneath this tilled layer contains an inherent large reservoir of water and mineral nutrients in bioavailable form (Kirkegaard et al. 2007;Kautz et al. 2013). Given the global mean maximum rooting depth of crops of 200 cm (Canadell et al. 1996) these resources are principally accessible to crops. However, the most common arable crops have shallower maximum rooting depths ranging from about 50 cm to 150 cm depth (Fan et al. 2016). Thus, projects like DeepFrontier (Thorup-Kristensen et al. 2020) aimed at developing deep rooting crops by breeding (Wasson et al. 2012), or by improving deep rooting through crop management such as early sowing (Rasmussen and Thorup-Kristensen 2016) and crop rotation (Thorup-Kristensen and Kirkegaard 2016). Further suggestions on the exploitation of deep soil layers by crops include a change from annual to perennial crop species (Glover et al. 2010;Crews and Cattani 2018), planting deep rooting pre-crops on arable soil Seidel et al. 2019), and agroforestry, where shallow rooted crops are intercropped with deep rooted species (Cardinael et al. 2015).
However, many soils have root restricting layers in the upper meter of soil, which affect, for example, up to 71% of Germany's agricultural soils (Schneider and Don 2019a). Effective means to reduce the penetration resistance for crops include physical melioration such as deep loosening (Schneider and Don 2019b). Deep loosening enhances not only root growth into the subsoil (Jakobs et al. 2019) but also increases water infiltration (Blanco-Canqui et al. 2017), which promotes crop yield especially in dry years (Marks and Soane 1987;Olesen and Munkholm 2007;Gaiser et al. 2012). Deep loosening combined with the incorporation of compost improves grain yield (Jakobs et al. 2019). Whether deep loosening also results in a deepening of the uptake depth of water and nutrients, and even causes nutrient uplift, remains elusive. Nutrient uplift, as a result of uptake of nutrients from the mineral soil and their return to the organic topsoil via multiple pathways, leaves behind a characteristic depth profile of increasing concentrations of the most plant essential nutrient elements in the exchangeable and extractable fraction of soil from depth to topsoil (Jobbágy and Jackson 2001). Thus, even though nutrient uplift can easily be detected by means of concentration measurements, it remains unknown from which depth uplifted nutrients are originally sourced. Given that root length measurements are highly demanding, an alternative potential analytical tool is a geochemical proxy such as the radiogenic strontium isotope ratio 87 Sr/ 86 Sr.
Strontium has four stable isotopes ( 84 Sr, 86 Sr, 87 Sr, 88 Sr), among which 87 Sr results from the radioactive β − decay of 87 Rb (half-life time = 4.88 × 10 10 yr −1 (Holden 1990)), eponymous for this 'radiogenic' isotope system. Unlike for non-traditional stable isotopes, any natural and instrumental mass-dependent isotope fractionation affecting the 87 Sr/ 86 Sr ratio is corrected for online during mass spectrometric measurements. Thus, the 87 Sr/ 86 Sr ratio of primary minerals only depends on its original ratios of 87 Rb/ 86 Sr and 87 Sr/ 86 Sr, and geological age. Following chemical weathering and biological uptake of Sr, inorganic and organic matter carries the 87 Sr/ 86 Sr ratio of its source. Comprehensively, the 87 Sr/ 86 Sr ratio has become the proxy of choice in provenance studies aiming at fingerprinting the geographical origin of terrestrial material across scientific disciplines including anthropology, archaeology, ecology, geochemistry, and geology. For example, in the agricultural sector 87 Sr/ 86 Sr ratios were applied to trace the origin of beverages such as wine (e.g., Almeida and Vasconcelos 2001;Marchionni et al. 2013;Coldwell et al. 2022), orange juice , and bottled water (e.g., Montgomery et al. 2006;Voerkelius et al. 2010); foods such as rice (Kawasaki et al. 2002;Ariyama et al. 2012), wheat (Liu et al. 2016(Liu et al. , 2017, olives (Techer et al. 2017), onions (Hiraoka et al. 2016), and other vegetables (e.g., Swoboda et al. 2008;Song et al. 2014;Aoyama et al. 2017); and products with stimulating effects such as coffee beans (Rodrigues et al. 2011) and marijuana (West et al. 2009). Whereas in these studies the 87 Sr/ 86 Sr ratio of the products is directly matched with the 87 Sr/ 86 Sr ratio of the source as a whole, the radiogenic Sr isotope ratio is also suited to fingerprint the depth of Sr uptake across an individual soil depth profile. In doing so, the 87 Sr/ 86 Sr ratio of individual plant organs pinpoints the 87 Sr/ 86 Sr ratio of the bioavailable fraction of Sr in soil -that is Sr in soil water and/ or the exchangeable fraction of Sr in soil (Song et al. 2014) -as demonstrated for grass, shrubs and trees (McCulley et al. 2004;Poszwa et al. 2004;Bélanger and Holmden 2010;Coble et al. 2015;Aguzzoni et al. 2019;Uhlig et al. 2020).
The objective of this study was to evaluate field experiments in which subsoil management was performed by deep loosening in different ways, and to explore the following questions by means of the elemental composition and the radiogenic Sr isotope ratio of shoots and their soil resources: 1.) To which depth does deep loosening deepen the nutrient uptake depth? 2.) Does deep loosening trigger nutrient uplift? 3.) If so, how fast is nutrient uplift?

Field site
The field study was conducted at the Experimental Station of the University of Bonn "Campus Klein-Altendorf" located in Rheinbach (6° 59′ 29'' E, 50° 37′ 21'' N). The climate is temperate and humid with a long-term (1956 to 2020) mean annual temperature of 9.6 °C, mean annual precipitation of 603 mm, and a pronounced water deficit of about -640 mm throughout the summer season. The soil is classified as Haplic Luvisol (hypereutric, siltic) developed from loess (IUSS Working Group WRB. 2015), and is characterized by a silty clay loam texture with clay accumulation in the subsoil between 45 and 95 cm soil depth. The CaCO 3 content is < 1 g kg −1 from 0 -1.27 m soil depth and rises to 127 g kg −1 below 1.27 m depth (Barej et al. 2014). The bulk soil density increases from 1.29 g cm −3 to 1.52 g cm −3 from topsoil to 1 m soil depth (Barej et al. 2014). For a more detailed description of the soil the reader is referred to Barej et al. (2014).

Experimental field trial and subsoil management
Field experiments were conducted at two central field trials (CF), namely CF 1-1 and CF 2, being ~2.1 km apart from each other. Each central field trial was subdivided into single plots allowing various subsoil management methods complemented by a control with three field replicates at CF 1-1, and four field replicates at CF 2. A detailed description on the technical realization of subsoil management is given in Jakobs et al. (2019), Hinzmann et al. (2021), Schmittmann et al. (2021). In brief, subsoil management was carried out by strip-wise melioration in September of the respective starting year (CF 1-1: 2016, CF 2: 2019) and included i) deep loosening (DL), ii) deep loosening with the incorporation of biowaste compost (DLB) from kitchen waste of private households, and iii) deep loosening with biowaste compost in combination with the deep rooting pre-crop lucerne (DLB luc.). Quantities of incorporated organic material amounted to 5 kg m −2 at CF 1-1 and 3 kg m −2 at CF 2. Deep loosening was applied to a depth of 60 cm. At CF 1-1 organic material was manually incorporated into subsoil by strip wise removing topsoil, mixing compost with subsoil from the E/B and Bt1 horizon at 30 cm to 60 cm depth and re-adding topsoil. To allow strip wise subsoil management of a larger area, organic material was mechanically incorporated at CF 2 into the E/B horizon ranging from 30 to 45 cm. Visual inspections did not provide indications that there was any unwanted admixture of compost into the Ap horizon during compost incorporation into subsoil. Prior to the subsoil management both central field trials were limed with converter lime (Lhoist Germany, Rheinkalk GmbH) in 2015 (0.3 kg m −2 ) and 2016 (0.4 kg m −2 ) and annually fertilized with ~30 g m −2 calcium ammonium nitrate. Preceding cultivars were wheat and barley at CF 1-1 and grass at CF 2. In addition, only at CF 2 lucerne was cultivated as pre-crop in 2016 in half of the field plots, regularly harvested, and mulched in 2018. After subsoil management mustard was sown at both field trials as catch crop and mulched prior to performing a crop rotation in the order of spring barley followed by winter wheat. Winter wheat at CF 1-1 was sown on 23 rd October 2017 and spring barley at CF 2 was sown on 25 th March 2020. No irrigation was performed from sowing to sampling.

Sampling and sample preparation
Winter wheat (Triticum aestivum L.) and spring barley (Hordeum vulgare L.) were sampled at flowering stage on 30 th May 2018 at CF 1-1 and on 22 nd June 2020 at CF 2, respectively. In doing so, bunches of aboveground shoots were collected from each field replicate, each of which contained a minimum of ten shoots. The reason for omitting roots from sampling was that only the topmost part of the roots was removable from the hardened topsoil. Thus, root biomass, accounting for about a third of the crop's total biomass under non nutrient deficient conditions (Lopez et al. 2023), would have been sampled incompletely. Representative aliquots of biowaste compost were taken in 2016 at CF 1-1 only.
Plant and soil samples were transferred from field to laboratory on the day of sampling and frozen at − 20 °C. Prior to soil sieving (< 2 mm fraction), or dissection of plant samples into whole ear, stem and leaves, samples were dried by lyophilization for a minimum of 24 h at − 55 °C using a Christ Beta 1-8LD plus freeze drier (Martin Christ Gefriertrocknungsanlagen GmbH, Germany). Aliquots of compost and converter lime were dried in an oven for 24 h at 60 °C. Plant organs were milled in 100 ml sealable HDPE bottles equipped with tungsten carbide milling balls using a shaker (Collomix Agia 330, Collomix GmbH, Germany). To minimize analytical effort and costs particularly on the analyses of Sr isotope ratios the field replicates of soil samples were pooled and homogenized using a rotator overhead shaker (BioSan Multi Bio RS-24, Germany).

Analytical methods
Soil digestion using Li 2 B 4 O 7 fusion and the analyses of element concentrations using inductively coupled plasma optical emission spectroscopy (ICP-OES) was performed at the Central Institute of Engineering, Electronics and Analytics (ZEA-3) of Forschungszentrum Jülich GmbH. The remaining analytics comprising microwave assisted sample digestion, element concentration analyses by quadrupole inductively coupled plasma mass spectrometry (ICP-MS), analyses of loss on ignition, and Sr isotope ratios by multicollector inductively coupled plasma mass spectrometry (MC-ICP-MS) were performed at the Institute of Bio-and Geosciences (IBG-3) at Forschungszentrum Jülich GmbH that is equipped with a self-made, low particulate clean laboratory.

Sample digestion and concentration analyses
Plant and compost samples. A maximum of 400 mg of powdered and homogenized whole ears, stems, leaves, and compost were digested using a pressurized microwave digestion system (turboWAVE, Milestone Srl, Italy) and a mixture comprising 3.5 ml 15 M HNO 3 , 2 ml 30% H 2 O 2 and 2.5 ml deionized water (18.2 MΩ cm, TOC < 3 ng g −1 , Merck Millipore, Germany). Plant digests were centrifuged at 8000 rpm for 5 min (Beckman-Coulter Allegra X30R, USA, with a rotor of fixed-angle (25°) and a maximum radius of 97 mm) to separate the dissolved sample from silicate residues and the supernatant was pipetted off. The residuum was rinsed three times each with 1 ml deionized water, centrifuged and the supernatants were combined to ensure the complete transfer of the dissolved sample. Element concentrations were analysed by ICP-OES (iCAP 6500, Thermo Scientific, Germany) following the Cs-HNO 3 method described in Schuessler et al. (2016). Relative uncertainties, given as the relative deviation of the measured element concentration of a standard reference material (SRM) to its certified value, were assessed by processing NIST SRM 1515 apple leaves with each sample digestion batch and by processing NIST SRM 1573a tomato leaves, NIST SRM 1575a pine needles and ERM-CD281 rye grass occasionally. Relative uncertainties were reported in Table S5 (Uhlig 2022) and typically amounted to better than ± 10%. Element concentrations of shoots ( [X] shoot ) were calculated from Eq. 1, where j refers to the crop organs ear, stem, or leaf, and f denotes the relative proportion of the dry biomass of the respective crop organ j to the total dry biomass of the shoot.
Bulk soil and converter lime samples. Prior to digestion a representative aliquot of 1 g of pooled and homogenized soil was milled using an agate mortar and pestle. Then, fusion digestion was performed by mixing 50 mg of soil with 250 mg of Li 2 B 4 O 7 followed by heating at 1050 °C for 30 min in a muffle furnace. Element concentrations were analysed by ICP-OES (iCAP 6500, Thermo Scientific, Germany) and concentrations of major elements were expressed as oxide concentrations. Loss on ignition (LOI) was determined by heating 0.5 g of soil at 550 °C for 2 h in a muffle furnace to rescale oxide concentrations to 100%. An aliquot of 50 mg of converter lime was dissolved in 6 M HCl, evaporated to dryness and re-dissolved in 0.3 M HNO 3 . Element concentrations of converter lime were analysed by ICP-MS (Agilent 7900, Agilent Technologies Inc., USA). Relative uncertainties of the fusion digestion were assessed by repeat digestion of NIST SRM 2709a San Joaquin soil, reported in Table S1b (Uhlig 2022), and typically amounted to better than ± 5%.

Exchangeable fraction in soil
The exchangeable fraction in soil was extracted by suspending two grams of dried and sieved soil (< 2 mm fraction) in 20 ml of 1 M NH 4 OAc and gently shaking this suspension for 2 h with a rotator overhead shaker (BioSan Multi Bio RS-24, Germany). The suspension was centrifuged at 8000 rpm for 5 min, and the supernatant was passed through a 0.20 µm acetate filter. Element concentrations were measured by ICP-MS (Agilent 7900, Agilent Technologies Inc., USA). Relative uncertainties were assessed by repeat analyses of NIST SRM 1640a natural water that is appropriate to resemble the sample matrix of the exchangeable fraction in soil, reported in Table S4 (Uhlig 2022), and typically amounted to better than ± 5% (Table S4).

Radiogenic Sr isotopes
Strontium was purified from matrix elements using converted disposable pipettes packed with 200 µl TrisKem SR-B50-S resin (50-100 µm). A sample aliquot containing 400 ng Sr was loaded onto the resin. Matrix elements were removed by elution with 5.5 ml 7.5 M HNO 3 and Sr was eluted with 2 ml deionized water. After evaporation to dryness, organic crown-ethers released from the Sr resin were oxidized in a 1:1 (v/v) mixture comprising 30% H 2 O 2 and concentrated HNO 3 , heated in closed PFA vials first for 1 h at 85° C and then for 3 h at 170 °C. After evaporation to dryness and redissolution in 0.3 M HNO 3 the purity of the Sr fraction was checked by ICP-MS analysis.
The 87 Sr/ 86 Sr ratio was measured with a MC-ICP-MS (Nu Plasma II, Nu Instruments Ltd, UK) equipped with an APEX-IR (Elemental Scientific Inc., USA) at low mass resolution in 1 block of 100 cycles with 10 s integration time. The intensity of the blank yielded ≤ 20 mV, which was negligible compared with the intensity of the sample yielding about 20 V on mass 88. The signals of the masses 82, 83 and 85 were measured together with the Sr masses (84,86,87,88) to correct for Kr and Rb interferences on mass 84 and 87, respectively. The measured 87 Sr/ 86 Sr ratio was normalized to the 86 Sr/ 88 Sr ratio of 0.1194 (Steiger and Jager 1977) using an exponential law to correct for natural and instrumental massdependent isotope fractionation. Repeated analyses of SRM 987 with and without Sr purification were used to determine the long-term accuracy of the method yielding mean values and two standard deviations of 0.71029 ± 0.00004 (N = 6, with Sr purification) and 0.71028 ± 0.00006 (N = 67, without Sr purification). Within two standard deviations the measured 87 Sr/ 86 Sr ratios of SRM 987 agree with its reference ratio (0.71025 ± 0.00005, level of confidence 95%, N = 474) reported in the GeoReM database (Jochum et al. (2005), http:// georem. mpch-mainz. gwdg. de) for MC-ICP-MS analyses.

Statistical analysis
Null hypothesis significance testing was carried out on the difference among the control and subsoil managed plots in the variables i) biomass, ii) element concentrations and iii) Sr isotope ratios using the data analysis software OriginPro 2022 (version 9.9.0.225; OriginLab Corporation; Northampton, MA, USA). Because field replicates of soil samples were pooled, statistical analyses were limited to field replicates of shoots. Prior to significance testing, datasets were checked for normal distribution (Shapiro-Wilk, p < 0.05) and variance equality (two sample F-test, p < 0.05). If both apply, the significance of difference was analysed using a two-sample t-test (p < 0.05). In case of variance inhomogeneity, a Welch-correction was performed. In case of non-normally distributed datasets, a non-parametric test was performed (Wilcoxon-Mann-Whitney test, p < 0.05).

Results
The elemental composition and the radiogenic Sr isotope ratio of individual crop organs, compost, converter lime and pooled soil samples including bulk soil (Li 2 B 4 O 7 digestion) and the exchangeable fraction in soil (1 M NH 4 OAc) is reported in Table S1-S5 (Uhlig 2022). Biomass and element concentrations of shoots averaged from individual field replicates are summarized in Table 1. Concentrations of selected elements, namely macronutrients (K, Ca, Mg, P), micronutrients (Fe, Mn, Zn), a plant beneficial element (Na), and non-nutritive elements (Ba, Sr) are displayed for the CF 1-1 trial in Fig. 1 and for the CF 2 trial in Fig. 2 for the compartments shoots (panels: a-e, p-t), the exchangeable fraction in soil (panels: f-j, u-y), and bulk soil (panels: k-o, z-ad). The 87 Sr/ 86 Sr ratios are shown in Fig. 3 for both trials.
Shoot biomass (Table 1) and bulk soil density (not shown, personal communication, Julien Guigue, September 2022) among control and managed plots were identical within uncertainty. Therefore, element concentrations rather than element stocks were considered hereafter. Overall, element concentrations of shoots, the exchangeable fraction in soil, and bulk soil were identical among control and deep loosening (DL) plots. Consequently, no management effects could be found with DL. In contrast, element concentrations and the radiogenic Sr isotope ratio of mainly the exchangeable fraction in soil differed from the control plots, when DL was combined with the incorporation of biowaste compost (DLB) and when lucerne (luc.) was cultivated as deep rooting precrop. Consequently, management effects were found with DLB and DLB luc., which are detailed next.

Elemental composition of shoots and soil compartments
Shoots barely responded to subsoil management with biowaste compost (DLB). Exceptions towards increased element concentrations at DLB plots of both trials were found for the plant-beneficial element Na. Its magnitude was much more pronounced than for micronutrients like Zn at CF 1-1, and Fe and Mn at CF 2. However, pre-cropping with the deep rooting lucerne resulted in an increase in element concentrations of Na, Mn, and P in shoots cultivated on the control plots. Beyond that, concentrations of the macronutrients Ca and Mg (both CF trials) and concentrations of the micronutrients Fe and Mn (CF 2) were lower than required for adequate growth (cf. data points with vertical dashed green lines in Fig. 1 and Fig. 2). Yet, the possibility of growth limitation was discounted because shoot biomass remained identical among the control and DLB plots.
The exchangeable fraction in soil showed management effects at both CF trials. At CF 1-1, concentrations of K, Ca and Na increased at the depth of compost incorporation and above (K, Ca, Mg, P) or below (Na). At CF 2, all elements responded to subsoil management, indicated by increased concentrations of K, Ca, Mg, P, Na, Fe, Mn and Sr, and decreased concentrations of Zn and Ba at the depth of compost incorporation.
Bulk soil revealed management effects towards increased concentrations of Ca, P, Zn, and Sr at the depth of compost incorporation at CF 2. These effects were consistent with the admixture of biowaste compost being three to five times higher in concentration of these elements than bulk soil. The relative contribution of elements in the exchangeable fraction to bulk soil was on average about 31% (Ca), 6.1% (Sr), 4.7% (Ba), 2.6% (Mg), < 1% (remaining elements). Thus, apart from Ca only a minor fraction of monoand bivalent elements was bioavailable in soil.
In the exchangeable fraction, the 87 Sr/ 86 Sr ratio ranged from 0.70896 to 0.71269 and featured four characteristics. First, 87 Sr/ 86 Sr increased from topsoil to the base of the sampled subsoil, suggesting the presence of distinct Sr sources in the exchangeable  Table S1b and Table S3 (Uhlig 2022). Vertical brown bars next to the most righthand panels illustrate soil horizons according to IUSS Working Group WRB. (2015). Capital letter A indicates significant (p < 0.05) concentration differences in shoots among management and the control fraction with depth. Second, the 87 Sr/ 86 Sr ratios of all topsoil samples were shifted from more radiogenic 87 Sr/ 86 Sr ratios towards the less radiogenic 87 Sr/ 86 Sr ratio of converter lime (0.70913), which is roughly identical with modern-day seawater (0.70918 (Veizer and Compston 1974;Palmer and Edmond 1989;Allègre et al. 2010)). Thus, the dissolution of lime and the subsequent entry of Sr and other elements  Table S1b and Table S3 (Uhlig 2022). Vertical brown bars next to the most right-hand panels illustrate soil horizons according to IUSS Working Group WRB. (2015) and white horizontal dashed lines subdivide the Ap horizon at CF 2 into upper, middle and lower Ap horizon. Capital letters refer to significant (p < 0.05) concentration differences in shoots among the control and managed field trial. Specifically, A refers to C -DLB, B refers to C luc. -DLB luc., C refers to C -DLB luc., D refers to C -C luc contained in converter lime (e.g., B, Ca, Co, Cu, Cr, Fe, Mg, Mn, Mo, P, Si, Zn) into the exchangeable fraction in soil exerts a major control on the elemental and isotope composition of the exchangeable fraction in topsoil. This finding is consistent with the fact that the application of agricultural lime on non-calcareous soils leads to a change in the elemental and isotope composition of the soil as demonstrated by Thomsen and Andreasen (2019). Third, at CF 2 the 87 Sr/ 86 Sr ratio was invariant in the lower two Ap horizons, which reflects soil homogenization through ploughing activities throughout past decades. Fourth, consistent with concentration measurements, management effects were evident at CF 2 by two lines of evidence. On the one hand, in both controls the 87 Sr/ 86 Sr ratios of the exchangeable fraction in the E/B horizon right below the ploughing depth, ranging from 0 -30 cm, differed from the 87 Sr/ 86 Sr ratios in the ploughing horizon. On the other hand, the invariance in 87 Sr/ 86 Sr ratios continued to the E/B horizon at the DLB managed plots. This finding suggests an additional and less radiogenic Sr source in the exchangeable fraction, namely organic-derived Sr, that cannot be found in the bulk digest of biowaste compost. Given that biowaste compost is sourced from multiple kitchen households, Sr is likely present in forms with different solubility. Thus, future studies should additionally focus on the bioavailable form of Sr in compost, which is speculated to be less radiogenic than its bulk digest.
A match of 87 Sr/ 86 Sr ratios among shoot organs and the exchangeable fraction occurred in different  Table S4 (Uhlig 2022)) at the central field trials CF 1-1 (top panels) and CF 2 (bottom panels) cultivated with winter wheat and spring barley under different subsoil managements and the control. Vertical dashed lines bracket the 87 Sr/ 86 Sr ratios of shoot tissue and filled vertical bars denote a match of the 87 Sr/ 86 Sr ratio among shoot tissue and the exchangeable fraction in soil in the respective soil horizon. Vertical brown bars next to the most right-hand panels illustrate soil horizons according to (USS Working Group WRB. (2015) and white horizontal dashed lines subdivide the Ap horizon at CF 2 into upper, middle and lower Ap horizon. Error bars of shoots from CF 2 denote the two-fold standard deviation of four field replicates. Error bars of the remaining samples are smaller than the symbol size and refer to the analytical uncertainty (two-fold standard error) soil horizons in the control and managed plots. For example, at CF 1-1 an isotope ratio match among ears and the exchangeable fraction occurred in the Bt1 horizon in the control plot, from the Ap to the Bt1 horizon in the DL plot, and from the Bt1 to the Bt3 horizon in the DLB plot. In contrast, at CF 2 an isotope ratio match among shoot organs and the exchangeable fraction occurred from the middle Ap to the E/B horizon in the control and in the control in combination with the deep-rooting pre-crop lucerne. Moreover, an isotope ratio match among crop organs and the exchangeable fraction extended to the Bt1 horizon in the DLB and DLB luc. plots at CF 2.

Discussion
Nutrient uptake depth inferred from 87 Sr/ 86 Sr ratios The 87 Sr/ 86 Sr ratio of crop organs serves as a fingerprint of the 87 Sr/ 86 Sr ratio of the exchangeable fraction in soil at the depth of Sr uptake. This depth is interpreted as the mean depth from which water and nutritive elements in bioavailable form are sourced. It also encompasses water and nutrient uptake from above this depth, for example, after germination during the crop's early growing stages, and below this depth, for example, at the time of sampling at which the present-day uptake depth is not yet well recorded in the 87 Sr/ 86 Sr ratio of shoots due to reservoir effects. Four major implications can be drawn from the results.
First, in the control plots the mean nutrient uptake depth of winter wheat lay in the Bt1 horizon, thus ranged from 50 to 60 cm depth (Fig. 3a), whereas the mean nutrient uptake depth of spring barley ranged from the top of the middle Ap horizon at 20 cm depth to the bottom of the E/B horizon at 40 cm depth (Fig. 3k). Consequently, winter wheat took up water and nutrients from greater depth than spring barley. This finding is consistent with the fact that winter wheat spent five months longer on the field upon the day of sampling than spring barley, allowing its roots to reach deeper soil regions than spring barley. According to Thorup-Kristensen et al. (2009) a cultivation sowed in winter (e.g., winter wheat) indeed grows two times deeper into the soil than a cultivation sowed in spring (e.g., spring wheat).
Second, subsoiling by deep loosening alone resulted into a mean nutrient uptake depth, which involved shallower soil horizons than in the control.
Specifically, the mean nutrient uptake depth in the DL plots ranged from the topmost Ap horizon at 0 cm depth down to the bottom of the Bt1 horizon at 60 cm depth (Fig. 3b). In contrast, the mean nutrient uptake depth in the control was limited to the Bt1 horizon, thus ranged from 50 to 60 cm depth (Fig. 3a). This finding was surprising because deep loosening was supposed to reduce soil penetration resistance to crop roots to enable roots for the access of deeper soil horizons for crop nutrition than in the non-loosened control. Thus, because deep loosening did not cause a shift of the mean nutrient uptake depth towards deeper soil horizons than in the control, root restricting layers apparently do not occur in the Luvisol at CF 1-1. Instead, deep loosening had the opposite of the supposed effect, namely an intensification of nutrient uptake in shallower soil horizons than in the control, which even hindered roots from foraging into deeper non-loosened soil for water and nutrient uptake.
Third, cultivating lucerne as deep rooting precrop resulted in the same mean nutrient uptake depth of the succeeding crops than without precropping with lucerne, which ranged from the top of the middle Ap horizon at 20 cm depth down to the bottom of the E/B horizon at 45 cm depth (Fig. 3l). On the one hand, this finding conflicts with the underlying hypothesis that lucerne creates deep biopores (Han et al. 2015), which allow roots of succeeding crops to develop a deeper rooting system (Han et al. 2017) than without pre-cropping with lucerne. On the other hand, this finding demonstrates that, although lucerne can significantly dry out the subsoil by transpiration, pre-cropping with lucerne had no negative effect on the mean nutrient uptake depth, even when the re-wetting period faces a severe drought with a water deficit amounting to about -640 mm.
Fourth, deep loosening in combination with the incorporation of biowaste compost into subsoil resulted not only in a deepening of the maximum rooting depth (Bauke et al., n.d.) but also in a deepening of the mean nutrient uptake depth (Fig. 3f, m-n), which extended to below the depth of compost incorporation (Fig. 3f, m-n). This deepening can likely be attributed to an improved water holding capacity of the organic matter added to subsoil incorporated into subsoil (Lal 2020).
Shifts in the mean nutrient uptake depth towards deeper soil horizons through DLB and DLB luc. were also accompanied by increases in element concentrations in shoots. The most affected elements were the redox-sensitive micronutrients Mn (DLB and DLB luc. t CF 2) and Fe (DLB luc. at CF 2), and the plant beneficial element Na (DLB at both trials, and DLB luc. at CF 2). The release of Mn and Fe might have been promoted by increased microbial utilization of the compost substrate, and thus lower redox potential in addition to possible releases of chelating agents. In contrast, the concentration increases of Na cannot solely be attributed to increased abiotic weathering because the soil pH at CF 2 increased from 6.9 in the control plot to 7.5 after compost incorporation into subsoil (Uhlig 2022). Even though some minor portions of nutrients might have been mobilized from subsoil by enhanced chemical weathering caused by microbial priming (Fang et al. 2023), this contribution is likely negligibly small compared with the supply of these elements in bioavailable form with biowaste compost. We thus attributed the increased element concentrations in shoots to peaks in the depth of compost incorporation (DLB and DLB luc. plots). This raises the question, whether subsoiling with compost results not only in nutrient uptake followed by harvest, but also in nutrient uplift, hence in the entry of these elements into the exchangeable and extractable fraction in topsoil as recently speculated by means of Mg stable isotopes at the same field site (Uhlig et al. 2022). A manifestation of which presents a sustainable side-effect of subsoiling with biowaste compost because the large geogenic reservoir of bioavailable nutrients situated below the depth of compost incorporation becomes an active component of the nutrient cycle in the soil-plant system that would otherwise remain unused by plants.

Nutrient uplift
The 'nutrient uplift hypothesis' from Jobbágy and Jackson (2001), formerly known as the 'base pump' (Vejre and Hoppe 1998), implies that mineral nutrients are taken up from the mineral soil, cycled through biota, returned via three potential pathways to topsoil to finally end up in the exchangeable and extractable fraction in topsoil. The three potential pathways include i) throughfall; ii) solubilised litterfall, harvest residues, and necrotic biomass; and iii) the so-called sink-regulated transport that is long-distance transport of nutrients through the phloem and its extrusion by roots back into soil, which is most pronounced for the phloem-mobile macronutrients K, Mg and P (White 2012) and which serves plants for biosensing nutrient availability in the growth medium (Ho and Tsay 2010;Krouk et al. 2010). No matter which pathway applies, soil profiles, strongly affected by nutrient uplift, are characterized by declining concentrations of mineral nutrients in the exchangeable and extractable fraction with depth. The magnitude of this decline is most pronounced for mineral nutrients that often limit plant growth like K and P because a high demand in plants is opposed to a low supply in bioavailable form in the soil (e.g., in the exchangeable and extractable fraction). The timescale over which nutrient uplift occurs was found hitherto to range from years (≤ 3 yr) after cultivating deep rooting precrops in cropping systems (Han et al. 2021) to more than a decade (≤ 15 yr) after afforesting a sand dune . However, the depth from which nutrient uplift occurs and the kinetics of nutrient uplift in cropping systems remain poorly constrained. Hence, Han et al. (2021) suggested the performance of tracer experiments to better constrain nutrient uplift in cropping systems. Given that the incorporation of biowaste compost into subsoil added nutrients in bioavailable form to the subsoil (cf. peaking concentrations in the exchangeable fraction at the depth of compost incorporation in Fig. 1f, h, j, Fig. 2f-h) the subsoil experiment performed in the present study can be considered as such a tracer experiment that operates at nutrient concentrations at natural abundance levels. It also circumvents the disadvantage that the 87 Sr/ 86 Sr ratio is unsuited to trace nutrient return to topsoil in the present study because the 87 Sr/ 86 Sr ratio in the exchangeable fraction of topsoil is dominated by converter-lime-derived Sr.
Along these lines, nutrient uplift was only found at CF 1-1. The magnitude of nutrient uplift, defined here as the proportional concentration increase in the exchangeable fraction in topsoil from C to DLB plots, was most pronounced for the macronutrients P (175%), K (120%), Mg (55%), and the micronutrient Mn (25%). Thus, the exchangeable fraction in topsoil indeed features a biological imprint, or in other words ecological stoichiometry. Provided that crop nutrition occurred under nutrient sufficient conditions, this form of biological pumping also explains the lack of increased concentrations of K, Mg, P, and Mn (only at CF 1-1) in shoots at the DLB plots. However, nutrient uplift was not evident at CF 2. The most likely reason is the six-months shorter cultivation period of crops at CF 2 because, in contrast to CF 1-1, winter wheat is not yet integrated into the crop rotation at CF 2. Thus, the interaction time of crops with the soil was longer at CF 1-1 allowing pathway iii) for causing nutrient uplift. Overall, nutrient uplift was evident as early as 21 months after subsoiling with biowaste compost at CF 1-1.
Role of Na in shifting the nutrient uptake depth below the depth of compost incorporation Although Na showed the most striking management effects towards increased concentrations in shoots and in the exchangeable fraction at the depth of compost incorporation, an uplift of Na into the exchangeable fraction of topsoil was not evident. Instead, Na concentrations peaked in the exchangeable fraction in the soil horizon below the depth of compost incorporation at both CF trials. This concentration peak is attributed to the solubilisation of compost-derived Na followed by its downward movement by gravitational water flow and adsorption at deeper soil horizons. Sodium, as a chemical twin of K in plants (Benito et al. 2014) can replace K (e.g., Subbarao et al. 2003;Wakeel et al. 2010) to act as an additional osmoticum to regulate the movement of stomata, which was found to be faster when Na was involved (Wakeel et al. 2011). Thus, Na reallocated towards deeper soil horizons could have been traced by roots biosensing for nutrient availability and mined to improve the water use efficiency of crops, thereby to cope with the water deficient years that the CF trials experienced. Hence, reallocated Na may have triggered the deepening of the mean nutrient uptake depth to below the depth of compost incorporation. Consequently, subsoiling with compost incorporation may not only improve the water use efficiency of crops but also shift the mean uptake depth of nutrients situated beyond the depth of compost incorporation as shown by means of 87 Sr/ 86 Sr ratios.

Conclusions
The performance of subsoil management by deep loosening or using lucerne as deep rooting pre-crop, and both of which in combination with compost incorporation was evaluated with respect to the nutrient budget in shoots, the exchangeable fraction in soil and in bulk soil, and the nutrient uptake depth by means of the elemental composition and the radiogenic Sr isotope ratio ( 87 Sr/ 86 Sr). Regardless of the applied subsoil management practice, negative effects on the nutrient budget in shoots and soil could not be found. On the contrary, Na contents in shoots consistently increased when subsoiling included compost incorporation. Moreover, using the 87 Sr/ 86 Sr ratio as a proxy for fingerprinting the mean nutrient uptake depth of crops revealed two key findings. First, subsoiling by deep loosening alone resulted into a more intense involvement of shallower soil horizons in nutrient uptake than in the control. Second, deep loosening with compost incorporation resulted into a deepening of the mean nutrient uptake depth below the depth of compost incorporation. This deepening could be triggered by biosensing roots, which not only traced increased levels of bioavailable Na below the depth of compost incorporation, but also utilized this compost-derived Na to improve the water use efficiency of crops during a summer drought. Finally, in settings, where agricultural lime is applied, the 87 Sr/ 86 Sr ratio is unsuited to trace nutrient uplift. However, an assessment of the elemental composition of the exchangeable fraction suggested nutrient uplift to occur within 21 months after subsoiling from the depth of compost incorporation into topsoil. Overall, subsoil management by deep loosening with compost incorporation causes a deepening of the mean nutrient uptake depth and nutrient uplift from depth to topsoil, thereby serves to be a more promising soil management practice to utilize soil resources beyond the hitherto intensively managed topsoil in the future.
Author contribution All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by David Uhlig. The first draft of the manuscript was written by David Uhlig and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.
Funding Open Access funding enabled and organized by Projekt DEAL. This work was supported by the German Federal Ministry of Education and Research (BMBF) for funding the project "BonaRes (Module A): Sustainable Subsoil Management-Soil3-II; subproject C" (grant number 031B0515C).
Data availability Research Data associated with this article can be accessed at GFZ Data Services under the reference Uhlig (2022). Tables S1-S5 include the dataset discussed in this publication along with further background data.

Declarations
Competing interest The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.