Startseite Characterization of soil permeability in the former Lake Texcoco, Mexico
Artikel Open Access

Characterization of soil permeability in the former Lake Texcoco, Mexico

  • Norma Patricia López-Acosta EMAIL logo , Alejandra Liliana Espinosa-Santiago und David Francisco Barba-Galdámez
Veröffentlicht/Copyright: 26. März 2019
Veröffentlichen auch Sie bei De Gruyter Brill

Abstract

The geotechnical subsoil conditions of the former Lake Texcoco represent a complex sequence of highly compressible lacustrine clays interbedded with layers and seams of harder and more permeable materials. Although the mechanical properties of these deposits have been extensively studied in the past, the information about their hydraulic properties is scarce. Currently, a comprehensive characterization of the hydraulic conductivity of this site has become necessary because of the construction of the New Mexico International Airport (NAIM). The present study describes a systematic evaluation of the hydraulic conductivity in the former Lake Texcoco through three different in-situ methods (well permeameter, LEFRANC and piezocone dissipation test). The measurements, taken from 155 locations, show a high spatial variability, with ranges spanning more than two orders of magnitude. The results also reveal that the estimated permeabilities vary significantly among methods. These discrepancies reflect the scale dependency of the hydraulic conductivity in the area caused by soil heterogeneities. A comparison of the presented results with previous studies demonstrates that piezocone tests provide representative results for the clayey formations, while LEFRANC tests better estimate the hydraulic conductivity of the permeable strata. Besides, CPTu tests yield more consistent values of hydraulic conductivity, with smaller dispersion than well permeameter and LEFRANC tests.

1 Introduction

The construction of the New Mexico International Airport (NAIM) has drawn renewed attention to the subsoil conditions of the former Lake Texcoco. This area overlies a series of soft clay and clayey silts, widely known for their high water content and large compressibility, interlayered by units of cemented clastic material [1]. In the past, several authors have studied the mechanical characteristics of the clayey deposits [1, 2, 3, 4, 5, 6]. However little research has been conducted to evaluate their hydraulic properties. In particular, the available information about the permeability of the deeper strata is scarce.

A comprehensive hydraulic characterization of the area is important for the design of airport and hydraulic structures under construction. Hydraulic conductivity k affects geotechnical analyses that involve water flow within the ground or earth structures, such as uplift pressure in excavations, regional subsidence and the effects of rain infiltrations on slope stability of lagoons and canals. It also exerts a strong influence on the design of pumping systems and soft ground improvement methods. Despite its great importance, the accurate evaluation of this parameter remains a complex task because of soil’s inherent heterogeneity. Many authors have documented the great spatial variability of k, finding coefficients of variation up to several orders of magnitude [7, 8, 9]. Furthermore, soil permeability can vary with the measurement scale. Some field features (e.g., fractures, cracks and pervious seams) create preferred flow paths that increase the value of k as the volume of tested material augments [10].

The methods to determine soil hydraulic properties include laboratory and in-situ tests. Although the former are economical, relatively quick and allow controlling different boundary conditions, they require an appropriate sample handling and provide just local k values. In-situ methods, instead, evaluate a more representative volume of soil under actual field conditions. However, most of them are costly, laborious and time-consuming. These drawbacks complicate the obtainment of large datasets required for studying the spatial variability of the hydraulic conductivity. As an alternative, soil permeability can be estimated indirectly from direct-push probings carried out to determine other geotechnical properties (e.g., piezocone and flat dilatometer tests). The simplicity, inexpensiveness and speed of these indirect in-situ methods enable the acquisition of a greater number of measurements. Nevertheless, most of the proposed interpretation theories are empirical and their applicability is restricted to specific soil types and field conditions.

This paper presents a systematic evaluation of soil permeability in the former Lake Texcoco through three different in-situ methods: well permeameter, LEFRANC and piezocone dissipation test. The characterization provides the permeability values of previously unexplored strata and applies univariate statistics to evaluate the influence of the various measurement techniques. Accordingly, the specific objectives of this study are: (a) to determine the hydraulic conductivity of clayey and permeable strata of the former Lake Texcoco using different in-situ methods; (b) to evaluate the ability of indirect in-situ methods to estimate accurate and reliable values of k; and (c) to quantify the scale effects in the determination of soil permeability in the study zone.

2 Methods

2.1 Geological and geotechnical context of the study area

The study area (bounded by latitude 1928’12” – 1933’36” N and longitude 9856’24” – 9901’12”W) lies at the lowest portion of the Basin of Mexico, formerly occupied by Lake Texcoco (State of Mexico, Mexico). The Basin of Mexico is a compound graben of 9600 km2 elongated in the NNE-SSW direction (about 125-km length and 75-km width). It is bounded on the north by the Pachuca Range, east by the Río Frío and the Nevadas Ranges, south by Chichinautzin Range and west by the Ajusco Volcano and Las Cruces Range [11]. The basin is an extensive plain with an average elevation of 2240 m above sea level and a subtropical climate [12]. It was originally open until 700,000 years ago when volcanic activity blocked its drainage resulting in the development of an extensive system of lakes: Zumpango, Xaltocan, Texcoco, Xochimilco and Chalco. During this period, lacustrine sediments, volcanic ash and other pyroclastic materials were deposited overlaying previous alluvial-fluvial sediments [13]. The basin remained closed until the completion of the Nochistongo cut in 1789. The construction of new drainage systems and land reclamation for urban expansion during the 20th century caused the gradual shrinking of the lakes, which today have practically disappeared [14]. Fig. 1 shows a geologic map of the Basin of Mexico.

Figure 1 Geologic map of the Basin of Mexico (adapted from [15]).
Figure 1

Geologic map of the Basin of Mexico (adapted from [15]).

The site stratigraphy corresponds to a typical Lacustrine Zone according to Mexico City Geotechnical Zoning [16]. Soil strata consist of soft clays and clayey silts interspersed with seams and layers of harder clayey silts with sands. As seen in Table 1, the soil strata are: a) Surface Crust (SC) of light brown fat clay of soft consistency with the presence of cracks; b) Upper Clay Formation (UCF), a thick layer of very soft lacustrine clay interspersed with relatively thin seams of volcanic origin and sandy silt; c) Hard Layer (HL) composed of thin layers of hard sandy silt with variable cementation; d) Lower Clay Formation (LCF), a green brown clay interspersed with gray fat clay of the same origin of the UCF; e) Deep Deposits (DD) formed by silts and sand interspersed with hard clays; f) Deep Clay Formation (DCF), a soft to very firm clay layer with similar characteristics to the others; and g) Deep Stratified Formation (DSF) composed of stratified deposits of clay, sand and silty sand.

