Selecting optimal calibration samples using proximal sensing EM induction and γ-ray spectrometry data: An application to managing lime and magnesium in sugarcane growing soil.

Calcium (Ca) and magnesium (Mg) are essential for growth of sugarcane leaves and roots, as well as respiration and nitrogen metabolism, respectively. To assist farmers decide suitable application rates of lime and Mg fertiliser, respectively, the Australian sugarcane industry established the Six-Easy-Steps nutrient management guidelines based on topsoil (0-0.3 m) Ca (cmol(+) kg-1) and Mg (cmol(+) kg-1). Given the heterogeneous nature of soil, digital soil mapping (DSM) methods can be employed to allow for the precise application rate to be determined. In this study, we examine statistical models (i.e., ordinary kriging [OK], linear mixed model [LMM], quantile regression forests [QRF], support vector machine [SVM], and Cubist regression kriging [CubistRK]) to predict topsoil and subsoil (0.6-0.9) Ca and Mg, employing digital data in combination (i.e., proximal sensing electromagnetic induction (EMI) [e.g., 1mPcon, 1mHcon, etc.], gamma-ray [γ-ray] spectrometry [i.e., TC, K, U and Th] and digital elevation model [DEM] derivatives). We also investigate various sampling designs (i.e., spatial coverage [SCS], feature space coverage [FSCS], conditioned Latin hypercube [cLHS] and simple random sampling [SRS]) and calibration sample size (i.e., n = 180, 150, 120, 90, 60 and 30). The predictions are assessed using Lin's concordance correlation coefficient (LCCC) and ratio of performance to interquartile distance (RPIQ) with an independent validation dataset (i.e., n = 40). The best results were for prediction of subsoil Mg, utilising CubistRK (LCCC = 0.82) with the largest calibration sample size (n = 180), followed by LMM (0.79), SVM (0.76), QRF (0.70) and OK (0.65). This was generally the case for topsoil and subsoil Ca. We also conclude that no single sampling design was universally better, and 180 samples were necessary for predicting topsoil Ca and Mg with moderate agreement (0.65 < LCCC < 0.80). However, with FSCS, a minimum of 120 samples were enough to calibrate a CubistRK model and achieve substantial (LCCC > 0.80) agreement for predicting subsoil Ca and Mg. With respect to soil use and management according to the Six-Easy-Steps, the sandy soil in the north and south require large application rate of lime (3.5 t/ha) and Mg (125 kg/ha), respectively. Conversely, varying amounts of fertiliser rates of lime (2.0, 1.5 and 1 t/ha) and Mg (50 kg/ha) were recommended where Vertosols were previously mapped.

Three-Dimensional Mapping of Clay and Cation Exchange Capacity of Sandy and Infertile Soil Using EM38 and Inversion Software.

Most cultivated upland areas of northeast Thailand are characterized by sandy and infertile soils, which are difficult to improve agriculturally. Information about the clay (%) and cation exchange capacity (CEC-cmol(+)/kg) are required. Because it is expensive to analyse these soil properties, electromagnetic (EM) induction instruments are increasingly being used. This is because the measured apparent soil electrical conductivity (ECa-mS/m), can often be correlated directly with measured topsoil (0-0.3 m), subsurface (0.3-0.6 m) and subsoil (0.6-0.9 m) clay and CEC. In this study, we explore the potential to use this approach and considering a linear regression (LR) between EM38 acquired ECa in horizontal (ECah) and vertical (ECav) modes of operation and the soil properties at each of these depths. We compare this approach with a universal LR relationship developed between calculated true electrical conductivity (σ-mS/m) and laboratory measured clay and CEC at various depths. We estimate σ by inverting ECah and ECav data, using a quasi-3D inversion algorithm (EM4Soil). The best LR between ECa and soil properties was between ECah and subsoil clay (R2 = 0.43) and subsoil CEC (R2 = 0.56). We concluded these LR were unsatisfactory to predict clay or CEC at any of the three depths, however. In comparison, we found that a universal LR could be established between σ with clay (R2 = 0.65) and CEC (R2 = 0.68). The LR model validation was tested using a leave-one-out-cross-validation. The results indicated that the universal LR between σ and clay at any depth was precise (RMSE = 2.17), unbiased (ME = 0.27) with good concordance (Lin's = 0.78). Similarly, satisfactory results were obtained by the LR between σ and CEC (Lin's = 0.80). We conclude that in a field where a direct LR relationship between clay or CEC and ECa cannot be established, can still potentially be mapped by developing a LR between estimates of σ with clay or CEC if they all vary with depth.

The location- and scale- specific correlation between temperature and soil carbon sequestration across the globe.

