Examination of Susceptibility to the Deﬁciency of Soil Water in a Forested Agricultural Area

: Mountainous regions present numerous obstacles to agriculture. These include the terrain, which is associated with surface erosion, as well as surface runoff, which washes away plant nutrients and weak soil. Spatial analysis is currently used in the study of various stochastic variables, especially those of high priority for soil water properties. Small watershed and basin-scale models were used to simulate the quantity of surface run-off, groundwater and predict the environmental impact of land use and land management practices. A new generation of the distributed hydrological models has greatly broadened simulation ﬁelds to soil and water diversiﬁed situations. The study also measured declines in slope and grain size distribution, factors impacting surface erosion and surface runoff. Multivariate statistics (canonical analysis) showed that soil moisture was most correlated both with agricultural land and forests, which is why it was used to create the model of spatial distribution. The model showed that salinity has the smallest forecast error in modeling, and thus best corresponds with the soil moisture. It is important to make a correct diagnosis of soil properties, and the degree of degradation. The assessment of the physiographic parameters of a basin will contribute to the development of proper usage and determine the quality of the water in the soil, which will be essential for forest resources and agricultural land in mountain areas exposed to surface erosion.


Introduction
In soil erosion research, it is important to determine the spatial variability (diversity) and spatial continuity (similarities). One uses the assumption that tested objects found closer together are more similar than objects found further apart [1]. A geographic information system (GIS) technique is commonly used in assessing the dynamics of rivers and the physiographic characteristics of a basin [2]. The evaluation of the quality of river basins is important in environmental monitoring [3]. The continuous spatial data of soil are difficult to present uniformly, and therefore new techniques and models are constantly being sought out, while results are verified using calibration [4]. Most frequently, multivariate statistics are used to assess the quality of river basins [5]. In environmental monitoring, the soil should be examined as well as other physical and chemical parameters.
Soil erosion negatively impacts plants by reducing the amount of available water [6][7][8]. It also has an impact on soil microorganisms. That is why the continuous control of the spatial distribution of soil parameters is substantial. Land use, next to soil moisture and precipitation, is one of the most valuable parameters in long-term studies of soil erosion [9][10][11]. For the spatial analysis of a basin, it is easiest to use a digital model of the terrain, on the basis of which rill erosion is assessed, among other things [12]. Monitoring the temperature in the soil contributes to a better understanding of the water content in the soil, especially for rainwater infiltration [13]. Changes in temperature affect porosity and Earth 2021, 2 533 soil particles. In clay soils predominantly consisting of loam, the temperature influences the interaction between the solid phase and liquid phase densities [14]. Temperature plays an important role in plant germination, root growth and respiration, and in the decomposition and mineralization of organic matter. It was also insightful in the weathering of cationic bases, the assimilation of nutrients, and the increase in crop yield times [15]. The spatial distribution for the mean temperature and volumetric soil moisture was noteworthy on arable land. The spatial monitoring and evaluation of EC (electric conductivity) are easiest to carry out using the method of telemetry and remote sensing of the environment. The spatial distribution of EC conducted by geostatistical analysis is very well suited for the modeling and detailed mapping of the soil [16].
The solution proposed in the article is a combination of multivariate statistical methods and a spatial model. The study of spatial properties is based on the interpolation of measuring points. Data interpolation should be performed for objects tested in multiple transects [17]. Hitherto, conventional kriging has been used to predict the content of organic matter in soil [18] and the spatial distribution of total carbon trapped in carbonate rocks [19]. The work uses geostatistical methods for the soil layer from 0 to 20 cm. Electrical conductivity showed greater recognition as high soil salinity. The salinity threshold of most cultivated vegetables is low, and salt tolerance is reduced when saltwater is used in crop irrigation techniques [20]. The latest research on a general predictive model based on the relationship between electrolyte concentration and negative clay soil particles is useful for determining the degree of soil dispersion. The study of soil salinity showed that monitoring of the top layers (0-15 and 15-30 cm) is possible with the use of GIS methods showing spatial differentiation [21].
The article presents the results of measuring soil resistance to penetration, which to some extent showed soil susceptibility to various forms of degradation. It is important from the point of view of irrigation because the hydraulic conductivity decreases due to the deterioration of the soil structure [22]. The aim of this study was to show the difference in soil water retention for small-scale and medium-size mountainous catchments.

