Vis-NIR Spectroscopy and Satellite Landsat-8 OLI Data to Map Soil Nutrients in Arid Conditions: A Case Study of the Northwest Coast of Egypt

The mapping of soil nutrients is a key issue for numerous applications and research fields ranging from global changes to environmental degradation, from sustainable soil management to the precision agriculture concept. The characterization, modeling and mapping of soil properties at diverse spatial and temporal scales are key factors required for different environments. This paper is focused on the use and comparison of soil chemical analyses, Visible near infrared and shortwave infrared VNIR-SWIR spectroscopy, partial least-squares regression (PLSR), Ordinary Kriging (OK), and Landsat-8 operational land imager (OLI) images, to inexpensively analyze and predict the content of different soil nutrients (nitrogen (N), phosphorus (P), and potassium (K)), pH, and soil organic matter (SOM) in arid conditions. To achieve this aim, 100 surface samples of soil were gathered to a depth of 25 cm in the Wadi El-Garawla area (the northwest coast of Egypt) using chemical analyses and reflectance spectroscopy in the wavelength range from 350 to 2500 nm. PLSR was used firstly to model the relationship between the averaged values from the ASD spectroradiometer and the available N, P, and K, pH and SOM contents in soils in order to map the predicted value using Ordinary Kriging (OK) and secondly to retrieve N, P, K, pH, and SOM values from OLI images. Thirty soil samples were selected to verify the validity of the results. The randomly selected samples included the spatial diversity and characteristics of the study area. The prediction of available of N, P, K pH and SOM in soils using VNIR-SWIR spectroscopy showed high performance (where R2 was 0.89, 0.72, 0.91, 0.65, and 0.75, respectively) and quite satisfactory results from Landsat-8 OLI images (correlation R2 values 0.71, 0.68, 0.55, 0.62 and 0.7, respectively). The results showed that about 84% of the soils of Wadi El-Garawla are characterized by low-to-moderate fertility, while about 16% of the area is characterized by high soil fertility.


