GGenLUZ produce SGM26v1.0: geoide 2026 de ultra-alta-resolución 30x30m para zonas afectadas por subsidencia en la COLM-Edo. Zulia

Investigación, Extensión, Datos/Soluciones/Modelos, Software – SGM26v1.0: Subsidence Geoid Model 2026, version 1.0, 30x30m ultra-high-resolution

Para mejorar la determinación de cotas de terreno con la técnica de la nivelación-GNSS en zonas de la COLM-Edo. Zulia afectadas por el fenómeno de subsidencia, para BMs PDVSA o vértices geodésicos cualesquiera, y con la finalidad ulterior de derivar tasas de hundimiento más precisas o de corregir rápida y rigurosamente la elevación de tales estaciones por los efectos tiempo-dependientes de la deformación vertical por subsidencia; GGenLUZ ha calculado un novedoso modelo geoidal regional, de extensión 1°x1° (~12100km²), de ultra-alta-resolución 30x30m y de alta-precisión ±3.3cm, denominado aquí SGM26v1.0: Subsidence Geoid Model 2026, version 1.0, ver Figuras 1-4.

Figura 1. DTM con zonas efectiva y de borde para SGM26v1.0: Subsidence Geoid Model 2026, version 1.0 [Acuña, 2026].

 

SGM26v1.0, como anteriores modelos de geoide calculados por GGenLUZ, ha sido determinado vía el cuasigeoide según la teoría de Molodensky [Heiskanen/Moritz, 1967]. Por primera vez para un modelo geoidal venezolano, en SGM26v1.0 se utiliza como modelo geopotencial de referencia una solución global de ultra-alta-resolución, a saber, G25X19W25 [Acuña, 2026], un modelo estático-compuesto formado por GOCO2025s [Oehlinger et al., 2025], XGM2019e [Zingerle et al., 2019] y WHU-CASM-UGM2025 [Liu et al., 2026], [https://icgem.gfz.de/], de muy alto grado y orden en su expansión armónico-esférica, i.e., n,m=11000 (resol. ~1.8x1.8km). Este detallado modelo geopotencial global comprende 60.5 millones de coeficientes tipo –C,S,sC,sS– (exac. 60516498 en 5.8 GB), dados originalmente en el sistema de marea terrestre permanente zero-tide [Mäkinen, 2021]. G25X19W25 es producto de la combinación en el dominio espectral de GOCO2025s(2-200) + XGM2019e(201-2159) + WHU-CASM-UGM2025(2160-11000).

Para calcular mediante síntesis cantidades geodésicas de referencia globales de superficie (i.e., anomalías de altura y anomalías de gravedad de 'aire-libre') durante el proceso de determinación del cuasigeoide SGM26v1.0 (Figura 3) vía la técnica 'remover-restituir', aquí se emplean los coeficientes geopotenciales G25X19W25 completos hasta el grado y orden máximo n,m=11000 de su expansión armónico-esférica, ahora convenientemente expresados en el sistema de marea terrestre permanente tide-free consistente con el elipsoide y campo de gravedad normal GRS80 [Moritz, 1980]. El nivel "0" del (cuasi)geoide lo establece el actual valor IAG del potencial de gravedad en el geoide W0=62636853.4 m²/s² [IAG, 2015]; ver Figura 2.

Figura 2. Cuasigeoide global G25X19W25 (n,m=11000, resol. ~1.8x1.8km) respecto al elipsoide GRS80 en SGM26v1.0.

Figura 3. Cuasigeoide gravimétrico SGM26v1.0 de resolución 30x30m, respecto al elipsoide GRS80.

Figura 4. Geoide híbrido final SGM26v1.0 (cuasigeoide gravimétrico ajustado localmente en 26 estaciones -BMs GPS/nivelación- PDVSA) de ultra-alta-resolución 30x30m, respecto al elipsoide GRS80/ITRF2020(2026.5), para áreas terrestres de la COLM-Edo. Zulia afectadas por subsidencia [Acuña, 2026].

Publicación en desarrollo . . .

Más información en:

Acuña, G. (2026): SGM26v1.0: geoide 2026 de ultra-alta-resolución 30x30m para zonas afectadas por subsidencia en la COLM-Edo. Zulia. Tópicos de Geodesia Geométrica. Junio 8-14, 2026. Cátedra Geodesia Geométrica (GGenLUZ). Dpto. de Geodesia Superior. Esc. de Ingeniería Geodésica. Fac. de Ingeniería. Universidad del Zulia (LUZ). Maracaibo, Venezuela.

Curso de extensión GGenLUZ: Altimetría satelital multimisión radar 'nadiral e interferométrica’ para monitoreo del espejo de agua en embalses venezolanos

Cursos/Talleres, Extensión – Curso de extensión 'on-line' GGenLUZ

GGenLUZ invita a estudiantes y profesionales de la ingeniería geodésica, geomática y geociencias afines, a participar en el curso de extensión, versión digital, titulado "Altimetría satelital multimisión radar 'nadiral e interferométrica' para monitoreo del espejo de agua en embalses venezolanos", que estará disponible para descarga desde este website a partir del lunes 8 de junio de 2026.

Poco tiempo después que la técnica satelital altimétrica entrara en su fase operacional permanente global con las misiones ERS-1 y TOPEX/Poseidon (1991/1992), comenzó su aplicación para el monitoreo de cuerpos de agua de superficie (ríos, lagos, embalses, tierras inundables) en zonas terrestres continentales. En la actualidad, hasta 9 misiones altimétricas (CryoSat-2, SARAL/AltiKa, Jason-3, Sentinel-3A/B, CFOSAT, Sentinel-6A(MF)/B y SWOT) registran global y permanentemente las variaciones de altura de esas superficies de agua continentales con precisiones de algunos centímetros, resolución espacial de hasta 100 m y tiempos de repetición de sus observaciones de pocos días. ... En este curso de 8 horas académicas, se describen los fundamentos de la técnica satelital altimétrica y su aplicación geodésica en el monitoreo de reservorios de agua dulce en Venezuela. Los tópicos a tratar incluyen: principios básicos y conceptos generales de la altimetría satelital radar 'nadiral e interferométrica'. Geometría de las observaciones. Superficies de referencia. Satélites y misiones altimétricas actuales. Generación y adquisición de datos altimétricos. Principales fuentes de error y su impacto sobre las observaciones. Procesamiento de datos altimétricos sobre zonas terrestres. Principales observables altimétricos. Correcciones geofísicas. Correcciones por errores instrumentales. Correcciones por el medio de propagación. Indicadores de calidad de las observaciones altimétricas. Selección y validación de datos. Combinación multimisión. Acceso a datos altimétricos multimisión. Principales bases globales de datos altimétricos multimisión. Descripción y acceso a productos altimétricos diversos: observaciones originales y pre-procesadas, cuadrículas y series de tiempo de variables altimétrico-dependientes obtenidas de procesamiento científico. Construcción de las observaciones corregidas de la superficie instantánea de aguas continentales. Generación de series de tiempo para el monitoreo del nivel instantáneo del agua en ríos, lagos y zonas inundables del país, y particularmente, del espejo de agua de los principales embalses venezolanos.

Para más información sobre el curso (p.ej., contenido, forma de pago, detalles para descarga, consultas, etc.) contacte al Prof. Gustavo Acuña a través del email gacuna@fing.luz.edu.ve, o del número WhatsApp +58-412-42.71.579.

Niveles actuales del espejo de agua en los embalses de Guri y Uribante según altimetría satelital 'radar-nadir' multimisión

Investigación, Extensión – Altimetría satelital en Guri y Uribante

La nota muestra los niveles actuales del espejo de agua en los embalses venezolanos de Guri (Edo. Bolívar) y Uribante (Edo. Táchira) para Mayo 19/20, 2026 como registrados por altimetría satelital 'radar-nadir' multimisión, ver Figuras 1 y 2.

Figura 1. Serie de tiempo de 34 años (actualizada al 19.05.2026) con las variaciones de elevación del espejo de agua en el embalse venezolano de Guri durante el periodo 1992-2026, derivada de observaciones altimétricas satelitales 'radar-nadir' multimisión WSE (TOPEX/Poseidon + Jason-1/2/3 + Sentinel-6A), y respecto al nivel medio del mar realizado por el geoide nacional VGM25v1.0 [Acuña, 2025]. La cota en m.s.n.m del espejo de agua del embalse más reciente registrada por altimetría multimisión, a saber al 19.05.2026, se estima sea 263.91 ± 0.05 m. El valor indica que el embalse actualmente presenta un 74% (aprox.) de su capacidad máxima, muy por encima de los niveles mínimos que registró durante las crisis eléctricas nacionales de los años 2003, 2010 y 2016.

Figura 2. Serie de tiempo de 18 años (actualizada al 20.05.2026) con las variaciones de elevación del espejo de agua en el embalse venezolano de Uribante durante el periodo 2008-2026, derivada de observaciones altimétricas satelitales 'radar-nadir' multimisión WSE (Jason-2/3 + Sentinel-6A), y respecto al nivel medio del mar realizado por el geoide nacional VGM25v1.0 [Acuña, 2025]. La cota en m.s.n.m del espejo de agua del embalse más reciente registrada por altimetría multimisión, a saber al 20.05.2026, se estima sea 288.53 ± 0.10 m. El valor prácticamente se corresponde a los niveles mínimos anuales que ha venido presentando el embalse en los últimos 13 años (desde 2013 aprox.); esto es cerca de un 22% por encima del nivel mínimo histórico que registró el embalse a mediados de 2009 (281 m).

Más información en:

Acuña, G. (2026): Niveles actuales del espejo de agua en los embalses de Guri y Uribante según altimetría satelital 'radar-nadir' multimisión. Tópicos de Geodesia Geométrica. Junio 1-7, 2026. Cátedra Geodesia Geométrica (GGenLUZ). Dpto. de Geodesia Superior. Esc. de Ingeniería Geodésica. Fac. de Ingeniería. Universidad del Zulia (LUZ). Maracaibo, Venezuela.

Altimetría satelital SWOT para monitoreo de subsidencia del dique costanero, periodo 2012.0-2026.1, región Lagunillas, COLM-Edo.Zulia

Investigación, Extensión – Monitoreo de subsidencia del dique costanero en Lagunillas con SWOT

La nota muestra resultados de la aplicación de observaciones altimétricas satelitales SWOT [https://podaac.jpl.nasa.gov/SWOT/] para el monitoreo de subsidencia de secciones seleccionadas del dique costanero en Lagunillas, COLM-Edo. Zulia, ver Figuras 1 y 2.

 

Como consecuencia del fenómeno progresivo de subsidencia, provocado por la extracción masiva de crudo en zonas de la COLM desde comienzos del siglo XX (1914 aprox.) [Jaeger et al., 1989], [Cardozo, 2016], la industria petrolera venezolana construyó a partir de 1923-1929, y desde esa época, monitorea, mantiene y sistemáticamente actualiza/mejora hasta el presente, un extenso muro costanero de contención de las aguas del Lago de Maracaibo. Esta mega-estructura proteje del ingreso de las aguas del lago -y de sus muy probables trágicas consecuencias- a las poblaciones costeras cuyas cotas de terreno, qué por la deformación de subsidencia, están por debajo del nivel lacustre (p.ej., hasta ca. de -10 m en áreas de máxima deformación), principalmente, zonas de Tía Juana, Ciudad Ojeda, Lagunillas, Bachaquero y San Timoteo / San Lorenzo.

 

El dique costanero de la COLM (Figura 1) es una de las más grandes obras de ingeniería hidráulica del país. El muro se extiende 48 km en longitud sobre la línea de costa regional, comprende 3 diques costaneros principales, 103 km de diques internos, 490 km de canales y 34 estaciones de drenaje.

Figura 1. Dique costanero, región Lagunillas, COLM-Edo. Zulia [Cardozo, 2016].

Aquí, aplicando mediciones altimétricas satelitales SWOT registradas en el primer trimestre de 2026 combinadas con mediciones convencionales de nivelación 2012 PDVSA, según el procedimiento descrito en [Acuña, 2026 ; https://ggenluz.blogspot.com/2026/04/determinacion-de-cotas-de-terreno-en.html], se derivan (y muestran) estimaciones de deformación vertical por subsidencia en secciones del dique costanero de Lagunillas (Figura 2).

 

El ensayo de esta técnica de monitoreo satelital, alternativa atractiva en términos de costos, tiempo de ejecución, horas de procesamiento y análisis, y precisión en resultados, respecto a la convencional nivelación geodésica, arrojó para el periodo 2012.0-2026.1 -aprox. 14 años sin mediciones de monitoreo geodésico formal por nivelación PDVSA-, valores de hundimiento por subsidencia en el dique costanero entre -38.8 y +3.1 cm, con valor medio de -11.5 cm, para 225 BMs (pernos de nivelación) sobre tal estructura. Los errores de la estimación se ubicaron entre ±0.1 y ±6.3 cm, con valor medio de ±2.5 cm. Para el 64% de los BMs, los errores mejoraron los ±5 cm; sólo 3 determinaciones fueron consideradas outliers.

 

En la Figura 2 se aprecia que la técnica aplicada identifica a las secciones norte y central del dique en Lagunillas como las más afectadas por la deformación vertical de subsidencia; mientras que su sección sur se muestra estable. También en la figura, parte de la sección norte del dique costanero en Bachaquero se presenta afectada significativamente por el fenómeno.

Figura 2. Deformación vertical por subsidencia en pernos de nivelación (BMs) PDVSA, periodo 2012.0-2026.1, registrada por altimetría satelital SWOT sobre el dique costanero, región Lagunillas, COLM-Edo. Zulia [Acuña, 2026].

Aunque la técnica satelital SWOT objeto de esta nota no supera en precisión a la nivelación geodésica convencional PDVSA 'in-situ' [Acuña, 2026], es una muy útil alternativa de monitoreo, rápida y confiable, de uso a gran escala sin necesidad de operaciones de campo, que puede aportar información valiosa en tiempo ligeramente diferido (5-11 días después de las observaciones SWOT) sobre el comportamiento del fenómeno de subsidencia en la región y su impacto en la infraestructura petrolera; especialmente en periodos de tiempo donde no fue posible la realización de formales (clásicas) mediciones geodésicas locales de deformación. 

Más información en:

Acuña, G. (2026): Altimetría satelital SWOT para monitoreo de subsidencia del dique costanero, periodo 2012.0-2026.1, región Lagunillas, COLM-Edo.Zulia. Tópicos de Geodesia Geométrica. Mayo 18-24, 2026. Cátedra Geodesia Geométrica (GGenLUZ). Dpto. de Geodesia Superior. Esc. de Ingeniería Geodésica. Fac. de Ingeniería. Universidad del Zulia (LUZ). Maracaibo, Venezuela.

GGenLUZ produce MDTSubCOLM0426v1.0: modelo de deformación de terreno por subsidencia 2004-2026 para la COLM-Edo. Zulia

Investigación, Extensión, Datos/Soluciones/Modelos, Software – MDTSubCOLM0426v1.0 ; _intp.yab/.exe

Las cotas PDVSA por nivelación geodésica convencional 2004, 2007 y 2012 [PDVSA, 2013], [Acuña, 2013] de 2089 BMs de control de subsidencia en la COLM-Edo. Zulia, y aquellas determinadas recientemente en esos sitios con altimetría satelital SWOT, época 2026.1 [Acuña, 2026], fueron combinadas óptimamente por GGenLUZ usando -fast least-squares prediction FLSP- para generar un modelo regional de alta-resolución 450x450m que representa la velocidad media vertical de la deformación del terreno por subsidencia en el periodo 2004-2026 y sus errores estimados, ver Figuras 1 y 2.

Figura 1. Velocidades medias verticales de deformación del terreno por subsidencia (cm/año) en BMs PDVSA según MDTSubCOLM0426v1.0 [Acuña, 2026]: Modelo GGenLUZ de Deformación de Terreno por Subsidencia en la COLM-Edo. Zulia, periodo 2004-2026, versión 1.0, de resolución 450x450m, con filtraje gaussiano de 900 m 'full-width'. 

Figura 2. Error en velocidades medias verticales de deformación por subsidencia (cm/año) según MDTSubCOLM0426v1.0 [Acuña, 2026].

En términos de magnitudes de velocidades medias verticales de deformación, MDTSubCOLM0426v1.0 muestra tasas de subsidencia extremas entre -12 cm/año y +1 cm/año, con valor medio en la región de -1.2 cm/año. Para el periodo 2004-2026, el modelo señala la zona de mayor deformación en Lagunillas y otras significantes de deformación moderada en Tía Juana, Bachaquero y San Lorenzo (ver Figura 1).

Respecto a los errores de las velocidades asociados a MDTSubCOLM0426v1.0, estos se acumulan con valores elevados en la zona de mayor deformación de Lagunillas, donde se presentan en forma de incertidumbres de hasta ±3.7 cm/año para el referido periodo. En general, como valor medio para la región, las velocidades muestran errores estándar de ±0.5 cm/año (ver Figura 2).

MDTSubCOLM0426v1.0 y su software de aplicación MDTSubCOLM0426v1.0_intp.yab/.exe [Acuña, 2026] permitirán estimar, o más precisamente, propagar en el tiempo, con calidad de pocos centímetros y en la zona de validez del modelo (ver Figura 1), la cota de cualquier BM o estación de control geodésico a una época especifica de interés en el periodo 2004-2026 (o próximo a éste), conocida su posición ITRF/GRS80 y cota en otra época, previa o posterior.

Sobre el uso de MDTSubCOLM0426v1.0:

El modelo de deformación es aplicable a través de su software de interpolación MDTSubCOLM0426v1.0_intp.yab/.exe [Acuña, 2026], ver Figura 3. Éste posibilita, de nuevo, para BMs o estaciones de control geodésico en la zona del modelo, con cotas niveladas o alturas elipsoidales GNSS, respectivamente, propagar en el tiempo 'corrigiendo por subsidencia' esos valores de elevación, conocida la posición ITRF/GRS80 de tales sitios. A continuación se presenta un ejemplo típico de la aplicación de MDTSubCOLM0426v1.0 con MDTSubCOLM0426v1.0_intp.yab/.exe.

Figura 3. Software GGenLUZ MDTSubCOLM0426v1.0_intp.yab/.exe [Acuña, 2026].

El ejemplo refiere al BM PDVSA 212A de control de subsidencia (también estación GPS/GNSS pasiva de medición episódica bienal), ubicado en Lagunillas, muy cerca de la zona de máxima deformación en la COLM por tal fenómeno, ver Figura 4.

El BM 212A fue medido con GPS de punto-preciso (PPP) durante varias sesiones en la campaña de subsidencia PDVSA de 2007 [Higuera/Suárez, 2007]; sus coordenadas resultantes refieren al ITRF2020/GRS80 en la época 2007.9 presentando calidad centimétrica (mejor a ±1-cm) en las 3 componentes de la posición [Acuña, 2025].

A partir de su altura elipsoidal observada (h=-18.7563 m), la cota del BM para la época de medición 2007.9 (H=-4.2508 msnl) fue determinada con nivelación-GPS y el geoide nacional VGM25v1.0 (N=-14.5055 m) desde los BMs profundos PDVSA 9201, 9202 y 9203. Ahora bien, aquí se emplea MDTSubCOLM0426v1.0_intp.yab/.exe para estimar su valor (cota) en la época actual, i.e., 2026.1 (aprox.).

El resultado de la determinación se muestra en detalle en las Figuras 4 y 5. Considerando un periodo de transformación de 18.2 años, y una velocidad media vertical anual de deformación por subsidencia de -1.76±0.5 cm/año (dada por el modelo MDTSubCOLM0426v1.0), la cota del BM fue propagada corrigiéndola por subsidencia a la época de interés 2026.1. Su valor estimado final, H=-4.5715 msnl, obtenido con un error de ±0.0919 m, refleja el impacto de -32.07 cm de hundimiento por subsidencia del sitio, entre las épocas 2007.9 y 2026.1. El valor propagado de la cota del BM 212A con MDTSubCOLM0426v1.0_intp.yab/.exe difiere sólo -5.99 cm del obtenido con altimetría satelital SWOT y -1.68 cm respecto al valor sintético predicho por regresión lineal según cuadrados mínimos utilizando el registro histórico de la cota del BM. Tales diferencias permanecen dentro del rango de error de la estimación con el software empleado. Esto es un indicador veraz del muy satisfactorio nivel de confiabilidad de la herramienta computacional objeto de esta nota técnica.

Figura 4. Aplicación de MDTSubCOLM0426v1.0_intp.yab/.exe en la estación GPS/BM PDVSA 212A para propagar su cota y corregirla por subsidencia desde la época de medición 2007.9 a la actualidad (2026.1 aprox.).

Figura 5. Resultados completos (archivo ASCII .out) de MDTSubCOLM0426v1.0_intp.yab/.exe para la estación GPS/BM PDVSA 212A [Acuña, 2026].

Una versión DEMO de MDTSubCOLM0426v1.0_intp.yab/.exe está disponible en el siguiente enlace: 

https://mega.nz/file/YIEXTA4b#j09q4mqtYqMApArcEkBj9qhEPUv8_vQ9GIz7vq69seo

Más información en:

Acuña, G. (2026): MDTSubCOLM0426v1.0 - modelo GGenLUZ de deformación de terreno por subsidencia 2004-2026 para la COLM-Edo. Zulia. Tópicos de Geodesia Geométrica. Mayo 10-16, 2026. Cátedra Geodesia Geométrica (GGenLUZ). Dpto. de Geodesia Superior. Esc. de Ingeniería Geodésica. Fac. de Ingeniería. Universidad del Zulia (LUZ). Maracaibo, Venezuela.

Solución total para la red de nivelación subsidencia PDVSA utilizando altimetría SWOT, época 2026.1, COLM-Edo. Zulia

Investigación, Extensión – altimetría SWOT para red de nivelación PDVSA de subsidencia

La nota muestra resultados definitivos -solución total para 2089 BMs- de la red terrestre de subsidencia PDVSA en la COLM-Edo. Zulia, época 2026.1, obtenidos por GGenLUZ a partir del uso de altimetría satelital SWOT y de su software de procesamiento elevCOLMswot_v3.0.yab/.exe [Acuña, 2026], ver Figura 1. 

Siguiendo el procedimiento descrito en publicación anterior de este blog (i.e., https://ggenluz.blogspot.com/2026/04/determinacion-de-cotas-de-terreno-en.html), han sido determinadas por GGenLUZ cotas de terreno época 2026.1 para la totalidad (2089) de los BMs PDVSA de control de subsidencia en la COLM-Edo. Zulia. Esto representa una 'solución-total' por altimetría SWOT para la señalada red. La estimación utilizó 9608132 observaciones WSE_raster registradas por altimetría satelital SAR-interferométrica SWOT en el periodo 01.01.2026-09.04.2026, siendo procesadas y analizadas con el software elevCOLMswot_v3.0.yab/.exe [Acuña, 2026]. Este cálculo empleó 522 horas de intenso trabajo computacional, repartidas durante 7.5 días en 3 laptops con procesador 11th Gen Intel(R) Core(TM) i7-11370H@3.30Hz.

 

Las cotas época 2026.1 para los 2089 BMs de subsidencia fueron determinadas como referidas al nivel medio PDVSA del Lago de Maracaibo de 1926, realizado por el geoide nacional VGM25v1.0 [Acuña, 2025] convenientemente ajustado a la región utilizando los BMs profundos PDVSA 9201, 9202 y 9203.

 

La Figura 1 muestra las alturas de terreno estimadas en los BMs PDVSA, regiones Tía Juana, Lagunillas y Bachaquero, 'drapeadas' mediante la técnica least-squares draping [Acuña, 2004] sobre un muy detallado DTM regional de ultra-alta-resolución 30x30m generado con datos SRTMV3 [https://www.earthdata.nasa.gov/data/instruments/srtm]. Claramente se observa en las imágenes la deformación del terreno por la subsidencia local para la época 2026.1, con cotas negativas de hasta -12 metros aprox. en Lagunillas, 'por debajo del nivel medio del Lago'. Las posiciones geodésicas de los BMs PDVSA en ITRF/GRS80 aparecen indicadas en la figura como puntos de color blanco.

 

En  toda la red de subsidencia PDVSA, el cálculo con altimetría SWOT arrojó valores de altura para los 2089 BMs entre 99.482 m y -12.002 m, con valor medio de 11.037 m, y las mayores magnitudes de deformación en Lagunillas. Sólo 13 determinaciones de altura con SWOT (0.6% del total) fueron consideradas insuficientemente precisas (outliers), y por tanto descartadas por sus elevados errores o anormal magnitud. Ver resultados en Figura 4.

 

(1a)

(1b)

(1c)

Figura 1. Elevaciones de terreno H (m.s.n.l, época 2026.1), respecto al geoide nacional VGM25v1.0 [Acuña, 2025], obtenidas del procesamiento de mediciones altimétricas satelitales SWOT con elevCOLMswot_v3.0.yab/.exe [Acuña, 2026], sobre 2089 BMs de la red terrestre de subsidencia PDVSA en la COLM-Edo. Zulia, regiones Tía Juana (1a), Lagunillas (1b) y Bachaquero (1c). Las imágenes muestran las elevaciones SWOT de los BMs 'drapeadas' mediante least-squares draping sobre un muy detallado DTM regional de ultra-alta-resolución 30x30m generado con datos SRTMV3 [https://www.earthdata.nasa.gov/data/instruments/srtm]. 

La distribución geográfica de los errores estimados en las cotas SWOT por elevCOLMswot_v3.0.yab/.exe para los BMs PDVSA, regiones Tía Juana, Lagunillas y Bachaquero, se presenta en la Figura 2.

 

El error medio de la estimación de cotas SWOT para toda la red fue sólo ±3.9 cm, variando entre ±0.1 cm y ±20.1 cm. En el 63% de los BMs el error de la estimación de cotas fue menor a ±4 cm, en el 94% menor a ±10 cm, mientras que en el restante 6% de los casos, los errores se ubicaron entre ±10 y ±20 cm. Estos valores son considerados aquí como excelentes resultados si se tiene en cuenta que la calidad media de las mediciones altimétricas SWOT (WSE) es ±10 cm (1-sigma) sobre zonas continentales.

(2a)

(2b)

(2c)

Figura 2. Error estimado por elevCOLMswot_v3.0.yab/.exe [Acuña, 2026] para elevaciones de terreno H (WSE_SWOT) en 2089 BMs PDVSA de control de subsidencia, época 2026.1, COLM-Edo. Zulia, regiones Tía Juana (2a), Lagunillas (2b) y Bachaquero (2c).

Utilizando las cotas SWOT época 2026.1 estimadas por altimetría para los BMs [Acuña, 2026], y sus correspondientes alturas por nivelación geodésica convencional PDVSA época 2012.0 [PDVSA, 2013], fueron calculados valores de deformación por subsidencia para tales sitios, correspondientes al periodo 2012.0-2026.1. La Figura 3 muestra la representación gráfica de esas cantidades en las regiones Tía Juana, Lagunillas y Bachaquero. Allí, en los últimos 14 años, la deformación por subsidencia según la altimetría SWOT alcanzó -1.7 m, aprox., en el foco máximo del fenómeno. En la figura, la deformación por subsidencia aparece 'suavizada' mediante filtraje gaussiano de 900 m ('full-width').

 

En toda la red, se verificaron valores de hundimiento por subsidencia en 1952 BMs (94% del total), variando entre -1.718 m y 0.182 m, con un valor medio de -0.114 m. 

 

En 124 BMs (6% del total) se obtuvieron deformaciones positivas (elevaciones o ascensos) cuyas causas serán motivo de estudio en futuras investigaciones. Se presume que este comportamiento se deba, en principio, a errores de medición combinados entre las determinaciones altimétrica y de la nivelación convencional.

(3a)

(3b)

(3c)

Figura 3. Deformación de terreno por subsidencia (altimetría SWOT 2026.1 'minus' nivelación geodésica convencional PDVSA 2012.0) en 2089 BMs, periodo 2012.0-2026.1, COLM-Edo. Zulia, regiones Tía Juana (3a), Lagunillas (3b) y Bachaquero (3c); con filtraje gaussiano de 900 m ('full-width').

Resultados seleccionados de la 'solución-total' por altimetría SWOT para la red de subsidencia PDVSA, obtenidos con elevCOLMswot_v3.0.yab/.exe, se listan en la Figura 4. Ahí, los datos y resultados SWOT más representativos para cada uno de los 2089 BMs aparecen indicados, aunque limitados en precisión.

 

El archivo completo de resultados elevCOLMswot_v3.0.yab/.exe está disponible en: https://mega.nz/file/gUMUmDRb#a0bY-60kSgDLyAz1GtNuFUIF5m1m5QDycsM5D2Xnirs.

 

.

.

.

Figura 4. Solución-total GGenLUZ, época 2026.1, para la red terrestre de nivelación PDVSA de control de subsidencia COLM-Edo. Zulia, de 2089 BMs, según altimetría SWOT [Acuña, 2026]. Resumen seleccionado de resultados elevCOLMswot_v3.0.yab/.exe.

Un ejemplo particular de la estimación de la cota de BMs 2026.1 y de su deformación por subsidencia 2012.0-2026.1 a través de la altimetría SWOT, como aplicada a la totalidad de la red de subsidencia PDVSA, se presenta en la Figura 5. El caso refiere al BM 1319 lugar donde la deformación por subsidencia en el periodo de estudio alcanzó el máximo valor (cota= -12.002 ± 0.058 msnl ; subsidencia= -1.718 m). 

Figura 5. Estimación de subsidencia 2012.0-2026.1 en el BM PDVSA (1319) de mayor deformación de terreno en Lagunillas, COLM-Edo. Zulia, según altimetría satelital 'interferométrica-SAR' SWOT. Resultados del software elevCOLMswot_v3.0.yab/.exe [Acuña, 2026].

Consideraciones sobre los resultados obtenidos:

 

i) Tal como ha sido establecido en publicación anterior, los resultados de la presente estimación de subsidencia en la COLM-Edo. Zulia con altimetría SWOT 'no son predicciones de un software de simulación', 'ni valores generados con IA', por el contrario, son determinaciones precisas -óptimas- basadas en mediciones geodésicas reales, masivas y de alta exactitud.

 

ii) Un error medio de ±4 cm (y menores a ±10 cm para el 94% de los 2089 BMs considerados) en la estimación de cotas 2026.1 para la red PDVSA de subsidencia mediante altimetría SWOT, es un 'muy satisfactorio' (excelente) resultado. La validez de tal afirmación se basa en que esa calidad de alturas se logra utilizando mediciones realizadas desde una plataforma espacial (satélite) orbitando a 891 km sobre la Tierra, midiendo la altura del terreno con errores medios de ±10 cm, sin ningún dispositivo o sensor en superficie, y sin tener que ejecutar labores de campo. Además, las observaciones empleadas son de libre acceso, y su procesamiento y resultados se realiza y obtienen de forma remota, con rapidez de pocos días, y aplicación posible para cualquier fecha (y lugar en la región) entre 2023.5 y la actualidad. 

 

iii) Es claro que el procedimiento de estimación de cotas de BMs con altimetría SWOT 'no supera' la calidad de ±1-3 cm de la nivelación geodésica convencional PDVSA en la COLM, que se mantiene como la técnica de mayor exactitud en la derivación de cotas para el monitoreo regional de subsidencia. Sin embargo, la altimetría SWOT ofrece resultados también con calidad geodésica (pocos cm), precisos y útiles a mayor escala, rápidos (disponibles en pocos días), sin necesidad de operaciones de campo, y a un costo que representa una pequeña fracción (aprox. un 8%) del requerido por la muy costosa, laboriosa y tardía (en meses) nivelación geodésica convencional.

 

Más información en:

Acuña, G. (2026): Solución total para la red de nivelación subsidencia PDVSA utilizando altimetría SWOT, época 2026.1, COLM-Edo. Zulia. Tópicos de Geodesia Geométrica. Mayo 3-9, 2026. Cátedra Geodesia Geométrica (GGenLUZ). Dpto. de Geodesia Superior. Esc. de Ingeniería Geodésica. Fac. de Ingeniería. Universidad del Zulia (LUZ). Maracaibo, Venezuela.

Curso de extensión GGenLUZ: Determinación de alturas de terreno (cotas) en BMs venezolanos utilizando altimetría satelital SAR-interferométrica SWOT

Cursos/Talleres, Extensión – Curso de extensión 'on-line' GGenLUZ

GGenLUZ invita a estudiantes y profesionales de la ingeniería geodésica, geomática y geociencias afines, a participar en el curso de extensión, versión digital, titulado "Determinación de alturas de terreno (cotas) en BMs venezolanos utilizando altimetría satelital SAR-interferométrica SWOT", que estará disponible para descarga desde este website a partir del lunes 4 de mayo de 2026.

Desde agosto 2023, datos altimétricos operacionales de la misión satelital SWOT – Surface Water and Ocean Topography [https://swot.jpl.nasa.gov/] de NASA/CNES, están disponibles libres, y sistemáticamente, accesibles menos de una semana después de su adquisición. En principio, tales datos se distribuyen en forma de valores corregidos de alturas WSE (referidas al geoide global EGM2008) para océanos y cuerpos de agua continentales (ríos, lagos, reservorios, tierras inundables y acumulaciones temporales de agua superficial). Esos valores de altura los registra SWOT a lo largo de una franja terrestre de 120-km de ancho, con una elevada densidad de observaciones de 2x2 km a 100x100 m, cuya huella (paso altimétrico) se repite en ciclos exactos de 21 días, y que debido a su porcentaje de solapamiento asegura la re-observación del satélite de sitios terrestres cada 11 días en promedio (2 veces cada 21 días). En regiones oceánicas, las alturas de superficie registradas por los altímetros radar-nadiral clase-Jason e interferométrico-SAR KaRIn de SWOT tienen una calidad media de ±3 cm, mientras que en zonas continentales la calidad de la misión es ±10 cm para aéreas cubiertas de agua, sin vegetación, de al menos 1-km². Seleccionando mediciones altimétricas SWOT, tipo raster, de 100x100 m de resolución espacial, colectadas en modo interferométrico, de aguas superficiales temporales localizadas en el entorno (p.ej., 0.1 a 0.3 km de radio) de estaciones de nivelación geodésica posicionadas con GNSS, es posible determinar en esos sitios correspondientes alturas físicas H (cotas o elevaciones de terreno) con alta precisión, mediante predicción por cuadrados mínimos con funciones de covarianza espacio-tiempo y el geoide nacional de alta-resolución VGM25v1.0 [Acuña, 2025]. ... En este curso (de 4h académicas) se ofrecen detalles de la misión SWOT, sus principales aplicaciones y productos, además de instruir a los participantes en el procedimiento de la determinación geodésica de cotas de terreno con datos altimétricos SWOT en Venezuela, técnica contenida y operativizada a través del software GGenLUZ elevCOLMswot_2.0 [Acuña, 2026] desarrollado inicialmente para la estimación de alturas SWOT en BMs de control de subsidencia PDVSA en la COLM-Edo. Zulia.

Para más información sobre el curso (p.ej., contenido, forma de pago, detalles para descarga, consultas, etc.) contacte al Prof. Gustavo Acuña a través del email gacuna@fing.luz.edu.ve, o del número WhatsApp +58-412-42.71.579.