Use of Pedotransfer Functions in the Rosetta Model to Determine Saturated Hydraulic Conductivity (Ks) of Arable Soils: A Case Study

: A key parameter for the design of soil drainage and irrigation facilities and for the modelling of surface runoff and erosion phenomena in land-formed areas is the saturated hydraulic conductivity (Ks). There are many methods for determining its value. In situ and laboratory measurements are commonly regarded as the most accurate and direct methods; however, they are costly and time-consuming. Alternatives can be found in the increasingly popular models of pedotransfer functions (PTFs), which can be used for rapid determination of soil hydrophysical parameters. This study presents an analysis of the Ks values obtained from in situ measurements conducted using a double-ring inﬁltrometer (DRI). The measurements were conducted using a laboratory permeability meter (LPM) and were estimated using ﬁve PTFs in the Rosetta program, based on easily accessible input data, i.e., the soil type, content of various grain sizes in %, density, and water content at 2.5 and 4.2 pF, respectively. The degrees of matching between the results from the PTF models and the values obtained from the in situ and laboratory measurements were investigated based on the root-mean-square deviation ( RMSD ), Nash–Sutcliffe efﬁciency ( NSE ), and determination coefﬁcient (R2). The statistical relationships between the tested variables tested were conﬁrmed using Spearman’s rank correlation coefﬁcient (rho). Data analysis showed that in situ measurements of Ks were only signiﬁcantly correlated with the laboratory tests conducted on intact samples; the values obtained in situ were much higher. The high sensitivity of Ks to biotic and abiotic factors, especially in the upper soil horizons, did not allow for a satisfactory match between the values from the in situ measurements and those obtained from the PTFs. In contrast, the laboratory measurements, showed a signiﬁcant correlation with the Ks values, as estimated by the models PTF-2 to PTF-5; the best match was found for PTF-2.


Introduction
Water permeability is a key property of soils, especially with regard to the design of soil irrigation and drainage facilities, modelling of surface runoff and erosion phenomena in land-formed areas, and environmental processes occurring in porous media [1][2][3]. The study of soil physics, and in particular the determination of the filtration coefficient (Ks) based on direct methods, is a very interesting issue; nevertheless, it is both labour-intensive and costly [4][5][6][7][8][9].
Climate change in a region or environment entail changes in the method of soil cultivation and plant production. Changes in the interactions between agriculture (i.e., soil compaction) and natural environment (i.e., weather, soil conditions) are a key feature of the transitions which scientists are trying to explain. Therefore, learning the landscape is an important tool in the proper management of water resources in rural areas. Agricultural In this work concentrated on the full depth soil pits (150 cm) and including an estimate of Ks for surface layers (0-50 cm) and additionally for a deeper levels of soil (50-150 cm). It was considered that for engineering purposes, such as irrigation or drainage planning, it is advisable to know the Ks parameter.
The aim of this research was to test the possibility of using PTFs in the Rosetta program to determine the Ks values of arable soils based on easily accessible input data, i.e., the soil type, granulometric composition, density, and water content. We hypothesize that the Ks values obtained from PTFs in the Rosetta program for arable soils correspond to the values obtained from (i) direct in situ and (ii) laboratory measurements. To test the hypothesis we examined the arable soils in Central Europe (district of Racibórz, south Poland) in a 150 cm depth soil pits.

Description of Study Area
The study area is located in the southwestern part of the Silesian Voivodeship (Poland) in the district of Racibórz (Figure 1). According to Kondracki's [39] division of Poland into physio-geographical regions, the study area is located in the Central European Plain (31), in the macroregion of the Silesian Lowlands (318. 5), and on the border of two mesoregions: the Głubczyce Plateau (318.58) and Racibórz Basin (318.59). The main watercourse flowing through the region is the Oder River, along with its left tributary, the Psina River.

Meteorological Conditions
In terms of climate, the study area is considered as one of the warmest areas in this region. In general, Poland has a mostly temperate climate, in transition between an oceanic climate dominating in the north and west of the country, and a continental climate in the south and east. In the multiannual period from 1971-2000, the mean annual air temperature was 8.5 °C, and total precipitation was 616 mm (according to the Institute of Meteorology and Water Management-station in Racibórz). Higher mean temperatures and a decrease in rainfall have been observed over the last decade. The characteristics of the meteorological conditions are listed in Table 1.
To assess the meteorological conditions of the region where the study was conducted, Selyaninov's hydrothermal coefficient (HTC) values were calculated, as a meas-

