**Soil Structure Stability under Di**ff**erent Land Uses in Association with Polyacrylamide E**ff**ects**

**Amrakh I. Mamedov 1,\*, Atsushi Tsunekawa 1, Nigussie Haregeweyn 2, Mitsuru Tsubo 1, Haruyuki Fujimaki 1, Takayuki Kawai 1, Birhanu Kebede 3, Temesgen Mulualem 4, Getu Abebe 3, Anteneh Wubet <sup>4</sup> and Guy J. Levy <sup>5</sup>**


#### Academic Editor: Alejandro Rescia

Received: 1 January 2021; Accepted: 26 January 2021; Published: 29 January 2021

**Abstract:** Soil structural stability is a vital aspect of soil quality and functions, and of maintaining sustainable land management. The objective of this study was to compare the contribution of four long-term land-use systems (crop, bush, grass, and forest) coupled with anionic polyacrylamide (PAM = 0, 25, and 200 mg L<sup>−</sup>1) application on the structural stability of soils in three watersheds of Ethiopia varying in elevation. Effect of treatments on soil structural stability indices were assessed using the high energy moisture characteristic (HEMC, 0–50 hPa) method, which provides (i) water retention model parameters α and *n*, and (ii) soil structure index (SI). Soil (watershed), land use and PAM treatments had significant effects on the shape of the water retention curves (α, *n*) and SI, with diverse changes in the macropore sizes (60–250; >250 μm). Soil organic carbon (SOC) content and SI were strongly related to soil pH, CaCO3 soil type-clay mineralogy, exchangeable Ca2+, and Na<sup>+</sup> (negatively). The order of soil SI (0.013–0.064 hPa−1) and SOC (1.4–8.1%) by land use was similar (forest > grass > bush > cropland). PAM effect on increasing soil SI (1.2–2.0 times), was inversely related to SOC content, being also pronounced in soils from watersheds of low (Vertisol) and medium (Luvisol) elevation, and the cropland soil from high (Acrisol) elevation. Treating cropland soils with a high PAM rate yielded greater SI (0.028–0.042 hPa−1) than untreated bush- and grassland soils (0.021–0.033 hPa−1). For sustainable management and faster improvement in soil physical quality, soil properties, and land-use history should be considered together with PAM application.

**Keywords:** land use; soil organic carbon; structure stability; soil type; polyacrylamide; dryland

#### **1. Introduction**

#### *1.1. Soil Structure Stability and Its Importance*

Soil structure and aggregation are central physical properties of soil that control a wide array of soil properties and functions including water retention and infiltration [1], susceptibility to erosion and the movement of associated contaminants [2], aeration, gaseous exchanges, and greenhouse gas emission [3], C sequestration, soil organic carbon (SOC) protection [4], soil organic matter mineralization [5], and biogeochemical cycling of essential elements such as macro- and micronutrients [6]. Hence, monitoring of soil structure and stability is vital in determining the sustainability of land use and management practices in both agricultural and natural ecosystems [7].

In Ethiopian highlands, like to similar regions, soils under long-term grassland and forest may contain higher levels of SOC content, which is characterized by persistent particulate and mineral-associated fractions, enhanced microbial activity, and a positive C balance [8]. All those, in turn, may improve soil hydraulic properties, and yield a better soil structure with a heterogenic aggregate and pore size distribution (PSD), compared with soils under long-term cropping, where soil structure is substantially modified by soil management [9,10]. Changing land use from grassland or forest to cropland significantly impacts soil quality, decreases SOC, and negatively affects soil functions such as soil aggregation and gas exchange, water and nutrients availability, and plant growth [7,11,12]. Restoring good soil structure and the proper functioning of degraded cropland soils, by using traditional conservation measures such as no-till practices, are highly dependent on climate and soil type, and texture, and may take over 20–40 years; yet significant positive trends could be noted already after 5–10 years of conversion of cropped land to grassland [12–16].

#### *1.2. Ethiopian Highland Soils: Land Use and Soil Degradation*

The Ethiopian highlands, an important region of natural biodiversity and agriculture production, are characterized by high erosion rates (interrill, rill, and gully, >40 t ha<sup>−</sup>1) and of CO2 efflux, and the deterioration of soil hydrological characteristics and cycles. The issue of land degradation is further complicated by global climate change associated with monsoon–type behavior, delivering significant June–September rainfall and flooding, that turn the soil resources vulnerable to physical and chemical degradation [17,18]. Sustainable watershed management and regreening projects (e.g., conservation systems, building terraces, the establishment of agroforestry), promoted by the Ethiopian government about 20 years ago, have led to satisfactorily integrated sustainable land management and livelihoods activities. Yet, erosion in many areas remains challenging. Moreover, the predicted changes in climate and land use indicate great variability and high levels of erosion rates in the upstream areas and increasing sediment loads in the Blue Nile [19].

Numerous studies in natural and farmed field plots in Ethiopian watersheds, on the impacts of short- and long-term land use and management (e.g., agroforest, grass, and croplands), and conservation measures (e.g., soil bunds by grass, use of enclosure, restrained grazing) on the improvement of soil quality, have recently been performed. The results for SOC, cation exchange capacity (CEC), exchangeable cations (Ca2<sup>+</sup>, Mg2+, Na+), pH, and available nutrients (P and K), soil microbial activity, aggregate stability, crop productivity, and erosion, showed that the efficiency of change in land use and implementation of conservation practices were related to the prevailing agro-ecology and the nature and duration of the used conservation measures [20,21]. Results of these local studies were in line with results from international ones [11,14–16,22]. Moreover, the decline in soil quality by land-use change and under high rain intensity, were related to CEC and clay depletion, translocation and leaching or losses of particulate and dissolved organic matter by runoff and erosion, and clogging of subsoil macropores by the suspended clay-size material. Heavy rains may also cause soil-saturation-related runoff and subsequent shallow lateral flow that increases the loss of fine particles from diverse land uses [8,23,24].