Introduction
Soil is a very complex ecosystem made up of biotic and abiotic factors that strongly differ from one environment to another. The characterization, modelling and mapping of soil properties are key factors for implementing good agricultural management practices [1][2][3][4] to maintain ecological balances and prevent land degradation in arid and semiarid environments. As an example, the accumulation of salts and soil nutrients in arid conditions is affected by many factors, such as topography, geology, climate, soil moisture, land use, agricultural activity, and local environmental conditions [5][6][7][8]. The traditional methods for estimating different soil properties typically involves extensive field work and laboratory analysis and, therefore, are not only expensive and time consuming but also may be affected by significant uncertainty. Therefore, over the last four decades, to model and map soil properties in a cost-effective manner at various scales, remotely sensed imagery has been proposed and used in combination with field measurements [9][10][11][12]. Important soil properties such as salinity, texture, minerals, and organic matter have been successfully characterized and investigated using multispectral scanner (MS), Landsat-8 operational land imager (OLI), Landsat-5 thematic mapper (TM), Landsat-7 enhanced thematic mapper plus (ETM+) [13].
Over the past two decades, scientists throughout the world have focused their interest on new technologies such as the visible-near-infrared (Vis-NIR) spectroscopy to identify and characterize soil in terms of (but not only) clay mineralogy, soil organic matter (SOM), soil composition, and soil texture [14][15][16][17]. It is well recognized that the absorption spectrum in the NIR zone (780-2500 nm) can be used for estimating H2O, CO2, OH, SO4, and CO3 groups [18]; furthermore, soil nutrients can be identified using NIR spectroscopy, particularly for estimating N, K, and P soil content (with expected satisfactory coefficients of correlation around 0.72 and 0.68 for N and K, respectively), and with higher value in the case of phosphorus (around 0.84) [19]. Moreover, additional independent studies have shown that calcium, potassium, magnesium, and sodium can be predicted using statistical models such as partial least-squares regression (PLSR) [20,21]. As a whole, today, Vis-NIR techniques are recognized to be effective for the quantitative retrieval of soil characteristics and usually provide good indications of soil quality [22]. Nevertheless, some critical issues have still to be faced, such as, for example, the estimation of carbonate and the gypsum contents that is still today a controversial issue. In fact, some studies highlighted that spectroscopic techniques cannot suitably predict the carbonate content (correlation lower than 0.52), whereas other studies pointed out that the joint use of spectroscopic techniques and PLSR improved the estimation with correlation values ranging from 0.86 to 0.91 [23][24][25].
From the methodological point of view, analytical methods based on changes in specific reflectance (in the visible range from 400 to 700 nm, and in the near-infrared range from 700 to 2500 nm [26,27]), enable the discrimination of different soil properties, such as pH, organic carbon, electrical conductivity, texture, nitrate-nitrogen, available phosphorus, exchangeable potassium, cation exchange capacity, exchangeable calcium, and exchangeable aluminum. Moreover, several prediction models have been used to assess soil properties based on reflectance spectroscopy, such as artificial neural networks (ANN), partial least square regression (PLSR), stepwise multiple linear regression (SIMR), multivariate adaptive regression splines (MARS), locally weighted regression (LWR), and principal components regression (PCR) [18,28].
As a whole, today, one of the major challenges to be faced is the need to develop low-cost methods for mapping soil properties over large areas and, on the other hand, it is important to consider that agricultural management needs a rapid analysis to identify the deficiency of elements in the soil and crops. To cope with this issue, Vis-NIR reflectance spectroscopy coupled with satellite data can suitably complement in situ analyses [29,30]. The timely availability of quantitative information on soil properties and their spatial distribution is extremely relevant for sustainable agricultural to achieve development, reducing the negative effects on soil and environment [31][32][33]. This is extremely important in arid and semi-arid areas which have several limiting factors for soil fertility, such as low nitrogen, phosphorus, scarcity of irrigation water, and low soil organic matter. Moreover, the mapping soil properties and fertility provides good indicators of land degradation [34][35][36] and/or evidence of land capacity.
An effort in this context is made in this paper, which is focused on the evaluation of soil nutrients (N, P, K), SOM, and pH in the arid area of Wadi El-Garawla (the northwest coast of Egypt) jointly using chemical analyses, Vis-NIR spectroscopy and satellite Landsat-8 data that are freely available from the NASA web site. In detail, the PLSR was used firstly to model the relationship between the averaged values from the analytical spectral devices (ASD) spectroradiometer and the soil's available nitrogen (N), phosphorus (P), and potassium (K) content, along with the pH and the soil organic matter (SOM); and (ii) secondly to retrieve N, P, K, pH, and SOM values from the OLI images. Thirty soil samples were selected to verify the validity of the results. The randomly selected samples included the spatial diversity and characteristics of the study area.
The approach herein proposed enabled us to (I) model the relationship between the Spectral reflectance by ASD spectroradiometer and the laboratory analysis of soil of soil properties (N, P, K), pH, and SOM; (II) map the predicted soil properties using OK; (III) map the predicted soil properties from Landsat OLI images; and (IV) map the soil fertility status.
Today the availability of open satellite data from national and international space agencies strongly facilitates the investigation of soil properties, and their timely availability enables a prompt update and spatial distribution over a large area as necessary to support soil management strategies and to update information on the input parameters of crop models.

Experimental Site
The investigated area is located on the northwestern side of the coastal zone in the western desert area of Egypt. Wadi El-Garawla is located about 18 km east of city of Marsa Matruh as shown in Figure 1. The river pours into the Mediterranean Sea and extends approximately 22 km from south to north with varying slope rates [37]. The study area covers approximately 65.02 km 2 and lies between longitudes 27°14′30" and 27°24′30″E and latitudes 31°3′30″ and 31°16′0″N. Wadi El-Garawla has many varieties of environmental conditions typical for that region [38,39]. The rainfall in the studied area ranges between 105.0 to 200 mm/y and the average temperature ranges between 8.1 and 18 °C in the winter and 20 and 29.2 °C in the summer.
The study area is characterized by the scarcity of vegetation cover during the summer and autumn seasons. The vegetation begins to increase at the end of winter and spring, when seasonal herbs and plantings grow depending on the winter precipitation [1,30]. The soil temperature regime of the area is thermic and the soil moisture regime is torrid. In addition, the soils were classified in two orders-Entisols and Aridisols-and divided into five subgroups: Typic Calcigypsids, Typic Haplogypsids, Typic Haplocalcids, Typic TorriPsamments, and Lithic Torriorthents [39].

