**Assessing Water Infiltration and Soil Water Repellency in Brazilian Atlantic Forest Soils**

**Sergio Esteban Lozano-Baez 1,\* , Miguel Cooper <sup>2</sup> , Silvio Frosini de Barros Ferraz <sup>3</sup> , Ricardo Ribeiro Rodrigues <sup>1</sup> , Laurent Lassabatere <sup>4</sup> , Mirko Castellini <sup>5</sup> and Simone Di Prima 4,6**


Received: 22 February 2020; Accepted: 7 March 2020; Published: 12 March 2020

**Abstract:** This study presents the results of the soil hydraulic characterization performed under three land covers, namely pasture, 9-year-old restored forest, and remnant forest, in the Brazilian Atlantic Forest. Two types of infiltration tests were performed, namely tension (Mini-Disk Infiltrometer, MDI) and ponding (Beerkan) tests. MDI and Beerkan tests provided complementary information, highlighting a clear increase of the hydraulic conductivity, especially at the remnant forest plots, when moving from near-saturated to saturated conditions. In addition, measuring the unsaturated soil hydraulic conductivity with different water pressure heads allowed the estimation of the macroscopic capillary length in the field. This approach, in conjunction with Beerkan measurements, allowed the design better estimates of the saturated soil hydraulic conductivity under challenging field conditions, such as soil water repellency (SWR). This research also reports, for the first time, evidence of SWR in the Atlantic Forest, which affected the early stage of the infiltration process with more frequency in the remnant forest.

**Keywords:** Beerkan method; infiltration; forest restoration; soil water repellency

#### **1. Introduction**

The United Nations has declared the period 2021–2030 the decade of restoration to scale up existing initiatives, such as the Bonn Challenge, to restore degraded ecosystems [1]. It is expected that restoration will not only help to slow climate change through carbon sequestration, provide food, and increase biodiversity [2], but will also have hydrological benefits because of the perceived association between forest cover and soil hydrological ecosystem services [3,4]. In this context, it is necessary to better understand the consequences of forest regrowth on soil hydrological processes, such as water infiltration, which is fundamental to maintain productive soil-water-plant interactions, and also to control soil erosion and runoff, soil moisture content, and groundwater recharge in the ecosystems [5–7].

Estimating saturated and unsaturated soil hydraulic conductivities is crucial for interpreting and modeling soil hydrological processes. In addition, knowledge of these properties may provide information on the impact of land use on soils characteristics [8], which are rarely considered in studies of forest restoration [9]. During the last years, many infiltration methods and devices have been developed to determine soil hydraulic properties [10]. Among them, the Beerkan method [11] is becoming very popular in soil science because it constitutes a simple and an inexpensive way to determine the saturated soil hydraulic conductivity, *K<sup>s</sup>* , in the field [12,13]. On the other hand, the mini-disk infiltrometer (MDI) is a routinely used method for measuring infiltration rates under negative pressure head in the field. The MDI is easily transportable and easy to use on hillslopes, thus, it substantially facilitates the replicability of the measurements [10].

Our previous study in the Atlantic Forest of Brazil [14] used the Beerkan protocol at three land covers, namely pasture, 9-year-old restored forest, and remnant forest. Our results showed that water repellency impacted water infiltration, yielding convex shaped cumulative infiltration curves. However, this observation was not carefully assessed. Similarly, many studies on tropical soils have reported some indirect effects of water repellency on water infiltration, such phenomenon is still poorly documented, especially in comparison with temperate regions [5,15–17]. Soil water repellency (SWR) or hydrophobicity is a transient soil property with which soils increase the resistance to wetting and infiltration. It is spatially and temporally very variable [18,19]. This is caused mainly by amphiphilic molecules produced by plants and organism, and generally occurs after forest fires or dry periods. Other factors that can be related to water repellency are the soil texture, soil temperature, pH, water content, soil organic carbon, land use, and plant cover [18–21]. In addition, recent research highlighted that climate change could increase the water repellency of soils, due to the increasing occurrence of extreme events such as droughts, which create the soil conditions (i.e., high temperatures and low soil water content) that promote the water repellency [22].

Currently SWR is receiving increased attention in the scientific literature, due to the important hydrological effects. For example, SWR reduces infiltration capacity, increases runoff rates as well as leaching of agrochemicals and soil erosion, also it can affect negatively the crop production, nutrients, and plant-available water [17,23,24]. On the other hand, SWR has positive impacts on soil aggregate stability and organic carbon sequestration [22]. Müller and Deurer [17] reported the benefit of SWR for the arid and semi-arid climates, considering that this soil property reduces the loss of soil water by evaporation and allows the rainwater to reach deeper depths. Despite these efforts, our understanding of SWR is still limited [18], especially when subcritical phenomena occur.

