Performance of Conservation Techniques for Semiarid Environments : Field Observations with Caatinga , Mulch , and Cactus Forage Palma

Understanding small-scale hydrologic processes and the impact of soil conservation techniques are crucial in reducing runoff and sediment losses in semi-arid regions. This study was conducted in the Alto Ipanema River Basin, in Pernambuco State (Brazil). Soil and water dynamics were intensely monitored in twelve experimental plots with different coverage conditions (plot with bare soil—Bare; plot with natural vegetation—Natur; plot with mulch—Mulch; plot with Cactus Palma—Palma). By far, bare soil conditions produced higher runoff and soil losses. Mulch cover was close to natural vegetation cover, but still presented higher runoff and sediment losses. Palma, which is a very popular spineless cactus for animal feed in the Brazilian semi-arid region, presented an intermediate hydrologic impact in controlling runoff, enhancing soil moisture, and also reducing soil losses. Experiments were conducted in one hydrologic year (2016/2017) at three different sites. They were intensely monitored and had the same number of plots. This enabled us to carry out a robust performance assessment of the two soil conservation practices adopted (Mulch and Palma), compared to natural vegetation cover and bare soil conditions. Such low-cost alternatives could be easily adopted by local farms in the region, and, hence, improve soil reclamation and regional resiliency in a water-scarce environment.


Introduction
Erosion is a natural process, leading to serious environmental consequences, reducing agricultural productivity, and increasing the sediment amount for the water bodies downstream.Soil erosion is probably one of the most relevant environmental degradation challenges, especially in hilly shallow soils typical in arid and semi-arid regions [1][2][3][4].
Reliable local-scale measurements of runoff and erosion under natural rainfall for different soil cover conditions are still limited, especially in semiarid environments, with sparse vegetation cover [5,6].
Mulch minimizes the impact rain has on soil surfaces [7], which might contribute to improved soil fertility, increase water availability through enhancing infiltration, and reduce evaporation, thus minimizing nutrient losses and also controlling soil temperature variations [4,8,9].Another potential alternative largely adopted in Brazilian semiarid areas is the cultivation of forage Palma, a cactus crop that presents a modified stem a with thin, flat vertical structure.Although it has limited rainfall interception, it has a high leaf area index [10], and if cropped along contour lines, then it tends to reduce surface runoff [11].
The use of field erosion plots under natural rainfall allows the study of many hydrological processes at on local scale.Well-designed experimental plots provide an opportunity to improve understanding of the local-scale water budget and plant cover effects on runoff.The knowledge to be acquired on these experimental small plots is valuable for some specific issues, such as evaluation of soil and water conservation practices, water retention, infiltration, and soil water dynamics [12].
Small-scale plots both in laboratory and in natural watersheds can also successfully contribute to the understanding of major hydrological processes during extreme rainfall events [13,14].In semiarid environments it is very challenging to collect runoff and soil loss data of a reasonable quality under natural conditions, as the number of rainfall runoff events is very limited; therefore, each event might have an important role in evaluating the impact of soil and water conservation alternatives [15].
Santos and Montenegro [16] analyzed high intensity rainfall events and their contributions to soil disaggregation, transport, and deposition in the Pernambuco semiarid environment.Based on a time series of 29 years, the authors verified that the first half of the year was characterized by rainfall events with high erosive potential, and complex rainfall patterns were observed with a higher occurrence where the peak level fell at the beginning of each rainfall event.Previously, Santos et al. [11] verified the potential use of mulching and Palma cactus in increasing soil moisture, over a 301-day study period, during a wet year.This was conducted in the same area as in our investigation, although no attempt was made regarding rainfall runoff analysis.
In studies of water and soil conservation in watersheds in Africa, Wenninger et al. [17] highlighted that hydrological processes in semi-arid regions usually presented high spatio-temporal variability.Hence, it is important to obtain local rainfall and runoff measurements, for soil moisture dynamics, in order to properly understand runoff generation.
Northeastern Brazil, particularly the semi-arid region of Pernambuco State, is usually subject to high intensity local rainfall events, known as thunderstorms [18].Such events cause high runoff rates and sediment losses, requiring conservation alternatives to be adopted in order to prevent irreversible degradation of the topsoil.
The objective of this study was to investigate small-scale hydrologic processes and the impact of soil conservation techniques on reducing runoff and sediment losses.This study used textural characterization of the associated transported sediments, using runoff plots with different soil covers (Bare, Natur, Mulch, and Palma) under natural rainfall, in the Caatinga biome of Brazil.

