Next Article in Journal
Water Budget, Biological Water Use, and the Soil Hydrological Cycle across Typical Ecosystems of the Heihe River Basin
Previous Article in Journal
Hydraulic Transient Simulation of Pipeline-Open Channel Coupling Systems and Its Applications in Hydropower Stations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soil-Gas Diffusivity-Based Characterization of Variably Saturated Agricultural Topsoils

1
Department of Civil Engineering, Faculty of Engineering, University of Peradeniya, Peradeniya 20400, Sri Lanka
2
College of Pastoral Agriculture Science and Technology, National Field Scientific Observation and Research Station of Grassland Agro-Ecosystems in Gansu Qingyang, The State Key Laboratory of Grassland Agro-Ecosystems of Lanzhou University, Lanzhou 730020, China
3
State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Northwest A&F University, Xianyang 712100, China
4
Department of Soil and Physical Sciences, Lincoln University, P.O. Box 85084, Lincoln 7647, New Zealand
5
Department of Geosciences and Natural Resource Management, Center for Permafrost (CENPERM), University of Copenhagen, Øster Voldgade 10, DK-1350 Copenhagen, Denmark
6
Department of Civil and Environmental Engineering, Southern Methodist University, Dallas, TX 75205, USA
*
Author to whom correspondence should be addressed.
Water 2022, 14(18), 2900; https://doi.org/10.3390/w14182900
Submission received: 26 August 2022 / Revised: 8 September 2022 / Accepted: 14 September 2022 / Published: 16 September 2022
(This article belongs to the Section Soil and Water)

Abstract

:
Soil-gas diffusivity and its variation with soil moisture plays a fundamental role in diffusion-controlled migration of climate-impact gases from different terrestrial agroecosystems including cultivated soils and managed pasture systems. The wide contrast in soil texture and structure (e.g., density, soil aggregation) in agriculture topsoils (0–10 cm) makes it challenging for soil-gas diffusivity predictive models to make accurate predictions across different moisture conditions. This study characterized gas diffusivity and gas-phase tortuosity in soils sampled from managed pasture and cultivated sites in Sri Lanka at 0–10 cm depth, together with selected soil-gas diffusivity data from the literature. Soil-gas diffusivity was measured using a one-chamber diffusion apparatus using N2 and O2 as experimental gases. The measured diffusivity, together with literature data representing both intact and repacked soils, were tested against five existing widely known gas diffusivity predictive models. The tested models tended to mischaracterize the two-region behavior in some of the aggregated soils, suggesting the need of soil-specific diffusivity models to better describe gas diffusivity in agricultural soils. We suggested a new parametric two-region model, developed in line with literature-based models, to represent both unimodal and bimodal/two-region behavior of selected soils. The new model statistically outperformed the existing predictive models for both intact and repacked soils and, hence, demonstrated its applicability to better characterize site-specific greenhouse gas emissions under different soil water regimes.

1. Introduction

