Spatial prediction of the exchangeable sodium percentage at multiple depths using electromagnetic inversion modelling

Agriculture in the semi‐arid and arid areas of the world requires irrigation. However, in these areas, soils naturally contain large amounts of sodium (sodic) which can cause amongst other things, surface crusting on the topsoil or structural instability in the subsoil. The exchangeable sodium percentage (ESP) needs to be mapped to guide the application of gypsum. Whilst geostatistical techniques, such as ordinary, co‐ and 3‐D kriging have been used, they have often been criticized because they are unable to take into account soil knowledge concerning distribution, processes and factors of formation. The use of digital soil mapping methods which couple remote or proximally sensed data with soil information is increasingly becoming useful because of the production of high‐resolution ancillary data. In this study, we first invert (using EM4Soil software) the electrical conductivity (σ_a –mS/m) of DUALEM‐421 data collected along a single transect. In doing this, we generate a 2‐dimensional electromagnetic conductivity image (EMCI). We couple the estimates of electrical conductivity (σ – mS/m) at 0.30 m depth increments down to 1.5 m with measured soil ESP. We compare the results of inversion using various possible coil array configurations of the DUALEM‐421 to determine a suitable set of data. We conclude that the use of the DUALEM‐41 is optimal (r² = 0.70). We use the calibration to estimate ESP along adjacent transects where we also generate EMCI. We are thus able to estimate ESP at various depths across a clay plain and an associated prior stream channel. We conclude that the collection of additional transects of DUALEM‐421 data as well as the use of a quasi‐3‐D inversion modelling approach would improve prediction.     

A seasonal behavior of surface soil moisture condition in a reclaimed tropical peatland

Soil moisture condition is essential to regulate the release of soil carbon from a drained peatland since aerobic microbial activities can be encouraged through oxygen supply associated with dewatering the soil layer while they may be discouraged under too dry conditions. Aiming to characterize the soil moisture condition in a reclaimed tropical peatland, we monitored the volumetric water content at 5?cm depth (θ _5?cm), groundwater level (GWL) and rainfall for 20 months from March 2010 to November 2011 in an oil palm field in Nakhon-Si-Thammarat, Thailand. We also measured the soil water retention curve and the unsaturated hydraulic conductivity (k) for a series of matric potential (h) to simulate the moisture condition monitored in the field by using the Buckingham-Darcy's flux law. During the dry season in 2010, the θ _5?cm consistently stayed lower than 0.35?m³?m^–3 with the GWL lower than a depth of 30?cm. In the transition from the dry season to the rainy season in 2010, the GWL rose to the land surface with peaks and dips across the time for about one month with the θ _5?cm increasing toward saturation. During the rainy season where the GWL stayed near or above the land surface, the θ _5?cm remained the field-saturated value of 0.58?m³?m^–3 on average, less than the laboratory-saturated value of 0.63?m³?m^–3, suggesting the development of a significant amount of entrapped air-phase. Hysteretic behavior in the measured θ _5?cm–GWL relation also supported that the top soil layer refuses to absorb water in wetting processes. The simulated θ _5?cm based on the measured k(h) and soil water retention curves demonstrated that the ease with which the top soil dries during a dry season was due mainly to the low k(h) value in the dried condition, while the slope of the θ(h) curve was so moderate that the soil layer could retain moisture for maintaining liquid water supply to the surface from the dropped GWL. Sensitivity analyses while varying the magnitude of both k(h) and evaporation rate (E) suggested that the k(h) function was more deterministic than the value of E in making the land surface easily dried. As the GWL stayed lower than 30?cm in depth for a total of 187 days out of the year monitored, while surface-ponding conditions took place for 120 days of the year, it was concluded that either the extremely dried condition or the saturated-moisture condition had dominantly occurred in the study site through a year and, thus, there may only be a limited time when soil organic matter near the land surface is in favorable moisture conditions for aerobic decomposition.     

Determination of chemical transport properties for different textures of undisturbed soils

