Soil-based Ecosystem Services Research Articles

Disentangling soil-based ecosystem services synergies, trade-offs, multifunctionality, and bundles: A case study at regional scale (NE Italy) to support environmental planning

The explicit use of ecosystem services (ESs) assessments has been called as a way to guide environmental decision making, yet the promise of the ES approach lies behind its potential. A way to consolidate the approach could be to introduce some aspects into the ESs assessments which might have been neglected so far. Such aspects are mainly: (1) a focus on the complex ESs relations (such as synergies and trade-offs) that can impact the supply of multiple SESs (soil ecosystem services), and (2) focus on potential drivers of SESs relations. We applied bivariate and multivariate approaches to SESs indicators derived from a solid pedological knowledge of the Emilia-Romagna study area in NE Italy. We focused on 7 SES: (1) habitat for soil organisms, (2) filtering and buffering capacity, (3) contribution to microclimate regulation, (4) carbon sequestration, (5) food provision potential, (6) water regulation, and (7) water storage capacity. These SESs were estimated through a combination of point observations, and pedotransfer functions (PTF) estimates spatialised over the area of interest with geostatistical simulation techniques. We found that SESs bivariate spatial relations could be categorised mainly in three types of patterns at regional scale, either: (1) synergistic SESs relations dominating at the region level, (2) trade-offs dominating, or (3) both kind of relations more or less equally frequent. Interestingly, in some cases the dominant regional SESs relation switched at a local level, and such switch was driven by soil properties. For the multivariate case (>2 SESs), two main results are highlighted. First, the combination of properties of some soils is so characteristic that they conform a single SESs bundle, as in the case of the rich SOM soils of alluvial origin in the NE of the region with low agricultural productivity, but high value in regulating SESs. Secondly, some SESs such as potential food provision and water regulation are more important than others to determine locations with high multi-services value at a regional level. This suggests that attention must be paid when ascribing high multi-services value locations as this is not independent of SESs relations. Overall, our results highlight the importance of soils in the potential supply of ESs and show that SESs relations are useful in the implementation of the concept in environmental assessments.

Spatial variability of clay minerals in a semi-arid region of Turkiye

Clay minerals are the source of many chemical and physical properties that influence the provision of soil-based ecosystem services. This study aimed to identify the most significant soil characteristics contributing to the spatial variability of clay minerals in a semi-arid region of Turkiye. Additionally, the study assessed the predictive capabilities of Classification and Regression Tree (CART), Random Forest Regression (RF) and eXtreme Gradient Boosting Regression (XGBoost) in estimating soil clay mineral content. Smectite+vermiculite (SMVR) was the most abundant clay mineral in the study area, followed by illite and kaolinite. Hyperparameter tuning significantly improved model accuracy, with root mean square error (RMSE) reductions ranging from 2.53% to 97.3%. The machine learning algorithms demonstrated varying performances in spatial prediction accuracy. The RF model achieved the lowest RMSE (8.587%) and the highest R2 values (0.796) for predicting SMVR. The XGBoost outperformed other models for kaolinite (RMSE: 4.814%, R2:0.713) and illite (RMSE:7.368%, R2:0.613). Exchangeable cations, particularly magnesium (Mg) and calcium (Ca), were identified as crucial factors influencing the spatial distribution of clay minerals. Among these, M concentration had the strongest influence on predicting both SMVR (38.1%) and illite (26.3%). Conversely, for kaolinite prediction, Ca concentration played the most significant role (38.7%), followed by Mg (19.93%). In conclusion, this study demonstrates the effectiveness of machine learning models, particularly XGBoost which achieved the lowest RMSE for all clay minerals investigated. These models offer a valuable tool for predicting clay mineral content in the Kazova Plain. The findings highlight the importance of parent material, weathering processes, and specific soil properties, such as exchangeable cations, in shaping clay mineral distribution. This knowledge not only contributes to a deeper understanding of soil formation in semi-arid environments but also practical applications. For instance, by predicting the abundance of SVMR, known for its high cation exchange capacity land managers can develop targeted strategies for optimizing fertilizer application in the Kazova Plain.

