Analyzing the Impact of Intensive Agriculture on Soil Quality: A Systematic Review and Global Meta-Analysis of Quality Indexes

: As the global population grows, the demand for food increases and puts a strain on food production systems and agricultural productivity, causing soil degradation. Soil quality indexes (SQIs) have been developed to maintain and improve soil quality. However, due to the variety of soils and SQIs, analyzing and comparing results has historically been difﬁcult. Therefore, in this study, we carried out a systematic review with meta-analysis focused on soil quality studies of agricultural soils under intensive agriculture using the uniﬁed weighted additive SQI methodology (SQI U ). Sixty-ﬁve observations of agricultural soil quality were analyzed. Chemical indicators were the most prevalent in the SQI U s, followed by physical and biological indicators. Conventional soil management had negative effects on soil quality (–7.55%). From the factors analyzed, the minimum database had a signiﬁcant effect on the soil quality results, but not the number of indicators that made up the SQI U . The SQI U made up of chemical–biological indicators (CB) presented negatively overestimated measurements of soil quality (–32.53%), exaggerating the degradation of the analyzed soils. The indicators that correlated most strongly with the size of the effect on agricultural soil quality were the cation exchange capacity (CEC) and microbial biomass carbon (MBC). The SQI U is a feasible tool to interpret the quality of agricultural soils around the world, as it makes it possible to obtain a simple and generalized view of soil conditions.


Introduction
The growing population and a shift towards high-consumption markets have increased pressure on the food production chain, resulting in a need for larger areas of cultivation to meet the population's nutritional needs [1]. This change in land use is a concern for government authorities who implement measures to improve and maintain agricultural productivity through the technification of conventional agricultural work in intensive agriculture. However, these measures have resulted in the emergence or escalation of soil-degrading processes [2,3]. Some intensive agricultural activities implement inadequate soil management practices, such as the cultivation of monocultures and the excessive and inadequate use of chemical fertilizers and pesticides, reducing soil biodiversity and provoking leaching processes, the eutrophication of aquifers due to runoff phenomena, and knowledge of trained personnel for the selection and degree of importance of the indicators to be included [11,12]; additive indexes (SQI A ), which analyze the indicators, considering them to have linear behavior [11,12]; weighted additive indexes (SQI W ), which use weights previously established by experts in the indicators to be included [11,12]; Nemoro indexes (SQI N ), used in soil quality studies around the world [11,12]; and unified weighted additive indexes (SQI U ), which establish the weights of the indicators based on statistical techniques. SQI W and SQI U are the most used in agricultural soils, and their application depends largely on those responsible for monitoring [11,12]. However, the SQI W methodology can result in the selection of indicators that do not correlate with the problem or phenomenon to be studied, since it is based on indicators previously established. In contrast, the SQI U methodology is not influenced by human subjectivity, and thus allows for the mathematical selection of those indicators that are most closely linked to the problem or phenomenon of interest. This guarantees a more precise selection of relevant indicators. For this reason, it was decided that the SQI U methodology would be used in this study.
To date, the only study focused on analyzing the methodologies and characteristics of the different SQIs developed on a global scale is a systematic review carried out by Sarmiento et al. (2018) [12]. They analyzed 70 research works published from 2001 to 2017, and identified 32 SQIs developed using various methodologies and study areas. Their analysis concluded that reliable indexes should include the following indicators: microbial biomass carbon (MBC), nitrogen (N) mineralization potential, soil pH, electrical conductivity (EC), CEC, porosity, BD, and hydraulic conductivity (HC). As previously stated, due to various factors involved in soil formation-e.g., parent material, topography, climate, and agricultural practices in the region-there is no consensus on a generalized SQI, or the number and type of indicators it should include. This heterogeneity makes it challenging to compare SQIs developed globally and establish standardized measurements or observations of soil quality. To address this issue, a systematic review with meta-analysis was conducted on a global scale, specifically focusing on studies that used the SQI U methodology for soil quality monitoring. The objectives of this study were to (i) determine the global impact of intensive agricultural practices on the quality of soils, and to (ii) investigate any potential effects on soil quality observations depending on factors such as the database used, and the number and type of indicators included in the SQI U . It was hypothesized that intensive agricultural practices have a negative impact on soil quality, and that there is a significant effect on quality observations based on the type of database used and the number and type of indicators included in the different SQI U s developed.