Agrichemicals usually contaminate groundwater via preferential flow, therefore determination of the preferential flow characteristics of soil is needed. One model that predicts solute transport due to preferential flow is the mobile–immobile (MIM) solute-transport model, which partitions total water content (θ; m³ m^?3) into mobile (θ_m) and immobile fractions (θ_im). In undisturbed soils, a method is proposed for determining the MIM model parameters, i.e. immobile water fraction (θ_im), mass transfer coefficient (α) and hydrodynamic dispersion coefficient (D _h). Breakthrough curves were obtained for five different soil textures in three replicates, by miscible displacement of Cl^? in undisturbed soil columns. Cl^? breakthrough curves were evaluated in terms of the MIM model. Analysis suggests that the values of D _h and α increased with lighter soil textures and θ_im increased with heavier soil textures. The values of θ_im ranged from 5.31 to 14.28% in different soil textures. Furthermore, values of θ_im were found to be related to soil clay content. Values of α ranged from 0.0257 to 0.32 h^?1 and values of D _h ranged from 0.36 to 11.2 cm² h^?1 in different soil textures. A significant linear correlation was obtained between α, θ_im, D _h and soil saturated hydraulic conductivity (K _s) and pore water velocity (v). A multivariate pedotransfer function was developed to estimate α, θ_im and D _h based on the geometric mean (d _g) and the standard deviation (σ_g) of the diameter of soil particles and soil organic matter content. The pedotransfer functions for D _h, θ_im and α were validated by independent data sets from other investigators.     

Modeling within field variation of the compaction layer in a paddy rice field using a proximal soil sensing system

A key characteristic of flooded paddy fields is the plough pan. This is a sub‐soil layer of greater compaction and bulk density, which restricts water losses through percolation. However, the thickness of this compacted layer can be inconsistent, with consequences such as variable percolation and leaching losses of nutrients, which therefore requires precision management of soil water. Our objective was to evaluate a methodology to model the thickness of the compacted soil layer using a non‐invasive electromagnetic induction sensor (EM38‐MK2). A 2.7 ha alluvial non‐saline paddy rice field was measured with a proximal soil sensing system using the EM38‐MK2 and the apparent electrical conductivity (EC_a) of the wet paddy soil was recorded at a high‐resolution (1.0 × 0.5 m). Soil bulk density (n = 10) was measured using undisturbed soil cores, which covered locations with large and small EC_a values. At the same locations (within 1 m²) the depth of the different soil layers was determined by penetrometer. Then a fitting procedure was used to model the EC_a – depth response functions of the EM38‐MK2, which involved solving a system of non‐linear equations and a R² value of 0.89 was found. These predictions were evaluated using independent observations (n = 18) where a Pearson correlation coefficient of 0.87 with an RMSEE value of 0.03 m was found. The EC_a measurements allowed the detail estimation of the compacted layer thickness. The link between water percolation losses and thickness of the compacted layer was confirmed by independent observations with an inverse relationship having a Pearson correlation coefficient of 0.89. This rapid, non‐invasive and cost‐effective technique offers new opportunities to measure differences in the thickness of compacted layers in water‐saturated soils. This has potential for site‐specific soil management in paddy rice fields.     

Pedo-transfer functions for estimating the hydraulic properties of paddy soils in subtropical central China

Pedo-transfer functions (PTFs) have been widely used to estimate soil hydraulic properties in the simulation of catchment eco-hydrological processes. However, the accuracy of existing PTFs is usually inadequate for use. To develop PTFs for local use, soil columns were collected from a double rice-cropped agricultural catchment in subtropical central China. The PTFs for saturated soil hydraulic conductivity (K_s) and parameters (θ_s, α, and n) of the van Genuchten model for the soil water retention curve (SWRC) were obtained based on soil’s basic properties, and compared with models developed by Li et al. in 2007 and Wösten et al. in 1999, respectively. Our results indicated that K_s in the range of 0.04–1087 cm d^?1 and θ_s in the range of 0.34–0.51 cm³ cm^?3 were both well estimated with the R²_adj of 0.72 and 0.87, respectively, but α (0.04–0.65 cm^?1) and n (1.05–1.21) were relatively poorly predicted with the respective R²_adj of 0.38 and 0.55, despite the use of more input parameters. Our local derived PTFs outperformed the other two existing models. However, if the local PTFs for paddy soils are not available, the Wösten et al. 1999 model can be proposed as a useful alternative. Therefore, this study can improve our understanding of the development and application of PTFs for predicting paddy soil hydraulic properties in China.     

