Aggregate Stability in Soil with Humic and Histic Horizons in a Toposequence under Araucaria Forest

Aggregate stability is one of the most important factors in soil conservation and maintenance of soil environmental functions. The objective of this study was to investigate the aggregate stability mechanisms related to chemical composition of organic matter in soil profiles with humic and histic horizons in a toposequence under Araucaria moist forest in southern Brazil. The soils sampled were classified as Humic Hapludox (highest position), Fluvaquentic Humaquepts (lowest slope position), and Typic Haplosaprists (floodplain). The C and N contents were determined in bulk soil samples. The chemical composition of soil organic matter was evaluated by infrared spectroscopy. Aggregate stability was determined by applying increasing levels of ultrasound energy. Carbon content increased from the top of the slope to the alluvial plain. Higher ultrasonic energy values for clay dispersion were observed in the C-rich soils in the lower landscape positions, indicating that organic compounds play an important role in the structural stabilization of these profiles. Both aliphatic and carbohydrate-like structures were pertinent to aggregate stability. In the Oxisol, organo-mineral interaction between carbohydrates and the clay mineral surface was the most important mechanism affecting aggregation. In soils with a higher C content (Humaquepts and Haplosaprists), stabilization is predominantly conferred by the aliphatic groups, which is probably due to the structural protection offered by these hydrophobic organic groups.


INTRODUCTION
Soil aggregate stability is one of the most important factors for soil conservation and maintenance of soil environmental functions.Aggregation affects soil capacity to store and stabilize organic C (Six et al., 2004;Kodešová et al., 2008), as well as soil water storage capacity and distribution in the landscape (Scheer et al., 2011;Schmidt et al., 2011;Berhe and Kleber, 2013).In addition, an increase in soil structural stability increases resistance against erosive agents and compaction (Schmidt et al., 2011;Novara et al., 2012;Berhe and Kleber, 2013;Chaplot and Cooper, 2015).
Aggregate stability is defined by several intrinsic soil properties.In soil profiles with lower soil organic matter (SOM) content, the main aggregation agents are Fe and Al oxides.The effect of these minerals on aggregate stability has been attributed to structural cementation, and is a subject extensively explored in the literature (Goldberg et al., 1988;Seta and Karathanasis, 1996;Pinheiro-Dick and Schwertmann, 1996;Barthès et al., 2008;Igwe et al., 2009).However, in soils with lower clay and Fe oxide content, aggregate stability is related to SOM (Six et al., 2002).
The relation between C content and aggregate stability was observed in soils in tropical climates in different degrees of weathering (Chaney and Swift, 1984;Gajic et al., 2006;Six et al., 2004;Inda Junior et al., 2007;Li et al., 2007).The SOM effect on aggregation has traditionally been attributed to the inner sphere complexes between the carboxyl groups and cations of the mineral structure by the ligand exchange mechanism (Chorover and Amistadi, 2001;Mikutta et al., 2011).Others types of organo-mineral interactions are also proposed: i) hydrophobic interactions; ii) C-O-alkyl groups and mineral hydroxyls; iii) cation bridges; iv) anion and cation exchange; and v) Van der Waals interactions (Stumm, 1992;Dick et al., 2009;Hanke et al., 2015).However, the effects of these organo-mineral interactions on aggregate stability have not been investigated in detail.
In mineral soils, the positive effect of SOM on aggregation is usually related to its content (Six et al., 2002;Wiseman and Püttmann, 2005;Inda Junior et al., 2007;Kodešová et al., 2009).However, the role of SOM chemical composition on aggregate stability has been less investigated.
Soil organic matter is composed of hydrophilic structures, such as carboxyl, N, and C-O-alkyl groups, as well as hydrophobic structures, such as aliphatic and aromatic groups.Recent theoretical models hold that hydrophilic groups, which are abundant in oxidizing environments, interact directly with the mineral surface; this occurs mainly with Fe and Al oxides, due to their high density of monocoordinated hydroxyls (Krull et al., 2003;Schöning et al., 2005;Wiseman and Püttmann, 2006;Mikutta et al., 2011).Some studies have indicated that in aerobic soils, the surface of Fe oxides can interact directly with carbohydrate-like structures that are more labile compounds, and not only via ligand exchange with carboxyl groups (Miltner and Zech, 1998;Schöning et al., 2005).
The SOM zonal model (Wershaw et al., 1996;Kleber et al., 2007) proposes that organo-mineral interactions occur between the mineral surface and the organic hydrophilic groups, while the aliphatic and aromatic structures constitute sites for hydrophobic interaction with other SOM micelles.In Brazil, hydrophobic structures have been observed in hydromorphic profiles and in soils where SOM is less oxidized due to hydric saturation in the pedoenvironment (Pérez et al., 1998).It is believed that water repellency can be caused by particles being coated by hydrophobic substances.These substances may be derived directly from plant residue decomposition or from specific products of microbial metabolism (DeBano, 2003;Buczko et al., 2005).Possibly, these hydrophobic zones can aid aggregate stabilization by covering organo-mineral complexes and larger structures.
Due to SOM self-association, hydrophobic protection can be a very important aggregation mechanism in C-rich soils.This is the case for humic and histic soil profiles from the high Rev Bras Cienc Solo 2017;41:e0160369 plain zones of southern Brazil.These studies have been used to evaluate the effect on water retention and environmental filters (Kadlec and Wallace, 2009).Aggregation is one of the most important factors for these soil functions.These C-rich soils under natural vegetation are an interesting environment to investigate the factors related to SOM that control structural stability since the aggregate stability mechanisms have reached their equilibrium stage.Studies that relate SOM chemical composition to aggregate stability in C-rich environments are rare, yet may contribute to knowledge regarding how SOM acts on the phenomenon of soil aggregation.
The hypothesis of this study was that, under the same environmental conditions, the structural stabilization of aerobic soils is mainly related to the organo-mineral interaction, while in C-rich profiles under hydromorphic conditions, hydrophobic protection is the most relevant mechanism.Thus, the objective of this study was to investigate the aggregate stability mechanisms related to SOM chemical composition in soil profiles with humic and histic horizons in a toposequence under an Araucaria moist forest in southern Brazil.

