Análisis geocientífico de colapsos estructurales en la ciudad de Juan C. Bonilla, Puebla, México (2021)

Geoscientific Analysis of Structural Collapses in the City of Juan C. Bonilla, Puebla, Mexico (2021)

Análise Geocientífica do Colapso Estrutural na Cidade De Juan C. Bonilla, Puebla, México (2021)

Rogelio Ramos Aguilarn1

Benemérita Universidad Autónoma de Puebla, Puebla,México

Patricia Máximo Romero2

Benemérita Universidad Autónoma de Puebla, Puebla, México

Máximo Ávila Cruz3

Benemérita Universidad Autónoma de Puebla, Puebla, México

Víctor Galindo López4

Benemérita Universidad Autónoma de Puebla, Puebla, México

María de los Dolores Ramírez Álvarez5

Benemérita Universidad Autónoma de Puebla, Puebla, México

Lorena Cárdenas López6

Benemérita Universidad Autónoma de Puebla, Puebla, México

Argenis Eduardo Morales González7

Benemérita Universidad Autónoma de Puebla, Puebla, México

María de la Cruz Vázquez García8

Imágenes Geográficas S.A. México

Victoriano Covarrubias Salvatori9

Consejo de Ciencia y Tecnología del Estado de Puebla, México

Agradecimientos

Se agradece el apoyo del Consejo de Ciencia y Tecnología del estado de Puebla para llevar a cabo el estudio, bajo convenio de colaboración CONCYTEP/214/2021.

Introducción

Los hundimientos geológicos existen en todo nuestro planeta, se originan en algunas ocasiones por excavación, filtración de agua o falta de ella, por un peso extra en la superficie, presentándose principalmente en zonas kársticas. La investigación se centra en el colapso estructural que se generó el 29 de mayo de 2021, en los campos de cultivo ubicados en la localidad de Santa María Zacatepec, perteneciente al municipio de Juan C. Bonilla, Puebla, México, dentro de la cuenca hidrológica del Río Alto Atoyac región hidrológica RH18 (INEGI, 2020).

Se realizaron estudios: geomorfológicos, geohidrológicos y geoestadísticos, con el objetivo de contar con datos e información estructural relacionada con la topografía, geología, geomorfología, geohidrología, cartografía y percepción remota del lugar; para comprender las posibles causas que originaron el colapso estructural (Máximo et al, 2017).

El trabajo se realizó con técnicas de prospección de campo y telemétricas, el cual se complicó, debido a que se desarrolló durante la temporada de lluvia, ya que el procesamiento de algunas de las imágenes satelitales, contaban con nubosidad arriba del 30%, lo que hacía imposible dar un seguimiento (Ramos et al, 2013).

Se realizaron cálculos de la prospección de campo y se compararon con los datos telemétricos para determinar si la precipitación de lluvia acumulada pudiera haber influido en el colapso.

Se encontró que la estructura colapsada pone en riesgo a la población que se encuentra cerca, ya que los resultados determinaron que probablemente exista un cambio paulatino de deformación, también se han generado pequeñas fisuras superficiales junto a la zona de estudio que suponen un mecanismo de falla continúa.

Los datos obtenidos, concluyen que la geodinámica terrestre y atmosférica de la zona, influyeron en el hundimiento geológico.

Antecedentes

El 29 de mayo de 2021 se presenta un colapso estructural en la localidad de Santa María Zacatepec del municipio de Juan C. Bonilla, ubicada a aproximadamente a 50 km de la capital del estado de Puebla. El trabajo de campo y gabinete cuenta con datos y salidas a la zona que abarcan hasta junio 26 de 2021; se determina el área de la superficie del hundimiento para realizar el desarrollo del cálculo del líquido que se desplaza dentro del mismo, se generaron modelos topográficos y geométricos y se obtuvo la media de la zona estudiada de 11,281.46 m2, cuyo dato sirvió para calcular el área interior del colapso en un modelo 3D, dando un valor promedio de 32,100.79 m2.

Para el análisis, se relaciona el modelo del geoide terrestre con el elipsoide terrestre, para poder desarrollar los modelos geoestadísticos y de riesgo geohidrológicos y geomorfológicos y graficar e interpretar los resultados.

