Development of mechanical soil stability in an initial homogeneous loam and sand planted with two maize (Zea mays L.) genotypes with contrasting root hair attributes under in-situ field conditions

Soil structure evolving from physical and biological processes is closely related to soil mechanical characteristics and texture. We studied the influence of substrate and genotype on the initial development of mechanical traits, differences between depths, and changes over the course of two years in the field. Plots were homogeneously filled with a loam and a sand and planted with two maize (Zea mays L.) genotypes (wild type (WT) and rth3 mutant) with contrasting root hair attributes. Undisturbed soil cores were taken in 2019 and 2020 at 14 and 34 cm depth. Confined uniaxial compression tests were performed to determine pre-compression stress (σpc), compressibility (Cc, Cs) and elasticity index (EI). Mechanical energy was calculated based on penetration resistance (PR) tests with a penetrometer needle resembling root geometries. σpc, Cc and Cs were significantly higher in loam as compared to sand, whereas the factor genotype proved to be negligible. Over time, σpc increased and Cc decreased in loam from 2019 to 2020 and Cs declined in both substrates. Higher mechanical energies were observed in loam and partially in WT. Required energy was higher at 14 cm than at 34 cm depth and decreased from 2019 to 2020 in sand. Airdry sand samples required four times as much energy than those at matric potential (Ψm) of -50 kPa. For the development of the mechanical traits examined texture proved to be the dominating factor and changes in soil stability could be observed within a short period of time.


INTRODUCTION
Soil mechanical characteristics evolve from the interplay of texture and structure of a soil, structure itself being influenced by a variety of biological, physical, and chemical factors (Bacq-Labreuil et al. 2019;Bryk et al. 2017;Grosbellet et al. 2011;Hallett et al. 2009;Naveed et al. 2014;Schon et al. 2017). Root growth constitutes a major biological factor for the formation of soil structure and depends on the plant's characteristics and the prevailing environmental conditions. Structural changes take place over time and are accompanied by the evolution of the mechanical properties of the soil. This development can be expected to be considerably different in the field compared to controlled lab conditions.
How root growth re-organises the pore system and aggregate size distribution, largely depends on soil texture. Bacq-Labreuil et al. (2019) showed that phacelia (Phacelia tanacetifolia Benth.) decreased soil porosity in a sandy loam as compared to the unplanted control but did not alter porosity of a clayey soil. The changes in bulk density (ρ b ) induced by root growth are mainly determined by soil structure and particle size. When roots grow in homogeneous soil, particles are rearranged, whilst in a more structured soil the existing pore space is utilized . Growth of coarse roots increase macropores, whereas finer roots display a high flexibility and improve micro-porosity (Bodner et al. 2014). Roots do not only greatly contribute to the creation of biopores but also improve connectivity of the pore system (Lucas et al. 2021). Furthermore, they have an influence on soil aggregation. In an agriculturally reclaimed mining area, macro-aggregate size distribution remained similar throughout the 24 years reclamation age (Pihlap et al. 2019), whereas aggregates <1000 μm increased in a clay through phacelia growth (Bacq-Labreuil et al. 2019).
Mucilage released at the root tip (Oleghe et al. 2017;Naveed et al. 2017) and an intact root tip (Iijima et al. 2004;Bengough and McKenzie 1997) ease the roots' way through the soil. The rhizosphere provides a habitat for soil microorganisms and root exudates have a major impact on bacterial community structure in the rhizosphere (Haichar et al. 2008). Roots contribute indirectly to soil structural evolution even after the plant's death, by providing nutrient and carbon to soil organisms (Li et al. 2015;Pett-Ridge and Firestone 2017), increasing both functional and taxonomic diversity though root detritus (Nuccio et al. 2020) and altering microbial activity through the spatial organisation of the pore network (Nunan et al. 2017). Furthermore, former root channels constitute an important root growth pathway in soil layers with a high mechanical impedance (Han et al. 2015). Both root exudates (Burak et al. 2021;Galloway et al. 2020) and existence of root hairs (Carminati et al. 2017) contribute to the formation of a rhizosheath, which in turn can stabilize soil (Wang et al. 2017) and affect both porosity and connectivity of the surrounding pores . In the light of these processes plant growth affects soil structure formation and stabilization by a wide range of interacting mechanisms, which in turn ultimately influences soil mechanical parameters. Roots may encounter layers which they are not able to penetrate, unfavourable mechanical conditions being a possible reason for that. For instance, Bengough et al. (2011) mentioned penetration resistances (PR) of 2 MPa reducing root elongation rates in maize by 50% compared to nonimpeded roots.
It is a combination of both abiotic and biotic factors that influences structure formation in a soil (Barto et al. 2010), oftentimes it is difficult to distinguish between them (Oades 1993). Physical factors such as climatic conditions and weather, influence soil structural development directly by freezing and thawing or swelling and shrinking cycles and indirectly through affecting plant growth. Root systems with finer roots can dry out the soil locally leading to cracks (Oades 1993), whereas larger roots can create macropores or use existing macropores to reach deeper soil regions  in search of water and nutrients. Leuther and Schlüter (2021) examined the impact of freezing and thawing cycles on soil structure and found that even only two of these cycles might already lead to a fragmentation of soil clods and an increase of unsaturated hydraulic conductivity. Wetting and drying leads to swelling and shrinking processes in the soil which affect macro-and mesoporosity and greatly depend on clay content and 1 3 Vol.: (0123456789) mineralogy (Diel et al. 2019). Stronger organic bonds leading to an increase in macroaggregates and their stability, can be expected between 2:1 clay minerals compared to a soil in which 1:1 clay minerals dominate (Denef and Six 2005). In managed grasslands Barto et al. (2010) showed that a higher percentage of sand in a soil interferes with soil aggregation.
The experimental site we examined consisted of plots that were homogeneously filled with two defined substrates and planted with two maize genotypes ). This set-up gave us the unique opportunity to study the development of soil mechanical properties following an initial homogeneous state under field conditions. The aim of this study was to quantify the effect of the factors substrate and genotype on soil stability, compressibility, elasticity, and mechanical energy required for root growth and to describe their spatial and temporal development.
The two genotypes employed in the experiment differ in their capability to form root hairs. The wild type (WT) forms root hairs, whereas in the rth3 mutant their elongation is suppressed. Root hairs play a central role in the cohesion between roots and soil particles (De Baets et al. 2020) and might therefore lead to higher stabilisation. However, differences in root length densities between genotypes might have a greater effect on soil mechanical properties than the actual root hairs. The samples in this study were taken between rows of maize due to spatial limitations. We hypothesized substrate to be the driving factor with higher stability and compressibility expected in loam as the higher clay content promotes soil aggregation thus enhancing soil stability. Bronick and Lal (2005) describe the important role of clay in the formation of compound particles that contribute greatly to soil aggregation. Shrinking and swelling processes in soils containing clay are crucial in soil structure formation (Dixon 1991). Higher clay contents and a greater level of aggregation in loam bring about a higher initial void ratio (ε) and its more pronounced decline when stresses exceed pre-compression stress (σ pc ) (Lebert and Horn 1991). The sand on the other hand is a single-grain substrate with weak cohesion between particles (Lebert and Horn 1991), less potential for aggregation and a pronounced primary pore system that leaves little room for compression.
Concerning the development over time and space, we expect σ pc to be higher closer to the soil surface, caused by the influence of environmental conditions, mainly drying and re-wetting. As structural development takes time, σ pc as a measure of stability is likely to increase from one year to the next and compressibility (compression index C c and swelling index C s ) is expected to decrease. Regarding PR, we predict higher values in loam compared to sand and a possible influence of genotype relating to their root system. We further hypothesize the heterogeneity to be greater at 14 cm depth and to increase with time due to the development of soil structure, all in all leading to higher values in the second year. Our overall hypothesis is, that in sand changes in mechanical characteristics caused by processes of structural formation will be less pronounced than in a loam.

