Subsidence Induced Faulting Hazard Zonation Using Persistent Scatterer Interferometry and Horizontal Gradient Mapping in Mexican Urban Areas

Subsidence and faulting have affected Mexico city for more than a century and the process is becoming widespread throughout larger urban areas in central Mexico. This process causes substantial damages to the urban infrastructure and housing structures and will certainly become a major factor to be considered when planning urban development, land use zoning and hazard mitigation strategies in the next decades. Subsidence is usually associated with aggressive groundwater extraction rates and a general decrease of aquifer static level that promotes soil consolidation, deformation and ultimately, surface faulting. However, local stratigraphic and structural conditions also play an important role in the development and extension of faults. In all studied cases stratigraphy of the uppermost sediment strata and the structure of the underlying volcanic rocks impose a much different subsidence pattern which is most suitable for imaging through satellite geodetic techniques. We present examples from several cities in central Mexico: a) Mexico-Chalco. Very high rates of subsidence, up to 370 mm/yr are observed within this lacustrine environment surrounded by Pliocene-Quaternary volcanic structures. b) Aguascalientes where rates up to 90 mm/yr in the past decade are observed, is controlled by a stair stepped N-S trending graben that induces nucleation of faults along the edges of contrasting sediment package thicknesses. c) Morelia presents subsidence rates as high as 80 mm/yr. Differential deformation is observed across major basin-bounding E-W trending faults and with higher subsidence rates on their hanging walls, where the thickest sequences of compressible Quaternary sediments crop out. Our subsidence and faulting study in urban areas of central Mexico is based on a horizontal gradient analysis using displacement maps from Persistent Scatterer InSAR that allows definition of areas with high vulnerability to surface faulting. Correlation of the surface subsidence pattern through satellite geodesy and surface faults show that the principal factor for defining these hazardous areas is best determined not by solely using the subsidence magnitude rates but rather by using a combined magnitude and horizontal subsidence gradient analysis. This approach is used as the basis for the generation of subsidence-induced surface faulting hazard maps for the studied urban areas.