Study Sites
The study area (Figure 1) was located at the Alto Ipanema River Basin (AIRB), a sub-basin of Ipanema River, one of the sub-basins monitored by the Hydrology Network of the Semi-Arid Region (REHISA).It is part of the municipalities of Arcoverde and Pesqueira, in the Pernambuco State (Brazil).The area has a complex landscape, characterized by a high spatial variability of elevation and climate, and deciduous vegetation, which constitutes the Caatinga Biome.Previous field investigations were adopted [10,16] as sources of information for the physical and chemical characteristics of soil in the experimental plots.These studies were located in the same pedological unit, and the soil properties were very similar among the three studied sites.The mean values from the data are presented in Table 1.The soil classification was Ultisol Eutrophic Typical and the infiltration capacity was 0.134 m h −1 .Four experimental plots with different cover conditions were considered at each of the three sites (corresponding to 12 plots in total).Soil and water monitoring was performed in experimental plots with different cover conditions (plot with bare soil-Bare; plot with natural cover-Natur; plot with mulch-Mulch; plot with Cactus Palma-Palma).Maintenance was carried out on the experimental plots before the beginning of the experiment, and minimal reworking was performed only to allow proper representativeness of the soil cover conditions (e.g., weed control procedures).
Figure 2 shows a general view of the investigated cover conditions (Bare, Natur, Mulch, and Palma): Bare-soil without any of natural or artificial cover on the plot; Natur-predominantly natural and/or spontaneous vegetation composed of small and medium-sized shallow caatinga, with predominant quince (Croton sonderianus) and jurema-preta (Mimosa hostilis Benth.);Mulch-dry grass The vegetation presents strong seasonality over time, characteristic of the Caatinga Biome.In the dry period, the native forest area presents loss of foliage (deciduous behaviour) and, in areas of sparse vegetation, large areas exhibit exposed soils.However, the cover conditions change considerably during the rainy season, which is characteristic of semi-arid regions, with rapid foliage regeneration [6,19].
Previous field investigations were adopted [10,16] as sources of information for the physical and chemical characteristics of soil in the experimental plots.These studies were located in the same pedological unit, and the soil properties were very similar among the three studied sites.The mean values from the data are presented in Table 1.The soil classification was Ultisol Eutrophic Typical and the infiltration capacity was 0.134 m h −1 .Four experimental plots with different cover conditions were considered at each of the three sites (corresponding to 12 plots in total).Soil and water monitoring was performed in experimental plots with different cover conditions (plot with bare soil-Bare; plot with natural cover-Natur; plot with mulch-Mulch; plot with Cactus Palma-Palma).Maintenance was carried out on the experimental plots before the beginning of the experiment, and minimal reworking was performed only to allow proper representativeness of the soil cover conditions (e.g., weed control procedures).
Figure 2 shows a general view of the investigated cover conditions (Bare, Natur, Mulch, and Palma): Bare-soil without any of natural or artificial cover on the plot; Natur-predominantly natural and/or spontaneous vegetation composed of small and medium-sized shallow caatinga,

Climate Data
According to Köppen's classification, the climate in the region was BSsh (extremely hot, semiarid).The mean annual precipitation, according to local historical data (1910 to 2012) was 671.90 mm.According to climate classification, total annual evaporation was ~1600 mm, and the mean temperature was 27.40 °C, varying from 18.5 to 29 °C.According to the aridity index, on average more than 60% of the region presented low susceptibility to desertification and exhibited medium susceptibility spots [20].The predominant vegetation was the hypoxerophilic caatinga and cactus [21].
Climatic data were recorded at an Automatic Weather Station installed near the experimental plots (Pesqueira, Brazil), as observed on the map (Figure 1).The station consisted of a set of sensors and a communication interface for data recording (CR1000 datalogger) and transfer including: an anemometer, a rain gauge, a temperature sensor, a relative humidity sensor, and a pyranometer, recording data hourly.

Data Processing
From October 2016 to October 2017, runoff and sediment yield were monitored and characterized in the 12 experimental plots.The experimental plots were already established some years ago [11], and they were delimited by brick walls 0.25 m above the soil surface and inserted 0.10 m into the soil.Downstream of the plots, a drain system collected runoff into two consecutive tanks 1 m³ in volume (Figure 3).During the study period, no runoff events reached the maximum storage capacity of the collection tanks, nor were significant losses of the stored water due to evaporation observed.Snapshots were taken by the end of each rainfall event using a camera.