Utilizing a DUALEM‐421 and inversion modelling to map baseline soil salinity along toposequences in the Hunter Valley Wine district

In the oldest commercial wine district of Australia, the Hunter Valley, there is the threat of soil salinization because marine sediments underlie the area. To understand the risk requires information about the spatial distribution of soil properties. Electromagnetic (EM) induction instruments have been used to identify and map the spatial variation of average soil salinity to a certain depth. However, soils vary with depth dependent on soil forming factors. We collected data from a single‐frequency and multiple‐coil DUALEM‐421 along a toposequence. We inverted this data using EM4Soil software and evaluated the resultant 2‐dimensional model of true electrical conductivity (σ – mS/m) with depth against electrical conductivity of saturated soil pastes (EC_p – dS/m). Using a fitted linear regression (LR) model calibration approach and by varying the forward model (cumulative function‐CF and full solution‐FS), inversion algorithm (S1 and S2), damping factor (λ) and number of arrays, we determined a suitable electromagnetic conductivity image (EMCI), which was optimal (R² = 0.82) when using the full solution, S2, λ = 3.6 and all six coil arrays. We conducted an uncertainty analysis of the LR model used to estimate the electrical conductivity of the saturated soil‐paste extract (EC_e – dS/m). Our interpretation based on estimates of EC_e suggests the approach can identify differences in salinity, how these vary with parent material and how topography influences salt distribution. The results provide information leading to insights into how soil forming factors and agricultural practices influence salinity down a toposequence and how this can guide soil management practices.     

Short-Term Impacts of Gray Water Irrigation on Soil Characteristics in Land-Treatment Vegetation Filters

This study on soil quality parameters was conducted from December 2008 to July 2010 to see the short-term impacts of gray water irrigation on physical and chemical properties of soil. The results envisaged slight positive impacts of wastewater irrigation on upper soil layers (0–30 cm) with an increase in soil organic matter (SOM) from 1.45% to 1.89% and nutrient content increases from 199 to 267 kg ha^?1 nitrogen (N), from 12.4 to 27.1 kg ha^?1 phosphorus (P), and from 131 to 136 kg ha^?1 potassium (K). Slight changes were also recorded in the soil pH and electrical conductivity (EC), which increased from 6.39 to 6.81 and 145 to 229 μS cm^?1 respectively, in the upper soil layers. The soil bulk density (BD) slightly decreased to 1.34 g cm³ in the treatment plot from 1.45 g cm³ at the commencement of the trial, which indicates good soil conditioning under the gray water irrigation trial.     

Laboratory calibration of time domain reflectometry to determine moisture content in undisturbed peat samples

Time domain reflectometry (TDR), while widely used to measure volumetric water content (θ) and bulk electrical conductivity (BEC) in unsaturated granular soils, remains less studied in peat than mineral soils. Empirical models commonly used in mineral soils are not applicable to peat for accurate determination of θ from measured apparent dielectric permittivity (?). Past studies for peat report highly variable calibrations, and suggest differences in origin of organic matter, degree of decomposition and bound water to explain such variability. This study shows that bound water appears to have minimal impact on calibration because of its negligible volumetric fraction at the low bulk densities of peat. Increased volumetric air fraction at the same θ values attributed to high porosity of peat makes the ?‐θ relationships of mineral soils inapplicable. Temperature effects on ? resulted in a correction factor for θ. The temperature correction factor decreased with decreasing θ and was determined experimentally to lie between ?0.0021 m³ m^?3 per °C for θ≥ 0.79 m³ m^?3 and ?0.0005 m³ m^?3 per °C for θ = 0.35 m³ m^?3. The decreasing value of the correction factor with θ can be explained by dependence of the ?‐θ relationship on properties of free water alone. Temperature dependence of BEC was close to that of soil solution. Maxwell‐De Loor's four‐phase mixing model (MDL) based on physical properties of the multiphase soil system can efficiently simulate the effect of increased air volume and varying soil temperature on the ?‐θ relationship in peat. In addition, linear ?‐θ calibration in peat can be improved when BEC is included in the calibration equation.     

Mapping soil salinity in 3‐dimensions using an EM38 and EM4Soil inversion modelling at the reconnaissance scale in central Morocco