Much research has been conducted to understand the spatial distribution of soil carbon stock and its temporal dynamics. However, an agreement has not been reached on whether increasing global temperature has a positive or negative feedback on soil carbon stocks. By analysing global maps of soil organic carbon (SOC) using a spherical wavelet analysis, it was found that the correlation between SOC and soil temperature at the regional scale was negative between 52° N and 40° S parallels and positive beyond this region. This was consistent with a few previous studies and it was assumed that the effect was most likely due to the temperature-dependent SOC formation (photosynthesis) and decomposition (microbial activities and substrate decomposability) processes. The results also suggested that the large SOC stocks distributed in the low-temperature areas might increase under global warming while the small SOC stocks found in the high-temperature areas might decrease accordingly. Although it remains unknown whether the potential increasing soil carbon stocks in the low-temperature areas can offset the loss of carbon stocks in the high-temperature areas, the location- and scale- specific correlations between SOC and temperature should be taken into account for modeling SOC dynamics and SOC sequestration management.

Evaluating a Bayesian modelling approach (INLA-SPDE) for environmental mapping.

Understanding the uncertainty in spatial modelling of environmental variables is important because it provides the end-users with the reliability of the maps. Over the past decades, Bayesian statistics has been successfully used. However, the conventional simulation-based Markov Chain Monte Carlo (MCMC) approaches are often computationally intensive. In this study, the performance of a novel Bayesian inference approach called Integrated Nested Laplace Approximation with Stochastic Partial Differential Equation (INLA-SPDE) was evaluated using independent calibration and validation datasets of various skewed and non-skewed soil properties and was compared with a linear mixed model estimated by residual maximum likelihood (REML-LMM). It was found that INLA-SPDE was equivalent to REML-LMM in terms of the model performance and was similarly robust with sparse datasets (i.e. 40-60 samples). In comparison, INLA-SPDE was able to estimate the posterior marginal distributions of the model parameters without extensive simulations. It was concluded that INLA-SPDE had the potential to map the spatial distribution of environmental variables along with their posterior marginal distributions for environmental management. Some drawbacks were identified with INLA-SPDE, including artefacts of model response due to the use of triangle meshes and a longer computational time when dealing with non-Gaussian likelihood families.

Irrigation salinity hazard assessment and risk mapping in the lower Macintyre Valley, Australia.

In the Murray-Darling Basin of Australia, secondary soil salinization occurs due to excessive deep drainage and the presence of shallow saline water tables. In order to understand the cause and best management, soil and vadose zone information is necessary. This type of information has been generated in the Toobeah district but owing to the state border an inconsistent methodology was used. This has led to much confusion from stakeholders who are unable to understand the ambiguity of the results in terms of final overall risk of salinization. In this research, a digital soil mapping method that employs various ancillary data is presented. Firstly, an electromagnetic induction survey using a Geonics EM34 and EM38 was used to characterise soil and vadose zone stratigraphy. From the apparent electrical conductivity (ECa) collected, soil sampling locations were selected and with laboratory analysis carried out to determine average (2-12m) clay and EC of a saturated soil-paste extract (ECe). EM34 ECa, land surface parameters derived from a digital elevation model and measured soil data were used to establish multiple linear regression models, which allowed for mapping of various hazard factors, including clay and ECe. EM38 ECa data were calibrated to deep drainage obtained from Salt and Leaching Fraction (SaLF) modelling of soil data. Expert knowledge and indicator kriging were used to determine critical values where the salinity hazard factors were likely to contribute to a shallow saline water table (i.e., clay ≤35%; ECe>2.5dS/m, and deep drainage >100mm/year). This information was combined to produce an overall salinity risk map for the Toobeah district using indicator kriging. The risk map shows potential salinization areas and where detailed information is required and where targeted research can be conducted to monitor soil conditions and water table heights and determine best management strategies.

Modeling Coastal Salinity in Quasi 2D and 3D Using a DUALEM-421 and Inversion Software.

Rising sea levels, owing to climate change, are a threat to fresh water coastal aquifers. This is because saline intrusions are caused by increases and intensification of medium-large scale influences including sea level rise, wave climate, tidal cycles, and shifts in beach morphology. Methods are therefore required to understand the dynamics of these interactions. While traditional borehole and galvanic contact resistivity (GCR) techniques have been successful they are time-consuming. Alternatively, frequency-domain electromagnetic (FEM) induction is potentially useful as physical contact with the ground is not required. A DUALEM-421 and EM4Soil inversion software package are used to develop a quasi two- (2D) and quasi three-dimensional (3D) electromagnetic conductivity images (EMCI) across Long Reef Beach located north of Sydney Harbour, New South Wales, Australia. The quasi 2D models discern: the dry sand (<10 mS/m) associated with the incipient dune; sand with fresh water (10 to 20 mS/m); mixing of fresh and saline water (20 to 500 mS/m), and; saline sand of varying moisture (more than 500 mS/m). The quasi 3D EMCIs generated for low and high tides suggest that daily tidal cycles do not have a significant effect on local groundwater salinity. Instead, the saline intrusion is most likely influenced by medium-large scale drivers including local wave climate and morphology along this wave-dominated beach. Further research is required to elucidate the influence of spring-neap tidal cycles, contrasting beach morphological states and sea level rise.