Table 1

Stratigraphic description of the study site.

StratumIDUSCS ClassificationThickness (m)Unit weight, γ (kN/m3)
Surface CrustSCCH0 – 1.615
Upper Clay FormationUCFCH18 – 3512
Hard LayerHLML2 – 316
Lower Clay FormationLCFCH5 – 1213
Deep DepositsDD-8 – 1018
Deep Clay FormationDCFCH9 – 1115
Deep Stratified FormationDSF-16 – 1917

The thickness of the strata remains relatively constant, with exception of the UCF that increases from northeast to southwest (Table 1). Clayey formations are characterized by their extremely high water content (varying from 300 to 400 %, but in some areas of the UCF values of up to 600 % have been registered) and void ratios (between 2-4 in the DCF, 3-6 in the LCF and above 6 in the UCF), high plasticity index and extraordinary compressibility [1, 2, 4]. Anisotropy is an important factor in Texcoco’s soils due to the presence of permeable seams and their intrinsic geological formation process. Generally, clays have greater horizontal permeability kh than vertical permeability kv, with ratios of kh/kv ranging from 1.1 to 10.0 [17, 18, 19, 20, 21, 22, 23]. Field or laboratory estimates of the anisotropy of the former Lake Texcoco clays have not been obtained yet, however, considering their layering, a value of kh/kv between 2 and 5 may be assumed [4, 18, 20].

The groundwater table is close to the surface. Previous pumping wells for brine extraction have caused a significant depletion of the hydrostatic pressure influencing the stress condition at the site (Fig. 2). In-situ piezometric measurements show that the drawdown is gradual in the study area, beginning at a lower depth in the north (between 5 and 10 m) and deepening to the south (between 20 and 24 m). In addition, the subsoil of the site has a high salt content, even the presence of gas has been detected (mainly methane).

Figure 2 Initial stress conditions at the study site: a) pore pressure, b) effective stress.
Figure 2

Initial stress conditions at the study site: a) pore pressure, b) effective stress.

2.2 Previous studies of hydraulic conductivity

Most prior research at the former Texcoco Lake has focused on evaluating the hydraulic conductivity of the uppermost strata, i.e. up to the DD. In 1974, [24] proposed a mathematical model to predict the settlement that would occur in the construction of four artificial lakes for the Lake Texcoco project. The authors estimated the permeability of the UCF and LCF from consolidation tests, whereas they obtained the permeability of the HL and DD from field pumping tests performed in the area between 1967 and 1968. In 1989, [13] investigated the hydraulic behavior of the Texcoco saline aquifer system by evaluating historical data, field studies, and numerical analyses of underground flow and solute transport. For the hydrogeological model, they determined the permeabilities of the HL and DD from a series of pumping tests and the hydraulic conductivity of the UCF and LCF from standard consolidation tests performed by [1, 2]. More recently, [4] estimated the horizontal hydraulic conductivity of the UCF and LCF fromexcess pore pressure dissipation records of a piezocone test (CPTu) done in 2000, while [25] performed a series of flexible wall permeameter tests. Table 2 presents a summary of hydraulic conductivities determined in the area according to the aforementioned studies.

Table 2

Summary of permeabilities previously estimated in the former Lake Texcoco.

Hydraulic conductivity, k (m/s)
Soil stratumHerrera et al

(1974) [24]
Rudolph et al

(1989) [13]
Alanís-González

(2003)* [4]
Lucero-Rivera

(2018) [25]
SC---1.31×10−08
UCF5.44×10−095.00×10−092.13×10−091.61×10−09
HL9.26×10−058.00×10−05--
LCF1.67×10−105.00×10−094.39×10−106.98×10−10
DD8.68×10−051.00×10−04-5.58×10−10
DCF---4.46×10−09
  1. Note: SC = Surface Crust; UCF = Upper Clay Formation; HL = Hard Layer; LCF = Lower Clay Formation; DD = Deep Deposits; and DCF = Deep Clay Formation. *Value of horizontal hydraulic conductivity.

2.3 Current measurement of saturated hydraulic conductivity k of soils

Three different in-situ measuring techniques were used to estimate the hydraulic conductivity in the study area: Well Permeameter Method (USBR 7300-89), LEFRANC (NF P94-132) and Piezocone (CPTu) test. The test campaign was part of the geotechnical survey performed for the construction of the New Mexico International Airport (NAIM) and included for the first time an extended study of the lower formations (DCF and DSF). Fig. 3 provides information about the location of the 155 test completed (21 USBR, 15 LEFRANC and 119 CPTu), including a location map of the study site.

Figure 3 Overview of the study area and location of the field tests.Note: R-# = Runway; C# = CPTu test; L# = LEFRANC test; U# = USBR test.
Figure 3

Overview of the study area and location of the field tests.

Note: R-# = Runway; C# = CPTu test; L# = LEFRANC test; U# = USBR test.

2.3.1 Well Permeameter Method (USBR 7300-89)

The well permeameter method, also known as USBR test, is an in-situ test originally developed for assessing the hydraulic conductivity of semipervious and pervious soils (k ≥ 1×10−7 m/s) along canal alignments and reservoir sites [26]. In general, it consists of injecting water into an uncased well and measuring the flow rate under a constant gravity head (Fig. 4a) Data analysis depends on the location of the groundwater table and impervious soil layers. For high groundwater tables, the hydraulic conductivity k can be estimated as:

Figure 4 In-situ tests: a) Well permeameter (USBR), b) LEFRANC, c) Piezocone dissipation test (CPTu) (photos courtesy of Federal Electricity Commission and GEOSOL).
Figure 4

In-situ tests: a) Well permeameter (USBR), b) LEFRANC, c) Piezocone dissipation test (CPTu) (photos courtesy of Federal Electricity Commission and GEOSOL).

(1)k=qV2πh2ln(h/r)(h/Tu)1+12(h/Tu)2