Large areas of Morocco require irrigation and although good quality water is available in dams, farmers augment river water with poorer quality ground water, resulting in salt build‐up without a sufficient leaching fraction. Implementation of management plans requires baseline reconnaissance maps of salinity. We developed a method to map the distribution of salinity profiles by establishing a linear regression (LR) between calculated true electrical conductivity (σ, mS/m) and electrical conductivity of the saturated soil‐paste extract (ECe, dS/m). Estimates of σ were obtained by inverting the apparent electrical conductivity (ECa, mS/m) collected from a 500‐m grid survey using an EM38. Spherical variograms were developed to interpolate ECa data onto a 100 m grid using residual maximum likelihood. Inversion was carried out on kriged ECa data using a quasi‐3d model (EM4Soil software), selecting the cumulative function (CF) forward modelling and S2 inversion algorithm with a damping factor of 3.0. Using a ‘leave‐one‐out cross‐validation' (LOOCV), of one in 12 of the calibration sites, the use of the q‐3d model yielded a high accuracy (RMSE = 0.42 dS/m), small bias (ME = ?0.02 dS/m) and Lin's concordance (0.91). Slightly worse results were obtained using individual LR established at each depth increment overall (i.e. RMSE = 0.45 dS/m; ME = 0.00 dS/m; Lin's = 0.89) with the raw EM38 ECa. Inversion required a single LR (ECe = 0.679 + 0.041 × σ), enabling efficiencies in estimating ECe at any depth across the irrigation district. Final maps of ECe, along with information on water used for irrigation (ECw) and the characterization of properties of the two main soil types, enabled better understanding of causes of secondary soil salinity. The approach can be applied to problematic saline areas with saline water tables.     

Detecting a leachate plume in an aeolian sand landscape using a DUALEM‐421 induction probe to measure electrical conductivity followed by inversion modelling

Solid waste poses a serious health risk when it is disposed of inadequately because water‐based solutions derived from the decomposition of solid waste products (leachate) can enter groundwater systems via plumes. To assess the public health risk and potential ecological impacts, we require knowledge on the pedological and hydrogeological settings in which waste is disposed. This is particularly the case in coarse textured highly permeable soil. To rapidly collect data, geophysical methods such as direct current (dc) resistivity techniques have been used. Moreover, non‐contact electromagnetic (EM) induction instruments have also been employed. The aim of this research was to demonstrate how the inversion using a 1‐dimensional inversion algorithm with lateral constraints of the apparent electrical conductivity (σ_a) measured in the horizontal coplanar (HCP) and perpendicular co‐planar arrays (PRP) of a DUALEM‐421 EM induction probe can be used to develop a two‐dimensional model of the true electrical conductivity (σ) within a Quaternary aeolian sand in the Tuggerah Soil Landscape southeast of Sydney in Australia. Our results from 2D models of σ accord with estimates of bulk electrical conductivity (σ_b) of a leachate plume and uncontaminated groundwater, the stratigraphy of the Tuggerah soil landscape unit and the depth of sand used to landscape the decommissioned landfill. Further research is needed to determine the origin of the plume and a quasi‐3D modelling approach is applicable.     

Compatibility of root growth and tuber production of potato cultivars with dynamic and static water‐saving irrigation managements

The important root characteristics of root length density (RLD) and root mass density (RMD) generally differ among irrigation managements and potato cultivars. The objective of this study was to investigate the RLD and RMD variations and their functional relationships with gross potato tuber yield for two commercial potato cultivars, Agria and Sante, under different irrigation strategies. Full irrigation and water‐saving irrigation strategies, deficit and partial root drying irrigations, were applied statically (S) and dynamically (D) based on daily crop evapotranspiration. Results showed that SPRD had significantly greater RLD (3.64 cm/cm³) and RMD (132.7 μg/cm³) than other irrigation treatments. Between the potato cultivars, Agria had significantly larger values of RLD (3.50 cm/cm³) and RMD (138.7 μg/cm³) than Sante. The functional relationship between the root growth characteristics and tuber yield showed that under water‐saving irrigations, Agria increased root mass at the expense of gross tuber yield but Sante increased root mass to maintain larger gross tuber yields. However, Agria produced more roots and gross tuber yield than Sante, and it is concluded that Agria is a more drought‐tolerant potato cultivar, which is recommended for tuber production in regions where water might be scarce. It was shown that larger root production in potatoes was associated with improved tolerance to water stress.     

The conversion of permittivity as measured by a PR2 capacitance probe into soil moisture values for Des Moines lobe soils in Iowa

Measurement of soil moisture is essential for irrigation scheduling and capacitance sensors have been widely used to monitor soil moisture at different depths. Two approaches for converting permittivity measures using the capacitance probe (PR2, Delta‐T Devices) to soil water content are to (a) use the default equation and parameters provided by the manufacturer, and (b) use site specific calibration equations. The objective of this study was to evaluate the performance of the manufacturer’s default equation and in‐situ calibrated equations for estimating soil water content. Permittivity measurement using the PR2 probe coincided with soil sampling during the growing seasons in 2006, 2007 and 2008 for Des Moines lobe soils in north‐central Iowa. The default equation provided by Delta‐T Devices for the PR2 probe estimated the soil water content for 3 years with an average root mean square error (RMSE) and index of agreement (IoA) values of 0.097 cm³/cm³ and 0.587, respectively. The default equation was calibrated by a 1‐year (2006) and a 2‐year (2006 + 2007) data set. The resultant statistics indicate that site specific calibration gives more accurate estimates of soil water content compared to the uncalibrated default equation. Three‐year averaged RMSE and IoA values were 0.049 cm³/cm³ and 0.742 for equations calibrated by the 1‐year data set, and 0.034 cm³/cm³ and 0.807 for equations calibrated by the 2‐year data set. The results from this study indicate that a site specific calibration is necessary for the PR2 probe, and equations calibrated by data from a longer period performed better than data from a shorter period. Where a site‐specific field calibration cannot be applied, coefficients are provided for various cropping systems in Des Moines Lobe soils of Iowa based on the results from this study.     

Coupled effects of biochar use and farming practice on physical properties of a salt-affected soil with wheat–maize rotation

Purpose

Being carbon-rich and porous, biochar has the potential to improve soil physical properties, so does conventional farming practice. Here, a field trial was conducted to investigate the combined effects of biochar use and farming practice on the physical properties of a salt-affected compact soil for wheat–maize rotation in the Yellow River Delta region.

Materials and methods

Salix fragilis L. was used as feedstock to produce biochar in the field via aerobic carbonization at an average temperature of 502 °C, terminated by a water mist spray, for use as a soil amendment at 0, 1, 2, and 4 g kg^?1 doses (CK, T1, T2, and T3, respectively). Farming practices included rotary tillage/straw returning for wheat sowing, spring irrigation, no-tillage seeding of maize, and autumn irrigation. Both cutting ring and composite samples of the soil were collected at four stages of wheat–maize rotation (22, 238, 321, and 382 d after the benchmark date of land preparation for wheat sowing) for the determination of soil properties by established methods.

Results and discussion

Rotary tillage/straw returning reduced soil bulk density (BD) from 1.48 to 1.27 g cm^?3 (CK) and 1.14 g cm^?3 (T3) and increased saturated hydraulic conductivity (Ks) from 0.05?×?10^?5 to 0.75?× 10^?5 cm s^?1 (CK) and 1.25?× 10^?5 cm s^?1 (T3). This tillage effect on BD and Ks gradually disappeared due to the disturbance from the subsequent farming practice. Biochar use lessened the disturbance. At maize harvest, BD was 1.47 (CK) vs. 1.34 g cm^?3 (T3), and Ks was 0.06?×?10^?5 (CK) vs. 0.28?×?10^?5 cm s^?1(T3); in comparison with CK, T3 increased Na⁺ leaching by 65%, Cl^? leaching by 98%, organic carbon content by 40.3%, and water-stable aggregates (0.25–2 mm) by 38%, indicating an improvement in soil properties.

Conclusions

Biochar use and rotary tillage improved soil physical properties (BD, Ks) and favored soil aeration, water filtration, and salt leaching, which further helped the accumulation of soil organic carbon, the formation of water-stable aggregates, and the amelioration of salt-affected compact soil.

     

Evaluating soil salinity status from bulk electrical conductivity and permittivity

In the range of volumetric water content, θ, from about 0.12 cm³ cm^–3 to saturation the relation between bulk electrical conductivity, C_b, and bulk electrical permittivity, ε, of mineral soils was observed to be linear. The partial derivative ?C_b/?ε appeared independent of the moisture content and directly proportional to soil salinity. We found that the variable X_s = ?C_b/?ε determined from in situ measurements of C_b(θ > 0.2) and ε(θ > 0.2) can be considered as an index of soil salinity, and we call it the ‘salinity index’. Knowing the index and sand content for a given soil we could calculate the electrical conductivity of the soil water, C_w, which is a widely accepted measure of soil salinity. The two variables from which the salinity index can be calculated, i.e. C_b and ε, can be read simultaneously from the same sensor by time-domain reflectometry. Quantities and symbols a constant /dS m^–1 b constant c constant /dS m^–1 C _b electrical conductivity of bulk soil /dS m^–1 C _b′ constant equal to 0.08 dS m^–1 C _s electrical conductivity of a solution used to moisten soil samples /dS m^–1 C _w electrical conductivity of soil water defined as the soil salinity /dS m^–1 C _wref reference salinity (that truly existing) resulting from the procedure of moistening samples, expressed as C_s + C_r/dS m^–1 C _r baseline value of C_s due to residual soluble salts present in the soil /dS m^–1 d constant D dry soil bulk density /g cm^–3 l slope r ratio S sand content /% by weight t time /s X _s salinity index /dS m^–1 X _si initial salinity index when distilled water is used to moisten soil samples /dS m^–1 Y a moisture-independent salinity-dependent variable /dS m^–1 z coordinate along direction of flow of the soil solution ε′ constant equal to 6.2 ε relative bulk electrical permittivity (dielectric constant) of the soil θ volumetric water content determined thermogravimetrically using oven-drying /cm³ cm^–3     

Development of an apparatus for measuring one‐dimensional steady‐state heat flux of soil under reduced air pressure

One of the best ways to evaluate the coupled heat and mass transfer in soil is to measure the heat flux and water distribution simultaneously. For this purpose, we developed an apparatus for measuring the one‐dimensional steady‐state heat flux and water distribution in unsaturated soil under reduced air pressure. The system was tested using four samples with known thermal conductivity (0.6–8.0 W m^?1 K^?1). We confirmed that the system could measure the one‐dimensional steady‐state heat flux under a fixed temperature difference between ends of the samples over a wide range of thermal conductivity values. Time domain reflectometry was used to measure the water distribution with a repeatability of less than ± 1.0%. We used the apparatus to measure the soil heat flux and distribution of water content and temperature under steady‐state conditions with reduced air pressure. The initial volumetric water content, θ_ini, of the soil samples was set at 0.20 and 0.40 m³m^?3. For a θ_ini of 0.20, the heat flux was not significantly affected by air pressure, and the water content on the hot side decreased whilst that on the cold side increased, i.e. a pronounced water content gradient was formed. For a θ_ini of 0.40, the heat flux increased sharply with reduced air pressure, and the water content did not change, i.e. a homogeneous water distribution was observed. The increase in the heat flux with air pressure reduction is caused by the vapour transfer in soil pores. We found that a large vapour transfer took place in the soil with the homogeneous water distribution, and that the vapour transfer was less in the soil with the pronounced water content gradient. These experimental facts were entirely different from the traditional knowledge of vapour transfer in soil under temperature gradients. A lack of data on heat flux must have resulted in the previously incorrect conclusions. The new apparatus will serve to clarify the intricate phenomena of thermally induced vapour transfer in unsaturated soil in further experiments.     

土壤垂向分层和均匀处理下水分差异的数值探讨

在现有众多的陆面过程模型中,对土壤水分的定量描述一般是假设垂向分布均匀,取表层土壤质地来表示整个垂向土壤质地。垂向分层和均匀处理下的土壤水分是存在差异的,这种差异有多大目前少有研究。设置3组不同饱和导水率组合的层状土壤代表不同区域的非均匀土壤,取3组层状土壤的上层土壤代表整个均匀土壤,通过建立一维土壤水分运动模型分析这种差异,同时分析饱和导水率、饱和含水量、残余含水率、孔隙大小分布参数和形状参数对层状土壤和均匀土壤的渗透量和储水量差异的敏感性,探讨垂向层状和均匀处理下土壤水分运动的差异。研究结果表明:1)建立的一维土壤水分运动模型模拟的土壤水分剖面与Yeh解析解和室内五水转化试验的土壤水分剖面一致,表明模型无论是考虑还是不考虑根系吸水都具有可靠性。2)采用垂向均匀方式处理,上下层饱和导水率相差越大的层状土壤,各水文变量的差异越大。当层状土壤上下层饱和导水率相差1.5倍时,层状土壤和均匀土壤的水分分布差别小于0.05 cm~3×cm~(-3);而当层状土壤上下层饱和导水率相差达3.3倍时,层状土壤和均匀土壤的水分分布差别达0.15 cm~3×cm~(-3),渗漏量相差20 cm以上,储水量相差5 cm左右。3)相对于层状土壤下层,均匀土壤下层的持水能力更差,水流速度更快,导致下层水分分布减小,渗漏量增加,储水量减小。4)形状参数n对渗透量的敏感性最强,土壤孔隙大小分布参数对储水量的敏感性最强,形状参数n其次。在实际应用中,如果一个区域的土壤上下层饱和导水率相差较大,那么垂向均匀处理可能会导致很大的误差,和实际土壤的水分分布相差很大,这会严重影响土壤水分的准确估计,在实际处理中需要认真考虑。     