Systematic Review
A literature search was conducted to find experimental studies that developed SQIs in intensive agricultural soils employing the SQI U methodology. The databases used are presented in Table 1. Google Analytics tools were used to identify the most frequently used keywords related to the topic of soil quality in agricultural soils, finally selecting "soil quality index", "agricultural soil", "principal component analysis", and "unified weighted additive index". With the keywords and synonyms, the search equation to be used in the databases from Table 1 was established as follows: (soil quality index OR SQI) AND (agricultural soil) AND (PCA OR principal components analysis OR unified weighted additive index). Boolean operators were also used to obtain more accurate results in the search process.
The articles were selected through a three-step filtering process, which was independently carried out by three members of the research team. For the first filter, every article had to meet several criteria, including containing at least one of the selected keywords in the title or abstract, being a research article (not a review), and having a DOI. The second filter required the articles to be published in English between 2011 and 2021. After applying the second filter, the articles were downloaded. If an article could not be downloaded, the corresponding author was contacted for assistance. Finally, for the third filter, articles had to meet specific requirements, such as including quantitative (numerical) soil indicators (physical, chemical, or biological), providing statistical information (including the standard deviation and mean), using experimental controls during the research, and presenting Agronomy 2023, 13, 2166 4 of 21 sufficient data to calculate these parameters, using Equation (1). By applying these three filters, we were able to select the most relevant and appropriate articles for our study [13].
where x i is the mean of the control group or treatments used, and SD is the standard deviation of the control group or treatments. Articles that did not meet the criteria during the filtering process were excluded. At the end of the filtering process, the members of the research team who participated in the selection compared their databases, and those articles that were present in the three databases were selected, thereby forming a single database.