Meteorological Conditions
In terms of climate, the study area is considered as one of the warmest areas in this region. In general, Poland has a mostly temperate climate, in transition between an oceanic climate dominating in the north and west of the country, and a continental climate in the south and east. In the multiannual period from 1971-2000, the mean annual air temperature was 8.5 • C, and total precipitation was 616 mm (according to the Institute of Meteorology and Water Management-station in Racibórz). Higher mean temperatures and a decrease in rainfall have been observed over the last decade. The characteristics of the meteorological conditions are listed in Table 1. To assess the meteorological conditions of the region where the study was conducted, Selyaninov's hydrothermal coefficient (HTC) values were calculated, as a measure of the precipitation efficiency in a given month. The calculation was performed as follows [40]: Land 2021, 10, 959 5 of 21 In the above, P is the monthly sum of precipitation (mm), and Σt is the sum of the mean daily air temperature values in a given month (in • C).
Based on the HTC values calculated for the months of the growing season from the 2010-2019, and as shown in Table 1, it was determined that 58.4% of the months were dry (8.3% extremely dry, 13.3% very dry, 28.5% dry, and 8.3% relatively dry); in 13.3% of the months, the conditions were optimal; and 28.3% of the months were humid (8.3% moderately humid, 6.7% humid, 8.3% very humid, and 5.0% extremely humid). In general, the variable meteorological conditions in the study area, characterised by periods with excesses and shortages of water, indicate the need to use drainage and irrigation facilities for agriculture. To determine their technical parameters, information on the soil permeability is needed. This confirms the advisability of conducting research on the possibility of using mathematical models to determine the soil Ks values, as an alternative to time-consuming in situ or laboratory tests.

Field Measurement and Soil Sampling
The research was conducted on arable lands during the agricultural season from 2012 to 2015. The field soil tests were conducted at 16 measurement-control points (up to a depth of 150 cm) in Wojnowice (2 points), Bojanów (2 points), Owsiszcze (6 points), Strzybnik (2 points), and Tworków (4 points). Three to five genetic levels were identified for each soil profile.
The soil infiltration was measured in situ, at depths of 10 cm (topsoil) and 35 cm (subsoil), using a double-ring infiltrometer (DRI) method ( Figure 2). The DRI comprised an inner ring (9.5 cm diameter) and outer ring (19.5 cm diameter) inserted into the ground at a depth of 10 cm. The DRI was inserted by using a falling weight-type hammer striking on a wooden plank placed uniformly on top of the ring, and without undue disturbance to the soil surface. Each ring of the DRI was filled with a constant head of water level, and the outer ring helped when checking the lateral flow from the inside ring, so as to better estimate the infiltration, reducing losses. The Ks value was estimated when the water flow rate inside the inner ring reached a steady state [1], which in the case of the studied soils, lasted approximately 3-4 h. The infiltration rate was calculated for the respective time intervals ∆T as follows [41]: where V is the volume of water (cm 3 ) added to the inner ring at time ∆T (s), and D r is the diameter of the inner ring (cm 2 ). For steady infiltration (as a constant value for the soil), the Land and Water Development Division [42] developed infiltration classes, as follows: very slow-<0.024 m·day −1 ; slow-0.024 ÷ 0.12 m·day −1 ; moderately slow-0.12 ÷ 0.48 m·day −1 ; moderate-0.48 ÷ 1.56 m·day −1 ; moderately rapid-1.56 ÷ 4.20 m·day −1 ; rapid-4.20 ÷ 5.81 m·day −1 ; and very rapid->5.81 m·day −1 .
Undisturbed soil samples were taken from each genetic horizon, using cylinders with a volume of 100 cm 3 (three replicates). In addition, approximately 1.0 kg of disturbed soil from each genetic horizon was used to determine the soil texture, and other laboratory analyses were performed.

Laboratory Analysis
The properties of the collected soil samples were determined as follows.