Is the organic carbon-to-clay ratio a reliable indicator of soil health?

Climate action plans under the Paris Climate Agreement and other national commitments aimed at improving soil-based ecosystem services require the operational monitoring of soil carbon (C). The European Union is aiming to enhance soil health, and as part of the proposed Soil Monitoring Law, the European Commission recommends the monitoring of the soil C loss indicator among other soil health indicators. In this study, we evaluate the feasibility of the proposed soil C loss indicator by assessing its performance using the EU-wide 2009 LUCAS soil survey data. The proposed indicator is the soil organic carbon (SOC) to clay ratio, with a threshold value of 1:13. The results are also compared with the C stock changes reported by countries to the climate convention (UNFCCC). Our results reveal that the variation in SOC and clay content at European scale exceeds that of the data used to develop the proposed indicator. We also found that the variation in the SOC content was influenced not only by clay content but also by climate and land-use reflecting C input levels. Therefore, the defined threshold is inadequate for detecting degraded soils if the SOC and clay content are beyond the conditions used to establish the criteria. Furthermore, major discrepancies were observed between the soil carbon stock changes reported by the national greenhouse gas (GHG) inventories and the proportions of degraded soils identified by using the soil C loss indicator. We conclude that employing a single indicator such as SOC:Clay ratio with one threshold value for all soils across various land covers, management practices, and climatic conditions, as defined by the European Commission for the Soil Monitoring Law, is inappropriate for monitoring soil C loss.

Linkages between plant functional diversity and soil-based ecosystem services in urban and peri-urban vacant lots

Linkages between plant functional diversity and soil-based ecosystem services in urban and peri-urban vacant lots

Comparative analysis of soil quality indexing techniques for various tree based land use systems in semi-arid India

Inappropriate management of land use systems is one of the main factors that leads to soil quality degradation and its quantification is crucial to their sustainable utilization planning. The objective of the research is to evaluate how various land-use systems, viz., natural forest, tree plantations of Tectona grandis, Terminalia bellirica, Swietenia macrophylla, Artocarpus hirsutus, Melia dubia based agroforestry system, horticulture (Mangifera indica) and agriculture systems impact the soil physicochemical and biological characteristics in semi-arid climatic conditions of India. Principal component analysis followed by linear and non-linear scoring methods was employed to compute the soil quality index (SQI). The soil attributes viz., dehydrogenase activity, acid phosphatase activity, soil available nitrogen, potassium, calcium, porosity, and soil available iron emerged as significant indicators for assessing the soil quality. Among different SQIs, non-linear weighted SQI can efficiently assess soil quality. Based on the non-linear weighted SQI, the order of the systems studied was natural forest (0.973) > Swietenia macrophylla (0.756) > agroforestry (0.737) > agriculture (0.556) > Tectona grandis (0.416) > Terminalia bellirica (0.373) > Artocarpus hirsutus (0.343) > Mangifera indica (0.208). The study concludes that converting natural forests into different land-use systems deteriorated the soil quality. Identifying soil indicators will help rapidly diagnose soil degradation, assess soil-based ecosystem services, and design appropriate land management practices in the future.

AMF-SporeChip provides new insights into arbuscular mycorrhizal fungal asymbiotic hyphal growth dynamics at the cellular level.

