Fluvisols Contribution to Water Retention Hydrological Ecosystem Services in Different Floodplain Ecosystems

Water retention is an important hydrological ecosystem service of active floodplain soils. The aim of the study was to evaluate the soil chemical, physical, and hydrological properties in Fluvisols in three different ecosystems that have an impact on water retention hydrological ecosystem services (WRHESs). We selected 16 localities along the Štiavnica River in Central Slovakia, 8 located in riparian zones (RZ), 5 in arable lands (AL), and 3 in grasslands (GL). Soil samples were collected from two layers (0–10 and 20–30 cm). In the laboratory, the soil physical (soil texture) and soil chemical properties (pH, soil organic carbon content, humic and fulvic acid ratio) were determined. Using undisturbed soil samples, the soil physical characteristics (particle density, bulk density, porosity, and actual soil moisture–SMa) were measured. With the help of pedotransfer functions, hydrological soil properties (field water capacity–FWC, wilting point–WP, available water capacity–AWC) were estimated. The recorded properties differed between the localities, ecosystems, and two layers. The SMa values showed a higher soil water retention potential of extensively used ecosystems, such as GL and RZ. However, the hydrological properties estimated by pedotransfer functions (FWC, WP, AWC) showed a higher soil water retention potential in AL localities. This indicated that for calculations, selected pedotransfer functions (particle size fractions, organic matter, and bulk density) and other soil or ecosystem properties (e.g., vegetation cover, meteorological conditions) have an impact on WRHESs. One such soil factor can be the quality of organic matter. On the basis of the results of the ANOVA, significant differences emerged between the different ecosystems for selected basic chemical, physical, and hydrological properties. The effect of the soil layer on the soil properties was revealed only in the case of SOC. The results indicated the effect of different ecosystems on soil WRHES and the importance of extensively managed ecosystems, such as RZ and GL. From this point of view, the reduction in the RZ and GL areas during a period of the last 70 years is negative. The findings should be taken into account in future sustainable floodplain management and landscape architecture.


Introduction
Ecosystems are being studied more in light of the concept of ecosystem services (ESs). ESs are defined as the contributions of ecosystems to human well-being [1]. ESs underpin all aspects of human wellbeing-from basic livelihood to moderate prosperity and sustainable development [2], directly or indirectly contributing to the achievement of sustainable development goals (SDGs) [3]. The 2030 Agenda, adopted by all UN member states in 2015, provides a shared blueprint for sustainability. Its 17 SDGs highlight the pathways to achieve the social, economic, and environmental dimensions of prosperity and sustainability [4]. Embedded in the SDGs are targets to protect, restore, and promote Land 2022, 11, 1510 3 of 20 information [22][23][24]. The most commonly used proxy is land use and land cover data [25]. Biophysical methods can be classified as direct, indirect, and modeling methods. Direct measurement methods include measurements of a state from ecosystem observations and surveys. Direct measurements deliver a biophysical value of ESs in physical units, which correspond to the units of the indicator, and quantify or measure a stock or a flow value. Direct measurements can be used as primary data input to other methods or directly as ES indicators [26]. The availability of accurate soil information for managers or data analysis for modeling is a crucial factor [27]. The use of indices that represent soil quality and their potential to provide ecological functions and ESs is desirable [28].
Soil ecosystem services and their assessment depend on the properties of the soil and their interaction. An assessment of soil HESs focused on soil water retention service can be based on direct measurements of soil hydrological properties, their estimation, or a combination of these two approaches. Although soil properties have been recognized for their importance, a long duration of field work is required to obtain relevant data [29]. Therefore, they are usually estimated through different models or with the help of pedotransfer functions. For WRHESs, assessment can be used as measured or estimated hydrological properties, such as actual soil moisture (SMa), field water capacity (FWC), wilting point (WP), available water capacity (AWC), and soil water deficit index (SWDI). These hydrological properties are mainly related to the soil water retention ecosystem service. Soil moisture forms and the capacity for water in the soil is an important factor for plant growth and against flooding. The projections of climate change indicate greater droughts, which can lead to reductions in soil moisture. Higher values of SMa, FWC, WP and AWC, in terms of volumetric soil water content, indicate a higher soil WRHES capacity. In addition, SWDI positive values mean an excess of water in the soil. The hydrological properties of the soil give a good indication of the potential to provide WRHESs. However, these hydrological properties are dynamic and depend on other soil physical, chemical and biological properties, on land use and management. Therefore, it is crucial to also know the basic soil properties that impact the hydrological properties and their interactions. Knowing these relationships can help with land use and ecosystem management practices.
We evaluated soil WRHES by hydrological soil properties in the floodplain along the Štiavnica River in central Slovakia. The floodplain is relatively small, but, with respect to geographical origin, local measurement of hydrological properties can offer better results than in a large territory [30]. Additionally, the character of the watercourse and its surrounding is diverse due to the great difference in altitude from spring to inflow, different geological and climate conditions, and land cover. The objectives of the study are: (a) to evaluate soil WRHESs with the help of hydrological soil properties measurement and estimation in Fluvisols in three different ecosystems along the Štiavnica River in central Slovakia (b) to measure basic soil chemical and physical properties, and to evaluate mutual interactions between them and the hydrological properties influencing WRHESs, (c) to analyze the effect of ecosystems, land uses, and soil layers on the soil hydrological properties that have impact on soil WRHESs, and (d) to assess changes in land use and land cover in the floodplain with the impact on WRHESs.