Se utilizaron datos telemétricos de plataformas satelitales como RapidEye, SkySat y Landsat, con software especializado para procesar la información y obtener modelados cartográficos y gráficos en 2D y 3D; las coordenadas utilizadas son geográficas y UTM. Se localizaron 18 puntos como referencia para generar el estudio gravimétrico, geoestadístico y de gasto.

Los datos de longitud, latitud y altura elipsoidal se convirtieron a una altura ortométrica que es la distancia entre la superficie física de la Tierra y del geoide en forma vertical, se obtuvo con el software libre del Instituto Nacional de Estadística y Geografía (INEGI) la curva tangente a la dirección de la gravedad en cualquier punto (Ramos et al, 2012).

La localidad de Santa María Zacatepec se ubica en el municipio de Juan C. Bonilla, en el estado de Puebla, se localiza en la parte centro-oeste del estado. Sus coordenadas geográficas son: 19º 05’ 30’’ y 19 º 08’ 36’’ de latitud norte y 98º 18’ 24’’ y 98º 25’ 36’’ de longitud oeste (Figura 1).

Figura 1. Ubicación de localidad de Santa María Zacatepec, Juan C. Bonilla, Puebla

Elaboración propia

Metodología

Se aplica la aceleración gravitacional normal para obtener la aceleración observada en el punto, el modelo utilizado es el EIGEN-6C4, a continuación, se calcularon las reducciones de Bouguer y aire libre (incluyendo la eliminación de la masa entre el punto y el geoide e igualmente la reducción de aire libre) para calcular sus respectivas gravedades y anomalías.

En el proceso de imágenes satelitales se combinaron los canales y bandas de los sensores satelitales para determinar algunas características de la región; se utilizaron datos multitemporales de las plataformas SkySat, Landsat y RapidEye para analizar la geología, topografía y uso de suelo.

El procesamiento de imágenes reveló algunas de las características que rodean a la zona colapsada, identificándose varias zonas de cultivo con condiciones de suelo blando y geología de roca volcánica.

Cálculo gravimétrico de las anomalías

Durante el desarrollo gravimétrico (Ramos et al, 2021) se realizó un análisis de los datos obtenidos para determinar la diferencia entre los campos gravitatorios terrestres reales y normales del área de estudio, la anomalía o perturbación potenciales (T), diferencia entre potenciales de gravedad (W) y energía de gravedad normal (U), el modelo matemático se muestra a continuación (1).

T(x,y,z)=W(x,y,z)-U(x,y,z) (1)

Para el desarrollo de los cálculos se utilizó el modelo de referencia geodésico GRS80, que es un sistema formado por el elipsoide de referencia global y del campo gravitatorio (el geoide), en la práctica, los elipsoides de referencia GRS80 y WGS84 son los mismos; se aplicó la Fórmula Internacional de Gravedad Normal adoptada por la Sociedad Internacional de Geodesia (2).

γ = 9.780318 (1 + 0.0053024 sen2 φ – 0.0000058 sen22 φ) (2)

Donde φ es la latitud y γ la gravedad normal dada en m/s2.

Se aplicaron para el estudio modelos prácticos matemáticos (3), donde F es la reducción por aire libre en mgal (milésima de gal, que equivale a 1 cm/s2) y H es la altura.

F=0.000003086*H (3)

La anomalía de aire libre (ΔgAL), utiliza la gravedad observada (gobs), la gravedad normal (γ) y la reducción de aire libre (F) (4).

ΔgAL=gobs+F-γ (4)

En este análisis se utilizó el valor de ρ=2.67 gr/cm3, modelo matemático (5).

AB=0.1119*H (5)

La anomalía de gravedad de Bouguer representa los cambios de grosor en la corteza terrestre y la densidad a nivel cortical, modelo (6).

ΔgB=gobs-γ-F+A+B (6)

Los modelos matemáticos generales para calcular las anomalías son (7,8, 9, 10, 11, 12):

● Atracción Gravitatoria de la capa de Bouguer:

AB=2πGρh=0.1119 h mgales (7)

● Gravedad reducida por Bouguer:

