Assessing Water Inﬁltration and Soil Water Repellency in Brazilian Atlantic Forest Soils

: 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 inﬁltration tests were performed, namely tension (Mini-Disk Inﬁltrometer, 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 di ﬀ erent water pressure heads allowed the estimation of the macroscopic capillary length in the ﬁeld. This approach, in conjunction with Beerkan measurements, allowed the design better estimates of the saturated soil hydraulic conductivity under challenging ﬁeld conditions, such as soil water repellency (SWR). This research also reports, for the ﬁrst time, evidence of SWR in the Atlantic Forest, which a ﬀ ected the early stage of the inﬁltration process with more frequency in the remnant forest. BEST-steady algorithm (%). The picture represents water droplets resisting inﬁltration into forest soil due to the water repellency.


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, 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, θ i (cm 3 cm −3 ), and the soil bulk density, ρ b (g cm −3 ). 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.

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].
Appl. Sci. 2020, 10, x FOR PEER REVIEW 3 of 14 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, θi (cm 3 cm −3 ), and the soil bulk density, ρb (g cm −3 ). 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.

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]. The saturated soil hydraulic conductivity, Ks (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, ts (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 ts value was determined as the first value for which: The saturated soil hydraulic conductivity, K s (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 s (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 s value was determined as the first value for which: where I reg (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 [31,32]. An illustrative example of t s 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 s (i.e., whenÊ > 2), while steady-state conditions establishes for all data points measured after time t s (i.e., whenÊ ≤ 2). 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 [31,32]. An illustrative example of ts 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 ts (i.e., when ̂ > 2), while steady-state conditions establishes for all data points measured after time ts (i.e., when ̂ < 2). 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, θs (cm 3 cm −3 ), considering the values of dry bulk density, ρb, previously determined.

Estimating the Saturated Soil Hydraulic Conductivity, Ks
We estimated Ks 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 Ks through an infiltration experiment of the Beerkan type [13] and an estimate of the macroscopic capillary length, λc (mm), expressing the relative importance of the capillary over gravity forces during water movement in unsaturated soil [33][34][35]. Firstly, the experimental steady-state infiltration rate, is (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]: where γ and γw are dimensionless constants [36,37] related to the infiltration front shape, that are commonly set at 0.75 and 1.818, and rd (mm) is the radius of the containment ring. Two different scenarios were considered to apply the SSBI method. The first scenario considered the MDI experiments, carried out with pressure heads of ℎ −20 = −20 mm and ℎ 0 = 0, to estimate λc by the following equation [38]: where ,ℎ −20 and ,ℎ 0 (mm 3 h -1 ) are the steady flow rates corresponding to ℎ −20 and ℎ 0 , respectively, and they were estimated as follows: For this scenario, we firstly averaged for each plot the individual is values, then plot-dependent λc values were estimated by Equation (3) (Table S1). 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, θ s (cm 3 cm −3 ), considering the values of dry bulk density, ρ b , previously determined.

Estimating the Saturated Soil Hydraulic Conductivity, K s
We estimated K s 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 s through an infiltration experiment of the Beerkan type [13] and an estimate of the macroscopic capillary length, λ c (mm), expressing the relative importance of the capillary over gravity forces during water movement in unsaturated soil [33][34][35]. Firstly, the experimental steady-state infiltration rate, i s (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, K sS (mm h −1 ) (the subscript S is used to indicate SSBI), as follows [33]: where γ and γ w are dimensionless constants [36,37] related to the infiltration front shape, that are commonly set at 0.75 and 1.818, and r d (mm) is the radius of the containment ring. Two different scenarios were considered to apply the SSBI method. The first scenario considered the MDI experiments, carried out with pressure heads of h −20 = −20 mm and h 0 = 0, to estimate λ c by the following equation [38]: where Q s,h −20 and Q s,h 0 (mm 3 h −1 ) are the steady flow rates corresponding to h −20 and h 0 , respectively, and they were estimated as follows: For this scenario, we firstly averaged for each plot the individual i s values, then plot-dependent λ c values were estimated by Equation (3) (Table S1).
The second K sS dataset was obtained considering λ c = 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, K sB (mm h −1 ) (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 K sB , by the following equation [41]: where b s (mm) is the intercept of the regression line fitted to the last data points of the I versus t plot.
The A (mm −1 ) 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]: 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 K sS and K sB in terms of factors of difference, FoD, calculated as the highest value between K sB and K sS divided by the lowest value between K sB and K sS . Differences between K sS and K sB not exceeding a factor of two were considered indicative of satisfactory K s predictions [33].

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].

Data Analysis
Following similar investigations [10,12,14], unique values of clay, silt, sand, OC, ρ b , θ i , and θ s 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 K sB , K sS , 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 K sS , K sB , K -20 , K -5 , K 0 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.

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). 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 ρb 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 ρb 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 θi ranged from 0.12 to 0.32 cm 3 cm −3 and the soil was significantly wetter in plots P1U, P2M, R4S, R4M, and F5I ( Figure S1).

Assessing SSBI Estimates
Both SSBI scenarios always yielded physically plausible estimates (i.e., positive Ks values). For the first scenario (i.e., λc 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 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 ρ b 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 ρ b 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 θ i ranged from 0.12 to 0.32 cm 3 cm −3 and the soil was significantly wetter in plots P1U, P2M, R4S, R4M, and F5I ( Figure S1).

Assessing SSBI Estimates
Both SSBI scenarios always yielded physically plausible estimates (i.e., positive K s values). For the first scenario (i.e., λ c estimated from multi tension experiments), the K sS 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., λ c = 83 mm), K sS 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 λ c values resulted in a slightly better estimation of K sS , yielding a lower mean FoD value, thus, only the first scenario was considered in the subsequent analysis.

Comparing BEST versus SSBI Estimates Under Soil Water Repellency Conditions
BEST-steady failed to estimate Ks in case of convex-shaped cumulative infiltration curves, which led to negative bs 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 Ks 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][53][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 Ks 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 and the intercept of regression line fitted to the last data points on the I vs. t plot. The magnitude of bs 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"

Comparing BEST versus SSBI Estimates Under Soil Water Repellency Conditions
BEST-steady failed to estimate K s in case of convex-shaped cumulative infiltration curves, which led to negative b s values and consequently to null K sB . The K sB 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 s 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][53][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 s 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 and the intercept of regression line fitted to the last data points on the I vs. t plot. The magnitude of b s 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 s 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). Lower SWR was detected in wetter soils. The correlation between θi 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. Ks 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 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 s values. Therefore, for all the experiments steady-state infiltration rates (i s ) were always reached before the end of the runs and after that the influence of hydrophobicity had ceased, so the K s 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.
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 s (mm) of regression line fitted to the last data points describing the steady-state conditions on the I vs. t plot, total duration, t end (s), total infiltrated depth, I end (mm), infiltrated depth at the equilibration time, I(t s ) (mm), equilibration time, t s (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 θ i 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 s data were positively correlated to ln(WDPT). This is logical, since both macropore flow (which affects the magnitude of K s ) 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, K sS 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, K sB and K sS 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.

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, K sS , 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 K sB 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 5 data, respectively. We also noticed a discrepancy between K sS and K 0 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). We also noticed a discrepancy between KsS and K0 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).

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 Ks 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 λc in the field. This approach, in conjunction with

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 s 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 λ c in the field. This approach, in conjunction with Beerkan measurements, allowed to generate better K s estimates based on field measurements, and also avoided any subjectivity caused by assuming a constant λ c value, which is often selected based on general descriptions of soil textural and structural characteristics when estimating K s from ponding infiltration experiments [34]. Nonetheless, developing alternative methods for estimating λ c is desirable for alleviating the amount of work necessary to accurately estimate K s .
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 −1 ), dry soil bulk density (ρ b in g cm −3 ), initial volumetric soil water content (θ i in cm 3 cm −3 ), and saturated volumetric soil water content (θ s in cm 3 cm −3 ), 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