Study Sites
The Smugawka ( Figure 1) and Mszanka ( Figure 2) creeks basins are situated in the Outer Western Carpathians and more specifically in Beskid Wyspowy mezoregion (Poland). The focus area is dominated by meadows, farmland, and forests. The dominant slopes in the region range from 10 to 20%. Four groups of soil texture were distinguished in the area of research. The farmland is made up of similar proportions of sandy clay loam and silt loam. Forest areas are dominated by loam, while clay loam appeared only on a small area of research. The bulk density throughout the both basins had a range from 1.42 to 1.72 g cm −1 .

Sampling and Data Handling
Measurements of soil resistance were carried out in three repetitions at one poi arable lands and grasslands at the 0-80 cm layer and in forest soil at the 0-20 cm l during subsequent field visits. In total, 27 samples were obtained for each type of use

Sampling and Data Handling
Measurements of soil resistance were carried out in three repetitions at one point in arable lands and grasslands at the 0-80 cm layer and in forest soil at the 0-20 cm layer, during subsequent field visits. In total, 27 samples were obtained for each type of use. The measurements were made using an electronic penetrologist with automatic control of the penetration speed by Eijkelkamp. The soil volume density with retained structure was calculated using the cylinder method with a cylinder volume of 100 cm 3 . The soil mass after drying at 105 • C was determined as well as the soil volume density (ρ o )-in the natural state according to the following relationship: where: m-soil dry mass (Mg), v-volume of soil in the natural state (m 3 ). Soil salinity (EC), moisture and temperature samples were performed at 80 spots by means of a test probe (type-HH2 moisture measure) using the TDR method. This method allows for noninvasive, accurate, and fully automatic measurements of electrical conductivity, moisture and temperature in the soil. For the analysis of particle size, averaged samples were taken (25 samples) at each test area. Each of the test samples was taken from the top layer of the soil (0-25 cm). The collected samples were dried, and then passed through a sieve with an aperture of 2 mm at room temperature. Determination of soil texture was made using Casagrande's aerometric method with Prószyński's modification. Multivariate statistics, including canonical analysis (CA) and nonmetric multidimensional scaling based on Euclidean distances (n-MDS), were performed with PAST software ver. 3.12.
A completely revised version of the hydrological model SWAT+ (Soil and Water Assessment Tool) was applied to show spatial distribution of the soil water and hydrological properties. Soil resistance, electrical conductivity, volumetric moisture, hydraulic conductivity and temperature in the soil layer from 0 to 20 cm were shown. Data collected in 2014-2018 were used. The results were averaged and presented in the form of spatial distribution for the entire catchment area. SWAT+ was aimed at optimizing spatial information by selecting the right algorithms. The method used was useful in forecasting measurement points with a similar spatial distribution of the examined features and in selecting those variables with the greatest strength of connection with surface runoff. The maps were generated using the QGIS version 3.16, intended for mapping and modeling the terrain surface. In order to cope with natural disasters, such as floods and droughts, and meet the growing needs of society, it is necessary to develop water resources for irrigation and drinking water supply. SWAT+ is a hydrological model that works over a period of time. The study used data from 1991-2020 to show the results. Validation was also performed for a six-year period (2014-2018) to verify the pool response using calibrated fitted values. In addition, the performance and evaluation of the model were analyzed using statistical parameters such as the correlation coefficient, Nash−Sutcliffe efficiency and mean square deviation coefficient. Additionally, SWAT+ was used to assess the impact of climate change and anthropogenic factors on the flow of water in the aquifer and the outflow of main waters to the canal and the main river.