where h is the height of water in the well, r is the radius of the well, q is the discharge rate under steady-state condition, V is the ratio of water viscosity at the temperature of the test to water viscosity at 20 C, and Tu the unsaturated distance between the water surface in the well and the groundwater table.

2.3.2 LEFRANC test

The LEFRANC test is a variable-head method developed to estimate the hydraulic conductivity of coarse and fine soils located below the groundwater table. The test consists of injecting water by gravity into a cavity of known dimensions and measuring the subsequent water level decline over time (Fig. 4b) [27]. Percolation takes place through an infiltration chamber (usually 1 m long) filled with clean gravel located at the bottom of a casing. By evaluating the mass-balance equilibrium, the hydraulic conductivity k can be calculated as:

(2)k=AiClnh1h2t2t1

where Ai is the internal cross-section of the borehole, h1 and h2 are the differences in the total head at times t1 and t2, and C is a shape factor that depends on the ratio of the injection chamber length L to its diameter D. For ratios between 1.2 and 10, [28] suggests using the ellipsoid formula:

(3)C=2πLDlnLD+LD2+1

2.3.3 Piezocone dissipation test (CPTu)

In the piezocone (CPTu) test an instrumented steel probe is driven into the ground at a constant speed (2 cm/s) to obtain a continuous record of tip resistance, sleeve friction and pore water pressure induced during cone penetration (Fig. 4c) These measurements provide valuable information for soil profiling and estimation of its mechanical and hydraulic properties [29]. As an alternative to continuous penetration, the dissipation of the interstitial water pressure at a certain depth can be studied [4]. A dissipation test involves interrupting the steady penetration of the cone and measuring over time the decay of induced excess pore water pressure until reaching the equilibrium pressure or some percentage of it. Fig. 5 presents typical dissipation records for different clayey formations of the former Lake Texcoco.

Figure 5 Pore pressure dissipation curves of the CPTu tests at different depths in the study site.
Figure 5

Pore pressure dissipation curves of the CPTu tests at different depths in the study site.

The dissipation process is mainly radial, and it is controlled by the horizontal coefficient of consolidation Ch, which in turn is a function of the soil horizontal permeability kh [30]. There are many methods to estimate the value of kh using dissipation test data [18, 31, 32, 33, 34]. The selection of a specific interpretation theory depends on the amount of information available about the equipment, stratigraphy and initial stress conditions at the site. [35] studied the applicability of four different methods in the former Lake Texcoco and concluded that both the [18] and [34] proposals provide sufficiently representative results of the soil strata. Accordingly, the theory of [18] was chosen for the present article.

The solution proposed by [18] gives a rational interpretation method based on linear uncoupled consolidation analyses (i.e., Ch does not change during consolidation). It assumes that the soil behaves as a homogeneous, isotropic, linear-elastic material. Thus, Ch can be calculated as:

(4)Ch=R2Txtx

where R is the radius of the cone shaft, tx is the measured time to reach a degree of consolidation x and Tx is its corresponding time factor.

Eq. (4) implicitly considers that the dissipation is mainly controlled by the soil properties within a cylindrical cavity of radius R, though the failure zone caused by cone penetration extends a distance at least four to six times R [36]. In linear analysis, Ch is proportional to the ratio of permeability to compressibility, therefore [18] proposed the following equation to determine kh:

(5)kh=γw2.3pvoRRCh

where p’vo is the initial vertical effective stress, γw is the unit weight of water, RR is the recompression ratio controlling dissipation around the piezocone.

It is important to note that the piezocone provides only parameters in the horizontal direction. Parameters in the vertical direction can be estimated from empirical correlations that consider the anisotropy and stratification of soils, as discussed in Section 2.1.

3 Results

The hypothesis of linear consolidation is not valid for all soils. According to [30], linear solutions are reasonably applicable to normally consolidated clays, however [18] recommends checking the validity of this assumption when dissipation tests are first conducted at a new site. For this, they suggest calculating Ch with Eq. (4) at different degrees of consolidation. Large discrepancies between the estimated values of Ch indicate that coupling effects are significant and this methodology is not applicable. Table 3 shows the horizontal coefficient of consolidation Ch and hydraulic conductivity kh calculated at different degrees of consolidation for the dissipation data of the study site presented in Fig. 5. Discussion of results is given in Section 4.

Table 3

Horizontal coefficient of consolidation Ch and hydraulic conductivity kh calculated at different degrees of consolidation.

Depth (m)Soil stratumParameterUnitDegree of consolidation (%)
3040506070
24.38UCFChm2/s1.74×10−052.19×10−052.51×10−052.88×10−053.15×10−05
khm/s5.46×10−096.88×10−097.89×10−099.07×10−099.91×10−09
29.73LCFChm2/s8.00×10−058.26×10−058.94×10−058.49×10−059.35×10−05
khm/s1.59×10−081.64×10−081.78×10−081.69×10−081.86×10−08
45.3DCFChm2/s6.12×10−066.74×10−067.91×10−068.69×10−061.02×10−05
khm/s6.83×10−107.52×10−108.84×10−109.71×10−101.13×10−09
50.06DCFChm2/s2.49×10−052.74×10−052.79×10−053.11×10−053.53×10−05
khm/s2.56×10−092.82×10−092.87×10−093.20×10−093.63×10−09
67.59DSFChm2/s1.33×10−051.74×10−052.08×10−052.41×10−052.95×10−05
khm/s8.16×10−101.07×10−091.28×10−091.48×10−091.82×10−09
  1. Note: UCF = Upper Clay Formation; LCF = Lower Clay Formation; DCF = Deep Clay Formation; and DSF = Deep Stratified Formation.

Table 4 presents summary statistics of the hydraulic conductivity k for each in-situ test. Several authors have suggested that the variability of the hydraulic conductivity can be described using the lognormal distribution [37, 38, 39, 40, 41, 42], thus Table 4 includes the geometric mean and the standard deviation of the decimal logarithms of k as a more reliable measure of central tendency and dispersion of the dataset, respectively.

Table 4

Summary of hydraulic conductivity statistics currently estimated in the former Lake Texcoco.