Soil plot experiment
The soil plot experiment was established in October and November 2018 in the framework of the DFG priority programme 2089 "Rhizosphere Spatiotemporal Organisation -a key to rhizosphere functions". It constitutes a new approach of studying the development of soil mechanical stability parameters. The site is located at the experimental station of the Helmholtz Centre for Environmental Research in Bad Lauchstädt in Germany (N 51° 23.425440, E 11° 52.555980). Field plots were excavated and filled with two homogeneous substrates, a loam and a sand (loam: 32.5% sand, 47.9% silt, 19.5% clay; sand: 91.8% sand, 5.6% silt, 2.6% clay ), and planted with maize genotypes WT and rth3. The randomized block design consisted of six plots for each combination of texture and genotype, resulting in six field replicates. Fertilization ensured an equal supply with nutrients in both substrates. Maize was first sown in April 2019. Details regarding the experimental set-up including filling of field plots, fertilization, and agronomic measurements are described in Vetterlein et al. (2021).

Sampling and sample preparation
Two sets of undisturbed soil cores were taken at growth stage BBCH83 (early dough) (Meier and Biologische Bundesanstalt für Land-und Forstwirtschaft 1997) in 2019 and 2020, one for a confined uniaxial compression test, the other for PR measurements. For the compression tests, cylinders with a diameter of 10 cm and a height of 3 cm were employed, resulting in a ratio of diameter to height of 3.33, which is in accordance with the International Organization for Standardization (ISO 2017), whereas for the remainder of the experiments cylinders with a diameter of 10 cm and a height of 6 cm were utilized, giving enough room for PR measurements. Samples were taken at a specific location in each field plot at two different depths (14 cm and 34 cm upper cylinder edge) and stored in the dark at 4°C until further use. The sample location was between rows of maize with a distance of 22 cm to the next plant. The position was in close proximity to the location used for root sampling in Vetterlein et al. (2022) who took soil cores 10 cm from the foot of the maize plant. In preparation for the measurements, samples were saturated with tap water and subsequently put on ceramic suction plates until weight consistency was reached to obtain the desired matric potentials (Ψ m ) at the start of the measurements. These were -50 kPa Ψ m for the uniaxial compression test and a sequence of -3, -12.5 and -50 kPa Ψ m for the PR with an additional airdry step for sand in 2020. Note that the unit kPa is used both for the stress and for Ψ m , for which we explicitly stated "kPa Ψ m " to avoid confusion.

