Combining Laser-Induced Breakdown Spectroscopy and Visible Near-Infrared Spectroscopy for Predicting Soil Organic Carbon and Texture: A Danish National-Scale Study.

Laser-induced breakdown spectroscopy (LIBS) and visible near-infrared spectroscopy (vis-NIRS) are spectroscopic techniques that offer promising alternatives to traditional laboratory methods for the rapid and cost-effective determination of soil properties on a large scale. Despite their individual limitations, combining LIBS and vis-NIRS has been shown to enhance the prediction accuracy for the determination of soil properties compared to single-sensor approaches. In this study, we used a comprehensive Danish national-scale soil dataset encompassing mostly sandy soils collected from various land uses and soil depths to evaluate the performance of LIBS and vis-NIRS, as well as their combined spectra, in predicting soil organic carbon (SOC) and texture. Firstly, partial least squares regression (PLSR) models were developed to correlate both LIBS and vis-NIRS spectra with the reference data. Subsequently, we merged LIBS and vis-NIRS data and developed PLSR models for the combined spectra. Finally, interval partial least squares regression (iPLSR) models were applied to assess the impact of variable selection on prediction accuracy for both LIBS and vis-NIRS. Despite being fundamentally different techniques, LIBS and vis-NIRS displayed comparable prediction performance for the investigated soil properties. LIBS achieved a root mean square error of prediction (RMSEP) of <7% for texture and 0.5% for SOC, while vis-NIRS achieved an RMSEP of <8% for texture and 0.5% for SOC. Combining LIBS and vis-NIRS spectra improved the prediction accuracy by 16% for clay, 6% for silt and sand, and 2% for SOC compared to single-sensor LIBS predictions. On the other hand, vis-NIRS single-sensor predictions were improved by 10% for clay, 17% for silt, 16% for sand, and 4% for SOC. Furthermore, applying iPLSR for variable selection improved prediction accuracy for both LIBS and vis-NIRS. Compared to LIBS PLSR predictions, iPLSR achieved reductions of 27% and 17% in RMSEP for clay and sand prediction, respectively, and an 8% reduction for silt and SOC prediction. Similarly, vis-NIRS iPLSR models demonstrated reductions of 6% and 4% in RMSEP for clay and SOC, respectively, and a 3% reduction for silt and sand. Interestingly, LIBS iPLSR models outperformed combined LIBS-vis-NIRS models in terms of prediction accuracy. Although combining LIBS and vis-NIRS improved the prediction accuracy of texture and SOC, LIBS coupled with variable selection had a greater benefit in terms of prediction accuracy. Future studies should investigate the influence of reference method uncertainty on prediction accuracy.

A systematic benchmarking framework for future assessments of soil health: An example from Denmark.

Based on current evidence and established critical thresholds for soil degradation indicators, it is concerning that over 60-70% of European soils are unhealthy due to unsustainable management and the impact of climate change. Despite European and national efforts to improve soil health, significant gaps remain. The proposal for a Soil Monitoring and Resilience Law, to be implemented by the European Union, seeks to establish a framework for soil monitoring and promote sustainable management practices to achieve healthy soils by 2050. This requires extensive data collection and soil monitoring systems to accurately estimate soil health across Europe, considering the diversity of soil types, climates, and land uses. To establish a framework for soil monitoring, we must understand the site-specific status of soil and the ranges of soil health indicators across specific pedoclimatic regions. In our study, we evaluated the soil status in agricultural areas in Denmark using soil health indicators and a site-specific benchmarking approach. We compiled nationally representative datasets, combining point and model-informed data of soil parameters such as organic carbon content, bulk density, pH, electrical conductivity, clay-to-soil organiccarbon ratio, water erosion, and nitrogen leaching. By categorizing Danish agricultural soils into monitoring units based on textural classes, landscape elements, and wetland types, we calculated benchmarks for these indicators, considering different cropping systems. Our approach provided detailed point-based results and a spatially explicit overview of the status of soil health indicators in Denmark. We identified areas where soil deviates from the benchmarks of different indicators. Such deviations might indicate soil functions operating outside the normal range, posing potential threats to soil health. This proposed framework could support the establishment of a baseline for assessing the directionality of future changes in soil health. Moreover, it is adaptable for implementation by other countries to support assessments of soil health.

Unleashing the sequestration potential of soil organic carbon under climate and land use change scenarios in Danish agroecosystems.

