Rheological Parameters as Affected by Water Tension in Subtropical Soils

: Rheological parameters have been used to study the interaction between particles and the structural strength of soils subjected to mechanical stresses, in which soil composition and water content most strongly affect soil resistance to deformation. Our objective was to evaluate the effect of water tension on rheological parameters of soils with different mineralogical, physical, and chemical composition. Surface and subsurface horizons of four Oxisols, two Ultisols, one Alfisol, and one Vertisol were physically and chemically characterized; their rheological parameters were obtained from amplitude sweep tests under oscillatory shear on disturbed soil samples that were saturated and subjected to water tension of 1, 3, 6, and 10 kPa. In these samples, the rheological parameters linear viscoelastic deformation limit (γ L ), maximum shear stress (τ max ), and integral z were determined. By simple regression analysis of the rheological parameters as a function of soil water tension, we observed increased mechanical strength with increasing water tension up to at least 6 kPa, primarily due to increased capillary forces in the soil. However, increased elasticity assessed by γ L was not as expressive as the increase in structural rigidity assessed by τ max and integral z. Elastic deformation of the soil (γ L ) increases with the increase in the number of bonds among particles, which depend on the clay, total carbon, expansive clay mineral, and cation contents; however, maximum shear resistance (τ max ) and structural stiffness (integral z) mainly increase with clay, kaolinite, and oxide content by increasing the strength of interparticle bonds. A decrease in mechanical strength occurs for water tension of 10 kPa (the lowest water content evaluated) in sandy horizons or in horizons with a high proportion of resistant microaggregates (pseudosand), when associated with low bulk density, due to fewer points of contact between soil particles and therefore lower capillary force.


INTRODUCTION
Soil physical degradation from external stresses increases with decreasing structural strength of soil (Horn, 2003).Traditional techniques such as direct shear and static compressibility are widely used for determining mechanical parameters by stress-strain relations correlated with structural properties on the scale of aggregates, such as soil density and porosity and aggregate stability.As these properties depend on soil composition, such as particle size distribution, mineralogy, and organic matter content (Horn and Peth, 2011), the traditional techniques are limited in their ability to clearly determine soil behavior caused by mechanisms acting on a particle scale (Markgraf et al., 2006).
Rheometry is a technique that has been used to access soil reaction to external stress derived from intergranular bonds, quantifying soil structural strength with stress-strain relationships (Markgraf andHorn, 2006, 2007;Khaydapova et al., 2015).Thus, rheometry allows soil deformation to be related to particle size and mineralogy, concentration and type of ions, soil organic matter (SOM), and water content (Markgraf and Horn, 2009).
Particle arrangement during soil deformation is affected primarily by water content.At high water content, friction between soil particles is reduced and the stresses applied first overcome cohesion among water molecules and then the stronger cohesion that exists among solid constituents (Ghezzehei and Or, 2001;Markgraf et al., 2006).A decrease in soil water content strengthens the bonds among soil particles by an increase in cohesive forces, friction between particles, and capillary forces (adhesion) due to water surface tension and the number and curvature of menisci (Kemper and Rosenau, 1984;Gallipoli et al., 2003;Mitchell and Soga, 2005;Lourenço et al., 2012).
Soil water tension (suction) is the manifestation of these capillary forces (Santamarina, 2003).As the soil dries, water recedes in capillaries between particles and aggregates, and adhesion forces and surface tension on water menisci pull adjacent particles at high force (Kemper and Rosenau, 1984;Lourenço et al., 2012).However, capillary force depends on soil composition.In clay soils, which have smaller diameter pores, adhesion forces and capillarity are higher (Horn and Peth, 2011;Lourenço et al., 2012), while in very coarse soils, the number of menisci, and thus capillary force, can decline sharply at water tension higher than 6 kPa (Holthusen et al., 2012b), due to larger pore size.
With reduction in soil water content, there is an increase in cohesion due to an increase in concentration of soluble constituents in the liquid phase, such as carbonates and organic molecules, which may precipitate between particles and aggregates as semi-crystalline inorganic compounds or organic amorphous compounds (Kemper and Rosenau, 1984;Kemper et al., 1987;Mitchell and Soga, 2005).Soils display cohesion because of connections of particles through cations bridges, electrostatic attraction among clay particles, and attraction between soil particles and between water and soil particles by van der Waals forces (Kemper and Rosenau, 1984;Libardi, 2010), as well as the cementing effects of organic substances and iron and aluminum oxides (Reichert et al., 2010;Kämpf et al., 2012).Therefore, soil cohesiveness increases mainly due to decreasing particle size and increasing content of expansive clay minerals (Kemper and Rosenau, 1984;Kemper et al., 1987), SOM content, andcontact area between particles (Reichert et al., 2010).
The effect of these factors on soil mechanical resistance can be evaluated by rheological parameters obtained in amplitude sweep tests under conditions of oscillatory shear as the deformation limit (γ L ), which represents elasticity (recoverable strain) (Mezger, 2014); the maximum shear stress (τ max ), which characterizes the maximum resistance (Holthusen et al., 2012a); and the integral z, which represents structural stiffness (Markgraf et al., 2012).Several recent studies have shown that increasing water tension increases the structural strength of soils with different textures and mineralogy, observed by the increase in these parameters (Markgraf andHorn, 2006, 2007;Holthusen et al., 2010Holthusen et al., , 2012c;;Markgraf et al., 2012;Baumgarten et al., 2012).The increase in water tension from 0 to 3 kPa is sufficient to significantly increase soil shear strength as a result of stabilization of menisci, salt precipitation, and increased friction between particles (Holthusen et al., 2010).
All these factors and processes intervene in the links between particles and aggregates (Mitchell and Soga, 2005), modifying soil structural strength.However, for different Brazilian soils, little has been investigated about microstructural resistance based on rheometry.Thus, we aimed to assess the effect of water tension on rheological parameters of soils from southern Brazil with different mineralogical, physical, and chemical composition.