Study Area
The study was carried out in Slovakia on the 55 km long Štiavnica River floodplain. We selected 16 localities of which 8 localities were agricultural land (3 of them used as grasslands and 5 used as arable land), and 8 localities were part of riparian zones ( Figure 1). The river flows through the Štiavnické vrchy Mountains of volcanic origin and the Podunajská pahorkatina Hills. The Štiavnické vrchy Mountains are managed as protected areas, the Podunajská pahorkatina Hills are used intensively for agriculture. The lowest locality is located at 123 m above sea level and the highest locality is 648 m above sea level. Due to such a difference in altitude, the watercourse is on the boundary of the Cfb/Dfb regions [31] (using the Köppen climate classification) characterized by a temperate/continental climate. The warm temperate and humid zone Cfb occurred in the Southwestern part of Slovakia with the lowest altitudes. The temperate and humid zone Dfb is located further north at higher altitudes. An average annual precipitation ranges from 584 (the lower part of the river) to 794 (the upper part of the river) mm. An average annual temperature ranges from 7 • C (the upper part of the river) to 10.2 • C (the lower part of the river). Soils are classified as Fluvisol according to the World Reference Base for Soil Resources. Fluvisols developed on fluvial deposits have weak horizon differentiation but distinct 10-25 cm thick topsoil horizons (A horizon) followed by the A/C and C horizon as parent material.
The study was carried out in Slovakia on the 55 km long Štiavnica River floodplain. We selected 16 localities of which 8 localities were agricultural land (3 of them used as grasslands and 5 used as arable land), and 8 localities were part of riparian zones ( Figure  1). The river flows through the Štiavnické vrchy Mountains of volcanic origin and the Podunajská pahorkatina Hills. The Štiavnické vrchy Mountains are managed as protected areas, the Podunajská pahorkatina Hills are used intensively for agriculture. The lowest locality is located at 123 m above sea level and the highest locality is 648 m above sea level. Due to such a difference in altitude, the watercourse is on the boundary of the Cfb/Dfb regions [31] (using the Köppen climate classification) characterized by a temperate/continental climate. The warm temperate and humid zone Cfb occurred in the Southwestern part of Slovakia with the lowest altitudes. The temperate and humid zone Dfb is located further north at higher altitudes. An average annual precipitation ranges from 584 (the lower part of the river) to 794 (the upper part of the river) mm. An average annual temperature ranges from 7 °C (the upper part of the river) to 10.2 °C (the lower part of the river). Soils are classified as Fluvisol according to the World Reference Base for Soil Resources. Fluvisols developed on fluvial deposits have weak horizon differentiation but distinct 10-25 cm thick topsoil horizons (A horizon) followed by the A/C and C horizon as parent material.