This investigation aims to broaden our previous work [14], using the same location in the Brazilian Atlantic Forest. In particular, the specific objective was to compare both unsaturated and saturated soil hydraulic conductivity determined with simple and low-cost field infiltration methods (MDI and Beerkan), for three land covers, namely pasture, 9-year-old restored forest, and remnant forest. This paper includes the first measurements of SWR in the Brazilian Atlantic Forest and its relevance regarding soil hydraulic properties, which had never been investigated so far, to the best of our knowledge.

#### **2. Materials and Methods**

#### *2.1. Field Sites and Soil Sampling*

The study area (22◦53′ S, 46◦54′ W) is located in the county of Campinas, São Paulo State, Southeast Brazil. The area is located inside the sub-basin of Atibaia River (2800 km<sup>2</sup> ), which belongs to the Piracicaba River basin. The vegetation is classified as seasonal semideciduous forest. The zone is characterized by a complex geology located at the transition between the Atlantic Plateau and the Peripheral Depression geomorphological provinces, with Ultisols and Entisols as main soils [25]. The elevation varies from 600–900 m a.s.l. The climate is classified as Cwa according to the Köppen classification, with annual rainfalls of 1700 mm and mean annual temperature of 20 ◦C [26].

The three investigated land covers (pasture, P, restored forest, R, and remnant forest, F) correspond with those of Lozano-Baez et al. [14], with the use of the same 18 plots (7 × 7 m in size). These plots represent two pasture (P1 and P2), two restored forest (R3 and R4), and two remnant forest (F5 and F6). Each of these sites is further divided into three blocks (i.e., upslope, U, midslope, M, and downslope, D). For a detailed description of the field sites, the reader may refer to our previous work [13]. In brief, for a given plot, three undisturbed soil cores (5 cm in height and 5 cm in diameter) were collected at the 0–5 cm depth. With these samples we determined in the laboratory the initial volumetric soil water content, θ<sup>i</sup> (cm<sup>3</sup> cm−<sup>3</sup> ), and the soil bulk density, ρ<sup>b</sup> (g cm−<sup>3</sup> ). Three disturbed soil samples (0–10 cm depth) were also collected to determine the soil texture and the soil organic carbon content (OC). The soil texture was determined by the hydrometer method [27] and the OC was determined by the Walkley-Black method. – θ <sup>−</sup> ρ <sup>−</sup> –

#### *2.2. Unsaturated and Saturated Soil Hydraulic Conductivity Measurements*

Unsaturated soil hydraulic conductivity was measured using the Mini-Disc Infiltrometer (MDI) [28]. A total of 108 MDI experiments were carried out in the study sites. At each plot, we randomly selected six points with a minimum distance between measurements of 2 m. At the same sampling point, we used three different water pressure head values, *h*, in the ascending sequences −20, −5, and 0 mm, in order to sample several subdomains of the pore size distribution. Unsaturated soil hydraulic conductivity was calculated according to the method proposed by Zhang [29]. Before the MDI tests started, we removed the litter and leaves, and the sampled soil surface was gently levelled and smoothed. We used a thin layer of fine sand to ensure the contact between the infiltrometer and the soil (Figure 1). The thickness of the layer of sand was negligible and did not modify the imposed pressure head of the MDI. Visual readings of the water level were taken every 30 s until steady infiltration was nearly reached. For further descriptions of the MDI, details of measurements and calculations of hydraulic conductivity see Decagon Devices Inc. (Washington, DC, USA) [28]. − −

**Figure 1.** Picture in the forest site showing the mini-disk infiltrometer and the steel ring used for the Beerkan infiltration test.

− The saturated soil hydraulic conductivity, *K<sup>s</sup>* (mm h −1 ), was measured with ponding infiltration experiments of the Beerkan type [13]. At each plot, we performed seven Beerkan tests, for a total of 126 experiments. We used a steel ring with an inner diameter of 16 cm inserted to a depth of about 1 cm into the soil surface (Figure 1). In each infiltration point, a known volume of water (150 mL) was repeatedly poured into the cylinder at a small height above soil surface (i.e., a few cm) and the energy of the water was dissipated with the hand fingers to minimize the soil disturbance. Then, the time needed for each poured volume to complete infiltration was logged. This procedure was repeated until the difference in infiltration time between three consecutives trials became negligible.

