Soil samples were collected on the Mugan plain of Azerbaijan. The study are has a arid climate with a mean annual precipitation, evaporation, temperature of 30 mm , 1000 mm , 17˚C, respectively. Study area (Figure 1) is mainly located under sea level and the dominant soilsare Calcisols, Solonchaks and Calcaric Fluvisols. Salinization, gleyfication and carbonatation processes are very appreciable and plays significant role in these soils genesis and morphology [10] .

For representative result of soil properties have been selected 46 sites and 194 samples were collected over genetic horizon of soil sites. Soil samples were air dried and 2 mm sieved. Particle size distribution were determined by Pipette method, Organic matter content (OM) and N contents were determined using Walkey-Black and Kjeldahl methods, respectively; CaCO₃ content was determined with the gas-volumetric method; CEC was determined through the extraction in ammonium acetate. Soil pH and electric conductivity (EC) were determined using a 1:5 soil-water suspension

Air died and 2 mm sieved soil samples were put in a glass Petri dish for spectral scanning. The samples were scanned using ASD FielSpec 3 radiometer (Analytical Spectral Devices, Boulder , Colo. ) at the range of 350 - 2500 nm wavelength, using a contact probe with artificial light. Each soil spectrum was obtained as the mean of 10 scans. Spectralon® white reflectance standard was scanned after every 10 soils samples for white reference.

Pre-processing is important for remove undesirable affects from spectral data. A spectral pre-treatment was performed to reduce the influence of particle size and optical path length variation in the reflectance spectra [11] . The head and the tail (350 - 400 nm and 2450 - 2500 nm) of the spectra showed very high noise and were removed manually from data matrix. The remaining spectra were subjected to Multiplicative Scatter Correction (MSC) for removing scatter effect and first derivative with Savinski Golay algorithms to highlight the peaks of the spectra [12] . All pre-processing data were carried out using Unscrambler 9.7 (CAMO) (See Figure 2).

The spectral modeling process was carried out by Partial Least Squareregression (PLSR), which is one of the most common methodin Vis-NIR chemometrics analysis. PLSR is a method for relating to data matrix X and Y through a linear multivariate model [13] . PLSR is closely related to Principal component regression (PCR), but PLSR algorithm integrates the compression and regression steps and it selects successive orthogonal factors that maximize the covariance between predictor and response variables [14] . Since the number of soil samples set was not high, but very representative of the soil typologies in Mugan plains, we decided to use a leave-one-out cross validation. Then, PLSR model was calculated on all the 194 samples, whereas the validation of the model was carried out on “total set-1 sample”, iteratively.

Summary statistic of measured laboratory data is provided in Table 1. The range of measured data have a variable respectively CaCO₃ (215 - 3.60 g/kg), P (49 - 1 mg/kg), SOM (67.10 - 3.80 g/kg) and for the variable pH (9.50 - 8), K (0.11 - 4.40), have narrow range in the data set opposite to other properties. Standard deviation (SD) for total N, CaCO₃, P, pH and SOM is respectively 0.59, 39.4, 0.97, 7.85, 0.25, 9.37.

The cross validation approach was used with PLSR model to calibrate soil variables based on entire soil data set. Method was estimated using cross validation R² and RMSE for the soil properties. The result of PLSR model applied to soil properties statistics are shown in Table 2 and Figure 3 for variable soil measured and validated properties. PLSR applied to total N, CaCO₃, total P, CEC, EC, SOM and pH. The results showed that there was strong correlation between Vis-NIR spectra and most of the measured soil properties. Only pH showed low correlation with soil spectra. The best calibration models were obtained for EC (R²-0.96), CaCO₃ (R²-0.94), Total N (R²-0.86), Total P (R²-0.86), CEC (R²-0.83), SOM (R²-0.78). The poor calibration result for pH (R²-0.65) but these properties inherently poorly related with Vis-NIR spectroscopy [1] .

The statistics of cross-validation and result of PLSR models for CaCO₃, CEC, SOM, Total N, Total P, and pH in soils of Mugan plain are presented in Table 2. Standard error for measured vs predicted CaCO₃ (4.22 - 4.11), CEC (0.88 - 0.79), SOM (0.99 - 0.70), Total N (0.06 - 0.05), Total P (0.83 - 0.66) and pH (0.03 - 0.02) respectively.

A correlation between soil spectra and carbonates were investigated by using calibration set and measured soil properties data and there is in 3 region (1800, 2350 and 2360 nm) of spectral wavelength to predict soil carbonates [8] . Prediction result demonstrated that in our prediction set is reliable. The predicted soil samples count are 194 and RMSE = 2.89. The best result for quantifying soil carbonate ranged between 10% - 60% and the

average value is 10.8 (Table 1) in our samples.

Soil organic matter is major constitutes for soil properties and wide spectral range for assess SOM over spectral region. There is strong correlation between soil color and SOM. For this reason there is suggested to use both vis and NIR to get best prediction result [1] . R-square is 0.83 in our sample set and RMSE is 0.54. Bands around 1400, 1890 and 2200 nm are very sensitive in spectra. Therefore, in study area soils have low content of SOM (Table 1.) for this reason Vis regions have high reflectance because of light color. Cation exchange capacity (CEC) prediction result is poor than calibration set. R-square is 0.62 but in calibration set R-square was 0.83 and RMSE = 0.43 prediction quality decreased when using more samples in validation set than calibration set [13] . Generally CEC capacity and SOM demonstrate same spectral similarity [16] .

Table 3 contains statistic result of predicted properties of soil vis-nir spectroscopy measurement. The low result recorded for Total N, pH and CEC (R² 0.44, 0.51 and 0.62). Several scientific works demonstrates that pH hasn’t direct spectral responses. The pH regulated several properties of soil. One of them is CEC, SOM and clay content. Generally study area soils undergo calcification and salinization process. The positive result is CaCO₃ and EC have a good relation with spectral reflection. As a prediction result, R-square of EC is 0.82 and RMSE 0.09.

This paper aimed to predict CaCO₃, CEC, SOM, EC, pH, Total N and Total P content in the arid zone of Azerbaijan using Vis-NIR reflectance Spectroscopy based on PLSR_CV models.

Considering of carbonate, CEC, SOM, EC, pH, Total N and Total P importance in soil also the cost of laboratory measurement of them Vis-NIR is a reliable alternative. It is possible to reduce laboratory costs using nis-NIR spectroscopy. The soil of the study area undergoes different ecological problems such as organic matter losses, salinisation and carbonatation. To study them and to prevent these negative processes the area soils

SE^*: Standart error, SD^#: Standart deviation.

should be analyzed with powerful and reliable methods. PLSR based prediction model with cross validation algorithms were achieved positive result. The result of this study shows that it is possible to predict soil properties with reasonable accuracy using PLSR methods.