#### *1.3. Polyacrylamide (PAM) as a Soil Stabilizing Agent*

The future of sustainable land management should be based on appropriate practices under different land-use types, which allow the existence of agricultural production in balance with crop and soil systems to support adequate drainage capacity and soil quality [12]. It should be noted that, despite its importance, the contribution of land use, conservation measures, and amendments (e.g., polyacrylamide, lime, manure, biochar) on soil hydraulic properties (water retention, saturated-unsaturated hydraulic conductivity, effective porosity, infiltration rate) and structure stability of Ethiopian soils, has received merely limited attention [25–27].

Application of polyacrylamide (PAM) can considerably (i) enhance soil structure, pore continuity, and aeration via binding of soil particles and eliminating loss of clay particles [28], (ii) sustain SOC accumulation through increasing microbial activity, protecting of SOC in macro- and micro-aggregates, and slowing down of residue decomposition and mineralization rate [29,30], (iii) increase SOC accumulation in subsurface soil layers by providing an appropriate medium for the growth of plants (root penetration, uptake of water and nutrients) [31,32], and (iv) mitigate on-site and off-site impacts of erosion by decreasing runoff generation and erosion rate. PAM addition can, therefore, control the quality of water in various land use areas and could be used for restoring marginal, grass, and forest land [33,34].

Recent studies also revealed that application of PAM could improve soil structure stability and protect SOC, and increase the diversity of soil bacterial communities by increasing soil moisture content and regulating the C to N ratio, regardless of soil with and without plants [35,36]. Therefore, it is expected that PAM application to the Ethiopian highland soils could improve soil structure stability that in turn may facilitate SOC protection [32] and possibly increase SOC storage. It is further anticipated that the use of PAM could replace, at least partly, traditional conservation measures, which often include high-cost conservation practices whose impact may emerge only after 3–4 decades [13,16]. The objective of this study was, therefore, to compare the contribution of four long-term land-use systems (crop, bush, grass, and forest) resulting in significantly varying soil properties and SOC content, with and without anionic PAM application, on structure stability of soils in three watersheds of Ethiopian highlands varying in elevation level.

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

#### *2.1. Study Area and Soil Sampling*

The Upper Blue Nile basin of Ethiopia is characterized by a monsoon climate with a dry (November to April) and a wet (May to October) season, with >80% of the annual precipitation occurring from June to September. Three watersheds in this basin, Guder, Abagerima, and Dibatie, differing in soils, elevation (2500–2900, 1900–2100, and 1500–1700 m), mean annual rainfall (2450, 1340, and 1020 mm), and temperature (19, 23, and 24 ◦C), respectively, represent three different agro-ecological zones, that should be considered for sustainable land management [37]. At these three watersheds, livestock types are similar, and croplands occupy a larger percentage of the area than grazing, bushlands, and forest [37]. In Guder, soil type was Acrisol (Alfisol/Ultisol by US taxonomy) with an acidic reaction, and mixed clay mineralogy with a high content of kaolinite (and Fe, Al oxides), and a small fraction of smectite. In Abagerima, the soil type was Luvisol (Cambisol/Alfisol by US taxonomy) with an acidic reaction and mixed kaolinite-smectite clay mineralogy. In Dibatie, the soil type was Vertisol with a slightly acidic reaction and predominantly smectite clay mineralogy. The parent material of the studied area is related to volcanic rocks (e.g., mainly basalt, then rhyolite, tuff, ignimbrites). Both Acrisol and Luvisol develop agric horizons resulting from clay dispersion, transport, and accumulation [23,38].

Soil samples (0–15 cm) were collected from experimental field plots (~0.1–0.2 ha) located in each of the three watersheds. The plots were placed in four long-term land-use systems (cropland, degraded bushland, grassland for grazing, and forest, hereafter referred to as crop, grass, bush, and forest land) that exhibited varying soil properties, including significantly differing SOC contents (Table 1). In each of the 12 sites (3 watersheds × 4 land uses), three randomized samples (used as replicates), each of ~0.5 kg were taken from the same experimental plot, using a soil probe sampler. The samples were brought to the laboratory, air dried, and ground to pass through a 2-mm sieve and then analyzed. The soils were characterized for (i) electrical conductivity (EC) and pH in a 1:2.5 soil:water extraction; (ii) particle size distribution using the hydrometer method, (iii) cation exchange capacity by sodium acetate, (iv) exchangeable cations by ammonium acetate, (v) calcium carbonate content using the

volumetric calcimeter method, and (vi) organic matter content by wet combustion [39]. The determined properties of the soils are presented in Table 1. Also, aggregate separation to group sizes of 0.5–1.0 mm by dry sieving was done for HEMC measurement, and <0.25 mm by wet sieving was performed for SOC determination in micro-aggregates to be compared with SOC of bulk soil (<2 mm) [39].

**Table 1.** Selected physical-chemical properties (mean ± standard deviation) of the soils sampled from 0–15 cm depth in the study sites. In each column, means labeled with the same letter are not significantly different at *p* < 0.05 based on the Tukey–Kramer HSD test.