() − () ()


 ̂ = |

The equilibration time, *t<sup>s</sup>* (s), namely the duration of the transient phase of the infiltration process, was estimated according to the suggested criterion by Bagarello et al. [30] for analyzing cumulative infiltration data. More specifically, the *t<sup>s</sup>* value was determined as the first value for which:

$$\hat{E} = \left| \frac{I(t) - I\_{reg}(t)}{I(t)} \right| \le E \tag{1}$$

where *Ireg*(*t*) is estimated from regression analysis considering the last points, and *E* defines a given threshold to check linearity. Equation (1) is applied from the end of the experiment until finding the first data point that fits the condition *E* ˆ ≤ *E* [31,32]. An illustrative example of *t<sup>s</sup>* estimation using the commonly used value of *E* = 2% [30] is shown in Figure 2a. Transient infiltration conditions therefore occur from time 0 until time *t<sup>s</sup>* (i.e., when *E* ˆ > 2), while steady-state conditions establishes for all data points measured after time *t<sup>s</sup>* (i.e., when *E* ˆ ≤ 2). ̂ ̂ θ <sup>−</sup> ρ

**Figure 2.** (**a**) Example of estimation of the equilibration time, *t<sup>s</sup>* (s), and infiltrated depth at the equilibration time, *I*(*ts*) (mm) from cumulative infiltration and (**b**) water repellency cessation time, WRCT (s), as the intersection point of two straight lines, representing the initial (hydrophobic) and the late (wettable) stages of the *I* vs. *t* 0.5 plot of a Beerkan infiltration run affected by soil water repellency (SWR).

At the end of each infiltration test, we collected a disturbed soil sample within the infiltration surface to determine the saturated gravimetric water content, and thus the saturated volumetric water content, θ<sup>s</sup> (cm<sup>3</sup> cm−<sup>3</sup> ), considering the values of dry bulk density, ρb, previously determined.

#### *2.3. Estimating the Saturated Soil Hydraulic Conductivity, Ks*

γ and γ

,ℎ−20

λ

,ℎ<sup>0</sup>

length, λ – − − = γγλ + 1 We estimated *K<sup>s</sup>* by the Simplified method based on the near Steady-state phase of a Beerkan Infiltration run (SSBI), recently proposed by Bagarello et al. [32]. This method estimates *K<sup>s</sup>* through an infiltration experiment of the Beerkan type [13] and an estimate of the macroscopic capillary length, λ<sup>c</sup> (mm), expressing the relative importance of the capillary over gravity forces during water movement in unsaturated soil [33–35]. Firstly, the experimental steady-state infiltration rate, *i<sup>s</sup>* (mm h −1 ), is estimated by linear regression analysis of the last data points of the cumulative infiltration, *I* (mm), versus time, *t* (h), plot, describing the near steady-state condition. Then, SSBI estimates the saturated soil hydraulic conductivity, *KsS* (mm h −1 ) (the subscript S is used to indicate SSBI), as follows [33]:

$$K\_{\rm sS} = \frac{\dot{i}\_{\rm s}}{\frac{\gamma \gamma \gamma\_w \lambda\_{\rm c}}{r\_d} + 1} \tag{2}$$

ℎ−20 ℎ<sup>0</sup>

ℎ−20 − ℎ<sup>0</sup> λ ℎ−20 − ℎ<sup>0</sup> where γ and γ<sup>w</sup> are dimensionless constants [36,37] related to the infiltration front shape, that are commonly set at 0.75 and 1.818, and *r<sup>d</sup>* (mm) is the radius of the containment ring. Two different scenarios

(,ℎ−20

,ℎ<sup>0</sup> ⁄ )

2

 = were considered to apply the SSBI method. The first scenario considered the MDI experiments, carried out with pressure heads of *h*−<sup>20</sup> = −20 mm and *h*<sup>0</sup> = 0, to estimate λ*<sup>c</sup>* by the following equation [38]:

$$
\lambda\_c = \frac{h\_{-20} - h\_0}{\ln \left( Q\_{s, h\_{-20}} / Q\_{s, h\_0} \right)} \tag{3}
$$

where *Qs*,*h*−<sup>20</sup> and *Qs*,*h*<sup>0</sup> (mm<sup>3</sup> h −1 ) are the steady flow rates corresponding to *h*−<sup>20</sup> and *h*0, respectively, and they were estimated as follows:

