Evidencia GNSS+DInSAR preliminar de deformación co-sísmica en Venezuela asociada a los terremotos del 24 de junio de 2026.

Resumen

Se presenta una caracterización preliminar del desplazamiento superficial asociado a los eventos sísmicos Mw 7,0–7,5 del 24 de junio de 2026, en Venezuela, mediante pares interferométricos Sentinel-1 procesados con HyP3/GAMMA y observaciones de la estación GNSS CCS1. Los desplazamientos en línea de vista (LOS) se combinaron en un ajuste ponderado regularizado con los desplazamientos medidos en la estacion GNSS CCS1 para estimar las componentes este, norte y vertical (dE, dN y dU). Debido a la baja sensibilidad norte–sur de Sentinel-1, dN depende principalmente de la restricción GNSS y no constituye una estimación DInSAR independiente. La solución indica un movimiento regional dominante hacia el oeste -0.482 m en la zona de estudio DInSAR. La concordancia DInSAR y GNSS en CCS1 de 0.018 m, mide consistencia interna, no validación externa. Los resultados requieren validación GNSS independiente, corrección atmosférica, control de la referencia de fase y propagación de covarianzas antes de una interpretación geodésica definitiva.

Palabras clave: DInSAR; GNSS; Sentinel-1; desplazamiento cosísmico; desplazamiento LOS; Venezuela.

1. Antecedentes

El 24 de junio de 2026, la región centro-norte de Venezuela fue afectada por varios terremotos de gran magnitud. El primer evento ocurrió a las 22:04:34 UTC, con magnitud Mw 7,2, epicentro en 10,409°N, 68,557°O y profundidad de 27,6 km. El segundo ocurrió a las 22:05:11 UTC, con magnitud Mw 7,5, epicentro en 10,509°N, 68,504°O y profundidad de 10,0 km. Un tercer evento ocurrió a las 22:54:00 UTC, con magnitud Mw 5.0, epicentro aproximado en 10.885°N, 66.973°W y profundidad cercana a 10.0 km. Estos parámetros corresponden al catálogo preliminar del USGS y pueden revisarse conforme avance el análisis sismológico (USGS, 2026a, 2026b). Los epicentros se localizaron al oeste de Caracas, en la zona de interacción entre la placa del Caribe (CP), el bloque de los Andes del Norte (NA) y la placa Suramericana (SA).

El procesamiento GNSS preliminar de CCS1, ubicada en Caracas, estimó un salto cosísmico predominantemente hacia el oeste. En el sistema topocéntrico local ENU —este (E), norte (N) y arriba (U)— se obtuvieron DE = −0,464 m, DN = -0,011 m y DU = +0,021 m (Arenas, I., Suárez H., 2026). Esta convención se utiliza para las componentes GNSS topocéntricas en todo el artículo. El estudio complementa la observación puntual GNSS con interferometría SAR diferencial (DInSAR) obtenida de imágenes Sentinel-1 anteriores y posteriores a los eventos sísmicos. DInSAR mide cambios relativos de distancia en la dirección de línea de vista del satélite. La combinación de geometrías ascendentes y descendentes del Sentinel-1 mejora la separación este–vertical, pero proporciona poca sensibilidad a la componente norte (Wright, Parsons y Lu, 2004; Brouwer y Hanssen, 2024).

2. Datos SAR utilizados

Se utilizaron imágenes Sentinel-1 de banda C (λ ≈ 0,0555 m), nivel 1 Single Look Complex (SLC), modo Interferometric Wide Swath y polarización VV. Los cambios de vegetación, humedad y condiciones superficiales pueden reducir la coherencia; por ello, la fase no debe interpretarse automáticamente como desplazamiento del terreno (ASF DAAC, 2026a). Los pares incluyen adquisiciones de Sentinel-1A y Sentinel-1D. Los datos de Sentinel-1D estuvieron disponibles para usuarios desde abril de 2026 (ESA, 2026).

