Next Article in Journal
Managing the Santa Cruz River: Prospects for a Future Binational Minute
Previous Article in Journal
Groundwater Contamination Source Recognition Based on a Two-Stage Inversion Framework with a Deep Learning Surrogate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Statistical Approach for the Assessment of Saturated Hydraulic Conductivity Values of Unsaturated Urban Soils Obtained by Field Infiltration Tests

Groundwater Engineering Research Centre, Technical University of Civil Engineering, bd. Lacul Tei, 124, 020396 Bucharest, Romania
*
Author to whom correspondence should be addressed.
Water 2024, 16(13), 1908; https://doi.org/10.3390/w16131908
Submission received: 9 April 2024 / Revised: 15 May 2024 / Accepted: 1 July 2024 / Published: 3 July 2024
(This article belongs to the Topic Urban Hydrogeology Research)

Abstract

:
An evaluation and interpretation of the obtained results focusing the hydraulic conductivity of anthropogenic saturated soil, k, has been performed on an urban area vadose zone. Four methods have been used to quantify the hydraulic conductivity: the tube infiltrometer (TI), the double ring infiltrometer (DRI), the minidisk infiltrometer (MDI) and the inversed auger (IA). This study comprises (a) a comparative analysis of the results obtained by each method between several trials performed at the same location and at distinct locations within the plot, (b) a comparative analysis of the results of all methods, and (c) a statistical analysis regarding the correlation between k as a dependent variable and the infiltration area A as the main independent variable. To select the k values close or corresponding to the saturation state for TI and IA methods, a domain of validity was defined. A new parameter, k* = k/A, was introduced which represents the hydraulic conductivity corresponding to an infiltration surface unit (1 cm2). An increase in this ratio with the increase in the infiltration area, within the same method or between different methods, indicates the heterogeneity of the terrain but especially the fact that the infiltration area no longer represents the main independent variable on which the hydraulic conductivity depends for the saturated state.

1. Introduction

Studies of city-scale urban aquifer water dynamics should be based on accurate urban groundwater balance analysis including natural and human-induced water sources, geological and anthropogenic strata, and the entire set of the urban infrastructure elements. Since the knowledge of deposits located in the shallow urban subsurface is increasingly important for urban planning [1], there is a need for better classification of anthropogenic materials and their hydraulic conductivity. To date, there are very few published hydraulic studies on anthropogenic strata in cities, and most of them are of a pioneering nature. Previous studies developing city-scale urban hydrogeological models [2,3,4,5] demonstrated the necessity of correctly quantifying the anthropogenic strata hydraulic conductivity to properly assess aquifer dynamics as well as groundwater recharge from precipitation.
Techniques to quantify accurately the saturated hydraulic conductivity k (LT−1) have yet to be developed [6]. Several studies have highlighted the spatial variability of the saturated hydraulic conductivity in soil characterized by uniform texture, even when applying a single test method. This is due to numerous and variable factors, not all of which are detectable (preferential flow, cracks, roots, structure, and others). When different techniques with different theoretical bases are considered, the topic becomes still more complex, even in homogeneous soil.
Comparative studies show the differences both between the saturated k values achieved when applying the same method while modifying the measurement device infiltration area for the same lithological unit of natural ground, and between k values obtained when applying distinct methods (using different infiltration areas corresponding to the testing device) on the same natural ground [7,8]. Some authors [9,10], applying the DRI method on natural undisturbed ground, obtained insignificant differences between the values of hydraulic conductivity when varying the infiltration area, but the values for increased infiltration area showed an appreciable decrease in the standard deviation.
Other studies [11] show differences of about one order of magnitude between the saturated k values obtained when applying the DRI and the tension disk infiltrometer (TDI) methods on clay. Also, large differences have been obtained when applying the TI and the TDI methods in situ, as well as when laboratory tests are performed [12]. The differences between the results obtained by applying TI and TDI arise mainly due to the water volume able to pass through the macroporous medium. This is much greater in the case of TI, where the water moves vertically, than for TDI, where there is also a horizontal component.
Other authors show that the differences in the estimation of the k value are due to TDI membrane’s resistance when it is crossed by air bubbles due to the water flow geometry [13,14], or even due to the spontaneous development of plant roots (grass) during the test run [15]. Very small differences between the k values obtained by TI and IA tests on clayey soils without microstratification were highlighted by van Hoorn [16]. Bagarello et al. [17] highlighted that by using the TI method on sandy silty soils, more accurate results could be obtained than on clayey soils. Most authors point out that some methods applied on so called “homogeneous” natural lands overestimate and others underestimate the hydraulic conductivity value, the process of water infiltration into unsaturated zone being extremely complex. It depends on a multitude of factors such as the presence of macropores (due to the structure and texture of the terrain, plant roots, etc.), terrain heterogeneity, the requirements for preparation and completion of the test, the differences between the algorithms of the methods, the predominant lithological nature of the terrain, and last, but not least, the experience of the operator and of the analyst.
Through time, substantial changes occur in the composition, structure, and texture of the soil, especially in urban areas as a result of the anthropogenic activity. Urban soil is defined as non-agricultural material resulting from anthropogenic activities and emerges through the processes of filling, extraction, and contamination of the natural surface with a minimum thickness of 50 cm [18] or as material that was modified, altered, and transported by human activities in the urban environment [19].
Urban soils are highly heterogeneous and may consist of mixtures of demolition materials, household waste, slag, allochthonous material, and others. These can be captured in an allochthonous or in a native matrix. In relation to the “homogeneity” or “quasi-homogeneity” of the natural terrain, two main groups of urban soils can be distinguished: (a) urban soil made of lithological homogeneous material, for example, allochthonous clayey material coming from a well-defined lithological sequence and used to fill a negative morphological surface; and (b) urban soil consisting of lithological mixtures (e.g., sand, gravel, clay) with materials originating from demolitions and other sources.
In the first group, the urban soil can be considered as a “homogeneous” layer, but different from the native terrain. Even if, from the lithological point of view, its characteristics could be similar, there are differences of porosity, compaction, mechanical properties, and other attributes. These differences can occur both vertically and horizontally, depending on the anthropogenic activity the soil has been exposed to. In this situation, the assessment of the saturated hydraulic conductivity k of the vadose zone can be conducted by methods also used for natural soils; however, the use of the data must take into consideration the allochthonous nature of this strata as well as the lithological nature of the indigenous substrata which may have distinct hydraulic properties. Thus, field tests performed on natural and compacted loess showed that the hydraulic conductivity was reduced by 1–2 orders of magnitude, leading to a low rate of water infiltration [20]. At a small scale, laboratory studies have identified useful methods for simulating water hydrodynamics in the vadose zone for layered soils [21].
In the second group, the lack of homogeneity of urban soil is obvious and no similarity is expected between the results obtained by applying several methods of determining saturated hydraulic conductivity of the vadose zone. The determination of the hydraulic conductivity using usual methods can lead to very dissimilar values, even when using the same method. This is due to the fact that each method provides point-based information. Consequently, interpretation and validation of any results are critical.
Especially for “homogeneous” soils, by using test devices with large infiltration area, values of hydraulic conductivity are obtained for the saturated state much closer to a supposed real value. From a theoretical point of view, in the case of group (b) soils, the best way to determine the hydraulic conductivity for the saturated state would be the flooding of the entire surface, which is completely impractical. The only possibility is to perform tests on distinct points. However, the urban soil being heterogeneous over relatively short distances, the values obtained may differ considerably for the same test method but also between different methods. In this situation, the following questions can be asked: (1) How can the achievement of the saturation or incipient saturation state be judged, for the tests which allow multiple infiltration runs in the same location as TI and IA methods and how can the selection of the results be made corresponding to this state? (2) If, for a relatively small representative elementary volume (REV) comprising two very close test locations where methods with distinct infiltration surfaces are applied, the amplitude of the selected values is smaller for the device with a smaller infiltration area, does this mean that these values better characterize the tested environment? (3) If the hydraulic conductivity for the saturated state is an intrinsic value that characterizes the type of urban soil (heterogeneous) in the test point and it is dependent on the infiltration area, then what is the influence of other factors on the results? (4) Is it possible to obtain a characteristic value of the saturated hydraulic conductivity for a heterogeneous anthropogenic terrain?
In order to obtain the saturation state for homogeneous soils, either the soil in the location of the test device has to be saturated for a time, which may depend on the initial saturation of the material, or more trials have to be performed. When the values for consecutive attempts of estimation of k are sufficiently similar, it can be considered that the saturated state has been reached. In fact, using an infiltration technique in an initially unsaturated soil under ponding conditions, the field saturated soil hydraulic conductivity is obtained due to entrapped air bubbles. Císlerová et al. [22] highlighted that, at higher moisture content, the air entrapment in large pores sealed off by water films will increase drastically and the saturated hydraulic conductivity will accordingly decrease. According with Sakaguchi et al. [23], the saturated hydraulic conductivity measured on a soil that contains entrapped air can be smaller than the unsaturated hydraulic conductivity close to saturation. On the other hand, repeating the same experiment many times at the same point can induce weakening of the particle bonds and migration of small particles [24], leading to the change of hydraulic conductivity. Therefore, air entrapment in the soil is a complex phenomenon that can have effects on “saturated soil hydraulic conductivity” obtained by field test methods, this term being often used for practical purposes instead of field saturated hydraulic conductivity, which is more or less close to the real state of the soil saturation. As will be seen in the following sections, in the case of urban soil, there are often large differences between the values of hydraulic conductivity obtained after consecutive trials in the same location for TI and IA methods. In this sense, a domain of validity of the saturated hydraulic conductivity values has been defined, within which are included the values selected for the calculation of the average value of the hydraulic conductivity when these are close to or reaching the saturation state.
In a homogeneous and isotropic ideal soil, the value of the hydraulic conductivity for the saturated state k, determined by any method, would be the same in any location and the amplitude between observed values would be Apl = 0. If two or more distinct test methods are applied, using the same or different infiltration areas A, theoretically, k1 ≠ k2 ≠… ≠ ki should be obtained (i represents the test method number) and Apl > 0. Differences could appear due only to the specific calculation algorithm and/or operator error, and the infiltration area has no influence. In reality, the soil is not homogeneous, the homogeneity representing an idealization and a simplification of the analyzed REV. And so, as we highlighted before, the infiltration area represents an independent variable which can influence k values because this surface can comprise other influencing factors (soil heterogeneity, cracks, roots, etc.) that may or may not be dominant. If the same test method uses distinct infiltration areas, then the differences between the k values, obtained as a dependent variable, will also include the influence of the area variation as an independent variable and the amplitude Apl > 0 between saturated k values. Also, between the distinct applied methods, the amplitude will be Apl > 0 and will contain the influences induced by the operator, calculation algorithm, and the size of the infiltration area.
From the statistical point of view, selecting values with the lowest amplitudes better characterizes the study environment. If the infiltration area changes for different locations when applying the same or distinct methods, the occurrence of extreme values with relatively low frequency can have a disproportionate effect on the amplitude of the selection and consequently lead to a misinterpretation of the value of the dependent variable [25], and comparing only the k values, this aspect cannot be highlighted. To minimize this effect, in order to compare results obtained by different test methods, we have defined a new variable, k* = k/A, which represents the saturated hydraulic conductivity corresponding to a unit infiltration surface (equal to 1 cm2). We assumed this unitary surface to be homogeneous. Therefore, with the increase in the infiltration area, the value of the k/A ratio must decrease (hypothesis “0”). The values of the amplitudes determined for the variable k* better describes the results obtained by the test method, and the infiltration area is the main variable which controls the k values. In reality, there are other variables, especially at the surface of the soil, which can influence the k values: cracks, roots, and others. In urban soils, in addition to the heterogeneity, undetected underground cavities may also be present. If the k/A ratio between the methods increases with an increase in the infiltration area (hypothesis ”1”), then other variables have a decisive influence on the k value (either as single or cumulative variables). A statistical analysis has been also made for the results obtained from different methods to test whether the infiltration area of the frequently used test devices can be considered the main independent variable, which controls the value of saturated hydraulic conductivity. This should help to determine the cumulative weight of other independent variables, having an influence on the dependent variables k and k*, and establish the significance of the k/A ratio.
Answers to questions (1) to (4) are provided by this study through (a) a comparative analysis of the results obtained by each method between several trials performed at the same location and at distinct locations within the same plot and defining the domain of validity and selection of saturated k values in the field for TI and IA methods, (b) a comparative analysis of the results obtained by used methods, and (c) a stochastic analysis regarding the correlation between the anthropogenic soil hydraulic conductivity in the saturated state k as a dependent variable and the infiltration area A as the main independent variable using a qualitative interpretation based on k/A ratio significance.
In urban areas with a wide spread of anthropogenic soil, a correct determination of the hydraulic conductivity corresponding to the saturated state of this type of material has multiple applications, such as optimization of the urban drainage systems, analysis of the hydraulic interaction between urban infrastructure elements and green spaces, optimal installation of green infrastructure, and many others. The correct quantification of the water recharge from liquid or solid rainfall (snow), required for urban hydrological and hydrogeological studies, is still a challenge.

