Search

Invasive plant species (IPS) pose a significant threat to natural ecosystems, causing substantial global biodiversity loss. Despite the fact that only a few species invade Israel's coastal dunes, their impact on landscape and ecosystems is notable. The invasive plant Heterotheca subaxillaris, introduced a few decades ago, has aggressively spread across these dunes, becoming a major invasive species along the Israeli coast. To provide essential spatial information for ecologists and land managers, this study developed a methodology to identify and map H. subaxillaris at various scales by integrating unmanned aerial vehicle (UAV) and satellite imagery. During the flowering period, drone images were employed to create multiresolution ground references using the random forest classification method, achieving 97.4% accuracy for H. subaxillaris identification. These UAV-based references validated mapping results from Vision-1 and PlanetScope satellite images, with accuracies of 96.9% and 83.3%, respectively. Vision-1 image outcomes were used as a reference for extensive probabilistic mapping using Sentinel-2 images. Given the species' small size and dispersed nature, no previous study has utilized the synergy of UAV and satellite images for multiscale mapping of H. subaxillaris. This research significantly advances detection across varying spatial resolutions, serving as a foundation for future monitoring systems with regularly acquired data. In addition, the study offers an annual phenology-based approach for identifying flowering dates through spectral indices, contributing to phenological change investigations. The present method also helps remote sensing track IPS expansion in similar species with sparse canopy cover, providing vital geographic data for conservation and management to prevent future growth.

Diffuse reflectance spectroscopy has been extensively employed to deliver timely and cost-effective predictions of a number of soil properties. However, although several soil spectral laboratories have been established worldwide, the distinct characteristics of instruments and operations still hamper further integration and interoperability across mid-infrared (MIR) soil spectral libraries. In this study, we conducted a large-scale ring trial experiment to understand the lab-to-lab variability of multiple MIR instruments. By developing a systematic evaluation of different mathematical treatments with modeling algorithms, including regular preprocessing and spectral standardization, we quantified and evaluated instruments' dissimilarity and how this impacts internal and shared model performance. We found that all instruments delivered good predictions when calibrated internally using the same instruments' characteristics and standard operating procedures by solely relying on regular spectral preprocessing that accounts for light scattering and multiplicative/additive effects, e.g., using standard normal variate (SNV). When performing model transfer from a large public library (the USDA NSSC-KSSL MIR library) to secondary instruments, good performance was also achieved by regular preprocessing (e.g., SNV) if both instruments shared the same manufacturer. However, significant differences between the KSSL MIR library and contrasting ring trial instruments responses were evident and confirmed by a semi-unsupervised spectral clustering. For heavily contrasting setups, spectral standardization was necessary before transferring prediction models. Non-linear model types like Cubist and memory-based learning delivered more precise estimates because they seemed to be less sensitive to spectral variations than global partial least square regression. In summary, the results from this study can assist new laboratories in building spectroscopy capacity utilizing existing MIR spectral libraries and support the recent global efforts to make soil spectroscopy universally accessible with centralized or shared operating procedures.

Forest soils have large contents of carbon (C) and total nitrogen (TN), which have significant spatial variability laterally across landscapes and vertically with depth due to decomposition, erosion and leaching. Therefore, the ratio of C to TN contents (C:N), a crucial indicator of soil quality and health, is also different depending on soil horizon. These attributes can cost-effectively and rapidly be estimated using visible–near infrared–shortwave infrared (VNIR–SWIR) spectroscopy. Nevertheless, the effect of different soil layers, particularly over large scales of highly heterogeneous forest soils, on the performance of the technique has rarely been attempted. This study evaluated the potential of VNIR–SWIR spectroscopy in quantification and variability analysis of C:N in soils from different organic and mineral layers of forested sites of the Czech Republic. At each site, we collected samples from the litter (L), fragmented (F) and humus (H) organic layers, and from the A1 (depth of 2–10 cm) and A2 (depth of 10–40 cm) mineral layers providing a total of 2505 samples. Support vector machine regression (SVMR) was used to train the prediction models of the selected attributes at each individual soil layer and the merged layer (profile). We further produced the spatial distribution maps of C:N as the target attribute at each soil layer. Results showed that the prediction accuracy based on the profile spectral data was adequate for all attributes. Moreover, F was the most accurately predicted layer, regardless of the soil attribute. C:N models and maps in the organic layers performed well although in mineral layers, models were poor and maps were reliable only in areas with low and moderate C:N. On the other hand, the study indicated that reflectance spectra could efficiently predict and map organic layers of the forested sites. Although, in mineral layers, high values of C:N (≥ 50) were not detectable in the map created based on the reflectance spectra. In general, the study suggests that VNIR–SWIR spectroscopy has the feasibility of modelling and mapping C:N in soil organic horizons based on national spectral data in the forests of the Czech Republic.

Mapping soil organic carbon (SOC) stock can serve as a resilience indicator for climate change. As part of the carbon dioxide (CO2) sink, soil has recently become an integral part of the global carbon agenda to mitigate climate change. We used hyperspectral remote sensing to model the SOC stock in the Sele River plain located in the Campania region in southern Italy. To this end, a soil spectral library (SSL) for the Campania region was combined with an aerial hyperspectral image acquired with the AVIRIS–NG sensor mounted on a Twin Otter aircraft at an altitude of 1433 m. The products of this study were four raster layers with a high spatial resolution (1 m), representing the SOC stocks and three other related soil attributes: SOC content, clay content, and bulk density (BD). We found that the clay minerals’ spectral absorption at 2200 nm has a significant impact on predicting the examined soil attributes. The predictions were performed by using AVIRIS–NG sensor data over a selected plot and generating a quantitative map which was validated with in situ observations showing high accuracies in the ground-truth stage (OC stocks [RPIQ = 2.19, R2 = 0.72, RMSE = 0.07]; OC content [RPIQ = 2.27, R2 = 0.80, RMSE = 1.78]; clay content [RPIQ = 1.6 R2 = 0.89, RMSE = 25.42]; bulk density [RPIQ = 1.97, R2 = 0.84, RMSE = 0.08]). The results demonstrated the potential of combining SSLs with remote sensing data of high spectral/spatial resolution to estimate soil attributes, including SOC stocks.