pH and EC were determined in the solution part of 1:2.5 soil:water extract using deionized water. EC: electrical conductivity, CEC: cation exchange capacity, SOC: soil organic matter.

#### *2.2. Preparation of PAM-Treated Soil Aggregates*

Treating aggregates with PAM solution (0, 25, and 200 mg L<sup>−</sup>1) was conducted in accordance with the procedure detailed in Mamedov et al. (2010) [40]. The tested PAM concentration were based on previous relevant studies, where various PAM rates were used to understand the mechanism of soil structure stabilization, and to evaluate a strategy of PAM application along with other conservation measures [32–34]. An anionic PAM of high molecular weight (~18 <sup>×</sup> 106 Da) and 30% hydrolysis with a trading name of Superfloc A–110 (Kemira, Lakeland, Florida, USA) was used. Solutions of 25 and 200 mg L−<sup>1</sup> PAM were prepared with tap water (electrical conductivity, [EC] = 0.4 dS m−1, sodium adsorption ratio, [SAR] of = 1.2 [mmolc/L]0.5, and pH = 6.5) under constant stirring and slow addition of PAM granules over 4 h. Plastic boxes (30 × 60 × 3 cm) were filled with very coarse sand to form a 5-mm thick layer that was then covered with a high porosity (>100 μm pore size) filter paper allowing PAM molecules to diffuse to the aggregates from the coarse sand layer. Aggregates (0.5–1.0 mm) from a given soil sample were gently spread on the filter paper to form a monolayer of aggregates, and then saturated from below with 25 or 200 mg L−<sup>1</sup> PAM solution for 1 h (at a rate of 4 mm h<sup>−</sup>1) using a peristaltic pump, and the boxes were then covered and kept in their respective solution for 24 h to reach equilibrium. Thereafter, the boxes were drained and the aggregates were placed in an oven to dry at 60 ◦C for 24 h, and then the aggregates were sieved to eliminate broken aggregates.

#### *2.3. Determination of Soil Structural Stability Indices*

Soil pore size distribution (PSD) and structure stability indices were determined for 108 samples (3 watersheds × 4 land uses × 3 PAM rates × 3 replicates), using the high–energy moisture characteristics (HEMC) method. In the HEMC method, aggregates' wetting rate is accurately controlled, and energy of hydration, differential swelling, and compression of entrapped air are the main forces responsible for breaking down of aggregates. This method enables the detection of small changes in soil structure and had been successfully applied in the determination of structure stability indices of soils from humid and arid regions with a wide range of stability [41,42].

Briefly, 15 g of aggregates (0.5–1.0 mm) were placed in a 60 mm I.D. funnel with a fritted disc (pore size: 20–40 μm) to form a bed ~5 mm thick with a bulk density of ~1.05 g cm<sup>−</sup>3. Saturation of the fritted disc was ensured prior to placing aggregates in the funnel. The aggregates were wetted (fast = 100 mm h<sup>−</sup>1) from the bottom with the deionized water (EC = 0.04 dS m<sup>−</sup>1) using a peristaltic pump. Then a soil water retention curve (0 to −50 hPa), corresponding to drainable pores of 60 to 2000 μm, using small steps (1–2 hPa), was performed by the hanging column and pipette. Consequently, structure stability indices were then inferred from differences among the water retention and specific water capacity (dθ/dψ) curves (i.e., differences in PSD) of the treatments (Figure 1) by using the modified van Genuchten model [43,44]:

$$\boldsymbol{\Theta} = \boldsymbol{\Theta}\_{\rm r} + (\boldsymbol{\Theta}\_{\rm s} - \boldsymbol{\Theta}\_{\rm r}) \left[ \mathbf{1} + (\boldsymbol{\alpha} \,\boldsymbol{\upmu})^{\rm n} \right]^{(1/n - 1)} + \mathbf{A} \boldsymbol{\upmu}^2 + \mathbf{B} \boldsymbol{\upmu} + \mathbf{C} \tag{1}$$

$$\mathbf{d}\boldsymbol{\theta}/\mathbf{d}\boldsymbol{\psi} = (\boldsymbol{\theta}\_{\text{s}} - \boldsymbol{\theta}\_{\text{r}}) \Big[ \mathbf{1} + (\boldsymbol{\alpha}\,\boldsymbol{\uppsi})^{\mathbf{n}} \Big]^{(1/\mathbf{n}-1)} (1/\mathbf{n}-1)(\boldsymbol{\alpha}\,\boldsymbol{\uppsi})^{\mathbf{n}} \mathbf{n} \Big/ \Big[ \boldsymbol{\uppsi} (1 + (\boldsymbol{\alpha}\,\boldsymbol{\uppsi})^{\mathbf{n}}) \Big] + 2\mathbf{A}\boldsymbol{\uppsi} + \mathbf{B} \tag{2}$$

$$\text{SI} = \text{VDP/MS} \tag{3}$$

In Equations (1)–(3), ψ is the matric potential (hPa), θs and θr are pseudo saturated and residual gravimetric water contents, respectively. The parameters α (hPa<sup>−</sup>1), and *n* (dimensionless) represent the location of the inflection point and the steepness of the water retention curve, and the reciprocal of α is often equated with the air entry suction; A, B, and C are quadratic terms used to improve the fitting of the model to the water retention curve [43]. Equation (2) is used to determine the volume of drainable pores (VDP, kg kg<sup>−</sup>1), defined as the integral of the area under the specific water capacity curve (dθ/dψ) and above its baseline; MS (hPa) is the modal suction corresponding to the matric potential at the peak of the specific water capacity curve (Figure 1d) and relates to the most frequent pore size. The SI in Equation (3) is used to evaluate the susceptibility of the tested samples to breakdown; the higher the value of SI the less susceptible the aggregates are to slaking. The coefficient of variation between replicates of water content (θ, kg kg<sup>−</sup>1) was <5%.