Climate Data
According to Köppen's classification, the climate in the region was BSsh (extremely hot, semi-arid).The mean annual precipitation, according to local historical data (1910 to 2012) was 671.90 mm.According to climate classification, total annual evaporation was ~1600 mm, and the mean temperature was 27.40 • C, varying from 18.5 to 29 • C. According to the aridity index, on average more than 60% of the region presented low susceptibility to desertification and exhibited medium susceptibility spots [20].The predominant vegetation was the hypoxerophilic caatinga and cactus [21].
Climatic data were recorded at an Automatic Weather Station installed near the experimental plots (Pesqueira, Brazil), as observed on the map (Figure 1).The station consisted of a set of sensors and a communication interface for data recording (CR1000 datalogger) and transfer including: an anemometer, a rain gauge, a temperature sensor, a relative humidity sensor, and a pyranometer, recording data hourly.

Data Processing
From October 2016 to October 2017, runoff and sediment yield were monitored and characterized in the 12 experimental plots.The experimental plots were already established some years ago [11], and they were delimited by brick walls 0.25 m above the soil surface and inserted 0.10 m into the soil.Downstream of the plots, a drain system collected runoff into two consecutive tanks 1 m 3 in volume (Figure 3).During the study period, no runoff events reached the maximum storage capacity of the collection tanks, nor were significant losses of the stored water due to evaporation observed.Snapshots were taken by the end of each rainfall event using a camera.Each plot had eight polyvinyl chloride (PVC) access tubes installed 2 m apart from each other.Due to the depth of the impervious layer of the soil profile, the access tubes were installed mainly up to 0.60 m in depth, and the depths for the soil water content measurement were 0.20, 0.40, and 0.60 m.The soil water content was monitored every 15 days using a Neutron Probe Device model CPN 503 DR.The neutron probe readings were converted into soil water volumetric content based on preestablished calibration curves that were available for the studied soil.

Water and Sediment Measurements
For all effective rainfall events, runoff and sediments were manually collected individually in the tanks within one day after each event.In addition, the particle size of eroded soil was monitored during the experiment for different soil covers in order to evaluate the impact of runoff on soil texture.
Collection of the stored water only occurred in the first tank, because there was no runoff volume that exceeded the tank's capacity.The stored water was stirred for uniformity (for 3 min), and 1 L random samples were collected in triplicate.In order to estimate the sediment concentration, samples were dried in an oven (105 °C).Soil sand, clay, and silt fractions were determined by the Boyoucus densimeter method, following the methodology proposed by Embrapa [22].
Total rainfall and rainfall intensity were measured with an automatic rain gauge (TB4-L) located at each experimental site, with an average distance of 2 m from the plots.This equipment was connected to a CR1000 datalogger, programmed to record events every 5 min.

Statistical Analysis
The experimental plots within the sites (Malaquias, Edivaldo, and João) were randomly located.The experimental data were submitted to variance analysis.For comparison of the mean behavior for the treatments, Tukey's test was applied at a 5% significance level.Correlation analysis was carried out, scored by the R 2 value and significance of the slope coefficient, when appropriate.Each plot had eight polyvinyl chloride (PVC) access tubes installed 2 m apart from each other.Due to the depth of the impervious layer of the soil profile, the access tubes were installed mainly up to 0.60 m in depth, and the depths for the soil water content measurement were 0.20, 0.40, and 0.60 m.The soil water content was monitored every 15 days using a Neutron Probe Device model CPN 503 DR.The neutron probe readings were converted into soil water volumetric content based on pre-established calibration curves that were available for the studied soil.

Water and Sediment Measurements
For all effective rainfall events, runoff and sediments were manually collected individually in the tanks within one day after each event.In addition, the particle size of eroded soil was monitored during the experiment for different soil covers in order to evaluate the impact of runoff on soil texture.
Collection of the stored water only occurred in the first tank, because there was no runoff volume that exceeded the tank's capacity.The stored water was stirred for uniformity (for 3 min), and 1 L random samples were collected in triplicate.In order to estimate the sediment concentration, samples were dried in an oven (105 • C).Soil sand, clay, and silt fractions were determined by the Boyoucus densimeter method, following the methodology proposed by Embrapa [22].
Total rainfall and rainfall intensity were measured with an automatic rain gauge (TB4-L) located at each experimental site, with an average distance of 2 m from the plots.This equipment was connected to a CR1000 datalogger, programmed to record events every 5 min.