Las adquisiciones corresponden al 6, 13, 18 y 25 de junio de 2026 y cubren sectores de Aragua, Carabobo, La Guaira, Miranda y Yaracuy, además del Distrito Capital. Se generaron cinco pares interferométricos: dos ascendentes y tres descendentes (Tablas 1 y 2). Estos pares no representan cinco geometrías independientes, porque comparten trayectorias similares y una única fecha posterior.

                                                   Tabla 1. Imágenes SAR Sentinel-1 (SLC) empleadas.

                              Tabla 2. Geometría de imágenes SAR Sentinel-1 (SLC) empleadas.

3. Procesamiento DInSAR

Los pares Sentinel-1 IW-SLC se procesaron en ASF HyP3 con hyp3_gamma 9.1.0 y GAMMA 20240627 (Hogenson et al., 2020; ASF DAAC, 2026a, 2026c). Se utilizaron polarización VV, órbitas restitutas, 20 looks en rango, 4 looks en azimut, filtro adaptativo ADF = 0,6, desenrollado Minimum Cost Flow, DEM Copernicus GLO-30 y geocodificación WGS 84 / UTM zona 19N. La configuración 20×4 produce una resolución nominal aproximada de 160 m y un espaciado de píxel de 80 m (Tablas 3a y 3b). El flujo de trabajo incluyó calibración de los SLC, preparación del DEM, mosaico de subfranjas y ráfagas, corregistro, refinamiento mediante diversidad espectral mejorada, remoción de las fases de tierra plana y topográfica, filtrado adaptativo, estimación de coherencia, desenrollado y geocodificación. Los interferogramas conservan contribuciones de deformación, atmósfera, errores orbitales y del DEM, decorrelación y posibles errores de desenrollado.

La fase topográfica simulada con el DEM y los vectores orbitales se sustrajo de cada interferograma. Los productos se generaron con órbitas restitutas. Un reprocesamiento posterior con órbitas precisas deberá evaluar cambios en rampas y residuos. El corregistro se refinó iterativamente antes de aplicar el filtro ADF y el desenrollado MCF. No se aplicaron filtros de banda en rango o azimut, filtro de moteado ni un umbral explícito de desenrollado. La calidad se evaluó mediante coherencia, continuidad espacial e inspección de residuos; el parámetro ADF no representa una incertidumbre.

                                                Tabla 3a. Parámetros de procesamiento DInSAR.

                       Tabla 3b. Metadatos específicos y referencia automática de los productos HyP3.

Los productos se geocodificaron en WGS 84 / UTM zona 19N. Los archivos analizados incluyeron fase envuelta y desenrollada, desplazamiento LOS, coherencia, amplitud, vectores de mirada, ángulo de incidencia, DEM, máscara de agua y metadatos del procesamiento. El producto vertical automático de HyP3 se conservó sólo como archivo auxiliar y no se utilizó como observación independiente. Los productos principales se distribuyeron en formatos GeoTIFF y KMZ, junto con los archivos de parámetros necesarios para reproducir y auditar el procesamiento.

Las Figuras 1 y 2 muestran ejemplos de fase envuelta para una geometría ascendente (UID 2450) y una descendente (UID BA14). Cada franja de 2π equivale aproximadamente a 2,8 cm de cambio LOS para Sentinel-1; la densidad de franjas indica gradientes de fase, pero no permite por sí sola separar deformación, atmósfera y errores de desenrollado (ASF DAAC, 2026a). En la Figura 1b se observa el patrón típico de deformación co-sísmica al norte del estado Yaracuy (Agua Negra), próximo a ambos epicentros; la alta concentración de franjas indica desplazamientos superficiales de mayor magnitud en esa área.

   Figura 1a. Interferograma co-sísmicos ascendente (2450) – wrapped phase (18-25/JUNIO/2026).

                 Figura 1b. Zona de epicentros.                          Figura 1c. Zona estados La Guaira, Miranda.