Future global climate changes are expected to increase soil organic carbon (SOC) decomposition. However, the combined effect of C inputs, land use changes, and climate on SOC turnover is still unclear. Exploring this SOC-climate-land use interaction allows us to understand the SOC stabilization mechanisms and examine whether the soil can act as a source or a sink for CO2. The current study estimates the SOC sequestration potential in the topsoil layer of Danish agricultural lands by 2038, considering the effect of land use change and future climate scenarios using the Rothamsted Carbon (RothC) model. Additionally, we quantified the loss vulnerability of existing and projected SOC based on the soil capacity to stabilize OC. We used the quantile random forest model to estimate the initial SOC stock by 2018, and we simulated the SOC sequestration potential with RothC for a business-as-usual (BAU) scenario and a crop rotation change (LUC) scenario under climate change conditions by 2038. We compared the projected SOC stocks with the carbon saturation deficit. The initial SOC stock ranged from 10 to 181 Mg C ha-1 in different parts of the country. The projections showed a SOC loss of 8.1 Mg C ha-1 for the BAU scenario and 6 Mg C ha-1 after the LUC adoption. This SOC loss was strongly influenced by warmer temperatures and clay content. The proposed crop rotation became a mitigation measure against the negative effect of climate change on SOC accumulation, especially in sandy soils with a high livestock density. A high C accumulation in C-saturated soils suggests an increase in non-complexed SOC, which is vulnerable to being lost into the atmosphere as CO2. With these results, we provide information to prioritize areas where different soil management practices can be adopted to enhance SOC sequestration in stable forms and preserve the labile-existing SOC stocks.

Improved description of terrestrial habitat types by including microbial communities as indicators.

Soils host diverse communities of microorganisms essential for ecosystem functions and soil health. Despite their importance, microorganisms are not covered by legislation protecting biodiversity or habitats, such as the Habitats Directive. Advances in molecular methods have caused breakthroughs in microbial community analysis, and recent studies have shown that parts of the communities are habitat-specific. If distinct microbial communities are present in the habitat types defined in the Habitats Directive, the Directive may be improved by including these communities. Thus, monitoring and reporting of biodiversity and conservation status of habitat types could be based not only on plant communities but also on microbial communities. In the present study, bacterial and plant communities were examined in six habitat types defined in the Habitats Directive by conducting botanical surveys and collecting soil samples for amplicon sequencing across 19 sites in Denmark. Furthermore, selected physico-chemical properties expected to differ between habitat types and explain variations in community composition of bacteria and vegetation were analysed (pH, electrical conductivity (EC), soil texture, soil water repellency, soil organic carbon content (OC), inorganic nitrogen, and in-situ water content (SWC)). Despite some variations within the same habitat type and overlaps between habitat types, habitat-specific communities were observed for both bacterial and plant communities, but no correlation was observed between the alpha diversity of vegetation and bacteria. PERMANOVA analysis was used to evaluate the variables best able to explain variation in the community composition of vegetation and bacteria. Habitat type alone could explain 46% and 47% of the variation in bacterial and plant communities, respectively. Excluding habitat type as a variable, the best model (pH, SWC, OC, fine silt, and Shannon's diversity index for vegetation) could explain 37% of the variation for bacteria. For vegetation, the best model (pH, EC, ammonium content and Shannon's diversity index for bacteria) could explain 25% of the variation. Based on these results, bacterial communities could be included in the Habitats Directive to improve the monitoring, as microorganisms are more sensitive to changes in the environment compared to vegetation, which the current monitoring is based on.

Mapping of Agricultural Subsurface Drainage Systems Using a Frequency-Domain Ground Penetrating Radar and Evaluating Its Performance Using a Single-Frequency Multi-Receiver Electromagnetic Induction Instrument.