The soil texture was determined using the Bouyoucose-Casagrande areometric method, based on a measurement of the density of the soil suspension during progressive sedimentation, and a sieving method to fractionate the sand. The contents of the particle size classes (sand, 2.0-0.05 mm; silt, 0.05-0.002 mm; and clay, <0.002 mm) were determined according to the Soil Taxonomy system from the United States Department of Agriculture [43].  The soil bulk density (BD) was determined based on the gravimetric method, based on cylinders (100 cm 3 ) for determining the mass of the dry soil per volume. The weight of each soil core was determined after drying in an oven at 105 °C for approximately 18-24 h. The dry bulk density for each core sample was then calculated as follows [44]: In the above, Ms is the mass of the dry soil weight (g), and Vt is the volume of the total soil sample (cm 3 ).


The soil water retention was investigated based on determining the soil suction using ceramic plates in a 5/15 bar pressure plate extractor. The pressure plate equipment used in this study was manufactured by the American Soil Moisture Equipment Corporation. In engineering practice, the soil suction is usually calculated in units of pF, as follows: Here, h is soil suction pressure (in cm H2O).
In the laboratory, the soil water potentials were measured at 300 cm H2O ≈ 33 kPa ≈ 2.5 pF (representing the field capacity) and 15,000 cm H2O ≈ 1500 kPa ≈ 4.2 pF (representing the permanent wilting point) [45]. The soil volumetric water content (θV) at pF = 2.5 and 4.2 was determined based on the gravimetric method. In particular, it was calculated as the ratio of the amount of water in the soil sample to the dry weight of the soil, after drying in an oven at 105 °C for approximately 18-24 h, as follows [44]:

Laboratory Analysis
The properties of the collected soil samples were determined as follows.

•
The soil texture was determined using the Bouyoucose-Casagrande areometric method, based on a measurement of the density of the soil suspension during progressive sedimentation, and a sieving method to fractionate the sand. The contents of the particle size classes (sand, 2.0-0.05 mm; silt, 0.05-0.002 mm; and clay, <0.002 mm) were determined according to the Soil Taxonomy system from the United States Department of Agriculture [43].

•
The soil bulk density (BD) was determined based on the gravimetric method, based on cylinders (100 cm 3 ) for determining the mass of the dry soil per volume. The weight of each soil core was determined after drying in an oven at 105 • C for approximately 18-24 h. The dry bulk density for each core sample was then calculated as follows [44]: In the above, M s is the mass of the dry soil weight (g), and V t is the volume of the total soil sample (cm 3 ).

•
The soil water retention was investigated based on determining the soil suction using ceramic plates in a 5/15 bar pressure plate extractor. The pressure plate equipment used in this study was manufactured by the American Soil Moisture Equipment Corporation. In engineering practice, the soil suction is usually calculated in units of pF, as follows: Here, h is soil suction pressure (in cm H 2 O).
In the laboratory, the soil water potentials were measured at 300 cm H 2 O ≈ 33 kPa ≈ 2.5 pF (representing the field capacity) and 15,000 cm H 2 O ≈ 1500 kPa ≈ 4.2 pF (representing the permanent wilting point) [45]. The soil volumetric water content (θ V ) at pF = 2.5 and 4.2 was determined based on the gravimetric method. In particular, it was calculated as the ratio of the amount of water in the soil sample to the dry weight of the soil, after drying in an oven at 105 • C for approximately 18-24 h, as follows [44]: Land 2021, 10, 959 7 of 21 In the above, V H2O is the volume of water in the soil sample (cm 3 ), and V t is the volume of the total soil sample (cm 3 ).