gB=gAl+At-AB (8)

● Anomalía de Bouguer:

ΔgtB ΔgAl+At-AB (9)

ΔgB ΔgAl+At-AB (10)

ΔgP=g+δgF+δgP-γ (11)

δgF+δgP=0.1967*h (12)

Con los datos de la Tabla 1 se aplicaron los modelos matemáticos anteriores (tabla utilizada en el reporte presentado al Consejo de Ciencia y Tecnología del Estado de Puebla)

Tabla 1. Datos calcular y obtener las anomalías del aire, aire libre y Bouguer

Elaboración propia.

Resultados y análisis

Cartografía Satelital

Para los datos e imágenes multitemporales se emplearon las plataformas SkySat, Landsat y SkySat, para los cálculos geométricos se utilizó el software Math. En la Figura 2 y Figura 3 se muestra la combinación de bandas 4, 3, 2 (RGB) con resolución de 50 cm., registrándose un área promedio inicial de 11,281.46 m2 y un perímetro de 328.09 m. en la zona del colapso (Ramos et al, 2022).

Las coordenadas UTM utilizadas para el procesamiento de las imágenes son las siguientes:

Punto de control	Coordenada X	Coordenada Y
1	565780.42 E	2114947.23 N	Esquina superior izquierda
2	565940.12 E	2114948.50 N	Esquina superior derecha
3	565943.70 E	2114812.18 N	Esquina inferior derecha
4	565780.19 E	2114809.72 N	Esquina inferior izquierda

Figura 2. Área procesada con MultiSpec de imágenes de Landsat, RapidEye y SkySat de zona de colapso

Elaboración propia.

Figura 3. Área procesada con MultiSpec de imágenes de Landsat, RapidEye y SkySat de combinación de bandas en la zona de colapso

Elaboración propia.

Durante el desarrollo del estudio, se encontró que la geomorfología estructural colapsada se fue modificando en un lapso aproximado de un mes, observándose con los registros de las imágenes satelitales procesadas durante ese tiempo. La geodinámica terrestre y atmosférica de la zona, contribuyó para un mecanismo de modificación del lugar (Ramos et al, 2005).

Se analizaron y procesaron las variables meteorológicas de precipitación como factor principal de alteración geohidrológica, las salidas a campo sirvieron para corroborar con los sensores meteorológicos instalados las lecturas obtenidas con el desarrollo matemático analizado.

En las gráficas 1, 2 y 3 se muestran los datos del desarrollo estructural (perímetro, área y volumen) obtenidos por análisis espectral satelital, del 17 de mayo al 26 de junio.

Gráfica 1. Crecimiento del área en m2

Elaboración propia.

Gráfica 2. Crecimiento del perímetro en m

Elaboración propia.

Gráfica 3. Crecimiento del volumen en m3

Elaboración propia.

Estudio geológico y geomorfológico

Se identificaron las características de la zona de estudio en cuanto a sus geoformas encontradas (relieves y topografía), que indican el proceso o grado de evolución generado por procesos endógenos y exógenos, además del proceso antropogénico.

Se identificó que la zona desplomada ha sido utilizada en exceso para el desarrollo agrícola, extracción de materiales y agua, durante los últimos 100 años; el relieve identificado, muestra un paisaje plano producto de meteorización, generado por lluvia y viento principalmente; la lluvia y la radiación solar han contribuido al cambio de la geología dentro del proceso evolutivo.

Las geoformas identificadas son pliegues anticlinales y sinclinales producto de fallamiento originados por procesos tectónicos y volcánicos, originados por contacto dinámico, describiendo sus características en la Figura 4.

Figura 4. Características geológicas y geomorfológicas de la zona colapsada con procesamiento ortofotogramétrico: 1.-Efecto de compresión de la roca, 2.-Existe una importante cantidad de material metamórfico intemperizado, 3.-Posible deslizamiento por flujo, 4.-El almacenamiento de sedimentos suaves es producto de la agricultura, 5.- Se aglomeran las formaciones geológicas y planicies creadas por los procesos fluviales del río Actiopa

Elaboración propia.

Fotogrametría terrestre