Subsurface drainage systems are commonly used to remove surplus water from the soil profile of a poorly drained farmland. Traditional methods for drainage mapping involve the use of tile probes and trenching equipment that are time-consuming, labor-intensive, and invasive, thereby entailing an inherent risk of damaging the drainpipes. Effective and efficient methods are needed in order to map the buried drain lines: (1) to comprehend the processes of leaching and offsite release of nutrients and pesticides and (2) for the installation of a new set of drain lines between the old ones to enhance the soil water removal. Non-invasive geophysical soil sensors provide a potential alternative solution. Previous research has mainly showcased the use of time-domain ground penetrating radar, with variable success, depending on local soil and hydrological conditions and the central frequency of the specific equipment used. The objectives of this study were: (1) to test the use of a stepped-frequency continuous wave three-dimensional ground penetrating radar (3D-GPR) with a wide antenna array for subsurface drainage mapping and (2) to evaluate its performance with the use of a single-frequency multi-receiver electromagnetic induction (EMI) sensor in-combination. This sensor combination was evaluated on twelve different study sites with various soil types with textures ranging from sand to clay till. While the 3D-GPR showed a high success rate in finding the drainpipes at five sites (sandy, sandy loam, loamy sand, and organic topsoils), the results at the other seven sites were less successful due to the limited penetration depth of the 3D-GPR signal. The results suggest that the electrical conductivity estimates produced by the inversion of apparent electrical conductivity data measured by the EMI sensor could be a useful proxy for explaining the success achieved by the 3D-GPR in finding the drain lines.

Visible-Near-Infrared Spectroscopy can predict Mass Transport of Dissolved Chemicals through Intact Soil.

The intensification of agricultural production to meet the growing demand for agricultural commodities is increasing the use of chemicals. The ability of soils to transport dissolved chemicals depends on both the soil's texture and structure. Assessment of the transport of dissolved chemicals (solutes) through soils is performed using breakthrough curves (BTCs) where the application of a solute at one site and its appearance over time at another are recorded. Obtaining BTCs from laboratory studies is extremely expensive and time- and labour-consuming. Visible-near-infrared (vis-NIR) spectroscopy is well recognized for its measurement speed and for its low data acquisition cost and can be used for quantitative estimation of basic soil properties such as clay and organic matter. In this study, for the first time ever, vis-NIR spectroscopy was used to predict dissolved chemical breakthrough curves obtained from tritium transport experiments on a large variety of intact soil columns. Averaged across the field, BTCs were estimated with a high degree of accuracy. So, with vis-NIR spectroscopy, the mass transport of dissolved chemicals can be measured, paving the way for next-generation measurements and monitoring of dissolved chemical transport by spectroscopy.

Soil legacy data rescue via GlobalSoilMap and other international and national initiatives.

Legacy soil data have been produced over 70 years in nearly all countries of the world. Unfortunately, data, information and knowledge are still currently fragmented and at risk of getting lost if they remain in a paper format. To process this legacy data into consistent, spatially explicit and continuous global soil information, data are being rescued and compiled into databases. Thousands of soil survey reports and maps have been scanned and made available online. The soil profile data reported by these data sources have been captured and compiled into databases. The total number of soil profiles rescued in the selected countries is about 800,000. Currently, data for 117, 000 profiles are compiled and harmonized according to GlobalSoilMap specifications in a world level database (WoSIS). The results presented at the country level are likely to be an underestimate. The majority of soil data is still not rescued and this effort should be pursued. The data have been used to produce soil property maps. We discuss the pro and cons of top-down and bottom-up approaches to produce such maps and we stress their complementarity. We give examples of success stories. The first global soil property maps using rescued data were produced by a top-down approach and were released at a limited resolution of 1km in 2014, followed by an update at a resolution of 250m in 2017. By the end of 2020, we aim to deliver the first worldwide product that fully meets the GlobalSoilMap specifications.

Water and solute transport in agricultural soils predicted by volumetric clay and silt contents.

Solute transport through the soil matrix is non-uniform and greatly affected by soil texture, soil structure, and macropore networks. Attempts have been made in previous studies to use infiltration experiments to identify the degree of preferential flow, but these attempts have often been based on small datasets or data collected from literature with differing initial and boundary conditions. This study examined the relationship between tracer breakthrough characteristics, soil hydraulic properties, and basic soil properties. From six agricultural fields in Denmark, 193 intact surface soil columns 20cm in height and 20cm in diameter were collected. The soils exhibited a wide range in texture, with clay and organic carbon (OC) contents ranging from 0.03 to 0.41 and 0.01 to 0.08kgkg(-1), respectively. All experiments were carried out under the same initial and boundary conditions using tritium as a conservative tracer. The breakthrough characteristics ranged from being near normally distributed to gradually skewed to the right along with an increase in the content of the mineral fines (particles ≤50µm). The results showed that the mineral fines content was strongly correlated to functional soil structure and the derived tracer breakthrough curves (BTCs), whereas the OC content appeared less important for the shape of the BTC. Organic carbon was believed to support the stability of the soil structure rather than the actual formation of macropores causing preferential flow. The arrival times of 5% and up to 50% of the tracer mass were found to be strongly correlated with volumetric fines content. Predicted tracer concentration breakthrough points as a function of time up to 50% of applied tracer mass could be well fitted to an analytical solution to the classical advection-dispersion equation. Both cumulative tracer mass and concentration as a function of time were well predicted from the simple inputs of bulk density, clay and silt contents, and applied tracer mass. The new concept seems promising as a platform towards more accurate proxy functions for dissolved contaminant transport in intact soil.