Potential to predict depth‐specific soil–water content beneath an olive tree using electromagnetic conductivity imaging

Efficient monitoring of soil moisture is becoming increasingly important. To understand soil–plant–water dynamics, we evaluate the potential of using a multiple‐coil‐array electromagnetic induction instrument and inversion software to map soil moisture beneath an olive tree. On twelve different days, we collected apparent electrical conductivity (EC _a) data using a DUALEM ‐21S and the volumetric soil moisture (θ ) using a bank of soil moisture sensors on opposite sides of the tree. Using EM 4Soil, we inverted the EC _a data on five of the days and established a site‐specific calibration between estimates of true electrical conductivity (σ ) and θ . The strongest calibration relationship between σ and θ (R ² = 0.65) was obtained for a full‐solution, S2 algorithm and damping factor of 1.2. A leave one out cross‐validation (LOOCV ) showed the calibration was robust, with a root mean square error (RMSE ) of 0.046 m³/m³, a mean error (ME ) of 0.001 m³/m³ and a Lin's concordance of 0.72. We subsequently evaluated the calibration relationship on the seven remaining days and over a drying period of 120 days. This approach provides information about the temporal evolution of θ by a LOOCV of validation with a RMSE of 0.037, ME of −0.003 and a Lin's concordance of 0.54. Improvement could be achieved by aligning the DUALEM ‐21S in the same orientation as the sensors, with time‐lapse inversion also being advantageous.     