Terrestrial agroecosystems are predominant anthropogenic sources of methane (CH4) and nitrous oxide (N2O), two high-potency greenhouse gases which entail inevitable climate footprints such as global warming and region climate shifts. Both CH4 and N2O are strong greenhouse gases having a global warming potential 25-fold and 298-fold higher than that of CO2, respectively, over a 100-year time horizon [1]. Among different cultivated systems, paddy ecosystems contribute to approximately 20% of human-induced CH4 emissions, due mainly to the typical submerged conditions practiced in paddy cultivation which favor methanogenesis in the paddy rhizosphere [2]. Essentially, intrinsic soil properties such as soil organic carbon, total nitrogen, organic matter, texture, etc. significantly affect the soil microbial functions and, thereby, can alter gas emissions to a great extent [3], while the diversity of the soil microbial community also plays a crucial role in gas production [4]. In addition, extensive application of nitrogen (N) fertilizers in paddy landscapes favors enhanced denitrification-driven N2O emissions. Though paddy emissions occur mainly by ebullition under submerged conditions, paddy fields typically experience seasonal water fluctuations which enable the produced gases to migrate and emit largely through the soil. In fact, alternate wetting and draining (AWD) is an emerging water-saving strategy practiced in countries for paddy cultivation in which the mid-season draining causes the release of produced greenhouse gases, largely into the atmosphere [5,6,7].
Grazed pasture constitutes another prominent agricultural system which is responsible for nearly 28% of anthropogenic N2O in the atmosphere [8]. Pastoral N2O is primarily derived from applied N fertilizers and also from the manure of the grazing livestock and is produced via nitrification and denitrification in the presence of anaerobic conditions in pasture topsoil. Scattered patches of ruminant urine and feces are prominent sources of N that potentially create N2O hotspots in pasture topsoils. The changes in water regimes in the pasture subsurface upon irrigation or rainfall make a major impact on fate and transport of N2O in the subsurface.
Once produced, the gases must be transported through the subsurface via the topsoil and emit across the soil–atmosphere continuum in order to make an impactful climate effect. If the topsoil is sufficiently aerated, N2O will be oxidized and emitted into atmosphere as dinitrogen (N2), an environmentally benign gas as compared to N2O, with no climate forcing effect. Similarly, the draining of paddy soils may aerate the topsoils and, therefore, result in the emission of CO2 in lieu of CH4. Thus, the agricultural topsoils (0–10 cm) become critically important in regulating the climate in local, regional, and global contexts.
Migration of greenhouse gases in agricultural topsoils and their emission into the atmosphere are predominantly diffusion-controlled, particularly in the absence of wind-induced pressure gradients. In addition, soil physical properties such as soil texture (soil type) and structure (e.g., soil density, aggregation) play undeniable roles in agricultural CH4 and N2O emissions. In particular, soil aggregation and associated hierarchical structural development in agricultural topsoils have been frequently discussed in the literature [9,10,11,12]. Extensive application of agricultural machinery, high livestock stocking rates, and frequent tillage of agricultural soils are among major controls affecting soil functional structure in pasture and cultivated soils. Mitigation measures, therefore, are highly dependent on how well the diffusion-controlled gas transport processes in agricultural topsoils are understood and how accurately they can be accounted for in predictive numerical tools.
Gas diffusion in soil is commonly characterized by soil-gas diffusivity, Dp/Do, where Dp (m3 soil air m−1 soil s−1) and Do (m2 air s−1) are the soil-gas diffusion coefficients in soil and in free air, respectively. Measuring Dp/Do is, however, experimentally intensive and instrumentally challenging due to the requirements of specific apparatus and tight control of initial and boundary conditions. The predictive gas diffusivity models offer a convenient alternative to predict soil-gas diffusivity from easily measurable properties such as air-filled porosity (ε) and total porosity (Φ). In complex geomedia such as well-structured agricultural soils, however, the soil physical heterogeneity makes the application of developed models questionable, since widely used models may have not been validated against soils with distinct structural contrast. Therefore, a careful revisiting of existing predictive models, together with measured gas diffusivity data on a wide range of agricultural topsoils, is an important prerequisite in investigating their applicability and to make reliable emission estimates.
Generally, agricultural topsoils are considered to be well-structured aggregated soils, having both inter-aggregate pores (i.e., pores between the aggregates) and intra-aggregate pores (i.e., pores within the aggregates), resulting in a distinctive bimodal pore structure. The two pore regions are generally considered to be functionally analogous with respect to soil-moisture and soil-gas dynamics and, hence, it is common to describe them in terms of two additive mathematical expressions in modeling soil-moisture retention [13] as well as soil-gas diffusivity [9]. However, compaction-induced changes due to animal treading, mechanical implements, and tilling operations on agricultural fields often alter soil pore structure [14]. Compaction essentially reduces the macropore domains and increases the micropore domains of a soil, thus shifting the bimodal nature of the structured soil. Although the predictive models developed for soil-gas diffusivity in unimodal soils are abound in the literature, these models cannot be used directly for aggregated soils as they tend to mischaracterize the two distinct pore regions. Therefore, several models have been modified, extended, and reintroduced to predict the soil-gas diffusivity in well-structured aggregated soils [14,15] with additional model parameters to numerically differentiate the different pore domains. Consequently, most two-region models are descriptive in nature, and use few a priori measurements to inversely estimate the additional model parameters. The classical numerical models also account for the two structural status of measured soil samples, intact or structurally disturbed, which also require close attention when agricultural soils are investigated.
The main objective of this study is to characterize gas transport behavior in cultivated and pasture topsoils (0–10 cm), including both structurally intact and disturbed soils, and also to investigate the applicability of widely known soil-gas diffusivity predictive models against these structurally diverse topsoils. A series of diffusivity measurements on selected cultivated and pasture soils as well as literature data were used to characterize the cultivated and pasture topsoils (0–10 cm). We used measured data from Sri Lankan cultivated and pasture sites, together with additional supporting data from the literature, to represent a wide geographic origin including soils from Japan (Nishi-Tokyo) and the United Kingdom (Lexington). Five widely applied soil-gas diffusivity models were tested against diffusivity data, along with a multi-parameter two-region descriptive gas diffusivity model to better characterize the two-region behavior of topsoils.

2. Materials and Methods

2.1. Soils and Data

In total, 12 soils were considered, with 6 soils representing each intact and repacked soil group. Each group included three pasture soils and three cultivated soils. Undisturbed soils were sampled from Ambewela pasture site (6.8693° N, 80.7957° E), Sri Lanka, at a 0–10 cm depth from five locations, along a downgradient transect, which also showed a natural organic matter gradient. Surface cracks, visible to the naked eye, could be seen on the topsoil of the sampling site and retrieved samples, likely due to the extensive presence of agricultural livestock. In addition, three Sri Lankan cultivated (paddy) soils, two intact soils from Kurunegala (7.5255° N, 80.4390° E) and Polonnaruwa (7.9400° N, 81.0174° E), and one repacked soil from Kandy (7.2794° N, 80.5910° E) were also used in this study. All paddy soils had been mechanically tilled before plantation, and samples were taken from a 0–10 cm depth.
In addition to the unpublished data mentioned above, literature data on pasture and cultivated soils, both in intact and disturbed status, were used in this study.
Table 1 shows all 12 considered soils from Sri Lanka (Ambewela, Kurunegala, Polonnaruwa, Kandy, and Peradeniya), Japan (Nishi-Tokyo), and the United Kingdom (Lexington), with the names denoting the sampling location. Textural contrast and important soil physical properties of selected intact and repacked soils are also given in Table 1. For further details on the data from the literature, the reader is referred to the related literature mentioned in Table 1.

2.2. Methods