TestSoil stratum
UCFHLLCFDDDCFDSF
USBRn 40-----
Min (m/s)4.60×10−06-----
Max (m/s)1.14×10−04-----
mk(m/s)2.08×10−06-----
sk1.18-----
LEFRANCn612-14-13
Min (m/s)1.07×10−098.46×10−07-8.38×10−10-5.42×10−06
Max (m/s)2.55×10−083.51×10−04-6.95×10−04-1.41×10−03
mk(m/s)6.69×10−092.39×10−05-5.57×10−06-9.11×10−05
sk0.460.88-1.84-0.76
CPTu*n46-30574650
Min (m/s)4.21×10−10-5.60×10−105.84×10−104.32×10−102.43×10−10
Max (m/s)2.40×10−08-3.42×10−083.41×10−088.18×10−094.17×10−08
mk(m/s)3.88×10−09-4.84×10−095.19×10−092.36×10−092.93×10−09
sk0.36-0.410.480.320.49
  1. Note: UCF = Upper Clay Formation; HL = Hard Layer; LCF = Lower Clay Formation; DD = Deep Deposits; DCF = Deep Clay Formation; DSF = Deep Stratified Formation; n = sample size; Min = minimum; Max = maximum; mk = geometric mean of k; sk = standard deviation of log (k). *Value of horizontal hydraulic conductivity.

4 Discussion

4.1 Applicability of piezocone dissipation test (CPTu) to estimate hydraulic conductivity in the former Lake Texcoco

The values presented in Table 3 show that Ch remains approximately constant during the consolidation process for the different clayey strata. Furthermore, the maximum difference between kh estimated with Eq. (5) is less than half order of magnitude. As stated in Section 2.1, anisotropy is an important factor of Texcoco’s clays, therefore a direct comparison of the estimated hydraulic conductivity with CPTu and previous studies [4, 13, 24, 25] is not possible. Assuming a ratio kh/kv=3 [4, 18, 20], the geometric means of kv of the UCF, LCF, DCF and DSF are 1.29×10−9, 1.61×10−9, 7.86×10−10 and 9.75×10−10 m/s, respectively. These values are similar to those reported in Table 4 and show that the selected interpretation method provide representative results for the clayey formations in the former Lake Texcoco. This confirms the applicability of linear uncoupled consolidation theory for the analysis of dissipation curves and suggests that indirect in-situ methods, such as CPTu, constitute a valid alternative to traditional field tests for determining k values in the study area. The application of such indirect techniques can help to reduce costs and time, optimizing the geotechnical site characterization process. It should be noted that the values of the estimated kv are practically in the same order of magnitude of kh. However, a comprehensive study of the in-situ anisotropy of Texcoco clayey formations is still necessary to provide reliable estimates of the vertical hydraulic conductivity.

4.2 Influence of the measurement technique

The results presented in Table 4 also suggest that the estimated hydraulic conductivity vary significantly among methods. The different locations of the tests hinder a point-to-point comparison of the measurement techniques. However, an analysis of the global summary statistics shows that both the USBR and the LEFRANC tests provide consistently higher permeabilities and exhibit greater dispersion than CPTu. In general, the values estimated with the CPTu are consistent with those reported in the literature for the clayey formations (UCF, LCF and DCF) (Table 4), but lower for the DD. Conversely, LEFRANC tests provide similar results to those reported for the more permeable and heterogeneous strata (HL and DD). Fig. 6 contrasts the distribution of the estimated k values in three strata (UCF, DD and DSF) and shows that this difference can be as many as four orders of magnitude.

Figure 6 Boxplot comparison of hydraulic conductivity obtained using CPTu, LEFRANC and USBR methods: a) UCF, b) DD, c) DSF.
Figure 6

Boxplot comparison of hydraulic conductivity obtained using CPTu, LEFRANC and USBR methods: a) UCF, b) DD, c) DSF.

Several authors have reported similar behaviors, which have been attributed to differences in flow geometry, measurement scales and inherent disturbances during testing [10, 42, 43, 44, 45]. Although an exhaustive analysis of the causes of these discrepancies is beyond the scope of the present article, scale effects are likely the most significant factors. Compared to CPTu, USBR and LEFRANC affect a larger volume of soil, which increases the influence of material heterogeneity. Fractures, cracks and pervious seams in the soil matrix create preferred pathways that allow faster flow and increase the measured value of k [10]. The presence of preferred flow paths explain the higher hydraulic conductivities estimated with the USBR test. These measurements were performed at shallow depths of the UCF (Fig. 7)which are characterized by the existence of extensive fracture networks [13]. Scale effects may also clarify the discrepancies obtained in the hydraulic conductivity of the DD and the DSF. Pumping tests used by [13, 24] have a large area of influence, slightly greater than the evaluated with the LEFRANC test. By contrast, CPTu and flexible well permeameter tests evaluate a considerably lower volume of soil, which limits the effects of the strata heterogeneities. Consequently, the permeabilities reported by the former authors are almost four orders of magnitude higher.

Figure 7 Profiles of the hydraulic conductivity variability for each stratum.
Figure 7

Profiles of the hydraulic conductivity variability for each stratum.

4.3 Spatial variability of the hydraulic conductivity

Fig. 7 evidences the variability of the hydraulic conductivity across the study area. The k values span more than two orders of magnitude regardless of the stratum or the measurement

technique. This variability is typical of natural and compacted soils, which can have a coefficient of variation CV up to 500 % [7, 8, 9]. Despite their great dispersion, the data also suggest that all clayey formations (UCF, LCF, DCF and DSF) have similar hydraulic characteristics. According to CPTu results, the horizontal permeabilities of these strata range from 10−10 to 10−8 m/s with a standard deviation between 0.3 and 0.5 (Fig. 8).

Figure 8 Boxplot comparison of hydraulic conductivity estimated with CPTu test for different clayey formations.
Figure 8

Boxplot comparison of hydraulic conductivity estimated with CPTu test for different clayey formations.

The results show the limitations of traditional univariate statistics to describe the hydraulic conductivity variability in the area. In this context, geostatistical analysis represents a suitable alternative to study the spatial correlation of the sample data and obtain a better characterization of this soil property.

5 Conclusions