Las Figuras 2 (a y b) muestran que en los estados La Guaira y Miranda también se observan gradientes de fase, aunque con una concentración de franjas interferométricas menor en comparación con la zona epicentral. Estas figuras se interpretan como evidencia cualitativa de una señal regional que debe confirmarse con fase desenrollada, coherencia, comparación entre pares y controles atmosféricos. Esta información DInSAR preliminar permitió identificar a priori las zonas de mayor desplazamiento superficial relativo causado por el terremoto del 24 de junio de 2026.

   Figura 2a. Interferograma co-sísmicos descendente (BA14) – wrapped phase (18-25/JUNIO/2026).

                                               Figura 2b. Zona estados La Guaira, Miranda.

4. Referencia común e inversión GNSS–DInSAR

Las mediciones DInSAR representan cambios relativos de distancia en la dirección de línea de vista (LOS). Antes de comparar o combinar los interferogramas, los cinco productos interferometricos se expresaron respecto de una referencia común. HyP3 selecciona automáticamente píxeles de alta coherencia para el desenrollado, pero alta coherencia no garantiza estabilidad. Las referencias originales se documentan en la Tabla 3b.

Los mapas LOS corregidos se presentan en las Figuras 3 y 4. En la geometría ascendente se registraron valores cercanos a −0,43 m en el sector Yaracuy–Falcón–Carabobo y a −0,51 m en La Guaira. En la geometría descendente se observaron valores próximos a −0,23 m en La Guaira y Miranda. Con la convención de HyP3, los valores negativos representan movimiento alejándose del sensor. Estos resultados derivan de la fase desenrollada y deben interpretarse junto con la coherencia y los errores residuales.

                        Figura 3. Mapa Line-of-Sight (LOS) co-sísmico ascendente (18-25/JUNIO/2026).

               Figura 4. Mapa Line-of-Sight (LOS) co-sísmico descendente (18-25/JUNIO/2026).

Los productos LOS ascendentes y descendentes se combinaron para estimar dE, dN y dU mediante mínimos cuadrados ponderados. Esta convención (dE, dN y dU) se utiliza para las componentes DInSAR en todo el artículo. Las observaciones con mayor coherencia recibieron mayor peso. La geometría casi polar de Sentinel-1 y su combinación orbital determina mejor dE y dU, pero ofrece poca sensibilidad al movimiento norte–sur; por tanto, dN depende principalmente de fuentes externas.

                                    Figura 5. Geometría LOS y sus componentes (dE, dN, dU).

Por esta razon, se realizaron dos cálculos por separado: un cálculo “libre” basado exclusivamente en datos DInSAR, combinando órbitas ascendentes y descendentes, donde se asumió dN = 0 m; y un segundo cálculo “regularizada”, donde DInSAR se combinó con los desplazamientos conocidos en la estacion GNSS CCS1 para lograr obtener información en la componente dN. El peso del GNSS fue mayor que el de una observación LOS individual para estabilizar la solución.

Para el primer caso “libre”, la solución DInSAR produjo dE = −0,496 m, dN = 0 m y dU = +0,094 m. En el segundo caso “regularizada” por CCS1, la solución DInSAR produjo dE = −0,475 m, dN = −0,008 m y dU = +0,043 m. La proximidad con CCS1 representa consistencia interna, no validación independiente, porque la estación participó en la inversión. Para efectos de este trabajo, los resultados que se consideraron para los pasos subsecuentes fueron los de la solución GNSS–DInSAR regularizada con CCS1.

5. Comparación DInSAR versus GNSS