The 100 cm3 undisturbed soil samples (height = 3.95 cm, diameter = 5.70 cm) were collected from Sri Lankan sites using metallic annular cores, ensuring minimum disturbance to the soil. The retrieved samples, leveled and kneaded at both ends to remove redundant soil, were end-capped and wrapped with polythene to preserve ambient moisture before being transferred to the Geotechnical Engineering laboratory for characterization. Repacked samples, on the other hand, were prepared from bulk soil sampled from a 0–10 cm depth, sieved to obtain 2 mm fraction, and repacked to the same bulk density observed in the corresponding layer.
The samples were first characterized for soil-moisture content, particle density, bulk density, particle size distribution (sieve analysis and hydrometer method), and organic matter content (loss on ignition method). The samples were then saturated for 72 h and sequentially drained to the intended moisture levels by stepwise evaporation. At each moisture level, the samples were kept closed for a sufficient amount of time to reach the hydraulic equilibrium and moisture redistribution before diffusivity measurements.
Soil-gas diffusivity was measured following the one-chamber method introduced by Taylor [18] and developed by Schjønning [19]. The custom-fabricated PVC chamber is 20 cm in height and 3.4 cm in internal diameter and is provisioned with two values as inlet and outlet for priming. The chamber was checked for airtightness before the measurement campaign was initiated. For the measurements, the sample was mounted on top of the diffusion chamber and made it airtight. Then, the chamber was flushed with 99.99% N2 gas to remove all O2 inside the chamber. The sample was then opened to the atmosphere by allowing the atmospheric O2 to diffuse through the sample into the chamber. The increase in O2 concentration inside the chamber was monitored continually with an O2 sensor attached to the chamber wall. Calculation of Dp/Do was performed following both Taylor (1949) [18] and Currie (1960) [20] methods as follows:
The Taylor (1949) method is founded on Fick’s first law and can be presented by the following equation:
ln Δ C t Δ C 0 = D p H s H c t  
where Δ C t is the change in concentration inside the chamber (gm−3) ( Δ C t = C t = t C t = 0 ),   H s is the sample height (m), and H c is the height of the chamber (m). Dp can be calculated from the gradient of the graph of ln ( Δ Ct/ Δ C0) vs. time (t).
The Currie (1960) method, on the other hand, is founded on both Fick’s first law and the second law and the basic equation can be expressed as follows:
ln Δ C t Δ C 0 = D p α 1 2 ε t + ln 2 h H s α 1 2 + h 2 + h
where Δ C t is the change in concentration inside the chamber (gm−3) ( Δ C t =   C t = t   C t = 0 ), H s is the sample height (m), h is ε / H c , ε is the air-filled porosity (m3 m−3), and H c is the height of the chamber (m). Dp (soil-gas diffusion) can be derived from the slope of the plot of ln ( Δ C t Δ C 0 ) versus time (t).

2.3. Soil-Gas Diffusivity Modeling

2.3.1. Unimodal Dp/Do Predictive Models

Table 2 shows an array of unimodal (one-region) predictive models for soil-gas diffusivity that were tested against the measured diffusivity data mentioned in Table 1.

2.3.2. Bimodal (Two-Region) Dp/Do Models

The two-region soil-gas diffusivity models are sparse in the literature, and typically consist of two mathematical expressions describing Dp/Do in two regions. For example, the recent two-region model from Jayaratne et al. [12] takes the form of:
For Region 1:
D p D o = α 1 w β 1 ε Φ β 1
For Region 2:
D p D o = α 1 + α 2 ( 1 w ) β 2 ε w Φ Φ β 2
where α1, β1, and α2, β2 (dimensionless) are parameters characterizing the Region 1 (interaggregate) and Region 2 (intraaggregate), respectively, while w is a weighing factor that numerically demarcates Region 1 and Region 2. In an extended formulation, Shanujah et al. [21] presented a two-region model incorporating an additional parameter εp (cm3 cm−3), referred to as a gas percolation threshold, an air-filled porosity below which the gas diffusivity is negligibly small due to extremely high moisture-induced tortuosity in the soil-gas phase in Region 1. Thus, the model for Region 1 takes the form:
For Region 1:
D p D o = α 1 w ε p β 1 ε ε p Φ β 1
while for the Region 2, it is the same as in Equation (4).
In order to make the formulation more physically meaningful in the two regions, we present here a slightly modified version of the above two models as follows:
For Region 1:
D p D o = α 1 ε ω o Φ ω Φ ω o Φ β 1
For Region 2:
D p D o = α 1 + α 2 ε ω Φ Φ ω Φ
where wo (dimensionless) is the fraction of total porosity that corresponds to the gas percolation threshold; other symbols denote the same as above. The model implicitly assumes that Dp/Do in the intraaggregate pore region exhibits (or can be approximated to) a linear increase with air-filled porosity as assumed in most of the previously developed two-region models. Note that (ε − ωoΦ)/(ωΦ − ωoΦ) in Equation (6) represents the effective air saturation at a given air-filled porosity (ε) within Region 1, whereas (ε − ωΦ)/(Φ − ωΦ) in Equation (7) represents the same at a given air-filled porosity within Region 2. Despite the different formulation, the new model (Equations (6) and (7)) is intended to perform at the same level as the model from the two previous models (Equations (3)–(5)); therefore, we only tested the new model performance in statistical model tests, as presented below.
The gas-phase tortuosity (τ, dimensionless) can be described as the roundabout distance a gas molecule will travel in the gaseous phase between two predefined points which are a unit Euclidean distance apart. Based on this definition, τ can be deduced from soil-gas diffusivity as presented by Ball [22] as follows:
τ = ε D p D o

2.4. Statistical Analysis

Two statistical indices, RMSE and bias, were used to statistically analyze model performance and compare model predictions. The RMSE evaluates the overall fit of a model to the measured data.
RMSE = 1 n i = 1 n d i 2
The bias, on the other hand, evaluates whether a model overestimates (positive bias) or underestimates (negative bias) the observations.
Bias = 1 n i = 1 n d i
where di is the difference between the observed and predicted diffusivity values and n is the number of diffusivity measurements in a data set.
The performance of the selected models and the two-region model against the measured gas diffusivity data for intact soils and repacked soils expressed in terms of RMSE and bias is shown in Table 3 and Table 4, respectively.

3. Results and Discussion

3.1. Soil-Gas Diffusivity