Soil Sampling and Chemical Analysis
The soil sample sites were determined based on the characteristics and the heterogeneity of the area because surface properties differ from south to north and were acquired on 15th December 2019. The amount of transported sediments is much deeper in the south. One hundred surface soil samples (0-25 cm) were gathered using a random sampling method. All geomorphic units were represented by several soil samples. The collected samples were dried in the laboratory at a normal temperature and then sifted by a 2 mm sieve. The collected soil samples were chemically analyzed in a laboratory where SOM was analyzed based on Walkley and Black and soil acidity (pH) in soil saturated paste by PH meter according to previous methods [40]. The soil's available N content was measured for each soil sample using conventional chemical analysis via the Kjeldahl method. The available phosphorus content and available potassium content were determined using flame photometry [41]. Table 1 shows the basic statistics of chemical analysis and shows that the soil of the study area is slightly to moderately alkaline with pH values from 6.56 to 8.97. Total soluble salts differed widely from one site to another and had a wide range, as the electrical conductivity of the soil-saturated water (ECe) ranged between 0.11 and 10.53 dS/m. The cation exchange capacity (CEC) also differed from one site to another due to the ratio of the fine fraction and soil organic matter percentage, which ranged between 0.86 and 5.66 cmol/kg. The calcium carbonate percentage of the soils had a wide range, between 2% and 37%. The soil organic matter percentage (SOM%) ranged from almost none (0.04%) to low (1.57%).

Digital Image Processing
Operational land imager (OLI) Landsat 8 images are characterized by 15 m panchromatic and 30 m multi-spectral spatial resolutions with nine spectral bands. Firstly, two OLI images acquired on 15 December 2019 were downloaded from the U.S. Geological Survey (USGS). In particular, the blue to short-wave infrared portion of the spectrum were used in this study. The thermal bands were excluded, and the images were geo-rectified according to UTM coordinates. All further digital image processing and analyses of the OLI satellite images were executed using the standard approaches provided by the ENVI software. Afterwards, all OLI images were atmospherically corrected using the FLAASH module, and the spatial resolution of the visible/NIR bands was resampled to 15 m depending on the panchromatic band. The data were represented by calibration to spectral radiance and then transformed to surface reflectance [41]. The images were mosaicked by combining multiple images into a single composite image within a dereferenced output mosaic. Finally, all satellite images were corrected and matched with the ground measurements of the study.

Spectral Measurements of the Soil Samples
Analytical spectral devices (ASDs; ASD-4 field spectroradiometer, Boulder, CO, USA) can record a complete range of 350-2500 nm spectrum of 0.1 s. Therefore, an ASD was used to collect spectra over the visible and near-infrared regions for each soil sample at 1.4-2 nm intervals with a spectral resolution of 3-10 nm. The readings were calibrated using the white reference panel. To avoid any change in radiation conditions, the white reference was checked. An ASD spectroradiometer measures the reflectance, transmission, radiance, and irradiance of an object. The recorded data are usually affected by surrounding factors, such as sources of illumination, scanning time, atmospheric conditions, and the field-of-view of the device. Therefore, a contact probe was used to control for those factors in the laboratory. Spectral data were recorded concerning an external white reference panel. Afterwards, five spectra for each sample were recorded, and the average values for the five spectral readings were calculated. Thus, one value was obtained to express the spectral characteristics of each sample [9,42,43].