#### *2.4. Statistical Analysis*

Analysis of variance (ANOVA) tests were conducted using the SAS Proc GLM procedure [45] to assess the effects of soil type (watershed elevation), land use and PAM treatments and their interactions on soil properties, near saturation water retention curve parameters (α, *n*), and HEMC stability indices (θs, θr, SI). Treatments mean comparisons were done using the Tukey-Kramer HSD test at *p* < 0.05 (Supplementary material: Table S1). Pearson pairwise correlation, regression analysis (*p* < 0.05), and stepwise regression analysis were used to examine the effects of soil properties (e.g., SOC, EC, pH, CEC, cations) on structure stability indices (Table S2).

**Figure 1.** Water retention curves of the soils form (**a**) Guder (Acrisol), (**b**) Abagerima (Luvisol), and (**c**) Dibatie (Vertisol) watersheds used under crop, grass, bush, and forest land, and (**d**) the specific water capacity curves of soils from Abagerima (Luvisol). The dashed baseline in the specific water capacity curve represents the soil shrinkage line.

#### **3. Results**

Results of the ANOVA tests revealed that all analyzed soil properties (Table 1), structural stability indices, and parameters (α, *n*) were significantly (*p* < 0.001) affected by soil type (watershed elevation), land use, and their interaction (Table S1). The VDP and SI, *n* and MS were highly correlated [41], and thus we opted to focus the discussion mostly on the SI and the model parameter α in the presented analyses.

#### *3.1. Land Use and Soil Properties*