Study area and sampling
Soil sampling was performed in a toposequence under an Araucaria moist forest in the sedimentary basin of Curitiba in the "Primeiro Planalto Paranaense" (roughly translated: "first plateau of Parana"), Brazil, located in the environmental protection area of Iraí (EPA-Iraí), a hydrographic microbasin of the Canguiri River.
The soils sampled were derived from sedimentary rocks (Argillite and Arcosian deposits) of the Guabirotuba Formation.This geologic formation belongs to the Cenozoic Era, with ages ranging from the Miocene to Pleistocene (Salamuni and Stellfeld, 2001;Salamuni et al., 2004).According to Köppen classification system, the climate is classified as Cfb -subtropical humid with moderate winters and mild summers, with annual average rainfall ranging from 1,500 to 1,600 mm and an average temperature of 19.1 °C (Maack, 2012).This region has a sub-dendritic drainage pattern and the slopes range from strongly rolling to rolling in a divergent convex topography.The slopes are attenuated at the top and at the entrance of the floodplain, where there is a flat slope.
During soil sampling, pH was determined in distilled water.In OX, the pH values ranged from 4.0 to 4.9.In FH, the same properties ranged from 4.5 to 4.9 and in TS, from 3.7 to 4.0.
To investigate aggregate stability, undisturbed soil monoliths in replicates of approximately 2 dm 3 were collected and transported to the lab in closed containers.Subsequently, the Rev Bras Cienc Solo 2017;41:e0160369 monoliths were gently disaggregated manually, maintaining field moisture, and sieved (9.51 mm).After air drying, the samples were separated into different aggregate classes using sieves with 4.00, 2.00, 1.00, and 0.25 mm mesh.

Organic carbon and nitrogen content
The C and N contents were determined by dry combustion ("Variol El" elemental analyzer) of the soil bulk samples.The samples had been milled and sieved through a 0.2 mm mesh.