Data Collection
In this study, the target variables (the variables on which the meta-analysis focused) were the soil quality observations obtained through the SQI U developed in each study. Additionally, all the studied indicators related to the condition of the soil (physical, chemical, and biological indicators) were selected as secondary variables. To gather data, mean and standard deviation information for both the target and secondary variables was collected from the tables of results and figures, or calculated from the articles. The information extracted from the figures was obtained using the software Engauge Digitizer 12.2.1 (https: //markummitchell.github.io/engauge-digitizer/) (accessed on 15 August 2023) [14]. Data collection and database building were carried out independently by three members of the research team. When creating the databases, the research team only considered scientific articles that focused on "agricultural" land uses. In the case of articles that involved soil samples taken at different depths, only those from "0 to 30 cm depth" were considered. Similarly, if the articles compared different methods for generating SQIs, only observations that utilized the "SQI U " methodology were considered. In addition, the team converted the units of quantitative variables to the most used units in publications related to the topic of analysis. At the end of the process, these databases were compared to create a final database.

Statistical Analysis
Analyses were performed using R statistical software version 4.0.5 [15]. First, we performed a descriptive analysis of the selected studies to determine the year of publication, geographic location, and frequency of indicators in the various SQI U s developed. For the systematic review, a publication bias test was performed with an Egger's symmetry test featuring a significance level of p ≤ 0.05. At the same time, the fail-safe number (N fs ) technique was applied, with a significance level of p ≤ 0.05 [16]. The N fs measures the robustness of meta-analysis results by determining the number of omitted observations required to refute the conclusions of the analysis. A meta-analysis was then performed using the statistical packages Metafor [17], Metaviz [18], Puniform [19], and Robumeta [20]. The natural logarithm of the response ratio (lnRR)-defined as the "effect size"-was used to measure the effect caused by intensive agricultural activity on soil quality, obtained through the various SQI U s developed. We also analyzed the factors of the type of database used (total database and minimum database), the number of indicators that made up the SQI U s (less than or equal to five and more than five), as well as the type of indicators (chemical, physical, and biological). The lnRR of each observation was calculated using Equation (2) [16,21].
where lnRR is the natural logarithm of the response ratio or effect size, and x t and x c are the mean values of the treatments and control group, respectively. The variance of the effect size (ν) was calculated using Equation (3) [16,22].
where n t and n c represent the sample size of the treatment and control groups, respectively; s t and s c represent the standard deviations SD = SE √ n of the treatment and control groups, respectively; and SE represents the standard error. The inverse of ν was used to obtain the weight factor (w) using Equation (4) [21].
The w was used to calculate the average effect size lnRR for all observations according to Equations (5) and (6) [22]: Using Equation (7), we obtained the percentage values of the lnRR , which were used to express the effects of the different observations [21].
In addition, a "random effects" model was used to determine whether the values of the lnRR were significant. To determine the 95% confidence intervals (CI), a bootstrapping test was performed with 999 iterations. If the CI of the values of the lnRR did not cross the zero-effect line, the effects were considered significant (p ≤ 0.05).

Metaregression
The relationship between the various indicators that made up the developed SQI U s was analyzed, selecting those indicators whose frequency of integration was five or more. A minimum number of at least six observations per indicator in the various studies was taken into consideration [13].

Metamodeling
To select models that best explain the variability of the values of the lnRR of the various observations, we used combinations of two indicators closely correlated with the lnRR , without considering interactions between them. To determine which of the developed models was the best to predict the lnRR , the corrected Akaike information criterion (AICc) and the Bayesian information criterion (BIC), as well as the marginal coefficient of determination R 2 , were used at a significance level of p ≤ 0.05 as a metric of the a measure of a model's relative quality compared to other models [13].

A Systematic Review and Descriptive Analysis
Of the 5142 articles analyzed, 22 were selected and added to the final database (Table S1). Figure 1 illustrates the articles compiled from several searches performed in different databases (Table 1)  from the Asian continent, with China being the country with the highest scientific output in the development of SQI U s ( Figure 3). This fact could be attributed to China's policies on food security for its growing population, which present a need for scientific knowledge to maintain and improve the quality of agricultural soils in certain areas and traditional crops [23]. The present study found that in Latin America-specifically in Mexico, Argentina, and Brazil-there was only one research study focused on the development of SQI U s [24][25][26], which suggests that there is a lack of technical knowledge and expertise in this field in these countries. During the selection process of articles for the final database, it was found that most articles were not included due to the absence of important statistical information-such as SD or SE-and the lack of controls for the treatments. This highlights the importance of using appropriate experimental designs and open science practices-e.g., creating open databases-to ensure the proper presentation and accessibility of research results. ronomy 2023, 13, x FOR PEER REVIEW 6 of 2 attributed to China's policies on food security for its growing population, which presen a need for scientific knowledge to maintain and improve the quality of agricultural soil in certain areas and traditional crops [23]. The present study found that in Latin Amer ica-specifically in Mexico, Argentina, and Brazil-there was only one research study fo cused on the development of SQIUs [24][25][26], which suggests that there is a lack of technica knowledge and expertise in this field in these countries. During the selection process o articles for the final database, it was found that most articles were not included due to th absence of important statistical information-such as or -and the lack of control for the treatments. This highlights the importance of using appropriate experimental de signs and open science practices-e.g., creating open databases-to ensure the prope presentation and accessibility of research results.  Of the indicators that made up the various SQI U s developed in the selected studies, chemical indicators were the most used, followed by physical indicators, and finally biological indicators (Figure 4). This finding is consistent with Sarmiento et al. (2018) [12] who through a systematic review reached the same conclusions on the use of the types of indicators in the SQIs developed worldwide. The indicators that most frequently composed the SQI U s were soil pH and soil organic C (SOC), both belonging to the type of chemical indicators.  Of the indicators that made up the various SQIUs developed in the se chemical indicators were the most used, followed by physical indicators, a    Researchers have often used chemical and physical indicators to assess soil conditions, as some of the physical indicators are simple to use and do not require extensive training or sophisticated equipment. These indicators can be easily applied in the field [24,27]. Physical indicators provide information on the structure and functioning of the soil, but with less variability over time and, consequently, less sensitivity to changes in soil condition, as changes become apparent after a long period of time. As for chemical indicators, the interpretation provided by their analysis is complex, since their application and interpretation require preparation and knowledge of technical and laboratory procedures, making their implementation in the field more difficult. On the other hand, chemical indicators are more sensitive, providing more detailed information on nutritional conditions and biogeochemical cycles of the soil, and showing possible soil-buffering conditions of a lesser magnitude than physical indicators [7]. Regarding biological indicators, these have historically been the least utilized, despite being the most sensitive indicators, as they are able to detect small changes in soil conditions. This low use could be attributed to the lack of recognition of the role of microbial processes in soil organic matter cycling, difficulties with measurement, and uncertainties associated with linking laboratory-scale measurements to field processes. Among the biological indicators, the MBC is the most widely used [12,28], probably because the implementation of other biological indicators requires solid knowledge of biology, chemistry, and the development of microorganisms, together with the requirement of a longer period to carry out their analysis and interpretation. However, the use of biological indicators has become more common in the last decade, presenting itself as an option for analysis with greater accuracy in the detection of specific or smaller-scale cases, and being promoted by more researchers in the field of soil science [7,9,12].

Meta-Analysis
One of the key tests used to verify the results of a meta-analysis is the publication bias test. This test ensures that no observations from relevant studies have been omitted, which could affect the outcome of the meta-analysis. The results of the publication bias test are presented in Figure 5, using Egger's symmetry test, where the symmetry of the points in the graph is measured [16]. The test showed a significant result (p ≤ 0.05), indicating that no observations were omitted during the systematic review process. In addition, the N fs calculation was performed, which resulted in a value of 47.04 × 10 6 quality observations, with a significance level of p ≤ 0.0001. This confirms that the conclusions obtained from the meta-analysis are significant (p ≤ 0.05), and suggests that more than 47 million quality observations must have been omitted to obtain inconclusive results.

Effect on the Quality of Intensive Agricultural Land
The lnRR for the quality of intensive agricultural soils measured by the SQI U methodology is shown in Figure 6. The lnRR obtained was significant [MD of −7.55%, 95% CI (−14.79 to 0.30%), p ≤ 0.05, with a number of observations k = 65], implying an overall loss of soil quality of 7.55% ( Figure 6). The negative effect of intensive agricultural activity has been observed and reported by many studies. Soil quality deterioration can be caused by physical factors such as degradation, water and wind erosion, compaction, and reduced water availability, as well as chemical factors such as loss of fertility, salinization, eutrophication, and pollution. These factors are often the result of intensive agriculture and poor soil management practices [6,8]. The excessive use of chemical fertilizers in intensive agriculture has caused soil acidification and salinization, which results in an increased leaching of N compounds into nearby aquifers. This continuous deterioration of the soil has reached a point where the soil is abandoned [29,30]. Inappropriate land use refers to both the "non-rotation" or "non-cyclicality" of crops and the application of inappropriate or excessive regimes of pesticides and chemical fertilizers, overgrazing, deforestation, and land use change [31]. In order to address soil deterioration, techniques such as zero or conservation tillage have been developed, alongside the use of organic amendments such as manure, biosolids, and composts. These techniques aim to increase soil OM and improve soil physicochemical and biological conditions [11,13,32,33], which allows for the recovery of the soil structure and its proper functioning, in order to carry out basic ecological functions such as nutrient cycling, water supply, C-CO 2 fixation, and OM mineralization [34]. To understand the impact of soil management techniques, it is important to have methods that accurately assess the initial conditions of the soil. The information obtained can then be used by soil managers and authorities to make informed decisions. The implementation of SQIs through research studies is crucial in developing these methods and obtaining reliable data [9,35].
Another methodological factor considered in the various studies analyzed for the development of SQIUs was the implementation of two types of databases: total database (TDB) and minimum database (MDB). It was observed that the use of the TDB did not present a significant 0.05, with a number of observations = 51] (Figure 7), which was very close to that obtained by the overall meta-analysis model. The use of the TDB approach has a disadvantage in terms of operation, as it requires a large amount of resources and time to sweep through multiple indicators [56]. To minimize this, it is recommended to employ personnel with knowledge of the study area to establish indicators that are specific to the area and the problems to be addressed [27]. The The use of the TDB approach has a disadvantage in terms of operation, as it requires a large amount of resources and time to sweep through multiple indicators [56]. To minimize this, it is recommended to employ personnel with knowledge of the study area to establish indicators that are specific to the area and the problems to be addressed [27]. The use of the MDB approach showed quality values similar to the general meta-analysis model; moreover, the use of an MDB with only a few indicators reduces the amount of resources and time and can be used for the implementation of processes to improve or maintain soil quality.
Subsequently, the methodological factor of the number of indicators that made up the SQI U developed in various studies was analyzed, as well as the possible existence of a minimum or maximum number of indicators, to obtain accurate observations on soil quality. This was based on the inherent complexity of the interpretation of the results obtained and the suggestions made by the "Soil Management Evaluation Framework", which suggests using a minimum of five indicators to analyze soil quality [57]. As the number of indicators that make up the SQI U increases, the interpretation of its results becomes more complex. For this reason, an analysis was conducted with the creation of two subgroups of observations: "equal to or less than five indicators SQI ≤5 " and "more than five indicators (SQI >5 )". The analysis showed that there was no significant difference (p > 0.05) lnRR , with the following results observed for the subgroups: SQI ≤5 [MD of −6.81%, 95% CI (−17.54 to 3.92%), p > 0.05, with a number of observations k = 40] and SQI >5 [MD of −8.63%, 95% CI (−18.75 to 1.50%), p > 0.05, with a number of observations k = 25] (Figure 8). It is important to note that using SQI U s with a smaller number of indicators can be more practical for field implementation because it reduces costs and the use of resources, and makes interpretation easier [12]. In contrast, SQI U s with a larger number of indicators can complicate implementation, increase costs, and make interpretation more difficult. The number of indicators used in an SQI U depends on the type of database used and, since there is no standardized database, it can be challenging to interpret the quality results of the lnRR . If an SQI U has a significant number of indicators, it is recommended to use the MDB.
Finally, the last methodological factor analyzed was the type of indicators that made up the SQI U s.  (Figure 9). The above result was found to be much lower than the one obtained in the general meta-analysis model. This suggests that using only chemical and biological indicators in an SQI U could potentially overestimate the negative impact of intensive agricultural practices on the soil, which can lead to incorrect actions being taken in terms of soil recovery, wherein priority is given to areas that are not severely damaged, while areas that require urgent action may not receive enough resources. Sarmiento et al. (2018) [12] concluded that the robustness of an SQI depends in part on the integration of the three types of indicators-chemical, physical, and biological-into the SQI. Regarding the other subgroups, they did not present a significant (p > 0.05) lnRR . This may be because the subgroups had a small number of observations, averaging at 13, and more studies would be needed to get a more precise observation of the effect of the different subgroups on the lnRR [21]. Therefore, it is necessary to implement a greater number of studies around the world focused on the establishment of SQIs to obtain a generalized view of soil quality. Agronomy 2023, 13, x FOR PEER REVIEW 14 of 22 , SQIU subgroup with less than or equal to five indicators; , SQIU subgroup with more than five indicators; MD, mean effect per observation; CI-95%, confidence interval at 95% ; df, degrees of freedom; p, probability value; , coefficient of heterogeneity; , variance; black diamond, ′ and overall model interval; black dotted line, no effect line; red dotted line, mean model effect line; black squares, mean value of effect plus variability per observation; ′ , average effect size. For references, please refer to Figure 6.