Study area
The study was conducted with samples of eight soils of different classes from southern Brazil, with a humid subtropical climate according to the Köppen classification system.Soils widely varying in particle size and mineralogy were chosen, belonging to the orders Latossolo/Oxisol (four soils), Argissolo/Ultisol (two soils), Planossolo/Alfisol (one soil) and Vertissolo/Vertisol (one soil), as described in table 1. (1) Classification by the Sistema Brasileiro de Classificação de Solos (Brazilian Soil Classification System) -SiBCS (Santos et al., 2013) and by Soil Taxonomy (USDA, 2010).

Sample collection and preparation
Soil samples were collected in profiles described in other studies.The samples were collected in the middle of the A and B horizons, air-dried, ground, and sieved through a 2 mm mesh.With this soil fraction, analyses were performed to characterize chemical and physical properties (Table 1) and execute rheological analysis, as described below.

Rheometry
Rheological analyses were performed by the amplitude sweep test with controlled deformation using a compact modular rheometer (Anton Paar MCR 102, Austria) fitted with a parallel plate measuring device: a fixed roughened bottom plate of 50 mm diameter and a top roughened rotating plate of 25 mm diameter (model PP25/P2).
The soil samples for testing were prepared by moistening the soil to between 10 and 30 % of gravimetric water content and keeping them in a sealed bag for about 16 h.Moist soil was compacted into rings of approximately 3.6 cm diameter and 1.0 cm height so that the bulk density (ρ s ) of the samples (15 for each horizon) equaled the average of ρ s observed in the field (Table 1).The samples were then saturated with destilled water by capillarity.Three samples from each horizon (group 1) remained saturated and the other samples formed groups 2, 3, 4, and 5, each with three samples, which were drained to the water tensions (w t ) of 1, 3, 6, and 10 kPa, respectively, on a sand tension table (Reinert and Reichert, 2006).Samples were only drained to tensions of up to 10 kPa because in rheological tests if the soil water content is small, the sample cannot deform homogeneously, which can introduce errors in the parameters obtained (Mezger, 2014).
The samples were subsequently placed on the lower plate of the rheometer, extracted from the ring, and cut horizontally and vertically to an approximate height of 4.5 mm and diameter of 25 mm.This procedure was performed immediately before the test to adjust the sample diameter to the diameter of the top plate and the gap (space between the measuring plates) and avoid possible effects from the surrounding soil (Holthusen et al., 2010).
The amplitude sweep tests were performed under the following conditions: constant temperature = 20 °C; sample height (gap) = 4 mm; rest period before beginning the test = 30 s; variation of amplitude of deformation = 0.0001 at 100 %; angular frequency = 0.5 Hz; and number of measuring points = 30.The normal force does not exceed 12 N at the beginning of the test (for obtaining a soil-to-plate contact) and tends to 0 N at the end of the test.The test duration is about 14 min, which leads to low water loss of the sample.Samples with water loss greater than 10 % gravimetric water content, measured as water loss from before to after testing, were disposed.Furthermore, the first three points of the 30 measurement points where the deflection angle (φ) is less than 1 μrad were excluded due to the lack of sensitivity of the rheometer used.
Rheological tests were controlled and executed automatically by the software Rheoplus/32 V3.62 Anton Paar.Deformation (γ, %) was determined by the relationship between the deflection (s, m) in the outer edge and height of the sample (h, m) (Equation 1); shear stress (τ, Pa) was calculated using the torque (M, m N) measured in each γ and the radius of the upper plate (r, m) (Equation 2); and phase shift angle (δ, °) was determined by displacement of the τ response curve in relation to the γ curve.
From these results, the storage modulus (G', Pa) and loss modulus (G", Pa) were calculated according to equations 3 and 4, respectively; and the loss factor (tan δ, dimensionless) was calculated by the ratio of G" and G' (Equation 5).
Rev Bras Cienc Solo 2016;40:e0150286 The G' and G" moduli, tan δ, and τ can be presented graphically as a function of γ (Figure 1), where the maximum τ during the test (τ max ) is obtained from the γ versus τ curve (Figure 1b).The deformation limit (γ L ) was calculated as the γ in which a difference of 5 % occurs in G' with respect to their values in the linear viscoelastic range (LVE range) (Figure 1a).The integral z (dimensionless) was calculated as the sum of the area defined, in the lower limit, by the tan δ curve from the lower applied deformation (here γ ≈ 0.0001%) until the yield point where tan δ = 1 and, in the upper limit, by the line parallel to the abscissa and with ordered tan δ = 1 (Figure 1a).
Other rheological parameters obtained from the amplitude sweep test, such as the τ at γ L , and the γ and the value of G'= G" at the yield point -described by Markgraf and Horn (2007) and Holthusen et al. (2010) -were not used in this study due to high collinearity with the parameters γ L , τ max , and integral z.

Statistical analyses
The set of rheological data was analyzed separately for each horizon, resulting in 16 experiments (16 horizons) in a completely randomized design, with water tension (w t ) as the quantitative factor with levels 0, 1, 3, 6, and 10 kPa.The results of these experiments were subjected to normality analysis by the Shapiro-Wilk test (p<0.05), to homogeneity of error variances by the Bartlett test (p<0.05),and to variance by the F test (p<0.05).When the F test was significant, simple regression analysis was performed.These analyses were performed using SAS 9.2 software (SAS, 2010).
Variables without normal distribution or homogeneity of variance were subjected to the Box-Cox transformation using Action software (Equipe Estatcamp, 2014).Spearman correlation analysis (p<0.05) between the rheological parameters and plastic limit (PL), liquid limit (LL), and plasticity index (PI) was also performed using the SAS 9.2 software (SAS, 2010).

RESULTS AND DISCUSSION
Soil water tension (w t ) did not affect the linear viscoelastic deformation limit (γ L ) of most horizons.There was a linear increase in γ L , with increasing w t only for the LVd2 A 1 , LBdf1 Ap, PVd1 Bt, SXe Btg, and VEo A horizons (Figure 2); i.e., the decrease in water content retained in the soil increased the range of elastic deformation of these horizons, where the most intense effect was for LBdf1 Ap.This result is partly due to the higher clay content of these horizons (Table 1).The greater the clay content and activity (such as smectite), the greater the expression of cohesion and adhesion forces (Reichert et al., 2010), increasing the contacts and linkages between particles with increasing w t (Kemper and Rosenau, 1984).Furthermore, the increase in clay-water attraction by van der Waals forces decreases the contact angle of water menisci with clay particles to near zero, thus increasing wetting and capillary forces (Lourenço et al., 2012;Amarasinghe et al., 2014).In sandy soils, such as in the top (A) horizons of LVd1, PVd1, PVd2, and SXe (Table 1), this phenomenon is less pronounced (Lourenço et al., 2012).
The highest values of γ L were observed in LVd2 A 1 , LBdf1 Ap, LBdf2 A 1 (without significant increase with increasing w t in the last horizon), and VEo A and Biv (Figure 2), which are associated with greater total carbon (TC) content (23 to 44 g kg -1 ) (Table 1).Organic compounds increase soil aggregation and resistance (Zhang and Hartge, 1990;Markgraf et al., 2012;Buchmann et al., 2015) by augmenting cohesion and number of menisci in unsaturated conditions (Zhang and Hartge, 1990;Bachmann and Zhang, 1991).Furthermore, with increasing w t , organic compounds and salts are precipitated in contact areas between particles and aggregates, increasing soil resistance to deformation (Soulié et al., 2007;Holthusen et al., 2010).The increase in soil mechanical stability with increasing SOM is in accordance with the findings of Markgraf and Horn (2007), who observed a decrease in γ L after removal of SOM by oxidation of Oxisols from RS.
In the same horizons (LVd2 A 1 , LBdf1 Ap, LBdf2 A 1 , and VEo A and Biv), higher γ L (Figure 2) are also associated with greater cation content (Table 1), corroborating results of Holthusen et al. (2010).An increase in polyvalent cations such as Al 3+ , Ca 2+ , and Mg 2+ decreases the repulsive forces between negatively charged particles by reducing the diffuse double layer (Singh and Uehara, 1998), bringing clay particles, oxides, and organic molecules to link and forming bridges that bond silt and sand particles (Bronick and Lal, 2005).Furthermore, with decreasing soil moisture, cations are strongly held in negatively charged clay particles and cations in excess or associated with anions may precipitate as salts (Baumgarten and Horn, 2013), which increases the number of bonds between soil particles.
Higher γ L also occurred in both horizons of VEo because of the significant presence of smectite (Albuquerque et al., 2000), an expansive 2:1 lattice type clay mineral that increases soil plasticity (Table 1), which allowed the soil to resist greater deformations (Khaydapova et al., 2015).In contrast, the kaolinite predominant in the clay fraction in the others soils of this study (445 to 638 g kg -1 of clay fraction, data not shown) collapses already at lower deformations (Khaydapova et al., 2015), although it establishes stronger links between particles (Kämpf et al., 2012).The results of Markgraf and Horn (2007), indicating higher γ L in a Vertisol compared to Oxisols, also show the collapse of structure at lower deformations of kaolinitic soils.In VEo, the adsorbed water and the cations in the soil solution, and on surfaces and in the interlayer of expansive clays increase soil elasticity (Mitchell and Soga, 2005;Markgraf et al., 2006), represented by γ L (Figure 2), due to the increase in capillary forces and electrostatic and molecular forces (Mitchell and Soga, 2005;Markgraf and Horn, 2006).According to Kemper et al. (1987), smectite particles saturated with Ca 2+ bind with sufficient strength to resist greater deformations.
For maximum shear stress (τ max ) there was a linear or quadratic increase with increasing w t in the surface horizons (except for PVd1) and subsurface soils (except for LVd2, LBdf2, and PVd2) (Figure 3), with greater shear strength in more clayey horizons   1).These results are similar to those obtained by Holthusen et al. (2010), who suggest that the increase in effective stress of the menisci, the cemented bonds due to precipitation of salts, and the number of contact areas between particles resulting from soil water drainage are mechanisms responsible for increase in τ max .Thus, besides the existence of adhesion forces on the particle-water-air interface, the water strongly adsorbed on the surfaces of minerals by hydrogen bonds or held by osmotic forces associated with cations in the diffuse layer tends to maintain the mineral particles in direct contact (Kemper et al., 1987), increasing the cohesive forces.This increase in cohesion with increasing w t tends to be higher in clayey soils than in silty and sandy soils, corroborating the results of Kemper and Rosenau (1984).As observed in our soils (Figure 3), increased mechanical strength with increasing depth in different soils is due to greater clay content in subsurface horizons compared to surface horizons (Table 1); this was also observed by Markgraf and Horn (2006) and Baumgarten et al. (2012).
For the LBdf1 Bw 1 and LBdf2 A 1 horizons, however, τ max increases until w t of 6 kPa and decreases for w t of 10 kPa.The point of maximum τ max at 6 kPa in these soils is possibly due to the combination of low ρ s (Table 1) with the occurrence of a strong and resistant granular microstructure, whose aggregates have a diameter similar to the diameter of the sand fraction particles, called pseudosand (visually observed and checked by tactile sensation according to Santos et al., 2005).This combination reduces the number of contact points within a given volume of soil and, consequently, reduces the number of menisci and the effective stress (Horn, 2003;Gallipoli et al., 2003;Reichert et al., 2010).The effective stress defines the forces that can stabilize the soil particles against deformations (Horn and Peth, 2011).In spherical particles, like sand and pseudosand aggregates, the capillary force is smaller than in colloidal particles (Chatterjee et al., 2012), because there is no continuous water film with a predominantly concave shape (negative water pressure) responsible for particles binding (Lourenço et al., 2012).A large formation of pseudosand in the Latossolos Brunos (Typic Hapludox) in this study is derived from high Fe oxide content, predominantly goethite (data not shown).
In structured samples of coarse texture (≈650 g kg -1 of sand), Holthusen et al. (2012b) observed an increase in τ max with an increase in w t up to 3 kPa and a subsequent decrease at 6 kPa, whereas in homogenized samples, the τ max of the same substrate increased from 0 to 6 kPa.According to the authors, the substrate aggregation formed pseudosand particles and, thus, greater drainage reduced the number of menisci more than the increase in the individual forces of the menisci.Similarly, Zhang and Hartge (1990) observed an increase in cohesion of a sandy soil with different SOM contents in the w t range of 4 to 7 kPa, and a decrease in cohesion occurred with the additional increase in w t up to 30 kPa.Greater resistance is observed when the soil is at intermediate water contents because the w t is sufficient to approximate particles and there is still water to displace ions and molecules for low-energy binding sites on the surface of the particles (Kemper et al., 1987).Among these forces, the stabilizing forces depend mainly on the water menisci between particles instead of electrostatic forces (Holthusen et al., 2010).
In fine-grained soils, therefore, increasing stability occurs with drainage because the phenomenon of adhesion occurs most intensely and the cohesion phenomena increase.However, in soils with a higher content of coarse particles, the adhesion and cohesion phenomena are less significant and the reduction of the menisci number is faster, due to the smaller number of contact points between particles (Mitchell and Soga, 2005) and, thus, a maximum point of resistance occurs in w t that is lower than in clayey soils (Holthusen et al., 2012b).The same behavior occurs in soils with lower density, as observed for Latossolos Brunos (Typic Hapludox), namely LBdf1 and LBdf2 (Table 1, Figure 3) because the increased spacing between soil particles increases the radius of w t (kPa) w t (kPa) w t (kPa) w t (kPa) w t (kPa) w t (kPa) Rev Bras Cienc Solo 2016;40:e0150286 curvature of menisci, which reduces the force of each menisci and, consequently, the capillary forces (Amarasinghe et al., 2014).
A linear increase in τ max up to the w t of 10 kPa occurs in the sandier horizons of LVd1 Ap and Bw 1 , PVd2 A 1 and SXe A 1 (Figure 3, Table 1); however, this seems contradictory to the process described before.A decline in shear strength would be expected with higher w t , but even sandy soils may exhibit increased cohesive force due to the surface tension of convex menisci formed in larger quantities in these soils (Lourenço et al., 2012).Higher ρ s (Table 1), predominance of fine sand followed by very fine to medium sand for all these horizons (data not shown), and some clay or silt content may justify the maintenance of capillary forces and the increase in τ max even with a suction of 10 kPa in these sandy horizons.LVd1 and PVd2 also have small amounts of Fe oxides (data not shown), which act as cementing agents (Six et al., 2004), increasing the rigidity of the links between soil particles.Furthermore, it is possible that increased friction between sand particles compensates the loss of capillary force in mechanical strength, although the effect of w t in sandy soil explains most of the variance in τ max (Holthusen et al., 2012b).Friction is an important part of shear strength for sandy soils due to high tensions at contact points between particles; whereas, for clay soils, the resulting cohesion of the interparticle attractive forces increases shear resistance (Reichert et al., 2010).
Integral z, like τ max , increased linearly or quadratically with increasing w t (Figure 4).In most horizons, including LBdf2 A 1 and Bw 1 (no statistical difference), structural rigidity increases up to at least the w t of 6 kPa.An increase in integral z in drained samples at w t up to 15 kPa, mainly in soils with finer texture, was reported by Baumgarten et al. (2012).Higher integral z values generally occur with increasing w t , and indicate a more elastic and rigid microstructure by considering the rheological curves as a whole, and is a parameter influenced by many soil properties (Baumgarten et al., 2012;Markgraf et al., 2012).
The decline in integral z in the quadratic regressions at w t of 10 kPa in the LVd1 Ap, LVd2 A 1 and Bw 1 , LBdf1 Ap and Bw 1 , and PVd2 A 1 horizons (Figure 4) seems to be related to particle size distribution and ρ s , similar to the observations made for τ max .Horizons with higher sand content (LVd1 Ap and PVd2 A 1 ) or pseudosand (LBdf1 Ap and Bw 1 , and LBdf2 A 1 and Bw 1 ), as well as the horizons with low ρ s -compared to the usual values of their textural class -(LVd1 Ap, LVd2 A 1 and Bw 1 , LBdf1 Ap and Bw 1 , LBdf2 A 1 , and Bw 1 ) (Table 1) show a decrease in structural rigidity when increasing the w t from 6 to 10 kPa, probably due to a decrease in the number of contact points between particles and capillary forces.
In both horizons of VEo, however, there was a decrease in integral z up to the w t of 6 kPa.This behavior, which is contrary to the results of γ L and τ max for VEo (Figures 2  and 3, respectively), do not show greater adhesion and cohesion forces expected in this soil with decreasing water content.Thus, additional studies to understand the causes of this behavior are necessary.The LBdf2 A 1 and Bw 1 , PVd2 Bt 1 , and SXe Btg horizons showed no significant regression of integral z with increasing w t (Figure 4), probably due to the sigmoidal behavior in the LBdf2 horizons and the small variations in water content between the w t evaluated for PVd2 Bt 1 and SXe Btg (data not shown).
Based on the consistency limits, most samples were in the plastic range (between LL and PL) (data not shown).There was a high positive correlation (Spearman r) between γ L (at all w t ) and PL, LL, and PI (from 0.63 to 0.85); and for τ max there was no significant correlation (data not shown).The correlation was positive with integral z and PL (0.71), LL (0.71), and PI (0.52) at saturated conditions, but was negative with PI at w t of 10 kPa (-0.69) (data not shown).
The positive correlation between γ L and IP indirectly indicates that soil elasticity assessed by rheometry (γ L ) increases with the increase in soil plasticity evaluated by the consistency test (IP), and both are associated with the clay and TC contents, presence of expansive clay minerals, and high cation content.Thus, in general, the evidence from this study indicates that the increased number of bonding points between particles and aggregates increases soil elasticity (γ L ); and shear strength (τ max ) and soil structural stiffness (integral z) increase with increasing strength or rigidity of these bonds, which are mainly dependent of the clay content and of the kaolinite and oxide content.
For our soils, rheological parameters indicating elasticity and structural strength, assessed with respect to w t , are largely in agreement with those observed in published studies.Increased soil structural strength with increasing w t occurs due to the decrease in the number of water molecule layers coating particles and aggregates, reducing their mobility and promoting additional resistance to deformation by increasing solid-solid interactions -i.e., greater cohesion (Ghezzehei and Or, 2001).Moreover, as a result of reduced water content, the menisci exhibit stronger contraction force -greater effective stress -even though the number of menisci decreases (Gallipoli et al., 2003;Peth and Horn, 2011;Holthusen et al., 2010Holthusen et al., , 2012b)).This is in accordance with the increase in structural strength in most horizons (Figures 2, 3, and 4), with increasing w t in this study, up to at least 6 kPa.
There is, however, some variation between different horizons, and soil resistance does not always increase with increasing w t , as was also observed by Holthusen et al. (2010).This suggests that clayey and sandy soils with high pseudosand content and, or low bulk density lose mechanical resistance through loss of capillary forces at w t above 10 kPa (field capacity).Given the importance of rheological parameters for assessment of changes in volume relationships and the functionality of soils under application of external loads as well as assessment of the risk of erosion and landslides, among other phenomena (Horn and Peth, 2011), more studies should be conducted to ascertain the decline in mechanical shear strength (τ max ) of clay soils at higher w t .

CONCLUSIONS
Soil mechanical strength increased with increasing water tension up to 6 kPa, primarily due to increased capillary forces in the soil.With increasing water tension, the increase in structural strength, as measured by τ max and integral z, was more pronounced than the increase in elasticity, measured by γ L .
Increased bonding between particles promotes soil elastic (or recoverable) deformation (γ L ), and γ L increased with higher amounts of clay, total carbon, expansive clay, and cations.Higher friction and strength of interparticle bonds increase shear strength (τ max ) and structural stiffness of the soil (integral z), and τ max and integral z mainly increased with increasing clay, kaolinite, and oxide content.
A decrease in mechanical strength of some soils at the highest water tension evaluated (10 kPa) was observed in the sandy horizons and horizons with a high proportion of very resistant microaggregates (pseudosand), particularly when such soils are associated with low bulk density.Thus, it is believed that this decrease in mechanical strength is caused by a smaller number of contact points between particles and consequent reduction in capillary forces.However, more studies are needed to evaluate the strength of granular soils with low bulk density and soils with expansive clay under the influence of water tension in the soil.

Table 1 .
Classification, symbol, coordinates, horizons, depth, and physical and chemical characterization of eight studied soils