Se realizó un segundo procesamiento GNSS con órbitas precisas de la estación CCS1, que permitió calcular y corroborar el desplazamiento co-sísmico en Caracas, Venezuela. El procesamiento incluyó datos de tres días previos al terremoto (21 al 23 de junio de 2026) y tres días posteriores (25 al 27 de junio de 2026). Para evaluar el contexto regional, se procesaron bajo la misma metodología seis estaciones GNSS adicionales: Barinas (BRNS), Coro (COR1), Maracaibo (MARA) y Puerto Ayacucho (PRAY), pertenecientes a la red REMOS del IGVSB; y Aruba (CN19) y Curazao (CN40), pertenecientes a la red COCONET del GAGE Facility Data Center. La Figura 6 muestra la ubicación de las siete estaciones GNSS, los tres epicentros (M7.5, M7.2 y M5.0) y la zona de estudio DInSAR.

                                               Figura 6. Mapa de estaciones GNSS procesadas.

Se realizó el procesamiento por dos vías independientes: Posicionamiento de Punto Preciso (PPP) mediante el servicio CSRS-PPP de Natural Resources Canada (NRCan), y mediante el servicio Online Positioning User Service (OPUS) del National Geodetic Survey (NOAA/NGS). Se generaron soluciones estáticas y las diferencias de posición se expresaron en sus componentes topocéntricas GNSS (DE, DN, DU). Todos los resultados están referidos al ITRF2020 (época 2026.5). Los desplazamientos GNSS calculados (DE, DN, DU) se presentan en la Tabla 4.

                      Tabla 4. Desplazamiento co-sísmico para estaciones GNSS en la región.

CCS1 registra el desplazamiento dominante: DE = −0,463 m, DN = −0,007 m y DU = +0,03 m. Los valores milimétricos de las demás estaciones son comparables con la incertidumbre reportada y no demuestran por sí solos una respuesta co-sísmica significativa. Su interpretación requiere series temporales, covarianzas y pruebas estadísticas de cambio.

                    Tabla 5. Desplazamiento co-sísmico GNSS y DInSAR para estación CCS1.

La Tabla 5 compara la solución regularizada con el valor de CCS1 utilizado como restricción. La diferencia tridimensional es aproximadamente 0,018 m. Este valor cuantifica consistencia de calibración; no constituye validación externa ni sustenta un radio fijo de confianza. La validación debe realizarse con estaciones excluidas del ajuste y con incertidumbres propagadas.

6. Resultados exploratorios en puntos de interés (POI)

A partir de los resultados DInSAR, se calcularon desplazamientos horizontales y verticales en varios sitios de interés: aeropuertos, autopistas, centros poblados, distribuidores viales, estadios, puentes, puertos, túneles, viaductos, entre otras obras de ingeniería. Se seleccionaron 20 sitios de interés (POI) distribuidos en los estados Aragua, Carabobo, La Guaira, Miranda y Yaracuy. La Figura 7 (a, b c) muestra el mapa con la ubicación de los veinte POI seleccionados.

                                                  Figura 7a. Mapa con puntos de interés (POI).

          Figura 7b. POI cercanos a epicentros.                          Figura 7c: POI en La Guaira y Miranda.

La Tabla 6 presenta los vectores topocéntricos de desplazamiento DInSAR para dichos POI, ordenados ascendentemente por el desplazamiento horizontal total (dHz) acumulado del 6 al 25 de junio de 2026. Los desplazamientos horizontales (dHz) varían entre 0.432 m y 0.522 m, con un RMS de ±0.482 m. Los desplazamientos verticales (dU) varían entre 0.029 m y 0.048 m, con un RMS de ±0.038 m. Esta tabla incluye los resultados DInSAR para la estación GNSS_INTU en Caracas, y cuyos valores GNSS reportados por DIGA-UCV son DE= -0,446 m, DN= 0,062 m, DU= 0,006 m, lo cual confirma el desplazamiento observado hacia el oeste (DIGA, 2026).

         Tabla 6. Desplazamiento co-sísmico DInSAR en puntos de interés (POI), periodo 6-25/Jun/2026.

7. Interpretación preliminar

La distribución espacial DInSAR es compatible con un campo regional de deformación asociados a los eventos sismicos. La señal dominante corresponde a movimiento hacia el oeste en el corredor Caracas–La Guaira y coincide con CCS1. Las variaciones locales pueden incluir retrasos atmosféricos, decorrelación, efectos de la máscara de agua y errores topográficos o de desenrollado.