•
The Ks was measured under laboratory conditions using a laboratory permeameter ( Figure 3) and using the Darcy's law [46] with constant head method (Equation (6)) and falling head method (Equation (7) on undisturbed soil samples for three replications. The constant method can be used with virtually any soil, apart from poorly permeable soils such as clay, whereas the falling-head method is used to measure low-permeability soils, such as f.i. clay or peat samples [47,48]. The soil samples were placed in the laboratory permeameter and then saturated in water for 2-3 days. For the purpose of this study, a permeameter produced by the Eijkelkamp with a closed or open system and 25 holders was used. The Ks was calculated using the Darcy's [46] equation, as follows: Here, explanations below; Determining Ks using the constant head method was calculated from the following equation: Here, V is the volume of water flowing through the sample (cm 3 ); K is the permeability coefficient or 'Ks' (cm·day −1 or m·day −1 ); h is the water level difference inside and outside ringholder or sample cylinder (cm); L is the length of the soil sample (cm); i is the permeability rise gradient or h/L (-); A is the cross-sectional area of the sample (m 2 ); and t is the time used for flow through of water volume (day). During measuring the following parameters have been determined: L and A-constants, depending on the type of sample ring used; V-volume measured in the burette (1 mL = 1 cm 3 ); t-length of time lapse; h-calculated with the water levels measured with the water level meter.
Determining Ks using the falling head method was calculated from the following equation: Here, Ks is the permeability coefficient (cm·day −1 or m·day −1 ); h is the water level difference inside and outside ringholder or sample cylinder (cm) h 1 and h 2 water level difference inside and outside the ringholder at respectively t 1 (start) and t 2 (end); A is the surface of a cross-section of the sample (cm 2 ); L is the length of the soil sample (cm); t is the time between beginning and end of the measuring t 2 − t 1 (day); a is the cross-section surface of a ringholder or sample cylinder (cm 2 ) for a sample cylinder applies A = a; x is the evaporation factor (literature value): 0.0864 cm·day −1 or 0.000864 m·day −1 .
The permeability of the soil is also determined by viscosity of the soil solution. Viscosity depends on the temperature. The laboratory water temperature varies from 18 to 22 • C, whereas the average groundwater temperature is 10 • C. Therefore, for certain applications the permeability will have to be corrected for the viscosity of the soil solution (usually water). The Ks for viscosity was corrected using the following equation [49]: Here, K 10 is the corrected Ks at 10 • C (cm·day −1 or m·day −1 ); K T is the Ks at the applied temperature (cm·day −1 or m·day −1 ); h 10 is the dynamic viscosity of water at 10 • C (Pa·s); h T is the dynamic viscosity of water at T • C (Pa·s).
The soil classification was established according to the Polish Soil Classification [50], World Reference Base for Soil Resources (FAO and International Union of Soil Sciences (IUSS) Working Group World Reference Base (WRB), [51]), and USDA soil taxonomy (Soil Survey Staff, [43]).

Rosetta Description
The Ks values were calculated using the published neural network program Rosetta, with hierarchical PTFs based on five levels of input data [5]. The first level (PTF-1) used soil textural classes based on a lookup table providing parameter means for each USDA textural class. The second level (PTF-2) used the sand, silt, and clay percentages as inputs, and in contrast to PTF-1, provided hydraulic parameters that varied continuously with the texture. The third level (PTF-3) included the predictors used in the level PTF-2, along with the soil dry-bulk density. The fourth level (PTF-4) used PTF-3 and the soil volumetric water content (θ) at a water suction of 33 kPa (2.5 pF). The last level (PTF-5) comprised all of the other parameters, plus the θ value at a water suction of 1500 kPa (4.2 pF) [38,52]. The necessary input data for the Rosetta program (e.g., % soil texture group and bulk density) were determined based on the laboratory analyses of soil samples from the different soil horizons (Tables 2-4).
In the Rosetta program, the relationship between θ and the water suction (h), i.e., the water retention θ(h), as well as Ks, are described using the well-known Mualem-van Genuchten equations [53] and are given as follows: In the above, θ(h) is the soil volumetric water content (cm 3 •cm −3 ) at suction h (cm H2O); θs and θr are the saturated and residual water contents (cm 3 •cm −3 ) at h = 0 and 15,000 cm H2O, respectively; α (>0 in cm −1 ) is related to the inverse of the air entry suction; n (>1) is a measure of the pore-size distribution, and m = 1 − 1/n; and Ks is the saturated hydraulic conductivity (m•day −1 ), as mentioned above.

Rosetta Description
The Ks values were calculated using the published neural network program Rosetta, with hierarchical PTFs based on five levels of input data [5]. The first level (PTF-1) used soil textural classes based on a lookup table providing parameter means for each USDA textural class. The second level (PTF-2) used the sand, silt, and clay percentages as inputs, and in contrast to PTF-1, provided hydraulic parameters that varied continuously with the texture. The third level (PTF-3) included the predictors used in the level PTF-2, along with the soil dry-bulk density. The fourth level (PTF-4) used PTF-3 and the soil volumetric water content (θ) at a water suction of 33 kPa (2.5 pF). The last level (PTF-5) comprised all of the other parameters, plus the θ value at a water suction of 1500 kPa (4.2 pF) [38,52]. The necessary input data for the Rosetta program (e.g., % soil texture group and bulk density) were determined based on the laboratory analyses of soil samples from the different soil horizons (Tables 2-4).
In the Rosetta program, the relationship between θ and the water suction (h), i.e., the water retention θ(h), as well as Ks, are described using the well-known Mualem-van Genuchten equations [53] and are given as follows: In the above, θ(h) is the soil volumetric water content (cm 3 ·cm −3 ) at suction h (cm H 2 O); θs and θr are the saturated and residual water contents (cm 3 ·cm −3 ) at h = 0 and 15,000 cm H 2 O, respectively; α (>0 in cm −1 ) is related to the inverse of the air entry suction; n (>1) is a measure of the pore-size distribution, and m = 1 − 1/n; and Ks is the saturated hydraulic conductivity (m·day −1 ), as mentioned above.