2. Materials and Methods

Four methods were used to assess the saturated hydraulic conductivity of the vadose zone of urban soil: the tube infiltrometer method (TI), the double ring infiltrometer method (DRI), the minidisk infiltration method (MDI), and the inversed auger method (IA). The water infiltration through the soil, for the first three methods, takes place only through the area corresponding to the base of the test device. For the IA method, the infiltration occurs both through the area corresponding to the base of the drilling hole as well as through its lateral surface. A comparative analysis of the results was undertaken for the all four methods. The first three (TI, DRI, and MDI) have been separately compared with the IA method, which tested the soil at depths of 30, 60, and 90 cm in order to highlight the possible influence of the surface factors. The details of the devices are shown in Figure 1. The tests were carried out on the same site (urban unsaturated soil) in similar weather conditions during a drought period.

2.1. Tube Infiltrometer/Single Ring Infiltrometer Method (TI)

The device consists of a metal ring or a cylinder with a wall thickness of less than 5 mm, sharpened at the bottom to penetrate the unsaturated ground, at a depth d (Figure 1a). When water is poured into the tube, a hydraulic load H is created and maintained constant throughout the duration of the test by means of a reservoir, in which the volume of water is equal to the volume of water poured into the infiltrometer. Then, the water level in the reservoir is measured for certain time steps and the cumulative infiltration I (L) and the infiltration rate I = dI/dt (LT−1) are determined and used to compute the saturated hydraulic conductivity k (LT−1) using the following relation [26,27]:
I = A w · t + B w · t
where t (T) is time, and Aw (LT−1) and Bw (LT−0.5) are coefficients related, respectively, to the hydraulic conductivity and sorptivity. Assuming that, in the last part of the infiltration process, a stationary regime is reached, then I=Aw∙t + c = k∙a’∙f∙t + c, where Aw represents the slope of the straight line and c the intersection with the y axis, k (LT−1) is the hydraulic conductivity, a’ = 0.9084 is a constant, and f is a correction factor depending on soil texture α* [28,29], hydraulic head H (L), and the ring geometry as in Figure 1a.

2.2. Double Ring Infiltrometer Method (DRI)

The infiltrometer is made up of two cylinders, open at both ends, concentrically placed on the surface of the ground: an inner and an outer ring (Figure 1b) which have to penetrate the terrain surface at a recommended depth of d = 0.05 … 0.10 m [30]. The test consists of achieving a constant H (L) water level in both cylinders and recording water volumes added at certain time steps to determine the cumulative infiltration I (L) and infiltration velocity, v = dI/dt (LT−1) where t (T) is the time. The saturated state hydraulic conductivity k is calculated using Philip’s equation [31,32]:
I t = S · t 0.5 + A · t
where S (LT−1/2) represents the sorptivity and A (LT−1) is a constant that represents the influence of the gravitational component at saturation, A = m∙k and m = 2/3 = ct. [33]. Using field data, two graphs are built: I = f(t) and v = f(t).

2.3. Minidisk Infiltrometer Method (MDI)

The infiltrometer is designed to measure the hydraulic conductivity of unsaturated ground at suction values ranging from −0.5 to −7 cm. It is composed of two chambers that are filled with water (Figure 1c). The bubble chamber controls the value of the suction. The water used for the infiltration test is stored in the “water reservoir”. At the bottom of the tank, there is a sintered steel disc with a thickness of 0.3 cm and a radius r = 2.25 cm, the total volume of the tank being V = 135 mL. Adjusting the suction (negative pressure) in the upper chamber prevents the water infiltration through macropores such as cracks, wormholes, and others. The infiltration of water through the soil is determined by the hydraulic forces that characterize the soil matrix. The hydraulic conductivity is dependent on pore geometry, soil water content, and capillary effect [34]. The value of hydraulic conductivity for the saturated state is obtained when all pores (including macropores) are full of water and the water flows. Even for close locations, quantification of the water flow through macropores can be challenging. Application of suction prevents the filling of macropores, and the results obtained by this method provide the hydraulic conductivity for the saturated state that characterizes the clay matrix. In this way, we can compare the values obtained for saturated k by applying the other methods (for REV corresponding to each location), which also include the effect of macropores or cracks, with the values obtained only for the clay matrix. In this sense, given the very small infiltration area of the device, the MDI locations chosen has no visible cracks or traces of vegetation, and at the start of tests, it was assumed that the presence of wormholes or underground roots was excluded.
For the calculation of the hydraulic conductivity coefficient, the pairs of values obtained during the test of cumulative infiltration I (L) and time t (T) are required, adjusted with the following relation [35]:
I = C 1 · t + C 2 · t
where C2 (LT−1/2) is a sorption-related parameter, C1 (LT−1) is related to the hydraulic conductivity parameter C1 = k∙A, and represents the gradient of the cumulative infiltration curve vs. the square root of time; A is a value that depends on van Genuchten parameters (n and α) corresponding to a particular texture of the soil, ray of infiltration disk, and the suction at the disk surface which were determined for 12 soil texture classes [36]. Starting from these values, the parameters are calculated for the minidisk infiltrometer, with a radius of 2.25 cm and suction between 0.5 and 7 cm, by the manufacturer, Decagon Devices, Inc., Plullman, WA, USA [37].

2.4. Inversed Auger Method (IA)

The method consists of drilling a hole at a given depth cased with a steel screen with the same diameter, filling it with water, and measuring the change of the water level at certain time intervals (Figure 1d). The saturated k value is given by the following equation [16]:
k = log h o + r 2 log h t + r 2 t = 1.15 · r · tan α
where t is time and other variables are as given in Figure 1d.

2.5. The Urban Area Experimental Site and Measurements Set Up

The test site area was an artificial green area located in Bucharest, Romania. This zone was initially a depression which, about 30 years ago, was filled with material about 1 m thick. Buildings and green areas have since been established. The lithological structure of the site is derived from an investigation borehole located at about 30 m NE of the test perimeter (Figure 2a). Anthropogenic material with a thickness of about 1 m, located near the ground surface, overlays a silty clay stratum. A slightly confined shallow sandy aquifer strata, with a depth of the piezometric level of 3.56 m, is located under the clayey layer.
This test site was selected for the following reasons:
The upper part of the ground is made up of urban soil;
The area is partially covered by vegetation (trees and grass) and could be considered a human-made “green area”;
The ratio between the perimeter length and width is about 5:1 and the space between different test locations is small, as is shown in Figure 2b, leading to a relatively balanced compensation of the assessed values’ variation due to the terrain heterogeneity, because smaller soil volumes are functionally more homogeneous than larger volumes [38];
The area is not affected by pedestrian traffic and therefore has not been subjected to compaction or other actions. During the tests, only the operator’s access, on established routes, was allowed.

2.5.1. Sampling and Laboratory Tests

In the vicinity of the DRI and TI test location (Figure 2b), samples of disturbed and undisturbed soil were collected down to depths of 20 cm. Also, disturbed samples, up to a depth of approximately 90 cm, were collected from the drill holes associated with the IA tests. Samples were collected after performing the infiltration tests to avoid terrain disturbance, with the exception of IA locations where samples were collected during the drilling works. Granulometric analyses were performed on all samples. On undisturbed samples, consistency indices, natural and dry density, porosity, and the degree of saturation were determined. Laboratory determinations revealed a filler material up to 0.90 m depth consisting of sand and gravel, brick debris, glass fragments, and roots, attached to a clayey matrix. The values of the indices and parameters determined in the laboratory, by groups of depth, are presented in Supplementary S1 as electronic Supplementary Material (ESM). The clayey matrix falls between the groups of medium- and high-plasticity (10 < Ip < 20), and in terms of consistency, it has a hard state (0.5 < Ic < 0.75). An exception is marked by the locations DRI 4/P01, that shows a plastic consistent state (0.5 < Ic < 0.75), and by DRI 5, that indicates a strong plastic state (0.75 < Ic < 1). The activity index IA of the clay material, which characterizes the intensity of the water interaction with fine particles (clay fraction), is determined as a ratio between the plasticity index Ip and the percentage of particles smaller than 2 μm. This indicates a terrain with medium activity (0.75 < IA < 1.25), values below 0.75 characterize a low active ground, and those over 1.25 a very active one.