Physical Methods and Hydrological Properties Calculations
Physical soil properties were determined using undisturbed soil samples of 100 cm 3 in volume that were collected in September 2018 from a depth of 0-10 and 20-30 cm using core extracting tubes (Eijkelkamp Equipment for Soil Research, The Netherlands). Undisturbed soil samples were hermetically sealed and taken to the laboratory for analyses using the gravimetric method according to Novák [32] of bulk density (BD), particle density (PD), actual soil moisture (SMa), and porosity (as total-P t , capillary-P c , noncapillary-P n , semicapillary-P s ). Raw samples were weighed on laboratory scales. The samples were then saturated with water through filter paper for 24 h. The samples were then placed on dry filter paper and reweighed at certain intervals to determine other characteristics. After half an hour of water desaturation, we determined the moisture for 30 min (SM 30m ). After 2 h of water desaturation, we determined by weighing and recalculating the 2 h of moisture (SM 2h ). After the next 22 h of further desaturation (after a total of 24 h of desaturation), by weighing and recalculating, we determined the 24 h of moisture (SM 24h ). The samples were dried and again weighed. PD as the mass per unit volume of solids in a soil was determined by the pycnometer method according to Blake and Hartage [33]. The following formulas were used for calculations of physical and hydrophysical properties: Bulk density (BD) (Equation (1)): where m d is the mass of dry soil in g, V is the sample volume in core-extracting tube in cm 3 . Particle density (PD) (Equation (2)): where m d is the mass of dry soil in g, m w is the mass of water expelled by soil in cm 3 . Actual soil moisture (SMa) (Equation (3)): where m a is the mass of the soil with actual soil moisture, m d is the mass of the dry soil, BD is bulk density of soil in g cm −3 . Total volume of pores (P t ) (Equation (4)): where PD is the particle density of soil in g cm −3 , BD is bulk density of soil in g cm −3 . Non-capillary pores (P n ) (Equation (5)): where P t is the total volume of pores in %, SM 30 is the moisture at 30 min in %. Capillary pores (P c ) (Equation (6)): where SM 24h is the 24 h moisture in %. Semi-capillary pores (P s ) (Equation (7)): where SM 30m is 30 min moisture in %, SM 24h is 24 h moisture in %. The pedotransfer functions are used to estimate the hydrological properties of the soil as an important water determinant. They are simple to complex knowledge rules that relate available soil information to soil properties and variables that are needed to parameterize soil processes [34]. Many authors have used a fine particle size fraction of less than 0.01 mm to calculate hydrological properties [35][36][37]. Gupta and Larson [38] added in the equation, in addition to particle size fractions, also soil organic matter content and bulk density. By this method, we estimated FWC, WP, and AWC. All these hydrological parameters are crucial for the evaluation of WRHESs focused on water retention. However, it is important to stress that we used a combined approach for hydrological parameter estimation. We used both so-called static (fine particles content and organic matter content) and dynamic (bulk density, actual soil moisture) soil properties. The FWC was introduced to express soil-plant-water relationships as found in field conditions. It is the capacity of the soil to retain water against the downward pull of the force of gravity. It is the point where the soil water holding capacity has reached its maximum for the entire field. At this stage, only micropores or capillary pores are filled with water, and plants absorb water for their use. Part of the water at the field capacity is readily available to plants and microorganisms. The WP is the minimum water content in the soil that the plant needs to not wilt. If the soil water content decreases to this or any lower point, a plant wilts and cannot longer recover its turgidity. The AWC is defined as the difference between the water content in the FWC and the WP. We used the following formulas for the calculations introduced by Gupta and Larson for temperate pedotransfer functions [39]: Field water capacity (FWC) (Equation (8)): Wilting point (WP) (Equation (9)): where Sa is the percentages of sand, Si is the percentage of silt particles, Cl is the percentage of clay particles, SOM is the percentage of soil organic matter, and BD is bulk density. Available water capacity (AWC) (Equation (10)): where FWC is field water capacity in %, WP is the wilting point in %.
The SWDI characterizes drought based on soil moisture and calculated hydrological soil properties. This index is capable of adequately identifying the main attributes that define a drought event. For the calculation, the following formula was used: Soil Water Deficit Index (SWDI) (Equation (11)): where SMa is actual soil moisture in %, FWC is field water capacity in %, AWC is available water capacity in %. The SWDI was proposed by Martinez-Fernandez et al. [40]. In the study, we used the SWDI as proxy indicator. Data used for the calculation does not fulfil standard requirements (e.g., SMa estimation for one day only). Thus, we used this proxy indicator mainly to see differences between ecosystems. Higher values of SMa, FWC, WP, and AWC, in terms of volumetric soil water content, indicate higher soil WRHES capacity. In addition, SWDI positive values mean an excess of water in the soil. When SWDI is positive, soils have excess water; when it equals zero, soil is in the field capacity of the water content (that is, without water deficit). Negative values indicate soil drought, and the water deficit is absolute (wilting point) when the SWDI reaches ≤ −10. At this point, the soil water content is below the lower limit of available water for the plant [41]. The SWDI in the range 0-−2 means that mind, −2-−5 moderate, −5-−10 severe, and ≤−10 extreme drought level.