Hydrofluoric acid treatment
Demineralization of the samples was performed according to Gonçalves et al. (2003).In brief, 1 g of soil was weighed and treated with 30 mL of hydrofluoric acid solution (10 %, v/v).The samples were shaken manually for 30 s, followed by horizontal mechanical agitation for 2 h.The material was centrifuged (10 min -relative centrifugal force = 1050) and the supernatant was discarded.This procedure was repeated eight times.At the end, the solid residue containing the concentrated organic matter (SOM HF ) was washed three times with distilled water and oven dried at 50 °C.

Fourier Transform Infrared Spectroscopy (FTIR)
The SOM HF samples were analyzed by Fourier transform infrared spectroscopy (FTIR) (Shimadzu 8300) in KBr pellets (1:100), using 32 scans and a resolution of 4 cm -1 in the spectral range of 4000 to 400 cm -1 .The absorption bands were attributed according to Farmer (1974) and Tan (1996).From the FTIR spectra we calculated the aromaticity index (AI) (Chefetz et al., 1996) in which I C=C is the absorption intensity around 1630 cm -1 and I C-H is the absorption intensity around 2920 cm -1 , after establishing the baseline between 1800 and 1500 cm -1 and between 3000 and 2800 cm -1 .

Aggregate stability and dispersed clay content by ultrasound energy application.
To investigate soil aggregate stability by ultrasonic energy application, approximately 10 g of soil were recomposed from the mass proportions of each aggregate class in the bulk sample.The soil aliquots were moistened by capillarity and allowed to stand for 24 h.These samples were placed in 100 mL centrifuge glass tubes, to which 80 mL of deionized water was added.The soil mass was converted into volume by particle density, which was determined by the fact that soils with high SOM content tend to show lower particle density values.The method used to determine particle density was described in Claessen (1997).Particle density values observed were higher in OX (mean of 2.63 mm), intermediate in FH (mean of 2.49 mm), and lower in TS (mean of 1.81 mm).
Subsequently, increasing levels of ultrasound energy (0,50,100,200,400,600,800,1000,1200, and 1500 J mL -1 ) were applied to the suspension containing soil and deionized water.The equipment used was a Vibracel VCX 700 previously calibrated according to the method described by Christensen (1992).
The concentration of dispersed clay after energy application was determined by the pipette method (Gee and Bauder, 1986).The data of dispersed clay as a function of the energy level applied were fitted to the equation y = a (1 -e -bx ), in which "y" is the Rev Bras Cienc Solo 2017;41:e0160369 dispersed clay concentration and "x" is the energy level.This exponential model shows that the clay dispersion rate is decreasing and that the dispersed clay content tends to the maximum, represented by parameter "a" of the function, in which the amount of ultrasound energy applied (variable "x") tends to infinity.The energy for total dispersion (E max ) of the aggregates was assumed to be required to obtain a dispersed clay content equivalent to 99 % of the parameter "a" value.The energy level for each sample was considered an indicator of aggregate stability.

Water holding capacity in disturbed soil samples
Water holding capacity (WHC) was determined to perform correlation analysis with SOM chemical composition in order to identify changes in the degree of sample hydrophobicity due to hydrophobic groups.These results were used to support the effect of SOM on aggregate stability.
Approximately 10 g of soil was recomposed from the mass proportions of each aggregate class in the bulk samples.These samples were milled, sieved (2 mm), and air dried.The soil was placed in filter paper in a glass funnel in a pre-weighed flask.
Then, about 100 g of distilled water (previously weighed on an analytical balance) was added in small portions to the sample.The funnel was covered with plastic film to avoid water loss by evaporation.After 12 h, the water flask was weighed, and WHC was determined by the quotient for the retained water mass over the dry soil mass.