En las prospecciones realizadas durante el análisis de la zona colapsada se realizó fotogrametría terrestre, algunas imágenes se presentan en la Figura 5.

Figura 5. Fotogrametría terrestre de las prospecciones realizadas durante el análisis de la zona colapsada

Elaboración propia.

Aceleración de la gravedad con respecto al agua, anomalías de aire-libre, Bouguer e incremento de la gravedad

En la Tabla 2 se muestran los resultados de la aceleración de la gravedad del agua (m/s2 y en mgal) y las anomalías de Bouguer y aire-libre, las dos últimas columnas son los resultados del incremento de la aceleración de la gravedad.

Tabla 2. Resultados de la aceleración de la gravedad del agua (m/s2 y en mgal) y las anomalías de Bouguer y aire-libre

Vértice	Latitud	Elevación	Gravedad (g) m/s2	Gravedad (g) (mgal)	Gravedad Reducida por Aire-Libre (mgal)	Módulo de la Gravedad Normal en el Elipsoide (m/s2)	Módulo de la Gravedad Normal en el Elipsoide (mgal)	Anomalía Aire-Libre (mgal)	Atracción Gravitatoria a la capa de Bouguer (mgal)	Gravedad Reducida por Bouguer (mgal)	Anomalía de Bouguer (mgal)	Anomalía de Bouguer (mgal)
1	19.125555	2224	9.77899995	977899.995	978586.3214	9.784544517	978454.4517	131.869651	248.8656	978337.4558	-116.995949	-116.995949
2	19.125447	2224	9.77899989	977899.989	978586.3154	9.784544397	978454.4397	131.8756848	248.8656	978337.4498	-116.9899152	-116.9899152
3	19.125363	2224	9.77899984	977899.984	978586.3104	9.7845443	978454.43	131.880378	248.8656	978337.4448	-116.985222	-116.985222
4	19.12523	2224	9.77899976	977899.976	978586.3024	9.784544146	978454.4146	131.8878088	248.8656	978337.4368	-116.9777912	-116.9777912
5	19.124913	2224	9.77899959	977899.959	978586.2854	9.784543799	978454.3799	131.9055185	248.8656	978337.4198	-116.9600815	-116.9600815
6	19.124813	2224	9.77899953	977899.953	978586.2794	9.784543683	978454.3683	131.9111055	248.8656	978337.4138	-116.9544945	-116.9544945
7	19.124924	2223	9.77900268	977900.268	978586.2858	9.784546897	978454.6897	131.5961291	248.7537	978337.5321	-117.1575709	-117.1575709
8	19.125297	2223	9.77900289	977900.289	978586.3068	9.784547315	978454.7315	131.5752903	248.7537	978337.5531	-117.1784097	-117.1784097
9	19.125563	2223	9.77900304	977900.304	978586.3218	9.784547614	978454.7614	131.5604291	248.7537	978337.5681	-117.1932709	-117.1932709
10	19.125647	2223	9.77900308	977900.308	978586.3258	9.784547701	978454.7701	131.5557365	248.7537	978337.5721	-117.1979635	-117.1979635
11	19.12553	2223	9.77900302	977900.302	978586.3198	9.784547575	978454.7575	131.5622729	248.7537	978337.5661	-117.1914271	-117.1914271
12	19.124866	2223	9.77900265	977900.265	978586.2828	9.784546834	978454.6834	131.5993693	248.7537	978337.5291	-117.1543307	-117.1543307
13	19.124952	2223	9.7790027	977900.27	978586.2878	9.784546932	978454.6932	131.5945646	248.7537	978337.5341	-117.1591354	-117.1591354
14	19.125097	2223	9.77900278	977900.278	978586.2958	9.784547093	978454.7093	131.5864638	248.7537	978337.5421	-117.1672362	-117.1672362
15	19.125277	2223	9.77900288	977900.288	978586.3058	9.784547294	978454.7294	131.5764076	248.7537	978337.5521	-117.1772924	-117.1772924
16	19.124835	2223	9.77900263	977900.263	978586.2808	9.784546797	978454.6797	131.6011013	248.7537	978337.5271	-117.1525987	-117.1525987
17	19.125266	2223	9.77900287	977900.287	978586.3048	9.784547278	978454.7278	131.5770223	248.7537	978337.5511	-117.1766777	-117.1766777
18	19.125091	2223	9.77900277	977900.277	978586.2948	9.78454708	978454.708	131.5867994	248.7537	978337.5411	-117.1669006	-117.1669006