Statistical Analysis
The experimental plots within the sites (Malaquias, Edivaldo, and João) were randomly located.The experimental data were submitted to variance analysis.For comparison of the mean behavior for the treatments, Tukey's test was applied at a 5% significance level.Correlation analysis was carried out, scored by the R 2 value and significance of the slope coefficient, when appropriate.

Rainfall Event Analyses
The studied hydrologic year of 2016/2017 can be classified as water-scarce.During the study period there were 76 daily rainfall events, which totaled 404.20 mm (corresponding to 60% of the mean rainfall depth in the 2017 hydrological year).
Figure 4 presents the observed rainfall, runoff, soil losses, and soil moisture temporal behaviors occurring from October 2016 to October 2017, for rainfall higher than 5 mm.

Rainfall Event Analyses
The studied hydrologic year of 2016/2017 can be classified as water-scarce.During the study period there were 76 daily rainfall events, which totaled 404.20 mm (corresponding to 60% of the mean rainfall depth in the 2017 hydrological year).
Figure 4 presents the observed rainfall, runoff, soil losses, and soil moisture temporal behaviors occurring from October 2016 to October 2017, for rainfall higher than 5 mm.
Only seven events produced runoff.The highest precipitation intensity occurred on 4 March 2017, with a value of 70 mm during 0.5 hours, corresponding to a return period of 300 years, according to Silva et al. [18].
The water and sediment production characteristics for each rainfall event are presented in Figure 4, as well as the mean soil moisture for the 0 to 0.6 m layer, as a function of time.
Runoff generation and soil losses were very limited for the Natur soil cover condition, and the conservation practices of Mulch and Palma efficiently reduced water and soil losses, as expected.More detailed rainfall runoff and rainfall soil loss analyses have been carried out in Section 3.2, exploring the correlations between the variables and also their statistical behaviors.The main variations observed for the moisture values could be explained by semi-arid characteristics, with precipitation events of high intensity and small duration, followed by high evapotranspiration rates.Lopes et al. [23], studying the rainfall spatial distribution for a region of the Only seven events produced runoff.The highest precipitation intensity occurred on 4 March 2017, with a value of 70 mm during 0.5 hours, corresponding to a return period of 300 years, according to Silva et al. [18].
The water and sediment production characteristics for each rainfall event are presented in Figure 4, as well as the mean soil moisture for the 0 to 0.6 m layer, as a function of time.
Runoff generation and soil losses were very limited for the Natur soil cover condition, and the conservation practices of Mulch and Palma efficiently reduced water and soil losses, as expected.More detailed rainfall runoff and rainfall soil loss analyses have been carried out in Section 3.2, exploring the correlations between the variables and also their statistical behaviors.
The main variations observed for the moisture values could be explained by semi-arid characteristics, with precipitation events of high intensity and small duration, followed by high Water 2019, 11, 792 7 of 15 evapotranspiration rates.Lopes et al. [23], studying the rainfall spatial distribution for a region of the Brazilian semi-arid region, observed that there was a high spatio-temporal variation caused by the previously mentioned rainfall characteristics.
The different surface covers influenced soil moisture dynamics.At the beginning of the hydrological year, soil moisture was below the soil water content at the permanent wilting point (0.071 m 3 m −3 ).Soil moisture for bare soil was the lowest among the studied soil cover conditions.
Table 2 gives an overview of the different sites and plots, namely the relative position of the different cover plots in each site.Although there were differences, plot characteristics were approximately the same, and rainfall had the same magnitude order.In addition, Table 2 also presents the statistical analysis results for runoff, soil losses, and soil moisture for the different sites and soil cover conditions.Statistical differences occurred between the distinct cover conditions and bare soil for the three studied sites (Malaquias, Edivaldo and João-Table 2).Higher soil loss was observed for the bare soil condition, followed by Palma.Natural Caatinga cover and mulch presented lower soil erosions.
Runoff had different behavior, with natural cover generating the smallest runoff depth and bare soil conditions producing the highest runoff.Mulch reduced soil loss because small barriers of accumulated mulch material formed, approximately following the elevation contours.These barriers (Figure 5) promoted the deposition of soil particles upslope.Statistical differences were also detected for the mean soil moisture for the whole period, being more evident for the João Site.Araújo et al. [24] studied the spatial distribution of soil moisture for the Ipanema River Basin, where the experimental plots were inserted, and also observed low water availability conditions for the whole area in the same period.
Observing Figure 4, it can be verified that some moisture differences were more evident for periods with rainfall spells as a result of the rainfall pattern, evapotranspiration, and the limited soil water holding capacity at the plots.
According to Table 2, mulch was the most suitable conservation practice for maintaining the highest soil moisture values.The observed soil moisture contents varied over time during the experimental period, and variations were related to the different cover types.
Santos et al. [11], studying the same soil cover conditions for a wet hydrologic year in the same region, observed that soil surface conditions had a high influence on soil moisture variation, both in the dry and the rainy periods.It was verified that natural vegetation cover (Caatinga biome) presented the largest water content in soil compared to the other treatments for the entire rainy period.Bezerra et al. [25] highlighted the importance of monitoring soil water dynamics for different cover conditions in semi-arid Brazil, especially in the Caatinga biome, aiming to improve soil water storage and to provide experimental in situ data for soil losses resulting from rainfall events.