Statistical analysis
Results were first analyzed by descriptive statistics (mean and standard deviation) to verify natural variability of samples with respect to three replicates.After verifying that, in general, the natural variability of the sample set was low, the means were analyzed by multivariate statistics.To support the multivariate studies, simple linear regressions were also applied.
The multivariate analysis was divided into two parts and the statistical software used were i) MULTIV (beta version 309 for Windows); ii) Statistica 7.0; and iii) BioEstat 5.0.The first step was to verify the numerical distribution of the variables and clusters between sample units in order to investigate the need for sample treatment as separate groups.To do so, principal component analysis (PCA) was applied.The data matrix involving all variables (soil chemical and physical properties and SOM variables) was subjected to vector transformation between variables by a standardized dispersion.Correlation was chosen as a similarity measure.
Subsequently, the dispersion stability of the ordination diagram was verified by a bootstrap test.After confirming data dispersion stability, cluster analysis was applied.Group stability was verified by reapplication of the bootstrap test with 5 % significance (p<0.05).As will be presented later, each of the three soils was considered a distinct group and was treated separately in the second step of multivariate analysis.
The objective of the second step was to investigate the effect of SOM chemical composition on aggregate stability.For that reason, canonical correlation analysis (CCA) was applied.The advantage of this technique was its ability to identify the strongest associations between two variable sets simultaneously, rather than simply calculating simple pairwise correlations (Gittins, 1985;Recio-Vazquez et al., 2014).A smaller number of variables than sample units is required for use of this technique.
As the soil organic C content showed a correlation with aggregate stability (results presented and discussed later), only the RI of the SOM groups (identified by FTIR) were chosen as matrix predictor variables; as predicted variables, the E max values were chosen.
Rev Bras Cienc Solo 2017;41:e0160369 The models examined were constructed from two data matrices that produced two canonical variables: i) a linear combination from the RI of SOM groups; and ii) a linear combination of E max values.Thus, for each soil, a canonical function (canonical variable) was obtained.This canonical variable describes the aggregate stability as a linear function of the SOM groups.
After this procedure, the RI values of the SOM group for each soil were replaced in the canonical function obtained.These values were correlated with the original E max values.The objective of this procedure was to evaluate the coefficient of determination (R 2 ) and the regression significance between the canonical function and aggregate stability.
Finally, the correlation coefficient (r) and the significance of the RI of each organic group with each canonical function were analyzed separately.The objectives of this procedure were to identify the most important chemical groups to obtain the canonical function, as well as to rank the groups according to the R 2 and its significance.

C and N content
The C and N content ranged from 0.9 to 176.2 g kg -1 and from 0.1 to 11.1 g kg -1 , respectively (Table 1), and both properties decreased in the order TS > FH > OX.In TS, the C content of all layers was above 80 g kg -1 , confirming the occurrence of histic horizons in this profile.In FH, the C content ranged from 67.4 to 90 g kg -1 , indicating that this soil is a transition between the mineral soils of the slope and the organic soils of the floodplain.
In OX, the C content ranged from 0.9 to 56 g kg -1 , decreasing progressively from the top to the deepest layers (Table 1).In the A horizon (0.00-0.42 m), the properties C content, horizon thickness, color, and base saturation characterize the surface horizon as A humic.The subsurface layers are related to the mineral horizons (Bw and C).As expected, the C and N contents were lower in the subsurface layers than the humic horizon (Table 1).The increase in C contents from the top of the slope to the floodplain confirms the conversion gradient of humic content in histic horizons, with increasing soil water saturation.
The C/N ratio was generally higher in FH (17-24), followed by TS (15-19) and OX (6-18) (Table 1).Changes in the values of the C/N ratio by depth suggest the occurrence of different SOM dynamics in these landscape positions.In OX, the C/N ratio increases up to 0.60 m, reaching 18, and then decreases in depth until C/N = 6 (1.60-1.80m).In FH, the C/N ratio is, in general, higher than in the other soils, and the C/N values are similar between the different layers.In TS, the C/N ratio tends to increase with depth, indicating relative enrichment of N compounds in the soil surface.This difference is probably related to water regime changes in different landscape locations.The SOM decomposition is more limited in hydromorphic terrains than the more oxidizing conditions at the top of the slope.