Statistical Analysis and Model Performance Evaluation
The data set comprised the analytical results from the soil samples collected from the 16 soil pits. For the statistical analysis, the procedures provided by the program Statistica PL Version 12.5 were used, with a 5% significance level. The minimum and maximum values were determined for each physical parameter of the soil, and the arithmetic mean, median, standard deviation (SD), and coefficient of variation (CV) were computed. Moreover, a Spearman correlation test was conducted for the dataset. The correlation strength and direction of the relationship between two variables were determined based on Spearman's rank correlation coefficient (rho), which takes values from −1 to 1. A positive sign of the coefficient indicated the existence of a positive correlation, and vice versa. The closer the values were to −1 and 1, the stronger the correlation. When rho = 0, there was no correlation between the examined variables [54]. This statistical method was chosen after finding that most data were not normally distributed (Shapiro-Wilk test). All of data are presented in the tables and graphs, i.e., to produce a visual image that is helpful in interpreting the results.
Based on the available soil data in this study, five widely used PTFs were selected to assess their respective performances in estimating Ks, via comparison with the measured Ks values (via the DRI) for the soil cores collected in the field (using a laboratory permeameter). The calculated values of Ks obtained using Rosetta were compared to the corresponding measured values and evaluated using two statistical parameters: rootmean-square deviation (RMSD) and Nash-Sutcliffe efficiency (NSE). The RMSD gives the mean difference between the measured and calculated values of Ks, and is calculated as follows [38,55]: where: m i is the measured value of Ks, and p i is the corresponding (predicted) value of Ks as obtained using the Rosetta program.
The RMSD values are always non-negative and should be as low as possible; RMSD = 0 indicates a perfect fit of the model to the measurement data.
The NSE compares measured and predicted values (here, for Ks), and is given as follows [56]: In the above, n is the number of observations, m i and p i are the measured and predicted values of Ks, respectively, and m i is the mean of the measured Ks values.
The NSE values range from −∞ to 1. A value of NSE = 1 corresponded to a perfect match of calculated values to the measured values of Ks. The 'efficiency NSE' = 0 indicated that the calculations of Ks using the Rosetta program were as accurate as the mean of the measured data; in contrast, an efficiency NSE < 0 occurred when the measured mean was a better predictor than the model or, in other words, when the residual variance, as described by the numerator in the equation above, was larger than the data variance described by the denominator.