Chemical Methods
Soil samples were collected in September 2018 from a depth of 0-10 and 20-30 cm to observe differences in soil properties and hydraulic parameters between the upper and lower layers of the soil. Disturbed soil samples were collected for soil analysis. The soil samples were air dried, homogenized, sieved in 2 mm sieves, and the soil texture was analyzed using the pipette method according to Novák [32]. The particle size fractions (sand, silt, clay) were classified according to the United States Department of Agriculture (USDA) system. Sand has a particle size ranging from 0.05 to 2.0 mm, silt from 0.002 to 0.05 mm, and clay less than 0.002 mm. As chemical soil properties, pH was determined in H 2 O in a 1:2.5 ratio and in 1 M KCl solution in a 1:2.5 ratio by the potentiometric method, total soil organic carbon (SOC) using the oxdimetric method according to Tyurin (a modification of Nikitin) [42] which is similar to the Walkley-Black oxidation method. The SOC was oxidized by 0.2 M K 2 Cr 2 O 7 with H 2 SO 4 . After oxidation, excess dichromate was determined by titration with a Mohr salt solution. The content of humic and fulvic Land 2022, 11, 1510 7 of 20 acids was determined using the short fractionation method according to Kononova and Belchikova [43] using extraction with 0.1 M Na 4 P 2 O 7 .10H 2 O. Soil organic matter quality we evaluated by the ratio of humic to fulvic acids (HA/FA). Humic acids are more stable organic matter compounds. The ratio of humic and fulvic acids as indicator of soil quality has been suggested [44]. It is generally accepted that higher values of the HA/FA ratio are typical for more fertile soils with higher stability of organic matter.

Statistical Analysis
To identify the relationships between soil parameters, correlation analyses were performed using Pearson s correlation coefficients and principal component analysis (PCA). Soil parameters were analyzed using analysis of variance (ANOVA). All statistical tests were performed in the SPSS Statistics software (version 19.0).

Remote Sensing Data
To describe land use and land cover changes (LULC) in this study area, two dates with their respective LULC maps were combined: 1950 and 2019. The LULC map for 1950 is a product of the Czechoslovakian military topography mapping. The 2019 map was generated from satellite images. The following categories were analyzed: (1) forest areas, (2) riparian vegetation, (3) other nonforest woody vegetation, (4) mosaic structures, (5) permanent grasslands, (6) arable land, (7) permanent crops, (8) industrial units and industrial parks, (9) mineral extraction sites, (10) other artificial surfaces, and (11) water bodies and areas.

Soil Chemical and Physical Properties
The basic statistical characteristics of the chemical and physical properties of soils at the RZ, GL, and AL localities in two depths are reported in Tables 1-3. The measured properties showed their heterogeneity.   In general, soil pH/H 2 O ranged from 5.04 to 7.65, and pH/KCl from 4.13 to 7.20. Mean pH values were higher in the AL localities in both layers compared to the RZ and GL localities. Overall, the SOC ranged from 9.00 to 30.60 g kg −1 with higher content in the upper soil layers and at the RZ localities followed by the GL and AL localities. However, the quality of the organic matter expressed by the HA/FA ratio was better in the AL localities (expressed by higher values of the HA/FA ratio). In AL localities, the mean HA/FA ratio was higher in the upper soil layer, while in the RZ and GL localities it was the opposite (Table 1).
In general, the average clay content ranged from 4.25 to 42.28%. The lowest clay content was observed in the GL localities in both layers (range of 4.25 to 7.76%), followed by the RZ (range of 4.35 to 15.22%) and AL (range of 8.50 to 42.28%) localities (Table 2).
In general, PD ranged from 2.04 to 2.67 g cm −3 and BD from 0.64 to 1.54 g cm −3 . Mean BD values were higher at AL localities at both depths compared to RZ and GL localities. The mean P t value was highest at the RZ localities in both soil layers followed by the GL and AL localities. The mean P t value was higher in the upper layer at RZ and GL localities, and, conversely, the mean P n value was lower in the upper layer at both ecosystem types. The mean value of P c was higher in the upper layer in all ecosystems. The mean P s value was lowest in the AL localities compared to the RZ and GL localities. (Table 3).