Confined uniaxial compression test
Undisturbed soil samples were measured at Ψ m -50 kPa with an oedometer ( fig. 1; Eijkelcamp 08.67 Compression test apparatus, Giesbeek, Netherlands) and the software Physical soil test Version 2.0.4 (Eijkelcamp, Giesbeek, Netherlands) to obtain information on pre-compression stress (σ pc ), compression indices (C c ), swelling indices (C s ), and elasticity indices (EI). The basic idea of the oedometer tests is to load a soil sample confined in a cylinder with a defined pressure and measure the settlement in vertical direction. Above and below the built-in soil sample are sintered metallic plates through which water can flow that is pressed out of the pores. The oedometer is furthermore equipped with a tensiometer to allow measurement of Ψ m throughout the experiment. In order to level out any irregularities at the sample surface, the device exerted a pre-pressure of 5 kPa for 20 s, which we defined as the actual start of the test. This was followed by 11 log-equidistant loading steps ranging from 10 to 575 kPa and a final unloading step of 5 kPa. The duration of each step was 20 min for loam and 5 min for sand to account for different consolidation times of these substrates and was established in a pre-test prior to the experiment. Applied stress, vertical settlement and Ψ m were logged every 2 s. For the necessary calculations bulk density (ρ b ) was determined by drying the soil cores at 105°C for 48 h after measurements; particle density (loam 2.46 g cm -3 and sand 2.62 g cm -3 , Rosskopf et al. 2022) had been quantified at an earlier stage.
A typical stress-strain curve (Fig. 1b) shows the settlement behaviour of a soil subjected to a series of successively increasing vertical loads and generally consists of a recompression range and a virgin compression range. Vertical stress is registered in kPa, where F is the vertical force applied and A the crosssectional area of the soil sample. The settlement is reflected in the void ratio ε, namely: The first part of the curve is the so-called recompression range, as the soil already had been exposed to stresses in this order of magnitude in the past. Its slope is defined as C s : Evaluation of σ pc was based on the methods of Casagrande (1936) and the logistic function suggested by Gregory et al. (2006). The slope of the virgin compression range (Fig. 1b) is called C c , which was calculated analogously to C s : The expression "virgin compression range" emphasizes the fact that the soil has not previously been subjected to pressures to this extent. The EI was calculated for the final loading step of 575 kPa with Δlog virgin compression range 1 3 Vol.: (0123456789) as used by Peth et al. (2010). The ratio ranges from 0 (completely plastic behaviour) to 1 (completely elastic behaviour).

Penetration resistance experiment
For the PR experiment a stainless-steel cover was devised to ensure correct positioning of the tests with sufficient distance between penetration sites and cylinder edge, and to prevent evaporation during the experiment, as only the hole in use was left uncovered. As the samples were sufficiently large to accommodate 12 penetrations, the experiment was performed at three different Ψ m (-3, -12.5, and -50 kPa Ψ m and an additional air-dry step in 2020) with three replicates per cylinder. For air-drying, samples were placed on a wire mesh and put underneath a fume hood at room temperature until no further weight loss occurred. A material testing machine (100 kN Allround Table Top Zwick/Roell, Ulm, Germany) equipped with a sensitive microsensor with a nominal force of 10 N (accuracy grade 1 according to ISO 7500-1 (ISO 2018) down to 0.02 N), was employed to push a penetration needle at a constant rate of 120 mm h -1 to a depth of 20 mm into the soil. The effect of the penetration rate on force measurements had been evaluated beforehand by comparing the forces resulting from real root growth velocities to higher ones which are more practicable to apply in the lab (Rosskopf et al. 2022). The results suggested that at a rate of 120 mm h -1 the resulting forces can be regarded as representative for those obtained at real root growth rates. The penetrometer probe was non-recessed with a diameter of 1 mm (Oleghe et al. 2017) and a 15° semi-angle resembling root geometries (Ruiz et al. 2017). The use of a non-recessed shaft allowed us to retract the PR needle at insertion speed which was performed for one of the three replicate measurements per cylinder. Displacement and forces were logged every 10 μm. Results of three measurements per cylinder at each Ψ m were averaged and corrected for shaft friction. PR was calculated as follows: with F Z,m being the measured axial penetration force and A the cross-sectional area of the PR needle. To calculate the mechanical energy demand according to Ruiz et al. (2017), the following formula was employed: with dz representing the incremental length, and the limits of 2 and 18 mm chosen to guarantee full contact between the cone and the soil, with values for all replicates being available. Results were subsequently related to the length of one metre. Samples were weighed after drying for 48 h at 105°C to calculate ρ b .