The hydraulic conductivity is an important parameter for different geotechnical and geohydrological analyses that involve the movement of water within the soil or earth structures. This paper described an extensive geotechnical survey performed to characterize the subsoil hydraulic conductivity k at the former Lake Texcoco. This test campaign was part of the preliminary studies for the construction of the New Mexico International Airport (NAIM) and included for the first time an extended evaluation of unexplored strata (the deepest formations DCF and DSF). The study used three different in-situ methods: well permeameter (USBR), LEFRANC and piezocone dissipation test. The assessment of the results showed that the estimated permeabilities vary significantly among methods with differences up to four orders of magnitude. Piezocone dissipation tests gave consistently lower values of k than the well permeameter (USBR) and LEFRANC methods. CPTu tests yielded more consistent values of the horizontal hydraulic conductivity of the clayey formations, with smaller dispersion than well permeameter and LEFRANC tests, and therefore, they are more reliable. The possible reasons for these discrepancies were mainly related to scale effects. In effect, compared to CPTu, USBR and LEFRANC affect a larger volume of soil, which increases the influence of material heterogeneity. Fractures, cracks and pervious seams in the soil matrix create preferred pathways that increase the measured value and dispersion of k. On the other hand, LEFRANC tests provided better estimates of the hydraulic conductivity of the permeable strata. Finally, the results also demonstrated that more work is necessary to assess the great variability of the hydraulic conductivity in the study zone, whose values span more than two orders of magnitude regardless of the stratum or the measurement technique. Further research should focus on implementing geostatistical analyses to investigate its spatial distribution and in the evaluation of in-situ anisotropy.

Acknowledgement

The authors acknowledge the Grupo Aeroportuario de la Ciudad de México for the financial support (GACM-IIUNAMAD-SRO-CONV-DCAGI-SC-13-17) to carry out this research. We also express gratitude to the companies INGEUM, Federal Electricity Commission (CFE) and GEOSOL, who carried out the in-situ tests of this study.

References

[1] Marsal, R.J., Mazari, M., El subsuelo de la ciudad de México [The subsoil of Mexico City]. Universidad Nacional Autónoma de México, México D.F, 1959 (in Spanish).Suche in Google Scholar

[2] Marsal, R.J., Graue, R., El subsuelo del lago de Texcoco [The subsoil of Lake Texcoco]. In: Carrillo, N. (Ed.), El hundimiento de la Ciudad de México y Proyecto Texcoco. Volumen Nabor Carrillo, México D.F., 1969, 167–202 (in Spanish).Suche in Google Scholar

[3] Romo, M.P., Auvinet, G., Ovando, E., Mendoza, M.J., Taboada, V.M., Lermo, J.Y. et al Ingeniería geotécnica para el nuevo aeropuerto internacional de la Ciudad de México: soluciones en el ex Lago de Texcoco y en Zapotlán de Juárez [Geotechnical Engineering in the New Mexico International Airport: solutions in the former Lake Texcoco and in Zapotlán de Juárez]. Informe del Instituto de Ingeniería, UNAM, a Aeropuertos y servicios Auxiliares, 2001 (in Spanish).Suche in Google Scholar

[4] Alanís-González, R., Caracterización geotécnica del ex Lago de Texcoco: Ciudad de México [Geotechnical characterization of the former Lake Texcoco: Mexico City]. Masters thesis, Universidad Nacional Autónoma de México, Mexico, 2003 (in Spanish).Suche in Google Scholar

[5] Mayoral, J.M., Romo, M.P., Osorio, L., Seismic parameters characterization at Texcoco lake, Mexico. Soil Dyn. Earthq. Eng., 2008, 28(7), 507–521.10.1016/j.soildyn.2007.08.004Suche in Google Scholar

[6] Flores-Guzmán, M., Ovando-Shelley, E., Valle-Molina, C., Small-strain dynamic characterization of clayey soil from the Texcoco Lake, Mexico. Soil Dyn. Earthq. Eng., 2014, 63, 1–7.10.1016/j.soildyn.2014.03.005Suche in Google Scholar

[7] Lumb, P., Application of Statistics in Soil Mechanics. Soil Mechanics: New Horizons, Lee, I. K., ed., London, 1974.Suche in Google Scholar

[8] Benson, C.H., Probability Distributions for Hydraulic Conductivity of Compacted Soil Liners. J. Geotech. Geoenviron., 1993, 119(3), 471-486.10.1061/(ASCE)0733-9410(1993)119:3(471)Suche in Google Scholar

[9] Baecher, G.B., Christian, J.T., Reliability and Statistics in Geotechnical Engineering. John Wiley & Sons, Chichester, 2003.Suche in Google Scholar

[10] Schulze-Makuch, D., Carlson, D.A., Cherkauer, D.S., Malik, P., Scale dependency of hydraulic conductivity in heterogeneous media. Ground Water, 1999, 37(6), 904-919.10.1111/j.1745-6584.1999.tb01190.xSuche in Google Scholar

[11] de Cserna, Z., de la Fuente-Duch, M., Palacios-Nieto, M., Triay, L., Mitre-Salazar, L.M., Mota-Palomino, R., Estructura geológica, gravimetría, sismicidad y relaciones neotectónicas regionales de la Cuenca de México [Geological structure, gravimetry, seismicity and regional neotectonic relations of the Basin of Mexico]. Universidad Nacional Autónoma de México, Instituto de Geología, Boletín 104, 1988 (in Spanish).Suche in Google Scholar

[12] Lozano-García, M.S., Ortega-Guerrero, B., Late Quaternary environmental changes of the central part of the Basin of Mexico; correlation between Texcoco and Chalco basins. Rev. Palaeobot. Palyno., 1998, 99, 77–93.10.1016/S0034-6667(97)00046-8Suche in Google Scholar

[13] Rudolph D., Herrera I., Yates R., Groundwater flow and solute transport in the industrial well fields of the Texcoco saline aquifer system near Mexico City. Geofis. Int., 1989, 28(2), 363-408.10.22201/igeof.00167169p.1989.28.2.1035Suche in Google Scholar

[14] Ovando-Shelley, E., Ossa, A., Santoyo, E., Effects of regional subsidence and earthquakes on architectural monuments in Mexico City. B. Soc. Geol. Mex., 2013, 65, 157–167.10.18268/BSGM2013v65n1a12Suche in Google Scholar

[15] Ramos-Leal, J., Noyola-Medrano, C., Tapia-Silva, F., Aquifer vulnerability and groundwater quality in mega cities: Case of the Mexico Basin. Environ. Earth Sci., 2010, 61, 1309-1320.10.1007/s12665-009-0434-5Suche in Google Scholar