Model Calibration and Validation
The spectral modeling of the soil data was achieved using PLSR, which is considered one of the most common approaches in Vis-NIR chemometrics analysis. This method depends on making the relation between the data matrix X and Y through a linear multivariate model [44]. PLSR algorithm integrates the compression and regression steps and selects successive orthogonal factors that maximize the covariance between predictor and response variables [44]. The advantage of PLS regression is that all available wavebands can be incorporated in the model, while earlier studies indicate that PLS models include redundant wavelengths and selecting specific wavebands can refine PLS analyses [45]. The soil samples were representative of the variation soil types in the Wadi El-Garawla basin. Using a leave-one-out cross validation, the dominant absorption features of each soil variable (N, P, K, pH and OM) were determined using PLSR. One hundred soil samples were randomly divided into a subset of 70 samples used for calibration of a subset of 30 samples for validation. Modeling was performed using the PLSR adopted because it usually provides promising results for Vis-NIR analysis [46,47]. The PLSR models (one for each soil parameter) were evaluated by the coefficient of determination (R2), the root means square error (RMSE), and the mean of response (MR). In addition, R2 was used to describe the model validation, where "x" represents the soil parameter values (N, P, K, pH, and SOM), which was measured using chemical laboratory analyses and used as the reference values for the calibration phase, "y" is the predicted value, and "n" is the number of soil samples used for the calibration [48,49]. MR, RMSE, and root means square standardized error (RMSSE) were calculated according to Equation (2) [50], and NRMSE was applied according to Equation (3) [51].
where is the total number of samples, and pi is the vector of predicted values of the variable being predicted, with oi being the observed values.
where n is the number of observations or samples; o is the osi is the standardized observed value at place i; psi is the standardized predicted/estimated value at place I where NRMSE is defined as the normalized root mean square error, RMSE as the root mean square error, and as the standard deviation of y, which is used in [51], where it is explained that the standard deviation (sd)-based NRMSE represents the ratio between the variation not explained by the regression vs. the overall variation in y. Thus, if the regression explains all of the variation in , nothing is unexplained, and the RMSE, and consequently the NRMSE, is zero. If the regression explains some parts and leaves other parts unexplained, which is at a scale similar to the overall variation, then the ratio will be around 1.

Mapping Soil Properties Using Ordinary Kriging
Ordinary kriging was used to interpolate the predicted soil values obtained from the PLSR model. OK is a geostatistical model that uses a set of statistical tools to predict the value of a given soil property (N, P, K, pH, and SOM) at a location that was not sampled. The normal distribution pattern of the data was checked using the histogram tool and normal QQPlots. The trend analysis was a check for each parameter.
The general equation of the kriging estimator method is as follows [51,52]: * where * is the estimated variable at the location, represents the values of an inspected variable at the location, is the statistical weight that is offered to the sample located near , and is the number of observations in the neighborhood of the inspected point. The semivariogram of the selected soil parameter was achieved using the average squared differences among all pairs of values according to Equation (4) [52].
where ℎ is the semivariance for the interval distance classh, ℎ is the number of pairs of the lag interval, is the measured sample value at point , and ℎ is the measured sample value at the position ℎ .
We applied the multiple semi-variogram models (linear plateau, circular, spherical, exponential, exponential, and Gaussian) for each parameter dataset. The validation and suitability of each model was tested via such parameters as the root mean square error (RMSE), the mean standardized error (MSE), and the root mean square standardized error (RMSSE) [50][51][52]. All data processing and analysis for OK were done in the ArcGIS software package, version 10.4.

Mapping N, P, K, SOM, and pH Using the Landsat 8 OLI
Spectral reflections obtained by the ASD spectroradiometer were used to determine the N, P, K, SOM, and pH values from the OLI images. The average of the spectral reflectance of the ASD was calculated at regions similar to that of the OLI images: blue (0.450-0.515 μm), green (0.525-0.600 μm), red (0.630-0.680 μm), NIR (0.845-0.885 μm), SWIR1 (1.560-1.660 μm), and SWIR2 (2.100-2.300 μm). PLSR was used for relationship modeling between the averaged values from the ASD spectroradiometer and the N, P, K, pH, and SOM values. Consequently, the models were applied to retrieve N, P, K, pH, and SOM values from OLI images. Thirty soil samples were selected to verify the validity of the results. The randomly selected samples included the spatial diversity of the topographic characteristics of the study area, where the element concentrations differed from one location to another, as they are usually related to the surrounding factors. Finally, the resulting maps were validated by comparing the predicted values with the laboratory values using the correlation coefficient (R2) and the root mean square error (RMSE). The stepwise linear regression model was used to conduct a regression analysis between the spectral band calculated from the ASD spectroradiometer and soil laboratory analysis in situ for N, P, K, PH, and SOM.