$$Q\_s = \mathbf{i}\_s \pi \mathbf{r}\_d^2 \tag{4}$$

For this scenario, we firstly averaged for each plot the individual *i<sup>s</sup>* values, then plot-dependent λ*<sup>c</sup>* values were estimated by Equation (3) (Table S1).

The second *KsS* dataset was obtained considering λ*<sup>c</sup>* = 83 mm, since it represents the suggested first approximation value for most soils types [37,39].

The Beerkan Estimation of Soil Transfer parameters (BEST) method [11] was also applied to estimate the saturated soil hydraulic conductivity, *KsB* (mm h−<sup>1</sup> ) (the subscript B is used to indicate BEST). More specifically, among the three existing BEST algorithms, we used the BEST-steady algorithm [40], that estimates *KsB*, by the following equation [41]:

$$K\_{\rm sB} = \frac{\mathbf{C} \text{ i}\_{\rm s}}{\mathbf{A} \,\mathbf{b}\_{\rm s} + \mathbf{C}} \tag{5}$$

where *b<sup>s</sup>* (mm) is the intercept of the regression line fitted to the last data points of the *I* versus *t* plot. The *A* (mm−<sup>1</sup> ) and *C* constants are defined for the specific case of the Brooks and Corey [42] relation and taking into account soil moisture initial conditions as follows [36]:

$$A = \frac{\gamma}{\mathbf{r}\_{\rm d}(\theta\_{\rm s} - \theta\_{\rm i})} \tag{6}$$

$$C = \frac{1}{2\left[1 - \left(\frac{\Theta\_i}{\Theta\_s}\right)^{\eta}\right](1-\beta)}\ln\left(\frac{1}{\beta}\right) \tag{7}$$

where β is a coefficient commonly set at 0.6, and η is a shape parameter that is estimated from the analysis of the particle size data with the pedotransfer function included in the BEST procedure [11].

Following Bagarello et al. [32], the BEST-steady algorithm was chosen to check the SSBI method, compering *KsS* and *KsB* in terms of factors of difference, FoD, calculated as the highest value between *KsB* and *KsS* divided by the lowest value between *KsB* and *KsS*. Differences between *KsS* and *KsB* not exceeding a factor of two were considered indicative of satisfactory *K<sup>s</sup>* predictions [33].

#### *2.4. Soil Water Repellency Carachterization*

Some of the Beerkan runs provided cumulative infiltration curves with convex shapes, signaling the occurrence of SWR phenomena [43,44].Then, the water repellency secession time, WRCT (s), was estimated from the intersection point of two straight lines, representing the initial and the late stages of *I* vs. *t* 0.5 relationship [45,46] (Figure 2b). The persistence of water repellency was measured using the water drop penetration time (WDPT) test. This test is widely used to determine the persistence of water repellency, it is easy to perform in field and presents the hydrological implications of hydrophobicity, because the amount of surface runoff is affected by the time required for the infiltration of droplets [47]. At each plot, we selected five sampling points. The WDPT was carry out by placing 10 drops (0.05 mL) of distilled water on to the soil surface and recording the time for their complete infiltration. Following other investigations [48,49] the infiltration recording was stopped after 3600 s. Moreover, if the drop did not infiltrate after this time interval, the value of 3600 s was assigned for the WDPT [47].

#### *2.5. Data Analysis*

Following similar investigations [10,12,14], unique values of clay, silt, sand, OC, ρb, θ<sup>i</sup> , and θ<sup>s</sup> were determined for each plot by averaging the measured values. For these soil parameters, we assumed a normal distribution, thus no transformation was performed on these data before statistical analysis. In addition, the *KsB*, *KsS*, *K–*20, *K*–5, *K*0, and WDPT data were assumed to be log-normally distributed since the statistical distribution of these data is generally log-normal [50]. Statistical comparison was conducted using two-tailed *t*-tests, whereas the Tukey Honestly Significant Difference test was applied to compare our data set. The ln-transformed *KsS*, *KsB*, *K–*20, *K–*5, *K*<sup>0</sup> and WDPT data were used for the statistical treatment. A probability level, α = 0.05, was used for all statistical analyses. It is reasonable to presume that infiltrometer data can also vary depending on the initial soil moisture and its effect on SWR [31], therefore the Spearman's rank correlation coefficients (*r*) were used to evaluate the relative influence of the soil properties on the infiltration process. For all the statistical analyses the Minitab© computer program (Minitab Inc., State College, PA, USA) was used.