[16] Gobierno de la Ciudad de México (GCDMX), Normas Técnicas Complementarias para Diseño y Construcción de Cimentaciones: Ciudad de México, México, Gaceta Oficial de la Ciudad México, 15 December 2017, 8-711 (in Spanish).Suche in Google Scholar

[17] Tavenas, F., Jean, P., Leblond, P., Leroueil, S., The permeability of natural soft clays. Part II: Permeability characteristics. Can. Geotech. J., 1983, 20(4), 645-660.10.1139/t83-073Suche in Google Scholar

[18] Baligh, M.M., Levadoux, J.N., Consolidation after undrained piezocone penetration II: interpretation. J. Geotech. Eng.-ASCE, 1986, 112(7), 727-745.10.1061/(ASCE)0733-9410(1986)112:7(727)Suche in Google Scholar

[19] Hansbo, S., Design aspects of vertical drains and lime column installations. In: Southeast Asian Geotechnical Conference, SEAGS: Bangkok, Thailand, 1987, 8.1-8.12.Suche in Google Scholar

[20] Onoue, A., Consolidation of multilayered anisotropic soils by vertical drains with well resistance. Soils Found., 1988, 28(3), 75-90.10.3208/sandf1972.28.3_75Suche in Google Scholar

[21] Bergado, D., Asakami, H., Alfaro, M., Balasubramaniam, A., Smear effects of vertical drains on soft Bangkok clay. J. Geotech. Eng.-ASCE, 1991, 117(10), 1509-1530.10.1061/(ASCE)0733-9410(1991)117:10(1509)Suche in Google Scholar

[22] Shogaki, T., Moro, H., Masaharu, M., Kaneko, M., Kogure, K., Sudho, T., Effect of sample disturbance on consolidation parameters of anisotropic clays. In: Int. Sym. on Compression and Consolidation of Clayey Soils, Hiroshima, Japan, 1995, 561-566.Suche in Google Scholar

[23] Indraratna, B., Redana, I.W., Laboratory determination of smear zone due to vertical drain installation. J. Geotech. Geoenviron., 1998, 2(124), 180-184.10.1061/(ASCE)1090-0241(1998)124:2(180)Suche in Google Scholar

[24] Herrera I., Alberro J., León J., Chen B., Análisis de asentamientos para la construcción de los lagos del Plan de Texcoco: Ciudad de México, México [Settlement analysis for the construction of Plan de Texcoco Lakes: Mexico City, Mexico]. Instituto de Ingeniería, Universidad Nacional Autónoma de México, 1974 (in Spanish).Suche in Google Scholar

[25] Lucero-Rivera, J. R., Caracterización de la permeabilidad del subsuelo del ex Lago de Texcoco mediante pruebas de laboratorio [Characterization of subsoil permeability of the former Lake Texcoco through laboratory tests]. Masters thesis, Universidad Nacional Autónoma de México, Mexico, 2018 (in Spanish).Suche in Google Scholar

[26] United States Bureau of Reclamation (USBR), Procedure for Performing Field Permeability Testing by the Well Permeameter Method (USBR 7300-89). Earth Manual, Part 1, Earth Sciences and Research Laboratory Geotechnical Research, 1998.Suche in Google Scholar

[27] Chapuis, R.P., Overdamped slug test in monitoring wells: Review of interpretation methods with mathematical, physical, and numerical analysis of storativity influence. Can. Geotech. J., 1998, 35, 697–719.10.1139/t98-041Suche in Google Scholar

[28] Association Française de Normalisation (AFNOR), Reconnaissance et essais, Essai d’eau LEFRANC [Recognition and tests, LEFRANC water test]. Norme Française NF P 94-132 Sols, 2000 (in French).Suche in Google Scholar

[29] Mayne, P.W., Synthesis 368 on Cone Penetration Testing: State-of-Practice. Transportation Research Board, NCHRP Project 20-05, Task 37-14, 2007Suche in Google Scholar

[30] Levadoux, J.N., Baligh, M.M., Consolidation after undrained piezocone penetration I: Prediction. J. Geotech. Eng.-ASCE, 1986, 112(7), 707-726.10.1061/(ASCE)0733-9410(1986)112:7(707)Suche in Google Scholar

[31] Schmertmann, J.H., Guidelines for Cone Penetration Test (Performance and Design). Federal Highway Administration (FHWA), Report No. FWHA-TS-78-209, 1978.Suche in Google Scholar

[32] Parez, L., Fauriel, R., Le piezocone-ameliorations apportees a la reconnaissance des sols [Piezocone – Improvements to soil recognition]. Rev. Franç. Géotech., 1988, 44, 13-27 (in French).10.1051/geotech/1988044013Suche in Google Scholar

[33] Robertson, P.K., Sully, J.P., Woeller, D.J., Lunne, T., Powell, J.J.M., Gillespie, D.G., Estimating coefficient of consolidation from piezocone tests. Can. Geotech. J., 1992, 29(4), 539-550.10.1139/t92-061Suche in Google Scholar

[34] Robertson, P. K., Estimating in-situ Soil Permeability from CPT y CPTu. In: Proceedings of the 2nd International Symposium on Cone Penetration Testing, Pomona, CA, 2010.Suche in Google Scholar

[35] López-Acosta, N.P., Espinosa, A.L., Estimación de la permeabilidad en materiales finos mediante pruebas de disipación de presión de poro [Estimation of permeability in fine materials through pore pressure dissipation tests]. In: Memorias de la XXVIII Reunión Nacional de Ingeniería Geotécnica, SMIG: Mérida, México, 2016 (in Spanish).Suche in Google Scholar

[36] Levadoux, J.N, Baligh, M.M., Pore pressure during cone penetration in Clays. Dept. of Civil Engineering, Massachusetts Institute of Technology, Cambridge, Research Report R80-15, 1980.Suche in Google Scholar

[37] Freeze R.A., A stochastic-conceptual analysis of one-dimensional groundwater flow in nonuniform homogeneous media. Water Res., 1975, 11(5), 725–41.10.1029/WR011i005p00725Suche in Google Scholar

[38] Hoeksema, R.J., Kitanidis, P.K., Analysis of the spatial structure of properties of the selected aquifers. Water Res., 1985, 21(4), 563–72.10.1029/WR021i004p00563Suche in Google Scholar

