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 cm
2). 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.
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, k
ss = 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
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.
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)).