Figure 1 shows the soil-gas diffusivity against air-filled porosity for six intact soils from Sri Lanka and the UK, representing pasture (Figure 1a–c) and cultivated (Figure 1d–f) soils, together with the two-region descriptive model (Equations (6) and (7)).
In selected topsoils, the aggregation seems to have been disrupted due to the frequent treading of livestock ruminants and application of agricultural machinery, leaving a non-aggregated structure, as we see for most of the soils in Figure 1. A distinct exception to this is Ambewela pasture soil, which showed a linear behavior in the Dp/Do vs. ε relation, as shown in Figure 1a, likely due to the cracked nature of the soil formed from extensive animal treading and trampling, which has resulted in an altered soil structure. In fact, large cracks could be seen on the soil surface during sampling. Similar linear variations in cracked soils/fractured geomedia in characteristic diffusivity curves have been observed in literature [23], while the gradient of the linear model provides a useful implication of the average fracture angle of the geosystem.
In the Peradeniya pasture site, on the other hand, the use of livestock or agricultural implements was not extensive over a couple of years prior to sampling and, therefore, the soil has gradually developed into a well-structured aggregated soil. Consequently, Peradeniya soil has exhibited a distinct two-region behavior (Figure 1b). Note the decline in the rate of increase in gas diffusivity when the gases enter the intraaggregate pore space (region 2), likely due to the encounter of a more tortuous pore network within the aggregates. In fact, in strongly aggregated porous media (e.g., manufactured porous substrates), a distinct threshold in gas diffusivity has been observed in the literature [24] due to the presence of highly tortuous pore network inside aggregates. The UK pasture soil also showed a slight two-region behavior, owing to its weakly aggregated nature, but was not so pronounced as the Peradeniya soil. On the other hand, the three cultivated soils, Kurunegala, Peradeniya, and Polonnaruwa (Sri Lanka), which have undergone a regular mechanical tillage before cultivation, showed unaggregated nature and, hence, yielded a typical unimodal behavior in gas diffusivity.
Figure 2 shows the measured soil-gas diffusivity against air-filled porosity, together with the descriptive two-region model (Equations (6) and (7)) for six repacked soils: three pasture soils (two Peradeniya soils and the Nishi-Tokyo soil) and three cultivated soils (Peradeniya, Nishi-Tokyo, and Kandy soils). Interestingly, all six repacked soils showed a two-region behavior. Typically, since the sieving and repacking involve de-structuring of soil aggregates in repacked soils, they are likely to show fewer two-region characteristics than intact soils. Note that the repacking did not cause the crushing of aggregates; hence, the aggregated behavior of sieved soils was preserved in repacked samples. In comparison, the three cultivated soils showed a higher gas percolation threshold compared to the three pasture soils. In Peradeniya repacked soils, the less-compacted soils (1.0 g cm−3, Figure 2a) have higher interaggregate pore volume than the highly compacted soils (1.30 g cm−3, Figure 2b), implying a decrease in interaggregate pore space due to high compaction, with an overall decrease in Dp/Do.
To examine the performance of existing widely used predictive models on agricultural topsoils, we tested five selected predictive models, which comprise both classical and recent models, against the measured Dp/Do for undisturbed topsoils (Figure 3). Evidently, excluding the Millington (1959) model [25], all other models tended to markedly underpredict diffusivity in undisturbed soils. In fact, predictions were better for cultivated soils than for pasture soils due to the altered structural properties in pasture soils, as discussed above. The Millington and Quirk (1961) model [26], the widely accepted model for predicting soil-gas diffusivity in gas transport models, tended to overpredict the Dp/Do in dry soils and markedly underpredict the Dp/Do in wet soils, as reported in many previous works in the literature [17,27,28].
The WLR model, developed specifically for sieved and repacked soils, overpredicted the cultivated intact soils as expected, but underpredicted the undisturbed pasture soils due to the enhanced soil structure leading to pronounced diffusion properties. In the SWLR model, in general, a complexity factor smaller (Cm = 1.5) than originally proposed (Cm = 2.1) was required to better characterize the pasture topsoils. In Figure 3f, the behavior of the two-region descriptive model (Equations (6) and (7)) shows the degree of improved prediction that can be achieved if a two-region model is used in place of the classical unimodal models.
A similar scatterplot comparison was also made for repacked soils, as depicted in Figure 4. As expected, both WLR [27] and SWLR (with Cm = 1.0) [28] models performed better due to their intended applicability for repacked soils, and the Buckingham (1904) model [29] also provided promising results. Furthermore, as shown in Figure 4, the behavior of the two-region descriptive model (Equations (6) and (7)) demonstrates the improved predictability when optimized model parameters are used.
Overall, it is clear that the unimodal models apparently mischaracterized the agriculture topsoils with altered soil structure, leading to a significant disparity in measured and modeled diffusivity data for both intact and repacked soils. Application of a two-region descriptive model, with a few a priori diffusivity measurements and calibrated model parameters, can be used to better characterize agricultural topsoils. Such a descriptive model calibrated with measured data specifically for a given soil can be particularly useful when site-specific agricultural management practices are adopted to enhance crop productivity and minimize agricultural emissions. When access to diffusivity measurements is lacking, our analysis shows that the Buckingham model, albeit the oldest and the simplest among the classical models, can be used for reasonably good predictions for agricultural topsoils.
The parameterization of the two-region model (Equations (6) and (7)) is given in Table 5 for both intact and repacked soils. Excluding Kurunegala soil, all other soils have a percolation threshold less than 10% of the total pore space (i.e., w0 < 0.1). Numerical analysis further shows that the interaggregate porosity in bimodal soils corresponds to 43–73% of the total porosity.

3.2. Gas-Phase Tortuosity

The calculated tortuosity of the functional soil gaseous phase from the measured soil-gas diffusivity and air-filled porosity data (Equation (8)) for the intact and repacked agricultural topsoils are shown in Figure 5, together with calculated tortuosity from the two-region model (Equations (6) and (7)). Note that the tortuosity is presented in decreasing order, which denotes the increasing air-filled porosity. Clearly, the tortuosity is high in the wet region due to the pronounced moisture-induced effects on pore network with constrained gas diffusion, which gradually decrease as the soil drains and becomes minimum in dry soil. The observed minimum tortuosity values (around 1.1 for both intact and repacked soils) are generally lower than those of typical mineral soils (around 1.5–2.0) due to the enhanced pore connectivity in aggregated soils. The uncertainty in tortuosity measurements and predictions is also high in the wet region due to the random connection and disconnection of water bridges across soil particles, yielding a large scatter in tortuosity values. As expected, the scatter in tortuosity is typically higher for intact soils than for the repacked soils due to the intrinsic soil structural effects in intact soils. Notably, the two-region model demonstrated its usefulness in describing the tortuosity in agricultural topsoils.
This study has provided useful experimental and numerical insight into characterization of agricultural topsoils based on measured and modeled soil-gas diffusivity data. It is worthy to be noted that, due to the wide structural diversity of agricultural topsoils, further measurements on both intact and repacked soils will enable a more detailed characterization of soils. It should be noted that all measurements involved have been carried out in laboratory-controlled environments where natural environmental complexities (e.g., dynamics in temperature, evaporation, wind, and humidity) were eliminated. Such additional environmental factors were out of the scope of this study but must be accounted for when making more realistic conclusions. Results, therefore, need to be compared against field-measured data with caution.

