Soil quality indicators for monitoring the short-term effects of mined soil rehabilitation strategies for bauxite

ABSTRACT Mining is a significant driver of soil mobilization, which impacts its physical, chemical, and biological properties. Changes in land-use affect the distribution of organic matter fractions in stable aggregates, a process that is still poorly understood, especially in drastically altered areas. Recovering and monitoring soil quality to ensure the sustainable development of agricultural crops in these areas after mining is challenging. This study aimed to evaluate the influence of agronomic practices in soil rehabilitation in a bauxite-mined area after three years of field experiment installations through an assessment of organic properties in soil and aggregate classes; an attempt was also made at proposing and elaborating a Soil Quality Index (SQI), which encompasses the soil’s physical, chemical, and organic properties. Different combinations of fertilization treatments and ground cover plants intercropping with coffee were evaluated as rehabilitation practices. The results showed that after three years of rehabilitation, when organic (OF), chemical (CF), and OF+CF fertilizers were applied to the areas of coffee intercropped with Brachiaria (B), they provided higher C and N contents to the soil and aggregates classes, as well as the compartments of soil organic matter (SOM). The minimum set of soil quality indicators for reclaimed bauxite-mined areas was composed of organic indicators: labile organic carbon (LOC) and mineral-associated organic matter (C-MOM); chemical indicators including pH and effective cation-exchange capacity (t), and physical indicators such as the bulk density (BD) and stable aggregates index in water (SAIW). The t and pH were the variables most sensitive to the management systems implanted during the rehabilitation of the mined area, and, therefore, were considered the best indicators of soil quality. Brachiaria was the cover plant that contributed most to improving the soil quality of mined bauxite areas by increasing the SQI, especially when fertilized. In general, when applied to the Brachiaria, the OF+CF fertilization presented a SQI of 0.78, differing statistically from that of the natural vegetation (1.00). Fertilizers and cover crops in association with coffee in the bauxite-mined areas improved the physical, chemical, and organic properties of the soil, thus representing a viable option for reconditioning mineral exploration areas.


INTRODUCTION
Mineral extraction represents an important activity contributing to the economic development of a country.Mining affects the sustainability of natural resources, especially that of soil and water.Studies about environmental impacts caused by mining activities are of increasing environmental relevance.
The Brazilian mineral industry has helped generate important surpluses in the country's trade balance, having exported more than 403 million tons of mineral goods in 2017.This generated foreign exchange on FOB (Free On Board) US$ 28,3 billion, representing 13 % of Brazil's total exports and 30.5 % of the trade balance (IBRAM, 2018).Tonietto and Silva (2011) reported that 67 % of the ore reserves measured and indicated in Brazil are located in the Minas Gerais State.According to statistics from the National Department of Mineral Production (NDMP), Minas Gerais is the Brazil's second largest producer of bauxite and, in 2016, produced the equivalent of R$ 112,199,137.00 (DNPM, 2018).
Soil degradation reduces agricultural production capacity resulting from its intensive use, causing changes in soil's physical, chemical, and biological properties.In addition to deforestation, the removal of topsoil exposes the remaining soil to erosive processes.Moreover, the withdrawal of consolidated material exposes the surrounding communities to dust and noise pollution.Therefore, mining activities in Brazil are positive for the municipalities, not only due to the taxes collected, but also through the direct and indirect jobs they generate.Rehabilitation efforts in mining areas have been conducted by companies, environmental organizations, universities, and research institutes that have been seeking effective procedures to reestablish the essential processes of the soil and affected ecosystems (Carneiro et al., 2008).
Rehabilitation involves the process of reversing the damage affecting degraded lands and turning the latter into productive and self-sustaining lands (Gripp and Nonato, 1993).The Federal Constitution of Brazil, Decree 97.632/89 -art.3, establishes that the degraded areas must aim to return to a form of use to achieve stability of the environment, according to a pre-established plan for land-use.Rehabilitation practices that provide a higher contribution of soil organic matter (SOM) are more successful because it can improve soil properties.They are also important for maintaining the soil quality of the ecosystems (Santos et al., 2011;Pillon et al., 2011;Cavalcante et al., 2019).
One of the indicators of the highest importance in the evolution of soil quality is its organic matter content (Doran and Parkin, 1994;Canals et al., 2007).Organic matter helps to indicate the quality of the management system that should be adopted for the ecological recovery of degraded areas (Bao et al., 2017).The decrease of C and N contents, which are considered essential propertes to the proper functioning of soils, has been reported in papers regarding the evaluation of areas under rehabilitation after mining (Trindade et al., 2000;Carneiro et al., 2008;Lunardi Neto et al., 2008;Cavalcante et al., 2019).Studies conducted on bauxite mining soils in Minas Gerais have shown a reduction of 99 % of the total C and N contents, microbial biomass, and enzymatic activity of the soil (Carneiro et al., 2008).Some strategies for the rehabilitation of mined areas have been used to increase the SOM through their revegetation by grasses and native plants (Vickers et al., 2012;Kneller et al., 2018;Iskandar et al., 2022), the use of leguminous plant species (Mukhopadhyay and Masto, 2016), that besides of C, improves N content in the soil (Roberts et al., 2015;Rodríguez-Vila et al., 2016).An increase or decrease in the SOM due to soil management practices implies benefits or damage to the system, respectively (Canals, 2007).
As a soil quality indicator (SQI), Reeves (1997) reported that the soil organic carbon (SOC) was a fundamental SQI in long-term studies.The SOC frequently does not present a high sensibility to the changes caused by determinate management systems in the short and medium terms; therefore, different SOM fractions can instead be used to verify the susceptibility of soils to the management types employed, according to the degree of lability that present (Gregorich et al., 1988;Pinheiro et al., 2004;Conceição et al., 2008Conceição et al., , 2014)).The study of soil organic compartments using SOM fractionation has provided a better understanding of nutrient releases and organic matter stabilization under changes in soil management and the environment (Lima et al., 2008;Alcântara Neto et al., 2011;Pessoa et al., 2012).
Different properties such as cation exchange capacity (CEC), soil biological activity, soil structure, pH, rainwater infiltration, soil erosion, salinization, etc., are all related to soil quality.However, the SOM compartments obtained by physical fractionation have not yet been used to evaluate soil quality, although these fractions are of great importance for the monitoring and evaluation of edaphic changes caused by the mining process.
Considering the reconditioning of the physical, chemical, and organic soil properties of mined areas, this study used organic and chemical fertilizers, as well as cover crops to: (i) determine the C and N contents in soil, soil organic fractions, and aggregate classes; (ii) elect the minimum number of indicators of soil quality encompassing physical, chemical, and organic properties; and (iii) use quality indicators to investigate the effects of different fertilization procedures and cover crops on coffee cultivation on the recovery of bauxite-mined areas.