Soil Fertility Status of Wadi El-Garawla
Soil fertility status (SFS) represents the nutrient content of soil, the available nitrogen, phosphorus, and potassium, the organic matter content, and the soil reaction (pH) (Figure 2) which indicates the degree of suitability for most crops for specific uses. In this study, we relied on the criteria for crop growth and needs, which were suggested in [28,31]. Five factors were selected in this study for evaluating the fertility degree for most crops, as shown in Table 2. The selected factors were the available N, P, and K, the SOM, and the pH, which were produced based on the spectroscopy techniques using the above methodology. Each factor was reclassified using the Arc GIS spatial model, and its weight was taken according to standard methods. The soil fertility status was evaluated using the GIS spatial model based on the following Equation (6)  S ava. N S ava. P S ava K S OM S pH (6) where S is the score factor and ava. N, ava P, ava K, OM, and pH are factors that express, respectively. Figure 2 shows the steps of mapping different soil nutrients based on the integration of spectral reflectance and satellite images.   Table 3 shows the basic statistical data of the predicted soil analysis (available N, P, and K, as well as the pH and SOM). The maximum values of the N, P, and K, the pH, and the SOM were 60

Spectral Characteristics of Studied Soil
The results showed the variance of the spectral reflections of the soil of the study area. Two dominant absorption features were observed in the ultraviolet and near-infrared wavelength ranges (355 and 1080 nm) in response to the N concentration. The change in the curve of the electromagnetic spectrum was associated with changes in the elements and the chemical composition, as the location of the response changed with K, where the response was in the NIR portion of the spectrum at 983 nm. Furthermore, the concentration of phosphorus affects several parts of the spectrum, in the red, SWIR1, and SWIR2 regions (Figure 3).

Prediction of N, P, K, pH, and SOM
The quantitative prediction of the N, P, K, pH, and SOM maps was produced using the PLSR models with an accuracy value (R 2 ) calibration of 0.89 (N), 0.72 (P), 0.91 (K), 0.65 (pH), and 0.75 (SOM). These models were successfully validated with 30% of the soil samples. The calibration and validation were evaluated by RMSE, MR, NMRSE, and coefficient of determination, as described in Table 4. The validation of the models provided reasonable results for N, P, K, pH, and SOM R 2 validation = 0.87, 0.87, 0.9, 0.69, and 0.84.  Table 5 shows the statistics of the measured and predicted values, including the maximum and minimum values and standard deviation. A variation between the estimated and predicted values was observed, the maximum values of N, P, K, pH, and SOM, and the difference between the measured and predicted values were recognized: they ranged from 60.41 to 56.33 ppm for available N, 1.32-1.81 ppm for available P, 152-151.5 ppm for available K, 8.67-9.79 ppm for pH, and 1.05%-1.23% for SOM. There was also better convergence of the mean values between the measured and predicted values, and the means ranged from 43.1 to 40.2 ppm for available N, 0.96-1.2 ppm for available P, 75.8-76.71 ppm for K, 8.01-8.03 for pH, and 0.39-0.46% for SOM. On the other hand, the minimum values exhibited weak convergence between the measured and predicted values for all factors except pH.