4. Conclusions

This study investigated new measurements and literature data on soil-gas diffusivity to characterize agricultural topsoils (0–10 cm) representing pasture and cultivated soil systems. The diffusivity data showed a wide structural diversity in both intact and structurally disturbed samples, with some soils showing distinct two-region behavior. A two-region descriptive soil-gas diffusivity model, with model parameters calibrated against a few a priori measured data, seemed to be useful in more accurately characterizing the structure-dependent behavior the of considered agricultural topsoils when compared to the widely used predictive gas diffusivity models. The gas-phase pore tortuosity analysis showed pronounced moisture-induced pore tortuosity with a large scatter in wet soils, while the two-region model demonstrated an adequate description of the calculated tortuosity. The overall numerical analysis, thus, highlighted the importance of site-specific characterization of gas diffusivity and pore tortuosity in order to make better predictions of soil-gas emissions as a prerequisite to introduce agricultural management strategies in relation to mitigation of greenhouse gases. The laboratory-controlled measurements of this study must be compared cautiously with the results of field-based studies, which provide more realistic observations in the presence of other environmental factors.

Author Contributions

Conceptualization, validation, writing—original draft preparation, A.M.S.N.A. and M.M.T.L.; conceptualization, validation, writing—original draft preparation, U.D.H.N.A.; visualization, writing—review and editing, Y.L.; conceptualization, writing—review and editing, funding acquisition, T.K.K.C.D.; visualization, validation, W.F., J.F., T.Y. and X.M.; visualization, supervision, T.C., B.E. and K.S. All authors have read and agreed to the published version of the manuscript.

Funding