Arbuscular mycorrhizal fungi (AMF) form symbiotic associations with the majority of land plants and deliver a wide range of soil-based ecosystem services. Due to their conspicuous belowground lifestyle in a dark environment surrounded by soil particles, much is still to be learned about the influence of environmental (i.e., physical) cues on spore germination, hyphal morphogenesis and anastomosis/hyphal healing mechanisms. To fill existing gaps in AMF knowledge, we developed a new microfluidic platform - the AMF-SporeChip - to visualise the foraging behaviour of germinating Rhizophagus and Gigaspora spores and confront asymbiotic hyphae with physical obstacles. In combination with timelapse microscopy, the fungi could be examined at the cellular level and in real-time. The AMF-SporeChip allowed us to acquire movies with unprecedented visual clarity and therefore identify various exploration strategies of AMF asymbiotic hyphae. We witnessed tip-to-tip and tip-to-side hyphal anastomosis formation. Anastomosis involved directed hyphal growth in a "stop-and-go" manner, yielding visual evidence of pre-anastomosis signalling and decision-making. Remarkably, we also revealed a so-far undescribed reversible cytoplasmic retraction, including the formation of up to 8 septa upon retraction, as part of a highly dynamic space navigation, probably evolved to optimise foraging efficiency. Our findings demonstrated how AMF employ an intricate mechanism of space searching, involving reversible cytoplasmic retraction, branching and directional changes. In turn, the AMF-SporeChip is expected to open many future frontiers for AMF research.

Increasing complexity of agroforestry systems benefits nutrient cycling and mineral-associated organic carbon storage, in south-eastern Brazil

Agroforestry systems are often promoted as solutions to address land degradation and climate change. However, agroforestry is an umbrella term for a large variety of systems and it is not clear how their degree of complexity influences their provision of soil-based ecosystem services, such as soil organic carbon (SOC) storage and nutrient cycling. Furthermore, a knowledge gap remains whether agroforestry systems perform equally well on all soil types. The objectives of this study were 1) to assess the links between agroforestry complexity, nutrient cycling and SOC fractions, and 2) to assess how soil texture influences these relationships in Brazilian agroforestry systems. We sampled 59 agroforestry plots across 30 sites in Sao Paulo state, Brazil, and 8 monocrop sites (6 pastures and 2 crop monocultures). The 38 sites represented a soil textural gradient, ranging from very sandy to very clayey (clay content range 25 – 620 g kg−1). An Agroforestry Complexity Index (ACI) was defined based on tree species richness, stem density and pruning management. Nutrient (N, P, K, Ca, Mg) and C contents were determined in litter and soil (0–30 cm depth) samples, and mineral-associated organic C (MAOC) and particulate organic C (POC) in soil samples were assessed as well. ACI was positively associated with C, N, P, Ca and Mg stocks in litter, and these litter nutrients were in turn positively associated with the corresponding soil nutrient stocks. Associations between soil nutrients and MAOC were stronger on sandy soils than on clayey soils, particularly for P, Ca and Cation Exchange Capacity (CEC). For POC, robust relationships with nutrients were only found on sandy soils. Structural Equation Models indicated causal relationships between agroforestry complexity, P and Ca cycling, and MAOC and POC stocks in topsoils. Our results indicate that nutrients effectively cycle from in situ mulch into plant-available soil pools and highlight the synergies between nutrient cycling and stable C stocks that can be achieved in complex agroforestry systems. These synergies seem to be particularly strong on sandy soils (<15 % clay).

Are preferences for soil-based ecosystem services driven by spatial phenomena?

Agricultural soils provide multiple ecosystem services that affect human well-being. Soils’ potential to provide these ecosystem services varies spatially. Socio-demographic and other drivers of environmental preferences are also spatially variable. Therefore, preferences for soil-based ecosystem services are likely to be spatially heterogeneous, which may result in different policy priorities across locations. Understanding this spatial heterogeneity of preferences is therefore essential to guide public policy to protect healthy soils. We present a study that combines explorative and hypothesis-driven approaches to understand the spatial heterogeneity of preferences for four soil-based ecosystem services: climate regulation, clean water provision, drought protection and flood protection. Based on the results of a discrete choice experiment conducted on a representative sample of the German public, we first use global and local spatial autocorrelation measures to test whether there are any obvious patterns in the spatial distribution of preferences. Second, we use spatial lag models to test a number of hypotheses to explain the observed preference heterogeneity. We particularly focus on the spatial variability of relevant phenomena such as floods, droughts or nitrate pollution of groundwater, and their effects on the studied preferences. Lastly, we compare the results from both approaches in order to see whether the identified patterns are consistent with each other. We find weak patterns of spatial heterogeneity, but our hypotheses are all rejected. This suggests that salience of relevant phenomena and individual affectedness do not have an effect on preferences for soil-based ecosystem services.