Los desplazamientos se concentran en los estados Aragua, Carabobo, Distrito Capital, Miranda, La Guaira y Yaracuy. Los centros poblados con mayores desplazamientos superficiales incluyen la Gran Caracas, Maiquetía, Catia La Mar, Caraballeda y Naiguatá.

El sistema de comunicación terrestre que conecta Caracas con La Guaira presenta los valores máximos de desplazamiento hacia el oeste, especialmente en la autopista Caracas–La Guaira, el Viaducto y los túneles Boquerón I y II, donde se observan deformaciones significativas. Estos resultados no implican necesariamente afectación estructural de dichas obras, pero deben considerarse en futuros estudios detallados. Esta coincidencia espacial justifica estudios geodésicos y geotécnicos específicos, pero no demuestra deformación estructural ni daño en obras individuales. Esa evaluación requiere mediciones in situ, resolución espacial adecuada y análisis de ingeniería.

8. Limitaciones

Los cinco interferogramas comparten el 25 de junio como única fecha posterior y utilizan fechas previas diferentes. En consecuencia, no separan completamente la deformación co-sísmica de posibles señales presísmicas, postsísmicas tempranas o atmosféricas. No se aplicó una corrección troposférica independiente ni se estimaron rampas orbitales residuales. Los pares presentan dependencia temporal y espacial, y sus errores no deben tratarse como independientes. La única estación GNSS utilizada en la inversión fue CCS1; por ello, el producto carece de validación externa.

9. Comentarios finales

Este estudio geodésico es el primero que integra observaciones GNSS y DInSAR para determinar desplazamientos co-sísmicos preliminares en la región Caracas–La Guaira asociados al terremoto del 24 de junio de 2026. Los cinco pares Sentinel-1 analizados, dos ascendentes y tres descendentes, muestran una señal regional compatible con la deformación asociada a los terremotos del 24 de junio de 2026. La combinación de geometrías orbitales mejora la separación de las componentes este–oeste y vertical; sin embargo, la fecha posterior común y la ausencia de corrección atmosférica impiden atribuir toda la señal exclusivamente al desplazamiento co-sísmico.

La solución preferida está regularizada únicamente con la estación CCS1 y muestra un movimiento dominante hacia el oeste. La componente norte no está resuelta de forma independiente por DInSAR; su estabilidad proviene principalmente del prior GNSS. Por ello, los productos deben denominarse “solución GNSS–DInSAR regularizada con CCS1” y deben acompañarse de un análisis de sensibilidad. Para ciertas interpretaciones, una solución bidimensional este–vertical (E–U) puede ser más defendible que una tridimensional. La diferencia de aproximadamente 0.018 m entre la solución regularizada y CCS1 refleja consistencia de calibración, no una validación independiente.

Es poco probable que la tendencia general de desplazamiento (–0.482 m en la componente horizontal y +0.032 m en la vertical) cambie significativamente tras un refinamiento riguroso. Estudios independientes sismogeodésicos corroboran tendencias muy similares.

Por consiguiente, la magnitud del desplazamiento obliga a una pronta actualización y modernización del marco de referencia geodésico nacional (geocéntrico) y del sistema de alturas (geopotencial), particularmente por la deformación en las inmediaciones del mareógrafo de La Guaira, que sirve de referencia para el datum vertical de Venezuela. Asimismo, se recomienda aumentar el número de estaciones GNSS de observación continua (REMOS) y desplegar una red de reflectores de esquina (Corner Reflectors) para maximizar la integración GNSS–DInSAR y posibilitar el monitoreo integral de la geodinámica del país.

10. Referencias

• Alaska Satellite Facility DAAC (ASF DAAC) (2026a) Sentinel-1 InSAR Product Guide. Available at: https://hyp3-docs.asf.alaska.edu/guides/insar_product_guide/ (Accessed: 12 July 2026).