2.5.2. Tube Infiltrometer Set Up

The tests were performed at six locations, as shown in Figure 2b, with the distances between locations varying between 1.1 and 2.5 m. In each location, several trials were carried out. Tubes with different heights and diameters were used. The maximum height of pouring water into the tubes H, the average of the hydraulic heads of the trials performed in each location Havg, the soil penetration depth of the infiltrometer d, and the internal radius of the tube r (see Figure 1a) are given on ESM as Supplementary S2. The different number of trials within each test location is attributed to the different periods of time in which water infiltration occurred, all locations being tested over the same time period. The value of the constant α* was chosen corresponding to the texture-structure category of the soil and is confirmed by the results of the granulometry tests.
At the same time, to reduce the influence of α* variation on the values obtained for k, we used tubes with a height over 20 cm, as the influence of the chosen α* on k decreases with the increase in the H value [39].
For the six tests, multiple trials for each location were performed for the following reasons: (a) the tests were carried out in all locations at the same time; (b) the relatively small diameter of the tubes and the heterogeneity of the soil led to the idea that it is quite unlikely to reach the saturated state of the soil (corresponding to the steady state) after a single attempt. This could happen, in the immediate vicinity of the surface, either due to high hydraulic conductivity, to the access paths being obstructing by the water-driven fine particles and subsequently transported by the water flow, or due to entrapped air; and (c) the physical effect of the soil–water interaction on time intervals can be better emphasized.
For each location, the following parameters were computed:
The hydraulic conductivity k of each trial Tn (T-trial, n-number of trials in each test location), according to the hydraulic head corresponding to the water height, on the basis of the accumulated S (cumulative drawdown) records of the hydraulic head Hn at certain time intervals (registered in the water reservoir);
The hydraulic conductivity kdcH for each TIi location using the average of the hydraulic heads (Havg) measured for n trials at each location i and, where the cumulative infiltration calculation was used, the average of the absolute drawdowns S of the hydraulic head Hn recorded at that location for the same cumulative time values;
The kavg hydraulic conductivity as the average of the values determined on n trials at each location.

2.5.3. Double Ring Infiltrometer Set Up

Five DRI tests were performed. Rings of 28 and 53.2 cm diameter were placed in the field at a depth of 5 cm, the area of the inner tube base being S = 615.75 cm2. The height of the water in both rings was H = 7 cm, corresponding to the position of the upper level marker. The lower marker was positioned 1.6 cm below, this height difference corresponding to a volume of V = 1000 cm3 of water, needed to be supplemented to maintain a constant level. The duration of the tests varied between 73 and 132 min. By maintaining the constant water level, a steady state flow regime was achieved when two successive additions of water showed no significant differences in volume.

2.5.4. Minidisk Infiltrometer Set Up

A total of 18 MDI tests were performed. The locations were chosen in the vicinity of the performed TI, DRI, and IA tests with a suction of −5 cm (drought conditions, relatively high content of sand and presence of macropores) and parameter A = 8.4.

2.5.5. Inversed Auger Method Set Up

The IA tests (depths between 0.30 and 0.83 m and with a diameter of 10 cm) were performed in six locations, away from the TI tests, on an alignment located at a distance of approximately 1.80 m from the TI test alignment. Because the ground is made up of fillings, stainless steel filters have been introduced in the holes to support them during the tests (avoiding the wall clogging). Several trials were made at each location. Their number, the average hydraulic head on each hoavg location, and the geometric data according to the notations of Figure 1d) are presented on ESM as Supplementary S4. The variable hydraulic head test has been chosen. The following parameters were determined for each location:
The hydraulic conductivity k of each test, according to the hydraulic head given by the initial height of the water column in the drilling ho at to = 0, and the hydraulic head values hi based on the cumulative drawdown Δh and corresponding to different times ti;
The hydraulic conductivity kdch for each location AIi using the average of the hydraulic heads measured for n tests at each location i; for the calculation of the dynamic level hi the mean of the absolute drawdown h of the hydraulic head Hn recorded during the n tests on the location for the same cumulative time values was used;
The hydraulic conductivity kavg as the average of the values obtained for n tests at each location.

3. Results

3.1. Tube Infiltrometer (TI)

The k values for each location and each trial are shown in Table 1. The results of the experimental records for each trial and test location, as well as the cumulative graphs of infiltration by time I = f(t) and infiltration rate by time I = f(t), are given on ESM as Supplementary S2. Saturated hydraulic conductivity was computed using Equation (1).
Theoretically, in a REV, with homogeneous soils, the hydraulic conductivity is minimum when the soil moisture is minimum and reaches a maximum value when the soil reaches the saturation state, the two parameters being directly correlated. In this experiment, the absolute values of the hydraulic conductivity obtained for the unsaturated and saturated state of the urban soil were not compared, but the relative values obtained for the hypothetical saturation or incipient saturation states for consecutive trials in the same location were compared. Moisture measurements could not be performed between successive attempts to avoid disturbing the ground. Thus, between two hypothetical saturation states corresponding to the i−1 and i trials, the infiltration velocity tends to become constant and the corresponding REV hydraulic gradient decreases by decreasing the pathway taken by the water particle as a result of the pore saturation. Consequently, the hydraulic conductivity, corresponding to the two hypothetical saturation states, decreases as the soil saturation increases until it becomes constant. To avoid any confusion, note that, in the following, the decrease in the hydraulic conductivity value with the increase in the soil saturation refers to the hydraulic conductivity values corresponding to the hypothetical saturation state between successive trials. It should be noted that the values of kdcH hydraulic conductivity (column 20—Table 1) determined using the average of the hydraulic head tests and the mean of the water level oscillations differ greatly from the arithmetic mean kavg (column 21—Table 1) and they are close to the measured hydraulic conductivity k values for the previous tests in each location. This is interpreted as meaning that the terrain was saturated due to water infiltration or that it was in the early stage of saturation. This is not valid for the TI6 location where, though there was a higher number of tests, there is no clear decrease in the hydraulic conductivity value as the number of tests increases (increasing the saturation state).
The variation of hydraulic conductivity could be due either to the ground heterogeneity or to a lower cohesion compared to the natural soil (through processes of entrainment of the soil particles by the infiltration water, which can clog different pathways and which are subsequently released). This phenomenon is more pronounced with the presence of anthropogenic elements such as brick, concrete, and glass fragments (identified in this location) which contribute to the hydraulic instability of the clayey matrix and of the permeable porous medium.
Hydraulic conductivity shows relatively low values for the TI2 location, which tend to stabilize due to the fact that the field is predominantly made up of fine material (clayey and silty, 82%) (Supplementary S1 on ESM—S02/P01). The questions, especially for heterogeneous areas containing anthropogenic products, are by which criterion can the correct value of the hydraulic conductivity be determined and how many trials/tests would be needed for a location to determine the final saturated value of k. In the case of TI5, as the values decrease for the first tests, the land tends to saturation, and then they increase, probably due to water movement following new pathways, and finally decrease again as the saturation increases (Figure 3).
At the TI2 and TI4 locations, the hydraulic conductivity decreases as the number of tests increases and the water infiltration becomes stable due to the ground saturation. In TI1 and TI5, the k values first decrease, then rise, and finally again decrease as the water movement tends to a stationary regime. The decrease in k values followed by an increase can be attributed to the obstruction of pores and superficial fissures, to the presence of air, and to the soil heterogeneity, the water pathways being subsequently unlocked.
In TI3 and TI6 in particular, the scatter of values is much higher between consecutive trials, indicating a general increase in k values as the number of trials increases. This means that reaching the steady state regime is more difficult as the water infiltration pathways close and re-open by mobilizing and further removing fine particles. The first trials highlight smaller k values, probably due to the entrapped air. Finally, both curves indicate a decrease in hydraulic conductivity and possibly a tendency towards stabilization of the flow. Due to the relatively small surface of the studied area and the short distance between the test points, it can be argued that the heterogeneity of the urban soils leads to large variations of hydraulic conductivity.
Figure 3 shows an area delineated by two horizontal lines. The upper one represents the arithmetic average of the hypothetical saturation state of ki values obtained for each test (kavg, column 21 of Table 1) and the lower represents the value of the hydraulic conductivity obtained by the mean of the hydraulic head values and of the mean level oscillation values (kdcH, column 20 of Table 1). This interval was defined as “the saturated hydraulic conductivity validity domain”.
If we accept that all k values above the upper line are overvalued and those below the bottom line are underestimated, it means that the value that could be considered for each location would be in the range delimited by the two parallel lines, being the average of the kavg and kdcH values. Since, in the boundary range, there are usually more values whose distribution is not central, the value of hydraulic conductivity was calculated using a linear regression for the selected values: y = ax + b, where a represents the slope, b the intersection with the ordinate, y the value of k, and x the number of the selected tests. The equations of the regression curves for selected values within the validity domain, as it was defined, are given in Figure 3.
The value of kss hydraulic conductivity determined at each TI location represents the average of the values determined by the regression curve for the selected trials (within the interval bordered by the two parallel lines).
The number of trials at each location is important. For example, if, at the TI1 and TI4 locations, where the values obtained for k between the first and the last trials show very large differences, 5–6 attempts would have been performed, then the error of the calculated kss values should have been quite high. The required number of trials to reach the steady state flow cannot be known at the beginning of the tests. This can be estimated only after performing the tests, computing k values and determining the validity domain.
As the k value decreases when the number of trials (and consequently the ground saturation) increases, there should be an inverse correlation between the values obtained for k and the number of trials. This correlation can be seen in Figure 3 for all the locations with the exception of TI3 and TI6. For these two, a direct correlation can be observed (the k values increase with the increase in the tests number). From the physical point of view, these data sets should be considered as not being complete and a larger number of trials should be made. Also, for TI3, the regression curve for the selected and simulated values (used to compute the kss) shows an inverse slope, indicating an inverse correlation. Due to the small number of available values, the kss value is expected to change if more tests are performed. For TI6, the regression line indicates a direct correlation, but its position is very close to the horizontal. A horizontal line should probably have been obtained with a few more attempts. From a strict statistical point of view, this shows a lack of correlation, but the physical interpretation indicates the achievement of incipient steady state for the infiltrated water. This is also valid for TI1 and TI4 locations, where the regression lines are close to the horizontal.