TDR estimation of electrical conductivity and saline solute concentration in a volcanic soil

Relative to montmorillonitic or kaolinitic soils, volcanic soils have atypical dielectric characteristics that interfere with the applicability of the Time Domain Reflectometry (TDR) technique for soil moisture (θ) determination when common, empirical calibration equations are used. This particular dielectric response affects estimation of salinity in volcanic soils. Six TDR-based methods to estimate bulk electrical conductivity (σ_a) on a range of KCl saline reference solutions were compared, with Nadler's method giving the best results (R_1:1²=0.988). Three models (linear, non-linear and empirical) for predicting soil solution electrical conductivity (σ_w) based on σ_a and θ, were experimentally tested on 24 hand-packed soil columns varying in salinity (Br⁻) from 0.2 to 4.0 dS m⁻¹, each in four θ levels (36–58%). Rhoades' linear model performed better, especially for large water contents, than the other two (R_1:1²=0.986 vs. 0.976 and 0.983, respectively). An interpretation in terms of mobile vs. immobile volumetric fractions of water present in volcanic soils is suggested as a possible explanation for these results. The empirical model resulted over-parameterized and an alternative equation with fewer non-correlated parameters, σ_a=(aθ²+bθ)σ_w+cθ², is proposed and tested with good results in volcanic soils from the Canary Islands and New Zealand. The equation encompasses both the relative dielectric dominance of the mobile water fraction at high water content typical of volcanic soils, and of the immobile fraction at low water contents. Simultaneous measurements made with a standard four-electrode probe and TDR gave good correlation (R²=0.964). A good linear correlation was also found between tracer concentration in the soil solution and σ_w (R²=0.960). Nadler's and the new empirical model also tested with good results under dynamic (flow) conditions during a miscible displacement experiment in a large monolith using bromide as a tracer. The method reveals itself as a robust tool for solute transport studies under controlled salinity conditions in a volcanic soil.     