Studying soil properties in different land uses is crucial for understanding the current prevailing conditions in soils from different agro-ecological zones from conservation measures strategy perspective [37,45]. Soil properties showed large variation in pH (5.2–6.7), SOC (1.4–8.1%), cation exchange capacity (CEC, 25–40 cmolckg<sup>−</sup>1), exchangeable Ca2<sup>+</sup> (9.4–16.1 cmolckg<sup>−</sup>1), soil texture (sandy loam to clay), electrical conductivity (EC, 0.2–0.9 dS m<sup>−</sup>1), and CaCO3 (0.3–5.8%) (Table 1). Some clear associations of certain soil properties (e.g., SOC, pH, EC), with elevation and land use were noted. Content of SOC was highest and pH was lowest in Guder watershed (highest elevation and rainfall), followed by Abagerima (medium elevation and rainfall) and Dibatie (lowest elevation and rainfall). Differences in soil clay content (texture) were also linked to watersheds: Guder (loam) < Abagerima (clay loam) < Dibatie (clay). The Clay/CEC ratio (i.e., indication of the clay mineralogy [40], also decreased with the increase in elevation: Guder (0.3–0.7) < Abagerima (0.7–0.9) < Dibatie (1.6–2.5)

(Table 1). In each watershed, soil pH and SOC were related to land use, and decreased in the following order: Forest > Grass <sup>≥</sup> Bush > Crop. Largely, soil CEC, exchangeable Ca2<sup>+</sup> and K+, and EC and CaCO3 were lowest under cropland and highest under forest, grass, or bushland. The changes in exchangeable Mg2<sup>+</sup> (3.2–4.5 cmolckg<sup>−</sup>1) and Na<sup>+</sup> (2.4–6.3 cmolckg<sup>−</sup>1) with variations in land use were not consistent, yet the latter mostly was higher under cropland or bushland (Table 1).

In general, only in 10 out of 85 cases, meaningful correlations (i.e., *r* ≥ 0.7) were noted among the various soil properties studied (Table 2a). However, in each watershed, strong correlations (50% of cases) between SI or SOC and the mean of soil properties were found (Table 2b). For instance, soil pH or EC were weakly related to SOC and SI (*r* < 0.3–05) when all samples were involved (Table 2a). Nevertheless, when separating correlation analysis by watersheds, pH was very strongly associated (*r* > 0.9) with SOC and SI in the acidic Acrisol and Luvisol, and moderately with slightly acidic Vertisol; soil EC was closely related to SOC and SI in Luvisol and Vertisol (*r* > 0.7–0.9), but not in Acrisol (*r* < 0.7) (Table 2b). Correlation analysis for all soils studied indicated that exchangeable Ca2<sup>+</sup> was closely related to SOC (*r* = 0.7), while it was not related to CaCO3 content (*r* = 0.3) (Table 2a). However, analysis for each watershed alone showed better relation between Ca2<sup>+</sup> and CaCO3 (*r* = 0.5–0.6) (Table S2).

**Table 2.** Pearson pair-wise correlation coefficients for properties and structural index (SI) of the soils used. The correlations with *r* ≥ 0.7 are in bold. Units of properties are in Table 1.


EC: electrical conductivity, CEC: cation exchange capacity, SOC: soil organic carbon, CCR: clay to CEC ratio, ESP: exchangeable sodium percentage; \*, \*\*: significant at the 0.01 (0.05) and 0.001 level accordingly.

#### *3.2. Water Retention and Structure Stability of Untreated and PAM-Treated Soils*

#### 3.2.1. Water Retention of Untreated and PAM-Treated Soils

An evaluation of the shape of the water retention curves indicates the existence of differences associated with the distribution of the studied pores (>60 μm), with pore diameter being calculated from the matric potential (d = −300/ψ, where d is the equivalent pore diameter, μm). Land use and PAM treatments had considerable effects on the shape of the near saturation water retention curves, and consequently on its model parameters (θs, θr, α, *n*), and the soil structural stability indices (VDP, MS, SI) (Figures 1–3, Table 3).

**Figure 2.** Water retention curves of untreated (PAM 0) and polyacrylamide (PAM) treated (25 and 200 mg L<sup>−</sup>1) soils used under crop, grass, and forest land for (**a**–**c**) Guder (Acrisol), (**d**–**f**) Abagerima (Luvisol), and (**g**–i) Dibatie (Vertisol) watersheds.

**Figure 3.** Structural index (SI) of the soils as affected by (**a**) soil type (watershed), (**b**) land use, and (**c**) polyacrylamide (PAM) treatments (0—untreated, 25 and 200 mg L<sup>−</sup>1). Columns labeled with the same letter are not significantly different at *p* < 0.05 level: (**a**), within the soil (watershed); (**b**), within land use; and (**c**), within each soil (watershed).


**Table 3.** Effects of the treatments on the HEMC structure stability indices (VDP, MS, and SI) and near saturation water retention model parameters (θs, θr, α, *n*). Within each factor, the columns labeled with the same letter are not significantly different at *p* < 0.05 level.

VDP: volume of drainable pores, MS: modal suction, SI: structural index, θs and θr: saturated and residual water contents, α and *n*: represent the location and steepness of the water retention curve.

In all watersheds, the water retention curves of untreated soils for the different land uses exhibited similar trends, but of different magnitudes (Figure 1, Table 3). In the relative wet (−ψ > 12 hPa; pore size > 250 μm) and dry (−ψ < 24 hPa; pore size 60–125 μm) parts of the HEMC curves, soil moisture content (θ) at a given matric potential (ψ) was positioned in the following order of land use: forest > grass > bush > crop (excluding bushland soil at the dry end which was similar to the crop one). In the medium matric potential range (−ψ =12–24 hPa; pore size 125–250 μm) such differences for the different land uses were inconsistent (Figure 1, Table 3).

Treating the soils with PAM modified the shape of water retention curves in the entire range of the matric potential studied (ψ = 0–50 hPa) compared with the untreated ones, with the effect of PAM being dependent on PAM rate, land use, and soil type (watershed) (Table 3). In all three watersheds, water retention curves of PAM-treated soils were mostly to the left side of (or below) those of the untreated soils in crop and grassland soils (results for bushland were similar to grassland—data not presented). Moreover, a trend was noted whereby, the difference between the water retention of untreated and PAM-treated soil was higher in the cropland soils compared with the grass or forest ones. In the forest soil, the water retention curves of the untreated and PAM-treated samples were comparable in Guder and Abagerima watersheds, whereas in the Dibatie watershed (with very sparse forest) water retention of samples was somewhat similar to those observed in the grassland soils (Figure 2).

The contribution of the PAM rate in Acrisol (Guder), yielded only small differences between the retention curves of the two PAM rates at the entire range of the matric potential (Figure 2a–c). By contrast, in Luvisol (Abagerima), a distinct difference between the retention curves for the two rates of PAM was noted, especially at ψ > 10–12 hPa in the crop and grassland, where water retention curves of high rate PAM-treated soils were on the left side of the low rate PAM-treated soils, revealing to more draining water at a given matric potential (Figure 2d–f). In Vertisol (Dibatie), the impact of PAM rate on the water retention curves was similar to that noted in the Luvisol (Abagerima), but it extended also to the forest soil (Figure 2g–i).

#### 3.2.2. Structure Stability of Untreated and PAM-Treated Soils

Structure stability indices and model parameters by soil type (watershed), land use, or PAM differed: the *n* (and MS) decreased, while θs, θr, α, and SI (and VDP), increased with the following order of land use: crop < bush < grass < forest (Table 3). The SI values of the individual treatments varied over a wide range, between 0.013 to 0.079 hPa−1, which could be related to the difference in soil type and properties among the watersheds and land use. However, within each watershed, and regardless of land use, PAM addition significantly increased soil SI. Moreover, similar to SOC order, untreated and PAM-treated soils had the same order of SI by land use: crop < bush < grass < forest (Figure 3).

In each watershed, the PAM effect on SI was highest in cropland and lowest in the forest (Figure 3). Relative to the untreated soils, PAM application increased soil SI by 1.2–1.3, 1.2–1.8, and 1.6–2.3 times in the soils of Guder (Acrisol), Abagerima (Luvisol), and Dibatie (Vertisol) watersheds respectively. At a high rate, PAM increased SI in the following order: Vertisol ≥ Luvisol > Acrisol for each type of land use (expect forest in Abagerima), while at a low rate of PAM and for untreated samples, the order of SI for same land use between the watersheds (or soil types) was inconsistent, although untreated soils under bush, grass, and forest land from Guder and Abagerima were more stable than from Dibatie watershed (Figure 3). It should be noted that treating cropland soils with a high PAM rate yielded greater or comparable SI (0.028–0.042 hPa−1) than untreated grass and bushland soils (0.021–0.033 hPa<sup>−</sup>1) in all watersheds (Figure 3). However, increase in PAM rate was (i) effective only in cropland soil with the lowest SOC (<2%) in Acrisol, (ii) not significant only in forest soil with highest SOC (~8.0%) in Luvisol, and (iii) significant for all land use samples in Vertisol (SOC = 1.7–5.1%). In the Guder watershed, the PAM rate has not changed the SI of bush, grass, and forest land samples (Figure 3).

A strong linear relation (R<sup>2</sup> = 0.72) between soil SI and SOC for the untreated soils, and for the low rate PAM-treated soil (R2 = 0.64), was noted; for the high rate PAM-treated soils, no meaningful relation (R2 = 0.26) between SI and SOC was obtained (Figure 4a). The relation between SOC of micro-aggregates (<0.25 mm) and that of bulk soil exhibited a strong linear relation that could be associated with the effect of micro-aggregate stabilization by SOC on soil SI (Figure 4a,b).

**Figure 4.** Relations (*p* < 0.001) between (**a**) soil structural index (SI) and soil organic carbon (SOC) for untreated (PAM 0) and polyacrylamide (PAM)-treated soils (25, 200 mg L<sup>−</sup>1), and (**b**) SOC of micro-aggregates (<0.25 mm) and bulk soil (<2 mm).

A stepwise regression analysis revealed that in untreated soils, 92% of the variation in SI was associated mainly with SOC (71%), and CaCO3 (+12%). For the low rate PAM-treated soils, variation in SI was also associated with SOC (62%) and EC (+17%). However, for the high PAM application rate, 96% of the variation was related to EC, silt, SOC, pH, and clay content (Table 4). There was an exponential relation (3 watersheds × 4 land use × 3 PAM rate = 36 treatments) between soil SI and water retention model parameter α; the strength of this relation was inversely related to SOC level (Figure 5).


**Table 4.** Stepwise regression analysis of the effect of soil properties on the structural index (SI), for the untreated (PAM 0) and PAM-treated (25, 200 mg L<sup>−</sup>1) samples.

EC: electrical conductivity, CEC: cation exchange capacity; SOC: soil organic carbon, CCR: clay to CEC ratio.

**Figure 5.** Relations (*p* < 0.001) between soil water retention model parameter α and soil structural index (SI) for all samples (SOC ≤ 8.1%, solid line) and samples with SOC ≤ 3% (dash line). Soil samples include untreated (PAM 0) and polyacrylamide (PAM)-treated ones (25, 200 mg L<sup>−</sup>1).

#### **4. Discussion**

#### *4.1. Land Use and Soil Type (Elevation) E*ff*ects*

The impact of land use on the soil properties (Table 1) was in agreement with other studies on similar soils or from those in Ethiopia, revealing that topsoil attributes, such as CEC and exchangeable cations (Ca2<sup>+</sup>, Mg2<sup>+</sup>, or Na+), SOC, pH, and aggregate stability were considerably higher, in soils under forest, grass or pasture lands than under cropland [8,20,21]. Results of the pair-wise correlation analysis for all the soils studied (Table 2a) highlight the considerable contribution of land use type coupled with soil type/elevation on the examined soil properties. These observations were in agreement with recent studies [37,46]. Moreover, the variation in the soil properties and structure stability indices under various land use (Tables 1–3) could be related to the sensitivity of characteristics of the studied soil types to structure-modifying processes [12,14,16], which is also associated with elevation or watersheds. Soil structure stabilization in the studied watersheds of the Ethiopian highland under long-term conservation practices leads to a reduction of runoff and soil loss, minimizing loss of crop residues and

organic components, and to SOC accumulation [37,46,47]. However, the role of soil CaCO3 on soil biogeochemistry (e.g., pH, EC, Ca2<sup>+</sup>, CEC) in general and in particular its contribution to improving soil structure stability as observed in our study (Table 4), is still greatly overlooked in acidic soils, including in soils from Ethiopian highland [48].

Aggregate-structure stability (expressed in terms of SI or α) was directly associated with SOC irrespective of elevation (Table 2). The distinct relation between SOC and aggregate stability has been recognized long ago and has been ever since verified in many studies [21,44,47]. Aggregate stability was directly related to exchangeable Ca2<sup>+</sup> in Luvisol and Vertisol (Abagerima and Dibatie) and CaCO3 content in Acrisol and Luvisol (Guder and Abagerima) (Table 2b). Evidentially, to maintain stable structure in soils of Ethiopian highlands, agricultural practices should include, in addition to means for increasing soils SOC, timely monitoring of soil Ca2<sup>+</sup> and, if needed, adding some external source of Ca2<sup>+</sup> (e.g., lime), especially at high altitudes, to balance Ca2<sup>+</sup> deficiency and storage [27,34,48]. Similar to soils from semi-arid and arid regions, the stability of the soils in the Ethiopian highlands was susceptible to exchangeable Na<sup>+</sup> or sodicity (i.e., ESP). Even though the negative role of sodicity on soil structure stability was evident in our study (Table 2), this topic, as well as the inverse relation between sodicity and SOC, have not been receiving their due attention in studies related to these acidic soils [25,26], although elevated exchangeable Na<sup>+</sup> was also recently noted by other studies performed under natural condition or conservation practices in the highland region [8,46].

#### *4.2. Soil Organic Carbon and Polyacrylamide E*ff*ect on Soil Structure Stability Indices*

The shape of the water retention curves and hence the range of SI in untreated soils is associated with aggregation and can be explained by the conceptual aggregate hierarchy model [49], postulating that aggregates are sequentially formed through the action of organic binding agents leading to the formation of micro-aggregates (20–250 μm) and then macro-aggregates (>250 μm). Micro-aggregates could be formed within macro-aggregates by roots and microbial activity; aggregates and pores provide physical protection of SOC and mineralization [4,5,50]. Aggregate size influences the emissions of CO2, and correlations exist between aggregate size and macroporosity, number of pores, and pore size [3,28]. It should be noted that in the untreated soils, SOC content (and thus SI) is controlled by soil type or clay mineralogy (adsorption of SOC on clay minerals) and it is linked to CEC, pH, and CaCO3, i.e., soil attributes affected by the consequences of land use [4,7,9,51].

The water retention curves for forest, grass or bushland soils differed from that of cropland in the range of macro-pores (>250 μm; −ψ > 12 hPa) and micropores (60–125 μm; 12 < −ψ < 24 hPa) (Figure 2). The deviation of the grass and forest soil water retention curves from that of the cropland one at the micropore size range (60–125 μm) was associated with micro-aggregate stabilization through avoiding soil disturbance by tillage, and higher SOC content [4,6]. The observed strong correlations between SI and SOC (*r* > 0.9) for each soil type, linear relations between SOC and SI (R<sup>2</sup> = 0.72), and between SOC of bulk soil and micro-aggregates (R<sup>2</sup> = 0.92) (Table 2, Figure 4) emphasize that good soil structure could support soil SOC accumulation at micro-and macro-aggregate level [7]. However, stable micro-aggregates play a critical role in the long-term stabilization of soil organic matter, whereas less stable macro-aggregates provide only a small physical protection [50]. In turn, good and stable soil structure facilitate the supply of food source for soil microbes and subsequently SOC development, and hence, positive feedback exists between microbial activity, SOC accumulation, and soil structure formation [36,42,51,52].

Treating soils with anionic PAM was effective in stabilizing existing aggregates and improving bonding between and aggregation of soil particles. The magnitude of the PAM effects, similar to the effects of SOC, were soil type-clay mineralogy and land use (e.g., SOC, EC, pH, Ca2+) dependent (Table 3, Figure 4a) and could be related to the integrated PAM–Soil (Ca2<sup>+</sup>, Mg2<sup>+</sup>)–SOC interaction [51,53]. Positively charged cations (adsorbed on negatively charged mineral surfaces) facilitate the adsorption of negatively charged long-chain organic molecules through cation bridging. Aggregation depends on the strength of bonding governed by the chain length of the organic molecules and the type of

cations on the exchange sites. Aggregates can be held together more strongly by PAM rather than by electrolytes because a single PAM chain is linked to multiple soil particles [54,55]. Thus, the higher SI from PAM treatments could be attributable to the adsorption of the long-chain PAM molecules to soil particles, which act as cementing material, and minimize aggregate slaking, and enhances soil aggregate and structure stability more than organic or inorganic amendments [40,56].

The deviation of the water retention curves of the PAM-treated cropland soils from the untreated one was less evident than the aforementioned deviations observed in the case of the changes in land use (Figures 1 and 2); yet, at times changes took place in the regions of the micro- and meso-size pores (60–250 μm; −ψ < 12 hPa). This phenomenon could possibly be linked to the fact that, in contrast to the non-tilled land uses (forest, grass, and bushland), application of PAM does not lead to the buildup of new soil aggregates with time, but to preserving and strengthening existing aggregates [29,32,41]. Moreover, in the aggregate size range studied (0.5–1.0 mm), adsorption of PAM of high molecular weight used on soils takes place mostly on exterior surfaces of soil material, because the narrow pores in small–size aggregates may not allow penetration of the large PAM molecules into the aggregates [57]. This pattern of PAM adsorption to the aggregates explains the significant impact of the PAM rate on the SI (Table 3); the larger the rate, the greater the SI [28,34,58].

The favorable impact of PAM, irrespective of its rate, on SI of the slightly acidic soils under cropland (Figure 3), was significantly greater than on the other land uses. This observation was also noted in soils from a semi-arid region [58,59]. Moreover, it is worthwhile noting that the contribution of PAM to improving SI (Acrisol < Luvisol < Vertisol) by watershed was opposite to the trend of SOC accumulation in these soils (Acrisol > Luvisol > Vertisol); yet all soil samples with SOC < 3–5% used under crop, grass, and bushland were significantly affected by PAM application (Figure 3). This fact could be seen as analogous to findings from former studies, where the positive effect of PAM on SI was related to (i) soil texture, with PAM being more effective in loam with initially much weaker soil structure than in clay soils [28,29], and (ii) soil clay mineralogy, with PAM being more effective in unstable smectitic soils than in the stable kaolinitic ones [40,59–61]. Our observations regarding PAM contribution to improving aggregate stability being more pronounced in soils with lower SOC than with higher SOC content, further highlight that PAM is more effective in soils of a priori lesser stability. Moreover, at a low PAM rate, SOC and PAM could jointly contribute (Table 4, Figure 4a) to soil quality and macropore stabilization for short or long-term PAM application strategy [32,40].

The exponential relation between soil SI and model parameter α (Figure 5), reveals that an increase in α by SOC or PAM treatment may significantly increase soil SI due to the increase in the mean aggregate and pore size, and resistance of aggregates to slaking [10,34]. Excluding samples with high SOC content (>5%) from the relation increased the coefficient of determination R<sup>2</sup> from 0.68 to 0.87 (Figure 5), yet regression curves with the close exponents were "parallel". Thus, exponential relations between SI and α (which include untreated and PAM-treated samples), could be considered as important and might suggest that along with SI, parameter α could also be used for evaluating changes in aggregate-structure stability following changes in soil intrinsic or extrinsic conditions associated with land use and management, including PAM application [59,60]. Finally, it transpires from our findings that PAM, in combination with other conservation methods, can integrally regulate soil aggregate-structure stability. Thus, PAM can be applied for restoration of crop, grass, bushland soil, or even forest soils, particularly in "hot spot" areas associated with enhanced translocation and leaching of the clay particles and cations in degraded lands from tropical, humid, and arid regions [8,23,32–34,59,62].

#### **5. Conclusions**

It emerged from our study that change of land use from cropland soils (annually tilled) to land uses that do not require tillage (bush, grass, and forest) has substantial positive effects on aggregate-structure stability indices in the studied acidic soils from the Ethiopian highlands. The observed positive effect of CaCO3, and the negative effects of sodicity on SOC accumulation and soil structure development in the acidic soils studied, which have previously been mostly overlooked, need greater attention. The positive effects of the non-tilled land uses were ascribed to SOC accumulation and a lesser extent to CaCO3 content; the magnitude of the effects depended on soil type. A viable alternative in the form of using small amounts of PAM (e.g., a soil amendment), was found effective in stabilizing aggregates and improving soil structure in all the studied soils. The contribution of PAM to improving aggregate stability was more pronounced in soils of lower SOC content, thus emphasizing PAM efficacy as a soil stabilizing agent in soils of a priori weak structure stability. The mechanisms by which SOC and PAM improve aggregate stability seems to differ as predominant changes were observed in the pore size range of macropores (>250 μm) for changes in land use and hence SOC effect, and in the pore size range of 60–250 μm in the case of PAM application.

The efficacy of PAM application suggests that improving soil aggregate-structure stability is obtained almost immediately without the need to go through the time-consuming process of changing land use or using long-term no-till practices, to achieve the Sustainable Development Goals of the United Nations Land Degradation Neutrality challenge [63]. The results of PAM addition to the Ethiopian highland soils met our expectations that PAM addition (i) improves soil structure stability, and (ii) could replace, at least partly, traditional conservation measures, which include high-cost cultivation practices whose impact may emerge only after several decades. As the soil types and type of clay minerals and soil texture cannot be changed by management practices, employing practices such as PAM addition and SOC driven aggregation is crucial for short- and long-term soil structure stabilization under various land uses.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2071-1050/13/3/1407/s1, Table S1: Analysis of variance (ANOVA) of the effect of the treatments on soil properties, and near saturation water retention model parameters (θs, θr, α, *n*), and the HEMC structure stability indices (VDP, MS, and SI). W: watershed, LU: land use, PAM: polyacrylamide, VDP: volume of drainable pores, MS: modal suction, SI: structural index, α and *n*, are the location of the inflection point and the steepness of the water retention curve. θr and θs, are the residual and saturated water content, EC: electrical conductivity (1:2.5), CEC: cation exchange capacity, SOC: soil organic carbon; CCR: the ratio of clay content to CEC (indication of clay mineralogy); \*: *p* < 0.001, Table S2: Pearson pairwise correlation coefficients for properties of the soils used. Units of the respective properties are as in Table S1. The coefficients higher than 0.5 are significant at *p* < 0.01–0.001. EC: electrical conductivity (1:2.5); CEC: cation exchange capacity, SOC: soil organic carbon.

**Author Contributions:** Conceptualization, A.I.M., A.T., N.H., and G.J.L.; methodology, A.I.M., G.J.L., T.K., B.K., T.M., G.A., and A.W.; software, A.I.M. and H.F.; validation, A.I.M., A.T., and G.J.L.; formal analysis, A.I.M., M.T., and N.H.; investigation, A.I.M., H.F., B.K., T.M., A.W., and M.T.; resources, A.T., H.F., and T.K.; data curation, M.T. and H.F.; writing-original draft preparation, A.I.M. and G.J.L.; writing—review and editing, A.I.M., G.J.L., and A.T.; visualization, A.I.M.; supervision, A.T., A.I.M., and N.H.; project administration, A.T.; funding acquisition, A.T. and H.F. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Science and Technology Research Partnership for Sustainable Development (SATREPS)—Development of a Next-Generation Sustainable Land Management (SLM) Framework to Combat Desertification project, Grant Number JPMJSA1601, Japan Science and Technology Agency (JST)/Japan International Cooperation Agency (JICA).

#### **Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Acknowledgments:** The authors thank Arid Land Research Center, Tottori University, The Science and Technology Research Partnership for Sustainable Development (SATREPS), and Japan Science and Technology Agency (JST)/Japan International Cooperation Agency (JICA). We also would like to thank the Kemira Company (Lakeland, Florida, USA) for providing anionic polyacrylamide, SUPERFLOC A-110.

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

#### **References**

1. Chandrasekhar, P.; Kreiselmeier, J.; Schwen, A.; Weninger, T.; Julich, S.; Feger, K.H.; Schwärzel, K. Why we should include soil structural dynamics of agricultural soils in hydrological models. *Water* **2018**, *10*, 1862. [CrossRef]


**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2021 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/).

MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel. +41 61 683 77 34 Fax +41 61 302 89 18 www.mdpi.com

*Sustainability* Editorial Office E-mail: sustainability@mdpi.com www.mdpi.com/journal/sustainability

MDPI St. Alban-Anlage 66 4052 Basel Switzerland

Tel: +41 61 683 77 34 Fax: +41 61 302 89 18