Metaregression
As there is no consensus on which indicators are most relevant to soil quality, a metaregression was conducted to analyze the correlation between the most commonly used indicators in the different SQIUs and the ′ . The indicators most closely related to the ′ are shown in Figure 10.

Metaregression
As there is no consensus on which indicators are most relevant to soil quality, a metaregression was conducted to analyze the correlation between the most commonly used indicators in the different SQI U s and the lnRR . The indicators most closely related to the lnRR are shown in Figure 10. The CEC indicator is related to the soil's ability to retain nutrients and make them available to crops. Soils with high CEC values typically have a good structure and high nutrient-holding capacity [58,59]. The C/N indicator provides information on the quality of OM present in the soils. High values of this indicator show a greater presence of C compounds, fibers, and cellulosic compounds, which could hinder the metabolization and The CEC indicator is related to the soil's ability to retain nutrients and make them available to crops. Soils with high CEC values typically have a good structure and high nutrient-holding capacity [58,59]. The C/N indicator provides information on the quality of OM present in the soils. High values of this indicator show a greater presence of C compounds, fibers, and cellulosic compounds, which could hinder the metabolization and mineralization of OM by the microorganisms present in the soil [60]. The Ca indicator plays an important role in soil structure, as it helps to stabilize aggregates and prevent degradation due to water and wind erosion, thus reducing the loss of OM [61]. However, high concentrations of Ca can have negative effects on microorganisms and crops, causing osmotic stress [62]. On the other hand, moderate levels of Ca are essential for plant growth and development, as it contributes to cellular processes such as enzymatic activation and photosynthesis [59,63]. The MBC indicator has been reported as an indicator that is very sensitive to the type of soil management and is related to the capacity of the microorganisms to metabolize, mineralize, and fix the OM added to the soil, thereby preventing its loss through leaching or erosion phenomena [24]. High values of MBC indicate a high concentration of OM in the soil that microorganisms can metabolize and retain, which helps improve nutrient cycling and create favorable conditions for the growth and development of microbial communities in the soil. Consequently, having a diverse and functional microbial community makes the soil more resistant to disturbances caused by human activities or natural phenomena, and maintains the biological function of the soil. The SND and CLY indicators are related to soil structure. High SND values may indicate nutrient leaching and lower OM concentration, which can lead to increased erosion risk [64,65]. On the other hand, high CLY values may indicate higher soil fertility due to the presence of OM and nutrients, but could also be a sign of poor drainage and water infiltration issues [66].
As shown in the results of the metaregression process, there are indicators belonging to the physical, chemical, and biological types, which highlights the use of all of them for the establishment of soil quality through the implementation of an SQI U .