SOM chemical composition by FTIR
The SOM HF samples showed the same FTIR spectrum pattern; in figure 1, the FH spectra are shown as an example.The absorption bands observed in the spectra were 3420 cm -1 , attributed to the O-H stretching in H-bonds; two bands at 2950 and 2840 cm -1 , attributed to the C-H stretching of aliphatic groups; 1710 cm -1 , attributed to the C=O stretching of carboxyls; 1630 cm -1 , attributed to C=C vibration of aromatic groups; 1540 cm -1 , attributed to N-H deformation and to C=N stretching; 1245 cm -1 , attributed to C-O stretching and OH deformation of the carboxyl group; and 1075 cm -1 , attributed to the C-O bond of primary and secondary alcohols.In SOM studies, this organic function is considered a part of the carbohydrate-like structures (C-O-alkyl).
In TS, the RI 1630 ranged from 24.8 to 30.4 %, with similar values throughout the profile (Table 1).In this profile, the AI ranged from 2.8 to 4.4, and the higher values at depth may be the result of older aromatic preservation by hydromorphy.
The RI of carboxylic groups (RI 1710 ) ranged from 10.3 to 27.3 %.In general, these values were higher in OX (21.5±2.2) and TS (24.0±1.7)than in FH (16.4±5.5)(Table 1).The OX profile showed a small tendency to increase RI 1710 at depth, which is probably related to carboxylated hydrophilic compound percolation (Kalbitz et al., 2000).In TS, hydromorphic conditions may limit the decarboxylation process, which may be considered one of the most important processes of SOM oxidative decomposition.The lower redox potential of hydromorphic soils may limit carboxyl breakdown from the structure and, consequently, its conversion to CO 2 .Under these conditions, other biochemical fermentation pathways are used by microorganisms (Nelson and Cox, 2002).
In FH, poor drainage and water table oscillation limit carboxylic compound percolation to deeper soil layers.
In OX, the RI of the N compounds (RI 1540 ) ranged from 6.0 to 12.2 % and increased gradually at depth (Table 1).This increase was associated with a decrease in the C/N ratio, confirming relative enrichment in N structures, which suggests an increasing degree of humification.In TS, the RI 1540 oscillated around 11±1.3 throughout the profile.The FH profile showed IR 1540 values much lower (2.6±0.4)than those observed in the other soils, indicating a lower proportion of N compounds.These results are consistent with higher C/N values compared to the other profiles.
The OX soil showed the highest values of RI 1075 , which decrease from the surface to depth (Table 1).The higher abundance of these groups, with higher biochemical lability in the SOM composition of the surface horizon, is related to the contribution of vegetation biomass.At depth, the proportion of this group decreased, due to its preferential degradation in the aerobic environment.In the FH and TS samples, a positive linear correlation was observed between RI 1075 and clay content (Figure 2).This result suggests that the C-O groups that are part of carboydrate-like structures (primary and secondary alcohols), which are more polar than other C-O groups, interact with the polar sites of mineral surfaces.This hypothesis is supported by the higher iron extracted by ammonium oxalate and iron extracted by dithionite-citrate-bicarbonate (Feo/Fed) ratio of FH and TS in comparison to OX (Hanke and Dick, 2017).In these two C-rich soils, the low crystallinity Fe oxide content is higher than in OX and, therefore, they have a higher density of monocoordinated hydroxyls.
The reactive sites may interact with C-O-alkyl groups of SOM.The interaction between the carbohydrate-like structures and the clay-mineral surface has been observed by other authors (Schöning et al., 2005;Hanke et al., 2015).This result corroborates the existence of a contact zone between SOM and the mineral surface through the action of hydrophilic organic groups (Kleber et al., 2007).