Runoff and Soil Loss Correlations
Figures 6 and 7 verify the relationships between the runoff coefficient (defined as the ratio between runoff depth and total rainfall) and soil loss as a function of rainfall depth and intensity.Rainfall intensity explains better the observed variations in the runoff coefficient than rainfall depth.It was shown that mulching successfully protected the soil surface for all three sites, which allowed a higher infiltration rate and, thus, greater soil moisture storage.For Palma, a similar behavior was observed both for runoff and soil loss.
These results can be explained by the following mechanisms and surface characteristics, presented by Montenegro et al. [4]: (1) soil cover protection from direct impact of rain drops; (2) higher hydraulic roughness based on mulch cover, retarding surface flow, and enhancing infiltration; Statistical differences were also detected for the mean soil moisture for the whole period, being more evident for the João Site.Araújo et al. [24] studied the spatial distribution of soil moisture for the Ipanema River Basin, where the experimental plots were inserted, and also observed low water availability conditions for the whole area in the same period.
Observing Figure 4, it can be verified that some moisture differences were more evident for periods with rainfall spells as a result of the rainfall pattern, evapotranspiration, and the limited soil water holding capacity at the plots.
According to Table 2, mulch was the most suitable conservation practice for maintaining the highest soil moisture values.The observed soil moisture contents varied over time during the experimental period, and variations were related to the different cover types.
Santos et al. [11], studying the same soil cover conditions for a wet hydrologic year in the same region, observed that soil surface conditions had a high influence on soil moisture variation, both in the dry and the rainy periods.It was verified that natural vegetation cover (Caatinga biome) presented the largest water content in soil compared to the other treatments for the entire rainy period.Bezerra et al. [25] highlighted the importance of monitoring soil water dynamics for different cover conditions in semi-arid Brazil, especially in the Caatinga biome, aiming to improve soil water storage and to provide experimental in situ data for soil losses resulting from rainfall events.