We gratefully acknowledge the financial support for this research from the Asia-Pacific Network for Global Change Research under the project reference number “CRRP2020-07MY-Deepagoda”. Any opinion, findings, and conclusions or recommendations expressed herein are those of the authors and do not necessarily reflect the views of those providing the financial support.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We gratefully acknowledge the financial support for this research from the Asia-Pacific Network for Global Change Research under the project reference number “CRRP2020-07MY-Deepagoda”. Any opinion, findings, and conclusions or recommendations expressed herein are those of the authors and do not necessarily reflect the views of those providing the financial support.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Myhre, G.; Shindell, D.; Breon, F.-M.; Collins, W.; Fuglestvedt, J.; Huang, J.; Koch, D.; Lamarque, J.-F.; Lee, D.; Mendoza, B.; et al. Anthropogenic and natural radiative forcing. In Climate change 2013: The Physical Science Basis; Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Stocker, T.F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P.M., Eds.; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  2. Qiu, J. China cuts methane emissions from rice fields. Nature 2009. [Google Scholar] [CrossRef]
  3. Arunrat, N.; Sansupa, C.; Kongsurakan, P.; Sereenonchai, S.; Hatano, R. Soil Microbial Diversity and Community Composition in Rice–Fish Co-Culture and Rice Monoculture Farming System. Biology 2022, 11, 1242. [Google Scholar] [CrossRef] [PubMed]
  4. Muhammad, I.; Lv, J.Z.; Wang, J.; Ahmad, S.; Farooq, S.; Ali, S.; Zhou, X.B. Regulation of Soil Microbial Community Structure and Biomass to Mitigate Soil Greenhouse Gas Emission. Front. Microbiol. 2022, 13, 868862. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, X.; Zhou, T.; Liu, Y.; Zhang, X.; Li, L.; Pan, G. Effect of mid-season drainage on CH4 and N2O emission and grain yield in rice ecosystem: A meta-analysis. Agric. Water Manag. 2018, 213, 1028–1035. [Google Scholar] [CrossRef]
  6. Liao, B.; Wu, X.; Yu, Y.; Luo, S.; Hu, R.; Lu, G. Effects of mild alternate wetting and drying irrigation and mid-season drainage on CH4 and N2O emissions in rice cultivation. Sci. Total Environ. 2019, 698, 134212. [Google Scholar] [CrossRef]
  7. Toma, Y.; Takechi, Y.; Inoue, A.; Nakaya, N.; Hosoya, K.; Yamashita, Y.; Adachi, M.; Kono, T.; Hideto, U. Early mid-season drainage can mitigate greenhouse gas emission from organic rice farming with green manure application. Soil Sci. Plant Nutr. 2021, 67, 482–492. [Google Scholar] [CrossRef]
  8. Abdalla, M.; Jones, M.; Smith, P.; Williams, M. Nitrous oxide fluxes and denitrification sensitivity to temperature in Irish pasture soils. Soil Use Manag. 2009, 25, 376–388. [Google Scholar] [CrossRef]
  9. Resurreccion, A.C.; Moldrup, P.; Kawamoto, K.; Yoshikawa, S.; Rolston, D.E.; Komatsu, T. Variable pore connectivity factor model for gas diffusivity in unsaturated, aggregated soil. Vadose Zone J. 2008, 7, 397–405. [Google Scholar] [CrossRef]
  10. Resurreccion, A.C.; Moldrup, P.; Kawamoto, K.; Hamamoto, S.; Rolston, D.E.; Komatsu, T. Hierarchical, bimodal model for gas diffusivity in aggregated, unsaturated soils. Soil Sci. Soc. Am. J. 2010, 74, 481–491. [Google Scholar] [CrossRef]
  11. Ghanbarian, B.; Hunt, A.G. Universal scaling of gas diffusion in porous media. Water Resour. Res. 2014, 50, 2242–2256. [Google Scholar] [CrossRef]
  12. Jayarathne, J.R.R.N.; Chamindu Deepagoda, T.K.K.; Clough, T.J.; Nasvi, M.; Thomas, S.; Elberling, B.; Smits, K.; Kankanamge, C.D.T. Gas-Diffusivity based characterization of aggregated agricultural soils. Soil Sci. Soc. Am. J. 2020, 84, 387–398. [Google Scholar] [CrossRef]
  13. Durner, W. Hydraulic conductivity estimation for soils with heterogeneous pore structure. Water Resour. Res. 1994, 30, 211–223. [Google Scholar] [CrossRef]
  14. Jayarathne, J.; Chamindu Deepagoda, T.K.K.; Clough, T.J.; Thomas, S.; Elberling, B.; Smits, K.M. Effect of aggregate size distribution on soil moisture, soil-gas diffusivity, and N2O emissions from a pasture soil. Geoderma 2020, 383, 114737. [Google Scholar] [CrossRef]
  15. Chamindu Deepagoda, T.K.K.; Smits, K.M.; Ramirez, J.; Moldrup, P. Characterization of thermal, hydraulic, and gas diffusion properties in variably saturated sand grades. Vadose Zone J. 2016, 15. [Google Scholar] [CrossRef]
  16. Kreba, S.; Coyne, M.S.; McCulley, R.L.; Wendroth, O.O. Spatial and Temporal Patterns of Carbon Dioxide Flux in Crop and Grass Land-Use Systems. Vadose Zone J. 2013, 12, 1–16. [Google Scholar] [CrossRef]
  17. Chamindu Deepagoda, T.K.K.; Moldrup, P.; Schjønning, P.; de Jonge, L.W.; Kawamoto, K.; Komatsu, T. Density-Corrected Models for Gas Diffusivity and Air Permeability in Unsaturated Soil. Vadose Zone J. 2011, 10, 226–238. [Google Scholar] [CrossRef]
  18. Taylor, S.A. Oxygen Diffusion in Porous Media as a Measure of Soil Aeration. Soil Sci. Soc. Am. J. 1950, 14, 55–61. [Google Scholar] [CrossRef]
  19. Schjønning, P. A Laboratory Method for Determination of Gas Diffusion in Soil; Report S1773; Danish Institute of Plant and Soil Sci Science: Tjele, Denmark, 1985; (In Danish with English Summary). [Google Scholar]
  20. Currie, J.A. Gaseous diffusion in porous media. Part 2. Dry granular materials. Br. J. Appl. Phys. 1960, 11, 318–324. [Google Scholar] [CrossRef]
  21. Shanujah, M.; Deepagoda, C.; Smits, K.M.; Shreedharan, V.; Paramrswaran, T.G.; Babu, G.L.S. Gas diffusivity-based characterization of aggregated soils linking to methane migration in shallow subsurface. Vadose Zone J. 2021, 20, e20135. [Google Scholar] [CrossRef]
  22. Ball, B.C. Modelling of soil pores as tubes using gas permeabilities, gas diffusivities and water release. Eur. J. Soil Sci. 1981, 32, 465–481. [Google Scholar] [CrossRef]
  23. Kristensen, A.H.; Thorbjørn, A.; Jensen, M.P.; Pedersen, M.; Moldrup, P. Gas-phase diffusivity and tortuosity of structured soils. J. Contam. Hydrol. 2010, 115, 26–33. [Google Scholar] [CrossRef] [PubMed]
  24. Deepagoda, T.C.; Lopez, J.C.C.; Møldrup, P.; De Jonge, L.W.; Tuller, M. Integral parameters for characterizing water, energy, and aeration properties of soilless plant growth media. J. Hydrol. 2013, 502, 120–127. [Google Scholar] [CrossRef]
  25. Millington, R.J. Gas Diffusion in Porous Media. Science 1959, 130, 100–102. [Google Scholar] [CrossRef] [PubMed]
  26. Millington, R.J.; Quirk, J.P. Permeability of porous solids. Trans. Faraday Soc. 1961, 57, 1200–1207. [Google Scholar] [CrossRef]
  27. Moldrup, P.; Olesen, T.; Gamst, J.; Schjønning, P.; Yamaguchi, T.; Rolston, D.E. Predicting the gas diffusion coefficient in repacked soil: Water induced linear reduction model. Soil Sci. Soc. Am. J. 2000, 64, 158. [Google Scholar] [CrossRef]
  28. Moldrup, P.; Deepagoda, T.C.; Hamamoto, S.; Komatsu, T.; Kawamoto, K.; Rolston, D.E.; de Jonge, L.W. Structure-Dependent Water-Induced Linear Reduction Model for Predicting Gas Diffusivity and Tortuosity in Repacked and Intact Soil. Vadose Zone J. 2013, 12, 1–11. [Google Scholar] [CrossRef]
  29. Buckingham, E. Contributions to Our Knowledge of the Aeration of Soils; U.S. Department of Agriculture, Bureau of Soils: Washington, DC, USA, 1904.