Study area and experimental details
The study was carried out in the bauxite-mined area located in the Minas Gerais Forest Zone, in the municipality of São Sebastião da Vargem Alegre, Minas Gerais, Brazil (W 42° 35' 02.18" and S 21° 01' 58.98") (Figure 1).The ore extraction concession belongs to the Brazilian Aluminum Company -Votorantim Metals (Companhia Brasileira de Alumínio -Votorantim Metais).The terrain of the region is strongly rolling and the soil is classified as Latossolos Vermelho-Amarelos Distróficos típicos, according to the Brazilian Soil Classification System (Santos et al., 2018), which corresponds to an Oxisols (Soil Survey Staff, 2014).The climate is defined as high-altitude tropical (Cwa) according to the Köppen-Geiger classification system (Köppen, 1936), with a mean annual temperature and rainfall of 23.5 °C and 1,564 mm, respectively.
After three years of experiment implementation (Figure 2), the rehabilitation of the mined area was first assessed by Cavalcante et al. (2019).A randomized block design (RBD) experiment with a split-plot arrangement was used with four replicates, in which the treatments consisted of coffee cultivation, different types of fertilization procedures, and different cover crops between the rows: no fertilization (NF); organic fertilization (OF); chemical fertilization (CF); organic + chemical fertilization (OF+CF); no plants (NP); stylosanthes (S); Brachiaria (B); and an area covered by a fragment of native vegetation (NV) was considered as reference area.The specifications of the treatments are described in table 1.Further details can be found in Cavalcante et al. (2019).

Soil sampling and physical analysis
All soil samples were collected from 0.00-0.20 m soil layer, and their physical, chemical, and organic properties were analyzed.Soil physical propertiess were determined following the methodology proposed by Donagemma et al. (2011).Bulk density (BD) was obtained with the volumetric ring method using undisturbed soil.All other soil analyses were determined on <2 mm soil fractions.Clay, silt, and coarse and fine sand were obtained Rev Bras Cienc Solo 2023;47:e0220126 by physical dispersion with slow agitation for 16 hs (50 rpm) and chemical dispersion with NaOH 0.1 mol L -1 ; the sand contents were determined by sieving and the amounts of clay and silt were quantified with the pipette method.The total porosity (TP) was determined by the following equation: Eq. 1 in which: ρ s is the particle density obtained by the volumetric balloon method; microporosity (Micro) was determined by weighing after equilibrium at -6 kPa in a tension table, and macroporosity (Macro) was defined as: TP -Micro.
Wet separation of aggregates was performed by a Yoder type vertical oscillating apparatus (26 cycles min -1 and 5.0 cm vertical amplitude for 15 min), resulting in six stable aggregate classes: 4.75 -2.00; 2.00 -1.00; 1.00 -0.50; 0.50 -0.25; 0.25 -0.105; and <0.105 mm (Kemper and Chepil, 1965).The proportions of water-stable aggregates (WSA) were used to calculate the mean weight-diameter (MWD W ) and the geometric mean diameter (GMD W ) using the following expressions:  For the calculation of indices (MWD W , GMD W , SAI W , and AGRI W ), the amount of particles >2.0 mm (grit) was discounted from the dry weight of aggregates retained in the class 4.75 -2.00 mm (Castro Filho et al., 1998).

Soil chemical analysis
Among the chemical properties, the pH in H 2 O was determined potentiometrically using a soil:water (1:2.5)suspension.Phosphorous and K were extracted by Mehlich-1 and determined by molecular absorption and flame photometry (Defelipo and Ribeiro, 1997).

Analysis of the organic properties
In the analysis of organic properties, the total organic C and N contents and the organic matter fractions in soil samples were measured alongside obtaining three classes of stable aggregates by dry sieving: macroaggregates (4.75 -2.00 mm), mesoaggregates (2.00 -0.25 mm), and microaggregates (0.25 -0.053 mm) (Costa Junior et al., 2011, 2012;Cavalcante et al., 2019).The total organic carbon (TOC) and total nitrogen (TN) of the soil were determined, respectively, by wet oxidation methods with external heating (Yeomans and Bremner, 1988) and sulfuric digestion followed by distillation (Kjeldahl) (Tedesco, 1985).Labile organic C (LOC) was determined according to the procedure proposed by Blair et al. (1995), which was adapted for tropical soils by Shang and Tiessen, (1997).The determination of labile N (LN) was performed by adapting the Sahrawat (1982) method.
In the physical fractionation of SOM (Sohi et al., 2001), the free light organic matter (LOM L ) fraction was obtained with the densimetric method; the coarse (POM C ) and fine (POM F ) particulate organic matter and mineral associated organic matter (MOM) fractions were determined with granulometric fractionation.The LOM L separation was performed OF 50 Mg ha −1 (dry basis) of poultry manure (1/2 dose applied to the planting furrow and 1/2 cast between the rows).
CF 8 Mg ha −1 of dolomitic limestone with 80 % PRNT, with 1/3 of the dose applied to the planting furrow and the remaining 2/3 between the rows.Additionally, 1.5 Mg ha −1 of Bayovar natural reactive phosphate was applied to the bottom of the furrow and 0.70 ha -1 between the rows.

OF+CF
The final treatment combined the aforementioned OF+CF.After planting the coffee, OF+CF and CF fertilization schemes were conducted, and from the second year, coffee production fertilizers were applied, both as blanket fertilizers in the planting rows using NPK fertilizers on 20-05-20 and in similar doses (600 g plant −1 , divided into two applications during the growing season) among the treatments.

NP
In which pre-emergent herbicide was applied when necessary.

NV
The adjacent area, covered by a fragment of NV defined as a Seasonal Semi-deciduous Montane Forest (IEF, 2016), was considered in this study as a reference area, although it was not previously part of the described design.
Rev Bras Cienc Solo 2023;47:e0220126 through the centrifugation of 5 g of fine air-dried soil in triplicate, preceded by manual agitation with NaI (d = 1.60 g cm -3 ); the supernatant with LOM L was slowly filtered with the aid of a vacuum pump, thus separating the LOM L from the NaI solution.The C and N contents of the LOM L fraction were determined by dry combustion in a CHNS elemental analyzer (Perkin Elmer).The remaining soil material (heavy fraction -HF) was sifted using 0.210-and 0.053-mm sieves, yielding coarse POM (POM C , 2.00 -0.210 mm) and fine POM (POM F , 0.210 -0.053 mm) fractions.After drying, the organic C and total N contents of these fractions were determined according to Yeomans and Bremner (1988) and Tedesco (1985), respectively.The C and N contents of the mineral associated organic matter (MOM) were calculated by the difference between the total soil C and N and the organic fractions: LOM L +POM C +POM F .

Soil quality index (SQI)
The SQI was developed based on principal component analysis (PCA), to evaluate the effect of different management systems on the rehabilitation of the bauxite-mined area.
Methodologies used by Karlen and Stott (1994) and Brejda et al. (2000) were adapted for this purpose, following the procedure proposed by Cavalcante et al. (2021).
The following steps were performed: a) 37 properties (original variables that play important roles in the soil) of the management systems and reference area were initially selected; b) a normality hypothesis test was conducted using the Lilliefors test (p<0.10);c) the data that presented a normal distribution were submitted to a multicollinearity test with the traditional trail analysis method (Cruz and Carneiro, 2006), with the exclusion of highly correlated variables, that is, to verify the presence or absence of dependencies between the variables.The multicollinearity diagnostic was made using the variance inflation factor (VIF) method using statistical software, in which VIF values greater than 10 were considered (Hair et al., 2009); d) PCA was then performed as a method of factor extraction, to only interpret the factors that presented eigenvalues >1.0.The retained factors were submitted to the varimax rotation to maximize the relation between independent attributes.In each principal component, only variables with high factor loadings were retained, i.e., those with absolute values within 10 % of the highest factor loading (Andrews et al., 2002), so that in the end, properties representative of the physical, chemical, and organic soil quality were selected; h) the selected properties were submitted a new varimax rotation and the respective rotated factor loadings were used in the calculation of the relative weights of the attributes on SQI, according to the following equation: in which: Wi is the relative weight of the i-th variable in SQI; F 1i and F 2i are the eigenvalues of the principal components; Pi is the rotational factorial loading of the i-th variable; Pj is the rotational factorial loading of the j-th variable; i and j are indices for the variables; and n is the number of variables involved in the PCA.After the selection of the indicators, the values of each variable were normalized by means of a relative standardization to be included in the SQI, being transformed into scores of the indicators (S), from 0 to 1, according to Liebig et al. (2001) and Bhardwaj et al. (2011).The SQI was determined as: in which: SQI is a number from 0 to 1; W i is the PC weighting factor, a number between 0 to 1; and S i is the indicator score for variable i, a number between 0 to 1.
Higher scores in the model indicated a better soil quality or greater performance of soil function.