Mapping of Soil Nutrients of Based on Ordinary Kriging
The mapping of soil nutrient properties was based on retrieving the selected nutrient values based on the spectral fingerprints of each characteristic. The prediction models were applied by previous statistical analysis of N, P, K, pH, and SOM. Thus, OK was used to map the soil nutrients. The performance of ordinary kriging interpolation and the efficiency of the geostatistical model for each soil parameter were checked by such parameters as the RMSE, the MSE, and the RMSSE, as illustrated in Table 6. Results showed that the spherical model was suitable for available N, pH, and SOM, and the Gaussian model was suitable for available P and available K.  The map of the organic matter showed low SOM content as it varied from 0.01% to 1.23%. The results showed that the predicted pH varied from 7.12 to 9.79 with a mean pH of 8.03. In the current study, the results show that the integration of soil properties gives an acceptable overview of the fertile soil condition distribution.

Mapping of Soil Nutrients Using Landsat-8 OLI Images
Five equations were obtained, as represented in Equations (6)- (10). The obtained accuracy values for R 2 were 0.923 (N), 0.907 (P), 0.957 (K), 0.978 (pH), and 0.952 (SOM). The five models were applied to OLI imagery to retrieve the spatial distribution values of N, P, K, pH, and SOM. The accuracy assessment was achieved by the NRMSE Each of the variable values (N, P, K, pH, and SOM) were obtained as based on the averages of ASD reflectance spectroscopy and compared to the reflections of the satellite image ranges. Figure  5a-e show the spatial distribution of N, P, K, pH, and SOM, where the available N ranges between 18 and 50 ppm, available P between 0.4 and 2.8 ppm, available K between 9 and 156 ppm, soil PH between 7.30 and 8.28, and SOM between 0.02% and 1.4%. The results of the validation models indicate acceptable outputs for all the elements studied, as the R 2 was 0.7 ± 3.5, 0.68 ± 0.06, 0.55 ± 4.3, 0.62 ± 0.07, and 0.7 ± 0.02 for N, P, K, pH, and SOM, respectively. Meanwhile, the NRMSE values were 0.39, 0.29, 0.076, 0.22, and 0.18 for N, P, K, pH, and SOM respectively, as shown in Table 7.

Soil Fertility Status of Wadi El-Garawla
Spatial distribution maps resulting from spectral measurements of soil nutrients were more significant than those produced using satellite imagery. Therefore, fertility status in the study area was evaluated using spectroscopic measurements directly as performed through the spatial modeling of all characteristics. Figure 6 shows the spatial distribution of the status fertility of the study area depending on the integration of available N, P, and K, the SOM, and the pH. Three degrees of soil fertility were recognized in the study area, namely high, moderate, and low. Their respective area of coverage was 1019.136 ha (about 16% of the total study area), 2709.02 ha (43%), and 2536.9 ha (41%).