Spatial Distribution of Soil Properties
The EC of the soil in the both basin areas ranged from 10-20 mS m −1 (outside the basin area up to 80 mS m −1 ). The lowest value was recorded in the meadows, while the highest value of salinity appeared in the soils beneath the agricultural area. Unfavorable climatic conditions in mountainous regions, e.g., temperature inversion, which contributes to delayed vegetation in the valleys was the reason to show the spatial variability of temperature. The salt's direction of migration depends on the amount of water contained in the soil, as well as weather and climate conditions. This was the reason for the designation of the spatial variability of electrical conductivity (EC; salinity), EC was selected for spatial analysis because salt reduces the number of nutrients and temperature, which inhibits the growth of plant roots. The EC level was compared between agricultural areas (average 14.54 mS/m for Smugawka catchment) and forested areas (average 1.64 mS/m for Smugawka catchment) in the basin of the mountain.
Statistical analysis showed that throughout the study area, the estimated salinity of the basin prevailed at the level of 6 mS m −1 . In general, agricultural soils were characterized by the highest rate of salinity. The calculations were performed for the basin area equal to 7.75 km 2 and 54.5 km 2 for the range of autocorrelation determined using the exponential semivariogram of 2.5 km. The mean values of the soil temperature ranged from 14.38 • C to 17.36 • C for Smugawka catchment ( Table 1). The average temperature value comes out at 15.84 • C for the Mszanka catchment ( Table 2). An estimate of the temperature for the entire area made using the multivariate statistics showed that the most likely temperature of the soil, which is to be expected throughout the basin area is 16.5 • C.

Impact of Soil Properties on Land Use
A selection of parameters descriptive of the soils which are related to both agricultural lands and forests was made by correspondence analysis (CA). Soil moisture demonstrated the closest relationship with salinity (EC). Soil temperature correlated positively with axis 1 and negatively with axis 2. Clay and particulate matter in the soils had the greatest impact on the soil temperature. Soil moisture correlated negatively in both the agricultural land and the forests (Figure 3). It should be clarified which factors most determine the outcome and make the choice of soil parameters, and which have the lowest forecast error in determining the moisture.
The total variance for the two axes was close to 89.5%. In order to determine the gradation of the effect of all studied factors on the distribution of moisture in the basin area, multivariate statistics were used. Correspondence analysis (CA) allows for the specification of which parameters affect the size of the test characteristics to the greatest extent, identify, and prioritize the order of the factors, from those with the greatest impact to those with the least. In the case of the basin, land use had the least impact on the amount of soil moisture. Then came bulk density. The parameters with the greatest impact on soil moisture proved to be soil texture, temperature, and EC. The examined soil parameters for meadows and forests were positively correlated for both catchments. The soil data collected for the agricultural areas were negatively correlated. The best type of correspondence to describe the basins proved to be the squares representing sets of points related to forests (Figure 4).