Clay content and aggregate stability
The clay content was similar between OX (303 to 497 g kg -1 ) and FH (320 to 487 g kg -1 ), and lower in TS (123 to 187 g kg -1 ) (Table 1).Sample dispersion, in ultrasonic energy, showed an inflection range of dispersion curves.This range was followed by practically constant values of dispersed clay, indicating the destruction of soil aggregates.Examples of the dispersion curves can be found in figure 3.
The E max values of samples ranged from 33 to 1535 J mL -1 (Figure 4), with lower values in OX and higher in FH and TS (Figure 4).The E max range was lower for soils with higher mineral content (687 and 555 J mL -1 for OX and FH, respectively) and higher for TS (952 J mL -1 ).However, in OX, the E max decreased progressively at depth (from 720 to 33 J mL -1 ), with low values in the deepest layer (1.60-1.80m) and an energy level expressively lower than that observed in the surface layer (0.00-0.05 m).This result highlights the structural fragility in subsurface horizons in comparison to the humic horizon (0.00-0.42 m), in which E max showed a mean value of 529 J mL -1 .
In contrast to OX, FH and TS showed E max values ranging randomly between the different layers (Figure 4).In general, these values were high, characterizing high soil structural stability by ultrasound energy.
The SOM effects on aggregate stabilization may be observed in the positive correlation between C content and E max (Figure 5).The slope of the regression line (angular coefficient) decreases in the order FH (24) > TS (15) > OX (13), indicating that C content affects aggregate stability differently in all three soils.In OX, the proportion of C to the clay content is low (C/clay ratio values ranging from 0.002 to 0.17), suggesting the occurrence of a direct interaction between the organic groups and mineral surface (contact zone).In TS, the C/clay ratio was higher, ranging from 0.54 to 1.34, which is due to the low clay content of this profile compared to the other soils.This result suggests C saturation of the mineral surface, as well as that the SOM effect on aggregation is different from that in OX.In FH, the C/clay ratio values were intermediate, ranging from 0.15 to 0.22, indicating that the SOM effect on aggregation may be a combination of the main processes that occur in OX and TS.In contrast, it may be observed that the aggregate stability of FH and TS is extremely dependent on C content, since the regression intercept values (parameter "a") are negative.This means that the "ŷ" value tends to zero when "x = 0".In OX, the intercept value is positive (83), which theoretically represents the action of other soil properties (chemical and mineralogical) on aggregate stability.
In OX, no correlation was observed between clay content and E max .This suggests that the difference in C content and SOM chemical composition between the layers was more substantial in explaining the differences in structural stability than the textural aspects.
In addition, lower aggregate stability was found in the layers with higher clay content in the subsurface.
The absence of a correlation between clay content and aggregate stability was also observed by other authors in studies on weathered soils.This has been attributed to the more pronounced effect of certain clay surfaces on aggregation (Neves et al., 2006;Inda Junior et al., 2007).
In FH, as opposed to OX and TS, the clay content was important in explaining the E max range between soil layers (positive correlation: R 2 = 0.88 and p<0.01).The uniformity of C content between FH layers may have contributed to the texture effect on aggregate stability.In TS, high C content overlaps with the clay effect on aggregation.

Water holding capacity in disturbed soil samples
The WHC of soil samples was higher in TS, with values ranging from 3 to 6 g g -1 (Table 1).
In FH, WHC ranged from 0.9 to 1.6 g g -1 , and in OX, it was from 0.40 to 0.50 (Table 1).These results confirm the high water storage capacity in histic horizons (TS), as well as in mineral-organic transitional soil (FH).
The values observed in the present study are consistent with those reported for organic horizons by Soil Survey Staff (1992), where the sapric material has a storage capacity of about five times its dry mass.This capacity is lower for pedogenic mineral materials.Therefore, in addition to the potential for C-stock, these soils have high water retention, which is important for water cycle regulation and water distribution along slopes.
Although higher WHC values (TS and FH) were observed in the hydromorphic C-rich soils, no direct correlation was observed between WHC and C content.This suggests that in soils with a high SOM content, C content is not the determining factor, rather SOM chemical composition is.In FH and TS, there is a negative correlation between the RI of the aliphatic groups (RI 2950-2840 ) and WHC (Figure 6), indicating that, while these C-rich soils showed higher water storage capacity, the increase of SOM hydrophobicity leads to a decrease in WHC.
The relative aliphatic enrichment of SOM possibly leads to hydrophobic niche formations that repel water entry.In turn, this higher water repellency would contribute to SOM preservation.This occurs by isolation of more labile groups by hydrophobic compounds.This hydrophobic protection has higher biochemical recalcitrance, limiting the action of heterotrophic organisms.Hydrophobic encapsulation also prevents structural collapse by the instantaneous exit of internal air from the aggregate in the presence of water.
For that reason, this hydrophobic mechanism could contribute to aggregate stability and lead to SOM preservation inside these structures (Ellerbrock et al., 2005).Thus, there is a synergism between aggregate protection by hydrophobic structures and SOM preservation by physical occlusion of organic material.
In OX, this correlation was not observed, possibly due to the greater importance of surface adsorption mechanisms to water retention in the clay fraction (Figure 6).In this case, contrary to that observed in FH and TS, the contribution of the clay fraction to WHC overlaps with the contribution of SOM chemical composition.In FH and TS, although the increase in hydrophobicity may contribute to SOM preservation through biochemical recalcitrance, as well as an increase in aggregate stability, this characteristic may reduce the WHC (Ellerbrock et al., 2005).
Our results showed that for low SOM content, the WHC is governed by clay content.Although no direct relationship was observed, the increase in SOM content leads to an increase in WHC, which decreased in the order TS > FH > OX.However, for the C-rich environments, the WHC range seems to be conditioned by the SOM chemical composition.That is, the relative enrichment of aliphatic groups leads to the formation of hydrophobic niches that repel water entry and decrease WHC.The high WHC of hydromorphic soils, associated with higher C contents, emphasizes its importance in the maintenance of very specific environmental conditions.These environmental conditions act as important support for the local fauna and flora that need to be preserved.