Discussion
The PLSR was herein applied to the reflectance spectra measured in the arid conditions of Wadi El-Garawla (on the northwest coast of Egypt) to model, predict, and map the available N, P, and K, pH, and SOM from in situ analysis, Vis-NEARNEAR spectroscopy, and satellite OLI data. Outputs from our analysis indicated that PLSR is reliable and effective for predicting soil nutrients and characteristics, as already found in other diverse areas in the world [16,21,52]. This confirms that the soil's available N, P, and K, the pH, and the SOM can be predicted using Vis-NEAR spectroscopy and shortwave infrared (Vis-NIR-SWIR). The variation in the spectral responses of N, P, and K throughout the wavelength range (350-2500) are associated with the diverse behavior of the diverse elements [18,53].
In addition, the absorption at 1400 and 1900 nm is referred to as the overtone bands related to water and hydroxyls [17,44,54]. On the other hand, the responses to pH were observed in three regions: blue, green, and SWIR2 [55]. Furthermore, the responses for SOM were evident in the blue-SWIR1-SWIR2 regions. As a whole, our findings are consistent with the results of a previous study [56,57]. Spectral reflections are influenced by different soil characteristics and the concentrations of different elements [58]. Even if the determination of the parts of the wavelength that respond to dynamic changes in the concentration of elements is a complicated process, the statistical analysis can overcome these problems.
The obtained models to predict the soil content of nitrogen, phosphorus, and potassium indicate an ability to measure the elements with acceptable accuracy, and these outputs are consistent with many scientists [59][60][61]. The PLSR model based on the resampled measured spectra as a result of the calibration models can be effectively used to predict N, P, K, pH and SOM values, as is evident by the coefficient of determination: R2 values were 0.89, 0.72, and 0.91 for N, P, and K and 0.65 and 0.75 for pH and SOM, respectively. The results show the capability of reflectance Vis-NEAR spectroscopy to predict soil pH with a correlation of validation of 0.69, and these findings are consistent with [62][63][64][65]. Furthermore, the validation of the models of SOM was 0.48. Hence, the results of the prediction of the SOM content have become acceptable and are consistent with other researchers studying arid and semi-arid areas [66,67]. The high values of pH in the study area refers to an increase in the carbonate contents, and the parent material was limestone. These calcareous lands are common in the western desert in Egypt and the north of Africa [68,69]. Soil pH is a critical factor because it affects the availability of soil nutrients to plant roots, and because it affects the biological activity in different soil environments and the activity of enzymes [70]. The map of available phosphorus shows that the soil of Wadi El-Garawla has low phosphorus content, which may be due to parent lime material. In addition, the available phosphorus is associated inversely with soil pH values, which is consistent with [69][70][71].
The spatial distribution of SOM in the study area shows that the area, in general, is poor in the proportion of organic matter, less than 0.5% in most of the region. In the northern parts, a relative increase in the percentage of organic matter was observed due to the presence of seasonal crops. Despite the increase of SOM in the north of the area, it was not observed that it had a significant effect on reducing the soil pH in the study area, except for some small areas where a decrease in soil pH with increasing SOM was observed. On the other hand, the results reflect the soil characteristics in Wadi El-Garawla, which strongly vary according to several factors ranging from topographic characteristics, climate conditions, human activities, and soil types [1,72,73]. The accuracy of the spatial distribution maps of soil characteristics using OK show the acceptability degree of the results and their compatibility with other study. The RMSSE values were 1, 0.98, 1.01, 0.99, and 0.97 for N, P, K, pH, and SOM, respectively; this finding agrees with [74]. The RMSSE values were close to one, and the MSE values were close to zero for all parameters. This indicates that OK was appropriate and reliable for predicting the spatial distribution of the studied soil properties.
The study area represents the lands of the northern coast, as it is similar in climatic and topographic conditions and characterized by a low-to-moderate content of available N and P, except in the north where agricultural activities occur, and these results are consistent with those of previous studies [73,75]. The OK map of N, P, and K showed that the predicted values were associated with the SOM content. This large difference in SOM may be due to the variations in topography and climate conditions [76,77].
The integration of Landsat-8 OLI images with spectral measurements provided satisfactory results on the distribution of soil nutrients in Wadi El-Garawla, even though the area is characterized by a low concentration of soil nutrients. These results demonstrate the effectiveness of satellite images (OLI) in predicting different soil properties, and these results are consistent with other independent investigations [18].
The R 2 and root mean square error (Table 6) confirmed the expected results of retrieved soil nutrients by satellite images, where the R 2 for nitrogen and soil organic matter were around 0.71 and 0.7, respectively. Furthermore, the R 2 for phosphorous and soil pH were 0.68 and 0.62, while 0.55 was recorded for available potassium. The accuracy results of OK interpolation were better than the results obtained by satellites, but when using high resolution images satellite imagery, the results may be better than the interpolation method [18]. Wadi El-Garawla is exposed to winter monsoon rains that cause the removal of soil nutrients from the surface layers by the slope effect, where the slope increases from south to north and can cause several environmental hazards, such as drought and desertification [72][73][74][75][76][77][78][79][80].
However, these low values are due to mismanagement of agricultural activities and the location in a semi-dry climate system. The areas that have significantly high fertility as demonstrated by agricultural activity and natural cover, as reflected in the levels of N, P, and K and organic matter. The results showed that most of the area of Wadi El-Garawla is characterized by low-to-moderately fertile soil, except for some scattered areas to the north that are characterized by high fertility. The south of the area is characterized by shallow-to-very-shallow depths and a coarse texture in addition to an undulated surface topography. The fertile soils are characterized by a deep-to-moderate soil profile, with a flat or almost flat and gently undulating surface to the north [39,[78][79][80].
Agricultural activities on Egypt's northwest coast depend on the availability of water, and that depends on monsoon rains during the winter. The results showed that the Saharan areas in Egypt are poor in their fertile content compared with the soil of Nile Delta [81][82][83][84]. This means that the soil of Wadi El-Garawla requires appropriate management that is consistent with the nature of the fertile condition. The types of crops should suit the soil's characteristics and water availability, as well as the climate of the region.

Conclusions
VNIR-SWIR spectroscopy is a very helpful technique for evaluating macronutrients, SOM, and soil pH. The current study was based on the modeling of the relationship between the spectral response and the concentrations of different elements. In detail, the PLSR was herein applied to the reflectance spectra measured in the arid conditions of Wadi El-Garawla (on the northwest coast of Egypt) to model, predict, and map the available N, P, and K, pH, and SOM from in situ analysis, Vis-NEAR spectroscopy, and satellite OLI data. Thirty soil samples were selected to verify the validity of the results. The randomly selected samples included the spatial diversity and characteristics of the study area.
As a whole, results from our investigations pointed out that the red and near-infrared regions are the most sensitive portions of the spectrum to N and K concentrations, while the red, SWIR1, and SWIR2 regions are the most sensitive to phosphorus concentrations. On the other hand, the responses to pH occur in three regions: blue, green, and SWIR2. Furthermore, the responses for SOM occur in the blue, SWIR1, and SWIR2 regions. The results indicated that spectroscopy could efficiently predict the concentrations of different elements, with R 2 values of 0.89, 0.72, 0.91, 0.65, and 0.75 for N, P, K, pH, and SOM, respectively. On the other hand, the validation of the models provided reasonable results, where the R 2 and RMSE for N, P, K, pH, and SOM were 0.87 ± 0.11, 0.87 ± 0.24, 0.9 ± 0.24, 0.69 ± 0.19, and 0.84 ± 0.12, respectively. Moreover, the use of Landsat-8 OLI images can produce acceptable results on the spatial distribution of soil nutrients, where R 2 and RMSE were 0.7 ± 3.5, 0.68 ± 0.06, 0.55 ± 4.3, 0.62 ± 0.07, and 0.7 ± 0.02 for N, P, K, pH, and SOM, respectively. Soil fertility in Wadi El-Garawla was classified into three classes: high, moderate, and low, which respectively represented about 16%, 43%, and 41% of the total space of the study area.
The results show that Wadi El-Garawla is characterized by increasing concentrations of soil nutrients in the northern parts due to the removal of these nutrients by active water erosion from the south to the north as a result of the natural slope. The results illustrated the importance of using Vis-NIR for the quantitative prediction of soil nutrients as well as an alternative to chemical analysis procedures. Wadi El-Garawla is characterized by calcareous soils and has a low availability of nutrients. Therefore, the area requires special management that takes into consideration the spatial distribution of nutrients. Suitable crops that can grow in such soils need to be selected. Moreover, it is necessary to rely on organic additives to improve the soil properties, as it helps to facilitate soil nutrients in the calcareous soil.
The methodological approach herein adopted does provide a "fast," low-cost tool that can be promptly applied widely to improve food production efficiency and significantly improve soil conservation and preservation. Our effort is also a contribution to supporting the sustainable management of soil and food production, so it provides a reference for the operational use of Earth Observation (EO)-based tools for addressing ecological problems and poverty alleviation, one of the major global goals of the 2030 Agenda for Sustainable Development (SDGs). As a whole, we proposed that the use of EO for efficient monitoring of the soil conditions can be further improved in the future on a multiscale and multi-parameter scale by integrating EO-based information with socioeconomic factors economy, population, etc., thus offering an integrated system that can be operationally adapted to suitably support and improve the efficiency of local government.