Elaboración propia.

Cálculo y desarrollo del gasto hidrológico y cartográfico

El aforo o suma es el gasto hidrológico de un fluido que se traslada por un canal de forma geométrica regular o irregular, calculado en metros cúbicos sobre segundo (m3/s).

Para el análisis y estudio, se establecieron varios modelos matemáticos y geométricos terrestres, examinando diversos escenarios como límites de agua bajas y superficiales.

Para determinar el gasto del agua en la zona colapsada se aplicó el modelo general de gasto (13).

Q=AIC/360 (13)

Donde A = superficie de la cuenca en metros cuadrados, I = intensidad de la lluvia en milímetros por hora y C = es el coeficiente de escorrentía sin dimensiones.

Los datos se obtuvieron del Instituto Nacional de Estadística y Geografía (INEGI) con el Simulador de Flujos de Agua de Cuencas Hidrográficas (SIATL, 2019), las variables utilizadas dificultaron el desarrollo debido a la falta de datos meteorológicos en la zona, para el procesamiento de los datos se tomó con un modelo ideal geohidrológico de planicie; los datos se manejaron para la precipitación en mm/hr, se multiplicaron por 24 horas y después por 30 días.

Se tomó como base longitudinal del cauce el diámetro promedio del hundimiento geológico; cabe mencionar que para el estudio se tomaron en cuenta modelos ideales geométricos para calcular el gasto, tomando al colapso dentro del área de la subcuenca para fines prácticos.

Los datos de campo se compararon con los obtenidos por telemetría y estaciones terrestres meteorológicas.

Este proceso se realizó para los 18 puntos, tomado los datos promedios; en la tabla 3 se muestra el cálculo superficial para una escorrentía con un modelo de textura geomorfológica plana, la intensidad de la lluvia se toma como punto de referencia media para realizar el proceso, se realizaron una serie de cálculos con diversos valores para la intensidad de lluvia, el criterio tomado es utilizar escurrimientos parciales ya que se carece de datos previos meteorológicos, por lo que se tuvo que tomar lecturas de estaciones meteorológicas portátiles instaladas en 2 escuelas en el municipio de Calpan, Puebla, muy cercanas a la zona de estudio.

En la Tabla 3 se muestra el resultado del promedio superficial para el mes de junio de 2021, con intensidad de lluvia baja y escorrentía con factor fijo de 0.30, el gasto obtenido es de 0.00000273575 m3/s.

Tabla 3. Resultado del promedio superficial para el mes de junio de 2021

Datos para cálculo de promedio superficial de gasto de zona colapsada
Intensidad lluvia	1.05 mm/hr es igual a 0.0000002910000 m/s
Área subcuenca	11,281.46 m2
Coeficiente escorrentía	0.30
Longitud cauce	126 m
Gasto promedio superficial	0.00000273575 m3/s

Elaboración propia.

En la Gráfica 4 se aprecia el comportamiento del gasto en los primeros 4 puntos calculados, el punto 1 representa el lugar más bajo por lo que el gasto es mínimo, el punto 4 representa una ligera elevación por lo que el gasto respecto a la pendiente del terreno aumenta.

Gráfica 4. Comportamiento del gasto respecto a los primero 4 puntos procesados en la zona de estudio

Elaboración propia.

Es importante señalar que el proceso de filtración del agua superficial ocasionado por las lluvias que se presentaron en el lapso de un mes se procesó para analizar la escorrentía con un modelo cartográfico en 2d y 3d, además de analizar el perfil para determinar los puntos de mayor riesgo, donde se aprecia una dirección de norte a sur como se muestra en la figura 6.

Figura 6. El modelo cartográfico representa la dirección del agua (pendiente y aguas abajo) con el software libre Surfer, la trayectoria del agua tiene una longitud de 126 metros desde el punto más alto al más bajo y una dirección de norte a sur

Elaboración propia.

En la Figura 7, se apreciar el comportamiento del gasto en forma laminar de dos perfiles (3 y 4), la razón por la que se realizó este análisis fue para tener más información del comportamiento hidrológico dentro del colapso de la estructura.

Figura 7. El modelo estructural de la zona de estudio muestra una profundidad de 55 metros, en una longitud de entre 90 y 40 metros. Para analizar mejor el comportamiento laminar en la zona colapsada

Elaboración propia.

El proceso exógeno generado en el lugar es muy importante, ya que su dinámica atmosférica ha influido en la geomorfología del lugar, es decir el proceso evolutivo es producto no solo de la dinámica interna terrestre sino de la climática y meteorológica también Sinem et al, (2019).

Cálculo de la pendiente y su relación con la aceleración de la gravedad

Para la obtención de todos los datos, fue muy importante obtener las coordenadas de los 18 puntos analizados, así como calcular la pendiente o inclinación del punto superficial al punto más bajo del hundimiento; el procedimiento es dividirlo en cuatro cuadrantes colocados al norte, sur, este y oeste, y luego calcular la pendiente promedio de la estructura analizada.

El modelo o relación matemática para calcular la pendiente (14) es:

(14)

En la Gráfica 5 se presenta un modelo general de la aceleración de la gravedad de la zona analizada con relación a las elevaciones con la que se corroboró el adecuado procesamiento (a mayor altura menor la aceleración gravimétrica y a menor altura mayor la aceleración de la gravedad).

Gráfica 5. Los datos obtenidos corroboran que a mayor altura menor la aceleración gravimétrica y a menor altura mayor la aceleración de la gravedad

Elaboración propia.

Resultados y análisis

El análisis geocientífico del desplome estructural en el Municipio de Juan C. Bonilla en el estado de Puebla, México, pudo establecer la estratigrafía de la zona colapsada, con una planicie de 1000 metros a la redonda como zona de riesgo, referente a la zona derrumbada se determinó un mecanismo dinámico del basamento, con escurrimiento del acuífero promedio a 50 metros de profundidad; la zona de estudio se encuentra dentro de la cuenca del Río Alto Atoyac, con una geodinámica de varias décadas de desarrollo activo del manto acuífero, que con el paso del tiempo y el peso extra superficial y sobre explotación de agua, pudieron contribuir a generar un colapso.

Los datos obtenidos del estudio con modelos geométricos, matemáticos y satelitales son los siguientes: 44.5, 130.1 y 125.3 metros de profundidad, diámetro mayor y diámetro menor, respectivamente.

Los datos del perímetro, área y volumen se obtuvieron tomando un modelo geométrico elipsoidal, cuyos datos son: área superficial de 11,281.46 m2, perímetro de 401.036 m y volumen de 759,270.929 m3.

Se obtuvieron los promedios de: aceleración de la gravedad (9.77900181 m/s2), anomalía de Bouguer (-117.10757041 mgal o 0.00117108 m/s2), anomalía al aire libre (131.68342959 mgal o 0.00131683 m/s2), fuerza del agua (97.78857950 N/m3), gasto (0.00000273575 m3/s), presión atmosférica (773.138 hPa o mbar), pendiente en el terreno de la superficie (-1 m) y pendiente promedio del punto alto al bajo (-0.052631579 m).

Conclusiones

El análisis desde las geociencias del colapso estructural en el Municipio de Juan C. Bonilla en el estado de Puebla, México, pudo identificar de forma precisa la geología, geomorfología, hidrología y desarrollar la cartografía y geoestadística de la zona; se estableció que la estructura derrumbada tenía una morfología plana, y cuyo proceso geodinámico endógeno interrelacionó con procesos de intemperismo y meteorización como el agua, precipitación, viento, temperatura y radiación solar. Se concluye que la región analizada es un receptor de agua sustentada por la fractura de un ducto profundo natural que favoreció al hundimiento, generando oquedades (pequeñas cavernas y ductos hidrológicos de entrada y salida de agua).

La zona podría generar un riesgo importante para la población, ya que el proceso dinámico de la zona puede continuar durante varios meses o años, lo que podría ocasionar un nuevo colapso, siendo importante continuar con observaciones del lugar e instrumentar la zona con sensores de deformación para tener una idea de su proceso.

Este estudio puede complementarse con otros métodos geofísicos como los eléctricos y magnéticos, además de sismicidad inducida para determinar las características de vibración del suelo; hubo limitaciones al no permitir las autoridades volar drones por la cercanía del aeropuerto de Huejotzingo.

El presente estudio pretende contribuir a la prevención y gestión de riegos naturales en el estado, gracias al apoyo del Consejo de Ciencia y Tecnología del Estado de Puebla.

Referencias

Instituto Nacional de Estadística y Geografía, INEGI. (2019). Cálculo de gasto o caudal. http://antares.inegi.org.mx/analisis/red_hidro/PDF/CalculoGasto.pdf

Instituto Nacional de Estadística y Geografía, INEGI. (2020). Mapa digital. http://gaia.inegi.org.mx/mdm6/

Máximo-Romero, P., Ramos-Aguilar, R., Martínez-Ortiz, A. y Hernández-Sánchez, A. (2017). Daños estructurales en Ciudad de México, producidos por los sismos del 28 de julio de 1957 y el 19 de septiembre de 1985. Revista UNIMAR, 35(2), 161-186. https://revistas.umariana.edu.co/index.php/unimar/article/view/1533/pdf

Ramos-Aguilar, R., Rivera-Domínguez, A. y Mayorga-Raposo, R. (2005). Aplicación de imágenes satelitales para determinar el clima y la radiación solar en el estado de Puebla y la geohidrología de su zona norte. Ciencia Ergo Sum, 10(003), 283-294. https://www.redalyc.org/pdf/104/10410306.pdf

Ramos-Aguilar, R., Máximo-Romero, P., Narciso-Hernández, J., Mirón-Morales, M. y Beltrán-Cruz, M. (2012). Estudio geoestadístico para obtener la gravedad local, pendiente y cálculo hidrológico de las barrancas Xaltelulco, Tepeloncocone, Tenepanco, Colorada y Quimichule del volcán Popocatépetl. Boletín de Ciencias de La Tierra, (31), 65–84. https://revistas.unal.edu.co/index.php/rbct/article/view/31254

Ramos-Aguilar, R., Máximo-Romero, P., Soto-Cruz, B., Alcántara-Iniesta, S. y Pérez-Luna, J. (2013). Analysis of the Meteorological Variables for Puebla City 2011-2012 Applying the Modeling Ion-Wavelets in a Hypothetical Manner. Atmospheric and Climate Sciences, 3(4), 450–458. https://doi.org/10.4236/acs.2013.34046

Ramos-Aguilar, R., Máximo-Romero, P., Galindo-López, V., Ávila-Cruz, M., Yáñez-Pérez, G., Márquez-Máximo, R., y Covarrubias-Salvatori, V. (2021). Análisis técnico geoestadístico de los ríos Atoyac y Alseseca que atraviesan la ciudad de Puebla, México. Contactos, Revista de Educación en Ciencias e Ingeniería, (119), 5-22. https://contactos.izt.uam.mx/index.php/contactos/article/view/99

Ramos-Aguilar, R., Máximo-Romero, P., Galindo-López, V., Ávila-Cruz, M., Yáñez-Pérez, G., Márquez-Máximo, R., Hernández-Andrade, D. y Covarrubias-Salvatori, V. (2022). Análisis gravimétrico y cálculo de la fuerza de un líquido en el volcán Malinche, Puebla-Tlaxcala, México. Boletín de Ciencias de la Tierra, (49), 9-14. https://doi.org/10.15446/rbct.n49.91906

SIATL. (2019). Institutito Nacional de Estadística y Geografía. http://antares.inegi.org.mx/analisis/red_hidro/siatl

Sinem, E., Barthelmes, F., Reißland, S., Elger, K., Förste, C., Flechtner, F. y Schuh, H. (2019). ICGEM – 15 years of successful collection and distribution of global gravitational models, associated services, and future plans. Earth System Science Data, 11(2), 647–674. https://doi.org/10.5194/essd-11-647-2019

---