Condition of the soil resource and provision of ecosystem services from a Brazilian Oxidic soil under conventional and integrated livestock-based systems

Integrated crop-livestock-forest systems promote soil health and deliver more ecosystem services (ES) compared to conventional livestock systems, although most studies on the subject poorly describe the soil component of these systems. This preliminary studyassessed the condition of the soil resource and its role in the provision of ecosystem services from a Brazilian Oxidic soil under conventional and integrated livestock-based systems. Five systems were studied from pasture only to partial and fully integrated croplivestock-forest systems. All are located on one of the Brazilian Agricultural Research Corporation (Embrapa) Research Farms. Data for the quantification of the soil resource and ES for each of these livestockbased systems were retrieved from previous studiesand used to quantify soil health, the provision of food and fibre and climate regulation, as well as report their impacts on receiving environment. Soil organic carbon content, a key component of soil health, was higher in the most integrated system. Soil-based grass yields were lower in integrated systems due to competition for resources from the trees or space taken by crops but had the highest overall provision of food and fibre. Carbon sequestration by trees in the integrated systems offsets enteric methane emissions from beef production, and this ES contributes to mitigating climate change. Future studies should include analysis of all the naturalresources and a wider range of soil-based ecosystem services, along with impacts on receiving environmentsto provide a more complete picture of the performance of integrated livestock-based systems.

A global indicator of soil macroinvertebrate community composition, abundance and diversity

Macroinvertebrate communities are highly sensitive indicators of physical and chemical soil qualities. Their evaluation in field conditions is rather simple and could serve as proxy for soil-based ecosystem services that farmers and field technicians could use. We tested the hypothesis that an indicator of soil macroinvertebrate communities, with 14 taxonomic groups characterized at the order level, could be used in any region of the world to assess these communities. A synthetic indicator was calculated using data from 9 reference sites from tropical, subtropical and temperate regions, a set of 3694 data of the open access Macrofauna database and a new site for validation. Invertebrates were extracted with the standard ISO/TSBF methodology and characterized with a set of 14 large taxonomic units, plus an index of taxonomic richness and total density. At each of the 9 reference sites, 27 to 252 sample points, representing different types of plant covers and/or soil management options, were considered and compared with their respective local synthetic GISQ indicator values. These indicators, elaborated from Principal Component Analysis of the sample points data, are set to vary from 0.1 to 1.0 according to the composition, diversity and abundance of the community. Analyses showed great similarities among sites, with Factors 1 (from 21.9 to 36.9 % variance explained) expressing the overall abundance and diversity of the communities and Factors 2 (8.9 to 15.83 %) opposing sites with dominant soil ecosystem engineer populations (Earthworms, Ants, Termites and some Coleoptera) to sites dominated by litter transformer populations (Diplopoda, Isopoda and others). Indicator formulae were designed based on PCA analyses of each data set and a global formula was established with the Macrofauna database. Very high correlations (>0.95) were obtained among values calculated with local data set and the general formula calculated with data of the Macrofauna database, in 8 of 9 sites and in the validation site. We discuss the importance of having a single formula to transform data obtained with a simple standard field method in the building of public policies for soil-based ecosystem services payment.

Distribution Characteristics and Influencing Factors of Soil Biological Indicators in Typical Farmland Soils

