Chaussard Et Al. 2021 (Esp)

Machine Translated by Google
ARTÍCULO DE INVESTIGACIÓN Más de un siglo de hundimiento en la Ciudad de México: sin esperanzas

10.1029/2020JB020648
Recuperación significativa de capacidad de elevación y almacenamiento
Puntos clave: E. Chaussard1 , E. Havazli2 , H. Fattahi2 , E. CabralCano3 , y D. SolanoRojas4
• La integración de 115 años de datos
de deformación revela un hundimiento de 1 2
Departamento de Ciencias de la Tierra, Universidad de Oregon, Eugene, OR, EE.UU., Laboratorio de propulsión a chorro, California
hasta 50 cm/año en la Ciudad de México 3
Instituto de Tecnología, Pasadena, CA, EE. UU., Instituto de Geofísica, Universidad Nacional Autónoma de México,
• La mayor parte del hundimiento se debe 4
México CDMX, México, División de Ingeniería en Ciencias de la Tierra, Facultad de Ingeniería, Universidad Nacional
irreversible y sus tasas están correlacionadas
linealmente con el espesor del acuitardo Autónoma de México, México CDMX, México
superior
• La migración de agua de baja calidad al
acuífero productivo prepara el escenario
Resumen Se sabe que muchas áreas en todo el mundo experimentan hundimientos del suelo debido a la extracción de
para una crisis hídrica, que influye en el
panorama socioeconómico.
agua subterránea. Tradicionalmente se supone que la extensión y las tasas de hundimiento están controladas por los
lugares y volúmenes de extracción de agua subterránea. Aquí, reevaluamos este supuesto en el área metropolitana de la
Ciudad de México combinando conjuntos de datos estratigráficos, hidrológicos, geodésicos y demográficos. La integración de
Información de soporte:
115 años de nivelación con 24 años de radar interferométrico de apertura sintética (InSAR) y 14 años de datos de GPS
La información de respaldo se puede encontrar
en la versión en línea de este artículo. revela que las tasas de subsidencia han sido prácticamente constantes en la Ciudad de México desde al menos 1950
y alcanzan los 50 cm/año. El análisis de los datos de InSAR y GPS muestra que no existe una deformación elástica
Correspondencia a: significativa, lo que demuestra que el hundimiento es casi completamente irreversible. En la Ciudad de México, no existen
E. Chaussard, relaciones directas entre las fluctuaciones del nivel del agua subterránea y las tasas de subsidencia o entre las tasas
echaussard@gmail.com de bombeo y las tasas de subsidencia. Por el contrario, se aísla una fuerte relación lineal positiva entre las tasas de
subsidencia y el espesor del acuitardo superior. A través de la integración de estos conjuntos de datos a largo plazo,
Citación: pronosticamos que tomará aproximadamente 150 años alcanzar la compactación total del acuitardo superior, lo que puede
Chaussard, E., Havazli, E., Fattahi, H., conducir a una subsidencia adicional de hasta 30 m. Ahora que la superficie potenciométrica es más profunda que la mayor
CabralCano, E. y SolanoRojas, D.
parte del acuitardo, el agua intersticial de la arcilla, rica en sales, componentes químicos y contaminantes, ahora fluye
(2021). Más de un siglo de hundimiento en la ciudad
de México: no hay esperanzas de una recuperación hacia abajo hacia el acuífero productivo, disminuyendo así la calidad del agua. Finalmente, nuestro trabajo muestra que las
significativa de la elevación y la capacidad de almacenamiento. consecuencias del hundimiento del suelo influyen en gran medida en el panorama socioeconómico en el área metropolitana de la Ciudad de
Revista de investigación geofísica: Tierra
sólida, 126, e2020JB020648. https://doi.
org/10.1029/2020JB020648
1. Introducción
Recibido el 24 de julio de 2020
La sostenibilidad del agua es un problema ambiental importante que afecta a muchas áreas metropolitanas, ya sea en forma de
Aceptado el 23 MAR 2021
riesgos de inundaciones o falta de acceso al agua potable, pero ninguna está tan afectada como la Ciudad de México. La Ciudad
de México (Ciudad de México o CDMX, antes llamada Distrito Federal) es la capital de México y la ciudad más poblada de
América del Norte. En la CDMX viven 9.2 millones de personas y más de 21 millones viven en su área metropolitana. El área
metropolitana del valle de México (Zona Metropolitana del Valle de México, ZMVM), abarca la CDMX, parte del Estado de México
y un municipio del Estado de Hidalgo. La Zona Metropolitana del Valle de México (ZMVM) es el centro económico, financiero y
político de México y alberga al 17% de la población nacional, con un promedio de 92 personas por hectárea, y aporta el 23% del
Producto Interno Bruto nacional. Producto (PIB) (Pojani & Stead, 2017). La ZMVM es una paradoja hídrica: la ciudad tiene una
larga historia de inundaciones que impulsaron un manejo hidráulico que se remonta a la época prehispánica, mientras que hoy
en día una gran cantidad del agua dulce de la ciudad necesita ser traída de fuentes a más de 100 kilómetros de distancia.
La ZMVM está ubicada en una cuenca endorreica en el lecho plano de lo que alguna vez fue un gran sistema lacustre sin salida
de drenaje natural a una elevación de la base de 2200 m sobre el nivel del mar. Esta cuenca está rodeada por varias cadenas
montañosas que forman parte del Cinturón Volcánico Mexicano con elevaciones que alcanzan más de 5,000 m. La escorrentía
de agua de lluvia y manantial de estas laderas elevadas hacia el oeste, este y sur de la cuenca se acumuló en esta depresión y
se volvió salina debido a la evaporación, formando el antiguo lago de Texcoco (Figura 1a). La ciudad prehispánica de Tenochtitlán
(ahora Ciudad de México) fue fundada por los aztecas en el lago Texcoco en 1325. En la década de 1520, después de la
conquista española del Imperio Azteca, los esfuerzos por controlar las inundaciones en las llanuras lacustres condujeron a un
drenaje progresivo del lago a través de canales y diques. En 1900, el valle de México fue completamente drenado artificialmente
mediante un canal de 40 km de longitud al noreste de la cuenca y un túnel a través de las montañas. Hoy en día, las aguas
© 2021. Unión Geofísica Americana. residuales y pluviales se bombean fuera de la cuenca a través de un gran sistema de canales y túneles, mientras que el 30% del
Reservados todos los derechos. agua potable de la ciudad se importa a través de una red similar. El restante
CHAUSSARD Y AL. 1 de 18
21699356, 2021, 4, Descargado de https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2020JB020648 por Cochrane Panamá, Wiley Online Library el [17/05/2024]. Consulte los Términos y condiciones (https://onlinelibrary.wiley.com/termsandconditions) en Wiley Online Library para conocer las reglas de uso; Los artículos OA se rigen por la licencia Creative Commons aplicable.
Revista de investigación geofísica: Tierra sólida 10.1029/2020JB020648
Figura 1. (a) Ubicaciones de los sitios de datos continuos de GPS (cGPS) (triángulos) y de los sitios de nivelación histórica (rectángulos) en datos de satélite óptico Landsat de la Zona
Metropolitana del Valle de México (ZMVM). La extensión del antiguo lago de Texcoco antes del drenaje se muestra con un contorno cian (Niederberger, 1987) y la extensión del acuífero y
acuitardo se muestran en naranja y amarillo, respectivamente (Gobierno del Distrito Federal, 2004). (b) Datos disponibles de 1975 a 2020 utilizados en este trabajo. Abajo: número y
cobertura temporal de adquisiciones de datos SAR disponibles gratuitamente. Arriba: datos cGPS disponibles abiertamente (ubicaciones que se muestran en a), datos de nivelación y datos de
agua. Los datos de nivelación existen desde 1900 a partir de mediciones de hitos de la ciudad (línea discontinua negra) reemplazadas después de 1985 por estudios de redes de nivelación
(cruces negras) (Mazari & Alberro, 1991; Sistema de Aguas de la Ciudad de México, 2012). Los datos sobre el agua de trabajos publicados (las referencias se muestran en la parte inferior del
panel b) existen en forma de mapas de cabezales de agua (cruces), series temporales (líneas completas) y tasas de bombeo (líneas discontinuas), codificados por colores según los trabajos publicados.
El 70% del agua potable de la ciudad proviene de la extracción de agua subterránea de pozos ubicados a lo largo de la
cuenca con una profundidad promedio de 600 m y una profundidad máxima de 3.000 m. La sobreexplotación del sistema
acuífero está provocando graves hundimientos del terreno en la ZMVM, reconocidos desde principios del siglo XX
(Gayol, 1929). Los hundimientos exacerban el riesgo de inundaciones y dañan tanto la infraestructura de distribución
de agua como las redes de alcantarillado, provocando fugas y contaminación. Como resultado, el acceso al agua potable
es un factor importante en la segregación urbana de la riqueza y genera condiciones de vida deficientes.
Trabajos anteriores se han centrado en hidrodinámica (p. ej., Calderhead et al., 2011; Rudolph et al., 2006), hidroquímica
(p. ej., Esteller & Andreu, 2005; Esteller et al., 2012), zonificación de peligros (FernándezTorres et al., 2020), la
vulnerabilidad del agua subterránea (HernándezEspriú et al., 2014), el agotamiento del agua subterránea (p. ej., Cal
derhead et al., 2012a, 2012b; DavilaHernandez et al., 2014) y la subsidencia relacionada en la ZMVM y más allá (por
ejemplo, CabralCano et al., 2008; Castellazzi et al., 2017; Chaussard, Wdowinski, et al., 2014). Aquí, superamos los
límites de la aplicación de la geodesia espacial con una integración de datos de todo el campo en el ZMVM para mejorar
nuestra comprensión de los procesos que controlan el hundimiento de la tierra y sus impactos sociales. En particular,
investigamos la posibilidad de que exista deformación elástica y revisamos el supuesto tradicional de que las tasas de
hundimiento reflejan directamente las ubicaciones y los volúmenes de extracción de agua subterránea.
1.1. Estratigrafía del Valle de México
Aunque sólo quedan secciones muy pequeñas del antiguo lago de Texcoco, un gran sector de la ZMVM descansa sobre
los depósitos cuaternarios ricos en arcilla del lecho del lago, que constituyen el acuitardo superior, inicialmente muy
saturado con agua de poro rica en sal. Este acuitardo superior tiene hasta ~100 m de espesor (Gobierno del Distrito
Federal, 2004) (ubicación mostrada en la Figura 1a) y se superpone a una gruesa secuencia de sedimentos aluviales
piroclásticos formando un acuífero superior semiconfinado altamente productivo (OrtegaGuerrero et al., 1993). Entre
400 y 800 m de profundidad existe otra porción del acuífero, el acuífero inferior, de donde hoy en día se capta agua
(Birkle et al., 1995; MoralesCasique et al., 2018). Los acuíferos superior e inferior están separados por el acuitardo
inferior con un espesor máximo de 300 m, compuesto por depósitos lacustres del Plioceno intercalados con lavas
basálticas (Mooser et al., 1996). Dada la edad y la carga litostática experimentada por el acuitardo inferior, la
transmisividad y conductividad hidráulica de esas arcillas son significativamente menores que las del acuitardo superior.
Debajo del acuífero inferior, calizas de la formación Morelos de hasta 1,500 m de espesor constituyen el basamento
del Cretácico medio con una alta permeabilidad debido a una intensa fracturación y canales de disolución generalizados
(MoralesCasique et al., 2018).
1.2. Agua subterránea en la ZMVM
Más de 680 pozos registrados están operativos en la CDMX y cerca de 1,500 pozos registrados operan en la ZMVM, que
van desde 100 a 3,000 m de profundidad (CONAGUA, 2009). Durante décadas, las tasas de extracción anual del acuífero
(45–54 m3 /s) han excedido significativamente la recarga natural (20 m3 /s) (CONAGUA 2009), lo que llevó a la
sobreexplotación y provocó una reducción potenciométrica gradual de aproximadamente 1,0– 1,5 m por año (Lesser y
Asociados 2003; Tortajada, 2006a). Antes de la extracción extensiva de agua subterránea, persistían gradientes hidráulicos
ascendentes en todo el Valle de México (Durazo y Farvolden, 1989; Ortega y Farvolden, 1989) prevaleciendo condiciones de
flujo artesiano. En contraste, hoy en día el agua subterránea fluye hacia abajo por encima del nivel potenciométrico reducido
(HuizarAlvarez et al., 2004), lo que potencialmente conduce al ingreso de cloruro y otros componentes químicos al acuífero
productivo.
1.3. Historia del hundimiento del terreno en la ZMVM
La disminución de los niveles de agua en el acuífero reduce la presión del agua intersticial en la base del acuitardo y la
difusión de la anomalía de presión de poro negativa dentro del acuitardo conduce a su compactación, causando hundimiento
superficial (Freeze & Cherry, 1979; Santoyo et al., 2005). Si bien las perforaciones en busca de aguas subterráneas
comenzaron en la década de 1850, Gayol (1929) proporcionó la primera referencia de hundimiento en la ZMVM, seguida por
los trabajos de Cuevas (1936) y Carrillo (1948). La historia del hundimiento en el centro de la ciudad está bien documentada
y marcada por tres períodos distintos: (1) hundimiento a 8 cm/año entre 1935 y 1948; (2) subsidencia de 29 cm/año entre
1947 y 1958 debido a la deshidratación de acuitardos; y (3) hundimiento de 5 a 9 cm/año después de 1959, lo que refleja el
taponamiento de los pozos en el centro de la ciudad que drenan directamente el acuitardo superior (FigueroaVega, 1984;
OrtegaGuerrero et al., 1993) (Ver Figura 2e ). Entre 1940 y 1985 se observaron hasta 7,5 m de hundimiento del terreno y
una caída de la superficie potenciométrica de 32 m (Ortega & Farvolden, 1989). En los últimos 20 años, la resolución espacial
y temporal mejorada de la detección de hundimientos ha revelado tasas consistentes de hasta 35 a 40 cm/año al este del
centro histórico de la ciudad y tasas de hundimiento que son altamente variables espacialmente.
1.4. Procesos de hundimiento
Los procesos elásticos (reversibles) e inelásticos (permanentes) a menudo ocurren en el mismo lugar y al mismo tiempo
relacionados con la extracción de agua subterránea (Polonia e Irlanda, 1988). La deformación inelástica del suelo ocurre
cuando las cabezas de agua caen por debajo del nivel más bajo anterior, lo que lleva a un reordenamiento de los granos en
los acuitardos (Polonia e Irlanda, 1988; Wilson y Gorelick, 1996). Cuando la tensión excede la tensión previa a la consolidación,
incluso si los niveles de agua subterránea no descienden a profundidades mayores que el espesor del acuitardo, la
disminución en la presión del agua de poro intersticial conduce al colapso de las láminas de arcilla (Santoyo et al., 2005). Una
vez drenadas las arcillas pueden contraerse hasta un 30% de su volumen original (AlbertoJaime & MéndezSánchez, 2010).
Mientras que la conductividad hidráulica de los acuíferos es lo suficientemente grande como para ajustarse casi
instantáneamente a los cambios de carga, los acuitardos tienen una conductividad hidráulica mucho menor, lo que resulta en
un lento equilibrio de las fluctuaciones de carga, lo que lleva a una compactación continua incluso si los niveles de agua se
estabilizan (Chaussard, Wdowinski, et al., 2014). Si toda la ZMVM estuviera construida sobre arcillas de espesor uniforme, la
ciudad se hundiría de manera homogénea. Sin embargo, como la ciudad está sustentada por depósitos lacustres de espesor
variable intercalados con depósitos volcánicos que hunden de manera desigual, especialmente en los bordes del antiguo
sistema lacustre y alrededor de estructuras volcánicas aisladas dentro del valle, el área experimenta fuertes gradientes de
deformación. Estos gradientes están asociados con el desarrollo de franjas de intensa fractura superficial, que dañan las
estructuras urbanas, incluidos los edificios históricos y las infraestructuras urbanas críticas, como los sistemas de tuberías de
suministro de agua y alcantarillado y la infraestructura de gas natural (CabralCano et al., 2010, 2015). ; Cigna et al., 2011; Santoyo et al.
1.5. Encuestas InSAR anteriores
InSAR (Radar Interferométrico de Apertura Sintética) es el método de última generación para medir la deformación del suelo
con una alta resolución (píxel de 10 a 30 m para las misiones SAR espaciales operativas actuales) en áreas vastas (de
cientos a miles de kilómetros cuadrados) con alta precisión (hasta mm/año). El desplazamiento en la dirección de la línea de
visión (LOS) del radar de un satélite SAR se estima interfiriendo la señal de fase de los pases del mismo satélite SAR sobre
un área en diferentes momentos (Bürgmann et al., 2000). La deformación en los sistemas acuíferos de todo el mundo se ha
utilizado para (1) restringir las áreas de hundimiento, (2) definir las propiedades de los sistemas acuíferos,
Figura 2. (a–c) Mapa de velocidad vertical media del suelo de cada plataforma satelital utilizando un enfoque de series temporales basado en correlación: (a) Teledetección
europea (ERS)Envisat, (b) Satélite avanzado de observación terrestre (ALOS)1, (c) Centinela1. (d) Mapa de velocidad vertical media combinado de 1996 a 2020. (e)
Comparación entre las tasas de nivelación y hundimiento del radar de apertura sintética interferométrica (InSAR) para puntos de referencia de la ciudad (ubicaciones
mostradas en la Figura 1a) (datos de Mazari y Alberro, 1991 y Sistema de Aguas de la Ciudad de México, 2012). (f) Vista 3D de la ZMVM. El contorno de las áreas que se
hunden rápidamente (>10 cm/año) se muestra con una línea discontinua negra en la topografía y debajo se muestra el mapa de velocidad vertical combinado 3D InSAR
interpolado. El panel inferior resalta el cambio de elevación acumulativo desde 1950 considerando tasas de subsidencia lineal.
y (3) mejorar las prácticas de gestión (Amelung et al., 1999; Bell et al., 2008; Chaussard et al., 2013, 2017; Chaussard,
Bürgmann, et al., 2014; Chaussard, Wdowinski, et al., 2014 ; Miller y Shirzaei, 2015 ; En el ZMVM, Strozzi y Wegmuller
(1999) utilizaron datos SAR de 19961997 de los satélites europeos de teledetección (ERS) para aislar 30 cm/año de
hundimiento. CabralCano et al. (2008) combinaron datos SAR de los satélites ERS y Envisat con datos de GPS y
encontraron hasta 35 cm/año de hundimiento entre 1996 y 2005. LópezQuiroz et al. (2009) utilizaron datos de Envisat
de 20022007 y encontraron 38 cm/año de hundimiento y Osmanoglu et al. (2011) encontraron tasas similares al
combinar datos de Envisat con datos de GPS. Chaussard, Wdowinski et al. (2014) utilizaron datos del satélite Advanced
Land Observing Satellite (ALOS)1 entre 2007 y 2011 y también encontraron tasas de subsidencia de ~35 cm/año.
Samsonov y d'Oreye (2017) utilizaron 2008–
2012 datos ascendentes y descendentes del satélite RADARSAT2 y encontraron tasas de subsidencia vertical de 35
cm/año y movimiento horizontal de hasta 5 cm/año. Castellazzi et al. (2017) combinado Envisat, Radarsat2,
y Sentinel1 y encontraron tasas similares entre 2003 y 2016 y Du et al. (2019) utilizaron múltiples sensores para derivar
un mapa de hundimiento acumulado entre 2004 y 2017. Cigna y Tapete (2021) utilizaron datos de Sentinel1 de 20142020
y encontraron tasas de hundimiento de hasta 39 cm/año. SolanoRojas, Wdowinski, et al. (2020) consideraron varias
técnicas de posprocesamiento para derivar gradientes de deformación, mientras que Fernández et al. (2020) generaron
mapas de vulnerabilidad estructural de riesgo de subsidencia y distorsión angular. Finalmente, SolanoRojas, CabralCano,
et al. (2020), documentaron una aceleración del hundimiento en la ZMVM luego de los terremotos de Puebla y Tehuantepec
de septiembre de 2017. En este trabajo, combinamos todos los datos InSAR disponibles gratuitamente (19962020) con
datos GPS continuos (cGPS), datos de nivelación, datos de agua (ubicación de pozos, tasas de bombeo, niveles de agua
a lo largo del tiempo), datos de espesor de sedimentos (acuitardo superior y acuífero), datos de uso de la tierra y datos de
densidad de población, para reconstruir la historia de subsidencia del siglo pasado en la ZMVM y mejorar nuestra
comprensión de (1) los procesos que controlan las tasas de subsidencia de la tierra y la extensión espacial, y (2) cómo se
relaciona la subsidencia al panorama socioeconómico.
2. Métodos y datos
Para investigar si las tasas y extensión de subsidencia están controladas por las ubicaciones y volúmenes de extracción de
agua subterránea en la ZMVM, como se supone generalmente, nos basamos en la integración de datos de cuatro campos
(Figura 1b): estratigrafía (espesores de acuíferos y acuitardos), hidrología (ubicaciones de pozos). , tasas de bombeo y
niveles de agua), geodesia (InSAR, nivelación histórica y GPS) y demografía (tipo de uso de la tierra y densidad de población).
El momento de los datos de deformación y agua subterránea existente se muestra en la Figura 1b. En la Figura 1a se muestran las
ubicaciones de los sitios cGPS existentes con datos disponibles abiertamente y los sitios de nivelación histórica, así como la extensión
de acuíferos y acuitardos. Los espesores de acuitardos y acuíferos utilizados se reportan como parte de las Normas Técnicas para la
Construcción en la Ciudad de México (Gobierno de la Ciudad de México, 2019; Gobierno del Distrito Federal, 2004) obtenidos
interpolando datos de pozos profundos, estaciones de movimientos fuertes y mediciones de microtemblores. (Auvinet & Juárez, 2002;
PérezRocha et al., 1995). Los datos del nivel del agua subterránea se compilan de CarreraHernández y Gaskin (2007), Lesser y
GonzálezPosadas (2005) y Ruiz (2015), y las ubicaciones de los pozos y las tasas de bombeo se toman de CarreraHernández y
Gaskin (2007). Utilizamos datos de uso de suelo del “Centro de Investigación en Ciencias de Información Geoespacial” (CONACYT,
2020) y datos de densidad de población del “Facebook Connectivity Lab and Center for International Earth Science Information Network
CIESIN Columbia University (2016)”. Primero, presentamos el método de procesamiento InSAR y los enfoques utilizados para evaluar
la posible existencia de deformación elástica. En segundo lugar, analizamos los resultados de la deformación del suelo combinados con
el análisis de todos los conjuntos de datos mencionados anteriormente en la sección Resultados.
2.1. InSAR
Para los datos de los satélites ALOS1 y Sentinel1, utilizamos el software ISCE para producir interferogramas corregistrados
en una imagen de referencia (Fattahi et al., 2017; Rosen et al., 2015), mientras que para los datos de ERS y Para los
satélites Envisat, utilizamos el software GMTSAR (ya que ISCE no admite datos de ERS y Envisat). Para todos los
interferogramas, eliminamos las contribuciones topográficas utilizando el modelo de elevación digital de 1 segundo de arco
de Shuttle Radar Topography Mission (SRTM), realizamos múltiples miradas a los interferogramas para lograr una
resolución espacial de aproximadamente 90 m de espaciado y utilizamos el algoritmo de flujo de red de costo estadístico
( SNAPHU) para desenvolver fases (Chen y Zebker, 2002). Para rastrear la deformación del suelo entre la primera y la
última adquisición SAR de cada plataforma satelital, utilizamos el enfoque de series temporales basado en correlación
(Chaussard et al., 2015a, 2015b) a través de MintPy (Yunjun et al., 2019) para mejorar la cobertura espacial de estudios
anteriores, que se basaron en métodos de subconjuntos de líneas de base pequeñas (SBAS) o dispersores permanentes
(PS). Utilizamos el modelo atmosférico global meteorológico ERA5 para corregir el retraso troposférico (Jolivet et al., 2014)
y eliminar la señal de longitud de onda larga en forma de rampa para disminuir el retraso ionosférico, el retraso troposférico
residual y eliminar posibles señal intersísmica o postsísmica (Chaussard et al., 2016). Los errores correlacionados con las
líneas de base (errores del modelo de elevación digital [DEM]) se eliminan siguiendo el método de Fattahi y Amelung
(2013), la raíz cuadrática media (RMS) de las series temporales de fase residual para cada adquisición se calcula siguiendo
a Yunjun et al. (2019), y los errores de desenvolvimiento se identifican mediante la técnica de cierre de fase y se corrigen manualmente
Para cada plataforma satelital, el análisis de series temporales de InSAR proporciona un mapa de velocidad LOS media y
una serie temporal en cada píxel que mantiene características similares del terreno entre pasadas sucesivas (tiempo
coherencia >0,7). Los datos de ERS y Envisat se procesan juntos en series temporales (Chaussard, Bürgmann, Fattahi,
Johnson, et al., 2015), mientras que los datos de ALOS1 y Sentinel1 se procesan individualmente. Suponemos que el
componente horizontal de la deformación es insignificante y convertimos el movimiento LOS (dLOS) en desplazamiento vertical
(dv) basándose en el ángulo de incidencia medio θ de cada satélite (dv = dLOS/cosθ). Esta suposición está justificada por
trabajos previos (p. ej., Osmanoglu et al., 2011 y Samsonov & d'Oreye, 2017) que mostraron que el movimiento horizontal es
pequeño en comparación con el movimiento vertical (<5 cm/año), de acuerdo con nuestras observaciones de que las diferencias
entre las tasas verticales proyectadas de varios períodos de tiempo son <5 cm / año (Figura S1cS1e).
La buena concordancia entre los mapas de velocidad vertical ALOS1 y Sentinel1 (Figuras 2b, 2c y S1e), que no se
superponen temporalmente, sugiere que las tasas de subsidencia son en su mayoría constantes. Datos GPS de acceso abierto
de los archivos TLALOCNet (CabralCano et al., 2018), UNAVCO y NOAANGS y sus soluciones de Blewitt et al. (2018)
(Figura S2) confirman que las tasas de hundimiento han sido prácticamente constantes desde 1996. Con base en las tasas
constantes observadas por los satélites SAR individuales, las estaciones cGPS y en la concordancia entre InSAR y las tasas
de nivelación, combinamos las series temporales de los satélites individuales en una series temporales únicas para limitar la
deformación vertical media de 19962020. Para hacerlo, (1) hacemos referencia a cada serie temporal al mismo píxel, (2)
convertimos cada serie temporal a vertical usando las geometrías satelitales respectivas, (3) remuestreamos cada serie
temporal a una cuadrícula geocodificada similar y (4) aplicamos desplazamientos equivalentes en píxeles. para ajustar una
pendiente lineal más un desplazamiento a la serie temporal de cada píxel (Figura S1a). La serie temporal ALOS1 está
compensada para coincidir con la serie temporal ERSEnvisat en su período de superposición, ya que las Figuras 2a y 2b
confirman que ambas plataformas observan tasas de deformación similares. La serie temporal ERSEnvisatALOS1 está
vinculada a la serie temporal Sentinel1 calculando compensaciones por píxeles basadas en (1) la brecha temporal entre la
última fecha ALOS1 y la primera fecha Sentinel1, y (2) la tasa media de desplazamiento vertical de la serie temporal
Sentinel1, que tiene más adquisiciones que la serie temporal ERSEnvisatALOS1 (Figura S1a). El mapa de velocidad vertical
media resultante de 19962020 se muestra en la Figura 2d.
2.2. Método para aislar la posible deformación elástica
Para evaluar si una parte de la deformación InSAR observada es reversible, seguimos el método de Chaussard y Farr (2019)
y aplicamos un Análisis de Componentes Independientes (ICA) para separar fuentes independientes mezcladas en la señal de
deformación basándose en el algoritmo FastICA (Hyvärinen & Oja, 1997). ICA es un método computacional en el que se utiliza
la independencia estadística para aislar diferentes historias de desplazamiento de superficies dependientes del tiempo en
componentes independientes (CI) con distintos patrones espaciales (puntuaciones) y temporales (series temporales de valores
propios) (Chaussard et al., 2017; Gualandi et al., 2017). La cantidad de circuitos integrados a retener, que refleja la cantidad
de patrones de deformación con diferentes propiedades espaciotemporales al ignorar el ruido, se selecciona en función del
porcentaje de valores propios distintos de cero retenidos. Como el muestreo temporal de plataformas más antiguas (ERS,
Envisat, ALOS) no permite aislar posibles deformaciones a corto plazo (tiempo de revisión > 30 días), aplicamos el ICA a la
serie temporal Sentinel1 (hasta 6 días de repetición).
Los resultados se presentan para cada componente como una serie temporal de valores propios que muestra la magnitud de
la contribución de IC en cada adquisición (Figura S3b) y un mapa de puntuación que muestra los píxeles que experimentan la
serie temporal asociada (Figura S3c). El producto de la puntuación escalada en píxeles por el valor propio en un momento
dado da la contribución en términos de unidades de deformación. Nuestro ICA considera 3.724.957 muestras (píxeles) por
época y 251 épocas (adquisiciones de Sentinel1). El valor propio de un componente (IC1) corresponde al 99,3% de la suma
de todos los valores propios distintos de cero, mientras que un segundo componente retiene sólo el 0,2% de los valores
propios. Por lo tanto, solo un patrón de deformación es significativo en el ZMVM y capturado por IC1. Se obtienen resultados
similares al eliminar la tendencia lineal media de la serie temporal antes de realizar el ICA (IC1 = 97,2% de los valores propios
distintos de cero retenidos; IC2 = 0,6%) y en un Análisis de Componentes Principales (PC1 explica el 97,2% de la varianza,
PC2 explica sólo el 0,7%). Los resultados del ICA se analizan en la Sección 3.1.
3. Resultados
3.1. Subsidencia rápida y constante desde 1950
Las Figuras 2a a 2c muestran las velocidades verticales medias de InSAR procesadas individualmente con el método basado
en correlación para las cuatro plataformas satelitales (ERS, Envisat, ALOS1, Sentinel1), que abarcan el período 19962020. El
La extensión espacial y la amplitud de la deformación son similares, lo que confirma un hundimiento rápido y continuo al este del
centro histórico de la ciudad (ubicación de los sitios de nivelación, Figura 1a). El hundimiento vertical máximo observado es de
50 cm/año a 19.445 N y −99.032 E (Figura S4), justo al oeste del lago Nabor Carrillo (ubicación mostrada en las Figuras 1a y
S4). Esta área permanece coherente sólo en el mapa de velocidad media Sentinel1, que tiene el tiempo de repetición más corto
(Figura 2c). Esta tasa es mayor que las tasas máximas observadas en trabajos anteriores de 35 a 40 cm/año, ya que el área que
experimenta este hundimiento muy rápido solía estar descorrelacionada en análisis SBAS y PS anteriores. Por el contrario, una
serie temporal basada en correlación de los datos repetidos de Sentinel1 de 12 a 6 días permite recuperar la rápida deformación
en esta área que fue parcialmente drenada en 2015 como parte del esfuerzo de construcción del nuevo aeropuerto (que luego
fue abandonado) ( Figura S4). Las áreas que son comunes con trabajos anteriores muestran tasas de subsidencia similares
(Figuras 2a2d). Otra característica de la subsidencia ZMVM es la existencia de altos gradientes de deformación. Tasas de
hundimiento superiores a 10 cm/
año se observan a menos de 200 m de áreas estables que no se hunden (Figuras 2a a 2c y 2f). Los contornos de las dos áreas
que muestran tasas de subsidencia superiores a 10 cm/año se muestran con una línea discontinua negra en la Figura 2f.
(y en las Figuras 4 y 5 y en las Figuras S1, S3S6 y S8, S9). Está claro que la extensión de las áreas que se hunden está
positivamente correlacionada con la extensión del antiguo lago de Texcoco y con la extensión del acuitardo superior cuaternario
(Figura 1a). Los gradientes se explorarán más a fondo en la Sección 4.2 en el contexto de las fracturas mapeadas en el panorama
socioeconómico.
Las cuatro plataformas de satélite tienen diferentes ángulos de incidencia y, por tanto, diferente sensibilidad a la deformación
horizontal; La concordancia general entre los mapas de velocidad vertical media individuales proyectados utilizando los ángulos
de incidencia nominales de los satélites confirma además que la mayor parte de la deformación registrada por InSAR es vertical.
En las Figuras 2a a 2d, las deformaciones verticales lineales de las estaciones cGPS están codificadas por colores de manera
similar a la velocidad vertical de InSAR y muestran una concordancia completa. La serie de tiempo en la ubicación del píxel
InSAR de la estación MMX1, ubicada en el área de hundimiento, se compara con el componente de desplazamiento vertical del
GPS en la Figura S2. Al considerar todo el lapso de tiempo, la regresión lineal que mejor se ajusta a la deformación vertical del
MXX1 es una tasa de hundimiento de 27 cm/año, en comparación con una tasa de 26 cm/año en las series temporales verticales
combinadas de InSAR. El análisis de los residuos del ajuste lineal a la deformación vertical de MMX1 revela una tendencia
positiva continua a 5 mm/año, con la excepción de un período de 19 meses después del terremoto del 19 de septiembre de 2017
(Figura S2b). Al considerar la deformación MMX1 sólo antes del terremoto, se observa una pendiente similar en la deformación
vertical de InSAR y GPS. Este acuerdo entre conjuntos de datos es consistente con estudios previos a más corto plazo (Cabral
Cano et al., 2008; LópezQuiroz et al., 2009; Osmanoglu et al., 2011). Además, como lo sugirieron previamente Osmanoglu et al.
(2011), pueden existir pequeñas diferencias entre los datos de InSAR y GPS (<1 cm/año) debido a la variabilidad de las tasas de
subsidencia en distancias cortas muestreadas por InSAR, la naturaleza inherentemente discreta de las mediciones de GPS y la
diferente resolución espacial y ruido de cada técnica. Los residuos GPS de regresiones lineales se evalúan más a fondo en el
contexto de cambios a corto plazo (estacionales) en la siguiente sección.
La comparación entre los datos de nivelación y los datos de InSAR muestra una buena concordancia en las ubicaciones de los
puntos de referencia de la ciudad dentro del centro histórico de la ciudad (Figura 2e). Además, los datos de nivelación muestran
que se han observado tasas constantes de hundimiento desde la década de 1950 después de que se cerraron los pozos de agua
dentro del centro de la ciudad y se trasladaron a las entonces afueras de la ciudad. Con base en las tasas constantes observadas
por los satélites SAR individuales, las estaciones cGPS y la concordancia entre InSAR y las tasas de nivelación, asumimos
tasas de hundimiento lineales entre 1950 y 2020, reconociendo que puede existir una deformación no lineal muy localizada, y
derivamos la pérdida de elevación acumulada. (Figura 2f), que revela hasta 38,7 m de pérdida de elevación y muchas áreas con
más de 20 m de pérdida de elevación. En sistemas controlados por bombeo de agua subterránea, la deformación elástica puede
quedar enmascarada por la deformación continua a largo plazo (Chaussard & Farr, 2019). Investigamos la posible existencia de
dicha deformación elástica en el ZMVM en la siguiente sección.
3.2. Ausencia de deformación elástica “oculta” a corto plazo
La deformación a corto plazo (estacional) que se desvía de una tendencia lineal sugeriría que el sistema acuífero en
compactación tiene el potencial de recuperar la capacidad de almacenamiento bajo una gestión adecuada (Amelung et al., 1999;
Freeze & Cherry, 1979; Terzaghi, 1996). Aquí, investigamos la existencia potencial de dicha deformación a corto plazo "oculta"
dentro del rápido hundimiento a largo plazo con (1) el método ICA (Chaussard & Farr, 2019) aplicado a los datos de Sentinel1
de 20142020, y ( 2) análisis de los datos del cGPS de 20042020.
El ICA muestra que solo un componente de deformación está presente en el ZMVM, y la magnitud de IC1 representa el 99,3%
de la suma de todos los valores propios distintos de cero (Figura S3a). Se evalúa un segundo componente para garantizar
que no se pierda ninguna señal (Figura S3b y S3c), pero como IC2 solo retiene el 0,2% de los valores propios, este componente
no es significativo. La serie temporal de valores propios de IC1 revela una señal lineal con una fuerte pendiente negativa
(Figura S3b, negra), y el mapa de puntuación de IC1 (Figura S3c arriba) muestra valores positivos en las áreas de hundimiento
identificadas en los mapas de velocidad vertical media de la Figura 2. IC1 captura así el rápido y continuo hundimiento que
afecta las áreas resaltadas por los mapas de velocidad vertical media (Figura 2). Por el contrario, la serie temporal de valores
propios de IC2 revela una pequeña pendiente positiva entre 2014 y 2017 (cuando las adquisiciones de Sentinel1 se realizaban
cada 12 días), nivelándose a una pendiente insignificante después de 2016 (repetición de 6 días) y “saltos” anuales
correspondientes a cambios en las líneas de base perpendiculares. La puntuación IC2 está dominada por una señal de longitud
de onda larga hacia el oeste y el sur con valores de puntuación positivos (colores de amarillo a naranja) en la esquina suroeste (Figura S3
Esto sugiere que IC2 captura el ruido incrustado en las series temporales de InSAR y confirma que un componente de la
deformación domina en el ZMVM: el hundimiento constante, rápido y a largo plazo. Puede existir una deformación no lineal
localizada a corto plazo, como se refleja en pequeñas “manchas” en el mapa de puntuación IC2, pero es muy pequeña.
La ausencia de una señal significativa variable en el tiempo a corto plazo en el ICA muestra que no se está produciendo
ninguna deformación elástica (reversible) importante en el ZMVM. Esto contrasta con las observaciones realizadas en otros
ambientes de rápido hundimiento, como los acuíferos (prevalentemente clásticos) del Valle Central y el Valle del Silicio, en
California, EE. UU., donde en ambos casos, el ICA aisló deformaciones elásticas de corto plazo incrustadas en los acuíferos
de más largo plazo. término deformación (Chaussard, Bürgmann, et al., 2014; Chaussard & Farr, 2019; Chaussard et al., 2017).
La falta de deformación significativa a corto plazo, estacional o semestral, en la ZMVM confirma una recarga insignificante en
comparación con la extracción e implica que se está produciendo una degradación inelástica (irreversible) continua del sistema
acuíferoacuitardo subyacente a la ZMVM.
Además, investigamos el potencial de deformación a corto plazo en la señal de deformación vertical de cGPS eliminando la
tendencia a largo plazo (Figura S2). A primera vista, los residuos de desplazamiento vertical del GPS sugieren deformación
estacional (Figura S2). Sin embargo, esta estacionalidad se observa en todas las estaciones independientemente de su
ubicación en rocas blandas o duras, lo que sugiere que puede ser una señal de modo común con una huella regional (longitud
de onda larga) asociada con la carga de marea, la carga de agua superficial o la carga tectónica que son independientes de
los procesos hidrológicos locales. Para centrarnos en la deformación local inducida hidrológicamente, eliminamos dichas
señales de modo común utilizando una estación cercana en terreno estable. Utilizamos la estación TNAL en Sierra Nevada ya
que no muestra una deformación vertical clara. Se obtuvieron resultados comparables utilizando UXAL, una estación ubicada
más al este en Xalapa, Veracruz. Los residuos de GPS reducidos en modo común, que capturan toda la señal local que se
desvía de una tendencia lineal, se muestran en las Figuras S2b y S2c. No se observa una deformación estacional clara en las
estaciones ubicadas en la ZMVM tanto dentro (MMX1) como fuera (UNIP) de las áreas de hundimiento. Estas observaciones
son consistentes con el análisis InSAR y con trabajos anteriores que tampoco detectaron fluctuaciones estacionales notables
(Osmanoglu et al., 2011). La ausencia de un componente estacional en los datos de cGPS confirma que (1) la gran mayoría
de la deformación observada es irreversible y (2) la deformación no lineal, si está presente, solo está muy localizada, lo que
respalda los resultados de InSARICA.
El análisis de los residuos de un ajuste lineal en MMX1 revela un aumento de 3,5 cm/año en la tasa de subsidencia en
comparación con la tasa de fondo durante un período de 19 meses después del terremoto de Mw7,1 del 19 de septiembre de
2017, que causó terreno fuerte. movimiento en el ZMVM (Sahakian et al., 2018) (observado tanto en los residuos brutos como
en los residuos reducidos de modo común, Figura S2b). Mayor hundimiento tras este terremoto, también detectado por Solano
Rojas, CabralCano, et al. (2020), sugiere una disminución en la rigidez cercana a la superficie debido a la reordenación de los
granos que permite que los sedimentos se estiren más fácilmente. Viens et al. hicieron observaciones similares del efecto del
fuerte movimiento del suelo sobre la rigidez cercana a la superficie. (2018) después del terremoto de TohokuOki de 9,0 Mw
de 2011, donde se observó una disminución de la velocidad sísmica en el área de Tokio. Allí, se modeló una recuperación
logarítmica en varias escalas de tiempo con una reología viscoelástica para ilustrar el regreso del medio a su estado previo al
evento. En la ZMVM, el efecto del fuerte movimiento del suelo provocado por el terremoto se observa durante 19 meses (hasta
abril de 2019), pero un análisis más detallado está más allá del alcance de este trabajo y se abordará mediante trabajos
futuros que combinen análisis de ruido geodésico y sísmico.
3.3. Controles sobre tasas de hundimiento y extensión espacial
Reevaluamos el supuesto tradicional de que las tasas y extensión de hundimiento están directamente controladas por los lugares y
volúmenes de extracción de agua subterránea. Para ello utilizamos un enfoque triple. Primero, comparamos la extensión espacial de la
deformación del suelo con la evolución de los niveles de agua subterránea a través del tiempo y con las tasas de bombeo.
En segundo lugar, utilizamos un análisis multivariable que integra todos los conjuntos de datos disponibles (Sección 2) y, finalmente,
empleamos una regresión del proceso gaussiano (GPR) (un enfoque de aprendizaje automático bayesiano no paramétrico) para cuantificar
qué tan bien se pueden predecir las tasas de subsidencia mediante métodos combinados e independientes. conjuntos de datos (niveles de
agua, tasas de bombeo, espesor de acuífero, espesor de acuitardo, densidad de población y tipo de uso de la tierra).
La Figura S5 muestra diferencias escalonadas de 5 años en superficies potenciométricas en la ZMVM entre 1975 y 2000 (Figura
S5aS5e; CarreraHernández & Gaskin, 2007), para 20002003 (Figura S5f; Lesser & GonzálezPosadas , 2005). , y para 20032009
(Figura S5g, Ruiz, 2015). Si bien la reducción de los niveles de agua subterránea ocurre en las áreas que se hunden (línea
discontinua negra), la mayor disminución se observa fuera de las áreas que se hunden, lo que sugiere que las fluctuaciones del nivel
del agua y el hundimiento son procesos concurrentes pero que pueden no tener una relación directa o causal. o, alternativamente,
ese muestreo de la superficie del potenciómetro está sesgado espacialmente o es insuficiente.
La Figura S5i muestra que la gran mayoría de los pozos están ubicados fuera de las áreas de hundimiento como consecuencia de la
decisión de 1950 de migrar los pozos de producción fuera de la CDMX. Por lo tanto, no se detecta una correlación espacial obvia entre las
tasas de hundimiento y la ubicación de los pozos. HernándezEspriú et al. (2014) y Solano Rojas et al. (2015) sugirieron de manera similar
que las tasas de deformación del suelo no pueden correlacionarse únicamente con la intensidad del bombeo de agua subterránea debido
a las heterogeneidades estratigráficas existentes en el subsuelo.
Encontramos que entre 1975 y 2009 la disminución promedio del nivel del agua es aproximadamente constante en 1–1,5 m/
año, de acuerdo con CarreraHernández y Gaskin (2007), quienes informaron tasas de reducción de alrededor de 1 m/
año (y hasta 2,5 m/año en algunas zonas), OrtegaGuerrero et al. (1993) y Ortiz Zamora y Ortega Guerrero (2010), quienes informaron
tasas de reducción de 1 a 1,5 m/año, y HernándezEspriú et al. (2017), quienes informaron una caída regional del nivel del agua del orden
de 0,5 a 2 m/año. Con base en la evolución lineal de los niveles de agua subterránea, confirmada por la serie temporal de aguas
subterráneas de Solano Rojas et al. (2015), extrapolamos los datos de 19752009 para estimar la reducción de los niveles de agua entre
1950 y 2020 (Figura S5h). Reconocemos que esta interpolación temporal lineal de los niveles de agua subterránea a lo largo de la ZMVM
es una simplificación, ya que pozos específicos pueden experimentar fluctuaciones más complejas. Sin embargo, las series temporales
publicadas de niveles de agua locales confirman fluctuaciones lineales (SolanoRojas et al., 2015), y la falta de informes sistemáticos sobre
los niveles de agua nos impide desarrollar un modelo más complejo. Encontramos que en estos últimos 70 años, la disminución neta del
nivel del agua es en promedio de 107 m, mientras que el volumen total de agua subterránea extraída es de 65 a 70 km3 (Figura S5h). Esta
estimación se encuentra en el extremo superior del volumen calculado utilizando la extracción de agua reportada de pozos registrados de
65 km3. , considerando
alcanzando tasas de extracción de 5 m3 /s entre 1950 y 1960, 8,5 m3 /s entre 1960 y 1970, 16 m3 /s entre 1970 y 1980, 27 m3 /s entre 1980
y 1990, 45 m3 /s entre 1990 y 2000, 50 m3 /s entre 2000 y 2010, y 55 m3 /s desde 2010 (CONAGUA, 2009; Del Campo et al., 2014; Mejía
et al., 2007; Ruiz, 2015). La existencia de pozos de agua clandestinos en la ZMVM probablemente explica el hecho de que el volumen total
extraído calculado a partir de los cambios en el nivel del agua se encuentra en el extremo superior del volumen calculado utilizando los
volúmenes de bombeo reportados.
Para limitar los controles sobre las tasas de subsidencia y la extensión espacial, realizamos un análisis multiparamétrico y evaluamos las
posibles relaciones entre las tasas de subsidencia por píxeles y los siguientes conjuntos de datos: cambios en la cabeza del agua durante
un período de tiempo similar, tasas de bombeo, espesor del acuífero, espesor del acuitardo, densidad de población y tipo de uso de suelo
(Figura 3). Utilizamos un análisis de regresión multilineal (MLR), que se ha demostrado que permite aislar las causas del hundimiento de la
tierra (Abdollahi et al., 2019; Tien Bui et al., 2018). El análisis MLR confirma la ausencia de relaciones entre las tasas de hundimiento y los
cambios en la cabeza del agua o las tasas de bombeo (Figuras 3c y 3e). De manera similar, no se observa ninguna relación entre las tasas
de hundimiento y el tipo de uso de la tierra o la densidad de población (Figuras 3d y 3f; consulte la Figura S6 para los conjuntos de datos),
lo cual es comprensible dadas las ubicaciones de los pozos de extracción de agua. En contraste, los espesores del acuífero y del acuitardo
(que tienen atributos espaciales similares) muestran fuertes correlaciones con las tasas de subsidencia (R2 = 0,80; Figuras 3b y 3a). La
relación entre el espesor del acuitardo y las tasas de subsidencia tiene una desviación menor (0,008) y un error cuadrático medio más bajo
(0,0005) que la relación entre el espesor del acuífero y las tasas de subsidencia (0,011 y 0,0007) (Tabla S1).
Figura 3. Análisis de regresión multilineal (MLR) (datos codificados por colores según las tasas de bombeo) que evalúa las relaciones potenciales
entre las tasas de subsidencia por píxeles y: (a) espesor del acuitardo; (b) espesor del acuífero; (c) cambios en la cabeza del agua; d)
densidad de población; (e) tasas de bombeo, y (f) tipo de uso del suelo.
Esta relación lineal entre las tasas de subsidencia y el espesor de los acuitardos concuerda con las observaciones de Cabral
Cano et al. (2008) quienes correlacionaron las tasas de subsidencia con el espesor de la arcilla lacustre derivada de imágenes
sísmicas y Solano Rojas et al. (2015), quienes se basaron en un subconjunto de tasas de subsidencia. Por lo tanto, nuestro
análisis sugiere que por cada metro adicional de espesor de arcilla, la tasa de subsidencia aumenta en 3 mm/año (valor de
pendiente).
Figura 4. Representación esquemática del flujo de agua en la Zona Metropolitana del Valle de México (ZMVM) en 1975 (izquierda) y en 2020 (derecha). (Izquierda) El flujo ascendente
persistió a lo largo del fondo de la cuenca antes de las décadas de importante extracción de agua subterránea debido a la recarga e infiltración de la escorrentía de las montañas en el
borde de la cuenca, lo que condujo a condiciones artesianas. (Derecha) Por el contrario, las décadas de extracción de agua subterránea han reducido en su mayor parte el nivel
potenciométrico a un nivel más profundo que la base del acuitardo superior (Figura S9), lo que lleva a un flujo descendente de agua contaminada desde el acuitardo hacia el acuífero productor.
3.4. Pronóstico de deformación futura
Además, demostramos la validez de nuestros resultados de MLR empleando un enfoque GPR, un enfoque de
aprendizaje automático bayesiano no paramétrico. Este método se basa en el MLR al considerar relaciones lineales
y no lineales entre múltiples parámetros dependientes. Entrenamos los datos sobre los diversos parámetros (como
combinaciones e independientemente) y predecimos las tasas de subsidencia (Tabla S1; Figura S7). De acuerdo con
los resultados de MLR, se obtienen R2 comparables (0,8) con el enfoque GPR cuando se consideran los espesores
del acuitardo y del acuífero como predictores individuales, pero se obtiene una desviación menor cuando se utiliza el
espesor del acuitardo para predecir las tasas de hundimiento (Figura S7 ). Además, considerar más variables que el
espesor del acuitardo no mejora la precisión de la predicción (Tabla S1). El análisis GPR también muestra que los
datos no garantizan una relación más compleja que una regresión lineal entre las tasas de hundimiento y el espesor
del acuitardo. La complejidad agregada por una ley de ajuste de poder utilizada por Cigna y Tapete (2021) conduce a una mayor
Figura 5. (a) gradiente de la deformación vertical media (contornos blanco y cian) en comparación con las ubicaciones de las fracturas
y fisuras mapeadas (negro, datos de CarreónFreyre et al., 2017), superponiendo el radar interferométrico de apertura sintética interpolado
( InSAR) mapas de velocidad media mostrados en transparencia sobre una imagen Landsat de la Zona Metropolitana del Valle de
México (ZMVM). (b) y (c) Extensión de las áreas de hundimiento (contorno blanco) en un contexto social: (b) niveles socioeconómicos
(datos de Tortajada, 2006a), (c) viviendas con suministro de agua potable deficiente (datos de Tortajada, 2006a) .
criterio de información de Akaike que un ajuste lineal y probablemente sea el resultado del uso de datos interpolados, lo que
conduce a un muestreo sesgado.
Utilizamos las relaciones lineales observadas entre las tasas de subsidencia y el espesor del acuitardo para estimar el nivel
actual de compactación de los acuitardos y pronosticar la compactación restante (Figura S8). Encontramos que los acuitardos
están actualmente compactados en un 17% en promedio, con la mayoría de los valores oscilando entre el 15% y el 20% de
compactación (Figura S8f) y que las arcillas tardarán en promedio ~150 años en reducirse al 30% de su densidad. su volumen
original (Figura S8h; reducción máxima del volumen de arcilla original por drenaje en la ZMVM según AlbertoJaime & Méndez
Sánchez, 2010), lo que provocará un hundimiento adicional de hasta 20 m en el área de la CDMX y hasta a 30 m en la cuenca
de Chalco (Figura S8d).
4. Discusión
4.1. Cambios en el flujo de agua e impacto en la calidad del agua
En condiciones naturales, la descarga de agua subterránea en la cuenca se produjo como manantiales en el borde del Valle
de México, utilizados como fuentes de agua tanto en la época prehispánica como colonial por obras hidráulicas como el
acueducto de Chapultepec. Así, antes de la extracción extensiva de agua subterránea del acuífero semiconfinado, el flujo
ascendente persistía a través del acuífero y el acuitardo (Durazo y Farvolden, 1989; Ortega y Farvolden, 1989), prevaleciendo
condiciones artesianas en todo el fondo de la cuenca (Figura 4 a la izquierda) . ). Como resultado del flujo ascendente a través
del acuitardo, el fondo de la zona salina había migrado a un nivel por encima de la base del acuitardo (HuizarAlvarez et al.,
2004; OrtegaGuerrero et al., 1993, 1999) y permaneció aislado del agua subterránea del acuífero que se encuentra debajo.
Además, la capa de acuitardo era impermeable y, por lo tanto, actuaba para limitar la llegada de contaminantes de la superficie
al acuífero subyacente.
La deformación inelástica del suelo ocurre cuando la tensión en el sistema acuífero excede la tensión previa a la consolidación,
lo que lleva a una reorganización de los granos en los acuitardos, incluso si los niveles de agua subterránea no son más
profundos que el espesor de los acuitardos. Sin embargo, en la ZMVM, más de cien años de extracción de agua subterránea
han reducido la posición de la superficie potenciométrica a un nivel cercano o inferior al fondo del antiguo lago (dependiendo
del espesor del acuitardo) (Figuras S9 y 4). bien). Los datos hidrológicos y sedimentarios muestran que para 1990, menos del
20% de la superficie de la base del acuitardo estaba todavía a 20 m más profundo que los niveles de agua subterránea, y sólo
el 1% estaba a 40 m más profundo que los niveles de agua subterránea (Figura S9). En el año 2000, sólo el 15% de la
superficie de la base del acuitardo estaba todavía 20 m por debajo del nivel freático. Considerando la evolución lineal de los
niveles freáticos mostrada por Solano Rojas et al. (2015), para 2030, los niveles de agua habrán caído por debajo de la casi
totalidad del espesor del acuitardo.
Esperamos un retraso entre la disminución de los niveles de agua subterránea y la compactación de los acuitardos debido a
la baja conductividad hidráulica del acuitardo, lo que resulta en un equilibrio lento en relación con las fluctuaciones de las
cabezas. Por lo tanto, el hundimiento lineal continuo sugiere que la compactación de la capa de acuitardo poco profundo aún
no se ha "alcanzado" con la tensión actual y en continua evolución en el sistema. El acuitardo drena naturalmente a tasas
constantes, lo que resulta en tasas de subsidencia mayoritariamente constantes que pueden predecirse a partir del espesor
del acuitardo, en contraste con el período 19351950 durante el cual el bombeo se llevaba a cabo directamente dentro de los
acuitardos, aumentando la velocidad de drenaje. y provocando el hundimiento acelerado que se observa en los datos de
nivelación histórica. Se espera que las tasas de subsidencia observadas actualmente persistan a menos que los niveles de
agua vuelvan a subir, y la compactación residual continuaría a pesar de la posible recuperación del nivel de agua debido a la
baja conductividad hidráulica de las arcillas. Incluso si los niveles de agua subterránea aumentaran, no hay esperanza de
recuperar la mayor parte de la elevación ya perdida y la capacidad de almacenamiento del acuitardo perdida porque la
compresibilidad inelástica de las arcillas es de 1 a 3 órdenes de magnitud mayor que su compresibilidad elástica ( Pavelko,
2004; Riley, 1998).
Si bien no es necesario que los niveles de agua subterránea estén por debajo de la base del acuitardo para que las arcillas
se compacten (ya que la difusión de la anomalía de presión de poro negativa se transmite a través del sistema), la profundidad
de los niveles de agua en relación con la base del acuitardo influye en el agua. migración. El gran contraste en la conductividad
hidráulica entre el acuífero y el acuitardo debido a la despresurización del acuífero da como resultado un flujo vertical (Neuman
& Witherspoon, 1969), y la disminución de la superficie potenciométrica ha llevado a una inversión del gradiente hidráulico
(Figura 4 a la derecha). El flujo descendente ahora domina a través del acuitardo, mientras que el flujo ascendente permanece
en el acuífero hasta el nivel potenciométrico actual (Figura 4 derecha, Figura S9). Las arcillas del acuitardo son
transmitiendo así agua verticalmente a través de un “proceso de goteo” hasta el acuífero subyacente (Cardona & Hernández ,
1995). El agua de poro en el acuitardo es salina o peor, está contaminada con contaminantes superficiales o material orgánico
proveniente de aguas residuales con fugas, lo que lleva a que los acuitardos liberen una mezcla de sales y químicos al acuífero
subyacente y deterioran progresivamente y cada vez más la calidad del agua bombeada. Además, a medida que se contraen,
las arcillas se fracturan, lo que permite que los contaminantes y las aguas residuales se filtren desde la superficie (Castellazzi et
al., 2018; Rudolph et al., 1991). Por lo tanto, si bien la reducción de las cargas hidráulicas en el acuífero no controla las tasas y
la extensión de la subsidencia, es responsable de inducir la migración de agua de baja calidad del acuitardo al acuífero, lo que
eventualmente prepara el escenario para una crisis de contaminación del agua.
4.2. Contexto socioeconómico: daños a la infraestructura, niveles de pobreza y acceso al agua
Las fallas y fracturas del terreno superficial ocurren como resultado de hundimientos diferenciales, causando daños a
infraestructuras urbanas críticas (Holzer & Johnson, 1985). Calculamos los gradientes horizontales de hundimiento vertical para
aislar áreas que experimentan un gran hundimiento diferencial (tasas de hundimiento con grandes cambios de amplitud en
distancias espaciales cortas). Los cálculos de gradiente se basan en un operador de gradiente de Sobel en el que el gradiente
de un píxel es una suma ponderada de píxeles en la vecindad de 3 por 3. El gradiente resultante se procesa a través de un filtro
gaussiano para reducir el ruido y se traza como intervalos de contorno resaltando una relación de 1:1.000 entre la velocidad
vertical media y la distancia horizontal, en m/año/km (contornos máximos de 0,1 y 0,2). Luego comparamos las ubicaciones de
alto gradiente de InSAR con fisuras y fracturas mapeadas (Figura 5a, CarreónFreyre et al., 2017). La Figura 5a confirma que
las fracturas y fisuras mapeadas (líneas negras) están ubicadas en áreas de alto gradiente (contornos blancos y cian). Sin
embargo, varias áreas de gradientes altos carecen de fracturas y fisuras cartografiadas probablemente debido a: (1) el catálogo
de fracturas y fisuras está incompleto y el trabajo de campo futuro podría guiarse por este mapa de gradiente; (2) algunas áreas
tienen características urbanas limitadas, lo que impide el mapeo (es decir, en el centro del área de hundimiento norte cerca del
lago Nabor Carrillo); y (3) en algunas áreas la tensión se alivia progresivamente o se reparan las fracturas.
En la ZMVM, más de 1,1 millones de casas carecen de acceso a agua potable y más de medio millón carecen de acceso a
servicios sanitarios adecuados (Tortajada, 2006b, Figura 5c). Comparamos cualitativamente el paisaje social con áreas en
hundimiento y de alto gradiente en las Figuras 5b y 5c. Se observan altos gradientes de subsidencia, que resultan en daños a
los edificios y la infraestructura, en áreas que exhiben niveles socioeconómicos bajos y muy bajos (Figura 5b). Esta correlación
espacial entre niveles socioeconómicos y áreas con alto gradiente de deformación puede explicarse por el hecho de que los
habitantes de ingresos altos o medios que viven en áreas que experimentan una deformación diferencial alta son más aptos para
superar los desafíos asociados y, por lo tanto, son más resilientes; por el contrario, los habitantes de ingresos bajos y muy bajos
no pueden permitirse el lujo de reubicarse ya que el valor de sus bienes inmuebles está cada vez más deprimido por las fisuras
y fracturas.
Además, la mayor concentración de viviendas con acceso deficiente a agua potable se observa dentro de las áreas que se
hunden (Figura 5c). Este acceso limitado al agua puede explicarse por las siguientes observaciones. Las áreas de hundimiento
se ubican en el centro de la cuenca (Figura 5a), alejadas de las fuentes de recarga hídrica (la vertiente oriental de la Sierra de
las Cruces, la vertiente norte de la Sierra de Chichinautzin y la vertiente occidental de la Sierra Nevada, ubicaciones que se
muestran en la Figura 1a) y donde solo existe un número limitado de pozos de agua activos (Figura S5i). Sin embargo, para
llegar a esta parte de la ciudad donde no existen pozos (Figura S5i), el sistema de distribución de agua tiene que cruzar los
límites de las áreas de hundimiento y, por lo tanto, está expuesto a altos gradientes de deformación. Dichos gradientes de
deformación se correlacionan con la aparición de fallas y fracturas, lo que aumenta la aparición de daños a la infraestructura
(Figura 5a). Como resultado, el sistema de distribución de agua que llega a las áreas que se hunden rápidamente se ve afectado
de manera desproporcionada por los daños, lo que lleva a que los residentes experimenten una disminución del acceso y la
calidad del agua. Estas observaciones sugieren que los niveles de vida dentro de la ZMVM son en parte consecuencias de la
crisis del agua, y la crisis del agua refuerza aún más las disparidades existentes (FernándezTorres et al., 2020).
5. Conclusiones
La integración de los datos de nivelación de 19002010 con la deformación vertical del suelo InSAR de 19962020 y la
deformación cGPS de 20022020 muestra que las tasas de subsidencia han sido aproximadamente constantes al menos desde 1950 en el
ZMVM, cuando se tomó la última decisión de gestión profunda para el suministro local de agua subterránea a la ciudad
lugar. Nuestro enfoque de series de tiempo InSAR basado en correlación revela tasas de subsidencia de hasta 50 cm/año en
el sector noreste de la ciudad, en un área que aún no está urbanizada y donde la deformación permaneció sin detectarse
debido a la descorrelación en análisis InSAR anteriores. El análisis de componentes independientes de los datos InSAR y la
determinación del error de modo común en los datos cGPS de código abierto revelan que no existe una deformación elástica
significativa en el ZMVM. La integración de conjuntos de datos hidrológicos, estratigráficos, geodésicos y demográficos muestra
que las ubicaciones y magnitudes de la extracción de agua subterránea no pueden explicar de manera única la extensión
espacial y las tasas de hundimiento del suelo observado. En cambio, se aísla una fuerte relación lineal positiva entre las tasas
de subsidencia y el espesor del acuitardo superior. Nuestros datos muestran que por cada metro adicional de arcilla, la tasa de
subsidencia aumenta en 3 mm/año. A partir de 2020, en promedio, el acuitardo superior está compactado en un 17% y
estimamos que, en las condiciones actuales, puede llevar aproximadamente 150 años alcanzar la compactación completa, lo
que conducirá a un hundimiento adicional de hasta 30 m. Estas tasas de hundimiento persistirán a menos que los niveles del
agua vuelvan a subir a profundidades poco profundas. Incluso si se elevaran los niveles del agua, no hay esperanzas de
recuperar la gran mayoría de la elevación perdida y la capacidad de almacenamiento perdida del acuitardo.
Desde 1950, los pozos de agua están ubicados sólo fuera del acuífero confinado, pero la disminución de la presión de los poros
se transmite por todo el valle, lo que hace que los niveles de agua sean ahora más profundos que la mayor parte del espesor
Expresiones de gratitud
Los datos GPS son parte de la Red de de la arcilla. Como consecuencia, el flujo de agua en el acuitardo ahora es descendente y, dado que el agua de los poros del
Observación y Colaboración a Largo Plazo acuitardo está salada y probablemente contaminada, esta inversión del gradiente conduce a una disminución en la calidad del
Transfronteriza, Terrestre y Atmosférica
agua del acuífero. Finalmente, las áreas con altos gradientes de subsidencia coinciden con una alta densidad de fallas y
(TLALOCNet) operada por el Servicio de
Geodesia Satelital (SGS) del Instituto de Geofísica fracturas superficiales, lo que (1) aumenta la percolación de agua de baja calidad a través del sistema acuífero, (2) afecta la
Universidad Nacional Autónoma de México
infraestructura hidráulica y (3) disminuye aún más el acceso. a los recursos hídricos para la población que vive en las zonas en hundimiento
( UNAM) en colaboración con UNAVCO y con
áreas.
el apoyo de la Fundación Nacional de Ciencias
y la Administración Nacional de Aeronáutica y del
Espacio bajo el Acuerdo Cooperativo NSF Estos resultados sugieren que el hundimiento influye en gran medida en el panorama socioeconómico e incluso en la viabilidad
EAR1724794 y por NOAANGS. Los autores urbana futura de la Ciudad de México. Dado que sólo un tercio del consumo de agua en la ZMVM proviene de cuencas y
agradecen a Luis Salazar Tlaczani de SGS escorrentías adyacentes y el suministro restante se extrae de su sistema acuífero subyacente, el escenario está preparado
y a todo el personal de UNAVCO Inc. por el
mantenimiento de la estación, la adquisición para una doble crisis hídrica y de subsidencia si no se implementan acciones drásticas de gestión del agua ( como la reutilización
de datos, el soporte de TI y la distribución de aguas residuales tratadas para la recarga de aguas subterráneas y el riego de paisajes agrícolas y urbanos, por ejemplo,
de datos para estas redes. TLALOCNet y las
Chelleri et al., 2015 ;
operaciones relacionadas del SGS cuentan
con el apoyo de los proyectos 253760, 256012
y 2017015955, UNAMPrograma de Apoyo
a Proyectos de la Fundación Nacional de Declaración de disponibilidad de datos
Ciencias Número EAR1338091 y el
Consejo Nacional de Ciencia y Tecnología Los datos de ERS y Envisat utilizados están disponibles gratuitamente a través del repositorio WInSAR y los datos de ALOS1
(CONACyT). Proyectos de Investigación e
y Sentinel1 a través del repositorio de Alaska Satellite Facility. Las soluciones de series temporales de cGPS utilizadas son de
Innovación Tecnológica (PAPIIT) IN104213,
IN111509, IN1093153, IN1048183, IN107321 y libre acceso en el repositorio del Laboratorio Geodésico de Nevada. Los datos de agua utilizados están incluidos en Carrera
apoyo complementario de la UNAMInstituto de Hernández y Gaskin (2007), Lesser y GonzálezPosadas (2005) y Ruiz (2015). La ubicación de los pozos y las tasas de bombeo
Geofísica.
se incluyen en CarreraHernández y Gaskin (2007). Los datos de fallas y fisuras se incluyen en CarreónFreyre et al. (2017).
Los datos socioeconómicos utilizados en este estudio están incluidos en Tortajada (2006a).
Los espesores de acuíferos y acuitardos son de Gobierno del Distrito Federal (2004). Los datos de uso del suelo son del “Centro
Este trabajo se benefició del acceso a la
de Investigación en Ciencias de Información Geoespacial” (CONACYT, 2020) y los datos de densidad de población son del
computadora de alto rendimiento de la Universidad
de Oregon, Talapas. El primer autor (E. “Facebook Connectivity Lab and Center for International Earth Science Information Network CIESIN Columbia University
Chaussard) cuenta con el apoyo de la (2016) .” El procesamiento ISCE y GMTSAR son software de código abierto disponibles gratuitamente en WInSAR y Github.
Administración Nacional de Aeronáutica y
MintPy es un software de código abierto distribuido en Github.
del Espacio a través de la subvención
número 80NSSC19K0741. E. CabralCano
agradece el apoyo del proyecto PAPIIT
IN1048183, IN107321 y los proyectos CONACyT Referencias
256012 y 2017015955. Parte de la investigación
se llevó a cabo en el Jet Propulsion Laboratory Abdollahi, S., Pourghasemi, HR, Ghanbarian, GA y Safaeian, R. (2019). Priorización de factores efectivos en la ocurrencia de hundimiento del terreno y su mapeo de
del Instituto de Tecnología de California, en susceptibilidad utilizando un modelo SVM y sus diferentes funciones kernel. Boletín de Ingeniería Geológica y Medio Ambiente, 78, 4017–4034. https://doi.org/10.1007/
virtud de un contrato con la Administración s1006401814036
Nacional de Aeronáutica y del Espacio. La AlbertoJaime, P., & MéndezSánchez, E. (2010). Evolución de las propiedades arcillosas de la Ciudad de México afectadas por hundimientos del terreno. En D. Carreón
contribución de Emre Havazli y Heresh Fattahi a Freyre, M. Cerca, DL Galloway y JJ SilvaCorona (Eds.), Land hundimiento, peligros asociados y el papel del desarrollo de recursos naturales (339, p. 522). Wallingford,
este manuscrito es el análisis de series Reino Unido: Asociación Internacional de Ciencias Hidrológicas.
temporales InSAR del conjunto de datos Amelung, F., Galloway, DL, Bell, JW, Zebker, HA y Laczniak, RJ (1999). Sintiendo los altibajos de Las Vegas: InSAR revela el control estructural del hundimiento del suelo
Sentinel1 como se describe en la sección y la deformación del sistema acuífero. Geología, 27, 483–486. https://doi.org/10.1130/00917613(1999)
2.1. Los autores agradecen al editor asociado 027<0483:estuado>2.3.co;2
y al editor por sus comentarios constructivos. Auvinet, G. y Juárez, M. (2002). Caracterización geoestadística y simulación del subsuelo del Valle de México. En Modelado y simulación aplicados.
lación (págs. 208212), IASTED.
Bell, JW, Amelung, F., Ferretti, A., Bianchi, M. y Novali, F. (2008). El InSAR de dispersión permanente revela la respuesta estacional y a largo plazo del
sistema acuífero al bombeo de agua subterránea y la recarga artificial. Investigación de recursos hídricos, 44. https://doi.org/10.1029/2007wr006152
Birkle, P., Torres, V. y González, E. (1995). Evaluación preliminar del potencial de acuíferos profundos en la Cuenca del Valle de México. Ingeniería Hidráulica
en México, 10, 47–53.
Blewitt, G., Hammond, WC y Kreemer, C. (2018). Aprovechar la explosión de datos GPS para la ciencia interdisciplinaria. Eos, 99. https://doi.
org/10.1029/2018eo104623
Bürgmann, R., Rosen, PA y Fielding, EJ (2000). Interferometría de radar de apertura sintética para medir la topografía de la superficie terrestre y su
deformación. Revista anual de ciencias planetarias y de la Tierra, 28, 169–209. https://doi.org/10.1146/annurev.earth.28.1.169
CabralCano, E., Dixon, TH, MirallesWilhelm, F., DíazMolina, O., SánchezZamora, O., & Carande, RE (2008). Imágenes geodésicas espaciales de un
rápido hundimiento del suelo en la Ciudad de México. Boletín de la Sociedad Geológica de América, 120, 15561566. https://doi.org/10.1130/b26001.1
CabralCano, E., Osmanoglu, B., Dixon, T., Wdowinski, S., DeMets, C., Cigna, F. y DíazMolina, O. (2010). Mapas de hundimiento y peligro de fallas utilizando
PSI y redes GPS permanentes en el centro de México. En Actas del octavo simposio internacional sobre subsidencia de tierras
(págs. 1722) Querétaro, México.
CabralCano, E., PérezCampos, X., MárquezAzúa, B., Sergeeva, MA, SalazarTlaczani, L., DeMets, C., et al. (2018). TLALOCNet: Una columna vertebral
continua de GPSMet en México para la investigación sismotectónica y atmosférica. Cartas de investigación sismológica, 89, 373–381. https://doi.org/
10.1785/0220170190
CabralCano, E., SolanoRojas, D., OliverCabrera, T., Wdowinski, S., Chaussard, E., SalazarTlaczani, L., et al. (2015). Herramientas de geodesia satelital
para la evaluación del riesgo de hundimiento del terreno y fallas superficiales asociadas en el centro de México. Actas de la asociación internacional de
ciencias hidrológicas (372, págs. 255260). https://doi.org/10.5194/piahs3722552015
Calderhead, AI, Martel, R., Garfias, J., Rivera, A. y Therrien, R. (2012a). Gestión sostenible para minimizar el hundimiento de la tierra de un sistema acuífero
volcánico sobrebombeado: herramientas para el diseño de políticas. Gestión de recursos hídricos, 26, 1847–1864. https://doi.org/10.1007/
s1126901299907
Calderhead, AI, Martel, R., Garfias, J., Rivera, A. y Therrien, R. (2012b). Bombeo en seco: un creciente déficit presupuestario de agua subterránea inducido
por la urbanización, la industrialización y el cambio climático en un acuífero volcánico sobreexplotado. Ciencias ambientales de la tierra, 66, 1753–1767.
https://doi.org/10.1007/s1266501113989
Calderhead, AI, Therrien, R., Rivera, A., Martel, R. y Garfias, J. (2011). Simulando el hundimiento regional del terreno inducido por el bombeo con el uso de
InSAR y datos de campo en el Valle de Toluca, México. Avances en recursos hídricos, 34, 83–97. https://doi.org/10.1016/j.
advwaters.2010.09.017
Cardona, A. y Hernández, N. (1995). Modelo geoquímico conceptual de la evolución del agua subterránea en el Valle de México. Tecno
logia y Ciencias del agua, 10, 71–90.
CarreónFreyre, D., Cerca, M., GutiérrezCalderón, R., LópezQuiroz, P., AlcántaraDurán, C., GonzálezHernández, M., & CentenoSalas, F. (2017) . Posible
influencia de la subsidencia y fracturamiento en la Ciudad de México en las construcciones dañadas por el sismo del 19 de septiembre de 2017 (p. 13).
Geociencias. Recuperado de http://www1.cenapred.unam.mx/SUBCUENTA/22a%20SESIÓN%20ORDINARIA/VIII.%20ASUNTOS%20GENERALES/
Reporte%20CDMX%20sismo%202017.pdf
CarreraHernández, JJ y Gaskin, SJ (2007). Análisis espaciotemporal de la precipitación y temperatura diaria en la Cuenca de México.
Revista de Hidrología, 336, 231–249. https://doi.org/10.1016/j.jhidrol.2006.12.021
Carrillo, N. (1948). Influencia de los pozos artesianos en el hundimiento de la Ciudad de México. Actas de la segunda conferencia internacional sobre
mecánica de suelos e ingeniería de cimientos (7, 156159). Rotterdam, Holanda: Sociedad Internacional de Mecánica de Suelos e Ingeniería Geotécnica.
Castellazzi, P., Garfias, J., Martel, R., Brouard, C. y Rivera, A. (2017). InSAR para apoyar la urbanización sostenible sobre acuíferos en proceso de
compactación: el caso del Valle de Toluca, México. Revista internacional de geoinformación y observación de la Tierra aplicada, 63, 33–44. https://doi.
org/10.1016/j.jag.2017.06.011
Castellazzi, P., Longuevergne, L., Martel, R., Rivera, A., Brouard, C. y Chaussard, E. (2018). Mapeo cuantitativo del agotamiento de las aguas subterráneas
a escala de gestión del agua utilizando un enfoque combinado GRACE/InSAR. Teledetección del medio ambiente, 205, 408–418. https://
doi.org/10.1016/j.rse.2017.11.025
Chaussard, E., Amelung, F., Abidin, H. y Hong, S.H. (2013). Ciudades que se hunden en Indonesia: ALOS PALSAR detecta un rápido hundimiento debido
a la extracción de agua subterránea y gas. Teledetección del medio ambiente, 128, 150161. https://doi.org/10.1016/j.rse.2012.10.015
Chaussard, E., Bürgmann, R., Fattahi, H., Johnson, CW, Nadeau, R., Taira, T. y Johanson, I. (2015a). Acoplamiento intersísmico y potencial sísmico refinado
en la zona de la falla de HaywardCalaveras. Revista de investigación geofísica: Tierra sólida, 120, 8570–8590. https://doi.
org/10.1002/2015jb012230
Chaussard, E., Bürgmann, R., Fattahi, H., Nadeau, RM, Taira, T., Johnson, CW y Johanson, I. (2015b). Potencial de terremotos más grandes en el área este
de la Bahía de San Francisco debido a la conexión directa entre las fallas de Hayward y Calaveras. Cartas de investigación geofísica, 42, 2734–2741.
https://doi.org/10.1002/2015gl063575
Chaussard, E., Bürgmann, R., Shirzaei, M., Fielding, EJ y Baker, B. (2014). Previsibilidad de los cambios de carga hidráulica y caracterización del sistema
acuífero y las propiedades de las fallas a partir de la deformación del suelo derivada de InSAR. Revista de investigación geofísica: Tierra sólida, 119,
6572–6590. https://doi.org/10.1002/2014jb011266
Chaussard, E. y Farr, TG (2019). Un nuevo método para aislar la deformación elástica de la inelástica en sistemas acuíferos: Aplicación al San
Valle de Joaquín, CA. Cartas de investigación geofísica, 46, 10800–10809. https://doi.org/10.1029/2019gl084418
Chaussard, E., Johnson, CW, Fattahi, H. y Bürgmann, R. (2016). Potencial y límites de InSAR para caracterizar la deformación intersísmica
independientemente de los datos del GPS: aplicación al sistema de falla de San Andrés del sur. Geoquímica, Geofísica, Geosistemas, 17, 12141229.
https://doi.org/10.1002/2015gc006246
Chaussard, E., Milillo, P., Bürgmann, R., Perissin, D., Fielding, EJ y Baker, B. (2017). Detección remota de la deformación del suelo para monitorear las
prácticas de gestión de aguas subterráneas: aplicación al Valle de Santa Clara durante la sequía de California de 20122015. Revista de investigación
geofísica: Tierra sólida, 122, 8566–8582. https://doi.org/10.1002/2017jb014676
Chaussard, E., Wdowinski, S., CabralCano, E. y Amelung, F. (2014). Hundimientos de terreno en el centro de México detectados por ALOS InSAR
series de tiempo. Teledetección del medio ambiente, 140, 94106. https://doi.org/10.1016/j.rse.2013.08.038
Chelleri, L., Schuetze, T. y Salvati, L. (2015). Integrando la resiliencia con la sostenibilidad urbana en barrios abandonados: Desafíos y oportunidades de la
transición hacia la gestión descentralizada del agua en la Ciudad de México. Hábitat Internacional, 48, 122–130. https://doi.
org/10.1016/j.habitatint.2015.03.016
Chen, CW y Zebker, HA (2002). Desenvolvimiento de fases para grandes interferogramas SAR: segmentación estadística y red generalizada
modelos. Transacciones IEEE sobre geociencia y teledetección, 40, 17091719. https://doi.org/10.1109/tgrs.2002.802453
Cigna, F., CabralCano, E., Osmanoğlu, B., Dixon, TH y Wdowinski, S. (2011). Detección de fallas inducidas por subsidencia en áreas urbanas mexicanas mediante
interferometría de dispersión persistente y mapeo de gradiente horizontal de subsidencia. En 2011, simposio internacional de geociencia y teledetección del IEEE
(págs. 21252128), IEEE.
Cigna, F. y Tapete, D. (2021). Tasas actuales de hundimiento del terreno, peligro y riesgo de fallas superficiales en la Ciudad de México con Sentinel1 20142020
IW InSAR. Teledetección del medio ambiente, 253, 112161. https://doi.org/10.1016/j.rse.2020.112161
CONACYT. (2020). Centro de Investigación en Ciencias de Información Geoespacial. Recuperado de https://centrosconacyt.mx/
centropublico/centrogeo/
CONAGUA. (2009). Situación del Subsector Agua Potable, Alcantarillado y Saneamiento. Comisión Nacional del Agua (CONAGUA). Recuperado de http://
www.conagua.gob.mx/CONAGUA07/Publicaciones/Publicaciones/LibroAnexosYTablasSituaci%C3%B3nSAPAS.
pdf
Cuevas, J. (1936). Condiciones de fundación en la Ciudad de México. Conferencia internacional sobre mecánica de suelos e ingeniería de cimentaciones, (3, pp.
233–237).
DávilaHernández, N., Madrigal, D., Expósito, JL y Antonio, X. (2014). Análisis multitemporal de la subsidencia del terreno en el Valle de Toluca (México) mediante una
combinación de interferometría de dispersión persistente (PSI) y datos piezométricos históricos. Avances en teledetección, 3, 49–60. https://doi.org/10.4236/
ars.2014.32005
Del Campo, MM, Esteller, MV, Expósito, JL, & Hirata, R. (2014). Impactos de la urbanización en la hidrodinámica e hidroquímica de las aguas subterráneas del acuífero
del Valle de Toluca (México). Monitoreo y evaluación ambiental, 186, 2979–2999. https://doi.org/10.1007/
s1066101335953
Du, Z., Ge, L., Ng, AHM, Zhu, Q., Zhang, Q., Kuang, J. y Dong, Y. (2019). Hundimientos de larga duración en la Ciudad de México de 2004 a 2018 revelados por cinco
sensores de radar de apertura sintética. Degradación y desarrollo de la tierra, 30, 1785–1801. https://doi.org/10.1002/ldr.3347
Durazo, J. y Farvolden, RN (1989). El régimen de aguas subterráneas del Valle de México a partir de evidencia histórica y observaciones de campo. diario
nal de Hidrología, 112, 171190. https://doi.org/10.1016/00221694(89)90187x
Esteller, MV y Andreu, JM (2005). Efectos antrópicos sobre las características hidroquímicas del acuífero del Valle de Toluca (centro de México).
Revista de hidrogeología, 13, 378–390. https://doi.org/10.1007/s1004000403954
Esteller, MV, Rodríguez, R., Cardona, A. y PadillaSánchez, L. (2012). Evaluación de cambios hidroquímicos debido a la explotación intensiva de acuíferos: estudios
de caso de México. Monitoreo y evaluación ambiental, 184, 5725–5741. https://doi.org/10.1007/
s1066101123760
Laboratorio de Conectividad de Facebook y Centro para la Red Internacional de Información sobre Ciencias de la Tierra CIESIN Universidad de Columbia. (2016).
Capa de liquidación de alta resolución (HRSL). Imágenes fuente para HRSL © 2016 DigitalGlobe. Obtenido de https://ciesin.columbia.edu/data/hrsl/
Fattahi, H., Agram, P. y Simons, M. (2017). Un enfoque de diversidad espectral mejorada basado en red para el análisis de series temporales TOPS. IEEE
Transacciones sobre geociencia y teledetección, 55, 777–786. https://doi.org/10.1109/tgrs.2016.2614925
Fattahi, H. y Amelung, F. (2013). Corrección de errores DEM en series temporales InSAR. Transacciones IEEE sobre geociencia y teledetección, 51,
4249–4259. https://doi.org/10.1109/tgrs.2012.2227761
FernándezTorres, E., CabralCano, E., SolanoRojas, D., Havazli, E. y SalazarTlaczani, L. (2020). Mapas de riesgo de hundimiento del terreno y evaluación de
vulnerabilidad estructural de distorsión angular basada en InSAR: un ejemplo en la Ciudad de México. Actas de la Asociación Internacional de Ciencias Hidrológicas,
382, 583–587. https://doi.org/10.5194/piahs3825832020
FigueroaVega, GE (1984). Historia Caso No. 9.8, México, DF, México. En J. Polonia (Ed.), Guía para estudios de subsidencia de tierras debido a la extracción de agua
subterránea: estudios e informes en hidrología (40, págs. 217–232), UNESCO.
Freeze, RA y Cherry, JA (1979). Agua subterránea (pág. 604). Englewood Cliffs, Nueva Jersey: PrenticeHall.
Gayol, R. (1929). Breves apuntes relativos a las obras de saneamiento y desague de la capital de la republica y de las que, del mismo genero
necesita con urgencia. Revista Mexicana de Ingeniería y Arquitectura, 8.
Gobierno de la Ciudad de México. (2019). Decreto por el que se reforman, añaden y derogan diversas disposiciones del Reglamento de
Construcciones para el Distrito Federal. Gaceta Oficial de la Ciudad de México, 21 época, 2 de abril de 2019, No. 63 Bis (p. 12).
Gobierno del Distrito Federal. (2004). Normas técnicas complementarias del reglamento de construcciones del distrito federal. Tomo II: Normas técnicas complementarias
sobre criterios y acciones para el diseño estructural de las edificaciones, cimentaciones, viento,sismo, hidráulicas, arquitec. Gaceta Oficial del Distrito Federal, 14
época, 6 de octubre de 2004, Tomo II, No. 103Bis (p. 302). Obtenido de http://www.
paot.org.mx/centro/gaceta/Octubre2004.html
Gualandi, A., Avouac, J.P., Galetzka, J., Genrich, JF, Blewitt, G., Adhikari, LB, et al. (2017). Deformación pre y post sísmica relacionada con el terremoto de Gorkha de
Mw7,8 de 2015, Nepal. Tectonofísica, 714–715, 90–106. https://doi.org/10.1016/j.tecto.2016.06.014
HernándezEspriú, A., ArangoGalván, C., ReyesPimentel, A., MartínezSantos, P., Pita de la Paz, C., MacíasMedrano, S., et al. (2017).
Evaluación de fuentes de suministro de agua en zonas de recarga de acuíferos no gestionados: estudio de caso del Valle del Mezquital (México). Agua, 9, 4. https://
doi.org/10.3390/w9010004
HernándezEspriú, A., ReynaGutiérrez, JA, SánchezLeón, E., CabralCano, E., CarreraHernández, J., MartínezSantos, P., et al. (2014).
El modelo DRASTICSg: Una extensión del enfoque DRASTIC para mapear la vulnerabilidad del agua subterránea en acuíferos sujetos a hundimiento diferencial
del terreno, con aplicación a la Ciudad de México. Revista de hidrogeología, 22, 1469–1485. https://doi.org/10.1007/s1004001411304
Holzer, TL y Johnson, AI (1985). Hundimiento del terreno causado por la extracción de agua subterránea en áreas urbanas. GeoDiario, 11,
245–255. https://doi.org/10.1007/bf00186338
HuizarÁlvarez, R., CarrilloRivera, JJ, ÁngelesSerrano, G., Hergt, T. y Cardona, A. (2004). Respuesta química a la extracción de aguas subterráneas.
ción al sureste de la Ciudad de México. Revista de hidrogeología, 12, 436–450. https://doi.org/10.1007/s1004000403433
Hyvärinen, A. y Oja, E. (1997). Un algoritmo rápido de punto fijo para el análisis de componentes independientes. Computación neuronal, 9, 14831492. https://doi.org/
10.1162/neco.1997.9.7.1483
Jolivet, R., Agram, PS, Lin, NY, Simons, M., Doin, MP, Peltzer, G. y Li, Z. (2014). Mejora de la geodesia InSAR utilizando la atmósfera global
modelos. Revista de investigación geofísica: Tierra sólida, 119, 2324–2341. https://doi.org/10.1002/2013jb010588
Lesser, JM y GonzálezPosadas, D. (2005). El agua subterránea de la Ciudad de México . En V congreso de Aguas subterráneas. Hermosilló, hijo. Lesser y
Asociados, SA.
Lesser y Asociados. (2003). Mapa de niveles estáticos con base en información del Sistema de Aguas de la Ciudad de México durante agosto
2002 (pág. 200).
LópezQuiroz, P., Doin, M.P., Tupin, F., Briole, P., & Nicolas, J.M. (2009). Análisis de series de tiempo de la subsidencia de la Ciudad de México limitada por
Interferometría de radar. Revista de Geofísica Aplicada, 69, 115. https://doi.org/10.1016/j.jappgeo.2009.02.006
MazariHiriart, M., TapiaPalacios, MA, ZarcoArista, AE y EspinosaGarcía, AC (2019). Retos y oportunidades en lo urbano
Calidad del agua en la Ciudad de México. Fronteras en las ciencias ambientales, 7, 169. https://doi.org/10.3389/fenvs.2019.00169
Mazari, M. y Alberro, J. (1991). Hundimento de la Ciudad de México. En M. Mazari, & J. Kumate (Eds.), Los problemas de la Cuenca de
México (pág. 403). Ciudad de México: El Colegio Nacional.
Mejía, L., Murillo, S., & Reyes, J. (2007b). Plan de manejo del Valle de Toluca. Cuenca y Acuífero. Problemática y Perspectivas. Subdirección
General de Programación; Dirección Local Estado de México. CONAGUA (Comisión Nacional del Agua).
Miller, MM y Shirzaei, M. (2015). Caracterización espaciotemporal del hundimiento y levantamiento del suelo en Phoenix utilizando series temporales InSAR y
transformadas wavelet. Revista de investigación geofísica: Tierra sólida, 120, 5822–5842. https://doi.org/10.1002/2015jb012017
Mooser, F., Montiel, A. y Zúñiga, A. (1996). Nuevo mapa geológico de las cuencas de México, toluca y Puebla, estratigrafía, tectónica regional y aspectos geotérmicos
(pág. 63). Comisión Federal de Electricidad.
MoralesCasique, E., Zapata Norberto, B., & Herrera, GS (2018). Estimación de la conductividad hidráulica vertical y el hundimiento del terreno durante el flujo
unidimensional no lineal de agua subterránea y la consolidación en acuitardos altamente compresibles mediante el filtro Ensemble Kalman. American Geophysical
Union, Resúmenes de las reuniones de otoño de 2018. Obtenido de https://fallmeeting.agu.org/2018/
Motagh, M., Walter, TR, Sharifi, MA, Fielding, E., Schenk, A., Anderssohn, J. y Zschau, J. (2008). El hundimiento del terreno en Irán causado por
sobreexplotación generalizada de los embalses de agua. Cartas de investigación geofísica, 35. https://doi.org/10.1029/2008gl033814
Consejo Nacional de Investigaciones (1995). Abastecimiento de agua en la Ciudad de México: mejorando las perspectivas para la sostenibilidad. Washington, DC: Academia Nacional
Prensa. Recuperado de http://lanic.utexax.edu/la/Mexico/water/book.html
Neuman, SP y Witherspoon, PA (1969). Teoría del flujo en un sistema de dos acuíferos confinados. Investigación de recursos hídricos, 5, 803–816. https://
doi.org/10.1029/wr005i004p00803
Niederberger, C. (1987). Paléopaysages et archéologie préurbaine du bassin de México (Méxique) (p. 357). Centro de estudios mexicanos y
centroaméricaines.
Ortega, AG y Farvolden, RN (1989). Análisis computacional del flujo regional de aguas subterráneas y condiciones de contorno en la Cuenca de México.
Revista de Hidrología, 110, 271–294. https://doi.org/10.1016/00221694(89)901923
OrtegaGuerrero, A., Cherry, JA y Rudolph, DL (1993). Consolidación de acuitardos a gran escala cerca de la Ciudad de México. Agua subterránea, 31, 708–718.
OrtegaGuerrero, A., Rudolph, DL y Cherry, JA (1999). Análisis de hundimientos de suelo a largo plazo cerca de la Ciudad de México: Investigaciones de campo
y modelado predictivo. Investigación de recursos hídricos, 35, 3327–3341. https://doi.org/10.1029/1999wr900148
Ortiz Zamora, D., & Ortega Guerrero, A. (2010). Evolución de la subsidencia de terrenos a largo plazo cerca de la Ciudad de México: revisión, investigaciones de campo y
simulaciones predictivas. Investigación de recursos hídricos, 46. https://doi.org/10.1029/2008WR007398
Osmanoğlu, B., Dixon, TH, Wdowinski, S., CabralCano, E. y Jiang, Y. (2011). Se observa subsidencia en la Ciudad de México con InSAR de dispersión persistente.
Revista Internacional de Geoinformación y Observación Aplicada de la Tierra, 13, 1–12. https://doi.org/10.1016/j.jag.2010.05.009
Pavelko, MT (2004). Estimaciones de propiedades hidráulicas a partir de un modelo numérico unidimensional de deformación vertical del sistema acuífero Sitio Lorenzi,
Las Vegas, Nevada (p. 52). Departamento del Interior de EE. UU., Servicio Geológico de EE. UU. https://doi.org/10.3133/wri034083
PérezRocha, LE, Ordaz, M., & SánchezSesma, FJ (1995). Interpolación espacial de datos sísmicos: El caso del Valle de México. En Actas de la Décima Conferencia
Panamericana sobre Mecánica de Suelos e Ingeniería de Cimentaciones. Obtenido de http://classify.oclc.org/
clasificar2/ClassifyDemo?owi=5612459551
Polonia, JF e Irlanda, RL (1988). Hundimiento del terreno en el valle de Santa Clara, California, a partir de 1982 (p. 93). Departamento del Interior,
Servicio Geológico de EE. UU. https://doi.org/10.3133/ofr84818
Riley, FS (1998). Mecánica de los sistemas acuíferos: el legado científico de Joseph F. Polonia. En JW Borchers (Ed.), Land subsidence: Case Studies and current
Research, actas del simposio del Dr. Joseph F. Polonia sobre subsidencia de terrenos (págs. 1327). Publicaciones especiales de la Asociación de Geólogos
Ingenieros.
Rosen, PA, Gurrola, EM, Agram, PS, Sacco, GF y Lavalle, M. (2015). El entorno de computación científica InSAR (ISCE): un marco de Python para las ciencias de la
Tierra. En American Geophysical Union, resúmenes de la reunión de otoño de 2015. Obtenido de https://agu.confex.com/
agu/fm15/meetingapp.cgi/Paper/60264
Rudolph, DL, Cherry, JA y Farvolden, RN (1991). Flujo de agua subterránea y transporte de solutos en arcilla lacustre fracturada cerca de México
Ciudad. Investigación de recursos hídricos, 27, 2187–2201. https://doi.org/10.1029/91wr01306
Rudolph, DL, Sultan, R., Garfias, J. y McLaren, RG (2006). Importancia de la mayor infiltración debido a la extracción de agua subterránea en la desaparición de un
sistema lagunar de cabecera: Cuenca de Toluca, México. Revista de hidrogeología, 14, 115130. https://doi.org/10.1007/
s1004000504634
Ruiz, G. (2015). Estimación de la recarga freática en el acuífero de la Ciudad de México. Ciencias ambientales de Procedia, 25, 220–226. https://doi.org/10.1016/
j.proenv.2015.04.030
Sahakian, VJ, Melgar, D., Quintanar, L., RamírezGuzmán, L., PérezCampos, X., & Baltay, A. (2018). Movimientos del suelo de los terremotos del 7 y 19 de septiembre
de 2017 en Tehuantepec y PueblaMorelos, México. Boletín de la Sociedad Sismológica de América, 108, 3300–
3312. https://doi.org/10.1785/0120180108
Samsonov, SV y d'Oreye, N. (2017). Subconjunto de línea de base pequeña multidimensional (MSBAS) para análisis de deformación bidimensional: caso
estudiar Ciudad de México. Revista Canadiense de Teledetección, 43, 318–329. https://doi.org/10.1080/07038992.2017.1344926
Santoyo, MA, Singh, SK y Mikumo, T. (2005). Proceso de fuente y cambio de tensión asociado al 11 de enero de 1997 (Mw= 7,1)
Michoácán, México, terremoto de losa. Geofísica Internacional, 44, 317–330.
Sistema de Aguas de la Ciudad de México. (2012). El gran reto del agua en la Ciudad de México. Pasado, presente y prospectivas de solución
para una de las ciudades más complejas del mundo. Sistema de Aguas de la Ciudad de México.
SolanoRojas, D., CabralCano, E., FernándezTorres, E., Havazli, E., Wdowinski, S. y SalazarTlaczani, L. (2020b). Aceleración de subsidencia desencadenada
remotamente en la Ciudad de México inducida por los terremotos de Mw 7.1 de Puebla de septiembre de 2017 y los terremotos de Mw de 8.2 de Tehuantepec de
septiembre de 2017. Actas de la Asociación Internacional de Ciencias Hidrológicas, 382, 683–687. https://doi.org/10.5194/
piahs3826832020
Solano Rojas, D., Cabral Cano, E., Hernández Espriú, A., Wdowinski, S., DeMets, C., Salazar Tlaczani, L., et al. (2015). La relación de hundimiento del terreno InSAR
GPS y el abatimiento del nivel estático en pozos de la zona Metropolitana de la Ciudad de México. Boletín de la Sociedad Geológica Mexicana, 67, 273–283. https://
doi.org/10.18268/bsgm2015v67n2a10
SolanoRojas, DE, Wdowinski, S., CabralCano, E., Osmanoglu, B., Havazli, E. y PachecoMartínez, J. (2020a). Un enfoque multiescala para la detección y mapeo de
hundimientos diferenciales utilizando productos InSAR multiplataforma. Actas de la Asociación Internacional de Ciencias Hidrológicas, 382, 173–177. https://doi.org/
10.5194/piahs3821732020
Stead, D. y Pojani, D. (2017). La crisis del transporte urbano en las economías emergentes: una visión comparada. La crisis del transporte urbano en
economías emergentes (págs. 283295). https://doi.org/10.1007/9783319438511_14
Strozzi, T. y Wegmüller, U. (1999). Hundimiento del terreno en la Ciudad de México mapeado mediante interferometría diferencial SAR ERS. IEEE 1999 internacional
Simposio internacional de geociencia y teledetección (págs. 19401942), IEEE. https://doi.org/10.1109/IGARSS.1999.774993
Terzaghi, K., Peck, RB y Mesri, G. (1996). La mecánica de suelos en la práctica de la ingeniería (p. 592). John Wiley e hijos.
Tien Bui, D., Shahabi, H., Shirzadi, A., Chapi, K., Pradhan, B., Chen, W., et al. (2018). Mapeo de susceptibilidad a hundimientos de tierras en el sur
Corea utiliza algoritmos de aprendizaje automático. Sensores, 18, 2464. https://doi.org/10.3390/s18082464
Tortajada, C. (2006a). La gestión del agua en el área metropolitana de la Ciudad de México. Desarrollo de recursos hídricos, 22, 353–376. https://doi.
org/10.1080/07900620600671367
Tortajada, C. (2006b). ¿Quién tiene acceso al agua? Estudio de caso del área metropolitana de la Ciudad de México (Nº HDOCPA200616) (p. 47). Desarrollo humano
Oficina de Informe de Desarrollo (HDRO), Programa de las Naciones Unidas para el Desarrollo (PNUD).
Viens, L., Denolle, MA, Hirata, N. y Nakagawa, S. (2018). Reología compleja cerca de la superficie inferida de la respuesta del gran Tokio a fuertes movimientos del suelo.
Revista de investigación geofísica: Tierra sólida, 123, 5710–5729. https://doi.org/10.1029/2018jb015697
Wilson, AM y Gorelick, S. (1996). Los efectos del bombeo pulsado sobre el hundimiento del terreno en el Valle de Santa Clara, California. Diario de
Hidrología, 174, 375–396. https://doi.org/10.1016/00221694(95)02722x
Yunjun, Z., Fattahi, H. y Amelung, F. (2019). Análisis de series temporales de InSAR de línea de base pequeña: corrección de errores de desenvolvimiento y reducción de
ruido. Computadoras y geociencias, 133, 104331. https://doi.org/10.1016/j.cageo.2019.104331

Chaussard Et Al. 2021 (Esp)

Cargado por

Copyright:

Formatos disponibles

Chaussard Et Al. 2021 (Esp)

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Chaussard Et Al. 2021 (Esp)

Cargado por

Copyright:

Formatos disponibles

Machine Translated by Google

ARTÍCULO DE INVESTIGACIÓN Más de un siglo de hundimiento en la Ciudad de México: sin esperanzas

1.1. Estratigrafía del Valle de México

1.2. Agua subterránea en la ZMVM

1.3. Historia del hundimiento del terreno en la ZMVM

1.4. Procesos de hundimiento

1.5. Encuestas InSAR anteriores

2.2. Método para aislar la posible deformación elástica

3.1. Subsidencia rápida y constante desde 1950

3.2. Ausencia de deformación elástica “oculta” a corto plazo

3.3. Controles sobre tasas de hundimiento y extensión espacial

3.4. Pronóstico de deformación futura

4.2. Contexto socioeconómico: daños a la infraestructura, niveles de pobreza y acceso al agua

También podría gustarte