Climate change is expected to impact the water cycle and severely affect
precipitation patterns across central Europe and in other parts of the
world, leading to more frequent and severe droughts. Usually when projecting
drought impacts on hydrological systems, it is assumed that system
properties, like soil properties, remain stable and will not be affected by
drought events. To study if this assumption is appropriate, we address the
effects of drought on the infiltration behavior of forest soils using dye
tracer experiments on six sites in three regions across Germany, which were
forced into drought conditions. The sites cover clayey-, loamy- and sandy-textured soils. In each region, we compared a deciduous and a coniferous
forest stand to address differences between the main tree species. The
results of the dye tracer experiments show clear evidence for changes in
infiltration behavior at the sites. The infiltration changed at the clayey
plots from regular and homogeneous flow to fast preferential flow. Similar
behavior was observed at the loamy plots, where large areas in the upper
layers remained dry, displaying signs of strong water repellency. This was
confirmed by water drop penetration time (WDPT) tests, which revealed, in
all except one plot, moderate to severe water repellency. Water repellency
was also accountable for the change of regular infiltration to fingered flow
in the sandy soils. The results of this study suggest that the
Soils moderate how water moves through the vadose zone and govern the percolation of water to groundwater and streamflow. Soils not only store water for plant growth, function as a habitat for different biota and as transition zone to groundwater, but also are important – especially the top layers – for sorption and degradation of contaminants and (agri-)chemicals (Hendrickx and Flury, 2001). The efficiency of this important ecosystem service for groundwater and surface water protection depends on the behavior of pollutants in the soil and the hydrological transport processes (Keesstra et al., 2012). How fast water passes the vadose zone depends on its hydraulic soil properties and distribution such as pore volume distribution, soil aggregation, water repellency and rooting pattern.
Due to climate change and increasing human intervention, the global water cycle is expected to change with summer dryness and winter wetness in many regions across the world probably increasing, including western and central Europe (IPCC, 2012; Prudhomme et al., 2014). In addition, droughts are expected to be more frequent and severe in the future (Prudhomme et al., 2014; Seneviratne et al., 2006). Drought conditions can alter the hydrological functions of soils, and soil structure responds to drought with shrinkage and fracturing of soil aggregates. These soil shrinkage cracks channel the infiltrating water and therefore foster the bypassing of the soil matrix (Hendrickx and Flury, 2001; Ritsema et al., 1997) and alter the infiltration patterns in soil. Thus, the infiltration and redistribution of water within the soil changes and hence also the proportion of water reaching the groundwater (Hendrickx and Flury, 2001).
Soils under drought conditions are prone to become water repellent, depending on soil properties and organic matter content (DeBano, 1981, 2000). Due to modifications of the three-dimensional distribution and dynamics of soil moisture, water repellency has far reaching consequences for infiltration processes (Doerr and Ritsema, 2006). Water repellency hinders infiltration and thus either increases overland flow (Doerr and Ritsema, 2006) or redirects the water into preferential flow paths and creates instable wetting fronts (fingered preferential flow; Ritsema et al., 1993; Ritsema and Dekker, 2000; Dekker and Ritsema, 2000).
To assess the impacts of drought, rainfall exclusion experiments are
valuable and often applied tools (e.g., English et al., 2005; Phillips et al.,
2009; Da Costa et al., 2010; Kopittke et al., 2014), in addition to elevated
CO
To study drought effects, often extreme short-term events equivalent to droughts with occurrence probabilities of up to 100 or even 1000 years have been introduced to the examined soils (e.g., Glaser et al., 2013). By introducing these extreme events, the question of transferability of the results to natural systems in respect to the expected behavior under predicted future drought conditions arises. Therefore, this study employs a moderate rainfall reduction equivalent to an annual drought with a 40-year return period, in accordance to climate predictions, thereby avoiding tentativeness due to an overreaction to an unnatural extreme drought (Gimbel et al., 2015).
Location of the study sites (red squares) within Germany (light blue); southwest: Schwäbische Alb; center: Hainich–Dün; northeast: Schorfheide–Chorin.
To monitor changes in soil hydraulic properties, the changes in infiltration patterns in the soil after 2 years of prolonged drought were observed in three regions across Germany. Infiltration patterns were chosen because they reflect the integrated changes of soil hydrological functions and directly show how water moves in the soil under altered conditions. In this paper, we present results of several dye tracer infiltration experiments before and after 2 years of prolonged artificial drought. The objectives of this study are (i) to investigate whether droughts predicted by climate projections affect the infiltration behavior of forest soils, and (ii) whether changes in infiltration patterns can be attributed to changes in the hydrologic properties of the soils. Three hypotheses will be tested: induced drought alters infiltration patterns due to changes in soil hydraulic properties, e.g., soil water repellency and the formation of shrinkage cracks, leading to preferential flow paths and faster infiltration; (2) the main tree species have an effect on the magnitude of the observed response and (3) the drought will increase water repellency depending on tree species and soil properties.
To identify the influence of drought on infiltration patterns of forest
soils, six plots in three different regions across Germany were selected.
The plots were located in Schwäbische Alb (southwest Germany),
Hainich–Dün (central Germany) and Schorfheide–Chorin (northeast
Germany) (Fig. 1). All plots are part of the Biodiversity Exploratories
framework that incorporates, in total, 150 sites on grassland and 150 sites
in forest (for more information on the Biodiversity Exploratories, refer to
Fischer et al., 2010). In each of the frameworks, two forest plots were
selected, which are – within each framework – similar with respect to
topography and soil texture type (Fig. 2) but differ in tree species
composition. In each site, one plot with a coniferous and one with a
deciduous main tree species was selected. At the Schwäbische Alb and
Hainich–Dün sites, beech (
Soil horizons, texture and rock fractions of the six experimental plots. Soil type classification according to the World Reference Base for soil (FAO, 2006).
The Schwäbische Alb soils are shallow (25 to 35 cm) Leptosols on
Jurassic shell limestone with a high stone content (Fig. 2, top panels). The mean
annual temperature at this site is 6.5
The experiments of this study are part of the interdisciplinary project
“Global Change Effects on Forest Understorey: interactions between Drought
and Land-use Intensity” (Gimbel et al. 2015). The artificially imposed drought
was created by a 10 m
To observe the impact of reduced precipitation input on soil moisture, soil
moisture probes were installed on the drought and on the control subplots of
every site. The probes (5TM and 5TE, Decagon Devices Inc.) were inserted at
5, 15 and 30 cm depths in three replicates on the plots at 2, 3, and
4 m distance from the central trees. The accuracy according to the technical
data sheets of the 5TE and 5TM probes is
Scheme for profile excavation
Hydrophobicity in soil was measured with the water drop penetration time (WDPT) test (e.g., Bisdom et al., 1993). This test determines how long water repellency persists on a porous surface. The tests were performed immediately before the dye tracer experiments in 2013, in the drought and control profiles of the deciduous plots and in the drought profiles of the coniferous plots. For the WDPT tests, a water droplet is placed on a planar soil surface with a pipette and the time is measured until the water drop is completely taken up by the soil. The observation was stopped after exceeding a time of 3600 s. Depending on the profile depths, WDPT tests were performed in several depths of the profile, covering the main soil horizontation. In each depth, five sampling locations were used to traverse the profile, and the tests were repeated 3 times per location, resulting in 15 WDPTs per depth (Fig. 3b). The mean and maximum values of the WDPT test were classified after Bisdom et al. (1993) (Table 1).
The dye tracer experiments were conducted in August 2011 before installation
of the roofs and in August 2013 after 2 years of drought. For each
experiment, an area of 80 cm
The next day (after waiting at least 12 h), three vertical soil profiles
per sub-area were prepared. Keeping a 10 cm buffer stripe at the beginning
and between the individual sub-areas, every sub-area was divided in three
sections, spaced 10 cm from each other (Fig. 3a). To obtain the dye
pattern, the surface of the excavated soil profiles was smoothed with a
spatula and loose particles were removed with a brush, avoiding smearing.
Stones were left in place and shaped into relief when needed. Roots were
trimmed. Pictures were taken from each profile with a standard digital
compact camera with a resolution of 10 megapixels (3648
Classification of water repellency by WDPT, after Bisdom et al. (1993).
To objectively compare the flow pattern of the different profiles, we used
the image analysis procedure developed by Weiler and Flühler (2004). We
provide a short description of the process here, for more detailed
information, refer to Weiler (2001) and Weiler and Flühler (2004). The
image processing consists of three main steps. In the first step, geometric
distortion of the image is corrected by establishing a relationship between
the image pixel location and the true location on the soil profile. During
this step, the image is also scaled such that, one pixel corresponds to a
square of 1 mm
To obtain objective measures to compare the dye patterns of the different
profiles and sites, we derived three depth-related variables of the binary
images: (1) volume density, (2) surface density and (3) stained path width
as basis for further delineation of flow processes. The volume density (VD)
is similar to the frequently used dye coverage. It is defined as the stained
volume divided by the reference space and originated from the methods of
stereology, which relates a three-dimensional parameter to two-dimensional
measurements (Weibel, 1979). Surface density (SD) is defined as surface area
of an object divided by the volume of the reference space. Surface density
provides information on the size and number of features: a high SD is caused
by a large number of small objects, whereas a low SD indicates less but
larger objects (Weiler, 2001). The stained path width (SPW) is derived by
measuring the width of every stained object at a certain depth. The SPWs of
every depth were classified into three classes of
Normalized cumulated sums of soil moisture of the drought versus the control subplots of the investigated soils. The computations begin with the start of the rainfall reduction and end at the date of the dye tracer experiment 2013.
Mean and maximum water drop penetration times (WDPTs) of the control (green) and drought (red) plots. Orange lines and numbers refer to the WDPT classes after Bisdom et al. (1993) (see Table 1).
Figure 4 shows the normalized cumulated sums of the soil moisture measurements of the control and the drought plots over the course of 2 years. All plots developed a soil moisture deficit compared to the control plots in the upper 5 cm of the soil, as shown by the black line below the 1 : 1 line. The water deficit is also transduced to 15 and 30 cm depths in both Schwäbische Alb plots and in the coniferous plot of Hainich–Dün, but is generally less pronounced. The plots at the Schorfheide–Chorin site show no deficit (deciduous plot) or even a small plus in soil moisture (coniferous plot) compared to the control plot. The sandy soils of Schorfheide–Chorin are already very dry without drought treatment. The reverse moisture effect might be caused by effects associated with the installation and placement of the probes (e.g., preferential flow along probe cables, concentration of throughfall – although installations were performed with great care), or root effects, e.g., hydraulic redistribution. However, we did not find any signs for hydraulic redistribution in the data. The deciduous plot of the Hainich–Dün site experienced major probe failures due to animal damage during the summer months of 2012 and again in 2013. Therefore, only the data taken during the winter months could be used for the comparison. For this reason, the data do not cover the months with the highest expected soil moisture deficits.
Figure 5 shows the results of the WDPT test in 2013. All drought-treated plots at all sites – coniferous and deciduous – exhibit water repellency (WDPT data from control plot under coniferous not available). All control plot soils are wettable (WDPT class 1) or feature at least lower water repellency than the drought-treated plots. The coniferous plots under drought of Hainich–Dün and Schwäbische Alb showed higher WDPTs than the deciduous plots. This is valid for both mean and maximum values. The Schorfheide–Chorin deciduous plot showed higher water repellency than the coniferous plot. In all soil profiles, water repellency is the highest in the topsoil and diminishes at a depth of about 20 cm. However, in the Schorfheide–Chorin deciduous plot, water repellency is present up to a depth of 50–60 cm. When present, strong to severe water repellency is dominant in the measured drought-treated plots. Only the Hainich–Dün deciduous plot soil is classified as wettable on average and the Schorfheide–Chorin coniferous plot as slightly water repellent. The highest values in mean and maximum water repellency were found in the coniferous plots Hainich–Dün (mean 941 s; max 3600 s) in about 10–15 cm depth and in the Schwäbische Alb (mean 990 s; max 2340 s) in the topsoil (Fig. 5).
Differences between pre-drought and control plots (without drought treatment) not only reflect differences in soil structure, texture and moisture due to a distance of 20–40 m between the drought and reference plot, but also may include time-dependent changes of the soil characteristics, which are independent from the drought treatment. To ensure validity of the dye pattern analyses, it is necessary to assure comparability among the plots. To exclude time-dependent changes as reasons for differences in pre-drought and drought-treated dye patterns, the pre-drought patterns were checked against the patterns of the control plots. Figure 6 compares the pre-drought pattern and the control pattern of the deciduous plots. In addition, Fig. 7 provides box plots of VD for different depths of the pre-drought and the control profiles for direct comparison.
The Schwäbische Alb pre-drought plot (Fig. 6, top left panel) shows high VD in the top 10 cm in all profiles. The 40 and 60 mm sprinkling volume profiles show a high SPW in the top 5 to 10 cm. On the control plot (Fig. 6, top right panel) large areas of the profiles top 10 cm are also stained, but the VD and SPW are not as high as the pre-drought profiles. This is especially evident in the VD box plots of the upper 0–10 cm (Fig. 7, top panels). All Schwäbische Alb profiles have high stone contents, in some cases exceeding 50 % of the profile width (Fig. 6). Below 10 cm depth, the control plot profiles are almost completely stained. This pattern is similar to the pre-drought profiles. In general, the patterns of the control profiles are similar in VD, SPW values and distribution to the 20 mm pre-drought profile. The 60 mm control profile reflects the high VD and SPW values in top layers, which are characteristic of the 40 and 60 mm pre-drought profiles.
The Hainich–Dün pre-drought profiles (Fig. 6, center left panel) show low to medium SPW in all depths. VD values are high in the top 5 cm in all profiles and between 10 to 30 cm in the 40 and 60 mm sprinkling amount profiles. The 20 mm profile displays only small VD values below 10 cm depth. All profiles have a medium to high stone content (30–60 %) below 30 cm depth. The control plot profiles (Fig. 6, center right panel) are very similar in VD and SPW to the pre-drought profile pattern, but with generally lower VD in the 60 mm sprinkling amount profile (Fig. 7, center right panel). Except for the 20 mm profile, which displays no stones, the control plot profiles have a medium to high stone content below 25 cm depth. In all profiles, large areas of the profile stayed unstained. However, although having a low VD in top layer, the 60 mm control plot profile is not following the pronounced drop in VD between 5 and 10 cm depths and the subsequent rise between 15 and 25 cm, which is characteristic for all other profiles (pre-drought and control). These distinct differences are apparent in the box plots (Fig. 7).
In the Schorfheide–Chorin pre-drought profiles (Fig. 6, bottom left panel), high VD and SPW values are present. The highest VD and SPW values can be found in the 60 mm sprinkling amount profile. Below 10 cm depth, the 20 mm pre-drought profile displays only small to medium SPWs and – in comparison to the 40 and 60 mm profiles – small VD values. The control plot profiles (Fig. 6, bottom right panel) show, in general, high VD and SPW values, but have lower values in the top 10 cm than the pre-drought profiles (Fig. 6, bottom right panel and Fig. 7, bottom panels). This is more apparent in the 20 and 40 mm profiles (Fig. 7). In the pre-drought and control plots, infiltration reached down to depths of over 70 cm and no stones are present.
To summarize, the comparison between the pre-drought and control plots showed a broad agreement. Differences, that need to be accounted for, are the lower VD in the profile top layers, especially at the Schwäbische Alb and Schorfheide–Chorin site. These differences might be due to spatial heterogeneities, e.g., slight differences in soil layer boundary depths, given the distance between the control and the pre-drought plots (15 to 30 m). In addition, the initial conditions (soil moisture) were also slightly different possibly resulting in the observed differences. Choosing 10 cm steps for statistical comparison of the VD may in addition introduce differences, if soil layer boundary depths differ. Therefore, not only the VD, but also the SPWs and the determined flow processes need to be taken into account for comparison. However, the pre-drought and drought experiment were performed in close vicinity (1 m). In the Hainich–Dün, the drop and rise of VD in all profiles points to a soil layer boundary effect on infiltration. This is not time dependent and present in both pre-drought and control profiles; therefore, the comparability between the pre-drought and drought pattern is not affected.
As can be seen in Fig. 8, all plots show marked differences between pre- and after-drought infiltration patterns. The clayey and loamy sites (Schwäbische Alb and Hainich–Dün) develop unstained (i.e., unwetted) areas in the topsoil layers. This is more pronounced in the coniferous plots, where unstained areas are already visible in the pre-drought infiltration pattern. Figure 9 compares VDs in box plots for different depths of the drought and pre-drought profiles including the statistical significance.
At the Schwäbische Alb site, medium to low volume densities (VDs) were found on the pre-drought coniferous plot throughout the whole profile for the 20 and 40 mm sprinkling depth and high VD for 60 mm sprinkling depth (Fig. 8, top left panel). The drought 40 and 60 mm profiles are lower in VD in the top layers (0–10 cm), than the pre-drought profiles (Figs. 8 and 9, top left panels); the 40 mm profile displays even unstained areas (no VD). The 20 mm pre-drought profile is already very low in VD; therefore, the differences to the after drought profile are not distinct (Fig. 8, top left panel). The drought coniferous plot shows a rise of VD culminating around 20 cm depth (Fig. 8, top panels) for all sprinkling amounts (20, 40 and 60 mm). Below 20 cm depth, the 20 and 40 mm profiles show (Fig. 9, top left panel) higher VD in the after drought profiles than in the pre-drought profiles, whereas the 60 mm profile show the same extent of VD in the drought and in the pre-drought profile.
Comparison between stained path width (SPW) of pre-drought (2011) and control (2013) plot. Blue shades indicate the SPW classes. The sum of SPW is the volume density (VD) per depth. Grey and black indicate the VD of stones.
VD box plots of the drought and the pre-drought pattern. Depth
ranges are omitted, where one of the profile is shorter than the other.
Statistically significant (
Comparison between before drought (2011) and after drought (2013) stained path widths (SPWs) and flow processes for coniferous and deciduous stand plots. The sum of SPW is the volume density (VD) per depth. Grey and black indicate the VD of stones.
VD box plots of the drought and the pre-drought pattern. Depth
ranges are omitted, where one of the profile is shorter than the other.
Statistically significant (
The stained path widths (SPWs) of the Schwäbische Alb coniferous
pre-drought profiles are small to medium in the 20 and 40 mm profiles and
high in the 60 mm profile (Fig. 8, top left panel). After drought, low to medium
SPWs are dominant in the 20 and 60 mm profiles; high SPW values occur in the 40 mm profile below 20 cm. The flow processes identified in
this depth as matrix flow, are caused by local saturation due to low
The Schwäbische Alb deciduous plot shows in the 40 and 60 mm pre-drought profiles high SPWs and in all infiltrating volumes high VDs in the top layer (0–10 cm; Figs. 8 and 9, top right panels). Medium to high VDs are maintained throughout the whole 40 and 60 mm profiles, and to lesser extend in the 20 mm profile. The drought profiles show lower VD in the top 10 cm, compared to the pre-drought profiles (Fig. 9, top right panel). Below 20–25 cm depths, the 20 and 40 cm drought profiles show higher VD than the pre-drought profiles. However, the drought profiles are more similar in shape to the VD pattern of the control than to the pre-drought profiles (Fig. 6, top panels). Also, the stone contents in the three pre-drought profiles are higher than in the drought profiles (Fig. 8, top right panel).
The dominating flow types in the Schwäbische Alb deciduous pre-drought
profiles are identified as macropore flow with low, mixed and high
interaction, and as matrix flow, depending on soil layer and infiltration
volume (Fig. 8, top right panel). The flow processes identified in the top
layers of the 40 and 60 mm pre-drought profiles as matrix flow are caused
by local saturation due to low
The Hainich–Dün coniferous pre-drought profiles show low VD for all sprinkling amounts, especially in the topsoil between 4 and 22 cm (Fig. 8, center left panel). The 20 and 40 mm pre-drought profiles show unstained areas (no VD). The small VD values are even more pronounced in the drought profiles (Figs. 8 and 9, center left panels), in which all profiles exhibit unstained areas. Below the unstained layer, the VD rises to a maximum in 15 to 20 cm depth and drops again around 30 cm depth. The 20 and 60 mm drought profiles show throughout all depths low VD (Fig. 8, center left panel).
In all Hainich–Dün coniferous pre-drought profiles, no large SPWs occur and flow types are classified as macropore flow with low, mixed and high interaction. This applies also for the 20 and 60 mm drought profiles. In contrast, the 40 mm drought profile exhibits high SPWs between 15 and 25 cm depth. Therefore, the flow types in these depths are identified as matrix flow (Fig. 8, center left panel). The main flow types in the coniferous drought profiles are macropore flow with low, mixed and high interaction. The pre-drought profiles are dominated by macropore flow with low and mixed interaction. In both, pre-drought and drought profiles, the stone content is comparable (Fig. 8, center left panel).
The Hainich–Dün deciduous drought profiles exhibit smaller VD in the top 5 cm compared to the pre-drought profiles (Figs. 8 and 9, center right panels). Unstained areas are present in the top 5–10 cm of the 20 mm drought profile. The 40 and 60 mm pre-drought profiles show high VD values between 10 and 25 cm. High VD values are also present in the drought profiles, maintaining high values throughout the whole profile. While no high SPW values are found in the pre-drought profiles, high SPW values can be found in the 40 mm drought profile between 10 and 30 cm and in the 60 mm drought profile between 10 and 40 cm (Figs. 8 and 9, center right panels). The flow types of the deciduous pre-drought profiles are classified as macropore flow with low, mixed and high interaction. The drought profiles are also classified as macropore flow with low, mixed and high interaction and, where high SPW values occur, as matrix flow (homogeneous and heterogeneous) (Fig. 8, center right panel). The stone contents of the pre-drought and drought profiles increase with depth below 25–30 cm; the drought profiles exhibit a slightly higher stone content than the pre-drought profiles (Fig. 8).
The pre-drought pattern of the Schorfheide–Chorin site shows high SPWs and VDs in the top layers (0–10 cm depth) decreasing with depth (Fig. 8, bottom left panel). While the 20 and 40 mm pre-drought profiles show a maximum infiltration depth of about 45 and 30 cm, respectively, the 60 mm pre-drought profile is stained below 70 cm, exhibiting medium VD values (Fig. 8, bottom left panel). High SPW values are found in the 20 mm pre-drought profile up to a depth of 15 cm and in the 40 and 60 mm profiles up to 10 and 30 cm, respectively. The drought profiles of the coniferous plots show far lower VD values in the top layers compared to the pre-drought profiles (Figs. 8 and 9, bottom left panels). The 40 mm drought profile exhibits even an unstained layer at about 5 cm depth. High SPW values can be found in the 40 and 60 mm drought profile, not in the top layers, but between 20 and 25 cm depth (40 mm profile), and between 10 and 25 cm depth (60 mm profile) (Fig. 8, bottom left panel). This is reflected in the flow-type classification. Whereas matrix flow dominates the top layers in pre-drought profiles (at least the top 10 cm), matrix flow occurs below 10 cm depth in the 40 and 60 mm drought profiles (Fig. 8, bottom left panel).
The Schorfheide–Chorin deciduous pre-drought and drought patterns do not exhibit much differences in shape and in VD values in the 20 and 40 mm profiles (Fig. 8, bottom right panel). The largest differences in VDs can be found in the top 10 cm of the 20 mm profiles and in the 60 mm profile (Fig. 9, bottom right panel). In addition, the 20 mm drought profile exhibits an unstained layer around 40 cm depth (Fig. 8, bottom right panel). The difference between pre-drought and drought is more evident in the SPW values: high SPW values are found in the 40 and 60 mm pre-drought profiles in the top and bottom half of the profile; high SPW values are found in the drought profile in the bottom half, plus a small layer of 2 cm of high SPWs around 10 cm depth in the 40 mm profile (Fig. 8, bottom right panel). In the Schorfheide–Chorin deciduous pre-drought profiles, flow types of the 40 and 60 mm are dominated by matrix flow (Fig. 8, bottom right panel). However, all profiles in the pre-drought plots have a proportion of macropore flow. In the drought profiles, matrix flow only occurs in the bottom half of the 40 and 60 mm profiles.
To summarize, compared with pre-drought infiltration pattern, the drought pattern of all plots reveal differences in infiltration processes. For example, over 90 % of the depths ranges show significant differences in VDs between the drought and pre-drought site. Clayey and loamy soils behave similarly, developing unwetted soil layers. High SPW values in 20 to 30 cm depth of the drought pattern indicate local saturation. In sandy soils, the change from high SPW values of the pre-drought pattern to medium and low in the drought pattern exhibit a change from front-like to a more scattered infiltration. In general, the effects were more pronounced at the coniferous plots. These findings correspond well with the results of the WDPT tests: in the clayey and loamy soils (except Hainich–Dün deciduous plot), the unstained topsoil layers coincide with the high WDPTs (Fig. 5). Hainich–Dün coniferous plots stay unwetted up to a depth of about 15 to 20 cm and Schwäbische Alb plots to a depth of about 10 cm, which corresponds to the depths where the highest WDPT values were observed (Hainich–Dün: WDPT class 4; Schwäbische Alb: WDPT classes 4 and 3). In the sandy soils of the Schorfheide–Chorin profiles, low SPW values correspond to high WDPTs (class 2 and 3). Below the water repellent zone, SPW values increase again (Fig. 5, bottom panel).
The comparison of pre-drought infiltration patterns of the drought plots with patterns of the control plots (without drought treatment) showed broad agreements. All control plot profiles are comparable to the pre-drought plot profiles, including differences that can be addressed to small-scale heterogeneities of soil properties. When interpreting the patterns, the differences in VDs in the top layers of all plots need to be taken into account. When doing this, at all sites, the dye experiments before and during drought conditions can be directly compared.
In this study, it was hypothesized that the induced drought alters infiltration patterns due to changes in soil hydraulic properties (e.g., soil water repellency and forming of shrinkage cracks) and the main tree species has an effect on the magnitude of the response. The results of the infiltration experiment show clear evidence for changes in infiltration pattern as well as the importance of tree species on infiltration pattern: Schwäbische Alb plots have clayey soils with a high stone content, and show, in pre-drought and control plots, a slow and even infiltration. The drought-treated plots developed large areas with small volume densities and SPWs in the topsoil, while for deeper layers, broad stains (large SPWs) were observed, which cover the profiles for the most part (high VDs). This is typical for preferential flow that follows the shrinkage cracks of clayey soils or biopores of roots or soil fauna (Dekker and Ritsema, 2000; Hendrickx and Flury, 2001; Hardie et al., 2011). Water infiltrates quickly to deeper layers, bypassing a large proportion of the soil matrix. In deeper soil layers where the cracks or biopores end, local saturation occurs, and lateral redistribution into the soil matrix due to the now lower infiltration capacity and velocity can be observed. This also explains the similar pattern in the loamy Hainich–Dün soils.
A trend to more preferential flow was also observed in the Hainich–Dün plots, where the dense and loamy soils are also prone to shrinkage. Furthermore, in the Hainich–Dün drought profiles unstained (i.e., unwetted) areas in the topsoil layers were observed. This is more pronounced in the coniferous plot, where unstained areas were already visible in the pre-drought experiment. Preferential flow originates not only from cracks and biopores, but also from textual boundaries and instable wetting fronts (Doerr and Ritsema, 2006; Hendrickx and Flury, 2001). Unstable wetting fronts can occur due to air entrapment or hydrophobicity, which effectively hinders infiltration and redirects the water to structural and textural preferential flow paths (Doerr and Ritsema, 2006). The unwetted topsoil layers of the Hainich–Dün coniferous plots can be explained by the combination of severe water repellency and shrinkage cracks acting as effective bypasses.
In contrast to the other sites, Schorfheide–Chorin soils are sandy and highly permeable with low stone content. In both Schorfheide–Chorin plots, the infiltration patterns changed from a regular front-like stable infiltration to unstable, more scattered and fingered infiltration patterns. Following the conceptualization of unstable flow in water repellent soils by Ritsema et al. (1993, 2000), water flows, after entering the soil, through preferential pathways through the water repellent layer and distributes laterally in the divergence layer underneath. In fact, such flow patterns were pronounced in the Schorfheide–Chorin deciduous drought plot: medium VDs and SPWs up to a depth of 50–60 cm and larger SPWs in the layer beneath. This fits with the results of the WDPT tests, which show a slight to strong water repellency in the top 50–60 cm of the profile.
In general, drought induced major changes on the infiltration behavior of the examined soils. Clayey and loamy soils developed preferential flow. In these soils, the bypassing of the top 10–20 cm is fostered by water repellency, leading to unwetted topsoil layers. Sandy soils developed fingered infiltration patterns, due to the formation of a water repellent layer. In all three sites, the effects of the drought treatment were more pronounced in soils with coniferous main tree species than with deciduous main tree species.
In this study, it was hypothesized that the artificially induced drought will increase soil water repellency depending on the main tree species and soil properties. The highest water repellency was found in the coniferous plots of Schwäbische Alb and Hainich–Dün. Soils under coniferous trees often feature acidic soil conditions, which promote water repellency (Orfánus et al., 2014). In a study by Orfánus et al. (2014) liming practices and associated rise of pH values significantly reduced water repellency of former pine forest soils. Jost et al. (2004) explained the difference in recharge under a beech and a spruce forest stand, with higher hydrophobicity, and therefore the hindering of infiltration, combined with higher surface runoff of the spruce stand. This is in contrast to the findings of Buczko et al. (2006), who found the highest proportion of water repellent soils in pure beech stands compared to pure pine and mixed stands on sandy soils. However, in our study the sandy Schorfheide–Chorin plots showed higher mean WDPTs in the deciduous (beech) plot, than on the coniferous (pine) plot.
The soil texture can also influence the water repellency: a study of Gonzalez-Penaloza et al. (2013) suggests that water repellency is related to soil particle size. They induced water repellency by using different concentrations of stearic acid on samples of fine, medium and coarse sand. Water repellency was extreme in coarse-textured samples. The authors explained that by the lower specific surface compared to fine-textured samples and therefore smaller area that has to be covered by water repellent agents. We could not observe this effect in our sandy soils. The coarser-textured coniferous plot was less water repellent than the finer-textured deciduous Schorfheide–Chorin plot. However, water repellency can originate from a broad range of factors. The degree of water repellency of a soil also depends on the amount and type of organic matter that is incorporated in it (DeBano, 1981; Bisdom et al., 1993; Buczko et al., 2006; Vogelmann et al., 2013a, b), as well as the age and type of forest and litter type (Neris et al., 2013).
Several studies detected a significant impact of spruce litter on
infiltration processes, either by hydrophobicity (Schume et al., 2004) or
interception (Neris et al., 2013). Schume et al. (2004) found that spruce
litter can intercept up to 5 mm of precipitation and Neris et al. (2013)
found infiltration rates of 20 mm h
Furthermore, the plants of the forest understory can also influence hydrophobicity of the soil; plants are covered with a cuticle composed of hydrophobic liquids, embedded in a polyester matrix and wax crystalloids (Holloway, 1994; Barthlott and Neinhus, 1997). Water repellent plant coatings can be found in all plant life forms with a clear dominance among herbs (Neinhus and Barthlott, 1997; Dekker and Ritsema, 2000). It has even been discussed that hydrophobic exudates might be a strategy for plants, microorganisms and fungi, to suppress germination and growth of competing vegetation by reducing evaporation and nutrient leaching (Doerr et al., 2007).
Hydrophobicity is dependent on the moisture status of the soil, which is
defined by Doerr and Thomas (2003) as critical moisture or transition zone.
Vogelmann et al. (2013a) found a critical water threshold of 0.36 to
0.57 cm
In our rainfall exclusion experiment drought stress was not intense enough to induce mortality or strong changes in aboveground biomass of a particular species (Gimbel et al., 2015). Nevertheless, drought and water repellency may promote the die-off of fine roots, which thereupon contribute to the total organic matter in the soil. The amount of soil organic matter and its composition has a strong influence on the strength of water repellency (e.g., Vogelmann et al., 2013a; Bisdom et al., 1993; DeBano, 1981). Therefore, the die-off of fine roots may lead to a self-reinforcing circle of water repellency.
Soil organic matter can form micro- and macro-aggregates by acting as
binding agent between soil components (e.g., Tisdall and Oades, 1982; Annabi
et al., 2011) or by covering soil particles (e.g., Vogelmann et al., 2013a).
Vogelmann et al. (2013b) concluded in their study, that water repellency
leads to slower wetting of soil aggregates. Therefore, cohesive forces hold
up longer, which increases the resistance to disaggregation and, thus,
indirectly aiding to maintaining soil structure. Terrestrial fungi are also
in the focus of research concerning soil water repellency and aggregation
(e.g., Tisdall and Oades, 1982; Rillig and Mummey, 2006; Chau et al., 2012).
Zheng et al. (2014) found in three of nine species of ectomycorrhizal fungi
associated with
All of our experimental plots showed clear response to the drought treatment, irrespective of their soil type and vegetation cover. Especially the fast bypassing of the topsoil layer and the development of unstained and hence unwetted areas may bear consequences in the upcoming climate change. Sorption and degradation of contaminants is strongest in the topsoil and decrease with soil depth (Hendrickx and Flury, 2001). Thus, bypassing of the topsoil soil matrix foster early arrival times and high concentrations of contaminants in the groundwater, which was shown by several tracer field studies (e.g., Hendrickx and Flury, 2001; Ritsema et al., 1997; Hardie et al., 2011). Once formed, dry zones persist further wetting and additional water infiltrates through already existing preferential pathways, further stabilizing established flow paths (Dekker and Ritsema, 2000; Hagedorn and Bundt, 2002). Under present climate conditions, soil water repellency is already a widespread phenomenon (Buczko et al., 2006). For the predicted climate conditions, where droughts will be more common, an even higher level of hydrophobicity is to be expected, according to the findings of our rainfall reduction experiments.
Two years of rainfall reduction equivalent to an annual drought with a 40-year return interval was sufficient to change the soil properties and hence the infiltration pathways of six forest soils independent of soil type and tree species. All drought-treated soils, except one, developed slight to severe water repellency. Main tree species had a particular effect on hydrophobicity, but this only accounts for minor differences in infiltration pattern. The drought history or, more generally, the climatic condition in the past had more of an effect on the observed hydrophobicity and infiltration behavior than the actual antecedent soil moisture conditions of the soils. The results of this study suggest that drought effects on infiltration processes need to be considered in hydrological models to obtain realistic predictions regarding water quality and quantity in runoff and groundwater recharge.
The authors thank Emil Blattmann, Delon Wagner and Lukas Neuhaus for their technical assistance. Field support by Stefan Seeger, Frank Lehmann, Christoph Ries and Ruth Magh is gratefully acknowledged. The work has been funded by the DFG Priority Program 1374 “Infrastructure-Biodiversity-Exploratories” (DFG-Refno. WE4598/3-1, GE1090/11-1 and BR1698/16-1). We thank the managers of the three Biodiversity Exploratories frameworks, Swen Renner, Sonja Gockel, Andreas Hemp, Martin Gorke and Simone Pfeiffer for their work in maintaining the plot and project infrastructure, and Markus Fischer, the late Elisabeth Kalko, Eduard Linsenmair, Dominik Hessenmöller, Jens Nieschulze, Daniel Prati, Ingo Schöning, François Buscot, Ernst-Detlef Schulze and Wolfgang W. Weisser for their role in setting up the Biodiversity Exploratories project. Field work permits were issued by the responsible state environmental offices of Baden-Württemberg, Thüringen, and Brandenburg (according to § 72 BbgNatSchG). The article processing charge was funded by the German Research Foundation (DFG) and the Albert Ludwigs University Freiburg in the funding programme Open-Access Publishing. Edited by: M. Hrachowitz