Effects of tillage intensity on pore system and physical quality of silt-textured soils detected by multiple methods
Thomas Weninger A E , Janis Kreiselmeier B C , Parvathy Chandrasekhar B C , Stefan Julich C , Karl-Heinz Feger C , Kai Schwärzel B , Gernot Bodner D and Andreas Schwen AA Institute for Soil Physics and Rural Water Management, University of Natural Resources and Life Sciences, Muthgasse 18, 1190 Vienna, Austria.
B UNU Institute for Integrated Management of Material Fluxes and of Resources (UNU-FLORES), United Nations University, 01067 Dresden, Germany.
C Institute of Soil Science and Site Ecology, Faculty of Environmental Sciences, TU Dresden, 01737 Tharandt, Germany.
D Division of Agronomy, Department of Crop Sciences, University of Natural Resources and Life Sciences, Vienna, Konrad Lorenz-Straße 24, 3430 Tulln, Austria.
E Corresponding author. Email: thomas.weninger@boku.ac.at
Soil Research 57(7) 703-711 https://doi.org/10.1071/SR18347
Submitted: 23 November 2018 Accepted: 29 April 2019 Published: 1 July 2019
Journal Compilation © CSIRO 2019 Open Access CC BY-NC-ND
Abstract
Understanding the effects of agricultural management practices on soil functionality is an ongoing challenge in environmental science and agricultural practice. In the present study we quantified the effects of changes in tillage intensity on soil physical quality and pore size distribution after 6, 10 and 23 years. At three long-term tillage experimental sites in central Europe we analysed soils under four different soil management systems: conventional mouldboard tillage; chiselling + rotary harrow; rotary harrow; and no till. These treatments differed in mechanical intensity and depth. Pore size distributions were calculated from soil water retention curves based on high-resolution measurements. Subsequently, fractions of functional pore size classes and indicators of soil physical quality were determined and compared between the treatments. In addition, we evaluated the performance of two calculation approaches for pore size distribution: (1) fitting of a smoothing cubic spline; and (2) a bimodal van Genuchten function. The parametric function yielded a higher proportion of storage pores by approximately 3–5%. The combination of multiple measurement and evaluation methods enabled detailed comparison of soil physical characteristics between different tillage treatments. No-till soils showed a distinct lack of transmissive pores and higher bulk density, but similar plant-available water capacity, compared with the other treatments. Under all soil management systems, aeration deficits were observed, emphasising the high vulnerability for compaction of silt-dominated arable soils with a low organic matter content. Hence, the design of agricultural soil management strategies on such soils needs to consider the risks of compaction as thoroughly as erosion or chemical degradation.
Additional keywords: high-resolution measurements, hydraulic soil properties, pore size distribution, soil degradation, soil management.
Introduction
Pores, with their considerable different sizes and shapes, are the reaction zone of soil (Gupta et al. 2008). The soil pore size distribution (PSD), which is the fraction of distinct pore sizes in the whole pore system, determines the functionality of soil. This functionality is essential for civilisation and ensures the production of food, the cleaning of waste water or the storage of water and carbon. Especially valuable in times of changing climate and increasing demands on soil and water resources is the ability of soil to store water. Consequently, the loss of functional soil pore space, more generally physical soil degradation, is an actual global threat (Lal 2000; Eswaran et al. 2001; Commission of the European Communities 2006). Agriculturally used soils are most vulnerable and need to be protected and regenerated by adapted soil management strategies (FAO 2017).
The soil pore system is a combination of textural and structural pores (Lal and Shukla 2004). The former are generic voids between primary soil particles, whereas the latter originate from secondary soil formation processes, like aggregation, shrinking–swelling, freezing–thawing or biological activity. Of these two subsystems, only structure can be changed by soil management in order to improve the functionality of the soil, especially with a focus on water (Rabot et al. 2018). In arable land, optimisation of the tillage regime is thought of as means of improving soil structure. Hydrological and physical soil properties that are affected by tillage include water storage capacity, hydraulic conductivity or aggregate stability. However, for a meaningful evaluation of the advantages and restrictions of different tillage intensities, appropriate methods are needed that allow detailed analysis of the soil pore system and its functionality. Commonly used methods are either direct measurements from images, obtained from computed tomography (Grevers et al. 1989; Anderson and Hopmans 2013) or thin sections (Kubiena 1938; Elliot and Heck 2007), or indirect derivation from the soil water retention curve (SWRC; the relationship between soil water head h and volumetric water content θ) via the law of capillarity (e.g. Lal and Shukla 2004). The direct methods allow only analysis of the macropore system on reasonably sized samples due to limits in resolution (Schlüter et al. 2018). Nevertheless, a broad understanding of processes can be derived from such images, including the shapes of pores, organic particles and similar characteristics. The informative value of quantitative image analysis is highly dependent on methodological aspects such as resolution and image analysis algorithms (Anderson and Hopmans 2013). Conversely, indirect methods are based on model simplifications, which may imply uncertainty in the outcomes of soil physical studies (Kutilek 2004; Tuller and Or 2004). Indirect methods are commonly used at larger scales because they do not require expensive instruments and are less time consuming than direct methods. The results of such measurements provide valuable quantitative information about the effects of certain processes on linked soil functions and can be used directly for soil physical modelling.
The effects of different soil and land management strategies on soil pore characteristics and soil hydraulic properties have been the subject of extensive research (Horel et al. 2015; Schjønning et al. 2017; Blanco-Canqui and Ruis 2018). Different results have been reported as a consequence of tremendous variability in climate, soil texture and differences in soil management systems. Pires et al. (2017) characterised changes in the pore system of no-till (NT) and conventionally tilled (CT) soils using three-dimensional microscale computed tomography, micromorphological analyses of thin sections of impregnated blocks and indirect determination of PSD via the SWRC and fitting of a cubic spline. On clay soil in subtropical climate, Pires et al. (2017) found a larger, more complex pore system after 26 years of NT, supposedly caused by higher biological activity than under CT. In another study, only minor beneficial effects of NT were shown by Wairiu and Lal (2006), who used mercury porosimetry to detect PSD in two 38-year-old tillage experiments on silt loam. In contrast, Peña-Sancho et al. (2017) found lower porosity and a lack of macropores under NT as well as considerable seasonal variability in PSD under CT and reduced tillage (RT) on loam soil (23-year tillage experiment; SWRC derived from pressure plate apparatus). Kodesova et al. (2011) characterised macropore structure (h > –70 cm) under CT and 30 years grass coverage by inverse simulation of a multistep outflow experiment and micromorphological images of thin sections. Grassland soils, as an extreme case of NT, had more capillary pores and matrix pores than soils under CT. Similarly, Schwen et al. (2011) analysed the macroporosity of a silt loam in eastern Austria by inversely simulated tension infiltrometer experiments, finding found higher conductivity after 11 years of NT due to highly connective and less tortuous macropores compared with CT and RT, whereas CT and RT showed high seasonal variability in soil physical properties.
More generally, the review by Horel et al. (2015) summarised changes in soil hydraulic properties after a change in land use or soil management. As a main outcome, Horel et al. (2015) stated that negative effects, such as a decrease in plant-available water, increased bulk density and loss of soil organic matter, may be expected with intensifications in soil disturbance. The transformation of soil properties after a change in soil management strategies may last for several years or decades, especially on heavily textured soils. Nevertheless, Horel et al. (2015) also concluded that robust statements about the effects of distinct management strategies on the physical properties of soils are hampered by inconsistencies in scientific and agricultural methodology or the heterogeneity of soil.
Similarly, the analytical methodology used in the studies cited above is often not consistent. To achieve continuous results over the whole soil moisture range, a mathematical function needs to be approximated to measurements on the SWRC. To this end, two different approaches are primarily used: (1) a parametric function with the aim of condensing the information into a preferably small number of parameters; or (2) an interpolation via more complex functions, like a cubic spline defined by a vast number of parameters (Othmer et al. 1991). The latter method has the advantage of yielding approximations closer to measurements, and consequently allows a more precise interpretation of the resulting PSD. In contrast, parametric functions are the standard way to input soil physical properties into simulation models for soil water dynamics. Commonly, simple S- or C-shaped functions are fitted to the data. However, most soils show a more complex PSD, and the lack of flexibility in the function used may cause low goodness of fit and hamper the interpretation of differences between certain treatments (Durner 1994; Lozano et al. 2016). Hence, bi- or multimodal parametric functions are increasingly being used (Romano and Nasta 2016; Reynolds 2017). In this study we compared the two approaches using data measured by the evaporation method (Schindler et al. 2010) combined with dewpoint hygrometry (Campbell et al. 1973). We examined whether the choice of approximation scheme affects statements about soil physical quality (SPQ; Reynolds et al. 2009) and fractions of functional pore sizes (Lal and Shukla 2004), both of which were derived from the PSD.
The main objective of the present study was to quantify the changes in PSD in arable fields after a conversion from mouldboard ploughing to conservation tillage or NT. Consequently, the findings of the study would add information to the incomplete picture about the implications of different-intensity tillage strategies on the physical constitution of soil. The experimental locations were representative of silt-dominated arable soils under a temperate climate, which account for the biggest part of crop production in central Europe. To obtain highly informative and valuable results we applied a unique combination of different high-resolution measurement methods and evaluation approaches.
Materials and methods
Sampling sites and procedure
Soil sampling and experiments took place at three long-term tillage trials in north-eastern Austria and Saxony, Germany (Table 1). Four different tillage treatments had been established on the fields for 6, 10 and 23 years: (1) CT with a mouldboard plough and rotary harrow; (2) RT with a chisel plough and rotary harrow; (3) minimal tillage (MT) using only a rotary harrow; and (4) NT with a direct seeder. Intensive tillage operations (mouldboard and chisel plough) were conducted after harvest; before seeding of the subsequent cash crop, only the rotary harrow or direct seeder were used. Undisturbed soil samples were collected in steel cores (250 cm3; inner diameter 8.4 cm, height 5 cm) at the soil surface, as were disturbed samples adjacent to the cores. Sampling was conducted between three and five times at the three sites throughout the vegetation period in 2016 (Table 2). Volumetric water content, θ, before sampling was between 0.20 and 0.35 cm3 cm–3; hence, the soil was neither near saturation nor were drying cracks present on the soil surface, and the sampling points were placed in inter-row spaces to avoid machine tracks.
Measurement methods
The data for the retention curves (data pairs of θ and h) were obtained using the evaporation method and a HYPROP device (METER Group, Munich, Germany; Schindler et al. 2010) and a dew point hygrometer (WP4C PotentiaMeter; METER Group). In addition, hood infiltrometer (Schwärzel and Punzel 2007) experiments were conducted to measure hydraulic conductivity in the near-saturated range, bulk density was measured by oven drying (105°C, 24 h) of 250-cm3 core samples and saturated hydraulic conductivity was measured in the laboratory using the falling head method (Reynolds and Elrick 2002). The sampling and measurement procedures have been described in detail elsewhere (Weninger et al. 2018). Failure in one of the methods used was inevitable in single cases, and in such cases data from the whole respective dataset (i.e. sampling point) had to be excluded from further analyses. This explains the discrepancy between maximum possible data extent calculated from Table 2 and actual data extent in Table 3.
Data processing and statistics
Retention curves for each experimental plot were approximated to measured data using two different approaches: a cubic spline and a bimodal van Genuchten (bVG) function (Othmer et al. 1991). First, values for soil water head h were transformed to pF values (pF = log10(h)). The density of data points over the pF range was very heterogeneous, hence data were classified with a class width of 0.2 pF to balance weights of measurements. The means of all data inside the respective classes were used as nodes for the fitting. The high flexibility of the cubic spline led to implausible oscillations in the resulting curve in the transition zone between the measurement ranges of the two combined methods. Consequently, we used a smoothing cubic spline (sCUB) to balance these irregularities (R Core Team 2016; smoothing parameter = 0.5). The fitting yielded two continuous functions (sCUB and bVG) for each sampled point, which were used to predict θ for the centres of pF classes (width 0.2, as above). Subsequently, the predicted class centre values of all corresponding replicates (same site and treatment) were pooled (the sample size at this point is given in Table 3, column n) and the mean, s.d. and CV were calculated and analysed.
Three different outcomes were used to evaluate the results: (1) graphical interpretation of curves; (2) classification into functional pore size classes after Greenland (1981); and (3) capacity-based indicators for soil physical quality (Reynolds et al. 2009). For the graphical interpretation of curves, water content data (predicted as described above) was normalised by the maximum observed water content, θmax (approximation for saturated water content, θs), yielding effective saturation, Se(h). Results for the same site and treatment were averaged and plotted together with the first derivative of Se(h), which corresponds to the PSD.
For the second and third outcomes, absolute values for θ were used and the conversion between soil water head h on the retention curve and the corresponding pore diameter was made using the simplified capillary equation h = 1490r–1 (where h is soil water head (cm) and r is the equivalent pore radius (µm)). Following Greenland (1981), pores with a calculated diameter >500 µm were classified as fissures, those with diameters between 50 and 500 µm were classified as transmissive pores, those with diameters between 0.5 and 50 µm were classified as storage pores, those with diameters between 0.005 and 0.5 µm were classified as residual pore, and those with a diameter <0.005 µm were classified as bonding pores. The predicted values for the borders of classes were interpolated, the fractions of pore volume in the different pore classes were calculated and the results for the same site and treatment were averaged and analysed statistically for differences between treatments. For the third evaluation approach, soil physical parameters and capacity-based indicators for SPQ (Reynolds et al. 2009) were calculated according to Table 4. The selected indicators are widely used for interpretation of SPQ and ensure comparability to similar studies.
All data was processed using R version 3.3.1 (R Core Team 2016). Normality or log-normality of selected results for comparison of treatment effects was tested by visual interpretation of Q–Q plots. Saturated hydraulic conductivity values derived by falling head laboratory method (Ks(FH)) and in the field by hood infiltrometer experiments (Ks(HI)) were represented best by a log-normal distribution; other metrics followed a normal distribution. The significance of differences between treatments was analysed by analysis of variance (ANOVA), with Tukey’s honestly significant difference (HSD) test used as post hoc test to define coherent groups (de Mendiburu 2016). A paired t-test was used to detect significant differences between results derived by the two approximation approaches.
Results and Discussion
Effects of tillage intensity on physical soil conditions
Results of measured soil physical properties are given in Table 3, and SPQ indicators calculated from the retention curve are given in Table 5. Ks(FH) and Ks(HI) enable inference of the presence of connective macropores (fissures and transmission pores according to the classification used) because the water flow measured predominantly occurs through these pores. No consistent differences between treatments were found, but most of the lowest conductivities were found under NT systems and highest conductivities under the RT and MT systems (Table 3). These findings are in agreement with the review of Blanco-Canqui and Ruis (2018), who did not find systematic effects of tillage on saturated hydraulic conductivity. Nevertheless, Blanco-Canqui and Ruis (2018) identified more distinct differences for infiltration capacity, which was highest under NT. Furthermore, NT soils showed higher values of bulk density, ρd, than other treatments at two of the three sites (Table 3), which is in accordance with 24 of 62 available studies reviewed by Blanco-Canqui and Ruis (2018). However, in medium-textured soils (the term used in comparison with sandy and clayey or fine-textured soil by Blanco-Canqui and Ruis 2018, hence comparable to the soils sampled herein), 18 of 24 studies reviewed by Blanco-Canqui and Ruis (2018) found significantly higher ρd under NT. This agreed with the fact that medium-textured or silt-dominated soils are especially vulnerable to compaction by clogging of pores (Horn et al. 1995). In the present study results, air capacity (AC) was also distinctly lower under NT, whereas plant-available water capacity (PAWC) was in the same range as for all other treatments. In contrast, no systematic limitations in AC under NT were found by Reynolds et al. (2009), who analysed 13 soils containing less silt, and Lozano et al. (2016) on an Argentinian loam with soil organic carbon (SOC) content of 40–56 g kg–1. Consequently, higher SOC content could decrease the vulnerability of the sampled soils to compaction, and an improvement may be expected after a longer period of NT management (Murphy 2015; Blanco-Canqui and Ruis 2018).
Classifying SPQ indicators in terms of their agricultural usability (Reynolds et al. 2009) showed that, regardless of treatment, all soils sampled except one (ρd = 1.14 g cm–3) had too-high ρd, with values ranging from 1.25 to 1.51 g cm–3 (Table 3; the optimal ρd range for loamy soils is between 0.9 and 1.2 g cm–3; Reynolds et al. 2009). Similarly, results for relative field capacity (RFC) were between 0.724 and 0.900 cm3 cm–3 (except for one soil, in which RFC was 0.673 cm3 cm–3), indicating potential yield losses due to a lack of aeration (Table 5; the optimal RFC range is between 0.6 and 0.7 cm3 cm–3). Results for AC followed the same trend even though the variability between treatments and sites was higher, and certain soils could be denoted as optimal (>0.14 m3 m–3), whereas in others aeration was poor (<0.10 m3 m–3). The PAWC was at least 0.203 m3 m–3, hence ideal in all soils (optimal range >0.20 m3 m–3). This may also be explained by the high fraction of silt together with low organic matter content (e.g. Horn et al. 1995).
Pore size distribution
All SWRCs and their derived PSDs are compared in Figs 1 and 2. Interpretations were based on sCUB (Fig. 1) because it represented measured data better than bVG. The most intensive treatment CT showed a distinct bimodal character for PSD at all sites. In all three tillage experiments there was an obvious difference between NT and all other treatments, which showed a mode (peak) in the range of transmissive or course storage pores. This mode was not present in the NT treatment, which was surprising because NT was expected to have more biologically built macropore, such as earthworm burrows or decayed roots (e.g. Wairiu and Lal 2006; Schwen et al. 2011; Alvarez et al. 2014). Fig. 3 shows quantitative comparison of pore volume belonging to functional pore size categories. On average, NT soils had 0.082 m3 m–3 less transmission pores than the other treatments (referring to total pore volume; s.d. = 0.095). In contrast, the NT treatment was richer in residual pores. Results from comparable studies using a coarser, size-based classification were variable without a uniform trend (Blanco-Canqui and Ruis 2018). In contrast with the present study, significantly lower fractions of storage and transmission pores were found in CT than in several types of conservation tillage by Pagliai et al. (2004) and Abdollahi and Munkholm (2017).
The variability between PSD curves of the same treatment was exceptionally low, as evidenced by the resulting indicators for SPQ (Table 5) and functional pore size classification (Fig. 3). Especially for CT, we expected comparably higher temporal variability due to seasonal changes (e.g. Tebrügge and Düring 1999; Kargas et al. 2016; Soracco et al. 2018). Most likely the low variability between single samples within the same treatment resulted from the normalisation with θmax. Consequently, in future studies measurements will be adapted and extended to analyse the seasonal variability of the characteristics examined for multiple vegetation periods. The focus of the present study was on long- and mid-term changes (6, 10, 23 years) of soil physical conditions due to different tillage intensity, and the low variability during the sampled period strengthens the results and interpretations.
Effects of fitting procedure
Soil physical parameters derived from retention curves can be affected by the type of fitting procedure applied to the measured θ(h) data points. We tested two different approaches to derive continuous retention curves and quantified the differences in the evaluation metrics (Table 5). The differences were significant for all indicators, as analysed by paired t-tests (α = 0.05). Using bVG rather than sCUB resulted in an overestimation of PAWC and an underestimation of AC of <2%, whereas RFC was overestimated by approximately 3%. This resulted in a higher proportion of samples where AC was classified as poor. In contrast, the resulting classification of PAWC (ideal) and RFC (potentially lacking aeration) did not differ between the two approaches. The fitting of the bVG function yielded a higher variability in SPQ metrics that was the result of lower goodness of fits compared with the more flexible and data-driven sCUB (Othmer et al. 1991; Kastanek and Nielsen 2001). Nevertheless, the effective differences between the two approaches after pooling replicate data were distinctively smaller than on single measurements as reported by Othmer et al. (1991). Hence, appropriate replication increases the possibility of detecting differences in PSD or SPQ between certain soil management regimes also using the less flexible parametric function (bVG). Bimodal functions should be used, because the lack of flexibility in unimodal functions may hamper the detection of such differences in certain soils (Lozano et al. 2016).
The most distinct differences between retention curves derived by the two different functions were visible at Site C (Obersiebenbrunn, Austria). There, bVG exhibited clear bimodality with a second mode in the range of residual pores (Fig. 2). The sCUB at Site C followed an irregular unimodal shape, whereas other sites were more similar to bimodality (Fig. 1). This may be interpreted as tendency of bVG to yield poor fits to irregularly shaped retention data because there is a certain probability that a second mode is detected even if it is not fully supported by measurements. This is in contrast with the results of Romano and Nasta (2016), who found benefits in using bimodal functions even for weakly bimodal PSD. Fig. 3b again shows that using the parametric bVG function overestimated storage pores, which determine PAWC, compared with sCUB. However, considerable differences regarding larger pores responsible for aeration were only found at Site C, where the overall pattern showed least bimodality in pore size distribution.
Conclusions and relevance
The use of a cubic spline function was found to be preferable to the parametric bVG model due to higher flexibility, especially for irregularly shaped data. We also found significant differences between the two approaches when used as a basis for a functional classification of the soil pore system.
The soils analysed were silt dominated and showed a lack of aeration, whereas water capacity was optimal. The higher bulk density of the NT system compared with the other tillage treatments was primarily related to a lower fraction of large, aerated pores. Graphical representation of PSD revealed this lack of transmissive pores in NT soils, whereas NT soils had higher PAWC. Additional analysis of indicators for SPQ confirmed these interpretations, and low organic matter content increased the vulnerability of the soil to compaction. Although silt-dominated soils are most endangered by erosion and thus target sites for RT systems, a loss of aeration has to be controlled to avoid adverse conditions for root aeration, plant growth and consequently yield. The results presented herein are representative of a significant portion of the agricultural soils in central Europe, and the methodology allows sound and detailed interpretations. Future advances in modelling of changes in soil pore systems based on such results will improve the opportunities for the development of sophisticated agricultural management strategies.
Conflicts of interests
The authors declare no conflicts of interest.
Acknowledgements
The authors thank the administration of agricultural schools in the province of Lower Austria for supporting the field experiments and providing meteorological data. The Austrian part of this project was funded by the Austrian Science Fund (FWF; I2122-B16) and the German part of the study was funded by the Deutsche Forschungsgemeinschaft DFG (SCHW 1448/6-1, FE 504/11-1).
References
Abdollahi L, Munkholm LJ (2017) Eleven years’ effect of conservation practices for temperate sandy loams: II. Soil pore characteristics. Soil Science Society of America Journal 81, 392–403.| Eleven years’ effect of conservation practices for temperate sandy loams: II. Soil pore characteristics.Crossref | GoogleScholarGoogle Scholar |
Alvarez CR, Taboada MA, Perelman S, Morrás HJM (2014) Topsoil structure in no-tilled soils in the Rolling Pampa, Argentina. Soil Research 52, 533–542.
| Topsoil structure in no-tilled soils in the Rolling Pampa, Argentina.Crossref | GoogleScholarGoogle Scholar |
Anderson SH, Hopmans JW (2013) ‘Soil–Water–Root Processes: Advances in Tomography and Imaging.’ SSSA Special Publication 61. (Soil Science Society of America (SSSA): Madison, WI) https://doi.org/10.2136/sssaspecpub61
Blanco-Canqui H, Ruis SJ (2018) No-tillage and soil physical environment. Geoderma 326, 164–200.
| No-tillage and soil physical environment.Crossref | GoogleScholarGoogle Scholar |
Campbell EC, Campbell GS, Barlow WK (1973) A dewpoint hygrometer for water potential measurement. Agricultural Meteorology 12, 113–121.
| A dewpoint hygrometer for water potential measurement.Crossref | GoogleScholarGoogle Scholar |
Commission of the European Communities (2006) Communication from the Commission to the Council, the European Parliament, the European Economic and Social Committee and the Committee of the Regions: thematic strategy for soil protection. Technical report. Available at https://eur-lex.europa.eu/legal-content/EN/TXT/?uri=CELEX:52006DC0231 [verified 8 November 2018].
de Mendiburu F (2016) agricolae: statistical procedures for agricultural research. R package version 1.2–4. Available at https://cran.r-project.org/web/packages/agricolae/index.html [verified 13 June 2019].
Durner W (1994) Hydraulic conductivity estimation for soils with heterogeneous pore structure. Water Resources Research 30, 211–223.
| Hydraulic conductivity estimation for soils with heterogeneous pore structure.Crossref | GoogleScholarGoogle Scholar |
Elliot TR, Heck RJ (2007) A comparison of optical and X-ray CT technique for void analysis in soil thin section. Geoderma 141, 60–70.
| A comparison of optical and X-ray CT technique for void analysis in soil thin section.Crossref | GoogleScholarGoogle Scholar |
British Standards Institution (2018) BS EN ISO 14688–1:2018. Geotechnical investigation and testing – identification and classification of soil‚ Part 1 – identification and description.(BSI: London)
Eswaran H, Lal R, Reich PF (2001) Land degradation: an overview. In ‘Responses to land degradation. Proceedings of the Second International Conference on Land Degradation and Desertification’, Khon Kaen, Thailand. (Eds EM Bridges, ID Hannam, LR Oldeman, FWT Pening de Vries, SJ Scherr, S Sompatpanit) (Oxford Press: New Delhi, India)
FAO (2017) Voluntary guidelines for sustainable soil management. Available at http://www.fao.org/3/a-bl813e.pdf [verified 6 June 2019].
Greenland DJ (1981) Soil management and soil degradation. Journal of Soil Science 32, 301–322.
| Soil management and soil degradation.Crossref | GoogleScholarGoogle Scholar |
Grevers MCJ, de Jong E, Arnaud RJS (1989) The characterization of soil macroporosity with CT scanning. Canadian Journal of Soil Science 69, 629–637.
| The characterization of soil macroporosity with CT scanning.Crossref | GoogleScholarGoogle Scholar |
Gupta RK, Abrol I, Finkl CW, Kirkham MB, Arbestain MC, Macías F, Chesworth W, Germida JJ, Loeppert RH, Cook MG, Schwag GO, Konstankiewicz K, Pytka J, Oertli JJ, Singer A, Edmonds WJ, Feng Y, Feldman SB, Shang C, Zalazny LW, Ford PW, Clothier BE (2008) Soil pores. In ‘Encyclopedia of Soil Science’. (Ed. W Chesworth) pp. 693–698. (Springer: Dordrecht, Netherlands)
Horel Á, Tóth E, Gelybó G, Kása I, Bakacsi Z, Farkas C (2015) Effects of land use and management on soil hydraulic properties. Open Geosciences 7, 1442–1454.
| Effects of land use and management on soil hydraulic properties.Crossref | GoogleScholarGoogle Scholar |
Horn R, Domzzał H, Słowinská-Jurkiewicz A, van Ouwerkerk C (1995) Soil compaction processes and their effects on the structure of arable soils and the environment. Soil & Tillage Research 35, 23–36.
| Soil compaction processes and their effects on the structure of arable soils and the environment.Crossref | GoogleScholarGoogle Scholar |
Kargas G, Kerkides P, Sotirakoglou K, Poulovassilis A (2016) Temporal variability of surface soil hydraulic properties under various tillage systems. Soil & Tillage Research 158, 22–31.
| Temporal variability of surface soil hydraulic properties under various tillage systems.Crossref | GoogleScholarGoogle Scholar |
Kastanek FJ, Nielsen DR (2001) Description of soil water characteristics using cubic spline interpolation. Soil Science Society of America Journal 65, 279–283.
| Description of soil water characteristics using cubic spline interpolation.Crossref | GoogleScholarGoogle Scholar |
Kodesova R, Jirku V, Kodes V, Muhlhanselova M, Nikodem A, Žigová A (2011) Soil structure and soil hydraulic properties of Haplic Luvisol used as arable land and grassland. Soil & Tillage Research 111, 154–161.
| Soil structure and soil hydraulic properties of Haplic Luvisol used as arable land and grassland.Crossref | GoogleScholarGoogle Scholar |
Kubiena W (1938) ‘Micropedology.’ (Collegiate Press: Ames, IA)
Kutilek M (2004) Soil hydraulic properties as related to soil structure. Soil & Tillage Research 79, 175–184.
| Soil hydraulic properties as related to soil structure.Crossref | GoogleScholarGoogle Scholar |
Lal R (2000) Physical management of soils of the tropics: priorities for the 21st century. Soil Science 165, 191–207.
| Physical management of soils of the tropics: priorities for the 21st century.Crossref | GoogleScholarGoogle Scholar |
Lal R, Shukla MK (2004) ‘Principles of Soil Physics.’ (Marcel Dekker: New York, NY)
Lozano LA, Soracco CG, Villarreal R, Ressia JM, Sarli GO, Filgueira RR (2016) Soil physical quality and soybean yield as affected by chiseling and subsoiling of a no-till soil. Revista Brasileira de Ciência do Solo 40, 1–12.
| Soil physical quality and soybean yield as affected by chiseling and subsoiling of a no-till soil.Crossref | GoogleScholarGoogle Scholar |
Murphy BW (2015) Impact of soil organic matter on soil properties – a review with emphasis on Australian soils. Soil Research 53, 605–635.
| Impact of soil organic matter on soil properties – a review with emphasis on Australian soils.Crossref | GoogleScholarGoogle Scholar |
Othmer H, Diekkrüger B, Kutilek M (1991) Bimodal porosity and unsaturated hydraulic conductivity. Soil Science 152, 139–150.
| Bimodal porosity and unsaturated hydraulic conductivity.Crossref | GoogleScholarGoogle Scholar |
Pagliai M, Vignozzi N, Pellegrini S (2004) Soil structure and the effect of management practices. Soil & Tillage Research 79, 131–143.
| Soil structure and the effect of management practices.Crossref | GoogleScholarGoogle Scholar |
Peña-Sancho C, López M, Gracia R, Moret-Fernández D (2017) Effects of tillage on the soil water retention curve during a fallow period of a semiarid dryland. Soil Research 55, 114–123.
| Effects of tillage on the soil water retention curve during a fallow period of a semiarid dryland.Crossref | GoogleScholarGoogle Scholar |
Pires LF, Borges JA, Rosa JA, Cooper M, Heck RJ, Passoni S, Roque WL (2017) Soil structure changes induced by tillage systems. Soil & Tillage Research 165, 66–79.
| Soil structure changes induced by tillage systems.Crossref | GoogleScholarGoogle Scholar |
R Core Team (2016) R: A language and environment for statistical computing. (R Foundation for Statistical Computing: Vienna, Austria.) Available at https://www.r-project.org/ [verified 13 June 2019].
Rabot E, Wiesmeier M, Schlüter S, Vogel HJ (2018) Soil structure as an indicator of soil functions: a review. Geoderma 314, 122–137.
| Soil structure as an indicator of soil functions: a review.Crossref | GoogleScholarGoogle Scholar |
Reynolds WD (2017) Use of bimodal hydraulic property relationships to characterize soil physical quality. Geoderma 294, 38–49.
| Use of bimodal hydraulic property relationships to characterize soil physical quality.Crossref | GoogleScholarGoogle Scholar |
Reynolds WD, Elrick DE (2002) Falling head soil core (tank) method. In ‘Methods of Soil Analysis, Part 4, Physical Methods’. (Eds JH Dane, GC Topp) pp. 809–812 (Soil Science Society of America: Madison, WI)
Reynolds WD, Drury CF, Tan CS, Fox CA, Yang XM (2009) Use of indicators and pore volume–function characteristics to quantify soil physical quality. Geoderma 152, 252–263.
| Use of indicators and pore volume–function characteristics to quantify soil physical quality.Crossref | GoogleScholarGoogle Scholar |
Romano N, Nasta P (2016) How effective is bimodal soil hydraulic characterization? Functional evaluations for predictions of soil water balance. European Journal of Soil Science 67, 523–535.
| How effective is bimodal soil hydraulic characterization? Functional evaluations for predictions of soil water balance.Crossref | GoogleScholarGoogle Scholar |
Schindler U, Durner W, von Unold G, Mueller L, Wieland R (2010) The evaporation method: extending the measurement range of soil hydraulic properties using the air-entry pressure of the ceramic cup. Journal of Plant Nutrition and Soil Science 173, 563–572.
| The evaporation method: extending the measurement range of soil hydraulic properties using the air-entry pressure of the ceramic cup.Crossref | GoogleScholarGoogle Scholar |
Schjønning P, Lamandé M, Crétin V, Nielsen JA (2017) Upper subsoil pore characteristics and functions as affected by field traffic and freeze–thaw and dry–wet treatments. Soil Research 55, 234–244.
| Upper subsoil pore characteristics and functions as affected by field traffic and freeze–thaw and dry–wet treatments.Crossref | GoogleScholarGoogle Scholar |
Schlüter S, Großmann C, Diel J, Wu GM, Tischer S, Deubel A, Rücknagel J (2018) Long-term effects of conventional and reduced tillage on soil structure, soil ecological and soil hydraulic properties. Geoderma 332, 10–19.
| Long-term effects of conventional and reduced tillage on soil structure, soil ecological and soil hydraulic properties.Crossref | GoogleScholarGoogle Scholar |
Schwärzel K, Punzel J (2007) Hood infiltrometer – a new type of tension infiltrometer. Soil Science Society of America Journal 71, 1438–1447.
| Hood infiltrometer – a new type of tension infiltrometer.Crossref | GoogleScholarGoogle Scholar |
Schwen A, Bodner G, Scholl P, Buchan GD, Loiskandl W (2011) Temporal dynamics of soil hydraulic properties and the water-conducting porosity under different tillage. Soil & Tillage Research 113, 89–98.
| Temporal dynamics of soil hydraulic properties and the water-conducting porosity under different tillage.Crossref | GoogleScholarGoogle Scholar |
Soracco CG, Lozano LA, Villarreal R, Melani E, Sarli GO (2018) Temporal variation of soil physical quality under conventional and no-till systems. Revista Brasileira de Ciência do Solo 42, e0170408
| Temporal variation of soil physical quality under conventional and no-till systems.Crossref | GoogleScholarGoogle Scholar |
Tebrügge F, Düring R-A (1999) Reducing tillage intensity – a review of results from a long-term study in Germany. Soil & Tillage Research 53, 15–28.
| Reducing tillage intensity – a review of results from a long-term study in Germany.Crossref | GoogleScholarGoogle Scholar |
Tuller M, Or D (2004) Retention of water in soil and the soil water characteristic curve. In ‘Encyclopedia of Soils in the Environment’. (Ed. D Hillel) pp. 278–289. (Elsevier: Amsterdam, Netherlands)
Wairiu M, Lal R (2006) Tillage and land use effects on soil microporosity in Ohio, USA and Kolombangara, Solomon Islands. Soil & Tillage Research 88, 80–84.
| Tillage and land use effects on soil microporosity in Ohio, USA and Kolombangara, Solomon Islands.Crossref | GoogleScholarGoogle Scholar |
Weninger T, Bodner G, Kreiselmeier J, Chandrasekhar P, Julich S, Feger K-H, Schwärzel K, Schwen A (2018) Combination of measurement methods for a wide-range description of hydraulic soil properties. Water (Basel) 10, 1021
| Combination of measurement methods for a wide-range description of hydraulic soil properties.Crossref | GoogleScholarGoogle Scholar |