3.2. Inversed Auger (IA)

The hydraulic conductivity values were computed using Equation (4) and the results are shown in Table 2. The field recorded data and the corresponding graphs are presented in Supplementary S5 on ESM.
Similar to the TI tests, the hydraulic conductivity values kdch (column 19 in Table 2) differ from the arithmetic mean kavg (column 20, Table 2) and they are closer to the k values determined for the last trials at each location. This leads to the interpretation that the terrain was saturated or is in the early stage of so-called saturation. On the other hand, the last determinations of k values in the row corresponding to each location (especially for locations with a greater number of trials) show oscillating values when compared to the neighboring ones. The variation of the values may be due to the change in the structure of the ground, either due to its heterogeneity or due to lower cohesion than of the natural terrains, and the entrainment of the soil particles that can clog different access paths which are subsequently released.
It should also be considered that the water also infiltrates through the lateral surface of the drilling hole, and this differs from one location to the other: the deeper boreholes open places with a different makeup of the clay matrix. Consequently, the hydraulic conductivity corresponds to the tested depth range for each location. Figure 4 illustrates the changes of hydraulic conductivity as a function of the number of tests for each location.
As with the TI tests, in Figure 4, a validity domain is shown of the hydraulic conductivity values determined between kdch and kavg (the lowest and highest values). An exception is the location IA1, where only two trials were performed (whose values are outside the limits of the validity domain). We point out that tests at all locations were made in the same period of time, and the cumulative durations of the two IA1 trials exceed the cumulative duration of the tests in other locations. Therefore, in the case of IA1, the calculation value can be chosen as the average of the two values that delineate the validity domain. Generally, an inverse correlation is found between the hydraulic conductivity and the number of tests (a decrease in the hydraulic conductivity values with the number of tests).
In the locations IA3 and IA4, for T11 and T9 trials, respectively, there is a sudden increase in the k value of more than one order of magnitude. The hydraulic head applied in T11 of IA3 (Table 2) represents approximately 50% of the hydraulic head applied to other trials. In the case of T9 of the IA4 location, the hydraulic head is very close to the hydraulic head applied to the other tests. A high hydraulic conductivity variation could be the result of the infiltration process progression due to the action of the water on the heterogeneous media by mobilization/demobilization of the fine particles (the depths of the drill holes are 30 cm and 29 cm, Figure 2b). Also, the lithological constitution of the clayey matrix has to be considered with the presence of brick debris, splinters, or thickening elements along the walls of the borehole. According to Supplementary S1, the samples taken around the depth of 30 cm indicate a matrix made up of silty clay, but the interpretation based on the macroscopic description of the excavated ground corroborated the granulometry of the adjacent locations, indicating the presence of a sandy clay on the surface with the silty clay located at the bottom. To lower the hydraulic head by 50%, the hydraulic conductivity should be characteristic of the lower horizon, but its variation is far too large, so the lithological differences are not so important.
The value of hydraulic conductivity kss corresponding to each IA location (Table 3) represents the average of the values determined by the regression line for the selected trials (trials within the boundary bounded by the two parallel lines) except for IA1, where the calculated value kss was determined as the mean of kavg and kdch values, kss = 0.114 m/day.
Figure 4 shows that some values (especially those related to the final IA trials) are below the lower limit of the validity domain (below kdch value), with a downward but oscillating trend. Apparently, the elimination of these values would lead to an erroneous evaluation of kss, but it should be considered that the kdch value depends on the number of trials and averages of the hydraulic head changes. Therefore, its position in the graph changes with the number of trials, and the values between the two limits of the range indicate the tendency toward the beginning of a steady state regime. A clear example is represented by IA3, where the value obtained at the last trial rises within the domain, while the previous four are below the lower limit of the validity domain with a partially downward trend. The regression line of the selected values is very close to the horizontal, indicating the entry into the steady state. In locations IA2, IA4, and IA6, the large angle between the regression line and the horizontal indicates that the steady state has not been reached. However, for IA6, the close and slightly oscillating values of the last trials indicate the approach of early stage of steady state flow.
For both test methods (TI and IA), using the values of Table 1 and Table 2, the following were computed for each location: the standard deviations so of the values obtained for k for all tests, the standard deviations sso for the k values selected within the validity domain, and the standard deviations sss for the k values computed/simulated with the regression curve equation for the selected trials together with the corresponding coefficients of variation Vo, Vso, Vss. These parameters, together with the values of the hydraulic conductivity, are presented in Table 3 as a simple average of the selected values ks and as an average of the values simulated by the kss curve in the validity domain.
It can be observed that for each location, the standard deviation and the variation coefficient values decrease from the obtained group of values to the selected ones. The smallest values are for the simulated group of values and signify that the magnitude, as well as the uncertainty, associated with the attempts to determine the k value, decrease if the values defined on the validity domain are considered.
The decrease in the variation coefficient suggests that the simulated k values are closer to real values as an expression of geological environment variability. By comparing the statistical parameters s and V, we see a large variation from one location to another. This mainly indicates a different behavior of infiltrated water related to the heterogeneity of the ground (lithological differences, clogging/opening of access ways, air presence with implication on the saturation degree, etc.).
This is also emphasized by the fact that for the selected group of values, the statistical parameters are much closer, showing a decrease in the variability of the geological factors as the infiltration process approaches the “steady state”. The situation found in the TI4 location is significant, where the calculated statistical parameters have values greater than one and the differences between the k values (obtained for each of the 14 trials) are between one and two orders of magnitude. TI graph 4 of Figure 3 shows a normal evolution as the number of tests increases, indicating a gradual decrease in the k value with an increase in saturation. It can be considered that the excessively high values of k obtained for the first attempts are related to the flow of water through a dry soil with cracks with larger openings as well to vegetation. As voids fill and water interacts with the clay fraction, which can lead to swelling, the infiltration velocity is attenuated, leading to a decrease in k. In locations IA3 and IA4, the coefficients of variation for the group of the observed values are very high due to the results obtained in the T11 and T9 trials for the respective locations. By selecting the k values from the validity domain, in the cases of IA2 and IA6, the values of ks and kss are equal because only two values have been identified in this field with which the regression line was built.
Structural and lithological heterogeneity can influence the seepage process through the unsaturated zone, but this process is also related to the relatively short time of the contact between the water and the lithological/mineralogical environment. As consequence, it is not appropriate to check if the data (in this case, the k values determined by n trials on each location) come from the same geological set on the basis of their variance, as the urban soil is a mixture of natural and anthropogenic materials, chaotically mixed.
For comparison, shown in Table 4 are the values of hydraulic conductivity that define the validity domain (kdcH and kavg) for each location, the average of these values km (as the central value of the domain), and the hydraulic conductivity value obtained on the basis of the selected values kss according to the regression line. The kss hydraulic conductivity value is considered as the closest to the real saturated value and is used below as k.
With the exception of TI4 and TI6, all the km values are quite close to kss values determined by linear regression. It is obvious that by applying a polynomial regression curve of second degree or even higher, with a better degree of approximation, these differences would change. Therefore, a rapid evaluation can be made by calculating the value of km. However, it is more correct to determine the kss values by applying the appropriate regression curve.

3.3. Double Ring Infiltrometer (DRI)

To assess the hydraulic conductivity from field records of cumulative infiltration I and infiltration rate v, Equation (2) was used. The field data and the corresponding graphs I = f (t) and v = f (t) are shown in Supplementary S3 on ESM. The saturated hydraulic conductivity values are shown in Table 5. These can be considered as average values of the locations since the tests performed with a constant level were completed at steady state, so that at each location marked in Figure 2b, a single test was performed.
The k values obtained by this method are greater than 1 (m/d) and have a relatively large variation from one location to another.

3.4. Minidisk Infiltrometer (MDI)

The hydraulic conductivity value k was determined by applying Equation (3) and using the field data. The graphs I = f t and the corresponding regression curves are shown in Supplementary S4 on ESM. The values obtained for the hydraulic conductivity k are shown by Table 5. Compared to those obtained by the above-mentioned methods (TI and DRI), the values differ by up to three orders of magnitude.
An explanation of the very low values obtained for hydraulic conductivity would be that the infiltration area, being very small (device feature), largely overlaps the clay matrix of the filler material. An analysis of the values in relation to infiltration area is presented in the following sections.

4. Discussion

4.1. Comparative Analysis of the Results Obtained by the Four Testing Methods

Figure 2b shows a diagrammatic map of the test area with locations of the different types of tests. Positions of these locations can be measured as distance along an axis from the lower left to lower right, with the origin on the left. A plot of these distances against k values of each test location is shown in Figure 5a.
It can be observed that the hydraulic conductivity values assessed by applying the DRI and IA methods have the largest amplitude between the test locations. The values determined by the MDI method have very small amplitude and its average can be considered as the characteristic value of the studied area.
When the measured values are represented in ascending order of the infiltration areas corresponding to each test method (Table 5), the k amplitude increases with the increase in the infiltration area, which supports the conclusion that the values obtained by the MDI method are the appropriate ones. However, from the physical point of view, an experiment performed on wider surfaces should lead to results that characterize the investigated zone better. If, apart from the geological variable k, we define the variable k* = k/A which represents the hydraulic conductivity corresponding to a unitary infiltration surface (equal to 1 cm2) and reconstruct the graph of Figure 5a with the k/A values of Table 5, the distribution in Figure 5b is obtained.
Comparing the distributions shown by Figure 5a,b, the pattern of values is inverted in the last figure. The hydraulic conductivity corresponding to the surface infiltration unit shows the smallest amplitudes in the case of the IA and DRI methods. The maximum amplitudes of k/A were obtained for the MDI and TI methods. This can be explained as due to the fact that by applying the DRI and IA methods, the water infiltrates through a much larger surface area compared to the other two methods, so the obtained values for the saturated hydraulic conductivity by the DRI and IA methods characterize the tested locations better by increasing the estimation confidence. This is particularly important in assessing the hydraulic conductivity for urban soils, as they are most often made of heterogeneous material. The small values of the saturated hydraulic conductivity obtained by the MDI method (where the infiltration surface is the smallest and the tests were carried out with a negative pressure) would rather characterize the clay matrix of the material that forms the urban soil and not the urban soil in its entirety.
If the terrain had been ideally “homogeneous” (hypothesis “0” stated above), the hydraulic conductivity for the saturated state, determined by any of the methods, should have been the same and the k/A ratio should have decreased with increasing infiltration surface area. The amplitude of this ratio should be the smallest for the methods with the largest infiltration surface, and the graphical representation of these values (according to Figure 5b) should be represented by parallel horizontal lines for equal infiltration surfaces for the same method. In the present case (heterogeneous soil), the k values differ significantly from one location to another within the same test method (Figure 5a) given that the respective locations are very close; the compared values of the amplitudes of the k/A ratios (Figure 5b) seem to confirm the hypothesis that their value decreases with increasing infiltration area. However, in particular, the variations of k/A values between close locations, determined by the method with large and equal infiltration areas (DRI), indicate the presence of other variables independent of the infiltration area that can influence the saturated k value. Moreover, in the case of MDI tests from close locations (MDI14 and MDI9), where the determined k value is a characteristic of the clay matrix, the amplitude between the k/A values (A = ct.) is very high. This indicates either a lithological/granulometric variation or the existence of other variables (roots, macropores below the ground surface, etc.), or even the cumulative effect of these variables that could not have been highlighted at the beginning of the tests. Therefore, in the case of performing several tests, considering that k is dependent on the size of the infiltration area, the k/A ratio provides qualitative information on the influence of other variables on the k values determined by the magnitude of the amplitudes of the k/A ratios. The influence of the infiltration area on the determined saturated k values, as well as of other cumulated variables, is presented in the following sections.