#### **3. Results and Discussion**

#### *3.1. Soil Properties*

−

λ −

The 18 plots showed appreciable differences in soil texture. Sandy loam (i.e., P1M, P1D, R3U, R3M, R3D, and F5U) and sandy clay loam (i.e., P1U, R4S, R4M, R4D, F5M, and F5D) were the dominant soil textures among the plots, followed by clay loam (i.e., P2M, F6U, and F6M) and loam (i.e., P2U, P2D, and F6D). The soil texture of the plots is presented in the USDA textural triangle (Figure 3).

**Figure 3.** Textural distribution of the 18 plots in the USDA textural triangle.

– ρ *ρ* θ As pointed out by our previous study [14], the OC at the pasture sites in the soil depth 0–0.10 m was similar to remnant forest, while restored forest sites presented the lowest OC values. The highest ρ<sup>b</sup> values were observed in the restored forest R4, where the exposure of the soil and trampling pressure during the land-use history was greater in comparison with restored forest R3. Forest soils were characterized by the lowest ρ<sup>b</sup> values, which can be related to the heterogeneous soil structure and higher soil macroporosity in this cover [51,52]. At the time of sampling, the θ<sup>i</sup> ranged from 0.12 to 0.32 cm<sup>3</sup> cm−<sup>3</sup> and the soil was significantly wetter in plots P1U, P2M, R4S, R4M, and F5I (Figure S1).

λ

λ <sup>−</sup>

#### *3.2. Assessing SSBI Estimates*

Both SSBI scenarios always yielded physically plausible estimates (i.e., positive *K<sup>s</sup>* values). For the first scenario (i.e., λ<sup>c</sup> estimated from multi tension experiments), the *KsS* values ranged between 5.9 and 1486.8 mm h −1 . The mean FoD was equal to 1.36 (maximum value = 2.74) and the individual values were less than 2 and 1.5 for 89% and 78% of the cases, respectively (Figure 4). For the second scenario (i.e., λ<sup>c</sup> = 83 mm), *KsS* data ranged between 3.7 and 934.5 mm h −1 . The mean FoD was equal to 1.51 (maximum value = 2.37) and the individual values were less than 2 and 1.5 in the 90% and 53% of the cases, respectively. Therefore, using the estimated λ<sup>c</sup> values resulted in a slightly better estimation of *KsS*, yielding a lower mean FoD value, thus, only the first scenario was considered in the subsequent analysis.

λ λ − **Figure 4.** Empirical cumulative distribution function plot of the factors of difference between the saturated soil hydraulic conductivity values estimated by the BEST -steady (*KsB*) and SSBI methods (*KsS*). *KsS* data were estimated considering λ<sup>c</sup> = 83 mm (blue solid line) and the mean λ<sup>c</sup> values estimated for each sampled plot from the MDI experiments carried out with a suction of 0 and −20 mm (red dashed line).

#### *3.3. Comparing BEST versus SSBI Estimates Under Soil Water Repellency Conditions*

− – BEST-steady failed to estimate *K<sup>s</sup>* in case of convex-shaped cumulative infiltration curves, which led to negative *b<sup>s</sup>* values and consequently to null *KsB*. The *KsB* data ranged between 4.5 and 1394 mm h −1 (almost three orders of magnitude). The BEST-steady algorithm yielded physically plausible estimates (i.e., positive *K<sup>s</sup>* values) for 108 of 126 infiltration runs (i.e., 85.7% of cases). The percentage of successful runs was of 95.2% both for the pasture and restored forest (40 of 42 runs). With reference to the remnant forest, BEST led to a failure rate value of 33.3%, leading to lacks of estimates for 14 of 42 infiltration runs. In these cases, cumulative infiltration curves had convex shapes, which are typical for hydrophobia i.e., [42,45,46]. Such hydrophobia may result from significant amounts of organic matter content i.e., [52–54], originating from fauna and flora activities [55]. Soil texture also plays a major role on SWR, in particular, SWR is expected to increase for decreasing clay content. In this sense, our plots (i.e., F5U, R3U, R3D, P1M) with more sand content exhibited higher WDPT values. On the other hand, for the forest plots (i.e., F6U, F6M, F6D) the significant amounts of organic matter had a main role in generating relevant WDPT values also on finer textured soils [23,56].