Statistical analysis
To test the effects of substrate and genotype and their interaction on the variables measured, data were sorted according to year and depth, and multi-factorial ANOVA was calculated after preconditions were verified. Where applicable, a subsequent Tukey HSDtest was performed. The same approach was applied to determine differences between depths and years. The energy data derived from PR measurements did not meet pre-requisites for ANOVA, so data transformations had to be carried out beforehand. For this purpose, data for the groups according to depth and year were log-transformed prior to ANOVA, subsequently significant differences for the factors substrate, genotype and their interaction within each Ψ m were quantified. Airdry samples were analysed separately, as only sand values in 2020 were available.
Here, one-factorial ANOVA was used to determine differences according to genotype and to substrate. Splitting the data into loam and sand we then calculated whether the factors depth and year caused any differences within energy values. For all tests, H 0 was specified as no differences occurring between observation groups and the significance threshold (

Effect of genotype and substrate on bulk soil mechanical parameters
Result tables for all parameters including mean values and standard errors for all groups can be found in the online resource (Table S1) alongside full ANOVA tables (Tables S2-S5). The impact of genotype and substrate on σ pc is graphically presented in Fig. 2. At 34 cm in 2019 the only differences in σ pc caused by the factor genotype ( Fig. 2a; p=0.033) with higher σ pc values for rth3 compared to WT and the only significant interaction (p=0.015) between the two factors (higher σ pc for WT in loam and for rth3 in sand; not shown). Values (omitting outliers) ranged from 19-52 kPa in loam and from 15-45 kPa in sand (Fig. 2c, d).
Higher values in loam compared to sand were found in all groups apart from 34 cm in 2019 (Fig. 2c). Compressibility was influenced by the factor substrate (Fig. 3), which led to strongly contrasting values with loam displaying a much higher compressibility than sand within the entire stress-strain relationship. C c were significantly higher in loam in all groups (Fig. 3a, b), the same being true for C s (Fig. 3c, d) (all p<0.001; full statistical tables are included in the supplementary information). Both ranges and inter-quartile ranges of C c in sand were remarkably low. In loam C s values differed considerably according to depth and year, whereas values in sand were more similar to each other. Contrarily, neither C c nor C s indices were affected by the factor genotype (not shown). No interactions occurred in any of the compressibility statistics.
EI were considerably higher in sand with mean values of 0.120-0.129 as opposed to loam with mean values between 0.045 and 0.052 (not shown). P-values were below 0.001, see supplementary information for statistical details. Again, no influence of genotype on EI could be discerned.
Significant correlations were observed in loam between ρ b and σ pc (Fig. 4a), C c (Fig. 4b), and C s (Fig. 4c). R 2 decreased and p-values increased from one year to the next, especially for σ pc the importance of ρ b as an explanatory variable was diminished with time. σ pc and C c (Fig. 5a) or C s (Fig. 5b), respectively, were negatively correlated (Fig. 5), whereas a positive relationship was found between C s and C c (Fig. 5c). Again, the same pattern occurred with strong relationships in 2019, and weaker and non-significant ones in 2020.
In sand the corresponding correlations were much weaker and mainly not significant (data not shown, see supplementary information S6).

Effect of depth and year on bulk soil parameters
The influence of the factors depth and year was analysed here more specifically and is graphically presented in Figures 6, 7, 8, along with indication of statistical significance. Full statistical tables are included in the online resource (Tables S7-S11). In sand, ρ b was significantly lower at 34 cm depth than 1 3 Vol.: (0123456789) at 14 cm (Fig. 6a), and the factor year did not result in significant differences (Fig. 6b). The opposite was the case in loam with similar values according to depth but significantly higher ρ b in 2020 compared to 2019. In none of the treatments σ pc (Fig. 7) were affected by depth, but it increased for both genotypes in loam from 2019 to 2020 (Fig. 7b).
Regarding the compressibility, depth did not influence C c (Fig. 8a, b), and higher C s values (Fig. 8c, d) were measured at 34 cm depth in loam, significantly for WT. Year turned out to be the dominant factor for changing soil stability leading to lower compressibility (tendency in C c in L rth3) in loam in 2020 compared to 2019. In sand, C c remained constant over depths and years with a remarkably low range. Lower C s values in 2020 could also be observed in sand, albeit not significant for rth3.
Effect of genotype and substrate on energy required for root growth Generally, the energy needed for one metre root growth increased with decreasing Ψ m . To better understand the effect of the factors genotype and substrate, data were grouped according to depth and year, as these groups displayed distinct patterns. The factor genotype proved to be of minor relevance with significant differences in only two groups out of twelve (Fig. 9). Significantly higher energy values for WT roots were observed in 14 cm 2019 at Ψ m -12.5 kPa  and in 34 cm 2020 at Ψ m -50 kPa. No interactions between substrate and genotype could be observed. Whereas in 2019 no significant differences occurred, in 2020 differences between energy values in loam and sand had developed in 14 cm at Ψ m -50 kPa and in 34 cm depth at Ψ m -12.5 and -50 kPa (Fig. 10b, d). Highest values of up to about 2 J (~2.5 MPa) were found at Ψ m -50 kPa in loam. The samples from 34 cm depth in 2019 displayed the lowest range of values (Fig. 10c). The respective statistical table is S12.
Additional graphs depicting the course of the PR curves are located in the online resource (Figs. S1 and S2). In 2019 at 34 cm depth all treatments displayed similar curves. The most pronounced differences occurred at -50 Ψ m kPa in 2020 with loam WT being distinctly higher than loam rth3 which in turn was higher than sand.

Effect of depth and year on energy required for root growth
When visually depicting the influence of the factors depth (Fig. 11a, c) and year on penetration energies (Fig. 11b, d) in loam and sand, a clear pattern becomes visible. At 14 cm depth values were higher in both substrates with significant differences for four out of six groups. Required energy in loam (Fig. 11b) remained constant over the years, whereas energy in sand (Fig. 11c) declined significantly from one year to the next. In loam an interaction between the factors occurred at Ψ m -12.5 kPa (p = 0.0207) with higher values at 14 cm depth compared to 34 cm in 2019 and vice versa in the following year. The inter-quartile range of both loam and sand values was reduced from one year to the next. ANOVA tables are available in the online resource (Table S13). Due to technical limitations, measurements at airdry conditions (not shown) could only be carried out in sand and only a few were done in loam, not enough to be analysed statistically. In sand, no differences caused by depth or genotype occurred. As the soil becomes drier, the plant must overcome much higher resistances and therefore invest more energy in root growth. In sand, PR of approximately 5 MPa were reached in both depths.
Regarding significant correlations between any of the other parameters with energy values, correlations could be observed between ρ b and energy in sand at different Ψ m (Fig. 12a, b, c). The respective correlations in loam were not significant and displayed no differences between years (data not shown). The air-dry sand samples of 2020 did not significantly correlate with ρ b .

Influence of genotype and substrate on bulk soil mechanical parameters
The bulk soil properties measured were not affected by genotype. Even examining the rhizosphere only, an influence of root hairs on ρ b could be neither detected in barley (Koebernick et al. 2018) nor in maize ). Roots play a crucial role in the formation of aggregates as they enmesh (c) (d) Fig. 9 Energy required for a root of a defined maize genotype (WT = wild type, rth3 = root hair-less mutant) with 1 mm diameter to grow 1 m. Data were grouped according to depth (14 and 34 cm) and year (2019 and 2020). Data were logtransformed for statistical analysis and subsequently rescaled to the original data range for plotting. Measurements were taken at three matric potentials (-3, 12.5, and -50 kPa). Different letters within each matric potential (n =24) signify differences between genotypes particles and release organic compounds into the rhizosphere. Aggregates are stabilized by bacterial colonization and wetting-and-drying cycles, especially if clay particles are present (Bronick and Lal 2005). Koebernick et al. (2017Koebernick et al. ( & 2018 concluded that barley root hairs have a stronger effect on the inter-aggregate pore space, creating a higher pore volume and increasing the percentage of small pores in the rhizosphere thus compensating for the compaction created by roots. These processes happen on a very local scale, possibly too small to be able to influence bulk soil measurements such as the ones in this study. Looking at root hairs and their influence on mechanical properties, especially tensile strength, on the aggregate scale is a topic we aim to study at this site over the next years. Geological history and mechanical loading are factors influencing σ pc which can be ruled out in the present setup, because the plots were artificially filled and all use of machinery for agronomic measures was avoided, and the soil was not tilled. This provided us with a unique set-up to measure the development of soil mechanical parameters. Even though the samples in the present experiment were taken as undisturbed samples, they are -technically speaking -remoulded samples, as the plots were artificially packed with homogenized substrate. We drained all oedometer samples to Ψ m -50 kPa to provide a comparable Energy required for a root with 1 mm diameter to grow 1 m in the respective substrate (loam and sand). Penetration resistance is displayed as secondary y-axis. Data were grouped according to depth (14 and 34 cm) and year (2019 and 2020). Data were log-transformed for statistical analysis and subsequently rescaled to the original data range for plotting. Measurements were taken at three matric potentials. Different letters within each matric potential (n = 24) signify differences between the substrates starting point for our measurements (Figures 2, 3 , 4, 5, 6, 7, 8). In the present experiment the factor substrate had a major impact on σ pc with higher values in loam compared to sand. This is in accordance with our hypothesis, as pronounced aggregation in loam leads to increased mechanical soil strength (Horn et al. 1994), especially at clay contents of more than 15% (Lebert and Horn 1991) which was the case for the loam used. Soil strength is a function of forces between particles, interparticle friction and cohesion as well as the number of particle contacts per volume. In a single grain structure as for sand, the number of particle contacts and the forces between particles are much lower compared to a loam with a variety of particles sizes (Lebert and Horn 1991). Thawingand-freezing as well as drying-and-wetting processes Loam data were log-transformed prior to statistical analysis and were rescaled to the original scale for plotting. Significant differences are indicated for each matric potential, n =24 (a) (b) (c) Fig. 12 Correlations between bulk density and energy at matric potentials of -3, -12.5, and -50 kPa in sand, data were split by year, n = 24. Coefficient of determination R 2 and respective p-values are given further increase soil stability (Dexter 1988), additional factors are the presence of organic matter, biological processes, and the consequences of anthropogenic impacts (Dexter 1988;Horn et al. 2019). Horn et al. 1994 stated that σ pc in a homogeneous substrate corresponds to the effective stress as a function of pore water pressure. In a structurally unstable soil, the previous highest hydraulic stress (most negative pore water pressure) affects σ pc (Mosaddeghi et al. 2003). In loam we found a positive correlation between ρ b and σ pc (Fig. 4a), which is in accordance with An et al. (2015) and Mosaddeghi et al. (2003). Rücknagel et al. (2007) used ρ b and aggregate density as input for multiple linear regression to estimate σ pc and revealed that an increase in ρ b results in higher σ pc . According to Lebert and Horn (1991) with an increment in aggregation caused by a higher clay content, the importance of ρ b as determining factor for soil strength decreases. Interestingly, in the present study no significant correlation between ρ b and σ pc existed for sand. We noted that the sand samples were in fact very unstable from the beginning. According to Dexter (1988), one of the critical issues of soil stability is its ability to resist the influence of water. During sample preparation in our study (saturation to a standard Ψ m ), sand samples consolidated considerably and lost on average 1.72 mm of their initial height. This led to an unavoidable stabilisation of these samples with a mean decline in ρ b of 0.11 g cm -3 , a decrease in C s and C c , and an increase of σ pc compared to actual field conditions. The compressibility of loam was higher compared to sand both along the virgin compression range (C c ) (Fig. 8a, b) and the recompression range (C s ) (Fig. 8c,  d). This resulted in more defined curves in loam, whereas low initial void ratios and their weak decline during compression dominated in sand, which is in agreement with Gregory et al. (2006). The sand with its single grain structure showed only little potential to be compressed as the majority of grains were of similar size and the pore system dominated by primary pores. This resulted in a relatively small pore volume that could be compressed, with few soil particles that were small enough to occupy these pore spaces. This is in contrast to the loam with a higher variety of pore sizes as well as of particle sizes leading to greater compressibility. The small ranges of compressibility in sand reflect the low structural development in this substrate as opposed to loam. Low initial void ratios resulted in higher resistance to compression and therefore decreasing C c and C s and higher σ pc as higher number of particle contacts lead to increased frictional forces, less available pore space and more energy necessary to relocate water during compression if the initial void ratio is low. This is in accordance with a study by Keller et al. (2011) that revealed initial void ratio as a determining factor for σ pc , C c , and C s. Keller et al. (2011) further argued that using a semi-logarithmic curve as basis for calculating soil compression properties has a major impact on these values. However, as it is considered and accepted as the common approach it was also applied here.
In loam, ρ b was negatively correlated with C c and C s ; the corresponding correlations in sand were weaker than we expected as compressibility had already been reduced due to consolidation (Fig. 4) as discussed above. A negative correlation between ρ b and C c for two agricultural soils was reported by An et al. (2015).
EI was higher in sand than in loam which at first glance seems to be an unexpected result. However, EI was based on a final un-loading step after maximum loading of 575 kPa. The EI relates the changes in void ratio of this final step to the overall changes in void ratio of the entire curve, which were much higher in loam compared to sand. Therefore, absolute void ratio changes in the unloading step were in fact higher in loam than in sand, but the resulting ratio, i.e., the EI, was lower.

Spatial and temporal development of bulk soil mechanical properties
We hypothesized structural development to be higher at 14 cm than in 34 cm depth, because closer to the soil surface the influence of environmental factors associated with more intensive wetting and drying can be expected to be stronger. Anyhow, in our study, differences between depths were not very pronounced and only occurred for sand, where ρ b was lower at 34 cm depth and C s in loam with higher values at lower depth. Rain was scarce in both years and as the bottom of the plots (1 m) was filled with a 25 cm gravel layer with a drainage textile placed on top, plants could not access water from below the plots but depended on precipitation alone. Both the effects of precipitation and drought start from the top thus increasing the variability of Ψ m in the topsoil layer 1 3 Vol.: (0123456789) resulting in differences in wetting-and-drying cycles according to depth (Jorda et al. submitted).
The temporal development was more prominent than the spatial ones. A general setting of soil could be observed with higher ρ b in 2020 compared to 2019 in loam, which is not surprising regarding the substrates had been filled-in homogeneously. In loam σ pc increased and C c decreased with time. σ pc and C c values in the same loam in a lab experiment employing remoulded samples (Rosskopf et al. 2022), which can be considered analogous to initial field conditions of the current experiment, fit in well with this timeline. The potential for further compaction, i.e., reduction of the pore space, had already been greatly reduced in sand in 2020, which becomes evident when looking at the extremely narrow ranges of C c . Together with declining C s in both substrates, this coincides with higher root length densities in all treatments in 2020 compared to 2019 (Vetterlein et al. 2022). Increased root length densities bring about a stabilization of structure, hence higher σ pc and lower C c and C s . This corresponds to the declining relevance of ρ b for the bulk soil mechanical parameters and the decreased correlations between them. The general development of correlations regarding bulk soil measurements getting weaker and less significant over time, points to a diversification of explanatory variables.

Required energy and penetration resistances
It is commonly known that decreasing Ψ m leads to an increase in forces needed for soil penetration (Quang et al. 2012;Elbanna and Witney 1987;Wang et al. 2016). The wettest samples displayed no differences between treatments for substrate and genotype, whereas differences according to substrate were more pronounced at lower Ψ m as PR responds stronger to Ψ m at higher clay contents (Costantini 1996submitted). With the current experimental set-up, it was not possible to me). Both sampling years were very dry, with cumulative precipitation amounts of 180 mm in 2019 and 210 mm in 2020 from sowing until sampling at BBCH83 (German Meteorological Service 2021). This resulted in low soil moistures dropping even below the permanent wilting point (Jorda et al. asure PR at soil moistures around the permanent wilting point in loam, as the penetrometer needle did not resist the occurring forces. In order to be able to tackle this problem in the future, we are presently working on a comparison of a 1 mm and a 2 mm diameter needle, which will allow us to measure PR at drier conditions and relate the results to the ones obtained so far. Nevertheless, the few measurements that were made in airdry loam give us an idea of the resistances the roots encountered in the field. PR values of over 50 MPa were reached in loam in five measurements out of 17. Even in sand average values around 5 MPa were around five times higher compared to measurements done at -50 kPa Ψ m . A lab experiment using the same loam yielded energies twice as high at the permanent wilting point compared to -50 kPa Ψ m (Rosskopf et al. 2022).
When transferring the results to natural root systems, it has to be taken into account that a real root does not grow straight into the soil as the PR needle moves. Natural occurring roots utilize the available pore system following paths of least resistance to avoid compacted regions. Nevertheless, root-soil contact is indispensable for obtaining nutrients and water. Our micro penetrometer approach pays attention to root diameter, and cone angle and employs penetration rates which are low enough to prevent any soil plastic viscous effects that could occur otherwise. The forces exerted at the penetrometer conus are recorded and corrected for the friction at the interface between penetrometer shaft and soil. A recessed shaft would ensure that the measured forces only result from the tip but using such a shaft oftentimes does not entirely prevent soil -depending on its texture and moisture -from getting re-attached to it, which then constitutes an error difficult to quantify. That is why a nonrecessed shaft was used, which can be pulled out at insertion speed, thus allowing us to estimate shaft friction at any point in time during the measurement and to correct the original forces with these values. The forces encountered during the tip-wise growth of roots can therefore be accurately represented in our very localized measurements. A major difference between an artificial penetrometer and a naturally elongating root remains: the reduction of friction through the secretion of mucilage and the sloughing of root tip cells. Iijima et al. (2004) measured a 30% decrease in the resistance a root experiences due to the presence of an intact root cap and the secretion of mucilage. Bengough and McKenzie (1997) quantified that a pushed maize root experiences about 40% of the resistance a metal probe does, and a growing root 50-100% of the resistance of a pushed root. Keeping 1 3 Vol:. (1234567890) in mind that differences due to frictional aspects exist, we can conclude that with our experimental setup we do get close to what roots potentially experience in the field.
Generally, we have to state that genotype has a more pronounced influence on PR with a tendency to higher PR in WT, as it has on the measured bulk soil mechanical parameters. Within the PR tests, failure is not solely attributable to compression as it is the case for σ pc , but also to tensile strength of aggregates. Factors that lead to aggregation and stabilisation of aggregates, such as deposition of organic matter in the rhizosphere leading to an increase in tensile strength, play a major role on this more local scale. Root hairs can improve a crop's ability to deal with drought stress as Marin et al. (2021) showed with barley genotypes, but their impact on soil mechanical parameters at the field scale has not been considered so far. A tendency to higher energy values could be observed in WT, which was also the genotype with higher root length densities at all depths at BBCH83 ( fig. 5 in Vetterlein et al. 2022). Our undisturbed samples were taken within the depth intervals of this project partner, which did not show a pronounced depth gradient, thus enabling us to refer to them directly. In particular, the combination of loam and WT resulted in higher forces compared to the other treatments at more negative Ψ m. WT had significantly higher shoot dry weight compared to rth3 ( fig. 3 in Vetterlein et al. 2022) in both years, resulting in higher water demand for WT and subsequently an earlier onset and a more severe drought stress. This process is reflected in the Ψ m measurements in sand made by Jorda et al. (submitted). A possible explanation for higher energy values in WT might therefore be both the direct impact of increased rooting and the indirect effect of water depletion through roots.
In the laboratory experiment mentioned above, which employed the same loam and sand (Rosskopf et al. 2022), PR and energy were measured at the water contents corresponding to the same Ψ m as used here, thus resembling a perfect homogeneous state that might come close to initial field conditions. In the present study, energy values in loam and sand were approximately an order of a magnitude higher than in the lab, indicating a process of structural formation in the field. A further explanation for this can be found in the respective sieving and filling procedures resulting in much more homogeneous samples in the lab, where substrates were sieved to 1 mm and mixed and filled in in small portions, and in slightly lower ρ b compared to the field. Reproducing relationships between PR, ρ b and soil moisture in the field with repacked samples in the lab was not possible due to structural variability as shown in a study by Costantini (1996).
In all treatments penetration energy was higher at 14 cm depth as opposed to 34 cm which is in line with our hypothesis of the evolution of structure being more pronounced closer to the soil surface due to environmental factors. We assume the higher soil moisture variability and therefore more frequent drying out processes up to the permanent wilting point (Jorda et al. submitted) to be of major importance in this context. A further aspect might have been higher root length densities at BBCH83 at 14 cm depth in sand in both years and in loam in 2020 (Vetterlein et al. 2022).
Our hypothesis that advancing structural development leads to an increase of required energy over time could not be affirmed, in fact penetration energies for sand decreased from one year to the next, while for loam no changes were observed. With time, differences between the substrates emerged, indicating a stronger structural development in loam, as hypothesized. It seems that in loam areas of higher and of lower resistance have balanced each other out, so that differences between years did not appear. In sand, higher root length densities and non-decomposed roots from the previous year ( fig. S2 in Vetterlein et al. 2022) might have contributed to the decrease observed, as roots have a smaller effect on aggregation in sand than they do in loam. Not only the adjusted Ψ m but also maximum pre-drying plays an important role for the mechanical strength of soils (Hartge 1986). Precipitation was altogether higher in 2020 and as much less water is needed to fill up the sand plots with their higher ρ b and lower porosity compared to the loam plots, this resulted in the sand being wetter in 2020 than in 2019 and compared to loam and brought about the very uniform force-displacement curves in theses samples. Wang et al. (2016) observed that with subsequent wetting-and-drying cycles soil structure evolves resulting in higher heterogeneity of strength and a simultaneous decrease of both maximum and overall PR values.

CONCLUSIONS
The overall effect of genotype on bulk soil mechanical properties considered proved to be negligible, whereas PR might be stronger impacted by roots. Substrate with its implications for soil structure was the determining factor for stability, compressibility, and elasticity in the two years following the establishment of the present field experiment. As expected, stronger changes of soil mechanical parameters could be found in loam as opposed to sand. After only a year, higher σ pc and lower compressibility, especially C s , were measured, indicating a stabilisation of structure with time in loam but only marginally in sand, in which measurements showed an altogether more uniform behaviour. The sand proved problematic for our measurements, as it consolidated during saturation, distorting the stress-strain relationship and the resulting parameters. Differences of bulk soil parameters between depths were less pronounced but indicated stronger structural development in the top layer. Regarding energy needed for root growth, higher values were found in loam and partially in WT with stronger differences between substrates in drier conditions. Higher PR in the top layer confirmed our hypothesis of stronger stabilization processes close to the soil surface. Temporal development did not change PR in loam and decreased it in sand. In a further field campaign, an adjusted technical set-up with stronger penetration needles will provide us with more detailed information about PR at drier conditions, also accounting for small-scale heterogeneity of values within samples and including a methodological comparison of 1 mm and 2 mm diameter PR needles. Furthermore, we are evaluating PR values in relation to root length densities in a column experiment using substrates and genotypes as in the present study to give us further insights on the effect of the actual root system on PR.