4.2. Considerations Regarding the Average Values Obtained for k and k/A Ratio as Geological Variables Related to Each Test Method for Both Data Selection Categories

The arithmetic average of the hydraulic conductivity values for the entire tested area, considering the values obtained by applying all the four methods on all the locations (not grouped by selection intervals), is k = 1.058 m/day (Table 5). This value is very close to the one obtained by applying the TI method (1.062 m/day). The one corresponding to the DRI method (3.687 m/day) is the highest value, and the one determined by the MDI method (0.160 m/day) shows the lowest value (almost seven times smaller as general average). The amplitude values confirm the situations presented by Figure 5. Significant differences between the averages obtained by each method are observed. Apparently, according to the values of the obtained averages, the simple arithmetic average of all the k values selected for the four methods or the average obtained by applying the TI method would best characterize the studied site. The smallest amplitude of the k values was obtained by the MDI method (0.467 m/day), with the one obtained by the DRI (4.094 m/day) being the highest. The amplitude resulting from the TI method (2.548 m/day) occupies an intermediate position, which means that the lowest spread of values is obtained by the MDI method.
Figure 6 shows the relationships between the calculated k and k/A averages and the variation coefficients for the grouped and non-grouped selection intervals and for each testing method, respectively.
The lowest variation coefficient was obtained for the DRI method, both for the selection mean of the hydraulic conductivity value kDRI and for the hydraulic conductivity value corresponding to an infiltration surface unit k* = kDRI/ADRI. At the same time, in the case of ungrouped elements, the variation coefficient exceeds the value of 0.5 by very little, considered statistically a maximum value for a normal variation of the analyzed parameter. In the case of grouped elements, the value of the variation coefficient decreases below this limit. The variation coefficient determined on the basis of the data provided by the MDI method is located in an intermediary position between those obtained for TI, IA (highest values), and DRI (lowest value). These values show that the comparative analysis of the results made only on the basis of the magnitude of the amplitude between the maximum and minimum values can lead to erroneous conclusions.
As a general observation, the graphs from Figure 5 and Figure 6 indicate a large variation of the values of k and k/A parameters obtained by different methods.

4.3. The Correlation Analysis between the Hydraulic Conductivity and the Infiltration Area

The influence of infiltration area A as independent variable on the dependent variables D (k and k/A ratio) was examined. Using the results of Table, 5 the general variances s2D and the conditioned variances s2D(A) of the dependent variables were calculated and are shown in Table 6. As a measure of the variation of the dependent variables under the influence of other factors, the residual variance s2Do was determined. The statistical parameters which indicate the strength and type of correlations, the correlation ratio ηD(A) and correlation coefficient r, were also calculated. The estimation of the certainty of the linear or non-linear regression was evaluated by the parabolic regression coefficient Crp, defined by the dispersion of theoretical mean values (simulated) D’avg(A) (Table 6) compared to the mean of the observed values Davg (Table 5). The parabolic regression is sufficiently certain if CrpηD(A) ≠ ηA(D) ≠ |r|. In this work, only ηD(A) and r were computed due to the fact that it should be incorrect to consider that the infiltration area of the device depends on the hydraulic conductivity.
For each dependent variable, the analysis was performed for two situations: for all four applied methods and only for methods that tested the terrain surface (without IA method). Related to the type of analyzed method and according to the values highlighted in Table 6, the following are noted:
Case (a)—The correlation coefficient and the correlation ratio indicate a direct and quite strong correlation between k and A. The residual variance s2ko represents about 25% of the general variance s2k. The fact that the conditioned variance represents 75% of the general variance indicates the weight attributed to the independent variable (infiltration area) in determining the fluctuations of the dependent variable (saturated hydraulic conductivity). The scatter diagram of Figure 7a suggests a parabolic correlation.
Although the above mentioned condition regarding the Crp test value is not respected, due to the very small difference between the compared values (Crp = 0.856 < ηD(A) = 0.864), the parabolic correlation can be accepted, taking in consideration that the residual dispersion represents 25% of the general dispersion and ηD(A) = 0.864 ≠ |r|= 0.470. The conditioned means values of the dependent variable and of the selection interval centers for the independent variable were used to draw the graph of Figure 7b (for the grouped observed values). By using the coefficients of the regression curve (parabolic form), the theoretical (simulated) average values of the dependent variable k’avg(A) were obtained (against the observed average values kavg). The curve of the simulated values was constructed, keeping the same values for the infiltration area as in the first curve.
It can be seen that the regression curve equation of Figure 7b has a maximum. If we set to zero the first derivative of the respective equation, we obtain the pair of values of the maximum, namely A = 636.364 cm2 and k = 2.627 m/day. Therefore, the value of the hydraulic conductivity increases to a maximum that corresponds to an infiltration area slightly larger than that used by the DRI method. Then, it decreases with the increase in the infiltration area, suggesting an inverse correlation. Confirmation of this can be made by analyzing the test procedures corresponding to the test methods. In the MDI, TI, and DRI methods, the water infiltration takes place through the base area of the respective devices. In the case of the IA method, infiltration is also taking place through the lateral surface and the measured values characterize the urban soil at a greater depth in contrast to the other three methods where the measured values characterize only the superficial zone. In this respect, the correlation between the hydraulic conductivity and the infiltration area will be analyzed in Case (b) only for the MDI, TI, and DRI methods.
Case (b)—By eliminating IA values, the scatter diagram of Figure 8a could suggest a nonlinear correlation as well as a linear one. The linearity or non-linearity of the correlation was estimated with the values of the statistical parameters (Table 6).
The results indicate a strong and direct correlation between hydraulic conductivity and infiltration area. The correlation coefficient r = 0.831 has twice the magnitude of the previous case, which indicates the strength of the direct correlation. The tiny difference between the correlation ratio and the correlation coefficient shows that it could be also considered as a linear correlation. However, the value of the parabolic regression coefficient being greater than the correlation ratio Crp = 0.991 ≥ ηD(A) = 0.844 shows that the parabolic interpretation approximates much better the regression curve of Figure 8b, where the graphs were prepared according to the methodology previously presented.
The nonlinear dependencies between variables appear when the change of the values of the independent variable determines the non-uniform modification (“accelerated” or “slowed down”) of the dependent variable. In situations where the coefficients r and η have equal or very close values (the distinction between linear and parabolic correlation being made on the basis of Crpη inequality as in the present case), adjusting the observation data (a parabola instead of a regression line) does not provide informational advantages. In the analyzed situation, the significance of the parabolic correlation should not be interpreted in a geological sense, as the infiltration area is not a geological variable. The parabolic correlation is due to the large variation of the infiltration areas of the three test methods (TI, DRI, and MDI) as well as to the variation of these areas within the TI method, combined with the relatively small number of DRI tests.
The regression curve of Figure 8 indicates that in the study situation, the hydraulic conductivity increases without a limit as the infiltration area increases. It should be kept in mind that the graph of Figure 8 depicts only this experiment and the results depend on the number of tests performed by each method and consequently by the variation of the infiltration area. In the case of the DRI method, only five tests with the constant infiltration area were performed. If more tests had been performed with variable infiltration areas, a decrease in the slope of the regression curve should have been expected (tending towards the horizontal) as the infiltration area increases. The residual variance represents about 28% of the general variance and apparently, as in case (a), the infiltration area seems to be the main factor influencing the k value. However, the conditioned variance of 72% of general variance is related to the devices’ infiltration areas which cover different surfaces with different characteristics (e.g., cracks). To estimate if the infiltration area is the main factor which influences the k value, it is necessary to refer to a unit infiltration area, as in cases (c) and (d).
Case (c)—The correlation analysis between the ratio of the hydraulic conductivity value and the infiltration area k/A and the size of the infiltration area A was performed by the same methodology for all methods (Table 6 case (c)). The graph of Figure 9a, constructed with the ungrouped values of the parameters k/A and A, suggests a possible linear correlation between these two variables. The correlation ratio indicates a connection between k/A and the infiltration area A. The correlation coefficient, by its negative value, shows an inverse relationship between the two variables. That means that with increasing infiltration area, the value of the hydraulic conductivity related to a surface unit decreases. There is a higher degree of confidence for the values obtained with the test devices having larger infiltration areas. The value of the parabolic regression coefficient being much higher than that of the correlation coefficient Crp = 0.839 ≥ ηD(A) = 0.577, signify that the parabolic regression (Figure 9b) better approximates the relation between the two variables. The conditioned and residual variances represent 33% and 67% of the general variance, respectively. This means that other factors which were not considered at the beginning of the tests have a strong influence on the infiltration process and the comparative analysis between different methods of the k/A ratio values better emphasizes this influence. In this way, the differences in amplitude of the graphs from Figure 5 can be better explained. On the other hand, the decreasing values on k/A ratio with increasing A values (Figure 9b) confirm hypothesis “0”, stated above, that the infiltration area is the main factor which influences k values and none of the other factors, taken separately, has a greater weight than 33%, and that the average value of the saturated hydraulic conductivity for all methods kavg = 1.058 m/d (Table 5) better approximates the analyzed REV.
Case (d)—Keeping only the values of the k/A ratio for the test methods that characterize the terrain surface (TI, DRI and MDI) and eliminating the values for the IA method, the statistical parameters are presented in Table 6 case (d). The value of the correlation coefficient, r = −0.237, is about half that where all test methods were considered and indicates a poor inverse correlation between the two variables. Since the r value can be influenced by other variables, the strength of a correlation is better highlighted by the correlation ratio ηD(A) = 0.434, which indicates a moderate correlation between k/A ratio and infiltration area A for these methods. In conjunction with the value of the residual variance s2Do= 0.0000554, which represents 81% of the general dispersion s2D= 0.0000682, it means that the TI, DRI, and MDI methods are strongly influenced by other variables (roots, superficial cracks, wormholes, etc.). In Figure 10a, the ungrouped mean values of k/A ratio suggest a parabolic correlation between independent and dependent variables, which is confirmed by Figure 10b. The value of the parabolic regression coefficient is Crp = 1.
If we compare the percentage values of unconditioned variances for the last two analyzed situations for the k/A ratio (67% for all methods and 81% for testing methods only on the land surface), the following observations can be made: (1) the lower value obtained in the first case is due to the fact that the influence of the other variables is attenuated by the values obtained by the IA method, where the drying cracks or roots have a lower weight with increasing depth; (2) starting from the hypothesis “0”, it is observed that even for heterogeneous soil, the value of the ratio k/A decreases with the increase in the infiltration area (Figure 9). Even if the cumulative influence of other variables is 67%, all of them still exceed the value of the conditioned variable (33%), the infiltration area being the independent variable with the most important weight. In the second case, where the share of unconditioned variables is very high (81%), in Figure 10b, it is observed that on the left side of the graph, the value of the k/A ratio increases with the increase in the infiltration area, which confirms hypothesis “1”. This means that for values lower than 350 cm2 (approximately), the infiltration area no longer represents the predominant independent variable. The right side of the curve described by the simulated values indicates a decrease in the k/A ratio with the increase in the infiltration area. The infiltration area becomes the main independent variable, the cumulative effect of the other unconditioned variables being attenuated by the increase in the infiltration surface, confirming hypothesis “0”. In this sense, the determination of the independent variable with the largest weight (on the left side of the graph) can be achieved by multiple regression analysis, which is not part of this study. Based on pumping tests on natural saturated and fissured (structurally heterogeneous) rocks, Rovey and Charles [40] showed that in order to obtain a characteristic value of the hydraulic conductivity, it is necessary that the surface on which the tests are performed must be representative from the point of view of heterogeneity, so it must be higher than that considered statistically “homogeneous”. In this study, we demonstrated this by analyzing the variation of k/A ratio values versus the variation of A values for anthropogenic soil starting from a unit surface (1 cm2) hypothetically considered homogeneous. When the k/A values decrease with increasing A values, then the variation of soil “homogeneity” from the microscale to the macroscale is insignificant; vice versa, when the values of the ratio k/A increase with the increase in A, then this variation is significant.
It should be noted that the graphs in Figure 8 and Figure 10 have an indicative and explanatory character because the number of values representing the infiltration area was diminished by extracting those corresponding to the IA method. Those values covered a range between 508 and 1303 cm2 (Table 5), the series of values being very narrow and represented by those methods with very small area (MDI) or much larger (DRI). This has led to a small number of grouping intervals with unequal dimensions, which does not ensure accuracy of the simulated values. At the same time, the domains defined by the graphs should not be extrapolated because, even in the case of more data where the values of the infiltration area would have a uniform distribution, extrapolation is limited by the physical meaning of the variables. Thus, in the graph of Figure 10b, for values greater than 800 cm2 of the infiltration area, the value of the k/A ratio will not take negative values as the value of k can never be less than zero.