Statistical analysis
Two-way ANOVA was used to check the effects of fertilization, ground covers, and their interactions on organic fractions in different aggregate-size classes of soils in bauxite mining areas undergoing rehabilitation.Tukey's HSD post hoc test at p<0.05 was used to evaluate the differences between the means.The PCA loaded variables were subject to Pearson correlation analysis at p<0.01, p<0.05, and p<0.10.The data were analyzed using the Microsoft Excel XLSTAT software (version 2021.3.1)(Addinsoft, 2021).

Effect of management systems on C and N pools
The interaction between fertilization and cover crops was not observed in the soil (Table 2).Significant effect of fertilization was observed in the TOC, TN, C-MOM, N-LOM L , N-POM C , N-MOM contents, and C-LOM L , C-POM C, C-MOM, N-LOM L , N-POM C contents in the cover crops.The OF+CF fertilization contributed to the increase of these contents in soil, differing from NF, although having been equivalent to OF and CF (Table 3).However, the OF and CF fertilization treatments contributed to the higher N contents in LOM L .Brachiaria grass (B) increased C and N contents in the soil and soil fractions regarding the subplots without plants (Table 3).While the use of Stylosanthes (S), despite contributing to improve of C and N contents, did not statistically differ from areas without plants (NP).

Selection of minimum data set
Among the 37 selected variables for the selection of the minimum data set (Table 4), 18 did not present deviations following the normal distribution (Organic: C-LOM L , C-POM C , LN, N-LOM L , N-POM C, N-POM F , and N-MOM; Chemical: P, K, Mg 2+ , Al 3+ , H+Al, and m; Physical: Sand C , Sand F , Silt, Clay and TP), these were eliminated from the process.
According to the trail analysis method performed through the multicollinearity test, the diagnosis for the explanatory variables revealed that of the 19 soil parameters, 11 did not present multicollinearity (Organic: LOC, C-POM F , and C-MOM; Chemical: pH, Ca 2+ , t, and P-rem; Physical: BD, Micro, AGRI W , and SAI W ) (Table 5).A high frequency of correlations indicates that organic, chemical, and physical soil properties can be organized into factors based on their correlation structures.However, not all parameters correlated with each other, namely in the case of C-POM F and BD (Table 6).In general, the C-MOM, pH, Ca 2+ , t, P-rem, Micro, and SAI W were positively correlated with most soil properties, while in contrast, AGRI W was negatively correlated with C-POM F .
The three first factors had eigenvalues exceeding one, and 89.61 % of the total variance were explained by three principal components (PCs), being retained for interpretation (Table 7).The factors explained more than 90 % of the variation in C-MOM, t, and SAI W ,   7).These properties are relevant to both agronomic and environmental studies, since they provide information about soil processes or its behavior (e.g., LOC, C-MOM, and pH).Other properties are equally important, namely those that indicate the capacity of the soil to resist cation exchange (e.g., t) and those that establish fundamental relationships with hydrological processes, as well as providing an essential function in the supply and storage of water, nutrients, and oxygen in the soil (e.g., BD).
During the PCA analysis, the indicators were selected and the weights (W i ) of each variable were calculated as a function of the eigenvalues and the explicability of the  indicator by the retained factor (Table 8).The weights of the selected properties ranged from 0.07 to 0.21.The soil properties that had higher W i values in SQI were those linked to organic and chemical indicators and those properties characterized by lower W i values were linked to physical indicators.