Modeling Soil Organic Carbon at Regional Scale by Combining Multi-Spectral Images with Laboratory Spectra.

There is a great challenge in combining soil proximal spectra and remote sensing spectra to improve the accuracy of soil organic carbon (SOC) models. This is primarily because mixing of spectral data from different sources and technologies to improve soil models is still in its infancy. The first objective of this study was to integrate information of SOC derived from visible near-infrared reflectance (Vis-NIR) spectra in the laboratory with remote sensing (RS) images to improve predictions of topsoil SOC in the Skjern river catchment, Denmark. The second objective was to improve SOC prediction results by separately modeling uplands and wetlands. A total of 328 topsoil samples were collected and analyzed for SOC. Satellite Pour l'Observation de la Terre (SPOT5), Landsat Data Continuity Mission (Landsat 8) images, laboratory Vis-NIR and other ancillary environmental data including terrain parameters and soil maps were compiled to predict topsoil SOC using Cubist regression and Bayesian kriging. The results showed that the model developed from RS data, ancillary environmental data and laboratory spectral data yielded a lower root mean square error (RMSE) (2.8%) and higher R2 (0.59) than the model developed from only RS data and ancillary environmental data (RMSE: 3.6%, R2: 0.46). Plant-available water (PAW) was the most important predictor for all the models because of its close relationship with soil organic matter content. Moreover, vegetation indices, such as the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI), were very important predictors in SOC spatial models. Furthermore, the 'upland model' was able to more accurately predict SOC compared with the 'upland & wetland model'. However, the separately calibrated 'upland and wetland model' did not improve the prediction accuracy for wetland sites, since it was not possible to adequately discriminate the vegetation in the RS summer images. We conclude that laboratory Vis-NIR spectroscopy adds critical information that significantly improves the prediction accuracy of SOC compared to using RS data alone. We recommend the incorporation of laboratory spectra with RS data and other environmental data to improve soil spatial modeling and digital soil mapping (DSM).

Digital mapping of soil organic carbon contents and stocks in Denmark.

Estimation of carbon contents and stocks are important for carbon sequestration, greenhouse gas emissions and national carbon balance inventories. For Denmark, we modeled the vertical distribution of soil organic carbon (SOC) and bulk density, and mapped its spatial distribution at five standard soil depth intervals (0-5, 5-15, 15-30, 30-60 and 60-100 cm) using 18 environmental variables as predictors. SOC distribution was influenced by precipitation, land use, soil type, wetland, elevation, wetness index, and multi-resolution index of valley bottom flatness. The highest average SOC content of 20 g kg(-1) was reported for 0-5 cm soil, whereas there was on average 2.2 g SOC kg(-1) at 60-100 cm depth. For SOC and bulk density prediction precision decreased with soil depth, and a standard error of 2.8 g kg(-1) was found at 60-100 cm soil depth. Average SOC stock for 0-30 cm was 72 t ha(-1) and in the top 1 m there was 120 t SOC ha(-1). In total, the soils stored approximately 570 Tg C within the top 1 m. The soils under agriculture had the highest amount of carbon (444 Tg) followed by forest and semi-natural vegetation that contributed 11% of the total SOC stock. More than 60% of the total SOC stock was present in Podzols and Luvisols. Compared to previous estimates, our approach is more reliable as we adopted a robust quantification technique and mapped the spatial distribution of SOC stock and prediction uncertainty. The estimation was validated using common statistical indices and the data and high-resolution maps could be used for future soil carbon assessment and inventories.

Comparative mapping of soil physical-chemical and structural parameters at field scale to identify zones of enhanced leaching risk.