5. Conclusions

Four methods (TI, DRI, MDI, and IA), with different infiltration areas A of the devices, were used to quantify the saturated hydraulic conductivity k for an anthropogenic unsaturated soil consisting of sand, gravel, boulder elements, brick residues, glass, and concrete, trapped in a clayey matrix with roots of plants at the terrain surface. The ratio between test perimeter length and width is 17.7 m/3.81 m, so that the close locations of the different test methods led to a relative compensation to the effect of soil heterogeneity [38]. The average values of the saturated hydraulic conductivity obtained by the four test methods show significant differences between them. For the TI and IA methods, the saturated hydraulic conductivity values k determined for ungrouped values are relatively close (kTI = 1.062 m/day and kIA = 1.56 m/day). Comparing them to DRI values, they are more than two times lower (kDRI = 3.687 m/day), and comparing to MDI values, they are higher by one order of magnitude (kMDI = 0.160 m/day). For the TI and IA tests, several attempts were performed in the same location. Since the test procedure does not allow the assessment of the degree of saturation of the soil without disturbing it, this assessment was carried out after performing the tests by defining the domain of validity. The lower limit of the domain is represented by the average values of hydraulic conductivity determined on the basis of the average value of the hydraulic head variation. The upper limit represents the simple arithmetic mean of the observed values. The values above the upper limits are either considered overvalued or correspond to the unsteady state flow. Those below the lower limit are undervalued and could be explained either by temporary blocking of the water flow paths or due to the interaction of clay particles with the water, which can lead to the swelling of the clay material, or by the presence of entrapped air. For the selected hydraulic conductivity values inside the domain of validity, the values of standard deviation and of the variation coefficient are considerably lower compared to those for the observed values. The domain of validity presents the advantage that the final obtained values express the evolution of the physical phenomenon itself and does not consider only the statistically determined confidence interval on a range of values where, sometimes, the lower limit can be negative, which physically is not possible. The regression line inside the validity area indicates a steady state water seepage when it is horizontal or a steady state initial stage when the gradient is close to zero (answer to question (1) raised in the Introduction section).
The amplitudes calculated between the average values of the saturated state hydraulic conductivity, obtained for each method, present the highest values for the devices with the largest infiltration surfaces. Apparently, this would contradict the hypothesis, stated by various authors, that the value determined for the hydraulic conductivity k (saturated state) approaches the real value as the infiltration surface A increases. In this sense, a new parameter k* = k/A was introduced which represents the hydraulic conductivity corresponding to a surface unit (1 cm2) ideally considered homogeneous. Thus, the lowest values of the amplitudes of the k/A ratio, determined for each test method, were obtained for the methods with the largest test area, DRI and IA, and the highest values for the methods with the smallest surface areas of infiltration, MDI and TI (answer to question (2)). The statistical significance of the amplitude of the k/A ratio is in accordance with the physical reality. The decrease in the k/A value with the increase in the A value, within the same method or between different methods, indicates that the infiltration area represents the main independent variable that controls the variability of the hydraulic conductivity. The increase in the k/A value with the increase in the A value indicates that the infiltration area no longer represents the main independent variable on which the k value depends for the saturated state. With other words, in the first case, the variation of the terrain homogeneity is not significant from a small scale to a large scale, while in the second case, this variation is significant and there are one or more variables with a greater influence on the k value, such as cracks, roots, etc. (answer to question (3)).
The results obtained in this study should be applied to select a method from the point of view of its applicability function of the urban soil composition and should not be interpreted as confirming or invalidating the utility of one or another method. Correct evaluation of the saturated hydraulic conductivity has multiple urban uses and a correct lithological recognition is needed to choose a test method on unsaturated urban soil. For example, to evaluate the saturated hydraulic conductivity for a supporting layer (considered homogeneous) below a pervious pavement, devices with small infiltration area could be used. When evaluating shallow aquifer recharge or hypodermic flow, devices with large infiltration surfaces must be used and the obtained results must be based on a large number of tests imposed by the heterogeneity of the anthropogenic terrain (answer to question (4)).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w16131908/s1.

Author Contributions

Conceptualization, T.G. and D.G.; Methodology, T.G.; Validation, T.G. and D.G.; Formal analysis, T.G.; Investigation, T.G. and D.G.; Resources, D.G.; Data curation, T.G.; Writing—original draft, T.G.; Writing—review & editing, T.G. and D.G.; Supervision, D.G.; Funding acquisition, D.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported through the HORIZON-CL6-2023-CLIMATE-01-2, AWARD research project “Alternative Water Resources and Deliberation process to renew Water Supply Strategic Planning” Grant agreement ID: 101136987.

Data Availability Statement

Data is contained within the article (and Supplementary Materials).

Acknowledgments