This study explores the potential of using smartphone cameras to estimate soil attributes as a field-based alternative to traditional sampling and lab analyses that are time-consuming, expensive, and non-friendly to the environment. It aims to investigate whether the color information captured by smartphone cameras could provide spectral information that can be utilized for soil attribute estimation using regression techniques. The DN values from smartphone camera photos were transformed into reflectance using colored plastic plated marker (PPM) reference panels. The spectral information of the PPM keys, obtained using a laboratory spectrometer, was convolved into 3 RGB smartphone bands. Using the legacy soil spectral library (SSL) of Israel on the 3-reflectance calibrated RGB smartphone bands, a good agreement with the reflectance of the laboratory spectrometer was obtained. The soil properties were modeled with the converted reflectance using the multiple linear regression (MLR), partial least square regression (PLSR), and ridge regression (RR) algorithms. The models created using the smartphone data (RGB and 6 synthetic bands) predicted six soil properties: i.e., soil organic matter (OM%), CaCO3, free iron oxides (Fe-dithionite), Al2O3, soil surface area (SSA), and hygroscopic moisture with the coefficient of determination (R2) values of 0.58, 0.54, 0.60, 0.55, 0.55 and 0.64, respectively. The quality of predictions was also evaluated using root mean square error (RMSE) and the ratio of prediction to interquartile distance (RPIQ). This research provides the results of an experiment to convert a personal smartphone camera to a spectrometer with the help of SSL. Despite the low spectral resolution of smartphone cameras, the models achieved fair results in predicting six soil properties. These findings could promote the adoption of citizen science methodologies for sustainable soil management as smartphone devices that are available to all can be simply calibrated to reflectance and used for proximal sensing of soils.

Soil spectral libraries (SSLs) are important big-data archives (spectra associated with soil properties) that are analyzed via machine-learning algorithms to estimate soil attributes. Since different spectral measurement protocols are applied when constructing SSLs, it is necessary to examine harmonization techniques to merge the data. In recent years, several techniques for harmonization have been proposed, among which the internal soil standard (ISS) protocol is the most largely applied and has demonstrated its capacity to rectify systematic effects during spectral measurements. Here, we postulate that a spectral transfer function (TF) can be extracted between existing (old) SSLs if a subset of samples from two (or more) different SSLs are remeasured using the ISS protocol. A machine-learning TF strategy was developed, assembling random forest (RF) spectral-based models to predict the ISS spectral condition using soil samples from two existing SSLs. These SSLs had already been measured using different protocols without any ISS treatment the Brazilian (BSSL, generated in 2019) and the European (LUCAS, generated in 2009–2012) SSLs. To verify the TF’s ability to improve the spectral assessment of soil attributes after harmonizing the different SSLs’ protocols, RF spectral-based models for estimating organic carbon (OC) in soil were developed. The results showed high spectral similarities between the ISS and the ISS–TF spectral observations, indicating that post-ISS rectification is possible. Furthermore, after merging the SSLs with the TFs, the spectral-based assessment of OC was considerably improved, from R2 = 0.61, RMSE (g/kg) = 12.46 to R2 = 0.69, RMSE (g/kg) = 11.13. Given our results, this paper enhances the importance of soil spectroscopy by contributing to analyses in remote sensing, soil surveys, and digital soil mapping.

The hyperspectral (HSR) sensors Earth Surface Mineral Dust Source Investigation (EMIT) of the National Aeronautics and Space Administration (NASA) and Environmental Mapping and Analysis Program (EnMAP) of the German Aerospace Center (DLR) were recently launched. These state-of-the-art sensors have joined the already operational HSR sensors DESIS (DLR), PRISMA (Italian Space Agency), and HISUI (developed by the Japanese Ministry of Economy, Trade, and Industry METI and Japan Aerospace Exploration Agency JAXA). The launching of more HSR sensors is being planned for the near future (e.g., SBG of NASA, and CHIME of the European Space Agency), and the challenge of monitoring and maintaining their calibration accuracy is becoming more relevant. We proposed two test sites: Amiaz Plain (AP) and Makhtesh Ramon (MR) for spectral, radiometric, and geometric calibration/validation (CAL/VAL). The sites are situated in the arid environment of southern Israel and are in the same overpass coverage. Both test sites have already demonstrated favorable results in assessing an HSR sensor’s performance and were chosen to participate in the EMIT and EnMAP validation stage. We first evaluated the feasibility of using AP and MR as CAL/VAL test sites with extensive datasets and sensors, such as the multispectral sensor Landsat (Landsat5 TM and Landsat8 OLI), the airborne HSR sensor AisaFENIX 1K, and the spaceborne HSR sensors DESIS and PRISMA. Field measurements were taken over time. The suggested methodology integrates reflectance and radiometric CAL/VAL test sites into one operational protocol. The method can highlight degradation in the spectral domain early on, help maintain quantitative applications, adjust the sensor’s radiometric calibration during its mission lifetime, and minimize uncertainties of calibration parameters. A PRISMA sensor case study demonstrates the complete operational protocol, i.e., performance evaluation, quality assessment, and cross-calibration between HSR sensors. These CAL/VAL sites are ready to serve as operational sites for other HSR sensors.

Soils are an essential factor contributing to the agricultural production of rainfed crops such as barley and triticale cereals. Changing environmental conditions and inadequate land management are endangering soil quality and productivity and, in turn, crop quality and productivity are affected. Advances in hyperspectral remote sensing are of great use for the spatial characterization and monitoring of the soil degradation status, as well as its impact on crop growth and agricultural productivity. In this study, hyperspectral airborne data covering the visible, near-infrared, short-wave infrared, and thermal infrared (VNIR–SWIR–TIR, 0.4–12 µm) were acquired in a Mediterranean agricultural area of central Spain and used to analyze the spatial differences in vegetation vitality and grain yield in relation to the soil degradation status. Specifically, leaf area index (LAI), crop water stress index (CWSI), and the biomass of the crop yield are derived from the remote sensing data and discussed regarding their spatial differences and relationship to a classification of erosion and accumulation stages (SEAS) based on previous remote sensing analyses during bare soil conditions. LAI and harvested crop biomass yield could be well estimated by PLS regression based on the hyperspectral and in situ reference data (R2 of 0.83, r of 0.91, and an RMSE of 0.2 m2 m−2 for LAI and an R2 of 0.85, r of 0.92, and an RMSE of 0.48 t ha−1 for grain yield). In addition, the soil erosion and accumulation stages (SEAS) were successfully predicted based on the canopy spectral signal of vegetated crop fields using a random forest machine learning approach. Overall accuracy was achieved above 71% by combining the VNIR–SWIR–TIR canopy reflectance and emissivity of the growing season with topographic information after reducing the redundancy in the spectral dataset. The results show that the estimated crop traits are spatially related to the soil’s degradation status, with shallow and highly eroded soils, as well as sandy accumulation zones being associated with areas of low LAI, crop yield, and high crop water stress. Overall, the results of this study illustrate the enormous potential of imaging spectroscopy for a combined analysis of the plant-soil system in the frame of land and soil degradation monitoring.