• Alaska Satellite Facility DAAC (ASF DAAC) (2026b) How to cite HyP3. Available at: https://hyp3-docs.asf.alaska.edu/usage_guidelines/ (Accessed: 12 July 2026).

• Alaska Satellite Facility DAAC (ASF DAAC) (2026c) Sentinel-1 On Demand InSAR product metadata for products 4C6B, D83D, 5AB8, BA14 and 2450 [data set metadata]. NASA Earthdata Cloud. Processing dates 3–4 July 2026.

• Arenas Bermúdez, I. and Suárez Acosta, H.D. (2026) Evidencia GNSS preliminar de deformación cosísmica en Caracas asociada al terremoto del 24 de junio de 2026. Unpublished technical report. Maracaibo: Laboratorio de Geodesia Física y Satelital Dr. Melvin Hoyer, Universidad del Zulia.

• Brouwer, W.S. and Hanssen, R.F. (2024) ‘Estimating three-dimensional displacements with InSAR: the strapdown approach’, Journal of Geodesy, 98(12), article 110. doi: 10.1007/s00190-024-01918-2.

• European Space Agency (ESA) (2026) Sentinel-1D user data opening from 17 April 2026. Available at: https://sentinels.copernicus.eu/-/sentinel-1d-user-data-opening-from-17-april-2026 (Accessed: 12 July 2026).

• DIGA (2026) Reporte Preliminar de la evaluación sismo-geodésica del evento del 24-jun-2026. Departamento de Ingeniería Geodésica y Agrimensura de la UCV. Universidad Central de Venezuela. Facultad de Ingeniería. Departamento de Ingeniería Geodésica y Agrimensura.

• Hogenson, K., Kristenson, H., Kennedy, J., Johnston, A., Rine, J., Logan, T., Zhu, J., Williams, F., Herrmann, J., Smale, J. and Meyer, F. (2020) Hybrid Pluggable Processing Pipeline (HyP3): A cloud-native infrastructure for generic processing of SAR data [computer software]. Zenodo. doi: 10.5281/zenodo.4646138.

• National Geodetic Survey (NGS) (2026) Online Positioning User Service (OPUS). Silver Spring, MD: National Oceanic and Atmospheric Administration. Available at: https://geodesy.noaa.gov/OPUS/ (Accessed: 12 July 2026).

• Natural Resources Canada (NRCan) (2026) CSRS-PPP updates and information. Available at: https://webapp.csrs.nrcan.gc.ca/geod/tools-outils/ppp-info.php?locale=en (Accessed: 12 July 2026).

• United States Geological Survey (USGS) (2026a) M 7.2 – 21 km ENE of San Felipe, Venezuela, 24 June 2026. Available at: https://earthquake.usgs.gov/earthquakes/eventpage/us6000t7zc/executive (Accessed: 12 July 2026).

• United States Geological Survey (USGS) (2026b) M 7.5 – 20 km ESE of Yumare, Venezuela, 24 June 2026. Available at: https://earthquake.usgs.gov/earthquakes/eventpage/us6000t7zp/executive (Accessed: 12 July 2026).

• Wright, T.J., Parsons, B.E. and Lu, Z. (2004) ‘Toward mapping surface deformation in three dimensions using InSAR’, Geophysical Research Letters, 31(1), L01607. doi: 10.1029/2003GL018827.

11. Agradecimientos

Los autores agradecen a ASF DAAC, NASA Earthdata Cloud y ESA/Copernicus por el acceso a los datos Sentinel-1 y al procesamiento HyP3. Los productos HyP3 contienen datos Copernicus Sentinel modificados de 2026, procesados por ESA. Al IGVSB/REMOS y EarthScope/GAGE por los datos GNSS. Al servicio CSRS-PPP de Natural Resources Canada (NRCan) y al servicio OPUS del National Geodetic Survey de USA (NOAA/NGS) por el procesamiento de las observaciones GNSS, y al United States Geological Survey (USGS) por el acceso al catálogo sísmico preliminar.