Soil biodiversity drives soil-based ecosystem services and is an important indicator of soil health. To understand the responses of important soil biological indicators to different farmland use contexts, 72 fields in three agricultural regions of China were used as research objects. The distribution characteristics and the factors influencing six indicators (carbon and nitrogen contents of soil microbial biomass (MBC, MBN, respectively), soil respiration (SR), soil catalase activity (CAT), soil acid phosphomonoesterase activity (APA), and soil earthworms) were investigated using field monitoring and indoor analysis. The MBC, SR, CAT, and APA indicators showed significant differences among the regions (p &lt; 0.05). Correlation and redundancy analyses indicated that the important factors affecting MBC, MBN, and soil respiration were cation exchange capacity, total N, organic matter, hydrolytic N, and soil bulk density, whereas the important factors affecting APA and earthworms were total phosphorus, available phosphorus, and available potassium. None of these factors had a significant effect on CAT. Climatic conditions, soil types, and farmland practices all have complex impacts on soil biodiversity. The results showed that attention should be paid to improving the physical conditions of the soil and to increasing soil fertility levels when establishing sustainable farm management patterns.

Biochar application differentially affects soil micro-, meso-macro-fauna and plant productivity within a nature restoration grassland

Biochar is proposed as an option to sequester carbon (C) in soils and promote other soil-based ecosystem services. However, its impact on soil biota from micro to macroscale remains poorly understood. We investigated biochar effects on the soil biota across the soil food web, on plant community composition and on biomass production. We conducted a field experiment in a nature restoration grassland testing four treatments: two biochar types (herbaceous feedstock pyrolyzed at 400 °C or 600 °C – hereafter B400 and B600), and a positive (i.e. unpyrolysed biochar feedstock, hereafter Hay) and negative (no addition) control. Responses of plants and soil biota were evaluated one and three years after establishing the treatments.Soil pH and K concentrations increased significantly in the B600 treatment. Mite abundances were significantly higher in B400 whereas nematode abundances were highest in Hay (1st year) and lowest in B400 (3rd year). Other soil fauna groups (enchytraeids and earthworms) varied more between years than between treatments. Legume cover increased significantly in the biochar treatments but this effect was transient. Legumes, grasses and primary productivity also showed a statistically significant Treatment x Year interaction due to transitory effects that were no longer present by the 3rd year.Our results suggest that biochar produced from meadow cuttings and applied at the 10 t/ha rate cause transitory impacts on soil biota abundance and plant communities over the 3-year timeframe used for this experiment. Therefore, this type of biochar could potentially be used for soil carbon sequestration, with minimal impacts on soil biota abundance or diversity, within the groups studied here, or plant biodiversity and productivity. Further research is required to investigate the longer-term impacts of this potential soil C storage sink.

Investigating preferences for soil-based ecosystem services

Abstract Soil provides multiple benefits for human well-being that are largely invisible to most beneficiaries. Here, we present the results of a discrete choice experiment on the preferences of Germans for soil-based ecosystem services. In an attempt to reduce complexity for respondents, we express soil-based ecosystem service attributes relative to the site-specific potential of soils to provide them. We investigate how knowledge about soils, awareness of their contributions to human well-being, and experience with droughts and floods affect preferences. We find substantial yet heterogeneous preferences for soil-based ecosystem services. Only some measures of familiarity exhibit significant effects on preferences.

Development and spatialization of a soil potential multifunctionality index for agriculture (Agri-SPMI) at the regional scale. Case study in the Occitanie region (France)