[39] Sitar, N., Cawfield, J.D., Der Kiureghian, A., First-order reliability approach to stochastic analysis of surface flow and contaminant transport. Water Res., 1987, 23(5), 794–804.10.1029/WR023i005p00794Suche in Google Scholar

[40] Griflths, D.V., Fenton, G.A., Probabilistic analysis of exit gradients due to steady seepage. J. Geotech. Geoenviron., 1998, 124(9), 789–97.10.1061/(ASCE)1090-0241(1998)124:9(789)Suche in Google Scholar

[41] López-Acosta, N. P., Auvinet, G., Uncertainty in analyses of one-dimensional steady-state seepage through random porous media. Probabilist. Eng. Mech., 2011, 26, 501–510.10.1016/j.probengmech.2011.04.004Suche in Google Scholar

[42] López-Acosta, N.P., Barba-Galdámez, D.F., Espinosa-Santiago, A.L., Choque-Mamani, P.I., Data on horizontal hydraulic conductivity of fine-grained soils of the former Lake Texcoco (Mexico). Data in Brief, 2018, 19, 1670–1682.10.1016/j.dib.2018.06.013Suche in Google Scholar PubMed PubMed Central

[43] Rovey, C.W., Cherkauer, D.S., Scale dependency of hydraulic conductivity measurements. Ground Water, 1995, 33(5), 769–780.10.1111/j.1745-6584.1995.tb00023.xSuche in Google Scholar

[44] Gwenzi, W., Hinz, C., Holmes, K., Philips, I.R.,Mullins, I.J., Field-scale spatial variability of saturated hydraulic conductivity on a recently constructed artificial ecosystem. Geoderma, 2011, 166, 43–56.10.1016/j.geoderma.2011.06.010Suche in Google Scholar

[45] Ronayne, M.J., Houghton, T.B., Stednick, J.D., Field characterization of hydraulic conductivity in a heterogeneous alpine glacial till. J. Hydrol., 2012, 458–459, 103–109.10.1016/j.jhydrol.2012.06.036Suche in Google Scholar

Received: 2018-09-20
Accepted: 2019-01-10
Published Online: 2019-03-26

© 2019 Norma Patricia López-Acosta et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 Public License.