Soil Hydrological Properties Evaluation in Relation to Soil Water Retention Hydrological Ecosystem Services
The soil hydrological properties were estimated along the Štiavnica River for three different ecosystems ( The differences in soil hydrological properties between two soil layers within one ecosystem type were not very distinctive. The hydrological properties estimated with the help of pedotransfer functions (FWC, WP, AWC) showed a higher soil water retention potential at AL localities. On the contrary, the SMa and SWDI values showed a higher soil water retention potential in extensively used ecosystems, such as GL and RZ. This finding indicates that in addition to the soil, other factors are involved in the fulfillment of WRHESs.

Relationships between Soil Properties, Ecosystems, Land Uuse, and Soil Layers
Many significant correlations between soil properties were established (Table 4). Estimated hydrological properties were significantly correlated mainly with particle size fractions, the HA/FA ratio as an indicator of soil organic matter quality, and semi-capillary porosity. Correlations with pH, BD, total porosity, and capillary porosity were also recorded. We recorded no correlations between hydrological properties and SOC or noncapillary porosity.
Principal component analysis (PCA) was used to observe relationships between variables (Table 5, Figure 3). The first two factors can explain more than 60% of the total variance in the all and RZ localities, more than 70% in the GL localities, and more than 80% in the AL localities. The FWC, WP, and AWC, clearly strongly positively correlated with factor no. 1 at all, RZ, GL localities (more than 0.9), and AL localities (besides AWC). In contrast, at AL localities, FWC and WP clearly strongly negatively correlated with factor no. 1. SMa was strongly positively correlated with factor no. 1 at the AL localities (0.989). At the RZ and all localities, SMa was positively correlated with factor no. 2. in the AL localities. The FWC, WP, and AWC, clearly strongly positively correlated with factor no. 1 at all, RZ, GL localities (more than 0.9), and AL localities (besides AWC). In contrast, at AL localities, FWC and WP clearly strongly negatively correlated with factor no. 1. SMa was strongly positively correlated with factor no. 1 at the AL localities (0.989).
At the RZ and all localities, SMa was positively correlated with factor no. 2.     Based on the results of the ANOVA, significant differences emerged between the different ecosystems for pH/H 2 O, clay, silt, sand, BD, P t , P s , FWC, WP, and SWDI ( Table 6). The effect of the soil layer on the soil properties was revealed only in the case of SOC. The p-values shown in bold are statistically significant at p < 0.05.

Land Use and Land Cover Changes in Floodplain and Their Impact on Hydrological Ecosystem Services
The LULC composition of the floodplain suffered changes during a period of 70 years (Figure 4, Table 7). In 1950, the landscape was predominantly devoted to agricultural activities. Arable land accounted for 37.2% and permanent grasslands for 26.3% of the study area. The forest areas covered 25.6% and riparian vegetation only 1.4% of the study area. In 2019, due to the management decision, forest areas increased to a cover of 32.5%. The increase was recorded for other non-forest woody vegetation, permanent crops, industrial, and artificial surfaces. Arable land and permanent grasslands decreased to 31.4% and 16.9%, respectively. Riparian vegetation decreased to cover 1.1%. increase was recorded for other non-forest woody vegetation, permanent crops, industrial, and artificial surfaces. Arable land and permanent grasslands decreased to 31.4% and 16.9%, respectively. Riparian vegetation decreased to cover 1.1%.

Soil Hydrological Properties and Their Relationships with Soil Properties
Fluvisols located on the Štiavnica River floodplain showed heterogeneity in soil properties, including hydrological properties, which is in line with many studies [45,46]. Soil hydrological characteristics belong to the most important factors for plant growth, microbial activity, and nutrient cycling [47]. The estimated hydrological properties are important parameters used to characterize the capacity of soil to provide water for plant growth. They are useful for many ecosystem and agronomic purposes and for measuring relative differences among soils. They represent a practical and easily understandable coefficient for measuring the ability of a soil to store and retain water, and thus contribute to HESs. Higher values of hydrological property mean higher soil contribution to water retention HESs.
The texture of the soil and the content of organic matter play an important role in the physical and hydrological properties of the soil [48]. We found correlations between hydrological parameters and selected soil characteristics. The proportion of sand, silt, and clay and the amount of organic matter present in the soil is important [49,50]. The textural composition differed strongly between the localities, and the clay content varied between 4.25 and 42.28%. The SOC content was generally higher for topsoil (from 10.65 to 30.60 g kg −1 ) than for subsoil (from 9.00 to 22.95 g kg −1 ). The combination of soil organic matter amount with soil texture affects the bound water held in soil matrix either via electrostatic forces at adsorption sites and/or surface tension (capillary forces) in soil pores.
Water storage and retention are directly related to pore volume and its categories [51]. In the study, expressed as mean values, the RZ localities showed a higher total porosity in both soil layers and capillary porosity in the second soil layer. In the case of AL localities, the average Pc value was higher in the first soil layer, indicating the positive impact of proper agricultural soil management on hydrophysical soil properties in the surface soil layer. The semi-capillary porosity expressed by the average value was significantly higher at the RZ localities compared to AL in both soil layers. The proportion of pores is also reflected in the actual soil moisture content. Higher values were observed at the RZ localities in both soil layers. The forms of extensive land use, such as forests, show higher capillary porosity compared to agricultural land [52], which has a positive impact on the soil water regime. Porosity and other soil properties reflected in the values of hydrological parameters.
In the study, FWC, WP, and AWC were significantly correlated with particle size fractions, which is a logical consequence of using the soil texture in the calculation. However, we did not find a correlation between hydrological properties and the amount of organic carbon in the soil. Soil organic carbon can improve soil quality by improving soil structure and improving soil moisture regime. Therefore, SOC is also suggested as a means of drought management. Yet, there is still no clear consequence on its quantitative effect [53]. Libohova et al. [54] stated that there have been overestimates of the amount of water that soil organic matter increase can contribute to the available water holding capacity. However, we found a correlation between hydrological properties and organic matter quality expressed by the HA/FA ratio. Humic acids are typically comprised of aliphatic and aromatic moieties with multiple functional groups (e.g., -OH, -COOH), which enables their adsorption onto solids [55]. They are able to withstand degradation for long periods. Humic acids play several important roles, such as an increase in soil physical and biochemical activities by improving soil texture and water holding capacity [56][57][58].
The results also indicated that physical properties and hydrological parameters are controlled not only by the soil constituents, such as mineral and organic particles but also by other factors. The pH values of soil in the natural environment have an enormous influence on the biogeochemical processes of the soil [59]. Therefore, soil pH is described as the master soil variable that influences many biological, chemical, and physical properties and processes of the soil that can affect the availability of water for plant growth and biomass yield [60].

Hydrological Properties and Water Retention Hydrological Ecosystem Services
Hydrological properties are the most important in monitoring the water cycle and drought [35]. Higher values of hydrological properties mean higher soil capacity for water storage and retention HESs. Using pedotransfer functions as soil texture, SOC, and BD we estimated higher average values of FWC, WP, and AWC at AL localities compared to RZ and GL localities. On the contrary, taking into account also the actual soil moisture content for the SWDI calculation, we found higher values for SWDI at GL and RZ localities. This indicates the importance of other factors, including organic matter quality or vegetation cover, having an impact on water retention in soil.
The average proxy SWDI values were negative in all three ecosystems which should indicate soil drought. The SWDI is considered a promising approach to use soil physical characteristics to measure available soil moisture in the root zone [61,62]. It has a greater biophysical meaning than other vegetation indices or climatic variables, and even some methods based on soil moisture [63].
In this regard, the study results can indicate further possibilities on how to improve the WRHESs capacity. The addition of inorganic or organic additives in different forms to the soil can improve WRHESs. Traditional materials (farmyard manure) have long been used in agriculture, having a positive effect on water holding and retention capacity [64]. In recent years, numerous approaches have been proposed to alleviate soil water holding capacity, including biochar application [65] and superabsorbent polymers [66]. Mi et al. [67] proved that bentonite application would have beneficial effects on the soil water holding properties. However, these practices can not only have a positive effect on hydrophysical properties, as in the case of biochar addition [68], but also a negative effect, as in the case of pig slurry use, which supports soil hydrophobicity [69]. Therefore, for optimal land management in relation to WRHESs, various aspects should be considered and taken into account based on the knowledge and evaluation of basic physical and hydrological properties. In particular, sustainable management based on nature-based solutions can simultaneously provide environmental, social, and economic benefit [70] and provide diverse ecosystem services.

Impact of Ecosystems and Land Use on Soil Hydrological Properties and Services
The ANOVA revealed an effect of different ecosystems on some physical and hydrological properties. There are studies that confirmed the effect of ecosystems and land use and cover on the textural composition of Fluvisols, the bulk density, porosity, and the quality of the humus and thus the effect on the fulfillment of WRHESs [71][72][73][74].
These results emphasized the importance of riparian ecosystems and grasslands in floodplains. Mostly RZ provide high values of evapotranspiration and contribute to the fulfillment of HESs [75]. In addition, Felipe-Lucia and Comin [76] stressed the importance of the mosaic character of the landscape, which can support HESs. The character of the mosaic and its heterogeneity are guarantees of stability and sustainability. Maintenance of agricultural land in floodplains, together with riparian and grasslands ecosystems, extensive or proper intensive management of agricultural land, could be an ecological solution to maintain a balanced water budget and HESs. Understanding the relationships between soil parameters that affect WRHESs is of importance for sustainable floodplain management [77]. Furthermore, WRHESs can be considered as suitable indicators for assessing ecological stability and resilience [78]. Therefore, management and regulation strategies regarding these ecosystems are fundamental to the provision of WRHESs on a regional scale [79]. However, these facts are often ignored [80], and severe WRHESs losses in floodplains were observed when provisioning services were strengthened [81], often at the expense of RZ localities. In our study we also confirmed the decrease of areas of riparian vegetation and grasslands. In this case, the extension of riparian zones and grasslands in the future would be desirable. The increase in forest area was positive because afforestation can have many positive impacts, including carbon sequestration, increasing soil organic carbon, improving soil structure, improving the landscape and biodiversity. Compared to an extensive form of management, intensive arable cropping has a large impact on soil properties and their dynamic due to frequent soil disturbance. Although the arable land area decreased during 70 years in the observed territory, its homogenous character is not suitable for sustainable management. In addition, land management practices must be adopted in agriculture to ensure both provisioning and hydrological ecosystem services.

Conclusions
We found heterogeneity in the chemical and physical soil properties of the soil in the floodplain of the Štiavnica River reflected in the hydrological properties. Higher values of hydrological properties meant higher soil capacity for WRHESs. Using pedotransfer functions (soil texture, SOC, and BD), we estimated higher average values of FWC, WP, and AWC at AL localities compared to RZ and GL localities. On the contrary, also taking into account the actual soil moisture content, we found higher values for proxy SWDI at GL and RZ localities. This indicates the importance of other soil and ecosystem factors, including organic matter quality or vegetation cover, that impact WRHESs. In addition, the ANOVA results confirmed an effect of different ecosystems on most observed soil chemical, physical, and hydrological properties. The results indicated the effect of different ecosystems on soil WRHES and the importance of extensively managed ecosystems, such as RZ and GL. Their areas in the floodplain of the Štiavnica River were reduced during a period of the last 70 years. The extensively used ecosystems are irreplaceable in providing WRHESs in floodplains, and their areas should be extended in the future. Such results should be taken into account in land use management and landscape architecture. Floodplains provide multiple ecosystem services, and WRHES should not be underestimated because of many indirect effects on all other services.