Impact of Soil Properties on Land Use
A selection of parameters descriptive of the soils which are related to both agricu tural lands and forests was made by correspondence analysis (CA). Soil moisture demon strated the closest relationship with salinity (EC). Soil temperature correlated positivel with axis 1 and negatively with axis 2. Clay and particulate matter in the soils had th greatest impact on the soil temperature. Soil moisture correlated negatively in both th agricultural land and the forests (Figure 3). It should be clarified which factors most de termine the outcome and make the choice of soil parameters, and which have the lowe forecast error in determining the moisture. The total variance for the two axes was close to 89.5%. In order to determine the gra dation of the effect of all studied factors on the distribution of moisture in the basin are  Earth 2021, 2, FOR PEER REVIEW multivariate statistics were used. Correspondence analysis (CA) allows for the spe tion of which parameters affect the size of the test characteristics to the greatest e identify, and prioritize the order of the factors, from those with the greatest imp those with the least. In the case of the basin, land use had the least impact on the am of soil moisture. Then came bulk density. The parameters with the greatest impact o moisture proved to be soil texture, temperature, and EC. The examined soil param for meadows and forests were positively correlated for both catchments. The soi collected for the agricultural areas were negatively correlated. The best type of spondence to describe the basins proved to be the squares representing sets of poi lated to forests ( Figure 4). The 1991-2020 simulation showed that the runoff of surface waters is very het neous ( Figure 5). This study synthesizes the methods for estimating water moveme storage measured in a set of forested agricultural areas. Moreover, the highest annu face runoff in the Mszanka catchment was on average from 5460 mm to 7640 mm ( 6). Problems are often found with parameter estimation of different space−time Finally, the visual inspection of the water level response to water scarcity greatl moted the hydrological model in this research. The research revealed that the level of the Smugawka aquifer accumulated rain in the range from 84 to 141 mm (Figure 7), and the level of aquifer of the Mszanka ment area was from 348.6 to 357.1 mm (Figure 8). The average distance from the table for Smugawka was 7.1-8.3 m (Figure 9). For the Mszanka catchment, the wate depth was at 2.82-2.99 m ( Figure 10).   The 1991-2020 simulation showed that the runoff of surface waters is very heterogeneous ( Figure 5). This study synthesizes the methods for estimating water movement and storage measured in a set of forested agricultural areas. Moreover, the highest annual surface runoff in the Mszanka catchment was on average from 5460 mm to 7640 mm ( Figure 6). Problems are often found with parameter estimation of different space−time scales. Finally, the visual inspection of the water level response to water scarcity greatly promoted the hydrological model in this research.  The research revealed that the level of the Smugawka aquifer accumulated rainwater in the range from 84 to 141 mm (Figure 7), and the level of aquifer of the Mszanka catchment area was from 348.6 to 357.1 mm (Figure 8). The average distance from the water table for Smugawka was 7.1-8.3 m (Figure 9). For the Mszanka catchment, the water table depth was at 2.82-2.99 m ( Figure 10).
The results clearly showed that the larger catchment area is characterized by useful hydrological parameters for soil water retention. In order to compare the type of soil use on the selected parameters of soil, nonmetric multidimensional scaling based on Euclidean distances were used. The analysis showed that the parameters of forest soil were similar for both grassland and agricultural land. However, meadows were more closely tied to forests, as evidenced by the positive correlation between the measured parameters of the soil (Figure 3). It can therefore be concluded that the assessment of the quality of soil for meadows and forests will be similar because they display similar properties. For agricultural land, other factors, environmental variables, and parameters in the assessment of soil should be considered as well. For all the presented geostatistical methods, a representative sample selection is required, especially for the runoff trends analysis and future projections of hydrological patterns in small forested catchments [23], mathematical modeling [24], land use [25], spatiotemporal prediction of root zone soil moisture [26], the spatial distribution of the measured features [27][28][29], yield measure [30], integrated monitoring system of water quality [31], agricultural water management [32,33], and an appropriate interpretation of spatial statistics [34].
Earth 2021, 2, FOR PEER REVIEW 9 . Figure 6. Surface runoff for Mszanka catchment in simulation period 1991-2020.     The results clearly showed that the larger catchment area is characterized by useful hydrological parameters for soil water retention. In order to compare the type of soil use on the selected parameters of soil, nonmetric multidimensional scaling based on Euclidean distances were used. The analysis showed that the parameters of forest soil were similar for both grassland and agricultural land. However, meadows were more closely tied to forests, as evidenced by the positive correlation between the measured parameters of

Conclusions
Multivariate analysis can be used to identify factors affecting the amount of land use. Soil parameters such as EC, temperature, moisture, bulk density, and soil texture are suitable for the evaluation of soil quality for land use. Both EC and soil temperature are important parameters that significantly affect the soil moisture, yield, and productivity of soils. By analyzing the level of EC of the soil in the studied basins, it was found that the smallest salinity content was present in the forest, and the largest in the agriculture areas. This can affect the formation of significant differences between forest areas, which in turn creates the need for additional agronomic work in the areas around the basin located near the forest. The differences between the factors affecting land use become apparent from the water analyzed in the soil. Our study applies a methodology for the soil water examination that provides useful information for water resource assessment at catchment scale. In order to demonstrate the productive suitability of the soils, and to improve water accessibility that is important for agriculture, a lot of research should be done in these areas to obtain the data on the variation in the soil properties.