Preferential flow and particle-facilitated transport through macropores contributes significantly to the transport of strongly sorbing substances such as pesticides and phosphorus. The aim of this study was to perform a field-scale characterization of basic soil physical properties like clay and organic carbon content and investigate whether it was possible to relate these to derived structural parameters such as bulk density and conservative tracer parameters and to actual particle and phosphorus leaching patterns obtained from laboratory leaching experiments. Sixty-five cylindrical soil columns of 20-cm height and 20-cm diameter and bulk soil were sampled from the topsoil in a 15-m × 15-m grid in an agricultural loamy field. Highest clay contents and highest bulk densities were found in the northern part of the field. Leaching experiments with a conservative tracer showed fast 5% tracer arrival times and high tracer recovery percentages from columns sampled from the northern part of the field, and the leached mass of particles and particulate phosphorus was also largest from this area. Strong correlations were obtained between 5% tracer arrival time, tracer recovery, and bulk density, indicating that a few well-aligned and better connected macropores might change the hydraulic conductivity between the macropores and the soil matrix, triggering an onset of preferential flow at lower rain intensities compared with less compacted soil. Overall, a comparison mapping of basic and structural characteristics including soil texture, bulk density, dissolved tracer, particle and phosphorus transport parameters identified the northern one-third of the field as a zone with higher leaching risk. This risk assessment based on parameter mapping from measurements on intact samples was in good agreement with 9 yr of pesticide detections in two horizontal wells and with particle and phosphorus leaching patterns from a distributed, shallow drainage pipe system across the field.

Predictive mapping of soil organic carbon in wet cultivated lands using classification-tree based models: the case study of Denmark.

Soil organic carbon (SOC) is one of the most important carbon stocks globally and has large potential to affect global climate. Distribution patterns of SOC in Denmark constitute a nation-wide baseline for studies on soil carbon changes (with respect to Kyoto protocol). This paper predicts and maps the geographic distribution of SOC across Denmark using remote sensing (RS), geographic information systems (GISs) and decision-tree modeling (un-pruned and pruned classification trees). Seventeen parameters, i.e. parent material, soil type, landscape type, elevation, slope gradient, slope aspect, mean curvature, plan curvature, profile curvature, flow accumulation, specific catchment area, tangent slope, tangent curvature, steady-state wetness index, Normalized Difference Vegetation Index (NDVI), Normalized Difference Wetness Index (NDWI) and Soil Color Index (SCI) were generated to statistically explain SOC field measurements in the area of interest (Denmark). A large number of tree-based classification models (588) were developed using (i) all of the parameters, (ii) all Digital Elevation Model (DEM) parameters only, (iii) the primary DEM parameters only, (iv), the remote sensing (RS) indices only, (v) selected pairs of parameters, (vi) soil type, parent material and landscape type only, and (vii) the parameters having a high impact on SOC distribution in built pruned trees. The best constructed classification tree models (in the number of three) with the lowest misclassification error (ME) and the lowest number of nodes (N) as well are: (i) the tree (T1) combining all of the parameters (ME=29.5%; N=54); (ii) the tree (T2) based on the parent material, soil type and landscape type (ME=31.5%; N=14); and (iii) the tree (T3) constructed using parent material, soil type, landscape type, elevation, tangent slope and SCI (ME=30%; N=39). The produced SOC maps at 1:50,000 cartographic scale using these trees are highly matching with coincidence values equal to 90.5% (Map T1/Map T2), 95% (Map T1/Map T3) and 91% (Map T2/Map T3). The overall accuracies of these maps once compared with field observations were estimated to be 69.54% (Map T1), 68.87% (Map T2) and 69.41% (Map T3). The proposed tree models are relatively simple, and may be also applied to other areas.

Spatial soil zinc content distribution from terrain parameters: a GIS-based decision-tree model in Lebanon.

Heavy metal contamination has been and continues to be a worldwide phenomenon that has attracted a great deal of attention from governments and regulatory bodies. In this context, our study proposes a regression-tree model to predict the concentration level of zinc in the soils of northern Lebanon (as a case study of Mediterranean landscapes) under a GIS environment. The developed tree-model explained 88% of variance in zinc concentration using pH (100% in relative importance), surroundings of waste areas (90%), proximity to roads (80%), nearness to cities (50%), distance to drainage line (25%), lithology (24%), land cover/use (14%), slope gradient (10%), conductivity (7%), soil type (7%), organic matter (5%), and soil depth (5%). The overall accuracy of the quantitative zinc map produced (at 1:50.000 scale) was estimated to be 78%. The proposed tree model is relatively simple and may also be applied to other areas.