Using time domain reflectometry for monitoring mineralization of nitrogen from soil organic matter

The mineralization of nitrogen from soil organic matter is important when one tries to optimize nitrogen fertilization and assess risks of N losses to the environment, but its measurement is laborious and expensive. We have explored the possibilities for monitoring N mineralization directly using time domain reflectometry (TDR). Net N and S mineralization were monitored over a 101‐day period in two layers (0–30 and 30–60 cm) of a loamy sand soil during aerobic incubation in a laboratory experiment. At the same time electrical conductivity of the bulk soil, σ_a, was measured by TDR. A series of calibration measurements with different amounts of KNO₃ at different soil moisture contents was made with the topsoil to calculate the electrical conductivity, σ_w, of the soil solution from σ_a and θ. The actual σ_w was determined from the conductivity of 1:2 soil:water extracts (σ_1:2) with a mass balance approach using measured NO₃^– concentrations, after correction for ions present prior to the addition of KNO₃. The average N mineralization rate in the topsoil was small (0.12 mg N kg^?1 day^?1), and, as expected, very small in the subsoil (0.023 mg N kg^?1 day^?1). In the top layer NO₃^– concentrations calculated from σ_a determined by TDR slightly underestimated measured concentrations in the first 4 weeks, and in the second half of the incubation there was a significant overestimation of measured NO₃^–. Using the sum of both measured NO₃^– and SO₄^2– reduced the overestimation. In the subsoil calculated NO₃^– concentrations strongly and consistently overestimated measured concentrations, although both followed the same trend. As S mineralization in the subsoil was very small, and initial SO₄^2– concentrations were largely taken into account in the calibration relations, SO₄^2– concentrations could not explain the overestimation. The very small NO₃^– and SO₄^2– concentrations in the B layer, at the lower limit of the concentrations used in the calibrations, are a possible explanation for the discrepancies. A separate calibration for the subsoil could also be required to improve estimates of NO₃^– concentrations.     