Europe is faced to a competing demand for areas suitable for agriculture and areas suitable for urban development. To help building the general guidelines necessary for a sustainable development at the regional scale, we developed a soil potential multifunctionality index for agriculture (Agri-SPMI), which used the concepts of soil-based ecosystem services and soil functions. It considered the provision of a physical and chemical habitat for plant growth functions, the retention and transfer of water and pollutants function, and the carbon storage function. Three agricultural land uses were considered in this index: arable lands, perennial crops, and pastures. The assessment of the Agri-SPMI was done over a large territory, using a small set of inherent soil properties available in the French soil survey database at 1:250,000 scale. The output was in the form of two maps, which can be viewed in a GIS software: a map of areas with high Agri-SPMI and a map of the service of food, biomass, and fiber production. The concept of soil multifunctionality allowed highlighting environmental issues related to agriculture and to encourage land planners to take them into account in addition to the agronomic potential of soils. The framework used is fully transferable to other scales, other soil contexts, or other ecosystem services, with adaptation of the soil properties considered.

Assessment of joint soil ecosystem services supply in urban green spaces: A case study in Northern Italy

The multiple functions of urban soils secure the supply of ecosystem services to the urban population, but they are seldom taken into account in current urban planning. The purpose of this study is to highlight the multifunctionality of the soils of urban green spaces and to assess the influence of different types of urban green (i.e. parks and gardens, roadside green and agricultural fields) on the joint supply of soil-based ecosystem services. In a case study area in the city of Carpi (NE Italy), we focused on a set of soil functions underpinning ecosystem services and on possible synergies and trade-offs between them. We surveyed and sampled 19 urban green areas to assess the following soil functions: biological fertility, potential habitat for organisms, water regulation and storage, soil buffering capacity and carbon stock. Results showed differences and trends in the bundles of ecosystem services provided by the soils of urban green areas, highlighting the relevance of soil disturbance and vegetation cover density in affecting soil functions. For biological indicators, results showed a negative significant correlation to CaCO3 content, which is associated with the degree of soil disturbance, and that soils in urban areas do not always have compromised soil fauna and may provide the same level of biological quality as agricultural soils or forests.

Soil properties after one year of interseeded cover cropping in topographically diverse agricultural landscape

Planting cover crops within or following a cash crop may improve soil-based ecosystem services due to increased plant diversity and a longer duration of live vegetation coverage. We examined the effect of three different cover cropping systems on soil properties after one year of a three-year organic transition rotation at four agricultural field sites with contrasting topographical positions, namely depressions, slopes, and summits. The four studied systems were (1) cereal rye (Secale cereal L.) planted after corn (Zea mays L.) harvest (Rye); (2) a mixture of cold susceptible cover crop species, namely, oat (Avena sativa), winter pea (Pisum sativum), and radish (Raphanus sativus), interseeded into corn (WK); (3) a mixture of cold tolerant cover crop species, namely, annual ryegrass (Lolium multiflorum), Dwarf Essex rapeseed (Brassica napus), and crimson clover (Trifolium incarnatum), interseeded into corn (WH); and (4) a no-cover control (NC). While soil moisture was affected by topography, interseeding cover crops into corn did not influence soil moisture levels at the 0–10 cm depth for the studied year. Soil NO3- content was markedly higher in the WK system compared to cereal rye and WH cover crop treatments. The difference was especially pronounced in depressions and summits. Soil N mineralization rates followed the pattern WH>WK>Rye>NC and the effects were most pronounced in slopes. Soil microbial biomass C was the highest in depressions followed by summits and slopes, and in depressions the WH had higher microbial biomass than the other systems. There were no effects of cover crops and topography on soil C mineralization one year after the organic transition was initiated. The WH system increased the fraction of 0.053–2 mm aggregates and decreased > 2 mm aggregates in depressions. The results suggest that the effect of cover cropping can become evident already one year after the organic transition begins. The WH mixture interseeded into the cash crop was an optimal cover crop choice for improving soil characteristics as well as decreasing soil N leaching risks during organic transition in undulating agricultural terrain. However, the magnitude of the benefit provided by WH was mediated by topography.

Editorial: Assessment and Modeling of Soil Functions or Soil-Based Ecosystem Services: Theory and Applications to Practical Problems