The authors are grateful to James Petch, Emeritus Professor of the University of Manchester, and to Carl Axness, former Senior member of Sandia National Laboratories, for the valuable comments they made on this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dijkstra, J.J.; Comans, R.N.; Schokker, J.; van der Meulen, M.J. The geological significance of novel anthropogenic materials: Deposits of industrial waste and by-products. Anthropocene 2019, 28, 100229. [Google Scholar] [CrossRef]
  2. Vasquez-Sune, E.; Sanchez-Vila, X.; Carrera, J. Introductory review of specific factors influencing urban groundwater, an emerging branch of hydrogeology with reference to Barcelona, Spain. Hydrogeol. J. 2005, 13, 522–533. [Google Scholar] [CrossRef]
  3. Boukhemacha, M.A.; Gogu, C.R.; Serpescu, I.; Gaitanaru, D.; Bica, I. A hydrogeological conceptual approach to stydy urban groundwater flow in Bucharest city, Romania. Hydrogeol. J. 2015, 23, 437–450. [Google Scholar] [CrossRef]
  4. Gogu, C.R.; Gaitanaru, D.; Boukhemacha, M.A.; Serpescu, I.; Litescu, L.; Zaharia, V.; Moldovan, A.; Mihailovici, M.J. Urban hydrogeology studies in Bucharest City, Romania. Procedia Eng. 2017, 209, 135–142. [Google Scholar] [CrossRef]
  5. Gogu, R.C.; Boukhemancha, M.A.; Gaitanaru, D.; Moraru, I. Interaction between City Subsurface Infrastructure and Groundwater. In Proceedings of the International Conference on Urban Drainage Modelling, Palermo, Italy, 23–26 September 2018. [Google Scholar]
  6. Morbidelli, R.; Saltalippi, C.; Flammini, A.; Cifrodelli, M.; Picciafuoco, T.; Corradini, C.; Govindaraju, R.S. In situ measurements of soil saturated hydraulic conductivity: Assessment of reliability through rainfall–runoff experiments. Hydrol. Process. 2017, 31, 3084–3094. [Google Scholar] [CrossRef]
  7. Youngs, E.G. Estimating hydraulic conductivity values from ring infiltrometer measurements. J. Soil Sci. 1987, 38, 623–632. [Google Scholar] [CrossRef]
  8. Reynolds, W.D.; Elrick, D.E. Chapter 6: Measurement and Characterization of Soil Hydraulic Properties. In Soil-Water-Solute Process Characterization: An Integrated Approach; CRC Press: Boca Raton, FL, USA, 2005. [Google Scholar]
  9. Lai, J.; Ren, L. Assessing the size dependency of measured hydraulic conductivity using double-ring infiltrometers and numerical simulation. Soil Sci. Soc. Am. J. 2007, 71, 1667–1675. [Google Scholar] [CrossRef]
  10. Lai, J.; Luo, Y.; Ren, L. Buffer index effects on hydraulic conductivity measurements using numerical simulations of double-ring infiltration. Soil Sci. Soc. Am. J. 2010, 74, 1526–1536. [Google Scholar] [CrossRef]
  11. Köhne, J.M.; Alves, J.J.; Köhne, S.; Tiemer, B.; Lennartz, B.; Kruse, J. Double-ring and tension infiltrometer measurements of hydraulic conductivity and mobile soil regions. Pesqui. Agropecuária Trop. 2011, 41, 336–347. [Google Scholar]
  12. Fallico, C.; Migliari, E.; Troisi, S. Comparison of three measurement methods of saturated hydraulic condutivity. Hydrol. Earth Syst. Sci. Discuss. 2006, 3, 987–1019. [Google Scholar]
  13. Lee, D.M.; Elrick, D.E.; Reynolds, W.D.; Clothier, B.E. A comparison of three field methods for measuring saturated hydraulic conductivity. Can. J. Soil Sci. 1985, 65, 563–573. [Google Scholar] [CrossRef]
  14. Reynolds, W.D.; Bowman, B.T.; Brunke, R.R.; Drury, C.F.; Tan, C.S. Comparison of tension infiltrometer, pressure infiltrometer, and soil core estimates of saturated hydraulic conductivity. Soil Sci. Soc. Am. J. 2000, 64, 478–484. [Google Scholar] [CrossRef]
  15. Beven, K.; Germann, P. Macropores and water flow in soils. Water Resour. Res. 1982, 18, 1311–1325. [Google Scholar] [CrossRef]
  16. Van Hoorn, J.W. Determining hydraulic conductivity with the inversed auger hole and infiltrometer methods. In Proceedings of the International Drainage Workshop, Wageningen, The Netherlands, 16–20 May 1978; Wesseling, J., Ed.; ILRI publication: Nairobi, Kenya, 1979; pp. 150–154. [Google Scholar]
  17. Bagarello, V.; Iovino, M.; Lai, J. Testing steady-state analysis of single-ring and square pressure infiltrometer data. Geoderma 2015, 261, 101–109. [Google Scholar] [CrossRef]
  18. Bockheim, J.G. Nature and Properties of Highly Disturbed Soils, Philadelphia, Pennsylvania; Paper Presented before Div. S-5; Soil Science Society of America: Chicago, IL, USA, 1974. [Google Scholar]
  19. Craul, P. A description of urban soils and their desired characteristics. J. Soc. Arboric. 1985, 11, 330–339. [Google Scholar] [CrossRef]
  20. Meng, X.; Liao, H.; Zhang, J. Infiltration law of water in undisturbed loess and backfill. Water 2020, 12, 2388. [Google Scholar] [CrossRef]
  21. Batsilas, I.; Angelaki, A.; Chalkidis, I. Hydrodynamics of the Vadose Zone of a Layered Soil Column. Water 2023, 15, 221. [Google Scholar] [CrossRef]
  22. Císlerová, M.; Šimůnek, J.; Vogel, T. Changes of steady-state infiltration rates in recurrent ponding infiltration experiments. J. Hydrol. 1988, 104, 1–16. [Google Scholar] [CrossRef]
  23. Sakaguchi, A.; Nishimura, T.; Kato, M. The effect of entrapped air on the quasi-saturated soil hydraulic conductivity and comparison with the unsaturated hydraulic conductivity. Vadose Zone J. 2005, 4, 139–144. [Google Scholar] [CrossRef]
  24. Dikinya, O.; Hinz, C.; Aylmore, G. Decrease in hydraulic conductivity and particle release associated with self-filtration in saturated soil columns. Geoderma 2008, 146, 192–200. [Google Scholar] [CrossRef]
  25. Bomboe, P. Geologie Matematica. Volumul I. Analiza Statistica a Datelor Geologice (Mathematical Geology. Volume I. Statistical Analysisis of the Geological Data); University of Bucharest: Bucharest, Romania, 1979. [Google Scholar]
  26. Wu, L.; Pan, L. A generalized solution to infiltration from single-ring infiltrometers by scaling. Soil Sci. Soc. Am. J. 1997, 61, 1318–1322. [Google Scholar] [CrossRef]
  27. Wu, L.; Pan, L.; Mitchell, J.; Sanden, B. Measuring saturated hydraulic conductivity using a generalized solution for single-ring infiltrometers. Soil Sci. Soc. Am. J. 1999, 63, 788–792. [Google Scholar] [CrossRef]
  28. Reynolds, W.D.; Lewis, J.K. A drive point application of the Guelph Permeameter method for coarse-textured soils. Geoderma 2012, 187–188, 59–66. [Google Scholar] [CrossRef]
  29. Angulo-Jaramillo, R.; Bagarello, V.; Iovino, M.; Lassabatere, L. Infiltration Measurements for Soil Hydraulic Characterization; Springer: Berlin, Germany, 2016. [Google Scholar]
  30. Reynolds, W.D.; Elrick, D.E.; Youngs, E.G. Part 4. Single-ring and double or concentric-ring infiltrometers. In Method of Soil Analysis: Physical Methods; Soil Science Society of America, Book Series 5; Soil Science Society of America, Inc.: Madison, WI, USA, 2002. [Google Scholar]
  31. Philip, J.R. Numerical Solution of Equation of the Diffusion Type with Diffusivity Concentration Dependent. II. Aust. J. Phys. 1957, 10, 29–42. [Google Scholar] [CrossRef]
  32. Philip, J.R. The theory of infiltration: 1. The infiltration equation and its solution. Soil Sci. 1957, 83, 345–358. [Google Scholar] [CrossRef]
  33. Gupta, R.K.; Rudra, R.P.; Dickinson, W.T.; Patni, N.K.; Wall, G.J. Comparison of saturated hydraulic conductivity measured by various field methods. Trans. ASAE 1993, 36, 51–55. [Google Scholar] [CrossRef]
  34. Rose, C.W.; Stern, W.R. Determination of withdrawal of water from soil by crop roots as a function of depth and time. Soil Res. 1967, 5, 11–19. [Google Scholar] [CrossRef]
  35. Zhang, R. Determination of soil sorptivity and hydraulic conductivity from the disk infiltrometer. Sci. Soc. Am. J. 1997, 61, 1024–1030. [Google Scholar] [CrossRef]
  36. Carsel, R.; Parrish, R.S. Developing joint probability distributions of soil water retention characteristics. Water Resour. Res. 1988, 24, 755–769. [Google Scholar] [CrossRef]
  37. Decagon Devices, Inc. Mini Disk Infiltrometer User’s Manual; Decagon Devices, Inc.: Pullman, WA, USA, 2006. [Google Scholar]
  38. Bagarello, V.; Iovino, M.; Lai, J.B. Field and numerical tests of the two-ponding depth procedure for analysis of single-ring pressure infiltrometer data. Pedosphere 2013, 23, 779–789. [Google Scholar] [CrossRef]
  39. Reynolds, W.D. Chapter 56. Saturated hydraulic conductivity: Field measurement. In Soil Sampling and Methods of Analysis, Canadian Society of Soil Science; Lewis Publishers: Boca Raton, FL, USA, 1993; pp. 599–613. [Google Scholar]
  40. Rovey, I.I.; Charles, W. Digital simulation of the scale effect in hydraulic conductivity. Hydrogeol. J. 1998, 6, 216–225. [Google Scholar] [CrossRef]