Runoff and Soil Loss Correlations
Figures 6 and 7 verify the relationships between the runoff coefficient (defined as the ratio between runoff depth and total rainfall) and soil loss as a function of rainfall depth and intensity.Rainfall intensity explains better the observed variations in the runoff coefficient than rainfall depth.It was shown that mulching successfully protected the soil surface for all three sites, which allowed a higher infiltration rate and, thus, greater soil moisture storage.For Palma, a similar behavior was observed both for runoff and soil loss.
These results can be explained by the following mechanisms and surface characteristics, presented by Montenegro et al. [4]: (1) soil cover protection from direct impact of rain drops; (2) higher hydraulic roughness based on mulch cover, retarding surface flow, and enhancing infiltration; and (3) water retention of the mulch cover.In addition, contour barriers provided by the Palma reduced overland flow velocities, favoring sediment deposition.Runoff generation was strongly related to antecedent moisture conditions of a specific event, as observed in Figure 8a.Soil loss, however, was more correlated to rainfall intensity (Figure 8b).Characteristic rainfall temporal regimes in the Brazilian semi-arid region, where higher intensities were verified at the beginning of the event, were also observed.
Water 2019, 11, x FOR PEER REVIEW 10 of 16 Runoff generation was strongly related to antecedent moisture conditions of a specific event, as observed in Figure 8a.Soil loss, however, was more correlated to rainfall intensity (Figure 8b).Characteristic rainfall temporal regimes in the Brazilian semi-arid region, where higher intensities were verified at the beginning of the event, were also observed.Figure 9 describes the basic statistics for runoff and soil loss by comparing the different surface covers.By far, the bare soil plots produced more total runoff and sediments.Mulch cover was close to natural vegetation cover, but still had a higher runoff and sediment loss.Palma presented an intermediate hydrologic response (see also Figures 6 and 7).According to the box plot in Figure 9, the natural Caatinga cover was efficient at increasing soil water storage; this result was also reported by Caloiero et al. [26].Although the presence of native forest increased water consumption as a result of transpiration, it was also verified that soil moisture of nearby soil surfaces increased in comparison to areas where vegetation was removed.
Brasil et al. [6] studied the importance of the Caatinga canopy in reducing rainfall kinetic energy and increasing soil water storage.In addition to the aforementioned contributions, natural cover enhances infiltration [27] and reduces evapotranspiration [2].
The study by Kiani-Harchegani et al. [28] might explain the presence of outliers and extremes for runoff and soil loss data, which also occurred in their study when rain intensity and surface slope were varied.Such values were associated with rainfall intensities that were outside of the 99% probability of occurrence for the studied dataset [18].Figure 9 describes the basic statistics for runoff and soil loss by comparing the different surface covers.By far, the bare soil plots produced more total runoff and sediments.Mulch cover was close to natural vegetation cover, but still had a higher runoff and sediment loss.Palma presented an intermediate hydrologic response (see also Figures 6 and 7).Runoff generation was strongly related to antecedent moisture conditions of a specific event, as observed in Figure 8a.Soil loss, however, was more correlated to rainfall intensity (Figure 8b).Characteristic rainfall temporal regimes in the Brazilian semi-arid region, where higher intensities were verified at the beginning of the event, were also observed.Figure 9 describes the basic statistics for runoff and soil loss by comparing the different surface covers.By far, the bare soil plots produced more total runoff and sediments.Mulch cover was close to natural vegetation cover, but still had a higher runoff and sediment loss.Palma presented an intermediate hydrologic response (see also Figures 6 and 7).According to the box plot in Figure 9, the natural Caatinga cover was efficient at increasing soil water storage; this result was also reported by Caloiero et al. [26].Although the presence of native forest increased water consumption as a result of transpiration, it was also verified that soil moisture of nearby soil surfaces increased in comparison to areas where vegetation was removed.
Brasil et al. [6] studied the importance of the Caatinga canopy in reducing rainfall kinetic energy and increasing soil water storage.In addition to the aforementioned contributions, natural cover enhances infiltration [27] and reduces evapotranspiration [2].
The study by Kiani-Harchegani et al. [28] might explain the presence of outliers and extremes for runoff and soil loss data, which also occurred in their study when rain intensity and surface slope were varied.Such values were associated with rainfall intensities that were outside of the 99% probability of occurrence for the studied dataset [18].According to the box plot in Figure 9, the natural Caatinga cover was efficient at increasing soil water storage; this result was also reported by Caloiero et al. [26].Although the presence of native forest increased water consumption as a result of transpiration, it was also verified that soil moisture of nearby soil surfaces increased in comparison to areas where vegetation was removed.
Brasil et al. [6] studied the importance of the Caatinga canopy in reducing rainfall kinetic energy and increasing soil water storage.In addition to the aforementioned contributions, natural cover enhances infiltration [27] and reduces evapotranspiration [2].
The study by Kiani-Harchegani et al. [28] might explain the presence of outliers and extremes for runoff and soil loss data, which also occurred in their study when rain intensity and surface slope were varied.Such values were associated with rainfall intensities that were outside of the 99% probability of occurrence for the studied dataset [18].

Aspects Related to Granulometry
Erosion in the bare soil plot provided the highest sand percentage (Figure 10), and the highest percentages of silt and clay were observed for the natural vegetation cover plot.Palma cover crop presented similar behavior to bare soil, since the soil was largely unprotected from direct drop impact, with sand percentage clearly above the silt and clay percentages for the transported sediments.
Water 2019, 11, x FOR PEER REVIEW 11 of 16