Editorial: Assessment and Modeling of Soil Functions or Soil-Based Ecosystem Services: Theory and Applications to Practical Problems

Impact of Environmental Conditions and Management on Soil Arthropod Communities in Vineyard Ecosystems

The importance of soil biodiversity and soil-based ecosystem services in the context of viticulture has recently been emphasized. Over 85% of soil fauna species richness is represented by edaphic arthropod communities. Edaphic arthropod responses to soil characteristics and management practices can be considered as good bioindicators of soil quality. Here, 168 soil samples that were collected from 2014 to 2019 in several vineyards of different Italian wine-growing areas were analyzed to explore how arthropod communities respond to several factors that are characteristic of vineyard ecosystems. The analysis of the combined effects of the primary abiotic variables (the chemical and physical characteristics of soil) and management practices (organic vs. conventional, soil inter-row management) on soil biological quality (assessed by QBS-ar index) identified soil temperature and soil texture as the abiotic factors exerting the most significant effect on the QBS-ar values. Organic vineyards exhibited higher QBS-ar values compared to those of conventionally managed vineyards, and subsoiling negatively influenced the soil biological quality.

Attitudes and Preferences towards Soil-Based Ecosystem Services: How Do They Vary across Space?

Soil ecosystem services (ES) provide multiple benefits to human well-being, but the failure to appreciate them has led to soil degradation issues across the globe. Despite an increasing interest in the threats to soil resources, economic valuation in this context is limited. Importantly, most of the existing valuation studies do not account for the spatial distribution of benefits that soil ES provide to the population. In this study, we present the results of a choice experiment (CE) aimed at investigating spatial heterogeneity of attitudes and preferences towards soil conservation and soil ES. We explored spatial heterogeneity of both attitudes and welfare measures via GIS techniques. We found that citizens of the Veneto Region (Northeast Italy) generally have positive attitudes towards soil conservation. We also find positive willingness-to-pay (WTP) values for soil ES in most of the study area and a considerable degree of heterogeneity in the spatial taste distribution. Finally, our results suggest that respondents with pro-environmental attitudes display a higher WTP based on the geographic pattern of the distribution of WTP values and attitudinal scores across the area.

Land-use change affects soil hydro-physical properties in Mollisols of semiarid Central Argentina

The objective of the present study was to describe the relationships between soil physical properties and to identify variables that determine the pore system and hydrological functions which are specifically important in drought prone semiarid marginal croplands. We compared Mollisols under natural vegetation (NV) with their agriculture (AG) pairs in the central Argentinean Dry Pampas to obtain the range of values for variables that characterize the soil pore system. These were total carbon (OC), bulk density (BD), structural instability index (SII), volumetric weight of aggregates (VAW), total porosity (TP), macropores (Ma), saturated hydraulic conductivity (K), infiltration rate (IR), Dexter's S index (S), and least limiting water range (LLWR). To a depth of more than 23 cm, TP and Ma were considerably lower (~ 35% and 50% less, respectively) in AG versus NV. The changes in porosity and pore size resulted in 74% smaller IR and 65% lower K. The S index had mean values of 0.16 and 0.14 for the upper and middle layer respectively in NV; while for AG, these values were significantly lower (0.09 and 0.10, respectively). The LLWR was also lower in AG (21 and 59% less in the upper and middle layer, respectively), due to high soil resistance even at non-limiting BDs (1.17 and 1.32 g cm−3 for AG2 and AG1, respectively). Soil OC was positively correlated with TP, Ma, and S, and negatively correlated with BD and SII. Multiple regression models for S included SII, TP, and clay contents, while for LLWR the variables were OC, TP, and clay plus silt. Thus, these complex soil quality indices can be predicted by using simpler physical parameters, but texture and organic matter also have a strong effect on them. Moreover, our results contribute to a better understanding of the role of texture and organic matter in controlling soil physical quality and for maintaining soil-based ecosystem services.