SOM and aggregate stability -principal component analysis and canonical correlation
In order to verify the sample unit distribution and group formation in relation to SOM composition and aggregate stability, principal component analysis (PCA) was performed in association with cluster analysis.The SOM variables (C and N content, C/N ratio, and the RI of SOM groups) and E max values were used for this purpose.By the bootstrap resampling test, only three main components (axes) were considered valid (p<0.05) and approximately 91 % of total variance could be explained by these components.Of this proportion, about 79 % of the variance is described by components 1 and 2, that is, 72 % of total variance (Figure 7).
Cluster analysis identified the group formation from the sample units for each soil.This result corroborates the sample dispersion pattern in the ordination diagram (Figure 7).However, a greater similarity was observed between the sample of the OX surface horizon and FH samples, which is mainly due to higher C and N content, higher C/N ratio, higher RI of aliphatic and nitrogen groups, and higher aggregate stability.In turn, the OX samples formed an isolated group, strongly described by the variables E max , C-alkyl, C and N, carboxyl, and C-O-alkyl.The bootstrap resampling test performed in cluster analysis confirmed the stability of these three groups (p< 0.05).Thus, each soil class of a toposequence characterizes a specific environment and, for the purpose of this study, was analyzed separately.
Considering that the C content showed a positive correlation with the aggregate stability parameter (Figure 5) and in the interest of studying only the relationship between the SOM chemical composition with this property, we used canonical correlation analysis (CCA).In this analysis, the RI of different functional SOM groups were considered as an independent variable set, and the E max as a dependent variable.
The description of canonical functions, their significance, and their correlation with aggregate stability for each soil and the statistical parameters of each regression model are presented in figure 8.The weight of each independent variable in the canonical function is presented in table 2.
The canonical predictor variables (canonical functions), composed of the RI of the SOM groups, were significant (p<0.05) for the three soils analyzed.The correlation between these variables and E max was likewise high.However, only some organic groups were considered valid considering the correlation coefficients between each group and the canonical function (r > 0.50).Based on the coefficient of determination (R 2 ), the significant organic groups were ranked for each canonical function (Table 2).
In OX, the two most important organic groups were C-O-alkyl and C-alkyl; it should be noted that the effect of carbohydrate-like structures was higher than aliphatic chains for soil aggregation (Table 2).A similar response was observed in FH, in which the same groups were considered relevant for aggregate stability.These results suggest that in mineral soils, the importance of more labile and polar structures, such as C-O-alkyl groups, is more pronounced due to the higher clay content.The reactive hydroxyls on the clay-mineral surface may interact with SOM polar groups forming organo-mineral complexes.In a second mechanism, the C-alkyl structures would self-associate, forming a second zone of hydrophobic interactions (Kleber et al., 2007).This hydrophobic zone may contribute to the increase of structural stability by the physical protection of aggregates.In FH, both C-O-alkyl and C-alkyl groups were important for aggregate stability.This was probably due to the fact that FH is a transitional condition between organic and mineral soils.However, the higher R 2 of aliphatic structures (Table 2) suggests there is a more significant contribution of a hydrophobic protection mechanism than in OX.Probably, the  higher C content and lower clay content compared to OX favor mineral surface saturation by C, as well as SOM self-association.
In TS, only the C-alkyl structures were significant for determining the canonical function.
In this soil, although a significant correlation was observed between the clay content and RI of C-O-alkyl groups (Table 2), the carbohydrate-like structures were not relevant in explaining aggregate stability.This result is due to the high C content and the lower density of mineral sites (mean of clay content = 156 g kg -1 ) for interaction with polar SOM groups, favoring mineral surface saturation by SOM and aliphatic association.Thus, the protection conferred by the hydrophobic structures (C-alkyl) was the predominant mechanism in soil structural stability.
Our results suggest that there is a transition between structural stability mechanisms according to toposequence position: i) in the non-hydromorphic positions, where soils are weathered and have lower C content, the organo-mineral interaction is more expressive; and ii) in hydromorphic positions (lower third and floodplain), the hydrophobic physical protection is the most important mechanism.However, this does not mean that in FH and TS strong interactions between SOM groups and lower crystallinity Fe oxides cannot occur.Ferridrite and Lepidocrocite are typically found in C-rich environments, as well as in flooded soils.In fact, the Feo/Fed ratio in FH and TS ranged from 0.46 to 1.0.The same ratio in OX showed very low values (Hanke and Dick, 2017), confirming the higher proportion of lower crystallinity oxides in C-rich soils.However, this mechanistic contribution may be less important than the protection conferred by the aliphatic zone, which was confirmed by the CCA.
The existence of a transition between the aggregate stability mechanisms in the toposequence, related to the polar and apolar SOM compounds, can be supported by the correlations between the WHC and the E max in each profile.Whereas in OX there was no significant correlation between these two attributes (R 2 = 0.02), in FH and OX, this correlation was significant (p<0.01).
One of the main mechanisms that break the aggregates is water entry, which causes instantaneous air expulsion, collapsing the structure (Sullivan, 1990).In addition, separation between the components by the increase in the water layer is also responsible for the rupture.This is due to cohesive force reduction and particle adhesion, causing the component dispersion.Therefore, the existence of stable organo-mineral complexes and the protection conferred by hydrophobic groups are important for aggregate stabilization.