Aspects Related to Granulometry
Erosion in the bare soil plot provided the highest sand percentage (Figure 10), and the highest percentages of silt and clay were observed for the natural vegetation cover plot.Palma cover crop presented similar behavior to bare soil, since the soil was largely unprotected from direct drop impact, with sand percentage clearly above the silt and clay percentages for the transported sediments.In general, runoff carried more sand fraction for the earlier events of the year, as noted by the maximum discharges that were observed.Later, more intermediate and finer materials were observed.Sand components in the eroded sediment resulted from both raindrop impact and the higher surface runoff, which resulted in higher velocity and, hence, in higher transportation capacity, observed mainly for bare soil cover conditions.
An experiment to help understand sediment loss was performed by Silveira et al. [29].In spite of the fact that the study was conducted in urban environment, the authors observed that the nature of sediment accumulation was complex, and it depended on the type of soil cover, such as green cover, compacted soil, construction, and slope activities.In addition, combinations of bare soil and steeper slopes have greater contributions to sediment loss.
Sediments for the natural cover (Caatinga) treatment consisted mainly of fine particles, because interception prevented the direct impact of raindrops reaching the soil surface.Observations by Kiani-Harchegani et al. [28] under laboratory conditions highlighted the relevance in considering all the key variables involved in sediment transport processes for a better understanding of grain size dynamics.
Bashari et al. [30], studying the effects of soil textures on soil losses, verified that erosion in soils is complex, and more research is still needed to fully understand the erosion mechanism.
Sediment evolution, in terms of percent of clay, silt, and sand in our study, can be observed for the sequence of rainfall events in Figure 11.Lines show the behavior of the textural displacement in the United States Department of Agriculture (USDA) textural chart.According to Figure 11, there were common textural displacement patterns for each soil cover (Bare, Natur, Mulch, and Palma).In general, runoff carried more sand fraction for the earlier events of the year, as noted by the maximum discharges that were observed.Later, more intermediate and finer materials were observed.Sand components in the eroded sediment resulted from both raindrop impact and the higher surface runoff, which resulted in higher velocity and, hence, in higher transportation capacity, observed mainly for bare soil cover conditions.
An experiment to help understand sediment loss was performed by Silveira et al. [29].In spite of the fact that the study was conducted in urban environment, the authors observed that the nature of sediment accumulation was complex, and it depended on the type of soil cover, such as green cover, compacted soil, construction, and slope activities.In addition, combinations of bare soil and steeper slopes have greater contributions to sediment loss.
Sediments for the natural cover (Caatinga) treatment consisted mainly of fine particles, because interception prevented the direct impact of raindrops reaching the soil surface.Observations by Kiani-Harchegani et al. [28] under laboratory conditions highlighted the relevance in considering all the key variables involved in sediment transport processes for a better understanding of grain size dynamics.
Bashari et al. [30], studying the effects of soil textures on soil losses, verified that erosion in soils is complex, and more research is still needed to fully understand the erosion mechanism.
Sediment evolution, in terms of percent of clay, silt, and sand in our study, can be observed for the sequence of rainfall events in  The flow takes thick material first, and when the soil loses the loose particles, as well as increasing wetness, the finer particles begin to be carried off.The more unprotected soil, the greater the amount of coarse material in the transported material.Bare soil plots, because of their unprotected soil surfaces and higher observed runoffs, are subjected to more energy and sediment transport capacity, especially for heavier sand particles.Such behavior was observed for all events and sites in our study.
The USDA textural charts (Figure 12) highlight the variation of the granulometric composition of the original soil and for runoff composition.For the uncovered soil, it was possible to observe that initially there was a loss of finer sediment, and later a predominance of sand fraction loss occurred.The flow takes thick material first, and when the soil loses the loose particles, as well as increasing wetness, the finer particles begin to be carried off.The more unprotected soil, the greater the amount of coarse material in the transported material.Bare soil plots, because of their unprotected soil surfaces and higher observed runoffs, are subjected to more energy and sediment transport capacity, especially for heavier sand particles.Such behavior was observed for all events and sites in our study.
The USDA textural charts (Figure 12) highlight the variation of the granulometric composition of the original soil and for runoff composition.For the uncovered soil, it was possible to observe that initially there was a loss of finer sediment, and later a predominance of sand fraction loss occurred.
The paths (position sequence, illustrating the evolution in time) in the USDA textural classification chart (triangle) associated with the rainfall events for the bare soil plots have different paths associated with the other coverages.Natural cover presented an unmodified pattern, with coarser grain sizes at the beginning of the runoff evolving to a finer size when soils had a higher water content.
For natural cover conditions, the transported sediment texture was always thinner than for the original soil.For Mulch cover, however, rainfall events produced dynamic soil texture losses around the original soil composition.