We conducted a systematic review and inventory of recent research achievements related to spaceborne and aerial Earth Observation (EO) data-driven monitoring in support of soil-related strategic goals for a three-year period (2019–2021). Scaling, resolution, data characteristics, and modelling approaches were summarized, after reviewing 46 peer-reviewed articles in international journals. Inherent limitations associated with an EO-based soil mapping approach that hinder its wider adoption were recognized and divided into four categories: (i) area covered and data to be shared; (ii) thresholds for bare soil detection; (iii) soil surface conditions; and (iv) infrastructure capabilities. Accordingly, we tried to redefine the meaning of what is expected in the next years for EO data-driven topsoil monitoring by performing a thorough analysis driven by the upcoming technological waves. The review concludes that the best practices for the advancement of an EO data-driven soil mapping include: (i) a further leverage of recent artificial intelligence techniques to achieve the desired representativeness and reliability; (ii) a continued effort to share harmonized labelled datasets; (iii) data fusion with in situ sensing systems; (iv) a continued effort to overcome the current limitations in terms of sensor resolution and processing limitations of this wealth of EO data; and (v) political and administrative issues (e.g., funding, sustainability). This paper may help to pave the way for further interdisciplinary research and multi-actor coordination activities and to generate EO-based benefits for policy and economy.

Marine sponges harbor a diverse array of microorganisms and the composition of the microbial community has been suggested to be linked to holo-biont health. Most of the attention concerning sponge mycobiomes has been given to sponges present in shallow depths. Here, we describe the presence of 146 culturable mycobiome taxa isolated from mesophotic niche (100 m depth)-inhabiting samples of Agelas oroides, in the Mediterranean Sea. We identify some potential in vitro interactions between several A. oroides-associated fungi and show that sponge meso-hyl extract, but not its predominantly collagen-rich part, is sufficient to support hyphal growth. We demonstrate that changes in the diversity of culturable mycobiome constituents occur following sponge transplantation from its original mesophotic habitat to shallow (10 m) waters, where historically (60 years ago) this species was found. We conclude that among the 30 fungal genera identified as associated with A. oroides, Aspergillus, Penicillium and Trichoderma constitute the core mycobiome of A. oroides, and that they persist even when the sponge is transplanted to a suboptimal environment, indicative of the presence of constant, as well as dynamic, components of the sponge mycobiome. Other genera seemed more depth-related and appeared or disappeared upon host’s transfer from 100 to 10 m.

Water infiltration rate (WIR) into the soil profile was investigated through a comprehensive study harnessing spectral information of the soil surface. As soil spectroscopy provides invaluable information on soil attributes, and as WIR is a soil surface-dependent property, field spectroscopy may model WIR better than traditional laboratory spectral measurements. This is because sampling for the latter disrupts the soil-surface status. A field soil spectral library (FSSL), consisting of 114 samples with different textures from six different sites over the Mediterranean basin, combined with traditional laboratory spectral measurements, was created. Next, partial least squares regression analysis was conducted on the spectral and WIR data in different soil texture groups, showing better performance of the field spectral observations compared to traditional laboratory spectroscopy. Moreover, several quantitative spectral properties were lost due to the sampling procedure, and separating the samples according to texture gave higher accuracies. Although the visible near-infrared–shortwave infrared (VNIR–SWIR) spectral region provided better accuracy, we resampled the spectral data to the resolution of a Cubert hyperspectral sensor (VNIR). This hyperspectral sensor was then assembled on an unmanned aerial vehicle (UAV) to apply one selected spectral-based model to the UAV data and map the WIR in a semi-vegetated area within the Alento catchment, Italy. Comprehensive spectral and WIR ground-truth measurements were carried out simultaneously with the UAV–Cubert sensor flight. The results were satisfactorily validated on the ground using field samples, followed by a spatial uncertainty analysis, concluding that the UAV with hyperspectral remote sensing can be used to map soil surface-related soil properties.