Artikel in diesem Heft

  1. Regular Articles
  2. 2D Seismic Interpretation of the Meyal Area, Northern Potwar Deform Zone, Potwar Basin, Pakistan
  3. A new method of lithologic identification and distribution characteristics of fine - grained sediments: A case study in southwest of Ordos Basin, China
  4. Modified Gompertz sigmoidal model removing fine-ending of grain-size distribution
  5. Diagenesis and its influence on reservoir quality and oil-water relative permeability: A case study in the Yanchang Formation Chang 8 tight sandstone oil reservoir, Ordos Basin, China
  6. Evaluation of AHRS algorithms for Foot-Mounted Inertial-based Indoor Navigation Systems
  7. Identification and evaluation of land use vulnerability in a coal mining area under the coupled human-environment
  8. Hydrocarbon Generation Potential of Chia Gara Formation in Three Selected Wells, Northern Iraq
  9. Source Analysis of Silicon and Uranium in uranium-rich shale in the Xiuwu Basin, Southern China
  10. Lithologic heterogeneity of lacustrine shale and its geological significance for shale hydrocarbon-a case study of Zhangjiatan Shale
  11. Characterization of soil permeability in the former Lake Texcoco, Mexico
  12. Detrital zircon trace elements from the Mesozoic Jiyuan Basin, central China and its implication on tectonic transition of the Qinling Orogenic Belt
  13. Turkey OpenStreetMap Dataset - Spatial Analysis of Development and Growth Proxies
  14. Morphological Changes of the Lower Ping and Chao Phraya Rivers, North and Central Thailand: Flood and Coastal Equilibrium Analyses
  15. Landscape Transformations in Rapidly Developing Peri-urban Areas of Accra, Ghana: Results of 30 years
  16. Division of shale sequences and prediction of the favorable shale gas intervals: an example of the Lower Cambrian of Yangtze Region in Xiuwu Basin
  17. Fractal characteristics of nanopores in lacustrine shales of the Triassic Yanchang Formation, Ordos Basin, NW China
  18. Selected components of geological structures and numerical modelling of slope stability
  19. Spatial data quality and uncertainty publication patterns and trends by bibliometric analysis
  20. Application of microstructure classification for the assessment of the variability of geological-engineering and pore space properties in clay soils
  21. Shear failure modes and AE characteristics of sandstone and marble fractures
  22. Ice Age theory: a correspondence between Milutin Milanković and Vojislav Mišković
  23. Are Serbian tourists worried? The effect of psychological factors on tourists’ behavior based on the perceived risk
  24. Real-Time Map Matching: A New Algorithm Integrating Spatio-Temporal Proximity and Improved Weighted Circle
  25. Characteristics and hysteresis of saturated-unsaturated seepage of soil landslides in the Three Gorges Reservoir Area, China
  26. Petrographical and geophysical investigation of the Ecca Group between Fort Beaufort and Grahamstown, in the Eastern Cape Province, South Africa
  27. Ecological risk assessment of geohazards in Natural World Heritage Sites: an empirical analysis of Bogda, Tianshan
  28. Integrated Subsurface Temperature Modeling beneath Mt. Lawu and Mt. Muriah in The Northeast Java Basin, Indonesia
  29. Go social for your own safety! Review of social networks use on natural disasters – case studies from worldwide
  30. Forestry Aridity Index in Vojvodina, North Serbia
  31. Natural Disasters vs Hotel Industry Resilience: An Exploratory Study among Hotel Managers from Europe
  32. Using Monarch Butterfly Optimization to Solve the Emergency Vehicle Routing Problem with Relief Materials in Sudden Disasters
  33. Potential influence of meteorological variables on forest fire risk in Serbia during the period 2000-2017
  34. Controlling factors on the geochemistry of Al-Shuaiba and Al-Mejarma coastal lagoons, Red Sea, Saudi Arabia
  35. The Influence of Kaolinite - Illite toward mechanical properties of Claystone
  36. Two critical books in the history of loess investigation: ‘Charakteristik der Felsarten’ by Karl Caesar von Leonhard and ‘Principles of Geology’ by Charles Lyell
  37. The Mechanism and Control Technology of Strong Strata Behavior in Extra-Thick Coal Seam Mining Influenced by Overlying Coal Pillar
  38. Shared Aerial Drone Videos — Prospects and Problems for Volunteered Geographic Information Research
  39. Stable isotopes of C and H in methane fermentation of agriculture substrates at different temperature conditions
  40. Prediction of Compression and Swelling Index Parameters of Quaternary Sediments from Index Tests at Mersin District
  41. Detection of old scattered windthrow using low cost resources. The case of Storm Xynthia in the Vosges Mountains, 28 February 2010
  42. Remediation of Copper and Zinc from wastewater by modified clay in Asir region southwest of Saudi Arabia
  43. Sedimentary facies of Paleogene lacustrine dolomicrite and implications for petroleum reservoirs in the southern Qianjiang Depression, China
  44. Correlation between ore particle flow pattern and velocity field through multiple drawpoints under the influence of a flexible barrier
  45. Atmospheric refractivity estimation from AIS signal power using the quantum-behaved particle swarm optimization algorithm
  46. A geophysical and hydro physico-chemical study of the contaminant impact of a solid waste landfill (swl) in King Williams’ Town, Eastern Cape, South Africa
  47. Landscape characterization using photographs from crowdsourced platforms: content analysis of social media photographs
  48. A Study on Transient Electromagnetic Interpretation Method Based on the Seismic Wave Impedance Inversion Model
  49. Stratigraphy of Architectural Elements of a Buried Monogenetic Volcanic System
  50. Variable secondary porosity modeling of carbonate rocks based on μ-CT images
  51. Traditional versus modern settlement on torrential alluvial fans considering the danger of debris flows: a case study of the Upper Sava Valley (NW Slovenia)
  52. The Influence of Gangue Particle size and Gangue Feeding Rate on Safety and Service Life of the Suspended Buffer’s Spring
  53. Research on the Transition Section Length of the Mixed Workface Using Gangue Backfilling Method and Caving Method
  54. Rainfall erosivity and extreme precipitation in the Pannonian basin
  55. Structure of the Sediment and Crust in the Northeast North China Craton from Improved Sequential H-k Stacking Method
  56. Planning Activities Improvements Responding Local Interests Change through Participatory Approach
  57. GIS-based landslide susceptibility mapping using bivariate statistical methods in North-western Tunisia
  58. Uncertainty based multi-step seismic analysis for near-surface imaging
  59. Deformation monitoring and prediction for residential areas in the Panji mining area based on an InSAR time series analysis and the GM-SVR model
  60. Statistical and expert-based landslide susceptibility modeling on a national scale applied to North Macedonia
  61. Natural hazards and their impact on rural settlements in NE Romania – A cartographical approach
  62. Rock fracture initiation and propagation by mechanical and hydraulic impact
  63. Influence of Rapid Transit on Accessibility Pattern and Economic Linkage at Urban Agglomeration Scale in China
  64. Near Infrared Spectroscopic Study of Trioctahedral Chlorites and Its Remote Sensing Application
  65. Problems with collapsible soils: Particle types and inter-particle bonding
  66. Unification of data from various seismic catalogues to study seismic activity in the Carpathians Mountain arc
  67. Quality assessment of DEM derived from topographic maps for geomorphometric purposes
  68. Remote Sensing Monitoring of Soil Moisture in the Daliuta Coal Mine Based on SPOT 5/6 and Worldview-2
  69. Utilizing Maximum Entropy Spectral Analysis (MESA) to identify Milankovitch cycles in Lower Member of Miocene Zhujiang Formation in north slope of Baiyun Sag, Pearl River Mouth Basin, South China Sea
  70. Stability Analysis of a Slurry Trench in Cohesive-Frictional Soils
  71. Integrating Landsat 7 and 8 data to improve basalt formation classification: A case study at Buon Ma Thuot region, Central Highland, Vietnam
  72. Assessment of the hydrocarbon potentiality of the Late Jurassic formations of NW Iraq: A case study based on TOC and Rock-Eval pyrolysis in selected oil-wells
  73. Rare earth element geochemistry of sediments from the southern Okinawa Trough since 3 ka: Implications for river-sea processes and sediment source
  74. Effect of gas adsorption-induced pore radius and effective stress on shale gas permeability in slip flow: New Insights
  75. Development of the Narva-Jõesuu beach, mineral composition of beach deposits and destruction of the pier, southeastern coast of the Gulf of Finland
  76. Selecting fracturing interval for the exploitation of tight oil reservoirs from logs: a case study
  77. A comprehensive scheme for lithological mapping using Sentinel-2A and ASTER GDEM in weathered and vegetated coastal zone, Southern China
  78. Sedimentary model of K-Successions Sandstones in H21 Area of Huizhou Depression, Pearl River Mouth Basin, South China Sea
  79. A non-uniform dip slip formula to calculate the coseismic deformation: Case study of Tohoku Mw9.0 Earthquake
  80. Decision trees in environmental justice research — a case study on the floods of 2001 and 2010 in Hungary
  81. The Impacts of Climate Change on Maximum Daily Discharge in the Payab Jamash Watershed, Iran
  82. Mass tourism in protected areas – underestimated threat? Polish National Parks case study
  83. Decadal variations of total organic carbon production in the inner-shelf of the South China Sea and East China Sea
  84. Hydrogeothermal potentials of Rogozna mountain and possibility of their valorization
  85. Postglacial talus slope development imaged by the ERT method: comparison of slopes from SW Spitsbergen, Norway and Tatra Mountains, Poland
  86. Seismotectonics of Malatya Fault, Eastern Turkey
  87. Investigating of soil features and landslide risk in Western-Atakent (İstanbul) using resistivity, MASW, Microtremor and boreholes methods
  88. Assessment of Aquifer Vulnerability Using Integrated Geophysical Approach in Weathered Terrains of South China
  89. An integrated analysis of mineralogical and microstructural characteristics and petrophysical properties of carbonate rocks in the lower Indus Basin, Pakistan
  90. Applicability of Hydrological Models for Flash Flood Simulation in Small Catchments of Hilly Area in China
  91. Heterogeneity analysis of shale reservoir based on multi-stage pumping data
Heruntergeladen am 11.9.2025 von https://www.degruyterbrill.com/document/doi/10.1515/geo-2019-0010/html
Button zum nach oben scrollen