In the eastern sector of Mexico City the sub soil consists of high contrasting sequences (lacustrine and volcanic inter bedded deposits) that favor the development of erratic fracturing in the surface causing damage to the urban infrastructure. The high-resolution geophysical prospecting are useful tools for the assessment of ground deformation and fracturing associated with land subsidence phenomena.
The GPR method allowed to evaluate the fracture propagation and deformation of vulcano-sedimentary sequences at different depths, the main electrical parameters are directly related with the gravimetric and volumetric water content and therefore with the plasticity of the near surface prospected sequences. The active seismology prospection consisted in a combination of Seismic Refraction (SR) and Multichannel Analysis of Surface Waves (MASW) for the estimation of the velocity of the mechanical compressive (P) and the shear (S) waves. The integration of both methods allowed to estimate the geomechanical parameters characterizing the studied sequence, the Poisson Ratio and the volumetric compressibility.
The obtained mechanical parameters were correlated with laboratory measured parameters such as plasticity index, density, shear strength and compressibility and, GPR and seismic profiles were correlated with the mapped fracture systems in the study area. Once calibrated, the profiles allowed to identify the lithological contact between lacustrine and volcanic sequences, their variations of thicknesses in depth and to assess the deformation area in the surface. An accurate determination of the geometry of fracturing was of the most importance for the assessment of the geological risk in the study area.
The ground-fracturing phenomenon related with land subsidence has increased recently in Mexico City due to natural and anthropogenic processes. One of the areas with the highest differential subsidence is the lacustrine plain of Iztapalapa (Carreón-Freyre, 2011). The high-resolution geophysical prospecting, such as Ground Penetrating Radar (GPR) and active seismology, are useful tools for the assessment of ground deformation. In the eastern sector of Mexico City the sub soil consists of high contrasting sequences (lacustrine and volcanic inter bedded deposits) that favor the development of erratic fracturing in the surface causing damage to the urban infrastructure.
The GPR method is based on the generation of short electromagnetic pulses
generated by a transmitter that penetrate in the subsoil, and are reflected
according with the properties of optics laws, toward the receiving antenna.
This method allowed to evaluate the fracture propagation and deformation of
vulcano-sedimentary sequences at different depths (Daniels et al., 1998). By the use of the Common
Mid Point (CMP) arrangement was possible to measure the propagation velocity
of electromagnetic waves through the contrasting layers and allowed to build
high-resolution geological sections (Grandjean et al., 2000). The energy is reflected to the surface
by the contrasting electrical properties of each lithological sequence; the
main electrical parameters are dielectric constant,
Location of the study area, east of Mexico City with spatial distribution of ground-fracturing and geophysical study.
The MASW seismic method allows the record of the surface waves generated by a mechanical source in the surface (hammer) by the analysis of the dispersion curve of the fundamental mode, and eventually the higher modes. This curve can be inverted in a vertical section of shear strength because the surface waves are directly related by a factor of 0.97 with shear wave moreover, these take usually 70 % of the total seismic energy. The active seismology prospection consisted in a combination of Seismic Refraction (SR) and Multichannel Analysis of Surface Waves (MASW) for the estimation of the velocity of the mechanical compressive (P-) and the shear (S-) waves. The variation in depth of the shear velocity of seismic waves is directly related with some mechanical parameters of the prospected sequence such as: bulk modulus, Young's modulus, shear modulus and Poisson's ratio, used for the estimation of their deformation and potential fracturing.
The study area is located in the Basin of Mexico, east of the city, within a transition zone between volcanic and lacustrine materials (Fig. 1). The volcanic deposits come from a volcanic cone called “Peñón del Marques” composed by pyroclastic rocks and lava. Lake sediments are mainly clays, with large amounts of water content retained in their structure and with high plasticity values. The geophysical profiles were acquired at the base of the volcanic structure, over a high deformed and fractured area. GPR and seismic profiles are parallel, separated by 50 m. In the same area a geotechnical borehole was drilled for the characterization of the lithological sequence. In Fig. 1 the location of the study area is shown.
Common Mid Point (CMP) acquisition technique, where “S” denotes the transmitter location and “R” denotes the receiver location (Huisman, 2003).
GPR fundamentals are based in the electromagnetic induction, where time is a
decisive variable in conjunction with the spatial coordinates
Data and results of the three CMP, on the left side the velocity,
the field data on central side, and the processing on right side.
This is a high resolution method because it is capable to get 64 samples per
meter, allowing to identify areas of deformation and fracture by reflection
profiles. The Common Mid Point (CMP) technique was used to estimate the
propagation velocity of the electromagnetic waves in different media and
allowed the estimation of other physical properties. The CMP technique
consists on separating the transmitter and receiver at different constant
distances, as shown in Fig. 2. Three CMP profiles were acquired to
characterize the spatial variation of the geological materials by the
analysis of their velocity. Figure 3 shows the field data, processing,
and interpretation of velocities from CMP. A GPR profile of 200 m length
(radargram) was acquired complement the punctual measurements. The radargram
shown in Fig. 4 was acquired with a 100 MHz antenna with 32 scans m
The radargram shows three main layers: (a) a layer of volcanic sand 2 m width; (b) a lacustrine layer corresponding to clayey materials of variable thickness; and (c) a pyroclastic deposit is the subspace corresponding to volcanic sand until 10 m depth. The radargram shows a distortion of the sequence because of the high plasticity of the clayey materials, and fracturing recorded in the lower sandy layer.
Reflection Profile of GPR.
Mechanical properties inferred from wave velocities of body.
S-wave model, acquired from MASW method and a stratigraphic correlation with the borehole.
P-wave model, acquired from the determination of the first arrivals method and a stratigraphic correlation with the borehole. SR Method.
The MASW method deals with surface waves in the low frequencies (1–30 Hz)
and uses a shallow range of investigation (approximately 20 m depth). These
superficial waves are related with S-wave by a factor of 0.97. Shear
velocity is directly linked to the material stiffness and is one of the most
critical mechanical parameters. The shear-wave velocity (
The SR method involves the analysis of the travel times of first wave arrivals through the prospected layer to surface. The interpretation of the seismic data involves resolving the number of layers of the sequence, according to the velocity variation of each layer, and the travel-time from a given refractor to the ground surface.
Distribution of mechanical parameters obtained from modeling of
the P- and S-wave (SR Method):
Velocities obtained from the CMP data.
MASW and SR profiles were acquired using 16 geophones of 12.5 Hz acquisition frequency, spaced every 2.5 m with a total length of 40 m. The S-wave velocity is determined by analyzing the Rayleigh wave and, P-wave velocity by calculating the first arrivals. Figures 5 and 6 show the S-wave model and P-wave model respectively. In Table 1 the relationship of P-wave and S-wave velocity is presented.
The rate of the two body waves (P and S) determine the spatial distribution
of the different mechanical parameters. Equation (1) describe the
compression waves and Eq. (2) the shear wave.
The parameter
The bulk modulus, rigidity, and Young modulus, can be determined either from
static or dynamic experiments, involving the travelling of seismic waves
through the soil. The bulk modulus of a soil measures the resistance to uniform
compression. It is defined as the ratio of the infinitesimal pressure
increase to the resulting decrease of the volume. The shear modulus is the coefficient of
elasticity for a shearing force. It is defined as the ratio of shear stress
to the unit displacement per sample length. Mathematically is the tangential
force exerted and the related strain change. The elastic constant, called
the Poisson's ratio (
If
The Young modulus,
If seismic data are available, the spatial distribution of mechanical
parameters of the prospected sequence can be modeled to create
pseudo-geologic sections as shown in Fig. 7. The distribution of the
obtained mechanical parameters (
The Poisson constant distribution, varying from 0.38 to 048, and S-wave
velocities, varying from 90 to 180 m s
The distribution of the mechanical parameters (Fig. 7 and Table 2) also allowed
identifying the transition zone by the increase the thickness of the
lacustrine sequence which contact correspond to the location of fracture F
The results of applying the GPR, MASW and SR methods simultaneously show the persistence of the clayey lacustrine material with variable thickness encased by volcanic materials. This material was characterized in the laboratory a present gravimetric large water contents varying from 200 to 300 %.
The contrasting mechanical properties between the two types of material may
cause slippage over lithological contacts, causing the propagation of
deformation and fracturing of the sequence, as shown in the GPR profile of
Fig. 4. This can also be correlated with the distribution of the S-wave
velocities presented in Fig. 5, where the lower values (about 100 m s
The application of seismology and electromagnetic high resolution methods allow to characterize deformation and fracturing in a subsidence area. GPR studies allowed the monitoring of the deformation zone and fracturing and to correlate lithological variations through CMP profiles. Seismology shows to be very useful to quantify the mechanical behavior of the studied sequences.
The physical and mechanical parameters of the high deformability of the lacustrine sequence showing large water content and low shear strength. Fracturing showed to be strongly related with the lithological contact between materials having contrasting mechanical parameters in heterogeneous subsidence areas such as the east zone of Mexico City.
The authors thank his colleagues at the Centre for Evaluation of Geological Risk and the Center of Geosciences for the field support and technical suggestions for improving this article. The authors acknowledge the support of the DGAPA-UNAM research project No. IN114714.