Metamodeling
With the emergence of six indicators related to the lnRR that belong to three different categories-physical, chemical, and biological-it was decided to model them in order to observe which of them represented to a greater extent the variability of the soil quality lnRR resulting from the intensive use of agricultural soils. Table 2 shows the models that presented significant correlations. lnRR , average effect size (%); CEC, cation exchange capacity (meq 100 g −1 dry soil); C/N, C and N ratio; Ca, calcium (meq 100 g −1 dry soil); MBC, microbial biomass C (mg C mic kg −1 soil); SND, sand (%); CLY, clay (%); AICc, corrected akaike information criterion; BIC, Bayesian information criterion; R 2 , linear correlation of heterogeneity; p, probability value for residual heterogeneity.
The metamodeling resulted in the combination of the CEC and MBC indicators representing 97% certainty of the variability of the lnRR for the quality of the analyzed soils ( Table 2). It is worth noting that the model that presented the highest correlation with the lnRR -with lower BIC and RSME parameters-was made up of chemical and biological indicators. However, it is necessary to include physical indicators in the database to be implemented in the development of the SQI U to avoid overestimating the quality observations obtained to negative values. The developed model itself would not be a generalized SQI U to be implemented; it simply provides an idea of the set of indicators to be analyzed for the configuration of the database used in the development of the SQI U .

Conclusions
Soil quality has been negatively impacted by intensive agricultural practices. The use of a minimum database for analyzing soil quality has advantages such as accuracy, ease of use, and ease of interpretation. However, experts in the region to be analyzed should be consulted to select indicators that focus on the specific problems being monitored. The number of indicators did not significantly affect the results obtained from the SQI U s. Therefore, it is suggested that the use of an SQI U should involve a smaller number of indicators, which would have operational advantages, such as lower costs and easier interpretation. The most important indicators for soil quality in intensive agricultural practices were CEC and MBC. These should be included in the database along with other indicators selected by local experts when developing regional SQI U s. The SQI U is a useful and easily interpretable tool for understanding soil quality worldwide, providing a simple and generalized overview of soil quality conditions. Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/agronomy13082166/s1, Table S1: Studies that integrated the database.