CONCLUSIONS
Both C-alkyl and C-O-alkyl groups are important for aggregate stability as determined by ultrasound energy.
In the more weathered soil (OX), organo-mineral interactions between the C-O-alkyl groups and the clay surface seem to be the most important aggregation mechanism.
In C-rich soils of lower toposequence positions (FH and TS), the predominant stability mechanism seems to be the physical protection conferred by hydrophobic structure self-association involving C-alkyl groups.
Higher ultrasonic energy values for clay dispersion are observed in the C-rich soils (TS and FH), indicating that organic compounds play an important role in structural stability.
Considering that changes in soil use and management may be reflected in changes in SOM chemical composition, a decrease in the structural stability and C-stock capacity of these soils is expected.Therefore, future studies comparing natural and anthropic systems under the same environmental conditions could contribute to understanding the influence of human activities on aggregation mechanisms related to SOM.

Figure 5 .
Figure 5. Linear regressions between the organic C content and dispersion energy levels in a) OX, b) FH, and c) TS.

Figure 6 .
Figure 6.Linear regressions between water holding capacity and a) clay content in OX, b) relative intensity of C-alkyl in FH, c) relative intensity of C-alkyl in TS.Clay (g kg -1 ) 0 300 400 500 Water holding capacity (g g -1)

Figure 8 .
Figure 8. Linear regressions between the canonical variable and the ultrasound energy levels in a) OX, b) FH, and c) TS.

Table 1 .
Carbon and N content (mean±standard deviation), FTIR spectral relative intensities, aromatic index, clay content, and water holding capacity of three soils in a toposequence under Araucaria moist forest -Pinhais, PR, Brazil

Table 2 .
Canonical variables produced by SOM chemical groups, canonical coefficient of determination, and significance and ranking of organic groups for each canonical variable *: significance (p=0.05).