The BEST-steady algorithm was unable to provide positive *K<sup>s</sup>* values, showing that BEST can only be used when the soil does not exhibit hydrophobic effect, as suggested by Lassabatere et al. [57]. As shown in Figure 5, at increasing failure rates of the BEST method corresponded higher WDPT values, suggesting that where hydrophobic condition occurred, mainly in the remnant forest plots, it was the main cause of failure of BEST-steady. More specifically, BEST-steady requires both the slope

"hockey like"

and the intercept of regression line fitted to the last data points on the *I* vs. *t* plot. The magnitude of b<sup>s</sup> depends on the entire cumulative infiltration curve (including the transient phase) [58,59], therefore that term is sensitive to SWR that impedes the early wetting phase of the infiltration process. When soil hydrophobicity occurred, the *I* vs. *t* 0.5 plot exhibited the characteristic "hockey-stick-like" shape [46], hiding the estimation of *K<sup>s</sup>* trough BEST-steady [42]. On the other hand, SSBI differs by the term expressing steady-state condition, considering exclusively the final infiltration rate [13]. The exclusive use of this term allowed to consider only the final stage of the infiltration process, i.e., when the hydrophobicity effect on infiltration was diminished. In this investigation eighteen Beerkan infiltration tests exhibited a clear hockey-stick-like shape, mainly at the remnant forest plots, that allowed calculation of WRCT as the intersection point of two straight lines, representing the initial and the late stages of *I* vs. *t* 0.5 relationships [45] (Figure 2b).

**Figure 5.** Comparison between the water drop penetration time, WDPT (s), and the failure rate of the BEST-steady algorithm (%). The picture represents water droplets resisting infiltration into forest soil due to the water repellency.

θ − Table 1 shows the results of the WRCT and the equilibration time calculations. Water repellency always affected the very early stage of the infiltration process since the WRCT values raged between 14 and 93 s, and they were always lower than the *t<sup>s</sup>* values. Therefore, for all the experiments steady-state infiltration rates (*is*) were always reached before the end of the runs and after that the influence of hydrophobicity had ceased, so the *K<sup>s</sup>* values estimated by the use of the SSBI method could be always properly estimated considering the last data points of the infiltration curves. Limiting the hydraulic characterization to the stabilized phase avoided the uncertainties due to specific shape of the cumulative infiltration and a no clear distinction between the early- and late-time infiltration process because soil hydrophobic phenomena [58]. In other words, the results presented in this study suggest that if hydrophobicity affects the first stage of a Beerkan infiltration test, the SSBI estimates should characterize the hydraulic property of the soil properly. We believe that this result has practical importance because the use the SSBI method allowed us to maintain the integrity of the dataset, and to compare the hydraulic behavior of different sites with different land uses, where soil hydrophobicity only occurs in some circumstance.

Spearman's rank correlation coefficients of the hydraulic measured properties. Moreover, maintaining a small water head on the soil surface may be useful to study the infiltration process in macroporous repellent soils [48]. The SSBI method, covering the soil surface with a practically null depth of water, hence, lower than the commonly water-entry values for repellent soils [59], could allow the operator to characterize water infiltration occurring through either structural or other gaps in the water repellent layer or as fingered flow through zones of hydrophilic or less water repellent soil [60]. On the contrary, establishing several cm of ponded head of water on the infiltration surface is expected to overwhelm SWR [61]. In this investigation, the detection of hockey-stick-like shapes

**θ** −

−

−

suggested that maintaining a small water head on the soil surface helped to prevent excessive positive pressure from overcoming SWR [62], and allowed the detection of water repellency.

**Table 1.** Values of the intercept, *b<sup>s</sup>* (mm) of regression line fitted to the last data points describing the steady-state conditions on the *I* vs. *t* plot, total duration, *tend* (s), total infiltrated depth, *Iend* (mm), infiltrated depth at the equilibration time, *I*(*ts*) (mm), equilibration time, *t<sup>s</sup>* (s), and water repellency cessation time, WRCT (s), for the eighteen Beerkan infiltration runs affected by hydrophobicity.