Effect of deep ploughing on the water status of highly and less compacted soils for coconut (Cocos nucifera L.) production in Sri Lanka

Soil compaction limits soil water availability which adversely affects coconut production in Sri Lanka. Field experiments were conducted in coconut (Cocos nucifera L.) plantations with highly and less compacted soils in the intermediate climatic zone of Sri Lanka. Soil physical properties of sixteen major soil series planted with coconut were evaluated to select the most suitable soil series to investigate the effect of deep ploughing on soil water conservation. Soil compaction and soil water retention with respect to deep ploughing were monitored during the dry and rainy seasons using cone penetrometer and neutron scattering techniques, respectively. Evaluation of soil physical properties showed that the range of mean values of bulk density (BD) and soil penetration resistance (SPR) in the surface soil (0–10 cm depth) of major soil series in coconut lands was from 1.38 ± 0.02 to 1.57 ± 0.07 g/cm³ and 55 ± 10 to 315 ± 16.4 N/cm² respectively. The total available water fraction increased with clay content of soil as a result of high micropores. However, due to soil compaction, ability of soils to conserve water and to remain aerated was low for those series. Deep ploughing during the rainy and dry periods in highly compacted soils (BD > 1.5 g/cm³ and SPR > 250 N/cm²) greatly increased conserved soil water in the profile, while in less compacted soils (BD < 1.5 g/cm³ and SPR < 250 N/cm²) conserved water content was adversely affected. Soil water retention in bare soils of both highly and less compacted soil series was higher than that of live grass-covered soil. Amount of water conserved in ploughed Andigama series with respect to bare soils and grass-covered treatments during the severe dry period was 10.4 and 16.9 cm/m, while water storage reduction in the same treatments with ploughed Madampe series was 6.55 and 5.45 cm/m respectively. In addition, deep ploughing even in the effective root zone with live grass-covered highly compacted soils around coconut tree was favorable for soil water retention compared to that of live grass-covered less compacted soils.