Mapping the Spectral Soil Quality Index ( SSQI ) Using Airborne Imaging Spectroscopy

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 OPEN ACCESS Remote Sens. 2015, 7 15749 (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 (R) and ratio of performance to deviation (RPD) were R = 0.84; RPD = 2.43 and R = 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 R = 0.71 and R = 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.


Introduction
Soil quality (SQ) is considered to be a key element of sustainable soil function.SQ can be defined as the capacity of a soil-specific property to sustain the productivity of plants and animals, within the limits of the ecosystem (natural or managed), while maintaining or improving the quality of water and air, and supporting human health and habitation [1,2].Soil formation or genesis is strongly dependent on the environmental conditions of both the atmosphere and the lithosphere and on long-term management [3].SQ involves physical, biological, and chemical attributes that are merged together to indicate soil functioning [4,5].Several studies have proposed conceptual frameworks for monitoring SQ [4,6].These frameworks usually share common steps, including the selection of physical, biological, and chemical soil properties that are essential in terms of soil functioning [7].Monitoring SQI at large spatial scales remains expensive, as well as time and labor consuming, when using the standard procedures due to the large number of soil analyses involved.In addition, SQ mapping requires the integrated consideration of key soil properties and their variations in space and time [8].The developments of SQ maps with complete area coverage are needed for addressing climate change, land degradation, and hydrological processes in global and regional models.However, currently only a few spatially exhaustive datasets are available [9].
Remote sensing may offer possibilities for monitoring SQ at spatial scales relevant for modeling or management.Spectroscopy techniques that include the visible (VIS, 400-700 nm), near-infrared (NIR, 700-1100 nm), shortwave infrared (SWIR, 1100-2500 nm), and thermal infra-red (TIR, 8000-12000 nm) spectral ranges are well-known tools for monitoring and studying different aspects of soil surface properties and SQ with various levels of prediction accuracy [10][11][12][13][14]. Spectroscopy is a rapid, non-destructive, reproducible, and cost-effective analytical method that is used in food science, medical science, and all geoscientific disciplines, in general, and in soil science, in particular [11,14,15].The spectral analysis procedure involves the preparation of soil samples, spectral acquisition, the preprocessing of spectral data, and the selection of an appropriate statistical model and each step can affect the accuracy of the model for each individual indicator [16].Different pre-processing transformations (PPTs) have been applied to improve prediction ability [17,18], such as mean normalization, baseline offset, maximum normalization [17], first derivatives [19], the Savitzky-Golay smoothing algorithm [20], second derivatives [21], and generalized least squares weighting (GLSW) [22].Such methods usually remove baseline effects and spectral noise, and reduce the impact of particle size [23].Choosing the best performing PPTs depends on the dataset and requires testing to ascertain an effective model [18].
Reflectance spectroscopy is used in chemometrics to construct spectral groups' classification and regression models to predict many soil attributes.While regression methods are used to model the spectral signature of a target based on specific physical, biological, or chemical soil properties, classification is used to group the spectral signatures of soil into categories [24][25][26][27].The statistical models have included parametrical methods such as partial least squares-regression (PLS-R), which is perhaps the most commonly used regression method technique (e.g., [12,14,28]).An example of a parametric classification method is the partial least squares-discriminant analysis (PLS-DA), which is a method for the supervised classification of spectral data.In the current study, we developed a conceptual framework for SQ maps based on the parametric methods (PLS-R and DA).
Recent studies have presented the ability of reflectance spectroscopy to monitor SQ.Cohen et al. [29] successfully predicted several soil properties using VIS-NIR spectroscopy as indicators of wetland SQ.They concluded that spectroscopy offers both cost and statistical power advantages for the integrative assessment of SQ in wetlands.Velasques et al. [27] successfully classified different tropical agro-ecosystem soils by using spectroscopy.Idowu et al. [30] evaluated SQ by 39 physical, biological, and chemical soil properties, using a combination of field, laboratory, and spectroscopic methods in New York State, USA, and stated that in order to advance the application of VIS-NIR spectroscopy, research should begin by establishing the spectral assessment of SQ across all soil parameters.Kinoshita et al. [23] discussed the potential of soil spectroscopy for predicting multiple soil properties simultaneously by multidimensional SQI in a western Kenyan chronosequence.They found that laboratory analyses for several physical, biological, and chemical soil properties could be replaced with a spectral analysis for SQ scores.The use of reflectance spectroscopy for the multidimensional assessment of the impact of land-use change on SQ for three land uses (afforestation, traditional grazing, and agro-pastoral) was presented by Paz-Kagan et al. [31].Their study evaluated the spectral predictability of 14 SQ indicators along with an integrated spectral soil quality index (SSQI).Askari et al [10] used spectroscopic and chemometric analyses to classify soil structural quality using 40 sites in Ireland under arable and grassland management systems.They found that spectroscopy is a suitable technique for the quantitative assessment of SQ.The above recent studies pave the way for developing an applied tool that can be used for land-use management in the context of soil security, soil health, and soil protection.
Although remote sensing and soil spectroscopy have been recognized as potentially effective techniques for monitoring SQ, they have been only rarely used with imaging spectroscopy (IS) since the first paper of Ben-Dor et al. [32].The ability to map soil with an airborne IS sensor has progressed with the development of relatively cheap and better signal-to-noise sensors that are easy to operate and available to all.IS presents a new spatial dimension for spectral information and provides a new capability for many existing soil (point) spectroscopy applications.However, moving from point measurements to the IS domain not only means a transition from micro-to macroscales, but also involves a whole new set of challenges [3,33].Among these are the low signal-to-noise ratio of most of the IS sensors (relative to the point spectrometer), large datasets, the effects of the bidirectional reflectance distribution function (BRDF), the spectral instability effect [34], the spectral mixing problem, optical shifts from one pixel to another, atmospheric interference, and crust and vegetation cover [3].Because most of the applications for SQ have been developed for point spectrometry, their immediate adaptation for the IS domain requires proper attention and adequate solutions to minimize the abovementioned problems [3].Although the potential of IS as a technique for the rapid and simultaneous prediction of soil properties is known, the challenge is to adapt the application of IS into a diagnostic screening tool that can assist in developing reliable SQ surface maps for land management at large spatial scales [35].In addition, there is a need to connect point soil spectroscopy and IS in order to relate them to functional soil properties and quality.Several studies have presented the ability to upscale laboratory models to airborne IS.These studies include, for example: soil degradation (salinity, erosion, and deposition) maps, soil mapping and classification, soil genesis and formation, soil contamination, soil water content, and more [3,33,[36][37][38][39].However, most of these studies used either single or several soil properties and not an integrated approach such as the SQI.
The primary aim of this research is to prove and demonstrate the ability of IS means to evaluate soil properties and quality across anthropogenically induced land-use changes.This objective was fulfilled by implementing the SSQI with IS data obtained by a laboratory and field point spectrometer, as well as by an airborne imaging spectrometer, in two experimental sites located in Israel and Germany.The proposed SSQI was compared with a commonly used soil quality index (SQI) determined by traditional physical, biological, and chemical soil analyses and accomplished by a series of mathematical/statistical procedures.

Study Sites
Two study areas, one located in Israel (a long-term ecological research (LTER) site) and one in Germany (an intensive research site of the TERENO Observatory "Harz/Central German Lowland, Zacharias et al. [40]), were selected for examining the abovementioned aims.The sites differ in their environmental conditions, soil properties, and management regimes and, therefore, well represent the diverse areas in which the suggested SSQI can be tested.In each research site, three management regimes were tested.

Agricultural Research Site, Northwestern Negev Desert, Israel
The Migda experimental farm belongs to a series of LTER sites in the northern Negev Desert, Israel (Figure 1A) (34°25′E, 31°22′N) [41].The total area of the Migda experimental farm is about 162 ha, and the elevation is 100-120 m a.m.s.l.The area is characterized by a mean annual rainfall of 230 mm that is concentrated during the rainy season between November and April.Average daily temperature ranges from 10 °C in the winter to 30 °C in the summer.The bedrock lithology is an Eocenian chalk.The soils are sandy-loamy loess, defined by the local Israeli definition system, and calcic haploxeralf according to the USDA's definition, and the texture composition is 20% clay, 38% silt, and 42% sand.
The farm was established in 1960 by the Agricultural Research Organization (ARO), supporting extensive agriculture that includes grazing under different grazing regimes.The main crop is spring wheat [42], growing annually during the rainy season from November to April.The farm has been grazed every year by a flock of about 800 Awassi sheep and 600 goats, starting in late February, when the field is fully covered, continuing until May (green pasture), and again from June to December (dry pasture).The grazing in the farm is managed by controlling the intensity, stocking density, flock size, and timing of herd introduction into the field.Three experimental systems were examined: (1) an abandoned agricultural field with natural vegetation in an area of 5 ha of mainly annual plants.In this field, no cultivation, irrigation, fertilization, or grazing has been performed for the last 40 years; (2) an abandoned field in an area of 9.6 ha, mainly with annual plants, with grazing, but with no cultivation, fertilization, or irrigation; and (3) a monocultural agro-pastoral field of wheat in an area of 9.5 ha, with moderate grazing and cultivation but with no fertilization or irrigation.

Harz/Central German Lowland Agricultural Site
The Schä fertal research catchment is located in the lower eastern Harz Mountains, Germany (Figure 1B) (11°02′E, 51°39′N).The study area is about 144 ha, and the elevation is 393 m a.m.s.l.The climate regime is moderate continental that is influenced by its leeward position.The average annual rainfall is 636 mm with a mean annual temperature of 6.8 °C, ranging from 1.8 °C in January to 15.5 °C in July.Luvisols and Cambisols soil groups, according to the FAO's classification [43], cover the slopes and are intensively used for agriculture [44].Lower relief positions near the channel are characterized by the Gleysols soil group that is mainly used for pasture.All of the soils have a loamy character with a texture composition of about 21% clay, 53% silt, and 24% sand [45].The entire catchment is divided into five fields, managed by one agricultural cooperative society that provides information such as fertilization and crop yield.Ploughing and harrowing is conducted in autumn.Major crops in the catchment are winter cereal and winter rapeseed, whereas 3% of the area is forested.Three experimental systems were examined: (1) a fertilized agricultural field that was fertilized several days before the flight campaign with urea and ammonium nitrate.In addition, the field was harvested and sown with winter rapeseed; (2) an unfertilized agricultural field that was harvested and sown with winter rapeseed.The latest fertilizer application (nitrogen and sulphur) on the field was in April 2013; and (3) a forest area.

Data Collection and Analysis
The current study involves several spectral datasets as illustrated in Figure 2 (laboratory, field, and airborne).Two flight campaigns were performed using Aisa (AisaDUAL and AisaFENIX) hyperspectral sensors onboard light aircrafts.One flight campaign took place on 6 September 2013 over the German site and the other on 19 September 2014 over the Israeli site.Both campaigns were carried out at the end of the dry season to minimize the effects of soil water content, annual vegetation, and cloud effects.During the campaigns, the skies were clear and the agricultural fields were exposed with bare soil.The flight campaigns were accompanied by field and laboratory spectral measurements, and by soil sampling in which the samples were transferred to laboratories for physical, biological, and chemical analyses.The soil sampling design in the two sites covered the different experimental systems, and the sampling points were randomly selected in each system.The spectral and soil data analyses were jointly subjected to a chain of mathematical/statistical procedures, namely PCA, PLS-R, and PLS-DA.

Field and Laboratory Spectral Measurements
In each of the study sites, the soil spectra were acquired under field conditions (undisturbed samples) and under laboratory conditions with the portable analytical spectral device (ASD) Field Spec® Pro spectrometer.Prior to the laboratory spectral measurements, the soil samples were air dried, ground, and passed through a 2-mm sieve.The ASD covers a spectral range of 350-2500 nm and has a 25° field of view.The instrument was repeatedly calibrated to spectral reflectance using a standard white reference panel (Spectralon Labsphere Inc., North Sutton, NH, USA).A total of 50 soil samples were collected in the Migda site (20 in the monocultural agro-pastoral field; 20 in the abandoned field with grazing; and 10 in the abandoned agricultural field).A total of 25 soil samples were collected in the Schä fertal site (10 in the fertilized agricultural field; 10 in the unfertilized agricultural field; and 5 in the forest).Reflectance data in the laboratory were measured under stable illumination from two directions while the spectrometer's fiber aperture was fixed at a constant height of 18.5 cm above the sample platform; thus, the diameter of the sample was 8 cm.The bidirectional illumination reduced the effects of micro-topography shadowing.In addition, each sample was measured four times, while rotating the sample 90° between each reading.These four readings were later averaged to a final value representing the spectral resolution of 1 nm.The spectral resolution of the ASD varies from 3 nm in the VIS-NIR range to 8-10 nm in the SWIR range.We resampled the ASD's spectral band to 1 nm uniformly along the entire spectral region.Finally, the data were spectrally resampled to the spectral configuration of the airborne sensors (see next section).The field spectral measurement was vertical in relation to the soil surface using the same ASD spectrometer in the Israeli site that was used in the laboratory, however with unexposed fiber (i.e., exposed to interior fixed lighting and not to the sun).To reduce spectral noise, four spectral readings for each soil sample were measured and averaged to a final value representing the field sample.In the Schä fertal site, undisturbed soil samples were measured in the same laboratory conditions as mentioned above.The soil samples in Schafertal site were measured in the same laboratory conditions since the weather did not allow ASD field measurements.PPTs of the spectral data included three types: (1) A second-order transformation of the polynomial Savitzky-Golay smoothing [20,46] was performed to minimize variance between samples caused by grinding and the optical setup.This transformation was found to be an optimal spectral pretreatment in similar studies [19,47,48].(2) An autoscale transformation was performed, which is an exceptionally common pre-processing method that uses mean-centering followed by dividing each variable by the standard deviation (SD) of the variable.This approach is a valid method with which to correct differing variable scaling and units if the predominant source of variance in each variable is signal rather than noise.Under these conditions, each variable is scaled such that its useful signal has an equal footing with the other variables' signals [49].Finally, (3) a GLSW transformation was performed to produce a filter matrix based on the differences between pairs or groups of samples that should otherwise be similar [50].The single adjustable parameter, α, that defines how strongly GLSW downweights interferences was set to 0.02.Adjusting α towards larger values (typically above 0.02) decreases the effect of the filter, while smaller α values (typically 0.001 and below) apply more filtering.The prediction process included randomly selecting the spectral sample dataset (25%) to calibrate the model for the prediction of known samples.

Airborne Imaging Spectroscopy
Flight campaigns were performed over the two study systems using two different SPECIM's AISA hyperspectral airborne sensors, AisaFENIX (in Israel) and AisaDUAL (in Germany).AisaFENIX is a hyperspectral sensor with a single optic that provides a continuous image, covering the VIS-NIR-SWIR spectral region (380-2500 nm) with 448 spectral bands and a FWHM of 3.5-5.5 nm.The sensor was mounted onboard a light aircraft that was flown at an altitude of 762 m, resulting in a spatial resolution of 1 m, under good weather conditions, over the study area at the Migda site.Another data cube was acquired over a nearby ground calibration site, located on the flight trajectory, and included net targets usually used for a supervised vicarious calibration (SVC) procedure [51].The AisaDUAL sensor is a hyperspectral sensor with a computable pushbroom airborne imaging spectrometer, covering the 420-2450 nm spectral range with 367 spectral bands.A bandwidth of 4.5 nm exists in the VIS and NIR regions, and of 6.3 nm in the SWIR region.This sensor combines the AisaEAGLE (400-970 nm) and the AisaHAWK (970-2500 nm) sensors with different optics, boresighted to each other, resulting in a 2-3 pixel offset.The AisaDUAL sensor was mounted onboard a light aircraft that was flown over the study area at the Schä fertal site, at an altitude of 768 m, resulting in a spatial resolution of 1 m, under good weather conditions.An inertial measurement unit (IMU) enabled global positioning system (GPS) device was used to record the geo-positional data of the aircrafts for geometric rectification in the Migda and Schä fertal sites.
AisaFENIX raw data were processed using the CALIGEO software (Spectral Imaging Ltd), resulting in georectified radiance images.The Migda site radiance image was atmospherically corrected using the ACORN-6 software (ImSpec LLC, Advanced Imaging and Spectroscopy) with a gain factor derived from the SVC site image processing.The resultant reflectance image was used henceforth.The AisaDUAL raw data were radiometrically corrected with the CALIGEO software by Spectral Imaging Ltd.Afterwards, the reduction of miscalibration effects (ROME) destriping algorithm was applied to reduce sensor miscalibration effects [52].The geometric correction of the AisaDUAL data was performed with the PARGE software (ReSe Applications).The BREFCOR method was used for a surface-cover-dependent correction of bidirectional effects [53].BREFCOR was implemented in the ATCOR-4 package version 6.3 (ReSe Applications) and was also used to calculate the reflectance values for the image.The positional accuracy was improved with geo-referencing and high resolution orthophotos.At the end, six flight lines were stitched into one mosaic.
To minimize the vegetation signal in the image, the normalized difference vegetation index (NDVI; Tucker [54]) was calculated in the images.A threshold of NDVI ≤ 0.3 was selected to represent bare soil [55].The areas with an NDVI that was higher than 0.3 were masked out from further analysis.Due to vegetation cover in the abandoned agricultural field with natural vegetation in Migda, and in the forest in the Schä fertal site, the analyses of the images were further conducted only in the fields with bare soil.Following the image pre-processing, the spectra of the selected soil samples were extracted according to the GPS points that were sampled in the field.The extraction of the pixels from the images was calculated by averaging the pixels in a 1-m area around the selected point.In addition, elimination of the spectral channels affected by atmospheric water absorption in the images was done (excluded were the 1340-1448 and 1792-2000 spectral bands).
Figure 3 presents the average spectral signature of the laboratory, field, and airborne data, in the different experimental systems in the two study sites.Spectral differences exist between the laboratory, field, and airborne datasets.The spectral soil signatures are affected by soil aggregation, particle size distribution, water content, and the presence of additional elements, such as vegetation litter, rock particles, and mineral deposits.Therefore, laboratory standardization of the sample preparation and measurement protocols was carried out to minimize these BRDF effects.The spectral soil signatures from the field data did not include preliminary preparation of the soil, and these signatures are affected by the soil structures and water content.The soil signature of the image is affected by the soil structure and other effects such as atmospheric interference, the BRDF effect, and more.In addition to these effects, there are a different number of spectral bands in each sensor.

Soil Sampling and Laboratory Analysis
Soil samples were collected during the flight campaigns, at a depth of 0-0.15 m.The sampling was conducted following a stratified random survey methodology.As mentioned before, in total, 50 soil samples were collected in the Migda site and 25 in the Schä fertal site.The soil sampling was accomplished by accurate GPS and surface hardness (SH) measurements.Hydraulic conductivity (HC) measurements were performed only in Israel under dry soil conditions in the field.All soil samples were transferred to the laboratory and were stored unopened at room temperature until analysis.
The Cornell Soil Health Test (CSHT) protocols were adopted for analyzing 13 physical, biological, and chemical soil properties [5,30,56].The physical properties included soil texture (fractions of clay, silt, and sand), available water capacity (AWC), surface hardness (SH), and hydraulic conductivity measured under saturated conditions (HC).The biological properties included potential active carbon (PAC) and soil organic matter (SOM).The chemical properties included pH, electrical conductivity (EC), extractable ammonium (NH4 + ), extractable nitrate (NO3 − ), extractable phosphorus (P), and extractable potassium (K).All laboratory measurements were performed with CSHT's standards [5]; however, minor modifications were introduced due to the specific management practices, climatic regions, and available tools.These included: (1) available water capacity (AWC) that was measured by the soil moisture as water retention characteristics [57]; (2) NH4 + and NO3 − that were measured by potassium chloride extracts [58]; and finally (3) the hydraulic conductivity property that was measured by a mini-disk infiltrometer in the field [59].Although the soil texture--the composition of the three fractions of clay, silt, and sand--is a main factor for the SQ, it is barely affected by soil management and, therefore, is not treated as a quality variable in itself.Thus, soil texture is not included in the SQI.

Soil Quality Index (SQI)
Evaluation of the SQ was carried out using the general approach of the SQI, involving scoring functions for each of the abovementioned soil properties [60].The scoring functions were defined in a simple nonlinear polynomial framework.Each soil property was transformed through a scoring algorithm into a unitless score (0 to 1) representing the associated level of function in that system so that the scores could be combined to form a single value [1,61].The distribution of each of the soil variables was assessed based on a Gaussian function [62,63].The interpretation of the scoring function was integrated into an index calculated by a PCA [64,65].The index values ranged from 0 to 1; low values indicated poor soils, while high values indicated healthy soils [5].
The soil property values were recorded by the different algorithms (scoring functions) to transform them into unitless scores (Si) for each soil property (i), using the following equations [23,63,64]: where x is the normally distributed soil property value,  is the baseline value of the soil property where the score equals 0.5 (inflection point) or the population mean, and b is the slope tangent of the baseline curve or 2 −2 of the population.The general shape of each function is characterized by the nature of the soil property with respect to SQ that can be either (1) more is better (mib): an upper asymptotic sigmoid curve (negative slope) that characterizes AWC, SOM, PAC, NH4 + , NO3 − , and K; (2) less is better (lib): a lower asymptote (positive slope) that characterizes SH; or (3) optimum function (op): a Gaussian function that characterizes pH, EC, P, and HC.The specific shapes and the values are determined according to prior knowledge as presented in the literature and by experts (e.g., [63,64]).All the soil measurement scores were integrated from the previous stage into a single additive index value termed an SQI (Equation ( 2)).This value is considered to be an overall assessment of SQ, reflecting management practice effects on soil function.To evaluate the index, we used the PCA procedure for pinpointing combinations of variables that describe major trends in the data: where PWi is the PCA weighting factor.Standardized PCAs of all (untransformed) data that differed significantly between treatments in the different land uses were performed using the MATLAB package (MathWorks, Natick, Massachusetts, USA) with a PLS toolbox (Eigenvector, Wenatchee, Washington, USA).The equations were normalized to obtain a maximum SQI with a score of one.
Principal components (PCs) with eigenvalues higher than one that explained at least 5% of the variations of the data were examined [4,63].Under a particular PC, only the variables with a high factor loading were retained for SQ indications.High factor loading was defined as having an absolute value within 10% of the highest factor loading.When more than one variable was retained under a single PC, a multivariate correlation was employed to determine whether the variables could be considered redundant and, therefore, eliminated from the SQI.If the highly loaded factors were not correlated, then each was considered important and, thus, retained in the SQI.Among well-correlated variables, the variable with the highest factor loading (absolute value) was chosen for the SQI.Each PC explained a certain amount of variation (percent) in the total dataset, and this percentage provided the weight for the variables chosen under a given PC.In order to evaluate the relative importance of each waveband in each of the PLS-R models, the variable importance in projection (VIP) was computed to reveal the score for each wavelength [66].

Correlation between Soil and Spectroscopy Analyses
A PLS-R cross-validation was used to correlate the spectral data with all laboratory soil measurements and the integrated SQI.A PLS-R is a predictive module technique to predict soil attributes, even when there are many predictor variables that are highly collinear [14,67].The technique is closely related to a principal components regression (PCR).However, unlike the PCR, the PLS-R algorithm selects orthogonal factors that maximize the covariance between predictor (X spectra) and response variables (Y soil laboratory data or scores).The PLS-R analysis was applied with the full cross-validation of the Venetian Blinds method [65] with a 70 to 30 split ratio between the calibration and the validation sets.The ratio between the calibration and validation datasets was 35 to 17 and 15 to 8, in the Migda site and the Schä fertal sites, respectively.The PPTs that were used include the Savitzky-Golay algorithm, auto scaling, and GLSW.The prediction was evaluated by the root mean square error of calibration and cross-validation (RMSEC and RMSECV), as well as by the coefficient of determination (R 2 ) values of the relation between the predicted and observed soil properties, and by the ratio of performance to deviation (RPD), calculated as RPD = SD/RMSECV.A particular focus was placed on RPD and R 2 because the differences in units and the SDs of the indicators made the interpretation of RMSE difficult.Using the interpretation criteria of [23], the prediction models were categorized as "excellent" with RPD ≥ 2.5 and R 2 ≥ 0.80, "good" with RPD between 2-2.5 and R 2 ≥ 0.70, "moderate" with RPD between 1.5-2 and R 2 ≥ 0.60, and "poor" with RPD ≤ 1.5 and R 2 ≤ 0.60.

Spectral Soil Quality Index (SSQI)
A PLS-DA was performed to quantify the changes in SQ in the different land uses in the two study systems.The PLS-DA is a variant of PLS modeling and aims to find the variables and directions in multivariate space that determine the known classes in a calibration set.The predictor (X spectra) and the response variables are the classification category variables (Y values) [68].The PLS-DA is a linear regression method whereby the multivariate variables, corresponding to the observations (spectral descriptors), are related to the class membership for each sample (land-use category or treatment).This method provides an understandable graphical means of identifying the spectral regions of difference between the classes and also allows for a statistical evaluation as to whether the differences between classes are significant.It takes into account the spectra and the soil attributes as a whole and classifies the spectral data.The strength of the PLS-DA model is defined by the kappa coefficient and the total accuracy calculated from the model's confusion matrix [50].
The discriminant analysis has been used in soil analysis to classify soil attributes according to the different land uses for each site.The PLS-DA model was performed on the laboratory, field, and image datasets and on a combined laboratory and image datasets.The PLS-DA classification was performed with the PLS-Toolbox (Eigenvector, Wenatchee, Washington, USA).The number of latent variables (LVs) was selected based on the recommendation of the PLS-Toolbox optimizing algorithm [50].The PLS-DA output was used to develop a scoring function in an attempt to evaluate the SQ only by spectral differences.A proportional odds logistic model was used for evaluating the scores of the spectral SQI from the output.The proportional odds model is based on the cumulative probabilities of the coefficient of variation (CV) and the (LV).Consequently, the proposed SSQI is a function of the cumulative probability scoring class, ranging from 0 to 1 as in the SQI: where T 2 refers to the Hotelling's T-squared distribution value [50] that represents a measure of the variation in each sample within the model.The T 2 indicates how far each sample is from the center of the model (score = 0) and represents the score distance as an SD within the PLS-DA; ß is the slope of the function and is calculated from a range of the minimum and the maximum values of the function.Under a particular CV, only the variables with a high factor loading were retained for the SSQI, and under a particular LV, only the variables with a high factor loading were retained.High factor loading was defined as having an absolute value within 10% of the highest factor loading [4].The cumulative variance of the model is the scores for the individual samples, and the coefficients of the LV are the weighting factors obtained from the PLS-DA model.Each LV explains a certain amount of variation in the total dataset; this percentage provides the weight for variables chosen under a given LV.

Statistical Analysis
The soil samples from each treatment (land use) were tested using a one-way analysis of variance (ANOVA) for all soil properties. .The separation of means was subjected to a Tukey Honest Significance Difference (HSD) test.Differences in soil properties were tested for significance at p ≤ 0.05 between changed land uses and by the results of the F-statistic test.The statistical analysis was performed with STATISTICA Version 11 software (StatSoft, Tulsa, OK, USA).

Soil Properties and Soil Quality
The mean value of all the soil properties from the Migda site is presented in Table 1, along with their SDs and significance values.The results of the three land uses show significant differences in all soil properties except for P. Significantly higher values of SH and SOM were found in the abandoned agricultural field with no grazing than in the other two land uses.In the agro-pastoral agricultural field, significantly higher values of HC, NH4, and K were found than in the other two land uses.Table 2 represents the respective Pearson correlation coefficients (R) for the measured soil properties.The results indicate soil properties with significant correlations (R ≥ 0.5, bold numbers) and with highly significant coefficient correlations (R > ±0.8 bold numbers with (*)).To avoid properties that could be considered as redundant, the multivariate correlation was tested for R ≥ ±0.8.Later, the soil properties with a high factor loading were eliminated from the SQI.High correlations were found between HC and SH and NH4 (R = −0.83and R = 0.79; p > 0.01, respectively) and between SH and NO3, and K (R = −0.92 and R = −0.88;p > 0.01, respectively).In addition, a high correlation was found between K and NH4 (R = 0.88; p = 0.01).
The results of the soil properties from the Schä fertal site are shown in Table 3.The three land uses show significant differences in all soil properties except for P. Significantly higher values of SH, SOM, PAC and NH4 + were found in the forest than in the agricultural land uses.In the agricultural land uses, significantly higher values of pH, EC, NO3 and K were found than in the forest.Significant differences were found between the two agricultural fields in AWC, SH, and pH.The higher pH value in the fertilized agricultural field than in the unfertilized field might be explained by the urea application several days before the campaign [69].In addition, significant differences in soil texture were observed between the three land uses.The soil in the forest had higher clay content, which is related to the ability to retain nutrients (higher cation exchange capacity) and to bind more organic matter.Table 4 represents the Pearson correlations.High negative correlations were found between SOM and pH and NO3 (R = −0.92;R = −0.79,p = 0.01, respectively).In addition, a positive correlation was found between EC and NO3 (R = 0.85; p = 0.01).

Table 1.
Soil quality properties for the Migda site, Israel, with the following treatments: (A) abandoned field no grazing; (B) agro-pastoral grazing; and (C) abandoned field with grazing.Statistics include: average value, standard deviation (SD), and indication of significant differences between treatments represented with small letters (a, b, c).

Soil Properties
Abandoned    The SQI was developed from the results of the transformed scoring of the soil properties from the two sites.Three PCs that explained 76.71% of the variance of the original data were chosen for an analysis of data redundancy from the Migda site.The PCA indicated three components with eigenvalues ≥1 (Table 5) that were selected, and the cut-off was placed at the third component.The loading values of the first PC1 (38.35% of variance) indicated HC, SH, pH, EC, and NH4, and had values within 10% of the highest value.The SH and HC were significantly correlated; therefore, they were associated with the highest loading value, and HC was selected for PC1 (38.35% of variation).For PC2 (28.13% of variation), the loading values were AWC and SOM.They were not significantly correlated, and both were selected for the SQI.For PC3 (10.23% of variation), PAC and NO3 were within 10% of the highest loading value.The total SQI values for the Migda agricultural site were, for the abandoned agricultural field with no grazing, a score of SQI = 0.48, for the agro-pastoral field, SQI = 0.61, and for the abandoned agricultural field with grazing, SQI = 0.35, with significant differences between all treatments (F = 137.5;p > 0.01).For the Schäfertal site, the PCA indicated four components with eigenvalues ≥1 (Table 6); thus the cut-off was placed at the fourth component.Four PCs that explained 81.95% of the variance of the original data were chosen for an analysis of data redundancy.The loading values of the first PC (38.9% of variance) indicated AWC, SH, SOM, PAC, and NH4, and had values within 10% of the highest value.They were not significantly correlated, and all were selected for the SQI.For the PC2 (20.38% of variation), the loading values were EC and NO3 (R = 0.85), which were correlated and associated with the highest loading value; therefore, the EC was selected for the SQI.For the PC3 (13.01% of variation), K and for the PC4 (9.66% of variation), pH were within 10% of the highest loading value.The total SQI values for the Schä fertal site were, for the fertilized agricultural field, a score of SQI = 0.51, for the unfertilized agricultural field, SQI = 0.39, and for the forest, SQI = 0.49, with significantly higher SQI in the fertilized agricultural field and the forest than in the unfertilized agricultural field (F = 13.1;p > 0.01).Figures 4 and 5 present the results of the SQI scores in the Migda site and in the Schä fertal site, respectively: (A) the SQI scores and (B) the scores of the physical, biological, and chemical components.In Figure 4A, the Migda site, significant SQI differences were found between all agricultural fields (F = 137.5;p > 0.01).In Figure 5A, the Schä fertal site, no significant SQI differences were found between the fertilized agricultural field and the forest; however, significant differences in the SQI were found in the fertilized agricultural field and the forest compared to the unfertilized agricultural field.Figures 4B and 5B show the results of the SQI representing significant differences between the physical, biological, and chemical soil components.

Spectral Correlation of Soil Quality Properties
Table 7 presents the results of the PLS-R analysis including the number of latent variables, the coefficient of determination (R 2 ), the RPD, and the VIP projection.The resultant spectral data that explain a good prediction for two sites are marked in bold in Table 7 and in Figure 6 and 7.The excellent and good results ("excellent" with RPD ≥ 2.5 and R 2 ≥ 0.80; "good" with RPD between 2-2.5 and R 2 ≥ 0.70) of the PLS-R model prediction in the Migda site include clay, sand and silt content, SH, pH, NH4, and NO3.In the Schä fertal site, clay content, SH, SOM, pH, EC, and K showed good results.In addition, the VIP was computed to reveal the score for each wavelength with excellent and good results.The sensitivity bands that were identified by the VIP for the different soil properties and SQI values are presented in Table 7.These bands are in agreement with those that were previously found in other studies [14,31,70].Figures 6 and 7 shows scatterplots of correlations between soil spectroscopy and laboratory-measured soil values, for the Migda site and Schä fertal site, respectively.The soil properties for the calibration dataset with a coefficient of determination range between 0.67 and 0.95.Figures 6 and 7 represents the results by RMSC, RMSCV, R 2 and the number of LVs for each soil property as an example.The results of the calculated SQI versus spectral measurements in the two study sites, along with their RPDs, are presented in Figure 8.The coefficients of determination in the Migda site were R 2 = 0.84 and RPD = 2.43, and in the Schä fertal site, they were R 2 = 0.78 and RPD = 2.10, both with good prediction accuracy (Figure 8).Table 7. Results of the partial least squares-regression (PLS-R) analysis in terms of spectral regions that are indicative for the soil properties in the two study systems: (1) Migda site, Israel, and (2) Schä fertal site, Germany.Statistics include the number of latent variables (LVs) selected in the PLS-R model, correlation of determination (R 2 ), and the ratio of performance to deviation (RPD).Bold numbers refer to the prediction models that were categorized as "excellent" with RPD ≥ 2.5 and R 2 ≥ 0.80, and "good" prediction models with RPD between 2-2.5 and R 2 ≥ 0.70.VIP is the variable importance in projection that presents that wavelength (in nanometers) of selected soil properties with excellent and good prediction results.

Spectral Soil Quality Index (SSQI)
Table 8 shows the total accuracy and the kappa coefficient values of the field and laboratory models, and the prediction model of the whole image in the two sites.The proportional odds in the PLS-DA classification of the spectral samples are presented in Figures 9 and 10 for the Migda and the Schä fertal sites, respectively.The results of the classification of the laboratory spectral data had a total accuracy of 1 and a kappa coefficient of 1 in the two sites.The same results were achieved after resampling the data to the AISA sensors (moving from 2000 to either 448 or 366 spectral bands).The results of the classification of the field spectral data had a total accuracy of 0.96 and 0.88, and a kappa coefficient of 0.93 and 0.88 in the Migda and Schä fertal sites, respectively.The results of the classification of the pixels that were extracted from the image showed the same results as the field spectral data, in the two sites.The results of the classification of the combined model (image and laboratory) had a total accuracy of 0.96 and 0.88, and a kappa coefficient of 0.94 and 0.82, in the Migda and Schä fertal sites, respectively.The results of the combined model (image and laboratory) improved the total accuracy.The PLS-DA provides an explicitly quantitative approach to predict the cumulative probability of soil spectral samples that belong to different soil conditions.The results show the ability to classify land uses according to soil function and the effect on the classification accuracy.As expected, moving from laboratory spectroscopy to field spectroscopy and further to IS reduces the classification accuracy.The classifications of the soil condition scores in the two systems were well separated, and each sample represents SSQI scores.The SSQI algorithm was implemented on the output results of the PLS-DA combined model (laboratory spectral data and the extracted pixels from the image).The predictions of the SSQI model in the two sites are shown in Figures 11 and 12.The classification predictions of the image had a total accuracy of 0.92 and 0.82, and a kappa coefficient of 0.91 and 0.80, in the Migda and Schä fertal sites, respectively.The SSQI indicates a higher value than the SQI; the model of the SSQI is a proportional model that is not based on the individual probability of each class but on the cumulative probabilities.Therefore, the proportions between classes that explain the changes caused by management are more essential than the actual values.The extracted pixels from the SSQIs predicted from the spectral data were significantly correlated to the SQIs calculated from the laboratory-measured data in the two sites (F = 9.75, p ≥ 0.01; F = 13.57,p ≥ 0.01, respectively).The correlations between the SSQI and the SQI were R 2 = 0.71 and R 2 = 0.7, in the Migda and Schä fertal sites, respectively.The spatial selection of the sampling points and the number of sampling affect the strength of the correlation coefficient between the SSQI and SQI.

Discussion
To the best of our knowledge, at present, most studies dealing with optical remote sensing technologies for SQ assessment focus on several soil properties as indicators of SQ for a specific soil type or land use [10], and much less research exists on the potential of IS for developing SQ surface maps.These studies mainly use point spectroscopy, either in the laboratory or in the field, and do not implement IS on a regional scale.In the move from point spectroscopy to IS, the transition is not only from micro-to macroscales, but a whole new set of challenges is encountered.In the current study, the prediction of SQIs was explored in the context of two different experimental agricultural sites under different climatic regimes using spectroscopic techniques and IS.It was found that the SSQIs from IS were not significantly different from traditional SQIs.The correlations between the extracted pixels from SSQI map and the SQI that were measured in the field were R 2 = 0.71 and R 2 = 0.7, in the Migda and Schä fertal sites, respectively.The analysis of the current study demonstrated the ability and reliability of SSQI and key soil properties to be efficiently predicted by VIS-NIR-SWIR spectroscopy.

Soil Properties and the Soil Quality Index
Many soil properties can be measured, depending on the expected soil function, land use, and management system, for assessing SQ.However, the process of measuring and analyzing a large number of soil properties is usually time-consuming and costly.In this study, the Cornell Soil Health Test (CSHT) protocols were adopted for analyzing 13 physical, biological, and chemical soil properties.Selecting indicators (soil properties) that are sensitive enough to detect the effects of management systems on SQ is essential.For example, in our study, soil texture was excluded from the analysis because it is a fixed soil property and cannot be readily influenced by management practices [4,30,71].We focused here on the effect of management on the SQ and soil properties.However, developing soil texture maps is of great importance as this property determines the water retention of the soil and can serve as an input for hydrological models and simulations [72].In addition, conventional methods for selecting proper indicators based on the PCA were used.Data redundancy was effectively reduced by using PCA and by reducing the soil properties with high multivariate correlations that were considered as redundant.Different weights were selected for each soil property according to the PCA results.These results showed that an appropriate number of soil properties can provide information for evaluating land-use impact on SQ in different study sites.
The SQI in the Migda site shows significant differences between all land uses and between their biological and chemical components.The SQI in the Schä fertal site shows no significant differences between the fertilized agricultural field and the forest land use.However, a significant reduction in the SQI value was found in the unfertilized agricultural field and the other two land uses.The differences between the two agricultural fields might be related to the urea fertilizer application on the fertilized agricultural field several days prior to the field campaign.It is known that urea fertilization leads to substantial increases in the pH value in the first three weeks after application [73].The fertilized agricultural field showed a significantly higher pH value than the unfertilized agricultural field.In addition, even if no significant differences in the SQI were found between the fertilized agricultural field and the forest in the Schä fertal site, the physical, biological, and chemical components differed significantly.This fact reflects an important advantage of the SQI model as an integrative approach that recognizes many processes in soils [74].The SQI equation is valid for establishing the degree of soil degradation as a function of soil properties [2,75,76].In addition, it enables land-use effects on soil properties to be identified and monitored by the proportions between the physical, biological, and chemical components.

Predictability of Indicators and Indices
Overall, the predictions of individual indicators and the SQI were "good" to "excellent" for several soil properties (Table 7).However, in each site, the correlation of the soil properties by spectroscopy was different and could be related to the variety of chromophores, which vary with the environmental conditions and the status of the five soil formation factors (climate, topography, parent material, organic matter, and time), and long-term management.In addition, the accuracy is related to the soil properties' concentrations and proportions, and to the number of sampling points.For example, the correlation of NO3 and NH4 were ranked as "good" to "excellent" in the Migda site, and "poor" to "medium" in the Schä fertal site (Table 7).The concentration of NO3 in the Schä fertal site was too low (>0.03 mg/100 gr) to be identified by conventional laboratory techniques, and therefore, it couldn't be identified by spectroscopic means.Spectroscopy is a reliable tool for assessing soil properties; however, it is site-specific and related to the soil function and structure.Previous studies have shown the ability of spectroscopy to predict several soil properties, such as texture, SOM, AWC, NH4, NO3, EC, and pH, with different prediction levels [15,25,31,66].These differences are probably associated with the energy of the absorbance and reflectance of molecular bonds in the NIR region due to the combination of C-H, N-H, C-O, C-N and O-H groups (chromophores) [77][78][79], minerals, water, and nutrients.
The spectral reflectance of a specific soil is a consequence of the entire chromophore interaction with the electromagnetic energy.This is related to the spectral curve that can serve as a footprint of the chromophore's overall existence in the examined matter.The ability to predict SQI is not only due to these primary soil properties, but also to a range of secondary properties that improve the prediction ability.The mechanism by which the prediction of these soil properties is possible depends on how well they are correlated with primary soil variables.The predictions of the integrated soil properties and the SQI predictions had "good" prediction accuracy (Figure 8).These results demonstrate that VIS-NIR-SWIR spectroscopy is suitable for the multidimensional evaluation of SQI.The SQIs were significantly correlated, in the Migda site (R 2 = 0.84; p < 0.01) and the Schä fertal site (R 2 = 0.78; p < 0.01).Furthermore, after finalizing the SQI, the PLS-R technique can be used rapidly, making it possible to analyze a large number of samples in a practical and timely manner as a prediction model.These qualities make spectroscopic analysis, combined with PLS-R, an attractive method for environmental monitoring, especially for modeling SQ in changed land uses [10].

Spectral Indices for Soil Quality Assessment
In both sites, no differences were found between the classification accuracy of the laboratory and the resampled laboratory model (reduction from 2000 to either 448 or 366 bands).The reduction in the spectral bands and the resampling of data to the AISA spectral bands did not affect the classification accuracy [80].Other studies showed the same trend in which no advantages were found with continuous spectral data [81].The classification accuracy depends on several factors, such as the spectral variance between classes, the number of classes, the number of sampling points, noise effects, the spectral resolution, and the PPTs.PPTs can improve the classification accuracy by increasing the variability between classes while decreasing the variability within classes.In this study, we used the second-order transformation of the polynomial Savitzky-Golay smoothing, an autoscale transformation, and GLSW.These PPTs were found to be very effective in data transformation, leading to a highly accurate classification.This conclusion is in agreement with a similar soil study [18].Pre-processing may be relevant for addressing some of the problems associated with IS.In this study, no significant differences were found between the classification accuracy of the field and the extracted pixels from the image model in both sites.This fact can be related to soil aggregation, particle size distribution, water content, and the presence of additional elements in the soil, such as vegetation litter, rock particles, and mineral deposits.However, differences between the classification accuracy of the laboratory and the field model (undisturbed soil sample) were found and can be related to the standardization of the sample preparation and the measurement procedure.The integrated model of extracted pixels from the image and the laboratory improved the classification accuracy (Table 8).In addition, the classification accuracy of the Migda site was higher than that of the Schä fertal site, which can be related to the number of sampling points that reduced the predication accuracy of the model.
IS has drawbacks relative to point spectrometry, such as a low signal-to-noise ratio, atmosphere attenuation, a varying field of view for every pixel, spectral instability, a low integration time for a given pixel, a spectral mixing problem, optical shifts from one pixel to another and more.To upscale the model to IS, additional corrections must be amplified, such as the BRDF effect, and the removal of water absorption bands, masking the vegetation signal, and more.In Figure 12, which shows the SSQI surface map of the Schä fertal site, the striping of the BRDF is shown; this is due to the mosaic of the six AISA flight lines that were stitched into one.The BRDF effect reduces the total accuracy of the SSQI model in the Schä fertal site.Because most of the applications for soil were developed for point spectrometry, their immediate adaptation for the IS domain requires proper attention and adequate solutions to minimize the above problems [3,36,82].The SSQI map (Figure 11) shows spatial pattern of soil quality, in the abandoned agricultural field with grazing.This pattern may be related to the small wadi that crosses the field and may affect the grazing density and distribution.This study provides a framework for assessing SSQI by the VIS-NIR-SWIR spectral range with IS as a preliminary tool for detecting hotspots of soil degradation at a regional scale.The PLS-DA model describes the possible separation of predefined soil conditions.Furthermore, these results demonstrate the sensitivity of reflectance spectroscopy to management changes in converted land-use soil conditions.

Monitoring Soil Quality with IS
IS is a promising tool for studying soil properties in large spatial domains.Most of the studies using IS addressed specific soil properties and not integrative SQ monitoring approaches [3].IS is a rapid, nondestructive, reproducible, and cost-effective analytical method, and therefore useful for SQ assessment.Although the IS approach is a cost-effective method, it hasn't been fully adopted because it is difficult to process and has not yet been recognized by many end-users [36].As the IS product is a geo-positioning mosaic comprising many spectral points, traditional (quantitative) approaches that successfully work for point spectrometry measurements may also be suitable for the imaging domain.The benefits of this technique include a reduction of the sampling processing time and an increase in the number of samples that can be analyzed within time and budget constraints, and hence, an improvement of the detection of changes in SQ in a given area.The SSQI is a diagnostic tool for assessing SQ in changed land uses and treatments.The correlations between the SSQI and the SQI were good in both sites.However, there is a need for the additional development of modeling approaches to create links between soil point spectroscopy and IS at a large spatial scale., Examining the cause of the soil degradation requires an extensive assessment of SQ by soil measurements and the use of models such as the SQI (using the physical, biological, and chemical components).Note that the use of IS as presented in this study can be relevant only in area with bare soil.In this study, the forests in Schä fertal site and abandoned agricultural fields in Migda site were excluded from the analysis since they were covered with vegetation.
SQ policies usually address specific management goals, such as productivity, waste recycling, and environmental protection [60].Thus, methodologies for SQ assessment should be able to measure the specific soil functions and the soil ecosystem services associated with these management goals.The ability of spectral analysis for evaluating SQ in point spectroscopy and IS will allow field implementation with IS and portable spectrometer equipment as new ways for monitoring and mapping SQ in order to address soil threats in a near real-time manner [10].This is an important implication for sustainable agricultural management, in particular, and land use, in general, because greater spatial and temporal resolution sampling will be possible at a lower cost.The SSQI maps can be used in various ways.Firstly, they may help in segmenting the landscape into soil-landscape units for each soil composition according to the land-use dynamics.Secondly, IS data can be analyzed using physically and chemically based methods to derive soil properties and soil function.Moreover, IS can be used as a data source supporting digital soil mapping [9].Finally, IS methods can facilitate the mapping of inaccessible areas by reducing the need for extensive time-consuming and costly field surveys in the future [9,83].

Conclusions
The current study underscores the potential application of reflectance spectroscopy and IS as reliable diagnostic screening tools for assessing SQ at a local scale, as well as the many soil attributes that are used to compute the SQI.The results demonstrate that: (1) correlations between the predicted spectral values and the calculated SQI were R 2 = 0.84; RPD = 2.43 and R 2 = 0.78; RPD = 2.10 in the Israeli and the German study sites, respectively; (2) the key soil properties of SQ can be predicted with "good" or "moderate" accuracy, but these are site-specific and related to the environmental conditions and the status of the five soil formation factors, as well as to long-term management; (3) the SQI for both agricultural experimental sites were predicted with "good" accuracy and were as good as those calculated from laboratory-measured data; (4) the utilized linear parametric models, PCA, PLS-R, and PLS-DA, improved the ability to predict soil properties and SQI, and they enabled us to overcome the problems of multivariant and highly collinear data such as that gained by spectroscopy; (5) the classification of soils using the PLS-DA statistical model into spectrally defined entities provides a basis for spatially explicit and quantitative definitions for developing the SSQI at a regional scale; and (6) the SSQI score related to land-use management could be differentiated by reflectance data, thus demonstrating that the SSQI is related to the productive function of the soil.The VIS-NIR-SWIR spectroscopy and IS can be used as a reproducible method for monitoring the effects of land use and for identifying hotspots of soil degradation at a regional scale.IS can provide a reliable approach for assessing SSQI and can be used for a rapid, low-cost analysis of soil conditions at appropriate scales under different management, land-use, and soil types.Thus, SSQI maps may help in evaluating the effects of management practices and in developing a better understanding of soil functions in the near future.

Figure 1 .
Figure 1.Aerial photograph of the study areas with their land-use categories: (A) Migda site, Israel; and (B) Schä fertal site, Germany.

Figure 2 .
Figure 2. Study flowchart for developing the soil quality index (SQI) and the spectral soil quality index (SSQI) by applying laboratory, point, and image spectroscopy data along with complementary soil measurements and analyses.

Figure 3 .
Figure 3.The average spectral signatures of the two study sites: Migda, Israel (A, B, C) and Schä fertal, Germany (D, E, F).The spectral dataset included: (A and D) laboratory spectral signatures; (B and E) field spectral signatures; and (C and F) airborne imaging spectral signatures.The land-use classes are characterized by the different colors.

Figure 4 .
Figure 4. (A) Scores of soil quality indices (SQIs) for the three land uses in the Migda site, Israel: abandoned field with no grazing (red color), agro-pastoral grazing (green color), and abandoned field with grazing (blue color); and (B) the SQI that was calculated by physical, biological, and chemical analyses.Capital letters above the error bars represent significant differences between land uses.Small letters within the columns represent significant differences between soil components.

Figure 5 .
Figure 5. (A) Scores of soil quality indices (SQIs) for the three land uses in the Schä fertal site, Germany: Fertilize agriculture field (red color), unfertilized agriculture field (green color), and forest (blue color); and (B) the SQI that was calculated by physical, biological, and chemical analyses.Capital letters above the error bars represent significant differences between land uses.Small letters within the columns represent significant differences between soil components.

Figure 6 .
Figure 6.Scatterplots of cross-validation (CV) predicted values versus measured values for several soil properties for the calibration dataset for all land uses in the Migda site, Israel.Calibration models were developed with a partial least squares-regression (PLS-R).RMSEC: root mean square error of calibration; RMSECV: root mean square error of cross-validation; P: phosphorus p; NH4: ammonium; NH3: nitrate, SH: surface hardness.The colors of the spots represent land-use types: abandoned field with grazing (blue), abandoned field with no grazing (red), agro-pastoral grazing (green).

Figure 7 .
Figure 7. Scatterplots of cross-validation (CV) predicted values versus measured values for several soil properties for the calibration dataset for all land uses in the Schä fertal site, Germany.Calibration models were developed with a partial least squares-regression (PLS-R).RMSEC: root mean square error of calibration; RMSECV: root mean square error of cross-validation; EC: electric conductivity; K: potassium; NO3: nitrate; SOM: soil organic matter; SH: surface hardness.The colors of the spots represent land-use types: forest (blue), fertilized agricultural field (red), unfertilized agricultural field (green).

Figure 8 .
Figure 8. Scatterplot correlation of soil quality indices (SQI) and reflectance spectroscopy values from the laboratory dataset for the changed land uses: (A) the Migda site, Israel; and (B) the Schä fertal site, Germany.Calibration models were developed with a partial least squares-regression (PLS-R).RMSEC: root mean square error of calibration; RMSECV: root mean square error of cross-validation.

Figure 9 .
Figure 9. Partial least squares-discriminant analysis (PLS-DA) classification of the different land uses in the Migda site, Israel, using data of: (A) field spectroscopy; (B) laboratory spectroscopy; (C) airborne imaging spectroscopy; and (D) merged image and laboratory spectroscopy.Dashed circles indicate the 95% confidence level.

Figure 10 .Table 8 .
Figure 10.Partial least squares-discriminant analysis (PLS-DA) classification of the different land uses in the Schä fertal site, Germany, using data of: (A) field spectroscopy; (B) laboratory spectroscopy; (C) airborne imaging spectroscopy; and (D) merged image and laboratory spectroscopy.Dashed circles indicate the 95% confidence level.

Figure 11 .
Figure 11.Hyperspectral imaging spectroscopy of the spectral soil quality index (SSQI) at Migda site, Israel.

Field No Grazing Agro-Pastoral Grazing Abandoned Field with Grazing
Note: AWC: available water content; SH: surface hardness (penetration); HC: hydraulic conductivity (infiltration), PAC: potential active carbon; SOM: soil organic matter; EC: electric conductivity; NH 4 : ammonium; NO 3 : nitrate; P: phosphorus; K: potassium; small letters indicate significant differences between treatments.Values in each column followed by the same letter do not differ significantly at p < α using an ANOVA Tukey test.a high values; b medium values; c low values.

Table 2 .
Pearson correlation coefficients for the measured soil quality properties in the Migda site, Israel.Bold numbers indicate significant differences with p ≤ 0.05 and bold numbers with (*) indicate highly significant differences with R ≥ 0.8 and p ≤ 0.01.AWC (%

Table 3 .
Soil quality properties for the Schä fertal site, Germany, with the following treatments: (A) agriculture 1; (B) agriculture 2; and (C) forest.Statistics include: average value, standard deviation, and indication of significant differences between treatments.
Note: AWC: available water content; SH: surface hardness (penetration); PAC: potential active carbon; SOM: soil organic matter; EC: electric conductivity; NH 4 : ammonium; NO 3 : nitrate; P: phosphorus; K: potassium; small superscript letters indicate significant differences between treatments.Values in each column followed by the same letter do not differ significantly at p < α using an ANOVA Tukey test.a high values; b medium values; c low values.

Table 4 .
Pearson correlation coefficients for the measured soil quality properties in the Schä fertal site, Germany.Bold numbers indicate significant differences with p ≤ 0.05 and bold numbers with (*) indicate highly significant differences with R ≥ 0.8 and p ≤ 0.01.

Table 5 .
Results of the principal component analysis (PCA) of soil properties in the Migda site, Israel.Bold and underlined values indicate underlined factors corresponding to the indicators included in the indices.

Table 6 .
Results of the principal component (PC) analysis of soil properties in the Schä fertal site, Germany.Bold and underlined values indicate underlined factors corresponding to the indicators included in the indices.