Lower SWR was detected in wetter soils. The correlation between θ<sup>i</sup> and ln(WDPT) was significant (*r* = −0.67, *p* = 0.002) (Table 2). This result was in line with the reasoning that the soil water content governs the interaction between soil particles and amphiphilic organic molecules, resulting from degradation of tree tissues, that coat soil particles and may be responsible for SWR [45]. The transition from wettable to hydrophobic status (and vice versa) is generally associated to a critical range of soil moisture [63]. The lower water content of this range defines the condition below which the medium is water repellent, the higher identifies the condition above which the medium is wettable. *K<sup>s</sup>* data were positively correlated to ln(WDPT). This is logical, since both macropore flow (which affects the magnitude of *Ks*) and water repellency phenomena were relevant at the remnant forest plots. In brief, the correlation between these two variables is not the result of a causal connection but the concomitancy of two processes: hydrophobia and macropore flow, which also lead to mainly subcritical water repellency. In addition, we conclude that hydrophobia had no effect on the estimation of the saturated hydraulic conductivity. Indeed, in opposite case, *KsS* and ln(WDPT) would have a negative correlation. Consequently, we assumed that the SSBI method proved efficient for detecting SWR and estimating properly the soil saturated hydraulic conductivity, at the same time. Lastly, *KsB* and *KsS* had a positive correlation with a value close to unity. The two estimators provide close estimates, as discussed above with the FoD. We then can conclude that soil hydrophobicity only affected the failure rate of the BEST-steady algorithm (Figure 5), without affecting the quality of its estimate when the method worked.

**Table 2.** Spearman's rank correlation coefficients of the hydraulic measured properties.


#### *3.4. Unsaturated versus Saturated Soil Hydraulic Conductivity*

The two types of infiltration tests, i.e., tension and ponding experiments, highlighted a clear increase of the hydraulic conductivity, especially at the remnant forest plots, when moving from near-saturated to saturated conditions (Figure 6). It is important to underscore that saturated and unsaturated conditions could be affected by the different soil texture. Our previous study [14] shows that soils with higher clay content (i.e., P2M, F6M, F6D) evidenced greater variations, by contrast, sandy soils (i.e., P1U, P1M, F5D) had lower variation. The mean values of the ratios between saturated, *KsS*, and near-saturated soil hydraulic conductivity, *K–*20, were 10.7, 21.5, and 118.3, for the pasture, restored forest and remnant forest, respectively. A similar trend was also detected when the *K*–5 values were considered, with the mean values of the ratios equal to 2.2, 5.6, and 23.7. Similar results were also obtained when *KsB* values were considered, with the values of the ratios equal to 10.6, 17.5, 92.0, and 2.2, 4.6, 17.4, for the *K*–20 and *K*<sup>5</sup> data, respectively. We also noticed a discrepancy between *KsS* and *K*<sup>0</sup> data, especially at the Forest site, because only under ponded conditions at the surface the macropores are activated [63]. The increase of the difference between saturated and unsaturated conditions can be explained by the activation of macroporosity at the forest plots [64]. Overall, the soil in the remnant forest is heterogeneous and characterized by a dominance of complex macropores. For example, a higher soil macroporosity and total porosity have been reported in the same forest soil by our previous work [14]. This soil macroporosity resulted from the better soil structure, which is caused by the high amount of biopores, roots, soil fauna activity and greater inputs of organic matter [52,63,64]. Moreover, soil variability at the scale of a few meters could have been less represented by the MDI, due to the small diameter of the infiltrometer (i.e., 4.5 cm). *–* – –

− − *– –* <sup>−</sup> − − **Figure 6.** Comparison of the mean saturated soil hydraulic conductivity values estimated with BEST-steady, *KsB* (mm h −1 ), and the SSBI method, *KsS* (mm h −1 ), and hydraulic conductivity, *K*<sup>0</sup> , *K–*<sup>5</sup> , and *K–*<sup>20</sup> (mm h −1 ), values measured with the minidisk infiltrometer under a tension of 0, −5, and −20 mm. For a given plot, means that do not share a letter are significantly different according to the Tukey honestly significant difference test (*p* < 0.05).

#### **4. Conclusions**

To improve the soil hydraulic characterization of the soils under different land uses in the Atlantic Forest of Brazil, we measured and compared the unsaturated and saturated soil hydraulic conductivity, using the MDI and Beerkan method for three land covers, namely pasture, 9-year-old restored forest, and remnant forest. This research reports, for the first time, provide evidence of SWR in the Atlantic Forest, especially in the remnant forest. Our measurements demonstrated that SWR affected the early stage of the infiltration process. The comparison between alternative methods to estimate *K<sup>s</sup>* allowed to account for the effect of SWR on water infiltration measurements. In particular, when there are evidences of SWR, our results suggest using the SSBI method instead of BEST-steady to avoid the failure of the analysis in case of string SWR. Indeed, the SSBI method allowed to maintain the integrity of the infiltration dataset, facilitating the hydraulic comparison between different land uses. Tension