Soil surveys are critical for maintaining sustainable use of natural resources while minimizing harmful impacts to the ecosystem. A key soil attribute for many environmental factors, such as CO2 budget, soil fertility and sustainability, is soil organic matter (SOM), as well as its sequestration. Soil spectroscopy is a popular method to assess SOM content rapidly in both field and laboratory domains. However, SOM source composition differs from soil to soil, and the use of spectral-based models for quantifying SOM may present limited accuracy when applying a generic approach to SOM assessment. We therefore examined the extent to which the generic approach can assess SOM contents of different origin using spectral-based models. We created an artificial big dataset composed of pure dune sand as a SOM-free background, which was artificially mixed with increasing amounts of different organic matter (OM) sources obtained from commercial compost of different origins. Dune sand has high albedo and yields optimal conditions for SOM detection. This study combined two methods: partial least squares regression for the prediction of SOM content from reflectance values across the 400–2500 nm region and a soil spectral detection limit (SSDL) to judge the prediction accuracy. Spectral-based models to assess SOM content were evaluated with each OM source as well as with a merged dataset that contained all of the generated samples (generic approach). The latter was concluded to have limitations for assessing low amounts of SOM (

Symptoms of root stress are hard to detect using non-invasive tools. This study reveals proof of concept for vegetation indices’ ability, usually used to sense canopy status, to detect root stress, and performance status. Pepper plants were grown under controlled greenhouse conditions under different potassium and salinity treatments. The plants’ spectral reflectance was measured on the last day of the experiment when more than half of the plants were already naturally infected by root disease. Vegetation indices were calculated for testing the capability to distinguish between healthy and root-damaged plants using spectral measurements. While no visible symptoms were observed in the leaves, the vegetation indices and red-edge position showed clear differences between the healthy and the root-infected plants. These results were achieved after a growth period of 32 days, indicating the ability to monitor root damage at an early growing stage using leaf spectral reflectance.

Soil contamination by potentially toxic elements (PTEs) is intensifying under increasing industrialization. Thus, the ability to efficiently delineate contaminated sites is crucial. Visible–near infrared (vis–NIR: 350–2500 nm) and X-ray fluorescence (XRF: 0.02–41.08 keV) spectroscopic techniques have attracted tremendous attention for the assessment of PTEs. Recently, the application of fused vis–NIR and XRF spectroscopy, which is based on the complementary effect of data fusion, is also increasing. Moreover, different data manipulation methods, including feature selection approaches, affect the prediction performance. This study investigated the feasibility of using single and fused vis–NIR and XRF spectra while exploring feature selection algorithms for the assessment of key soil PTEs. The soil samples were collected from one of the most heavily polluted areas of the Czech Republic and scanned using laboratory vis–NIR and XRF spectrometers. Univariate filter (UF) and genetic algorithm (GA) were used to select the bands of greater importance for the PTE prediction. Support vector machine (SVM) was then used to train the models using the full-range and feature-selected spectra of single sensors and their fusion. It was found that XRF spectra alone (primarily GA-selected) performed better than single vis–NIR and fused spectral data for predictions of PTEs. Moreover, the prediction models that were derived from the fused data set (particularly the GA-selected) enhanced the models’ accuracies as compared with the single vis–NIR spectra. In general, the results suggest that the GA-selected spectra obtained from the single XRF spectrometer (for As and Pb) and from the fusion of vis–NIR and XRF (for Pb) are promising for accurate quantitative estimation detection of the mentioned PTEs.

Potassium is a macro element in plants that is typically supplied to crops in excess throughout the season to avoid a deficit leading to reduced crop yield. Transpiration rate is a momentary physiological attribute that is indicative of soil water content, the plant’s water requirements, and abiotic stress factors. In this study, two systems were combined to create a hyperspectral–physiological plant database for classification of potassium treatments (low, medium, and high) and estimation of momentary transpiration rate from hyperspectral images. PlantArray 3.0 was used to control fertigation, log ambient conditions, and calculate transpiration rates. In addition, a semi-automated platform carrying a hyperspectral camera was triggered every hour to capture images of a large array of pepper plants. The combined attributes and spectral information on an hourly basis were used to classify plants into their given potassium treatments (average accuracy = 80%) and to estimate transpiration rate (RMSE = 0.025 g/min, R2 = 0.75) using the advanced ensemble learning algorithm XGBoost (extreme gradient boosting algorithm). Although potassium has no direct spectral absorption features, the classification results demonstrated the ability to label plants according to potassium treatments based on a remotely measured hyperspectral signal. The ability to estimate transpiration rates for different potassium applications using spectral information can aid in irrigation management and crop yield optimization. These combined results are important for decision-making during the growing season, and particularly at the early stages when potassium levels can still be corrected to prevent yield loss.

The investigation of iron oxides in soil using spectral reflectance is very common. Their spectral signal is significant across the visible–near infrared (VIS–NIR) spectral range (400–1000 nm). However, this range overlaps with other soil chromophores, such as those for water and soil organic matter (SOM). This study aimed to investigate the effect of different SOM species on red soil from Israel, which is rich in hematite iron oxide, under air-dried conditions. We constructed datasets of artificially mixed soil and organic matter (OM) with different percentages of added compost from two sources (referred to as A2 and A5). Eighty subsamples of mixed soil–OM were prepared for each of the OM (compost) types. To investigate the effect of OM on the strong iron-oxide absorbance at 880 nm, we generated two indices: CRDC, the absorbance spectral depth change at 880 nm after continuous removal, and NRIR, the normalized red index ratio using 880 and 780 nm wavelengths. The different OM types influenced the soil reflectance differently. At low %SOM, up to 1.5%, the OM types behaved more similarly, but as the OM content increased, their effect on the iron-oxide signal was greater, enhancing the significant differences between the two OM sources. Moreover, as the SOM content increased, the iron-oxide signal decreased until it was completely masked out from the reflectance spectrum. The masking point was observed at different SOM contents: 4% for A5 and 8% for A2. A mechanism that explains the indirect chromophore activity of SOM in the visible region, which is related to the iron-oxide spectral features, was provided. We also compared the use of synthetic linear-mixing practices (soil–OM) to the authentic mixed samples. The synthetic mixture could not imitate the authentic soil reflectance status, especially across the overlapping spectral position of the iron oxides and OM, and hence may hinder real conditions.

A novel hyperspectral-physiological system that monitors plants dynamic response to abiotic alterations was developed. The system is a sensor-to-plant platform which can determine the optimal time of day during which physiological traits can be successfully identified via spectral means. The directly measured traits include momentary and daily transpiration rates throughout the daytime and daily and periodical plant weight loss and gain. The system monitored and evaluated pepper plants response to varying levels of potassium fertilization. Significant momentary transpiration rates differences were found between the treatments during 07:00–10:00 and 14:00–17:00. The simultaneous frequently measured high-resolution spectral data provided the means to correlate the two measured data sets. Significant correlation coefficients between the spectra and momentary transpiration rates resulted with a selection of three bands (ρ523, ρ697 and ρ818nm) that were used to capture transpiration rate differences using a normalized difference formula during the morning, noon and the afternoon. These differences also indicated that the best results are not always obtained when spectral (remote or proximal) measurements are typically preformed around noon (when solar illumination is the highest). Valuable information can be obtained when the spectral measurements are timed according to the plants’ dynamic physiological status throughout the day, which may vary among plant species and should be considered when planning remote sensing data acquisition.

Reflectance of light across the visible, near-infrared and shortwave infrared (VIS-NIR-SWIR, 0.4–2.5 µm) spectral region is very useful for investigating mineralogical, physical and chemical properties of soils, which can reduce the need for traditional wet chemistry analyses. As many collections of multispectral satellite data are available for environmental studies, a large extent with medium resolution mapping could be benefited from the spectral measurements made from remote sensors. In this paper, we explored the use of bare soil composites generated from the large historical collections of Landsat images for mapping cropland topsoil attributes across the European extent. For this task, we used the Geospatial Soil Sensing System (GEOS3) for generating two bare soil composites of 30 m resolution (named synthetic soil images, SYSI), which were employed to represent the median topsoil reflectance of bare fields. The first (framed SYSI) was made with multitemporal images (2006–2012) framed to the survey time of the Land-Use/Land-Cover Area Frame Survey (LUCAS) soil dataset (2009), seeking to be more compatible to the soil condition upon the sampling campaign. The second (full SYSI) was generated from the full collection of Landsat images (1982–2018), which although displaced to the field survey, yields a higher proportion of bare areas for soil mapping. For evaluating the two SYSIs, we used the laboratory spectral data as a reference of topsoil reflectance to calculate the Spearman correlation coefficient. Furthermore, both SYSIs employed machine learning for calibrating prediction models of clay, sand, soil organic carbon (SOC), calcium carbonates (CaCO3), cation exchange capacity (CEC), and pH determined in water, using the gradient boosting regression algorithm. The original LUCAS laboratory spectra and a version of the data resampled to the Landsat multispectral bands were also used as reference of prediction performance using VIS-NIR-SWIR multispectral data. Our results suggest that generating a bare soil composite displaced to the survey time of soil observations did not improve the quality of topsoil reflectance, and consequently, the prediction performance of soil attributes. Despite the lower spectral resolution and the variability of soils in Europe, a SYSI calculated from the full collection of Landsat images can be employed for topsoil prediction of clay and CaCO3 contents with a moderate performance (testing R2, root mean square error (RMSE) and ratio of performance to interquartile range (RPIQ) of 0.44, 9.59, 1.77, and 0.36, 13.99, 1.54, respectively). Thus, this study shows that although there exist some constraints due to the spatial and temporal variation of soil exposures and among the Landsat sensors, it is possible to use bare soil composites for mapping key soil attributes of croplands across the European extent.

Soil mineralogy can be used to study changes in the environment affecting the soil surface, such as dust from the desert through Aeolian processes, which is one of the sources that determine the mineral nature of the soil. Ground- and field-based hyperspectral longwave infrared images, acquired before and after dust dispersion on the soil surface, were processed and analyzed by applying a procedure for determining soil surface mineralogy from the emissivity spectrum, using two indices―SQCMI (the Soil Quartz Clay Mineral Index) and SCI (the Soil Carbonate Index)―to identify changes in the abundance of quartz, clay minerals and carbonates on the surface, caused by the settling dust particles. Mineralogical changes were identified, depending on the mineral composition of the dust compared to the soil surface mineralogy.

Earth observation (EO) has an immense potential as being an enabling tool for mapping spatial characteristics of the topsoil layer. Recently, deep learning based algorithms and cloud computing infrastructure have become available with a great potential to revolutionize the processing of EO data. This paper aims to present a novel EO-based soil monitoring approach leveraging open-access Copernicus Sentinel data and Google Earth Engine platform. Building on key results from existing data mining approaches to extract bare soil reflectance values the current study delivers valuable insights on the synergistic use of open access optical and radar images. The proposed framework is driven by the need to eliminate the influence of ambient factors and evaluate the efficiency of a convolutional neural network (CNN) to effectively combine the complimentary information contained in the pool of both optical and radar spectral information and those form auxiliary geographical coordinates mainly for soil. We developed and calibrated our multi-input CNN model based on soil samples (calibration = 80% and validation 20%) of the LUCAS database and then applied this approach to predict soil clay content. A promising prediction performance (R2 = 0.60, ratio of performance to the interquartile range (RPIQ) = 2.02, n = 6136) was achieved by the inclusion of both types (synthetic aperture radar (SAR) and laboratory visible near infrared–short wave infrared (VNIR-SWIR) multispectral) of observations using the CNN model, demonstrating an improvement of more than 5.5% in RMSE using the multi-year median optical composite and current state-of-the-art non linear machine learning methods such as random forest (RF; R2 = 0.55, RPIQ = 1.91, n = 6136) and artificial neural network (ANN; R2 = 0.44, RPIQ = 1.71, n = 6136). Moreover, we examined post-hoc techniques to interpret the CNN model and thus acquire an understanding of the relationships between spectral information and the soil target identified by the model. Looking to the future, the proposed approach can be adopted on the forthcoming hyperspectral orbital sensors to expand the current capabilities of the EO component by estimating more soil attributes with higher predictive performance.

With the increasing role that unmanned aerial systems (UAS) are playing in data collection for environmental studies, two key challenges relate to harmonizing and providing standardized guidance for data collection, and also establishing protocols that are applicable across a broad range of environments and conditions. In this context, a network of scientists are cooperating within the framework of the Harmonious Project to develop and promote harmonized mapping strategies and disseminate operational guidance to ensure best practice for data collection and interpretation. The culmination of these efforts is summarized in the present manuscript. Through this synthesis study, we identify the many interdependencies of each step in the collection and processing chain, and outline approaches to formalize and ensure a successful workflow and product development. Given the number of environmental conditions, constraints, and variables that could possibly be explored from UAS platforms, it is impractical to provide protocols that can be applied universally under all scenarios. However, it is possible to collate and systematically order the fragmented knowledge on UAS collection and analysis to identify the best practices that can best ensure the streamlined and rigorous development of scientific products.

Soil quality (SQ) assessment has numerous applications for managing sustainable soil function. Airborne imaging spectroscopy (IS) is an advanced tool for studying natural and artificial materials, in general, and soil properties, in particular. The primary goal of this research was to prove and demonstrate the ability of IS to evaluate soil properties and quality across anthropogenically induced land-use changes. This aim was fulfilled by developing and implementing a spectral soil quality index (SSQI) using IS obtained by a laboratory and field spectrometer (point scale) as well as by airborne hyperspectral imaging (local scale), in two experimental sites located in Israel and Germany. In this regard, 13 soil physical, biological, and chemical properties and their derived soil quality index (SQI) were measured. Several mathematical/statistical procedures, consisting of a series of operations, including a principal component analysis (PCA), a partial least squares-regression (PLS-R), and a partial least squares-discriminate analysis (PLS-DA), were used. Correlations between the laboratory spectral values and the calculated SQI coefficient of determination (R2) and ratio of performance to deviation (RPD) were R2 = 0.84; RPD = 2.43 and R2 = 0.78; RPD = 2.10 in the Israeli and the German study sites, respectively. The PLS-DA model that was used to develop the SSQI showed high classification accuracy in both sites (from laboratory, field, and imaging spectroscopy). The correlations between the SSQI and the SQI were R2 = 0.71 and R2 = 0.7, in the Israeli and the German study sites, respectively. It is concluded that soil quality can be effectively monitored using the spectral-spatial information provided by the IS technology. IS-based classification of soils can provide the basis for a spatially explicit and quantitative approach for monitoring SQ and function at a local scale.

Hyperspectral remote-sensing techniques offer an efficient procedure for mineral mapping, with a unique hyperspectral remote-sensing fingerprint in the longwave infrared spectral region enabling identification of the most abundant minerals in the continental crust—quartz and feldspars. This ability was examined by acquiring airborne data with the AisaOWL sensor over the Makhtesh Ramon area in Israel. The at-sensor radiance measured from each pixel in a longwave infrared image represents the emissivity, expressing chemical and physical properties such as surface mineralogy, and the atmospheric contribution which is expressed differently during the day and at night. Therefore, identifying similar features in day and night radiance enabled identifying the major minerals in the surface—quartz, silicates (feldspars and clay minerals), gypsum and carbonates—and mapping their spatial distribution. Mineral identification was improved by applying the radiance of an in situ surface that is featureless for minerals but distinctive for the atmospheric contribution as a gain spectrum to each pixel in the image, reducing the atmospheric contribution and emphasizing the mineralogical features. The results were in agreement with the mineralogy of selected rock samples collected from the study area as derived from laboratory X-ray diffraction analysis. The resulting mineral map of the major minerals in the surface was in agreement with the geological map of the area.

Introduced in 2011, the supervised vicarious calibration (SVC) approach is a promising approach to radiometric calibration and atmospheric correction of airborne hyperspectral (HRS) data. This paper presents a comprehensive study by which the SVC method has been systematically examined and a complete protocol for its practical execution has been established—along with possible limitations encountered during the campaign. The technique was applied to multi-sourced HRS data in order to: (1) verify the at-sensor radiometric calibration and (2) obtain radiometric and atmospheric correction coefficients. Spanning two select study sites along the southeast coast of France, data were collected simultaneously by three airborne sensors (AisaDUAL, AHS and CASI-1500i) aboard two aircrafts (CASA of National Institute for Aerospace Technology INTA ES and DORNIER 228 of NERC-ARSF Centre UK). The SVC ground calibration site was assembled along sand dunes near Montpellier and the thematic data were acquired from other areas in the south of France (Salon-de-Provence, Marseille, Avignon and Montpellier) on 28 October 2010 between 12:00 and 16:00 UTC. The results of this study confirm that the SVC method enables reliable inspection and, if necessary, in-situ fine radiometric recalibration of airborne hyperspectral data. Independent of sensor or platform quality, the SVC approach allows users to improve at-sensor data to obtain more accurate physical units and subsequently improved reflectance information. Flight direction was found to be important, whereas the flight altitude posed very low impact. The numerous rules and major outcomes of this experiment enable a new standard of atmospherically corrected data based on better radiometric output. Future research should examine the potential of SVC to be applied to super-and-hyperspectral data obtained from on-orbit sensors.

Soil mineralogy is an important factor affecting chemical and physical processes in the soil. Most common minerals in soils—quartz, clay minerals and carbonates—present fundamental spectral features in the longwave infrared (LWIR) region. The current study presents a procedure for determining the soil mineralogy from the surface emissivity spectrum. Ground-based hyperspectral LWIR images of 90 Israeli soil samples were acquired with the Telops Hyper-Cam sensor, and the emissivity spectrum of each sample was calculated. Mineral-related emissivity features were identified and used to create indicants and indices to determine the content of quartz, clay minerals, and carbonates in the soil in a semi-quantitative manner—from more to less abundant minerals. The resultant mineral content was in good agreement with the mineralogy derived from chemical analyses.

Manmade crude oil contamination, which has negative impacts on the environment and human health, can be found in various ecosystems all over the globe. Hyperspectral remote sensing (HRS) is an efficient tool to investigate this crude oil contamination where its electromagnetic spectrum is analyzed. This exploratory study used an innovative HRS imagery sensor to study the effect of petroleum hydrocarbon (PHC), found in crude oil, on the spectrum of soils across the longwave infrared (LWIR 8–12 μm) spectral region. This contrasts with previous studies that focused on shortwave and midwave infrared (SWIR 1–2.5 and MWIR 3–8 μm, respectively) regions. An outdoor HRS image of three different types of soils, contaminated with 11 PHC concentrations, was processed and analyzed. Since PHC is spectrally featureless in the LWIR region, the analysis focused on the spectral alteration of the dominant minerals in the soils. Good evaluation metrics of R2 > 0.83 and a root-mean-squared-error (RMSE) between 1.06 and 1.33 wt % showed that the PHC level can be predicted with relatively good accuracy, even without direct spectral features of crude oil PHC, using an airborne LWIR camera in field conditions. This study can be used as a proof of concept for future airborne remote sensing of PHC-contaminated soils.

Remote-sensing platforms are often comprised of a cluster of different spectral range detectors or sensors to benefit from the spectral identification capabilities of each range. Missing data from these platforms, caused by problematic weather conditions, such as clouds, sensor failure, low temporal coverage or a narrow field of view (FOV), is one of the problems preventing proper monitoring of the Earth. One of the possible solutions is predicting a detector or sensor’s missing data using another detector/sensor. In this paper, we propose a new method of predicting spectral emissivity in the long-wave infrared (LWIR) spectral region using the visible (VIS) spectral region. The proposed method is suitable for two main scenarios of missing data: sensor malfunctions and narrow FOV. We demonstrate the usefulness and limitations of this prediction scheme using the airborne hyperspectral scanner (AHS) sensor, which consists of both VIS and LWIR spectral regions, in a case study over the Sokolov area, Czech Republic.

Dual pushbroom hyperspectral sensors consist of two different instruments (covering different wavelengths) that are usually mounted on the same optical bench. This configuration leads to problems, such as co-registration of pixels and squint of the field of view, known as the boresight effect. Determination of image-orientation parameters is due to the combination of an inertial measurement system (IMU) and global position system (GPS). The different positions of the IMU, the GPS antenna and the imaging sensors cause the orientation and boresight effect. Any small change in the correction of the internal orientation affects the co-registration between images extracted from the two instruments. Correcting the boresight effect is a key and almost automatic task performed by all dual-system users to better analyze the full spectral information of a given pixel. Thus, the boresight effect is considered to be noise in the system and a problem that needs to be corrected prior to any (thematic) data analysis. We propose using the boresight effect, prior to its correction, as a tool to monitor and detect spectral phenomena that can provide additional information not present in the corrected images. The advantage of using this effect was investigated with the AISA-Dual sensor, composed of an EAGLE sensor for the VIS-NIR (VNIR) region (400–970 nm) and HAWK for the SWIR region (980–2,450 nm). During the course of more than six years of operating this sensor, we have found that the boresight effect provides a new capacity to analyze hyperspectral data, reported herein. Accordingly, we generated a protocol to use this effect for three applications: (1) enhancing the shadow effect; (2) generating a 3-D view; and (3) better detecting spectral/spatial anomalies based on sub-pixel edge detection. This paper provides examples of these applications and suggests possible uses from an airborne platform.

Image registration is widely used in remote-sensing applications. The existing automatic image registration techniques fall into two categories: Intensity-based and feature-based; the latter (which extracts structures from both images) being more suitable for multi-sensor fusion, detection of temporal changes and image mosaicking. Conventional image registration algorithms have proven to be inaccurate, time-consuming, and unfeasible due to image complexity which makes it cumbersome or even impossible to discern the appropriate control points. In this study, we propose a novel method for automatic image registration based on topology (AIRTop) for change detection and multi‑sensor (airborne and spaceborne) fusion. In this algorithm, we first apply image‑processing methods (SURF—Speeded-Up Robust Features) to extract the landmark structures (roads and buildings) and convert them to a features (vector) map. The following stages are applied in GIS (Geographic Information System), where topology rules, which define the permissible spatial relationships between features, are defined. The relationships between features are established by weight-based topological map-matching algorithm (tMM). The suggested algorithm presents a robust method for image registration. The main focus in this study is on scale and image rotation, when the quality of the scanning system is constant. These seem to offer a good compromise between feature complexity and robustness to commonly occurring deformations. The skew and the anisotropic scaling are assumed to be second-order effects that are covered to some degree by the overall robustness of the sensor.

The monitoring and quantification of soil carbon provide a better understanding of soil and atmosphere dynamics. Visible-near-infrared-short-wave infrared (VIS-NIR-SWIR) reflectance spectroscopy can quantitatively estimate soil carbon content more rapidly and cost-effectively compared to traditional laboratory analysis. However, effective estimation of soil carbon using reflectance spectroscopy to a great extent depends on the selection of a suitable preprocessing sequence and data-mining algorithm. Many efforts have been dedicated to the comparison of conventional chemometric techniques and their optimization for soil properties prediction. Instead, the current study focuses on the potential of the new data-mining engine PARACUDA-II®, recently developed at Tel-Aviv University (TAU), by comparing its performance in predicting soil oxidizable carbon (Cox) against common data-mining algorithms including partial least squares regression (PLSR), random forests (RF), boosted regression trees (BRT), support vector machine regression (SVMR), and memory based learning (MBL). To this end, 103 soil samples from the Pokrok dumpsite in the Czech Republic were scanned with an ASD FieldSpec III Pro FR spectroradiometer in the laboratory under a strict protocol. Spectra preprocessing for conventional data-mining techniques was conducted using Savitzky-Golay smoothing and the first derivative method. PARACUDA-II®, on the other hand, operates based on the all possibilities approach (APA) concept, a conditional Latin hypercube sampling (cLHs) algorithm and parallel programming, to evaluate all of the potential combinations of eight different spectral preprocessing techniques against the original reflectance and chemical data prior to the model development. The comparison of results was made in terms of the coefficient of determination (R2) and root-mean-square error of prediction (RMSEp). Results showed that the PARACUDA-II® engine performed better than the other selected regular schemes with R2 value of 0.80 and RMSEp of 0.12; the PLSR was less predictive compared to other techniques with R2 = 0.63 and RMSEp = 0.29. This can be attributed to its capability to assess all the available options in an automatic way, which enables the hidden models to rise up and yield the best available model.

The albedo of bare soil depends on its organic matter, iron oxide, carbonate contents, and reflectance geometry, features considered stable over time, and also depends on salinity, moisture and roughness, which change dynamically due to agricultural practices. This paper deals with the quantitative estimation of the amount of shortwave radiation that could be reflected by air-dried bare soils in clear-sky conditions within arable lands in Israel throughout the year, assuming that they were shaped by a plough, a disk harrow, or a smoothing harrow. An area of bare soils was extracted from Landsat 8 images, within the contours of arable lands. The radiation reflected from the bare soils was calculated by equations predicting variations in their half-diurnal albedo as the solar zenith angle function. Accordingly, laboratory reflectance data of Israeli soil samples were used. The results clearly showed annual variation in the amount of short-wave radiation reflected from all bare soils within arable lands. The minimum radiation occurred in the winter, between the 1st and 70th day of the year (DOY), and the maximum was identified in the summer between 200th and 250th DOY. This could reach about 3–5 PJ/day and 16–23 PJ/day, respectively.

The purpose of this study is to evaluate a realistic feasibility of using hyperspectral remote sensing (also termed imaging spectroscopy) airborne data for mapping asphaltic roads’ transportation safety. This is done by quantifying the road-tire friction, an attribute responsible for vehicle control and emergency stopping. We engaged in a real-life operational scenario, where the roads’ friction was modeled against the reflectance information extracted directly from the image. The asphalt pavement’s dynamic friction coefficient was measured by a standardized technique using a Dynatest 6875H (Dynatest Consulting Inc., Westland, MI, USA) Friction Measuring System, which uses the common test-wheel retardation method. The hyperspectral data was acquired by the SPECIM AisaFenix 1K (Specim, Spectral Imaging Ltd., Oulu, Finland) airborne system, covering the entire optical range (350–2500 nm), over a selected study site, with roads characterized by different aging conditions. The spectral radiance data was processed to provide apparent surface reflectance using ground calibration targets and the ACORN-6 atmospheric correction package. Our final dataset was comprised of 1370 clean asphalt pixels coupled with geo-rectified in situ friction measurement points. We developed a partial least squares regression model using PARACUDA-II spectral data mining engine, which uses an automated outlier detection procedure and dual validation routines—a full cross-validation and an iterative internal validation based on a Latin Hypercube sampling algorithm. Our results show prediction capabilities of R2 = 0.632 for full cross-validation and R2 = 0.702 for the best available model in internal validation, both with significant results (p < 0.0001). Using spectral assignment analysis, we located the spectral bands with the highest weight in the model and discussed their possible physical and chemical assignments. The derived model was applied back on the hyperspectral image to predict and map the friction values of every road pixel in the scene. Combining the standard method with imaging spectroscopy may provide the required expansion of the available data to furnish decision makers with a full picture of the roads’ status. This technique’s limitations originate mainly in compositional variations between different roads, and the requirement for the application of multiple calibrations between scenes. Possible improvements could be achieved by using more spectral regions (e.g., thermal) and additional remote sensing techniques (e.g., LIDAR) as well as new platforms (e.g., UAV).

The overarching goal of this paper was to espouse methods and protocols for water productivity mapping (WPM) using high spatial resolution Landsat remote sensing data. In a world where land and water for agriculture are becoming increasingly scarce, growing “more crop per drop†(increasing water productivity) becomes crucial for food security of future generations. The study used time-series Landsat ETM+ data to produce WPMs of irrigated crops, with emphasis on cotton in the Galaba study area in the Syrdarya river basin of Central Asia. The WPM methods and protocols using remote sensing data consisted of: (1) crop productivity (ton/ha) maps (CPMs) involvingcrop type classification, crop yield and biophysical modeling, and extrapolating yield models to larger areas using remotely sensed data; (2) crop water use (m3/ha) maps (WUMs) (or actual seasonal evapotranspiration or actual ET) developed through Simplified Surface Energy Balance (SSEB) model; and (3) water productivity (kg/m3) maps (WPMs) produced by dividing raster layers of CPMs by WUMs. The SSEB model calculated WUMs (actual ET) by multiplying the ET fractionby reference ET. The ETfraction was determined using Landsat thermal imagery by selecting the “hot†pixels (zero ET) and “cold†pixels (maximum ET). The grass reference ET was calculated by FAO Penman-Monteith method using meteorological data. The WPMs for the Galaba study area demonstrated a wide variations (0-0.54 kg/m3) in water productivity of cotton fields with overwhelming proportion (87%) of the area having WP less than 0.30 kg/m3, 11% of the area having WP in range of 0.30-0.36 kg/m3, and only 2% of the area with WP greater than 0.36 kg/m3. These results clearly imply that there are opportunities for significant WP increases in overwhelming proportion of the existing croplands. The areas of low WP are spatially pin-pointed and can be used as focus for WP improvements through better land and water management practices.

Soil spectroscopy has shown to be a fast, cost-effective, environmentally friendly, non-destructive, reproducible and repeatable analytical technique. Soil components, as well as types of instruments, protocols, sampling methods, sample preparation, spectral acquisition techniques and analytical algorithms have a combined influence on the final performance. Therefore, it is important to characterize these differences and to introduce an effective approach in order to minimize the technical factors that alter reflectance spectra and consequent prediction. To quantify this alteration, a joint project between Czech University of Life Sciences Prague (CULS) and Tel-Aviv University (TAU) was conducted to estimate Cox, pH-H2O, pH-KCl and selected forms of Fe and Mn. Two different soil spectral measurement protocols and two data mining techniques were used to examine seventy-eight soil samples from five agricultural areas in different parts of the Czech Republic. Spectral measurements at both laboratories were made using different ASD spectroradiometers. The CULS protocol was based on employing a contact probe (CP) spectral measurement scheme, while the TAU protocol was carried out using a CP measurement method, accompanied with the internal soil standard (ISS) procedure. Two spectral datasets, acquired from different protocols, were both analyzed using partial least square regression (PLSR) technique as well as the PARACUDA II®, a new data mining engine for optimizing PLSR models. The results showed that spectra based on the CULS setup (non-ISS) demonstrated significantly higher albedo intensity and reflectance values relative to the TAU setup with ISS. However, the majority of statistics using the TAU protocol was not noticeably better than the CULS spectra. The paper also highlighted that under both measurement protocols, the PARACUDA II® engine proved to be a powerful tool for providing better results than PLSR. Such initiative is not only a way to unlock current limitations of soil spectroscopy, but also offers considerable efficiency and cost- and time-saving possibilities, which lead to further improvements in prediction performance of spectral models.

The study tested a data mining engine (PARACUDA®) to predict various soil attributes (BC, CEC, BS, pH, Corg, Pb, Hg, As, Zn and Cu) using reflectance data acquired for both optical and thermal infrared regions. The engine was designed to utilize large data in parallel and automatic processing to build and process hundreds of diverse models in a unified manner while avoiding bias and deviations caused by the operator(s). The system is able to systematically assess the effect of diverse preprocessing techniques; additionally, it analyses other parameters, such as different spectral resolutions and spectral coverages that affect soil properties. Accordingly, the system was used to extract models across both optical and thermal infrared spectral regions, which holds significant chromophores. In total, 2880 models were evaluated where each model was generated with a different preprocessing scheme of the input spectral data. The models were assessed using statistical parameters such as coefficient of determination (R2), square error of prediction (SEP), relative percentage difference (RPD) and by physical explanation (spectral assignments). It was found that the smoothing procedure is the most beneficial preprocessing stage, especially when combined with spectral derivation (1st or 2nd derivatives). Automatically and without the need of an operator, the data mining engine enabled the best prediction models to be found from all the combinations tested. Furthermore, the data mining approach used in this study and its processing scheme proved to be efficient tools for getting a better understanding of the geochemical properties of the samples studied (e.g., mineral associations).

Hyperspectral remote-sensing in the reflected infrared and thermal infrared regions offers a unique and efficient alternative for mineral mapping, as most minerals exhibit spectral features in these regions, mainly in the shortwave and longwave infrared. Airborne hyperspectral data in both spectral regions, acquired with the AisaFENIX and AisaOWL (Specim) sensors over Makhtesh Ramon in Israel, were analyzed. Calculating the reflectance and emissivity spectra of each pixel in the shortwave infrared and longwave infrared region images, respectively, and determining mineral indices enabled identifying the dominant minerals in this area—kaolinite, calcite, dolomite, quartz, feldspars and gypsum—and mapping their spatial distribution in the surface. The benefit of using hyperspectral data from both reflected infrared and thermal infrared regions to improve mineral identification was demonstrated.

We explored the effect of raindrop energy on both water infiltration into soil and the soil's NIR-SWIR spectral reflectance (1200–2400 nm). Seven soils with different physical and morphological properties from Israel and the US were subjected to an artificial rainstorm. The spectral properties of the crust formed on the soil surface were analyzed using an artificial neural network (ANN). Results were compared to a study with the same population in which partial least-squares (PLS) regression was applied. It was concluded that both models (PLS regression and ANN) are generic as they are based on properties that correlate with the physical crust, such as clay content, water content and organic matter. Nonetheless, better results for the connection between infiltration rate and spectral properties were achieved with the non-linear ANN technique in terms of statistical values (RMSE of 17.3% for PLS regression and 10% for ANN). Furthermore, although both models were run at the selected wavelengths and their accuracy was assessed with an independent external group of samples, no pre-processing procedure was applied to the reflectance data when using ANN. As the relationship between infiltration rate and soil reflectance is not linear, ANN methods have the advantage for examining this relationship when many soils are being analyzed.

The commonly used analytic method for assessing total petroleum hydrocarbons (TPH) in soil, EPA method 418.1, is usually based on extraction with 1,1,2-trichlorotrifluoroethane (Freon 113) and FTIR spectroscopy of the extracted solvent. This method is widely used for initial site investigation, due to the relative low price per sample. It is known that the extraction efficiency varies depending on the extracting solvent and other sample properties. This study’s main goal was to evaluate reflectance spectroscopy as a tool for TPH assessment, as compared with three commercial certified laboratories using traditional methods. Large variations were found between the results of the three commercial laboratories, both internally (average deviation up to 20%), and between laboratories (average deviation up to 103%). Reflectance spectroscopy method was found be as good as the commercial laboratories in terms of accuracy and could be a viable field-screening tool that is rapid, environmental friendly, and cost effective.