Results and Discussion
According to the PTG [50], FAO and IUSS Working Group WRB [52], and USDA soil taxonomy [43], the examined soils were classified as follows.  As shown in Table 2 and Figure 4, the examined soils are heterogeneous in terms of their textures. Seven textural classes are observed in the examined soil profiles: silt loam (SiL; n = 42), sandy loam (SL; n = 10), loam (L; n = 6), clay (C; n = 3), sandy clay loam (SCL; n = 3), clay loam (CL; n = 2), and silty clay loam (SiCL; n = 2). The silt loam is the predominant soil texture. The SD values (Tables 3 and 4) for most soil properties are large, indicating moderate to strong variability; this is not uncommon for soils [12,57,58]. Typically, the CV values are < 45%, indicating the mean variability of the measurement data. The highest CV values are observed for sand (57.0%) and clay (57.2%). The sand content in the soil profiles is between 11% and 76%, the silt content is between 16% and 72%, and the clay content is between 8% and 65% (Tables 2 and 4). The high swelling-clay content could be a cause of the poor water flow through the samples, as well as a consequence of abnormal Ks results [59]. As shown in Table 2 and Figure 4, the examined soils are heterogeneous in terms of their textures. Seven textural classes are observed in the examined soil profiles: silt loam (SiL; n = 42), sandy loam (SL; n = 10), loam (L; n = 6), clay (C; n = 3), sandy clay loam (SCL; n = 3), clay loam (CL; n = 2), and silty clay loam (SiCL; n = 2). The silt loam is the predominant soil texture. The SD values (Tables 3 and 4) for most soil properties are large, indicating moderate to strong variability; this is not uncommon for soils [12,57,58]. Typically, the CV values are < 45%, indicating the mean variability of the measurement data. The highest CV values are observed for sand (57.0%) and clay (57.2%). The sand content in the soil profiles is between 11% and 76%, the silt content is between 16% and 72%, and the clay content is between 8% and 65% (Tables 2 and 4). The high swelling-clay content could be a cause of the poor water flow through the samples, as well as a consequence of abnormal Ks results [59]. With regard to the surface soil horizons (n = 32), the mean soil bulk density (ρb) is 1.55 g•cm −3 , whereas that for overall depths is 1.60 g•cm −3 (Tables 3 and 4). The soil volumetric water contents at the matric potentials of 33 kPa (2.5 pF) and 1500 kPa (4.2 pF) in both cases have similar mean values (0.350/0.351 cm 3 •cm −3 and 0.152/0.155 cm 3 •cm −3 , respectively).
Higher values of Ks are found in the upper horizons of the tested soils, owing to the significant macroporosity caused by plant roots, reconsolidation [60], and annual loosening from agrotechnical works and alternating soil freezing and thawing [61]. In addi- With regard to the surface soil horizons (n = 32), the mean soil bulk density (ρ b ) is 1.55 g·cm −3 , whereas that for overall depths is 1.60 g·cm −3 (Tables 3 and 4). The soil volumetric water contents at the matric potentials of 33 kPa (2.5 pF) and 1500 kPa (4.2 pF) in both cases have similar mean values (0.350/0.351 cm 3 ·cm −3 and 0.152/0.155 cm 3 ·cm −3 , respectively).
Higher values of Ks are found in the upper horizons of the tested soils, owing to the significant macroporosity caused by plant roots, reconsolidation [60], and annual loosening from agrotechnical works and alternating soil freezing and thawing [61]. In addition, zoogenic channels (especially earthworm channels) and plant residues increase soil permeability [62]. For these reasons, the samples taken from arable and sub-arable horizons are less dense than those from lower horizons; the latter are also subject to self-compacting (Tables 2-4). Table 5 presents descriptive statistics concerning the values of Ks for the topsoils and subsoils (n = 32), as obtained from in situ and laboratory measurements, and as estimated with the use of the PTF models in the Rosetta program (variant 1). In contrast, Table 6 presents statistical measures considering the Ks values of all of the soil horizons (n = 68), but excludes the in situ tests not conducted on larger and deeper horizons, for practical reasons (variant 2).  In variant 1, the characteristic values of Ks (maximum, mean, and median values) as obtained from the in situ measurements are approximately four times higher than the Ks values obtained from the laboratory measurements. Such large discrepancies between the results obtained from both direct methods result from the manner of conducting the measurements; in particular, small-volume soil samples, which do not allow for the inclusion of all of the factors influencing the soil permeability, were taken for the laboratory tests. The mean Ks value obtained from the in situ measurements, i.e., 0.85 m·day −1 , is four times higher than those estimated by PTF-1 and PTF-2, almost half lower than that obtained from PTF-3, and twice as high as the values obtained from PTF-4 and PTF-5. In addition, the median value from the in situ tests (0.40 m·day −1 ) is from 1.2 to 2.2 times higher than the values obtained from the four PTF models-the exception is model PTF-3, whose median value is four times higher than that from the in situ tests. The maximum and mean values of the Ks values obtained from laboratory tests, compared to those from the in situ method, are usually more similar to those obtained from the PTF models (Table 5). The mean Ks values obtained from the in situ measurements and PTF-3 indicate the mean class, whereas the results obtained from laboratory measurements and the other four PTF models indicate a moderately slow basic infiltration class [42].
Most of the subsurface horizons are characterised by low Ks values, usually not exceeding 0.50 m·day −1 . In variant 2, covering the analysis of all soil horizons, the mean value of Ks from the laboratory measurements (0.21 m·day −1 ) is very close to the results obtained from the PTF-1 and PTF-2 models, almost six times lower than the Ks value from PTF-3, and three times lower than the values estimated by PTF-4 and PTF-5 ( Table 6). The mean Ks values obtained from the laboratory measurements and PTF-1, PTF-2, and PTF-3 indicate the moderate class, whereas the results obtained from the other three PTF models indicate a moderately slow basic infiltration class [42].
Regardless of the amount of data taken for analysis (n = 32 or 68), the values of Ks are usually characterised by high or very high random variability (CV = 26.1-230.7%), not only in comparison to the other soil hydrophysical properties, but also compared to the methods used to determine Ks, i.e., field, laboratory, and PTF tests (Tables 5 and 6, Figure 5). This regularity has been described in various other studies, for example, Rezaei et al. [63]. As reported by Merdun et al. [64] and García-Gutiérrez et al. [65], the high spatial variability between the soil horizons with regard to Ks are the result of the soil heterogeneity, which hampers the prediction of Ks using mathematical models.
value of Ks from the laboratory measurements (0.21 m·day −1 ) is very close to the results obtained from the PTF-1 and PTF-2 models, almost six times lower than the Ks value from PTF-3, and three times lower than the values estimated by PTF-4 and PTF-5 ( Table  6). The mean Ks values obtained from the laboratory measurements and PTF-1, PTF-2, and PTF-3 indicate the moderate class, whereas the results obtained from the other three PTF models indicate a moderately slow basic infiltration class [42].
Regardless of the amount of data taken for analysis (n = 32 or 68), the values of Ks are usually characterised by high or very high random variability (CV = 26.1-230.7%), not only in comparison to the other soil hydrophysical properties, but also compared to the methods used to determine Ks, i.e., field, laboratory, and PTF tests (Tables 5 and 6, Figure  5). This regularity has been described in various other studies, for example, Rezaei et al. [63]. As reported by Merdun et al. [64] and García-Gutiérrez et al. [65], the high spatial variability between the soil horizons with regard to Ks are the result of the soil heterogeneity, which hampers the prediction of Ks using mathematical models. The Ks values obtained from the laboratory measurements are statistically significant (p < 0.05), negatively correlated with the clay fraction and density, and positively correlated with the sand fraction (Table 7). Similar observations were noted by Zhao et al. [13], who conducted studies on the Loess Plateau of China. The interdependence of Ks values and soil textures were well documented by Hillel [66], García-Gutiérrez et al. [65], and other researchers [3]. In the case of in situ measurements, the soil hydraulic conductivity is only significantly positively correlated with the sand fraction (rho = 0.29). The PTF-1 model is significantly positively correlated with the soil density, but negatively correlated with the clay fraction and soil volumetric water content (θ) at water suctions of 33 and 1500 kPa (2.5 pF and 4.2 pF). The PTF-2 and PTF-3 models show a significantly positive correlation with the sand fraction and negative correlations with the clay fraction and soil volumetric water content (θ) at 2.5 and 4.2 pF. The PTF-4 and PTF-5 models are significantly positively correlated with the sand fraction, and are negatively correlated with the dust fraction (Table 7). The Ks values obtained from the laboratory measurements are statistically significant (p < 0.05), negatively correlated with the clay fraction and density, and positively correlated with the sand fraction (Table 7). Similar observations were noted by Zhao et al. [13], who conducted studies on the Loess Plateau of China. The interdependence of Ks values and soil textures were well documented by Hillel [66], García-Gutiérrez et al. [65], and other researchers [3]. In the case of in situ measurements, the soil hydraulic conductivity is only significantly positively correlated with the sand fraction (rho = 0.29). The PTF-1 model is significantly positively correlated with the soil density, but negatively correlated with the clay fraction and soil volumetric water content (θ) at water suctions of 33 and 1500 kPa (2.5 pF and 4.2 pF). The PTF-2 and PTF-3 models show a significantly positive correlation with the sand fraction and negative correlations with the clay fraction and soil volumetric water content (θ) at 2.5 and 4.2 pF. The PTF-4 and PTF-5 models are significantly positively correlated with the sand fraction, and are negatively correlated with the dust fraction (Table 7). Red color means, that the determined correlation coefficient is significant (for the significance level α = 0.05).
No significant correlations were detected between the Ks values obtained from the in situ measurements and the PTF models. However, with the exception of the PTF-1 model, the Ks values from the laboratory tests as estimated by the four PTF models are positively correlated (Table 7). For these reasons, only the linear correlations between the values obtained from the laboratory tests and PTF models are presented in Figure 6a-e.
In general, the matching of the Ks results obtained from the direct in situ measurements with the results estimated by the Rosetta PTFs is poorer than that obtained from the laboratory measurements and estimated by the Rosetta PTFs. This is confirmed by the RMSD values, which range from 1.40 to 1.64 in the first case, whereas in the second case, they range from 0.19 to 1.58. In addition, the values of the coefficient of determination (R 2 ), which indicates the extent to which the model explains the gathered measurement data, are usually much higher when comparing the data obtained from the laboratory tests and that from the PTF models (Table 8). The RMSD values of the Ks results obtained from the laboratory measurements and estimated with PTFs reveal the best match for the case of the PTF-2 model. The NSE and R 2 values are the best in this case (Table 8), indicating that the granulometric composition (% of sand, silt, and clay) in the PTF-2 model has the greatest influence on the water permeability of the soils.