Figure 1. The schemes of the test devices. (a) tube infiltrometer, TI: H = ct.—poured water column height, d—soil penetration depth, r—infiltration radius; (b) double ring infiltrometer, DRI: H = ct.—poured water height; (c) minidisk infiltrometer, MDI; (d) Inversed auger, IA: H-borehole depth, Ho—poured water level depth, Ht—water level depth at time t, ho—water column height at time t = 0, Δh—drawdown between step time ti and ti+1, ht—water column height at time ti+1.
Figure 1. The schemes of the test devices. (a) tube infiltrometer, TI: H = ct.—poured water column height, d—soil penetration depth, r—infiltration radius; (b) double ring infiltrometer, DRI: H = ct.—poured water height; (c) minidisk infiltrometer, MDI; (d) Inversed auger, IA: H-borehole depth, Ho—poured water level depth, Ht—water level depth at time t, ho—water column height at time t = 0, Δh—drawdown between step time ti and ti+1, ht—water column height at time ti+1.
Water 16 01908 g001
Figure 2. The experimental site: (a) the lithological column; (b) the distribution of the test locations and the applied methods.
Figure 2. The experimental site: (a) the lithological column; (b) the distribution of the test locations and the applied methods.
Water 16 01908 g002
Figure 3. The variation of hydraulic conductivity value vs. the number of trials at each TI test location.
Figure 3. The variation of hydraulic conductivity value vs. the number of trials at each TI test location.
Water 16 01908 g003
Figure 4. The variation of the hydraulic conductivity value in relationship to the number of trials on each IA location.
Figure 4. The variation of the hydraulic conductivity value in relationship to the number of trials on each IA location.
Water 16 01908 g004
Figure 5. The analyzed parameters distribution along the test area length: (a)—distribution of k values, (b)—distribution of k/A values.
Figure 5. The analyzed parameters distribution along the test area length: (a)—distribution of k values, (b)—distribution of k/A values.
Water 16 01908 g005
Figure 6. The relationship between the calculated means of k and k/A vs. variation coefficients for ungrouped (a,b) and grouped selections (c,d).
Figure 6. The relationship between the calculated means of k and k/A vs. variation coefficients for ungrouped (a,b) and grouped selections (c,d).
Water 16 01908 g006
Figure 7. The correlation between k and A values for all methods: (a)—ungrouped observed mean values; (b)—grouped observed and simulated mean values.
Figure 7. The correlation between k and A values for all methods: (a)—ungrouped observed mean values; (b)—grouped observed and simulated mean values.
Water 16 01908 g007
Figure 8. The correlation between k and A values for TI, DRI, and MDI methods: (a)—ungrouped observed mean values; (b)—grouped observed and simulated mean values.
Figure 8. The correlation between k and A values for TI, DRI, and MDI methods: (a)—ungrouped observed mean values; (b)—grouped observed and simulated mean values.
Water 16 01908 g008
Figure 9. The correlation between k/A ratio and A values for all methods: (a)—ungrouped observed mean values; (b)—grouped observed and simulated mean values.
Figure 9. The correlation between k/A ratio and A values for all methods: (a)—ungrouped observed mean values; (b)—grouped observed and simulated mean values.
Water 16 01908 g009
Figure 10. The correlation between k/A ratio and A values for the TI, DRI, and MDI methods: (a)—ungrouped observed mean values; (b)—grouped observed and simulated mean values.
Figure 10. The correlation between k/A ratio and A values for the TI, DRI, and MDI methods: (a)—ungrouped observed mean values; (b)—grouped observed and simulated mean values.
Water 16 01908 g010
Table 1. The k values obtained by the TI method.
Table 1. The k values obtained by the TI method.
LocationParameterTrialskdcH, m/dkavg,
m/d
T1T2T3T4T5T6T7T8T9T10T11T12T13T14T15T16T17T18T19
123456789101112131415161718192021
TI1Hn, cm30.442.6--47.346.442.547.346.745.845.744.445.945.647.846.1--------44.6--
k, m/d4.631.53--3.183.512.562.792.041.571.771.61.471.180.960.88--------0.9962.119
TI2Hn, cm55.45957.656.5------------------------------57.1--
k, m/d0.380.240.160.17------------------------------0.160.238
TI3Hn, cm26.330.13529.821.232.13332.729.626.731.1----------------29.7--
k, m/d0.350.340.370.820.711.200.940.981.891.291.25----------------0.3950.923
TI4Hn, cm23.825.428.927.930.730.227.130.232.430.130.131.831.930.7----------29.4--
k, m/d17.610.32.531.480.970.600.390.220.190.460.470.180.140.12----------0.0492.546
TI5Hn, cm49.950.356.15452.554.453.452.4----------------------52.9--
k, m/d0.690.610.460.980.770.650.420.23----------------------0.2370.600
TI6Hn, cm44.620.545.649.545.245.249.648.446.54249.248.848.445.848.850.648.448.542.845.7--
k, m/d1.384.762.442.813.533.732.393.271.836.810.475.63.37.144.036.596.145.962.990.543.946
Table 2. The k values obtained by the IA method.
Table 2. The k values obtained by the IA method.
LocationParameterTrialskdch
m/d
kavg,
m/d
T1T2T3T4T5T6T7T8T9T10T11T12T13T14T15T16T17T18
1234567891011121314151617181920
IA1h0, cm50.272.7------------------------------------
k, m/d0.1450.094--------------------------------0.1080.120
IA2h0, cm43.944.648.145.450.2------------------------------
k, m/d1.311.120.830.420.36--------------------------0.3920.808
IA3h0, cm27.829.829.929.927.429.919.429.929.927.614.228.929.829.728.729.929.929.9----
k, m/d8.583.323.923.363.352.324.152.651.788.4122.72.372.761.561.641.471.123.372.3714.382
IA4h0, cm17.12131.927.128.627.127.230.627.531.131.331.731.931.330.731.831.9------
k, m/d12.244.604.813.124.103.763.912.8713.474.872.212.851.602.401.301.680.89--1.6434.156
IA5h0, cm2636.450.548.850.137.847.939.438.6----------------------
k, m/d1.921.541.591.340.701.060.891.010.61------------------0.9481.184
IA6h0, cm75.462.868.573.374.870.172.669.1------------------------
k, m/d3.442.022.041.371.601.431.191.36--------------------1.4661.805
Table 3. The statistical parameters values of k variable for the TI and IA test locations.
Table 3. The statistical parameters values of k variable for the TI and IA test locations.
Test TypeLocationObtained/Observed ValuesValidity Domain
Selected ValuesSimulated/Calculated
Selected Values
No. of Trialskavg, m/dsoVoNo. of Trialsks, m/dssoVsokss, m/dsssVss
TITI1142.1191.0840.51271.5940.2620.1651.5960.0960.060
TI240.2380.0990.41930.1910.0420.2180.1920.0320.167
TI3110.9230.4790.51920.7650.0790.1040.7650.0790.104
TI4142.5465.081.997120.6460.7151.1060.6460.5610.868
TI580.6000.2300.38340.4290.1520.3550.4320.1320.307
TI6193.9461.9650.498102.7470.7700.2802.7400.2570.094
IAIA120.1200.0360.302--------0.114----
IA250.8080.4190.51920.6250.2880.4600.6250.2880.460
IA3184.3825.0351.14993.2500.5790.1783.2540.2350.072
IA4174.1563.5030.843102.8480.8850.3112.8570.7410.260
IA591.1840.4400.37230.9850.0850.0860.9850.0260.026
IA681.8050.7290.40421.5160.1160.0761.5160.1160.076
Table 4. The k mean values determined by the TI and IA tests.
Table 4. The k mean values determined by the TI and IA tests.
Test TypeLocationk, [m/d]
kdcHkavgkm = (kdcHkavg)/2kss
TITI10.9962.1191.5571.596
TI20.160.2380.1990.192
TI30.3950.9230.6590.765
TI40.0492.5461.2970.646
TI50.2370.6000.4180.432
TI60.543.9462.2432.740
IAIA10.1080.1200.1140.114
IA20.3920.8080.6000.625
IA32.3714.3823.3763.254
IA41.6434.1562.8992.857
IA50.9481.1841.0660.985
IA61.4661.8051.6601.516
Table 5. The k, A, and k/A values for each type of test in ascending order of the infiltration area values.
Table 5. The k, A, and k/A values for each type of test in ascending order of the infiltration area values.
Test Typek,
(m/d)
A_area,
(cm2)
k/A,
(m/d/cm2)
k Amplitude,
(m/d)
k/A Amplitude,
(m/d/cm2)
A Amplitude,
(cm2)
TI30.76541.8540.018282.5480.0334333.268
TI40.64641.8540.01543
TI50.43257.0120.00758
TI20.19261.5840.00312
TI62.7474.9690.03655
TI11.59675.1220.02125
IA33.254507.8700.006413.3150.00632795.230
IA42.857531.0800.00538
IA50.985733.5000.00134
IA20.626807.3000.00078
IA61.5161188.8000.00128
IA10.1191303.1000.00009
DRI32.226615.7520.003624.0940.006650
DRI21.945615.7520.00316
DRI12.711615.7520.00440
DRI55.512615.7520.00895
DRI46.039615.7520.00981
MDI10.17815.9000.011190.4670.029370
MDI20.11515.9000.00723
MDI40.24115.9000.01516
MDI30.13615.9000.00855
MDI50.10515.9000.00660
MDI60.24115.9000.01516
MDI70.17815.9000.01119
MDI180.06315.9000.00396
MDI80.26215.9000.01648
MDI160.20915.9000.01314
MDI150.10515.9000.00660
MDI90.47115.9000.02962
MDI140.08415.9000.00528
MDI100.12615.9000.00792
MDI130.26215.9000.01648
MDI110.04215.9000.00264
MDI120.00415.9000.00025
MDI170.05215.9000.00327
Average1.058251.1140.00938------
TI average1.06258.7330.01703------
IA average1.560845.2750.00255------
DRI average3.687615.7520.00599------
MDI average0.16015.9000.01004------
Table 6. The correlation parameters between dependent and independent variables.
Table 6. The correlation parameters between dependent and independent variables.
DTypes of Analyzed MethodsRankArea Interval Limits,
cm2
Ai Interval Center,
cm2
fDi FrequencyDavg(A)s2Ds2D(A)s2DoηD(A)rCalculated D’avg(A)Crp
k,
m/d
(a) All1106035210.2252.2981.7170.5810.8640.4700.8450.856
26016011031.5091.279
316066041073.5062.419
4660116091020.8062.200
511601260121011.5160.796
612601310128510.1190.296
Σ=351.058--------------
(b) All without IA1104025180.1602.4451.7420.7040.8440.8310.1760.991
24060032061.0621.289
360063061553.6874.125
Σ=290.954--------------
k/A,
m/d/cm2
(c) All1106035210.010570.0000670.0000220.0000450.577−0.4450.0155840.839
26016011030.020300.013731
316066041070.005960.007560
4660116091020.001060.001675
511601260121010.001280.000784
612601310128510.000090.000871
Σ=350.00938--------------
(d) All without IA1104025180.010040.0000680.0000130.0000550.434−0.2370.0100421.000
24060032060.017030.016995
360063061550.005990.005846
Σ=290.01079--------------
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ghibus, T.; Gaitanaru, D. A Statistical Approach for the Assessment of Saturated Hydraulic Conductivity Values of Unsaturated Urban Soils Obtained by Field Infiltration Tests. Water 2024, 16, 1908. https://doi.org/10.3390/w16131908

AMA Style

Ghibus T, Gaitanaru D. A Statistical Approach for the Assessment of Saturated Hydraulic Conductivity Values of Unsaturated Urban Soils Obtained by Field Infiltration Tests. Water. 2024; 16(13):1908. https://doi.org/10.3390/w16131908

Chicago/Turabian Style

Ghibus, Traian, and Dragos Gaitanaru. 2024. "A Statistical Approach for the Assessment of Saturated Hydraulic Conductivity Values of Unsaturated Urban Soils Obtained by Field Infiltration Tests" Water 16, no. 13: 1908. https://doi.org/10.3390/w16131908

APA Style

Ghibus, T., & Gaitanaru, D. (2024). A Statistical Approach for the Assessment of Saturated Hydraulic Conductivity Values of Unsaturated Urban Soils Obtained by Field Infiltration Tests. Water, 16(13), 1908. https://doi.org/10.3390/w16131908

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop