Resumen
El trabajo se basa en un análisis general desde las geociencias del colapso en los campos de cultivo de Santa María Zacatepec, de la ciudad Juan C. Bonilla. En el desarrollo se aplica la técnica de Percepción Remota, permitiendo interpretar combinaciones de bandas y determinar la geometría, geología, geomorfología y cartografía del sitio. Se presenta un análisis resumido del trabajo que se realizó, además de mostrar modelos de simulación del gasto hidrológico y la aceleración del agua en las pendientes estructurales de escorrentía, datos que se usaron para la cartografía y modelado gráfico finales. Se realizó un estudio geomorfológico, geológico y topográfico del lugar, con técnicas telemétricas satelitales y aéreas; también se utilizó software especializado que se aplicó para levantamientos geodésicos, realizando cartografía en modelos 2D y 3D del área de análisis.
Las prospecciones del lugar se realizaron con el objetivo de conocer la forma de la estructura e identificar áreas potencialmente riesgosas para la población cercana al colapso. Se obtuvieron ortofotos para calcular y determinar la pendiente por distanciaelevación, la aceleración-altura, el gasto y la fuerza del agua; que sirvió para concluir las posibles causas que ocasionaron el colapso, y contar con datos técnicos suficientes que se compartieron con las autoridades de protección civil.
Palabras clave: colapso; geomorfología; geoestadística; aceleración; modelos
Abstract
This paper presents a general analysis from the geosciences of structural collapses. The analysis was conducted in the crop fields of Santa María Zacatepec, in the city of Juan C. Bonilla, Mexico. The analysis implemented the Remote Sensing technique, facilitating the interpretation of band combinations and determining the geometry, geology, geomorphology and cartography of the site under study. A summary analysis of the work performed is presented, in addition to showing simulation models of the hydrological outflow and water acceleration in the runoff structural slopes, data that were used for the final mapping and graphic modeling. A geomorphological, geological, and topographical study of the site was conducted, using satellite and aerial telemetric techniques. Specialized software was also used for geodetic surveys, mapping in 2D and 3D models of the analysis area. Surveys of the site were carried out to determine the shape of the structure and identify potentially risk areas for the population near the collapse events. Orthophotos were obtained to calculate and determine the slope by distance-elevation, acceleration-height, water flow and force, which served to conclude the possible causes that led to the collapses and have sufficient technical data that were shared with the civil protection authorities.
Keywords: collapse; geomorphology; geostatistics; acceleration; models
Resumo
O estudo fundamenta-se em uma análise geral das geociências do colapso nos campos de cultivo de Santa María Zacatepec, na cidade de Juan C. Bonilla, México. A técnica de sensoriamento remoto é aplicada no desenvolvimento, permitindo a interpretação de combinações de bandas e a determinação da geometria, geologia, geomorfologia e cartografia do local. É apresentada uma análise resumida do trabalho realizado, além de modelos de simulação do gasto hidrológico e da aceleração da água nos taludes estruturais de escoamento. Esses dados foram utilizados para o mapeamento final e a modelagem gráfica. Foi realizado um estudo geomorfológico, geológico e topográfico do local, utilizando técnicas telemétricas aéreas e de satélite. Também foi utilizado um software especializado para levantamentos geodésicos e mapeamento em modelos 2D e 3D da área de análise. Os levantamentos do local foram realizados com o objetivo de conhecer o formato da estrutura e identificar áreas de risco potencial para a população próxima ao colapso. Foram obtidas ortofotos para calcular e determinar a inclinação por distância-elevação, aceleração-altura, fluxo de água e força, o que serviu para concluir as possíveis causas que levaram ao colapso e para ter dados técnicos suficientes que foram compartilhados com as autoridades de proteção civil.
Palavras-chave: colapso; geomorfologia; geoestatística; aceleração; modelagem
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).
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 sen 2 φ - 0.0000058 sen 22 φ) (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
Gravedad reducida por Bouguer: (7)
gB=gAl+At-AB
Anomalía de Bouguer: (8)
Δ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)
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
565780.42 E
2114947.23 N
Esquina superior izquierda
565940.12 E
2114948.50 N
Esquina superior derecha
565943.70 E
2114812.18 N
Esquina inferior derecha
565780.19 E
2114809.72 N
Esquina inferior izquierda
Área procesada con MultiSpec de imágenes de Landsat, RapidEye y SkySat de combinación de bandas en la zona de colapso
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.
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.
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
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.
Fotogrametría terrestre de las prospecciones realizadas durante el análisis de la zona colapsada
Aceleración de la gravedad con respecto al agua, anomalías de airelibre, 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.
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.
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.
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.
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
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.
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
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).
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
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
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
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
» 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/
» http://gaia.inegi.org.mx/mdm6/ -
Máximo-Romero, P., Ramos-Aguilar, R., Martínez-Ortiz, A. y HernándezSá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
» 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
» https://www.redalyc. org/pdf/104/10410306.pdf -
Ramos-Aguilar, R., Máximo-Romero, P., Narciso-Hernández, J., MirónMorales, 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
» 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
» 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
» 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
» 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
» 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
» https://doi.org/10.5194/essd-11-647-2019
Fechas de Publicación
-
Fecha del número
Jan-June 2024
Histórico
-
Recibido
01 Mar 2023 -
Acepto
30 Ago 2023