Figure 1 .
Figure 1.Location of the study sites of Malaquias, Edivaldo, and João in Alto Ipanema River Basin, Pernambuco State, Brazil.

Figure 1 .
Figure 1.Location of the study sites of Malaquias, Edivaldo, and João in Alto Ipanema River Basin, Pernambuco State, Brazil.

Figure 2 .
Figure 2. Photographs of the experimental plots for the different cover conditions (plot with bare soil-Bare; plot with natural cover-Natur; plot with mulch-Mulch; and plot with Cactus Palma-Palma).

Figure 2 .
Figure 2. Photographs of the experimental plots for the different cover conditions (plot with bare soil-Bare; plot with natural cover-Natur; plot with mulch-Mulch; and plot with Cactus Palma-Palma).

Figure 3 .
Figure 3. Sketch of the experimental layout.

Figure 3 .
Figure 3. Sketch of the experimental layout.

Figure 4 .
Figure 4. Time series of rainfall, runoff, sediment loss, and soil moisture during the study period.The cover condition identification was inserted in the charts of João, thus representing all the other sites.

Figure 4 .
Figure 4. Time series of rainfall, runoff, sediment loss, and soil moisture during the study period.The cover condition identification was inserted in the charts of João, thus representing all the other sites.
by the same letter do not significantly differ in the same column, according to the Tukey test (p < 0.05).

Figure 5 .
Figure 5. Photograph detailing the soil surface after a rainfall event (plot with mulch cover).

Figure 5 .
Figure 5. Photograph detailing the soil surface after a rainfall event (plot with mulch cover).

Water 2019 ,
11, x FOR PEER REVIEW 9 of 16 and (3) water retention of the mulch cover.In addition, contour barriers provided by the Palma reduced overland flow velocities, favoring sediment deposition.

Figure 6 .
Figure 6.Runoff coefficient as a function of rainfall volume and intensity, for all sites and plots.

Figure 7 .
Figure 7. Soil loss as a function of rainfall volume and intensity, for all sites and plots.

Figure 6 .
Figure 6.Runoff coefficient as a function of rainfall volume and intensity, for all sites and plots.

Figure 6 .
Figure 6.Runoff coefficient as a function of rainfall volume and intensity, for all sites and plots.

Figure 7 .
Figure 7. Soil loss as a function of rainfall volume and intensity, for all sites and plots.Figure 7. Soil loss as a function of rainfall volume and intensity, for all sites and plots.

Figure 7 .
Figure 7. Soil loss as a function of rainfall volume and intensity, for all sites and plots.Figure 7. Soil loss as a function of rainfall volume and intensity, for all sites and plots.

Figure 8 .
Figure 8. Runoff coefficient (a) and soil loss (b) for all sites and all rainfall events.Antecedent soil moisture conditions are indicated by triangles, representing high antecedent soil moisture (moist conditions), and circles, representing low antecedent soil moisture (dry conditions).

Figure 8 .
Figure 8. Runoff coefficient (a) and soil loss (b) for all sites and all rainfall events.Antecedent soil moisture conditions are indicated by triangles, representing high antecedent soil moisture (moist conditions), and circles, representing low antecedent soil moisture (dry conditions).

Water 2019 ,
11, x FOR PEER REVIEW 10 of 16

Figure 8 .
Figure 8. Runoff coefficient (a) and soil loss (b) for all sites and all rainfall events.Antecedent soil moisture conditions are indicated by triangles, representing high antecedent soil moisture (moist conditions), and circles, representing low antecedent soil moisture (dry conditions).

Figure 11 .
Lines show the behavior of the textural displacement in the United States Department of Agriculture (USDA) textural chart.According to Figure 11, there were common textural displacement patterns for each soil cover (Bare, Natur, Mulch, and Palma).

Figure 11 .
Figure 11.Tracking percentages of sand, silt, and clay for the three sites in all soil covers, together with the seven rainfall events that produced erosion in the plots.

Figure 11 .
Figure 11.Tracking percentages of sand, silt, and clay for the three sites in all soil covers, together with the seven rainfall events that produced erosion in the plots.

Table 2 .
General characteristics (mean values and Tukey test for all seven events in each site) for the plots with surface covers: n-Natural; b-Bare soil; m-Mulch; and p-Palma.