Effect of management systems on soil organic properties
The use of inputs, whether of mineral or organic nature, can favor the increase of organic properties through different management practices, thus representing an important tool for the acquisition and maintenance of soil quality (Loss et al., 2009(Loss et al., , 2011;;Valadão et al., 2011;Cavalcante et al., 2019).
An increase in SOM due to herbaceous plants, principally those of the genus Brachiaria, has been observed in several papers (Pillon et al., 2011;Wendling et al., 2011;Carmo et al., 2012;Rossi et al., 2012;Cavalcante et al., 2019).The large number of roots incorporated into the soil by fast renovation and incorporating this grass, reflect the SOM contents and soil aggregation.It has been shown that stylosanthes, which present annual and bi-annual features and depend on restocking from natural sowing to remain in soils, presented a smaller contribution of C and N to soils through their root systems, thus limiting the recovery capacity of the soils.Thus, the effect of organic and chemical fertilization treatments on these variables is more strongly related to the development of herbaceous species, such as those that acquire their nutrient uptake through rich root systems and N 2 fixation (Silva and Corrêa, 2010;Barbosa et al., 2019), helping increase or maintain these properties in soils.
In general, the soil organic properties in areas with management systems were still lower than those observed in reference areas, that is, below the ideal condition of the reference considered in this research.However, the results showed that the properties essential for the proper functioning of the soil could be recovered through the adopted management practices.

Organic fraction distribution and its relationship with aggregate size
There was no significant interaction (p>0.05) between the fertilization treatments and the ground covers and the C, N, and SOM fractions in the evaluated aggregate classes, allowing for discussion of only the main effects of the treatments (Table 2).At 0.00-0.20 m soil layer, the rehabilitation alternatives influenced the distribution of organic compounds in the aggregate classes, showing that the restoration of vegetation promoted an increase in the organic contents and stability of the soil aggregates (Jiao et al., 2012;Raiesi, 2012;Deng and Shangguan, 2017).However, regardless of the treatment, the fractions obtained in the three structural classes were not equivalent or superior to those of NV, which had the highest values (Table 3).The greater contribution of plant waste and the fact that bauxite mining activities did not change the original condition of soil aggregation may be related to the results obtained in this situation.Changes in land management compromise soil quality (Deng and Shangguan, 2017).The results of management alternatives on the distribution of TOC and TN in the organic soil fractions were unique in each aggregate class (Table 2).
In situations where there were significant differences (p<0.05), it was confirmed that the fertilization treatments alternated with regard to those that provided the highest C or N values in the studied aggregate classes.As for ground cover, in general, the Brachiaria (B) effectively (p<0.05)contributed to an increase in C and N contents in the fractions for which significant differences were observed in macro-, meso-, or microaggregates (Table 3).Stylosanthes (S) practically did not alter the fractions of organic matter in which they had significant differences (p<0.05),like the absence of plants.
Free light fraction (LOM L ) comprises the organic compartments most sensitive to soil degradation by farming (Freixo et al., 2002) and can be considered a premature indicator of the decline of SOM contents (Six et al., 2002;Wu et al., 2004;Lima et al., 2008;Pulrolnik et al., 2009;Santos et al., 2013).This fraction contains root fragments, seeds, leaves, and twigs in the least advanced stages of decomposition, and it presents a significant seasonal and spatial variability since it is directly linked to the soil residue supply (Christensen, 2000).Among all the fractions, the LOM L presented the lowest proportion of C and N (Figure 3).The highest C and N contents in the LOM L fraction were found in the mesoaggregates, especially those that had received treatments consisting of fertilization and Brachiaria (OF/B, CF/B, and OF+CF/B).Loss et al. (2011) also observed a similar result.The incorporation of plant residues, either as a source of N or combined with the application of nitrogen fertilizers, contributed to higher percentages of C in mesoaggregates (Huang et al., 2015).In macroaggregates, there was a greater effect of N in all fractions.Better results of TN on macroaggregation have been observed (Hontoria et al., 2016).Among all the cases, the effect of CF fertilization was observed only in soil macroaggregates.
The OF led to higher values of C-POM C in the mesoaggregates and C-POM C , C-POM F , and N-POM F in the microaggregates, however, C-POM C and N-POM F in the meso-and microaggregates, respectively, which OF was like NF.These results confirm organic fertilization as an important management strategy to improve soil quality and highlight the importance of particulate fractions in evaluating soil quality.According to Rossi et al. (2012), the storage of C in particulate and light fractions is closely related to the recent input of vegetable residues, which supports the action of Brachiaria in maintaining these contents.According to the same authors, the particulate fraction may effectively indicate disturbances to land-use.Higher proportions of C and N of POM C and POM F were here found in the macro-and microaggregates, respectively (Figure 3).Huang et al. (2010) observed that the use of organic fertilization alone or in combination with inorganic fertilizer favored higher stocks of C-POM in soil microaggregates.The addition of nitrogen fertilizers may contribute to higher P contents in microaggregates (Wang et al., 2016).The difference in the proportions of free and particulate fractions observed in this study may have been related to mechanisms of the chemical stabilization of the SOM or to the degree of humification of the SOM in these fractions (Figure 2).Campos (2003) also revealed that the bigger the specific surface of the particle, the bigger its capacity to interact with organic matter.Accumulation of C in the particulate fraction and in the Rev Bras Cienc Solo 2023;47:e0220126 free fraction is closely associated with the recent supply of plant material (Rossi et al., 2012), which corroborates the action of Brachiaria in maintaining the soil C contents (Table 3).According to this same author, the particulate fraction may effectively indicate the differences between the use systems applied.
Among all the fractions, the highest C and N contents were found in MOM, in all aggregates classes, and with a contribution from all management systems (Figure 3).In the MOM fractions, the macro-, meso-, and microaggregates exhibited a similar behavior regarding the effect of fertilization, with higher values of N occurring in both classes of aggregates when a OF+CF fertilization procedure was used.This effect also was observed in C-MOM in the mesoaggregates.The dynamics of the MOM fraction would be closely related to the texture (Feller and Beare, 1997;Freixo et al., 2002) and the surface of the mineral particles (Campos, 2003).The largest proportion of C and N in this fraction could also be associated with the supply of organic matter, mainly by Brachiaria (Table 3), as suggested by Carmo et al. (2012).Santos et al. (2011) observed that the largest C contents in a native field area were associated with minerals, relating this effect to chemical protection mechanisms of organic matter with clay contents.Indeed, soils with low clay contents generally have smaller protection of organic matter and low capacity of the mineral fraction to keep a relative stock of C in MOM, leading to the vulnerability of the management system (Santos et al., 2013).However, the higher proportion of C and N in MOM indicates that the complexation mechanism had a higher relevance than the occlusion mechanism.Thus, MOM is the fraction most stable to management changes, being of great importance for the storage and availability of C and N in the soil.

Organic, chemical, and physical soil quality indicators
The SQI was significantly higher (p<0.05) for the NV (1.0 ± 0.03) than the other management systems (Figure 4).Higher SQI scores for preserved sites were previously found by Bhardwaj et al. (2011) and Chaves et al. (2017), who noted the alteration of soil properties when replaced by another management system with soil use intensity.
The mean contribution of the variables to the SQI (0.00-0.20 m soil layer) followed this order: t (0.23 ± 0.01) > pH (0.20 ± 0.01) > SAI W (0.08 ± 0.002) > C-MOM (0.07 ± 0.01) > LOC (0.06 ± 0.01) = BD (0.06 ± 0.01).The chemical variables showed a greater contribution to the SQI (Figure 4).The cation exchange capacity (CEC) is an indicator of soil fertility (Brady and Weil, 2007).The pH showed a significant positive correlation with t (r = 0.820; p<0.01) (Table 6), and it had a direct relationship with the availability of nutrients (Bohn et al., 2001).Mukhopadhyay et al. (2016) found a significant positive correlation between pH and CEC (r = 0.937; p<0.01) in the elaboration of the SQI for the assessment of recovered coal mine residues.The OF+CF/B management system contributed the most to the SQI values of t (0.27) and pH (0.22).OF/B and NF/NP treatments had significantly higher SAI W values than all other similar rehabilitation systems studied (Cavalcante et al., 2019).Mbagwu and Schwertmann (2006) also showed that soils with high levels of Al oxides than Fe oxides favored soil Rev Bras Cienc Solo 2023;47:e0220126 macroaggregation.However, the alteration of management systems may reduce SOM levels and, consequently, the stability of aggregates.
In general, the proportions of C and N in different compartments of organic matter (LOM L , POM C , POM F , and MOM) were smaller in the most labile fractions (LOM L : C = 0.1 %, N = 0.1 %; POM C : C = 1.5 %, N = 1.3 %; and POM F : C = 8.2 %, N = 2.2 %), reflecting on the dynamics and availability of these nutrients in the soil (t x C-MOM: r = 0.877; p<0.01); comparatively the highest levels in MOM (C = 90.2%; N = 96.4%), whose dynamics would be closely related to texture (Feller and Beare, 1997;Freixo et al., 2002).The interaction between the organo-mineral complexes would also be directly related to the surface of the mineral particles, as mentioned earlier.Santos et al. (2011) noted higher C contents in soils associated with minerals, relating this effect to the mechanisms of chemical protection of organic matter with the levels of clay.
According to the results obtained by Cavalcante et al. (2019), there is a great positive correlation between the carbon management index (CMI) and LOC (r = 0.998; p<0.01).Thus, short-term LOC values can provide a sensitive indicator for changes in soil C dynamics through the adoption of management practices for the rehabilitation of mined areas.Bulk density has also been used as an indicator of soil quality, as this is a dynamic property that is susceptible to use and easy to determine, since it is related to compaction and the relative restriction to root growth (Arshad et al., 2002).

CONCLUSIONS
When organic (OF), chemical (CF), and OF+CF fertilizers were applied to areas cultivated with coffee intercropped with Brachiaria (B), they increased C and N contents in the soils and aggregates classes, as well as in compartments of SOM.The minimum set of soil  quality indicators of reclaimed bauxite-mined areas could consist of organic indicators: labile organic carbon (LOC) and mineral-associated organic matter (C-MOM); chemical indicators: pH and effective cation-exchange capacity (t); and physical indicators: bulk density (BD) and stable aggregates index in water (SAI W ). The t and pH were the variables most sensitive to the management systems used in the rehabilitation of the mined area, and were thus considered the best indicators of soil quality.Brachiaria was the cover plant that contributed most to increasing the soil quality of mined bauxite areas by increasing the SQI, especially when fertilized.Thus, the use of fertilizers and ground cover plants intercropped with coffee in the bauxite-mined areas improved the soil's physical, chemical, and organic properties, providing a viable option for reconditioning mineral exploration areas.

Figure 2 .
Figure 2. Area preliminarily cultivated with coffee (last ten years) (a); reconfiguration of the mined area and replacement of the A horizon soil stored prior to mining (b); experimental area comprising zones of coffee cultivation, different types of fertilization treatments, and intercropping cover crops species (c); and experimental layout with ground cover treatments implanted between the crop rows (d) (Cavalcante et al., 2019).

Figure 3 .
Figure3.Carbon and nitrogen ratios of free light organic matter (LOM L ), coarse particulate organic matter (POM C ), fine particulate organic matter (POM F ), and mineral-associated organic matter (MOM) in soils and in aggregate size classes of soils from natural forest and a rehabilitated bauxite mining area subjected to different fertilization treatments and with different cover crops.

Figure 4 .
Figure 4. Soil quality index (SQI) means and contributions of each indicator in the SQI at 0.10-0.20 m soil layer, in areas covered with native forest and after three years of recovery with coffee using different fertilization treatments and ground cover plants.NF: no fertilization; OF: organic; CF: chemical; OF+CF: organic + chemical; NP: no plants; B: Brachiaria grass; S: stylosanthes; NV: native vegetation; LOC: labile organic carbon; MOM: mineral-associated organic matter; t: effective cation-exchange capacity; BD: bulk density; SAI W : stable aggregates index in water.

Table 1 .
Details of the treatments used in the rehabilitated area after bauxite mining with coffee cultivation and intercropping cover crops subjected to different fertilization procedures

Table 2 .
ANOVA results for the effects of fertilization, cover crops and their interaction on C and N fractions in different aggregate size classes of soils in bauxite mining area subjected to rehabilitation

Table 3 .
Carbon and N fractions in soil and in different aggregate size classes of soil under natural forest and in a rehabilitated bauxite mining areas with fertilization and cover crops NF: no fertilization; OF: organic; CF: chemical; OF+CF: organic + chemical; NP: no plants; B: Brachiaria grass; S: stylosanthes; TOC: total organic carbon; TN: total nitrogen; LOM L : free light organic matter; POM C : coarse particulate organic matter; POM F : fine particulate organic matter; MOM: mineralassociated organic matter.The means are followed by the same letter(s) within each line, which do not significantly differ from p<0.05 by Tukey's test.

Table 4 .
Means of the physical, chemical, and organic properties of the soil (0.00-0.20 m) in native forests and in a rehabilitated area after bauxite mining cultivated with coffee and intercropping cover crops subjected to different fertilization procedures * The variable from which the sample was extracted does not follow the normal distribution according to the Lilliefors test at <0.05; ns: not significant; NV: native vegetation; NP: no plants; B: Brachiaria; S: stylosanthes; TOC: total organic carbon; LOC: labile organic carbon; TN: total nitrogen; LN: labile nitrogen; LOM L : free light organic matter; POM C : coarse particulate organic matter; POM F : fine particulate organic matter; MOM: mineral-associated organic matter; BS: base sum; CEC: cation-exchange capacity; t: effective cation-exchange capacity; V: base saturation; m: saturation by aluminium; P-rem: P remaining; Sand C : coarse sand; Sand F : fine sand; BD: bulk density; Macro: macroporosity; Micro: microporosity; TP: total porosity; MWD W : mean weight diameter obtained by wet sieving; GMD W : geometric mean diameter obtained by wet sieving; AGRI W : percentage of water stable aggregates >2 mm; SAI W : stable aggregates index in water.

Table 5 .
Multicollinearity diagnostic for the explanatory variables, based on the variance inflation factor (VIF), of management systems evaluated at 0.00-0.20 m soil layer Bold numbers indicate variables that presented VIF >10.TOC: total organic carbon; LOC: labile organic carbon; TN: total nitrogen; POM F : fine particulate organic matter; MOM: mineral-associated organic matter; BS: base sum; CEC: cation-exchange capacity; t: effective cation-exchange capacity; V: base saturation; P-rem: remaining; BD: bulk density; Macro: macroporosity; Micro: microporosity; MWD W : mean weight diameter obtained by wet sieving; GMD W : geometric mean diameter obtained by wet sieving; AGRI W : percentage of water stable aggregates >2 mm; SAI W : stable aggregates index in water.

Table 6 .
Pearson correlation coefficients between soil properties, evaluated for 0.00-0.20 m soil layer, in areas of natural forest and bauxite mining undergoing rehabilitation procedures and being cultivated with coffee and intercropping cover crop with different treatments of fertilization

Table 7 .
Rotated factor loadings and communalities of a three-factor model of organic, chemical, and physical soil properties in a bauxite mining area under rehabilitation and cultivated with coffee and intercropping cover crops with fertilization treatments LOC: labile organic carbon; POM F : fine particulate organic matter; MOM: mineral-associated organic matter; t: effective cation-exchange capacity; P-rem: P remaining; BD: bulk density; Micro: microporosity; AGRI W : percentage of water stable aggregates >2 mm; SAI W : stable aggregates index in water.The bold factors are considered highly loaded and included in the index.

Table 8 .
Organic and chemical properties selected as soil quality indicators of soils from 0.00-0.20 m soil layer, their rotational factor loads and weights in the soil quality index (SQI) of areas covered with native forests and after three years of recovery with coffee using different fertilization treatments and cover crops LOC: labile organic carbon; MOM: mineral-associated organic matter; t: effective cation-exchange capacity; BD: bulk density; SAI W : stable aggregates index in water.Bold factors are considered highly loaded and included in the index.Rev Bras Cienc Solo 2023;47:e0220126