Conclusions
The search for alternatives to labour-intensive and cost-intensive in situ tests and laboratory methods for determining the value of Ks, at a time of global climate change, seems to be a priority. The increased frequency of periods with excesses and scarcities of water will require the greater use of soil irrigation and drainage facilities, and the modelling of surface runoff and erosion phenomena in land-formed areas. Correspondingly, knowledge regarding soil permeability is required for all technical and non-technical activities.
In general, when using a descriptive comparative analysis, it is difficult to show a significant relationship between Ks values from in situ tests and the diverse hydrophysical properties of arable soils, which are subject to cyclic loosening, reconsolidation, and/or freezing and thawing. The high sensitivity of this physical parameter to biotic and abiotic factors, especially in the soil upper horizons, hampered a satisfactory adjustment to the Ks values obtained from the PTF models. Therefore, in situ tests, although difficult

Conclusions
The search for alternatives to labour-intensive and cost-intensive in situ tests and laboratory methods for determining the value of Ks, at a time of global climate change, seems to be a priority. The increased frequency of periods with excesses and scarcities of water will require the greater use of soil irrigation and drainage facilities, and the modelling of surface runoff and erosion phenomena in land-formed areas. Correspondingly, knowledge regarding soil permeability is required for all technical and non-technical activities.
In general, when using a descriptive comparative analysis, it is difficult to show a significant relationship between Ks values from in situ tests and the diverse hydrophysical properties of arable soils, which are subject to cyclic loosening, reconsolidation, and/or freezing and thawing. The high sensitivity of this physical parameter to biotic and abiotic factors, especially in the soil upper horizons, hampered a satisfactory adjustment to the Ks values obtained from the PTF models. Therefore, in situ tests, although difficult and timeconsuming, should continue to be conducted to obtain a reliable measurement database that is free of errors resulting from the soil type variability, and that considers the influences of fauna and flora on the physical properties of soil. This is of particular importance for the modelling of surface runoff and erosion processes and design of surface irrigation systems (e.g., sprinklers), where only 40-50 cm of the upper horizon of the soil is wetted, and the technical parameters of the corresponding drainage facilities are highly dependent on soil permeability.
The in situ measurements of Ks are significantly correlated only with the laboratory tests; nevertheless, the values obtained in situ are much higher, and are influenced by a number of factors, including those of a methodological nature (e.g., the larger diameter of the measuring cylinder). In contrast, the laboratory measurements show a significant correlation with the Ks values estimated using the models PTF-2 to PTF-5. The best match is found (based on the RMSD, NSE, and R 2 values) when using the PTF-2 model in the Rosetta program. This means that the most reliable Ks values can be obtained from the percentages of the sand, dust, and clay fractions. Therefore, there is a rational alternative, in the form of the PTF, for the costly and labour-intensive determinations of Ks conducted in the laboratory. The use of the PTF-2 model allows for the determination of the hydraulic conductivity of deeper soil horizons, where the variability of the hydrophysical properties is lower and there are large methodological limitations in conducting field work or difficulties in taking intact structural samples for laboratory tests. This type of solution can be used to collect data for designing subsurface drainage and irrigation facilities, some of the technical parameters of which are determined by the Ks values.
In conclusion, the research results confirmed the remarkable difference between the performance of PTFs in the Rosetta and the Ks results obtained from field or laboratory measurements. The study thus provided evidence that the Ks values obtained from local dataset give better results performances in predicting the soil saturated hydraulic conductivity, compared to PTFs in the Rosetta derived from large datasets, which in turn confirm the limitation of applying PTFs developed from one region to other regions.