λ

(MDI) and ponding (Beerkan) infiltration tests provided a complementary information, highlighting a clear increase of the hydraulic conductivity, especially at the remnant forest plots, when moving from near-saturated to saturated conditions. This information is relevant to assess the infiltration recovery after forest restoration, as it signals soil structure heterogeneity and higher soil macroporosity. In addition, measuring the unsaturated soil hydraulic conductivity with different water pressure heads also allowed to estimate λ<sup>c</sup> in the field. This approach, in conjunction with Beerkan measurements, allowed to generate better *K<sup>s</sup>* estimates based on field measurements, and also avoided any subjectivity caused by assuming a constant λ<sup>c</sup> value, which is often selected based on general descriptions of soil textural and structural characteristics when estimating *K<sup>s</sup>* from ponding infiltration experiments [34]. Nonetheless, developing alternative methods for estimating λ<sup>c</sup> is desirable for alleviating the amount of work necessary to accurately estimate *K<sup>s</sup>* .

In this investigation we used the water repellency secession time (WRCT) and water drop penetration time (WDPT) to assess SWR. The SWR was observed in pasture and forest soils with higher sand content. As expected, SWR phenomena were less severe for increasing soil moisture conditions and more common on remnant forest soils. SWR has important hydrological effects, including water supply in forest ecosystems. Thus, SWR cannot be neglected in forest soil hydraulic studies and it must be accounted when developing hydrological models. Future research will focus on understanding the interactions between vegetation, soil biology and soil properties (i.e., including physical, chemical, and mineralogical properties) that are promoting SWR in Atlantic forest soils. Details on the effects of forest restoration on water repellency are severely lacking. It is also important to consider the temporal and spatial dynamic of the soil infiltration and water repellency. For example, future studies could quantify the spatial extent at the larger scale, upscaling the measurement from point to catchment scale.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2076-3417/10/6/1950/s1, Figure S1: Comparison between the mean soil organic carbon content (OC in g KgPress release−<sup>1</sup> ), dry soil bulk density (ρ<sup>b</sup> in g cm−<sup>3</sup> ), initial volumetric soil water content (θ<sup>i</sup> in cm<sup>3</sup> cm−<sup>3</sup> ), and saturated volumetric soil water content (θ<sup>s</sup> in cm<sup>3</sup> cm−<sup>3</sup> ), values for the 18 sampled plots. Bars indicate standard deviation. For a given variable and plot, means that do not share a letter are significantly different according to the Tukey honestly significant difference test (*P* < 0.05). The subscript letter refers to the landscape position (Upslope, Middleslope and Downslope) in each site, Table S1: Mean values for each sampled plot (P1, P2, R3, R4, F5 and F6) of the steady-state infiltration rates, *is*,*h*<sup>0</sup> and *<sup>i</sup>s*,*h*−<sup>20</sup> (mm h−<sup>1</sup> ), flow rates, *Qs*,*h*<sup>0</sup> and *Qs*,*h*−<sup>20</sup> (mm<sup>3</sup> h −1 ), obtained from the MDI experiments carried out with a pressure head *h*<sup>0</sup> = 0 and *h*−<sup>20</sup> = −20 mm, and macroscopic capillary length, λ*<sup>c</sup>* (mm), estimated by Equation (3).

**Author Contributions:** S.E.L.-B. carried out the data collection and wrote the initial draft; M.C., S.F.d.B.F. and R.R.R. conceived and designed the experiment; L.L., M.C. and S.D.P. revised, discussed, modified, and supplemented the ideas for the final draft. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was supported by the Fundação de Amparo à Pesquisa do Estado de São Paulo (BIOTA/FAPESP Program: 2013/50718-5 and 1999/09635-0), and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq 561897/2010-7 and Miguel Cooper's, Ricardo Ribeiro Rodrigues's and Silvio Frosini de Barros Ferraz's scientific productivity fellowships). This work was also supported through the "Programma Operativo Nazionale (PON) Ricerca e Innovazione 2014—2020 (Linea 1—Mobilità dei ricercatori, AIM1853149, CUP: J54I18000120001) funded by the European Regional Development Fund (ERDF) and the Italian Ministry of Education, University and Research (MIUR).

**Acknowledgments:** We are grateful to Luiz Felippe Salemi who gave their valuable suggestions. S.E.L.-B. would like to thank Julia Gardies, Daigard Ricardo Ortega, Monica Borda and Miller Ruiz for their support during field work.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

*Article*