Figure 1. Measured soil-gas diffusivity against air-filled porosity for six intact soils (af), together with the descriptive two-region model (Equations (6) and (7)). AMB; Ambewela (Sri Lanka), PDY; Peradeniya (Sri Lanka), LXT; Lexington (United Kingdom), KRN; Kurunegala (Sri Lanka); PLN; Polonnaruwa (Sri Lanka). P and C in parentheses denote Pasture and Cultivated, respectively.
Figure 1. Measured soil-gas diffusivity against air-filled porosity for six intact soils (af), together with the descriptive two-region model (Equations (6) and (7)). AMB; Ambewela (Sri Lanka), PDY; Peradeniya (Sri Lanka), LXT; Lexington (United Kingdom), KRN; Kurunegala (Sri Lanka); PLN; Polonnaruwa (Sri Lanka). P and C in parentheses denote Pasture and Cultivated, respectively.
Water 14 02900 g001
Figure 2. Measured soil-gas diffusivity against air-filled porosity for six repacked soils ((af), together with proposed descriptive two-region model (Equations (6) and (7)). PDY; Peradeniya (Sri Lanka), NTK; Nishi-Tokyo (Japan), KDY; Kandy (Sri Lanka). P and C in parentheses denote Pasture and Cultivated, respectively.
Figure 2. Measured soil-gas diffusivity against air-filled porosity for six repacked soils ((af), together with proposed descriptive two-region model (Equations (6) and (7)). PDY; Peradeniya (Sri Lanka), NTK; Nishi-Tokyo (Japan), KDY; Kandy (Sri Lanka). P and C in parentheses denote Pasture and Cultivated, respectively.
Water 14 02900 g002
Figure 3. Scatterplot showing the measured and modeled Dp/Do for intact soils using six predictive-descriptive models (af). AMB; Ambewela (Sri Lanka), PDY; Peradeniya (Sri Lanka), LXT; Lexington (United Kingdom), KRN; Kurunegala (Sri Lanka); PLN; Polonnaruwa (Sri Lanka). P and C in parentheses denote Pasture and Cultivated, respectively.
Figure 3. Scatterplot showing the measured and modeled Dp/Do for intact soils using six predictive-descriptive models (af). AMB; Ambewela (Sri Lanka), PDY; Peradeniya (Sri Lanka), LXT; Lexington (United Kingdom), KRN; Kurunegala (Sri Lanka); PLN; Polonnaruwa (Sri Lanka). P and C in parentheses denote Pasture and Cultivated, respectively.
Water 14 02900 g003
Figure 4. Scatterplot showing the measured and modeled Dp/Do for repacked soils, using six predictive-descriptive models (af). PDY; Peradeniya (Sri Lanka), NTK; Nishi-Tokyo (Japan), KDY; Kandy (Sri Lanka). P and C in parentheses denote Pasture and Cultivated, respectively.
Figure 4. Scatterplot showing the measured and modeled Dp/Do for repacked soils, using six predictive-descriptive models (af). PDY; Peradeniya (Sri Lanka), NTK; Nishi-Tokyo (Japan), KDY; Kandy (Sri Lanka). P and C in parentheses denote Pasture and Cultivated, respectively.
Water 14 02900 g004
Figure 5. Pore tortuosity analysis using the two-region descriptive model (Equations (6) and (7)) for (a) intact and (b) repacked soils. AMB; Ambewela (Sri Lanka), PDY; Peradeniya (Sri Lanka), LXT; Lexington (United Kingdom), KRN; Kurunegala (Sri Lanka); PLN; Polonnaruwa (Sri Lanka), NTK; Nishi-Tokyo (Japan), KDY; Kandy (Sri Lanka). P and C in parentheses denote Pasture and Cultivated, respectively.
Figure 5. Pore tortuosity analysis using the two-region descriptive model (Equations (6) and (7)) for (a) intact and (b) repacked soils. AMB; Ambewela (Sri Lanka), PDY; Peradeniya (Sri Lanka), LXT; Lexington (United Kingdom), KRN; Kurunegala (Sri Lanka); PLN; Polonnaruwa (Sri Lanka), NTK; Nishi-Tokyo (Japan), KDY; Kandy (Sri Lanka). P and C in parentheses denote Pasture and Cultivated, respectively.
Water 14 02900 g005
Table 1. Soils and data from the literature and their physical properties.
Table 1. Soils and data from the literature and their physical properties.
SoilVegetationSoil TypeSoil TextureBulk
Density
Organic MatterTotal
Porosity
Reference
Sand
(%)
Silt
(%)
Clay
(%)
Mg m−3(kg kg−1)(cm3 cm−3)
Intact Soils
Ambewela, SLpastureSandy loam54.139.56.30.72 (0.02)0.070.70This study
Peradeniya, SLpastureSandy loam72.125.12.81.33 (0.05)0.0960.46[14]
Lexington, UKpastureSilt loam7.367.325.41.51(0.06)0.050.43[16]
Peradeniya, SLcultivatedSandy loam82.06.511.51.420.0460.47[14]
Kurunegala, SLcultivatedSandy clay loam70.69.220.21.20 (0.018)0.0270.54This study
Polonnaruwa, SLcultivatedSandy clay loam61.216.921.91.20 (0.02)0.0430.53This study
Repacked Soils
Peradeniya-1, SLpastureSandy loam72.125.12.81.10 (0.05)0.100.57[14]
Peradeniya-2, SLpastureSandy loam72.125.12.81.30 (0.05)0.100.57[14]
Nishi-Tokyo, JPpastureSilt loamNA±NANA0.62 (0.05)NA0.74[17]
Peradeniya, SLcultivatedSandy loam82.06.511.51.420.0460.47[14]
Nishi-Tokyo, JPcultivatedSilty loamNA±NANA0.62 (0.05)NA0.74[17]
Kandy, SLcultivatedSandy loam61.520.518.01.05 (0.01)0.0370.57This study
Notes: SL: Sri Lanka; UK: United Kingdom; JP: Japan. ± NA: Not available.
Table 2. Selected predictive soil-gas diffusivity models used in this study.
Table 2. Selected predictive soil-gas diffusivity models used in this study.
Dp/Do ModelEquation
Buckingham (1904) D p D o = ε 2
Millington (1959) D p D o = ε 4 / 3
MQ (1961) D p D o = ε 10 / 3 2
WLR-Marshall
(Moldrup et al., 2000)
D p D o = ε 1.5 ε
SWLR
(Moldrup et al., 2000)
D p D o = ε 1 + C m ε
Notes: MQ, Millington and Quirk; WLR, Water-induced Linear Reduction; SWLR, Structure-dependent Water-induced Linear Reduction.
Table 3. Statistical analysis on model performance using RMSE (Equation (9)) and Bias (Equation (10)) for intact soils.
Table 3. Statistical analysis on model performance using RMSE (Equation (9)) and Bias (Equation (10)) for intact soils.
Dp/Do ModelPeradeniya-1, SLPeradeniya-2, SLNishi-Tokyo, JPTemuka, NZWakanui, NZTempleton, NZ
d = 1.0 g cm−3d = 1.3 g cm−3
RMSEBiasRMSEBiasRMSEBiasRMSEBiasRMSEBiasRMSEBias
Buckingham (1904)0.0437−0.0290.0594−0.04680.08040.00830.011−0.00880.0179−0.01370.00384−0.0027
Penman (1940)0.05340.04620.04180.03650.08150.07490.03510.02880.01730.00850.039790.0355
Marshal (1959)0.04540.03240.02930.01270.1160.07510.01210.00520.0083−0.00570.014240.01064
Millington (1959)0.07380.06150.05220.04210.14010.10630.0230.01490.0098−0.00010.025170.02033
MQ (1960)0.05690.02440.0392−0.00140.12280.04310.0085−0.00580.0155−0.01220.00352−0.00001
MQ (1961)0.0755−0.01330.0697−0.04680.1489−0.01220.0189−0.01460.0234−0.01690.00997−0.00761
WLR-Marshall0.0504−0.00720.0525−0.03370.11950.00720.0145−0.01160.0206−0.01540.0065−0.00524
SWLR0.0499−0.01220.0545−0.03670.1127−0.01230.0156−0.01240.0214−0.01580.00715−0.00572
Two Region0.01650.00230.01150.0010.01680.00360.00340.00090.0046−0.00210.004190.00108
Notes: MQ, Millington and Quirk; WLR, Water-induced Linear Reduction; SWLR, Structure-dependent Water-induced Linear Reduction.
Table 4. Statistical analysis on model performance using RMSE (Equation (5)) and Bias (Equation (6)) for repacked soils.
Table 4. Statistical analysis on model performance using RMSE (Equation (5)) and Bias (Equation (6)) for repacked soils.
Dp/Do ModelPasture SoilsCultivated Soils
Peradeniya-1, SLPeradeniya-2,Nishi-Tokyo, JPPeradeniya, SLNishi-Tokyo, JPKandy, SL
d = 1.0 g cm−3SL
d = 1.3 g cm−3
RMSEBiasRMSEBiasRMSEBiasRMSEBiasRMSEBiasRMSEBias
Buckingham (1904)0.0437−0.0290.0594−0.04680.08040.00830.011−0.00880.0179−0.01370.00384−0.0027
Millington (1959)0.07380.06150.05220.04210.14010.10630.0230.01490.0098−0.00010.025170.02033
MQ (1961)0.0755−0.01330.0697−0.04680.1489−0.01220.0189−0.01460.0234−0.01690.00997−0.00761
WLR-Marshall0.0504−0.00720.0525−0.03370.11950.00720.0145−0.01160.0206−0.01540.0065−0.00524
SWLR0.0499−0.01220.0545−0.03670.1127−0.01230.0156−0.01240.0214−0.01580.00715−0.00572
Notes: MQ, Millington and Quirk; WLR, Water-induced Linear Reduction; SWLR, Structure-dependent Water-induced Linear Reduction.
Table 5. Numerical parameterization of two-region soil-gas diffusivity model (Equations (6) and (7)) for pasture and cultivated soils.
Table 5. Numerical parameterization of two-region soil-gas diffusivity model (Equations (6) and (7)) for pasture and cultivated soils.
SoilSample StatusVegetation TypeCharacterization of Two-Region Soil-Gas Diffusivity Model(Equations (6) and (7))
woα1β1w1α2
Ambewela, SLIntactpasture00.410.991.00
Peradeniya, SLpasture0.050.141.100.430.09
Lexington, UKpasture0.0050.351.501.00
Peradeniya, SLcultivated0.0020.261.401.00
Kurunegala, SLcultivated0.20.473.011.00
Polonnaruwa, SLcultivated0.010.202.131.00
Peradeniya-1, SLRepackedpasture0.050.262.050.730.06
Peradeniya-2, SLpasture0.100.211.600.550.056
Nishi-Tokyo, JPpasture0.060.202.130.530.15
Peradeniya, SLcultivated0.050.232.500.700.08
Nishi-Tokyo, JPcultivated0.060.221.360.480.11
Kandy, SLcultivated0.050.133.740.580.18
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abeysinghe, A.M.S.N.; Lakshani, M.M.T.; Amarasinghe, U.D.H.N.; Li, Y.; Deepagoda, T.K.K.C.; Fu, W.; Fan, J.; Yang, T.; Ma, X.; Clough, T.; et al. Soil-Gas Diffusivity-Based Characterization of Variably Saturated Agricultural Topsoils. Water 2022, 14, 2900. https://doi.org/10.3390/w14182900

AMA Style

Abeysinghe AMSN, Lakshani MMT, Amarasinghe UDHN, Li Y, Deepagoda TKKC, Fu W, Fan J, Yang T, Ma X, Clough T, et al. Soil-Gas Diffusivity-Based Characterization of Variably Saturated Agricultural Topsoils. Water. 2022; 14(18):2900. https://doi.org/10.3390/w14182900

Chicago/Turabian Style

Abeysinghe, A. M. S. N., M. M. T. Lakshani, U. D. H. N. Amarasinghe, Yuan Li, T. K. K. Chamindu Deepagoda, Wei Fu, Jun Fan, Ting Yang, Xiaoyi Ma, Tim Clough, and et al. 2022. "Soil-Gas Diffusivity-Based Characterization of Variably Saturated Agricultural Topsoils" Water 14, no. 18: 2900. https://doi.org/10.3390/w14182900

APA Style

Abeysinghe, A. M. S. N., Lakshani, M. M. T., Amarasinghe, U. D. H. N., Li, Y., Deepagoda, T. K. K. C., Fu, W., Fan, J., Yang, T., Ma, X., Clough, T., Elberling, B., & Smits, K. (2022). Soil-Gas Diffusivity-Based Characterization of Variably Saturated Agricultural Topsoils. Water, 14(18), 2900. https://doi.org/10.3390/w14182900

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop