viernes, 18 de noviembre de 2022

Expresando en ITRF2020 las densificaciones nacionales de 3er. orden de SIRGAS-REGVEN

Investigación, Extensión - Densificaciones REGVEN en ITRF2020

Para actualizar y expresar en la reciente solución ITRF2020, las coordenadas de 742 vértices geodésicos de densificación REGVEN de orden C (±5cm), ver Figura 1, fueron transformadas desde la realización SIRGAS-REGVEN(1995)ITRF94,1995.4/GRS80 a SIRGAS-REGVEN(2022)ITRF2020,2022.5/GRS80.

Figura 1. Estaciones de densificación REGVEN orden C (742) vs. estaciones primarias
REGVEN (176) y REMOS (31) de órdenes A y B.

Este proceso ya fue ensayado en las estaciones SIRGAS, REGVEN y REMOS (órdenes A y B) de Venezuela, y descrito en una publicación anterior (Acuña, 2022a; https://ggenluz.blogspot.com/2022/05/adaptando-la-realizacion-2015-de-sirgas.html).

Aquí, las coordenadas objeto de transformación corresponden a posiciones 3D validadas por el Instituto Geográfico de Venezuela Simón Bolívar (IGVSB; http://www.igvsb.gob.ve/) -disponibles desde su portal web- para estaciones REGVEN de diversos proyectos de densificación ejecutados a nivel nacional. En el proceso se utilizó además el modelo venezolano de velocidades geodésicas VVM20, los parámetros de transformación entre soluciones ITRF del IERS, y el software GGenLUZ para transformaciones ITRF, Trn_ITRF,v4.0.

La transformación SR1995 -> SR2022 introdujo cambios medios en las coordenadas de las 742 estaciones de 37.3cm, 17.1cm y 1.2cm (rms) en latitud y longitud geodésicas, y altura elipsoidal, resp. En posición, más del 95% de la magnitud de las variaciones es causada por 27 años de desplazamiento acumulado en cada estación entre las épocas 1995.4 y 2022.5 debido al efecto de la tectónica regional de placas.

Pruebas de validación en estaciones REGVEN con coordenadas SR2022 “observadas” (obtenidas de mediciones GNSS directas), indicaron incertidumbres medias de 4.1cm, 3.2cm y 2.8cm (rms) en las coordenadas φλh SR2022 de las restantes estaciones que fueron derivadas por transformación. Estos valores cuantifican la calidad general del proceso de transformación aplicado aquí, ubicándose dentro del umbral de precisión tolerado en Venezuela para estaciones REGVEN de orden C (±5cm). 

Adicionalmente, para expresar las coordenadas iniciales SR1995 de los vértices de densificación REGVEN en las realizaciones 2000 y 2015 de SIRGAS-REGVEN, transformaciones respectivas fueron también incluidas.

La información generada en este trabajo para los vértices de densificación REGVEN será incluida en la nueva versión de la base de datos GGenLUZ para posiciones de estaciones geodésicas SIRGAS-REGVEN, SR95001522P_v2.0 (Acuña, 2022b; https://ggenluz.blogspot.com/2022/06/nuevo-software-ggenluz-sr950015p-base.html).

Más detalles sobre esta nota son dados en la siguiente publicación:

Acuña, G. (2022): "Expresando en ITRF2020 las densificaciones nacionales de 3er. orden de SIRGAS-REGVEN". Tópicos de Geodesia Geométrica. Noviembre 13-19, 2022. Cátedra Geodesia Geométrica (GGenLUZ). Lab. de Geodesia Física y Satelital. Dpto. de Geodesia Superior. Facultad de Ingeniería. Universidad del Zulia. Maracaibo, Venezuela.


sábado, 22 de octubre de 2022

Transportando alturas físicas (cotas) de vértices geodésicos sobre largas distancias en Venezuela

Investigación, Extensión - ... cont. la validación de VGM19v3.0

En geodesia, una de las principales aplicaciones de la determinación regional de alta-resolución del geoide, es producir detallados modelos de esa superficie equipotencial para su uso combinado con el posicionamiento preciso GNSS y entonces poder generar alturas físicas de terreno (cotas) vía el procedimiento diferencial de la nivelación-GNSS, o más comúnmente, vía la aplicación absoluta de la conocida fórmula de la geodesia geométrica: h = H + N (donde h es la altura elipsoidal de la estación proveniente del posicionamiento GNSS, H la altura física de terreno -a derivar- sobre el geoide, y N la ondulación del geoide proveniente de un modelo regional de alta-resolución y precisión).

La determinación de elevaciones de terreno de esta forma, es una alternativa muy atractiva en términos de costos operativos, tiempo de ejecución y calidad geodésica, a la en extremo exigente técnica de la nivelación geodésica convencional. Ésta última, es uno de los procedimientos de medición mas costosos en tiempo, dinero y trabajo en la ingeniería geodésica, especialmente cuando se desarrolla sobre largas distancias (i.e., decenas a cientos de kms.). Por ejemplo, sobre un trayecto carretero de aprox. 100 km, transportar una cota de terreno de un punto A a otro punto B, mediante nivelación geodésica convencional de 1er. orden, podría tomar ca. de  50 días de trabajo de campo para un grupo de medición (con todos los costos que esto implica); mientras que esa misma tarea empleando nivelación-GNSS podría tomar un máx. de 3 días de trabajo campo, y sobre todo sin los riesgos que implica el trabajo diario en carretera asociado a la nivelación convencional.

Figura 1. VGM19v3.0 - Venezuelan Geoid Model 2019, versión 3.0 (Acuña, 2021),
geoide gravimétrico ajustado a 383 estaciones nacionales GPS/nivelación.

En esta nota técnica se evalúa para Venezuela el potencial de la nivelación-GNSS cuando se aplica sobre largas distancias (entre 400 a 900 km aprox.) y se utilizan las posiciones de vértices geodésicos SIRGAS-REGVEN de orden A medidos durante la campaña GNSS nacional REGVEN2015 (Acuña et al., 2017), sus respectivas alturas niveladas, y la última versión (gravimétrica) del modelo geoidal venezolano VGM19v3.0 (Acuña, 2021), ver Figura 1.

A manera de ejemplo, en la Figuras 2 y 3 se muestran las estaciones LA_GUAIRA (GUA2) y LA_CANOA (CANO) midiendo durante la campaña REGVEN2015. Tales vértices son estaciones SIRGAS-REGVEN de orden A (±1cm), la primera, es excéntrica al BM principal del mareógrafo LaGuaira (costa centro-norte de Venezuela), lugar donde se define el datum vertical del sistema de alturas del país (LGVD1962); y la segunda, es el vértice geodésico que realiza el punto datum horizontal del sistema nacional de control geodésico convencional PSAD56 (Datum Provisional para Suramérica de 1956, ó LaCanoa/Hayford para Venezuela), en la zona centro-oriental del país. Además de sus respectivas coordenadas LAT,LON,ALT.ELIP conocidas en la realización SIRGAS-REGVEN(2015), ambas estaciones tienen asignados valores de ondulación del geoide provenientes del modelo VGM19v3.0, y cotas niveladas por el IGVSB que servirán para ejecutar el transporte de la altura física GUA2 -> CANO, y validar su resultado.

Figura 2. Estación SIRGAS-REGVEN orden A LA_GUAIRA (GUA2)
 durante la campaña REGVEN2015.

Figura 3. Estación SIRGAS-REGVEN orden A LA_CANOA (CANO)
 durante la campaña REGVEN2015.

La posición relativa de las estaciones GUA2 y CANO, y del resto de los vértices SIRGAS-REGVEN orden A en Venezuela, a saber, AMUAY (AMUA), MARACAIBO (MARA), SAN_ANTONIO (SNTN), AGUA_LINDA (AGUA), STA_ELENA (SELE), KAMÁ (KAMA) y CARÚPANO (CARU), se muestra en la Figura 4. El vértice JUNQUITO (JUNQ) no se considera en esta nota por su cercanía (aprox. 22 km) a la estación GUA2. La Figura 4 muestra también las conexiones entre la estación GUA2 y cada uno de los vértices antes señalados, representando el transporte de la cota de esa estación hacia cada vertice cuyos resultados se presentan en la Tabla 1. 

Figura 4. Conexiones entre estaciones seleccionadas SIRGAS-REGVEN(2015) orden A
para el transporte de alturas físicas (cotas) en Venezuela por nivelación-GNSS.

Para el transporte GUA2 -> CANO, su perfil de terreno de 406 km de longitud (ver Figura 5) indica diferencias de altura topográficas de hasta 2300 m a lo largo de la trayectoria, mostrando en la región alturas elipsoidales sistemáticamente inferiores a las elevaciones (alturas físicas) de terreno. Esto se debe a que en la zona el geoide VGM19 se ubica 22 m -en promedio- por debajo del elipsoide GRS80, una consecuencia del fuerte gradiente negativo en dirección noreste que observa el geoide en territorio venezolano (ver Figura 1). La aplicación de la nivelación-GNSS en este caso arrojó como resultado un error de aprox. 10 cm en la determinación de la altura física de la estación CANO desde GUA2, lo que representa un error relativo de sólo 0.3 ppm en distancia (ver Tabla 1). Este valor es aproximadamente el mismo al obtenido en promedio para el resto de los transportes ensayados en esta nota, ver Tabla 1. 

Figura 5. Perfil de terreno correspondiente al transporte de cota entre las estaciones
GUA2 y CANO mediante nivelación-GNSS.

Tabla 1. Resultados del transporte de la cota de la estación LA_GUAIRA (GUA2) hacia 8
vértices SIRGAS-REGVEN orden A en Venezuela, sobre largas distancias,
 empleando la técnica nivelación-GNSS.

Continuando con las pruebas de validación del más reciente modelo nacional del geoide, y de acuerdo a los resultados aquí señalados, es posible establecer que para un rango de distancia inter-estación entre 400 y 900 km, el procedimiento de nivelación-GNSS empleando el modelo geoidal VGM19v3.0, ofrece una calidad relativa de ± 0.3 ppm (i.e., ± 0.3 mm/km) en la determinación de alturas físicas (cotas en m.s.n.m, datum LGVD1962) de vértices geodésicos en Venezuela. Esto se corresponde cercanamente con las precisiones que podría ofrecer la nivelación geodésica convencional de 2do. orden, aunque necesitando solo una pequeña fracción del costo operativo que el último procedimiento implica.

Para más detalles sobre esta nota, puede consultar la siguiente publicación:

Acuña, G. (2022): "Transportando alturas físicas (cotas) de vértices geodésicos sobre largas distancias en Venezuela". Tópicos de Geodesia Geométrica. Octubre 23-29, 2022. Cátedra Geodesia Geométrica (GGenLUZ). Lab. de Geodesia Física y Satelital. Dpto. de Geodesia Superior. Facultad de Ingeniería. Universidad del Zulia. Maracaibo, Venezuela.


martes, 18 de octubre de 2022

PTReGVen3_v1.0: Software GGenLUZ para la transformación de coordenadas geodésicas entre las realizaciones SIRGAS-REGVEN 1995, 2000, 2015 y 2022

InvestigaciónExtensiónSoftware - PTReGVen3_v1.0.yab/.exe

PTReGVen3_v1.0 (Acuña, 2022a) es un software desarrollado durante Octubre 2022 por la cátedra Geodesia Geométrica de LUZ (GGenLUZ) para ejecutar transformaciones de coordenadas geodésicas entre las realizaciones 1995, 2000, 2015 y 2022 del sistema de referencia geodésico SIRGAS-REGVEN (SR) de Venezuela, ver Figura 1.

Figura 1. Software PTReGVen3_v1.0.yab/.exe.

Escrito en YaBasic_v2.90.2 para su ejecución en Linux/Unix o Windows, PTReGVen3_v1.0 representa la 3era. generación de la conocida herramienta computacional GGenLUZ/LGFS-LUZ PTReGVen (Acuña et al., 1998), PTReGVen2 (Acuña, 2011), inicialmente concebida para operativizar la transformación de coordenadas geodésicas entre los datum PSAD56 [LaCanoa/Hayford] y SIRGAS-REGVEN(1995) [ITRF94, 1995.4/GRS80] en Venezuela.

PTReGVen3_v1.0 aplica con extrema rigurosidad matemática los últimos parámetros de transformación del datum (PTD) calculados en el país que relacionan en todas las combinaciones posibles las 4 realizaciones SR hoy disponibles, ver Acuña (2022b). 

El software determina y aplica los cambios de coordenadas utilizando el modelo de similaridad clásico Molodensky-Badekas (Badekas, 1969) para la transformación de coordenadas geodésicas 3D. 

PTReGVen_v1.0 trabaja primero aproximando las diferencias sistemáticas entre las realizaciones SR dadas por los PTD, para luego modelar (y adicionar) en las estaciones a transformar, mediante colocación por cuadrados mínimos -LSC-, las diferencias residuales entre las soluciones SR previamente estimadas en 207 vértices geodésicos del país al momento de calcular los PTD.

PTReGVen3_v1.0 transforma las coordenadas de puntos simples o lista de n-estaciones, ingresando las coordenadas en una realización SR específica para obtener las coordenadas transformadas en las 3 realizaciones restantes.

Dependiendo del tipo de transformación ejecutada por PTReGVen3_v1.0, el software producirá coordenadas transformadas con calidad de ± 1-3 cm.

Una versión DEMO de PTReGVen3_v1.0 está disponible para los usuarios nacionales de GGenLUZ. Para descargarla contacte a su autor, Prof. Gustavo Acuña, a través del email gacuna@fing.luz.edu.ve.

Para más detalles sobre PTReGVen3_v1.0, puede consultar la siguiente publicación:

Acuña, G. (2022): PTReGVen3_v1.0: Software para la transformación de coordenadas geodésicas entre las realizaciones SIRGAS-REGVEN 1995, 2000, 2015 y 2022. Tópicos de Geodesia Geométrica. Octubre 16-22, 2022. Cátedra Geodesia Geométrica (GGenLUZ). Laboratorio de Geodesia Física y Satelital. Dpto. de Geodesia Superior. Facultad de Ingeniería. Universidad del Zulia. Maracaibo, Venezuela.


miércoles, 12 de octubre de 2022

Parámetros de transformación del datum (PTD) entre las realizaciones 1995, 2000, 2015 y 2022 del sistema de referencia geodésico nacional SIRGAS-REGVEN

Investigación, Extensión - PTD_SIRGAS-REGVEN

Utilizando 207 vértices geodésicos SIRGAS/REGVEN/REMOS en Venezuela, con coordenadas conocidas en las realizaciones SIRGAS-REGVEN de 1995 (Hérnandez et al., 2000), 2000 (Hernández et al., 2003), 2015 (Acuña, 2017) y 2022 (Acuña, 2022), se determinaron 6 grupos de 7+3 parámetros de transformación del datum (PTD) para relacionar las soluciones anteriores en todas las combinaciones posibles y habilitar una transformación simple entre sus coordenadas asociadas.

La Figura 1 muestra las posiciones geográficas de las 207 estaciones utilizadas para la determinación de los grupos de PTD; 176 vértices corresponden a estaciones pasivas tipo SIRGAS/REGVEN, y 31 a estaciones activas tipo REMOS.

Figura 1. Red nacional de estaciones SIRGAS/REGVEN/REMOS.

Una típica estacion REGVEN, con coordenadas conocidas en las 4 realizaciones hasta hoy disponibles del sistema de control geodésico venezolano SIRGAS-REGVEN, se muestra en la Figura 2.

Figura 2. Estación REGVEN Chiguará (CHIG), Edo. Mérida, durante la campaña REGVEN2015.

Los grupos de PTD fueron determinados mediante cálculo y ajuste por cuadrados mínimos, empleando el clásico modelo de similaridad Molodensky-Badekas (Badekas, 1969) como implementado en el software GGenLUZ TransDAT_v4.0 (Acuña, 1997-2022). El proceso de estimación utilizó las coordenadas geodésicas curvilíneas (finales observadas + sintéticas) de los vértices en sus respectivas soluciones además de sus estimaciones de error.

La Tabla 1 presenta los 6 grupos de 7+3 PTD calculados. Se indica para cada grupo, a) las realizaciones SIRGAS-REGVEN inicial y final consideradas, b) el número de puntos comunes utilizados en la estimación, c) 3 parámetros de translación (DZ,DY,DZ), d) 3 parámetros de rotación (EX,EY,EZ), e) 1 parámetro de escala (DM), f) las coordenadas geocéntricas cartesianas del centroide para el área de cálculo (XM,YM,ZM), y g) el error medio (s0) de la transformación Molodensky-Badekas en los puntos comunes cuando se utiliza tal grupo de parámetros.

Tabla 1. PTD entre realizaciones SIRGAS-REGVEN.

Al cambiar de una realización SIRGAS-REGVEN a otra, cada grupo de PTD aproxima las variaciones de coordenadas que producen, por un lado, el cambio en la solución global ITRF a la cual refieren las realizaciones, y por otro, los desplazamientos tiempo-dependientes (velocidades) en las estaciones por efecto de la deriva de placas téctonicas.

Así, por ejemplo, los parámetros que relacionan las realizaciones SIRGAS-REGVEN de los años 1995 (SR1995) y 2015 (SR2015) aproximan las variaciones en las coordenadas geodésicas de las estaciones REGVEN que se producen al cambiar del marco global ITRF1994 al ITRF2014, y al considerar los desplazamientos en las estaciones por deriva de placas entre las épocas 1995.4 y 2015.5.

Observando la magnitud de los valores de este 2do. grupo de parámetros en la Tabla 1, es posible identificar que la transformación SR1995 -> SR2015 introduce variaciones medias en las coordenadas iniciales de aprox. 28cm y 6cm en latitud y longitud geodésicas (ver parámetros DZ y DY), y -3cm en altura elipsoidal (ver parámetro DX), respectivamente.

Finalmente cabe señalar, los grupos de PTD aquí calculados fueron determinados con la finalidad de facilitar la actualización "sencilla, rápida, directa y masiva" de productos geodésicos y cartográficos nacionales, desde la realización SIRGAS-REGVEN de 1995 a una versión más reciente y mejorada, por ejemplo, aquella del año 2015 ó 2022; ésto sin necesidad del uso combinado de grupos de hasta 14 parámetros de transformación entre soluciones ITRF, modelos regionales de velocidad de placas tectónicas, y complejos modelos matemáticos para transformaciones de similaridad.

Los PTD listados en la Tabla 1 podrán ser aplicados con rigurosidad mediante el software GGenLUZ PTReGVen3_v1.0 (Acuña, 2022),  próximamente disponible a través de este blog.

Para más detalles sobre esta nota técnica, puede consultar la siguiente publicación:

Acuña, G. (2022): Parámetros de transformación del datum entre las realizaciones 1995, 2000, 2015 y 2022 del sistema de referencia geodésico nacional SIRGAS-REGVEN de Venezuela. Tópicos de Geodesia Geométrica. Octubre 9-15, 2022. Cátedra Geodesia Geométrica (GGenLUZ). Laboratorio de Geodesia Física y Satelital. Dpto. de Geodesia Superior. Facultad de Ingeniería. Universidad del Zulia. Maracaibo, Venezuela.




OBSERVACIÓN: la anterior información (i.e., PTD entre realizaciones SIRGAS-REGVEN) se ofrece libre a la comunidad geodésico/cartográfica del país pero sin ningún compromiso de garantía. Su autor, Prof. G. Acuña, y la cátedra GGenLUZ, no se hacen responsables por la aplicación, adecuada o nó, que los usuarios hagan de tales datos, ni de las consecuencias que sus resultados puedan generar.

lunes, 18 de julio de 2022

REGVEN2015v3.0: solución 3.0 de GGenLUZ para la red GNSS nacional de Venezuela, campaña 2015

InvestigaciónExtensión - Solución REGVEN2015v3.0

Aprovechando la disponibilidad de las nuevas versiones de los software de procesamiento GNSS GGenLUZ/GNSS-PPP (Acuña, 2022) y GipsyX (Bertiger et al., 2020; https://gipsy-oasis.jpl.nasa.gov/) en enero 2022, la campaña 2015 de REGVEN fue reprocesada por GGenLUZ, dando como resultado la solución REGVEN2015v3.0 de esa red GNSS nacional.

REGVEN2015 (ver Figura 1) representa la segunda campaña GNSS de remedición de la REd Geocéntrica de VENenezuela (Acuña et al., 2017; Acuña, 2019). Llevada a cabo entre junio 08 y julio 14 de 2015, REGVEN2015 fue coordinada y ejecutada por el Instituto Geográfico de Venezuela Simón Bolívar (IGVSB), contando además con la participación de instituciones nacionales públicas y privadas, p.ej., GGenLUZ/LGFS como centro de procesamiento y análisis, y PDVSA(occidente/oriente), INTEVEP, DIGECAFA, MPPRE, USB, UDO, LUZ, DelMonteOT, entre otras, en la realización de las observaciones. REGVEN2015 fue concebida para extender y actualizar a la solución global ITRF2014, el marco de referencia geodésico de Venezuela; además de posibilitar: a) la determinación de un modelo de deformación nacional para la red de control geodésico básico del país, b) la estimación de parámetros de transformación del datum entre REGVEN2015, PSAD56 y las soluciones REGVEN anteriores (1995 y 2000), y c) en conjunto con el modelo geoidal venezolano de ultra-alta-resolución 30x30m VGM17, construir (definir y realizar) un novedoso sistema de referencia vertical para la determinación de alturas físicas (normales) en el país. Durante REGVEN2015 se ocuparon en definitiva 145 estaciones pasivas tipo REGVEN1995-2000 (ver Figura 2), considerándose además 31 estaciones activas tipo REMOS (ver Figura 3). En términos generales, las estaciones pasivas fueron medidas durante 2 días en sesiones de 8 horas, empleando receptores geodésicos GPS y GNSS de doble/multi-frecuencia. En las estaciones activas REMOS, mediciones GNSS permanentes registradas durante 3 a 5 días fueron utilizadas. Más detalles en (Acuña et al., 2017).

Figura 1. Estaciones REGVEN y REMOS ocupadas durante la campaña REGVEN2015.

Figura 2. Medición GNSS en estación REGVEN (orden B) EL_RECREO (RECR), Edo. Zulia,
durante la campaña REGVEN2015.

Figura 3. Medición GNSS en estación REMOS (orden A) CARACAS (CRCS), Caracas D.C.,
durante la campaña REGVEN2015.

En 2017 y 2019, soluciones para REGVEN2015, versiones 1.0 (Acuña et al., 2017) y 2.0 (Acuña, 2019), resp., fueron generadas por GGenLUZ. Ambas soluciones son "solo-GPS", comprenden 172 estaciones (142 pasivas REGVEN y 30 activas REMOS) y refieren al ITRF2014, época 2015.5, elipsoide GRS80.

REGVEN2015v1.0 es una solución combinada "preliminar" que integra con COLUMBUSv4.5 varias soluciones independientes obtenidas con sistemas de procesamiento GNSS on-line vía web como CSRS-PPP, AUSPOS, APPS-GIPSY, OPUS, TrimbleRTX, SCOUT y GGenLUZ/GNSS-PPP. Las incertidumbres medias obtenidas (precisiones 1-sigma) en coordenadas latitud, longitud y altura elipsoidal para las posiciones de vértices REGVEN2015v1.0, órdenes A, B y C, fueron ± 2.7 mm, ± 3.5 mm y ± 8.0 mm, resp.

REGVEN2015v2.0 es una solución definitiva para la época de su estimación (marzo 2019); ésta fue obtenida a partir de un exclusivo procesamiento científico GPS-PPP con el software GIPSY-OASISv6.4. En esta versión 2.0, fueron corregidos algunos problemas de altura y tipo de antenas GNSS, e identificación de estaciones REGVEN reconstruidas presentes en la versión 1.0; además se logró incorporar días adicionales de medición en algunas estaciones permanentes REMOS. Como resultado, las incertidumbres medias obtenidas en coordenadas latitud, longitud y altura elipsoidal para las posiciones de vértices REGVEN2015v2.0 fueron ± 2.5 mm, ± 3.3 mm y ± 4.9 mm, resp. En la versión 2.0, el delicado modelaje troposférico y de mareas terrestres/oceánicas, y el uso de productos orbitales JPL que posibilitan la resolución de ambigüedades en el procesamiento GPS-PPP que aplica el software GIPSY-OASISv6.4, resultó en un significativo mejoramiento en la estimación (repetibilidad) de las alturas de las estaciones, casi por un factor 2, respecto a los resultados de la versión preliminar 1.0 de REGVEN2015.

En este trabajo, una nueva solución (final en la actualidad) fue obtenida del reprocesamiento de los datos REGVEN2015 empleando las recientes versiones 2.0 de GipsyX, y 3.0 de GGenLUZ/GNSS-PPP, como disponibles a comienzos de 2022. Esta solución, REGVEN2015v3.0, incluye ahora todos los datos disponibles GPS + GLONASS registrados en la referida campaña. Se adicionaron 4 nuevas estaciones para un total de 176 (145 REGVEN y 31 REMOS). En 136 estaciones se registraron datos GPS/GLN y en 40 datos solo GPS. Debido a la cercana correspondencia (a nivel del sub-centímetro) entre los resultados de ambos software de procesamiento GNSS para las posiciones 3D de estaciones REGVEN2015, REGVEN2015v3.0 fue calculada como una solución combinada, tambien referida al ITRF2014_2015.5/GRS80.

A diferencia de las realizaciones previas, REGVEN2015v3.0 es una solución basada en procesamiento GNSS-PPP estático que utiliza observaciones GPS+GLN de alta-resolución temporal cada 5 segundos, y productos orbitales precisos finales multi-constelación de 30 segundos de resolución. El procesamiento PPP con GipsyX_v2.0 emplea productos JPL exclusivos para GipsyX, mientras que GGenLUZ/GNSS-PPP_v3.0 trabaja con productos operacionales GFZ/CODE de estándar IGS.

En términos de resultados obtenidos, la solución combinada REGVEN2015v3.0 ha sido la de mayor calidad hasta ahora alcanzada; ésta mostró incertidumbres medias (1s) en coordenadas latitud, longitud y altura elipsoidal para vértices geodésicos de ± 2.2 mm, ± 2.8 mm y ± 3.9 mm, resp. Del total de 176 estaciones, 142 (80.5%) presentaron errores 3D inferiores a ± 1 cm, 26 (14.9%) < ± 2 cm, y 8 (4.6%) < ± 2.5 cm; ver Figura 4. En estaciones REMOS, las incertidumbres medias en latitud, longitud y altura elipsoidal fueron aun menores, a saber, ± 1.1 mm, ± 1.7 mm y ± 2.9 mm, resp.

Figura 4. Incertidumbres para posiciones 3D en solución GGenLUZ REGVEN2015v3.0.

Alturas físicas (elevaciones de terreno en m.s.n.m, datum La_Guaira) fueron estimadas para los vértices geodésicos REGVEN2015v3.0 mediante nivelación-GNSS empleando el geoide nacional VGM19v3.0. Comparaciones entre este tipo de alturas, y aquellas obtenidas por nivelación convencional en 55 estaciones REGVEN mostraron diferencias con RMS de ± 3.8 cm. Información geodésica complementaria como deflexiones de la vertical, gravedad total, ondulación del geoide, velocidades geodésicas, potencial de gravedad y número geopotencial, además de coordenadas geocéntricas cartesianas y planas globales proyectadas (UTM y LCC), fue también generada para cada vértice.

Considerando la excelente precisión de los resultados logrados con REGVEN2015v3.0, esta solución será la base de la próxima realización GGenLUZ para 2022 del sistema de referencia geodésico de Venezuela, a saber, SIRGAS-REGVEN(2022).

Más detalles sobre esta nota son dados en la siguiente publicación:

Acuña, G. (2022): "REGVEN2015v3.0: solución 3.0 de GGenLUZ para la red GNSS nacional de Venezuela, campaña 2015". Tópicos de Geodesia Geométrica. Julio 17-23, 2022. Cátedra Geodesia Geométrica (GGenLUZ). Lab. de Geodesia Física y Satelital. Dpto. de Geodesia Superior. Escuela de Ingeniería Geodésica. Facultad de Ingeniería. Universidad del Zulia. Maracaibo, Venezuela.



domingo, 3 de julio de 2022

Actualizando la determinación de la tasa de ascenso del nivel medio del mar en el Caribe venezolano

Investigación - Aumento del MSL en Venezuela (1992-2022)

El aumento del nivel medio global del mar y de la temperatura media terrestre son los principales indicadores del cambio climático.

Según el último informe del IPCC (Panel Intergubernamental sobre Cambio Climático; ONU, 2021), la temperatura registrada en el periodo 2015-2019 fue la más cálida observada de los periodos equivalentes anteriores, situándose 1.1°C por encima de los valores de la época pre-industrial (1850-1900), y 0.2°C por encima de los valores del periodo 2011-2015.

Por su parte, observaciones de satélites altimétricos desde la década de los años 90's evidencian que el nivel medio global del mar no solamente muestra un aumento sostenido durante los últimos 30 años, sino que además está acelerándose de +2.2 mm/año en 1993 a +5 mm/año hoy en día.

Ambos indicadores globales correlacionan positivamente con el aumento de la concentración de gases de efecto invernadero en la atmósfera, C02 principalmente, un 146% mayor a los niveles de la época pre-industrial, y aumentando aprox. un 2% anual en la actualidad.

Esta nota técnica trata sobre la última estimación GGenLUZ de la tasa de aumento del nivel del mar en el Caribe venezolano, periodo 1992-2022, como determinada en la costa central de país, para la posición del mareógrafo La Guaira, lugar donde se define el datum vertical de Venezuela (LGVD1962), origen del sistema nacional de alturas. Además, aquí se verifica si la referida tasa experimenta una aceleración local, y se proyecta su posible impacto futuro sobre la línea y tierras-bajas costeras de Venezuela.

La actual determinación GGenLUZ del aumento del nivel del mar en La Guaria se basa en observaciones satelitales multimisión de alturas corregidas de la superficie del mar instantáneas, registradas en el periodo 1985.5-2022.5, ver Figura 1. Se consideran las fases geodésicas y de repetición exacta de 12 misiones satelitales altimétricas, i.e., GEOSAT, ERS-1, TOPEX/Poseidon, ERS-2, GFO, JASON-1, ENVISAT, JASON-2, CRYOSAT-2, JASON-3 (ver Figura 2), SARAL/AltiKa y SENTINEL-3A. En total, cerca de 145500 mediciones satelitales altimétricas localizadas hasta 165 km alrededor del mareógrafo La Guaira son utilizadas.   

Figura 1. Observaciones satelitales altimétricas multimisión alrededor del mareógrafo La Guaira registradas en el periodo 1985.5-2022.5; tales mediciones son utilizadas por GGenLUZ en su determinación más reciente de la tasa de aumento del nivel del mar en el Caribe venezolano.


Figura 2. Satélite Jason-3 (https://www.jpl.nasa.gov/missions/jason-3), último hasta ahora de la serie TOPEX/Poseidon + Jason-1/2/3. Esta serie de satélites altimétricos NASA/CNES ha sido la más precisa y de mayor tiempo continuo de operación en el monitoreo global de los cuerpos de agua terrestres.

Los datos altimétricos mostrados en la Figura 1 corresponden a cantidades CorSSHs (alturas elipsoidales corregidas instantáneas de la superficie del mar) obtenidas de la base global de datos altimétricos RADS/DEOS (http://rads.tudelft.nl/rads/rads.shtml). Estas observaciones son tratadas con los últimos modelos para correcciones orbitales, instrumentales y por el medio de propagación, además de ser intercalibradas (mediante ajuste crossover) respecto a los datos de mayor calidad en órbitas precisas, i.e., la serie T/P+J1+J2+J3; esto permite el uso combinado y preciso de los datos multimisión. La incertidumbre media de las mediciones altimétricas fue ± 2.3 cm.

Luego de reducir de las observaciones CorSSHs la contribución de la superficie media del mar (MSSHs) proveniente del reciente modelo DTUMSS21 (Andersen et al, 2021), los valores de anomalías del nivel del mar (SLAs) obtenidos se utilizaron para generar por LSC en la locación del mareógrafo "virtual" La Guaira, valores medios zonales del nivel del mar (MSLs), ponderados en distancia y tiempo, con resolución temporal de 10 días. Para el periodo 1992-2022, un total de 10951 valores MSLs fueron estimados con una precisión media de ± 1.7 cm.

A la serie altimétrica satelital de valores MSLs se le ajustó por cuadrados mínimos un modelo de regresión lineal para estimar la tasa de ascenso local del nivel del mar, promedio para el periodo 1992-2022. El resultado fue +3.4 ± 0.4 mm/año, ver Figura 3.

Figura 3. Serie de tiempo 1992-2022 del nivel medio del mar (MSL) obtenida por altimetría satelital multimisión en la costa central de Venezuela (mareógrafo La Guaira), utilizando las observaciones indicadas en la Figura 1. La evolución en el tiempo del MSL en La Guaria (valores medios zonales cada 10 días) confirma una tasa de ascenso de +3.4 ± 0.4 mm/año durante los últimos 30 años de monitoreo satelital.

La tasa estimada aquí, +3.4 ± 0.4 mm/año con época media 2007.0, es significativamente mayor (+10%) a la determinada en 2011 por GGenLUZ, a saber, +3.1 ± 0.3 mm/año con época media 1997.0 (Acuña, 2011). Esto confirma que, en efecto, la tasa de aumento del nivel medio del mar en la costa central de Venezuela se incrementó (o aceleró) entre la última década del siglo pasado y la primera de este siglo, a razón de aprox. +1% anual, equivalente a, +0.03 mm/año.

Una visión más precisa de la aceleración del aumento del MSL en La Guaira aparece representada en la Figura 4.

Figura 4. Estimaciones para las 3 últimas décadas de la tasa de aumento del nivel medio del mar en el mareógrafo La Guaira, costa central de Venezuela, i.e., periodos 1992-2002 (+1.8 ± 0.2 mm/año), 2002-2012 (+2.9 ± 0.3 mm/año), y 2012-2022 ( +5.7 ± 0.4 mm/año).

Al comparar las tasas de aumento del nivel del mar en La Guaria obtenidas en periodos de tiempo similares correspondientes a las 3 últimas décadas (ver Figura 4), la aceleración del MSL se hace aun más evidente y dramática, mostrándose con total magnitud, sin la atenuación causada por los diferentes periodos de las muestras empleadas en las determinaciones generales de 2011 (Acuña, 2011) y 2022 (Acuña, 2022; esta publicación). La Figura 4 indica como el aumento del MSL en la costa central del Venezuela ha experimentado durante los últimos 30 años de monitoreo satelital, un incremento sostenido, variando desde +1.8 mm/año en la última década del siglo pasado, a los +2.9 mm/año y +5.7 mm/año, en las 2 primeras décadas de este siglo, respectivamente. Este incremento en la velocidad del aumento del MSL (o aceleración) se estima que sea en la actualidad cercano a +0.2 mm/año.

Entonces, para hablar de cifras definitivas podríamos establecer las siguientes cantidades: en La Guaira, el aumento del nivel medio del mar, como determinado por altimetría satelital multimisión durante el periodo 1992-2022, presenta magnitud de +3.4 ± 0.4 mm/año para la época media 2007.0, y está acelerándose a razón de +0.19 ± 0.07 mm/año.

Estos resultados se corresponden con las estimaciones publicadas por el IPCC en su último informe (IPCC/ONU, 2021); siendo además ratificados por determinaciones independientes de las misiones altimétricas de mayor precisión TOPEX/Poseidon, Jason-1, Jason-2 y Jason-3; ver Figura 5.

Figura 5. Estimaciones de la tasa de aumento del nivel medio del mar en La Guaira obtenidas independientemente a partir de datos registrados por los satélites TOPEX/Poseidon (TP), Jason-1 (J1), Jason-2 (J2) y Jason-3 (J3), en el periodo 1992-2022.

De nuevo, las estimaciones ahora independientes de las 4 misiones altimétricas indicadas en la Figura 5, vuelven a confirmar la significativa aceleración de la tasa de aumento del nivel medio del mar en el Caribe venezolano en los últimos 30 años; las observaciones satelitales indican que desde comienzos de la década de los 90's a la actualidad, la tasa se ha incrementado hasta por un factor 3, desde aprox. los +2 mm/año en 1992 a los +6 mm/año en 2022.

Tales números tienen efectos muy serios sobre las predicciones para 2050 y 2100 en cuanto a la posible pérdida (o retroceso) de la línea de costa venezolana en sus fachadas lacustre, caribeña y atlántica; y en lo que se refiere al aumento de la superficie de tierras costeras inundables por la penetración del mar en los territorios continentales e insulares del país.

En las zonas costeras de tierras de bajas y de pendientes suaves negativas hacia el mar, y por tanto, fácilmente erosionables e inundables, como las ubicadas en muchas zonas del litoral central, occidente y oriente de Venezuela, una tasa de aumento del nivel del mar como la arriba reportada (+3.4 mm/año) representará un posible ascenso de ese nivel para el 2050 de +15 cm, y para el 2100 de +32 cm.

Teniendo en cuenta que por cada 1-centímetro de ascenso en el MSL, la línea de costa podría retroceder aprox. 1-metro (Acuña, 2011), las anteriores estimaciones reprensentarían retrocesos de 15 m y 32 m para las fechas indicadas.

Ahora bién, tales valores ya habian sido proyectados en (Acuña, 2011) bajo la premisa que esa tasa de ascenso del MSL se mantuviese constante en el tiempo. Sin embargo, como hemos señalado antes en esta nota, esa no es la situación actual. Considerando la aceleración del MSL calculada aquí para el país (+0.2 mm por año), los valores proyectados antes para 2050 y 2100 posiblemente se multiplicarían por factores 3 y 8, respectivamente, ó quizás en una mayor proporción. Esto es dramático.

En (Acuña, 2011) se proyectaba que Venezuela tendría para el 2100 aprox. 1685 km² de tierras costeras bajo el nivel del mar (un 0.19% del territorio nacional continental), hoy para esa misma fecha, y si los factores que agudizan el calentamiento global (p.ej., el aumento anual en la concentración atmosférica de gases con efecto invernadero) no son ralentizados por políticas globales de protección medio ambientales, Venezuela tendrá bajo aguas marinas y lacustres casi 14 mil km² de superficies costeras (un 1.59% del total de tierras continentales del país); esto sin contar la muy posible y significativa reducción de la superficie de los territorios insulares nacionales, que en casos como Isla de Aves, pueda producir su desaparición definitiva, y en consecuencia, la pérdida de inmensas porciones de mar territorial venezolano y zona económica exclusiva.

Para más detalles sobre esta nota, puede consultar la siguiente publicación:

Acuña, G. (2022): "Actualizando la determinación de la tasa de ascenso del nivel medio del mar en el Caribe venezolano". Tópicos de Geodesia Geométrica. Julio 3-9, 2022. Cátedra Geodesia Geométrica (GGenLUZ). Lab. de Geodesia Física y Satelital. Dpto. de Geodesia Superior. Facultad de Ingeniería. Universidad del Zulia. Maracaibo, Venezuela.



domingo, 26 de junio de 2022

Evaluando la efectividad de ETHpredictor_v4.1 en la generación de rentabilidades diarias

InvestigaciónExtensiónSoftware ETHpredictor_v4.1.yab/.exe

ETHpredictor_v4.1 [Acuña, 2022] es una herramienta computacional de predicción del precio actual del Ethereum - ETH (para valores máximo y mínimo diarios) desarrollada recientemente por GGenLUZ; ver (https://ggenluz.blogspot.com/ 2022/05/ethpredictorv41-software-para-la.html).

A continuación se muestra una evaluación de la efectividad de ETHpredictor_v4.1 cuando se utiliza para generar rentabilidades diarias en dólares con ETH a través de operaciones de trading manual en el mercado de futuros de criptoactivos, i.e., operaciones financieras de compra y venta del ETH usando apalancamiento, en plataformas como Binance, Bitget, FTX, BybitCoinbase, etc.

La predicción para hoy domingo 26.06.2022 de ETHPredictor_v4.1 es mostrada en la Figura 1. Ésta fue calculada mediante colocación por cuadrados mínimos (LSC) con los valores del precio de apertura, cierre, máximo y mínimo del ETH para el día anterior (ayer, sábado 25.06.2022), además de datos similares de hasta un mes atrás. Estos datos los publica regularmente el portal Investing.com (https://es.investing.com/crypto/ethereum/historical-data).

Figura 1. Predicción de valores diarios máximo y mínimo del Ethereum 
con ETHpredictor_v4.1 para el día domingo 26.06.2022.

Tal como se verifica en la figura, ETHpredictor_v4.1 predijo por LSC que el valor del ETH alcanzaría una cotización máxima de 1270 ± 11 USD, y una cotización mínima de 1196 ± 13 USD.

Ahora bien, la Figura 2 muestra los valores reales del precio del ETH con resolución temporal de 5 minutos para el día de hoy, hasta aprox. el mediodía. El gráfico ha sido generado en la plataforma TradingView (https://es.tradingview.com/) y se presenta en forma de velas japonesas, un estándar para la visualización de precios en operaciones de trading.

Figura 2. Precios reales del Ethereum (resolución temporal de 5 minutos)
 para el día domingo 26.06.2022.

En la Figura 2 puede verse que, hasta el mediodía de hoy domingo, el precio del ETH alcanzó primero un valor máximo de 1278.75 USD a las 08:10 HLV, para luego caer y registrar un valor mínimo de 1210.00 USD a las 11:05 HLV.

Comparando la información en ambas figuras se puede verificar el alto nivel de acierto de la predicción de ETHpredictor_v4.1. Los valores  máximo y mínimo del precio del ETH estimados por ETHpredictor_v4.1 sólo se desvían 0.7% (8.75 USD) y 1.2% (14 USD) de las cotizaciones reales alcanzadas hoy por ese criptoactivo, respectivamente.

Si los valores de la predicción con ETHpredictor_v4.1 se utilizaran en el trading de ETH, el resultado sería el siguiente: en una plataforma de trading profesional operando en el mercado de futuros como Binance (https://www.binance.com/es), y usando un apalancamiento moderado de 10x, al momento en que el ETH alcanzara el valor máximo de predicción, se abriría una operación short (apuesta a la baja) con un monto de inversión de, p.ej., 100 USD. Esta operación se mantendría activa hasta que el ETH presentara un valor próximo al mínimo de predicción. Al cerrar en este punto la operación de trading, se habría obtenido una ganancia/día de 4.96%, esto con el apalancamiento a 10x significa una rentabilidad efectiva de 49.6%, es decir, 49.6 USD adicionales al monto de 100 USD utilizado como inversión !!!, y en un sólo día de trading ejecutando una única operación !!!

Para más información sobre ETHpredictor_v4.1.yab/.exe, y su forma de utilización en el trading manual o automático del Ethereum, Bitcoin y/o cualquier otra criptomoneda, contacte a su autor, Prof. Gustavo Acuña, a través del email: gacuna@fing.luz.edu.ve, o vía WhatsApp al +58-412-4271579.

Más detalles sobre el rendimiento de ETHpredictor_v4.1.yab/.exe pueden consultarse en la siguiente publicación:

Acuña, G. (2022): Evaluando la efectividad de ETHpredictor_v4.1 en la generación de rentabilidades diariasTópicos de Geodesia Geométrica. Julio 10-16, 2022. Cátedra Geodesia Geométrica (GGenLUZ). Laboratorio de Geodesia Física y Satelital. Dpto. de Geodesia Superior. Facultad de Ingeniería. Universidad del Zulia. Maracaibo, Venezuela.


miércoles, 22 de junio de 2022

SR950015P_v1.0: Base de datos para posiciones de vértices geodésicos en Venezuela, datum SIRGAS-REGVEN, realizaciones 1995, 2000, 2015 + 2022

 InvestigaciónExtensión, Software - SR950015P_v1.0.yab/.exe

GGenLUZ pone a disposición de la comunidad geodésica nacional el software SR950015P_v1.0.yab/.exe [Acuña, 2022a], el cual permite acceder a las posiciones finales de los vértices geodésicos de mayor importancia en el país, ver Figura 1.

Figura 1. Estaciones SIRGAS-REGVEN (185) en SR950015P_v1.0.yab/.exe.

SR950015Pv1.0.yab/.exe [Acuña, 2022a] materializa una simple base de datos que almacena y suministra a los usuarios nacionales de GGenLUZ, las posiciones definitivas de los vértices geodésicos SIRGAS, REGVEN y REMOS en el país, correspondientes a las realizaciones SIRGAS-REGVEN de 1995 (ITRF94_1995.4/GRS80), 2000 (ITRF2000_2000.4/ GRS80) y 2015 (ITRF2014_2015.4/GRS80).

SR950015Pv1.0.yab/.exe contiene hasta ahora 185 vértices geodésicos medidos en Venezuela entre Mayo,1995 y Junio,2022 empleando técnicas GNSS de precisión. Esta versión sólo considera vértices de los órdenes A (±1cm) y B (±2cm). En próximas versiones, SR950015P  incluirá monografías digitales y vértices SIRGAS-REGVEN de orden C (±5cm), correspondientes a las diferentes densificaciones del marco geodésico nacional realizadas hasta ahora (p.ej., por el IGVSB, PDVSA, GGenLUZ/LGFS-LUZ, alcaldías, gobernaciones, universidades, otros organismos del estado venezolano, empresa privada, etc.).

Además de brindar las posiciones finales de los vértices en las realizaciones SR1995, SR2000 y SR2015, en términos de coordenadas geodésicas curvilíneas LAT,LON,ALT.ELIP; cartesianas geocéntricas XYZ; y planas proyectadas globales N,E_UTM; con sus respectivas estimaciones de error; SR950015Pv1.0.yab/.exe genera también información geodésica relevante para cada locación que complementa las posiciones, por ejemplo, velocidades geodésicas en LAT,LON,ALT.ELIP, elevación por nivelación, ondulación del geoide, elevación GNSS/geoide, gravedad total, anomalía de gravedad de aire-libre, componentes de la deflexión de la vertical, potencial de gravedad y número geopotencial. La mayoría de estas cantidades provienen de observaciones y/o de los modelos nacionales de velocidad VVM20v1.0 [Acuña, 2022b], y de geoide/gravedad VGM19v3.0 [Acuña, 2021].

Finalmente, SR950015Pv1.0.yab/.exe muestra las posiciones de los vértices geodésicos en la solución GGenLUZ SIRGAS-REGVEN(2022) [Acuña, 2022c], expresada en el nuevo ITRF2020 (https://itrf.ign.fr/en/solutions/ITRF2020), y correspondiente a la época actual 2022.5, elipsoide GRS80. Esta solución proviene del reprocesamiento GPS/GLN hecho por GGenLUZ de REGVEN2015 empleando GipsyX (https://gipsy-oasis.jpl.nasa.gov/), y luego propagando las coordenadas obtenidas en ITRF2014_2015.5 al ITRF2020_2022.5. Para esto se empleó el grupo oficial de 14 parámetros de transformación del IERS/ITRF2020 entre ambos marcos (https://itrf.ign.fr/docs/solutions/itrf2020/Transfo-ITRF2020_TRFs.txt), el modelo de velocidad venezolano VVM20v1.0 [Acuña, 2022b] y el software Trn_ITRF, v4.0 [Acuña, 2022d].

SR950015Pv1.0.yab/.exe, y la mayor parte de la información geodésica contenida en la aplicación, fue preparada por el Prof. Gustavo Acuña, cátedra Geodesia Geométrica de LUZ (GGenLUZ; https://ggenluz.blogspot.com), durante el transcurso de 2022. Un ejemplo de sus resultados se presenta a continuación. Éstos refieren a la estación SIRGAS-REGVEN Maracaibo (MARA) como obtenidos de la versión DEMO del software.



Figura 2. Resultados típicos de SR950015Pv1.0.yab/.exe, archivo de texto MARA.out.

La versión DEMO de SR950015Pv1.0.yab/.exe está disponible para los usuarios nacionales de GGenLUZ. Para descargarla contacte a su autor, Prof. Gustavo Acuña, a través del email gacuna@fing.luz.edu.ve, o acceda directamente al siguiente enlace:


Para más detalles sobre SR950015Pv1.0.yab/.exe, puede consultar la siguiente publicación:

Acuña, G. (2022): SR950015P_v1.0: Base de datos para posiciones de vértices geodésicos en Venezuela, datum SIRGAS-REGVEN, realizaciones 1995, 2000, 2015 + 2022Tópicos de Geodesia Geométrica. Junio 19-25, 2022. Cátedra Geodesia Geométrica (GGenLUZ). Laboratorio de Geodesia Física y Satelital. Dpto. de Geodesia Superior. Facultad de Ingeniería. Universidad del Zulia. Maracaibo, Venezuela.


sábado, 18 de junio de 2022

GGenLUZ produce VVM20, el más reciente modelo de velocidad para estaciones de control geodésico en Venezuela

 InvestigaciónExtensiónDatos/Soluciones/Modelos VVM20, Venezuelan geodetic Velocity Model 2020

Continuando la serie de modelos nacionales VVM (VVM12, VVM15, VVM17; Acuña, 2015-2020), en esta oportunidad GGenLUZ produce VVM20 - Venezuelan geodetic Velocity Model 2020, el más reciente modelo de velocidad para estaciones de control geodésico en Venezuela (Acuña, 2022), ver Figura 1.

Figura 1. VVM20 (Acuña, 2022) - modelo de velocidad para estaciones
 de control geodésico en Venezuela. 

VVM20 aproxima las variaciones respecto al tiempo en la posición horizontal de estaciones de control geodésico en Venezuela, por efecto del movimiento permanente de las placas tectónicas que se identifican en el territorio nacional (i.e., CA: Caribe, SA: Suramérica, ND: NorAndes).

VVM20 utiliza los modelos geológico-geofísicos NNR NUVEL-1A (DeMets, Gondon, Argus y Stein, 1994), PB2002 (Bird, 2003) y MORVEL56 (DeMets, Gordon y Argus, 2010), y los geodésicos APKIM2005 (Drewes, 2009), VEMOS2009 (Drewes y Heidbach, 2009), ITRF2008 (Altamimi, Métivier y Collilieux, 2012), VEMOS2015 (Sánchez y Drewes, 2016) y VEMOS2017 (Drewes y Sánchez, 2020) para estimar una solución preliminar, media ponderada, que luego es ajustada mediante colocación por cuadrados mínimos (LSC) a las velocidades "reales(medidas)", como determinadas con GNSS, en 98 estaciones de control geodésico del país. Tales velocidades reales corresponden a estimaciones de la red continental SIRGAS-CON (www.sirgas.org) y aquellas derivadas de la comparación de las posiciones de la red geocéntrica nacional REGVEN entre sus realizaciones de 1995, 2000 y 2015 (www.igvsb.gob.ve), y de la red GNSS nacional activa REMOS, previamente transformadas a la última realización del ITRS, el ITRF2020 (https://itrf.ign.fr/en/solutions/ITRF2020; Altamimi, 2022).

Al igual que los modelos antes señalados, VVM20 sólo aproxima velocidades horizontales, posibilitando calcular las componentes del vector velocidad en latitud, longitud y XYZ como referidas al ITRF2020, para cualquier estación de control geodésico ubicada en territorio venezolano y zonas adyacentes.

VVM20 fue calculado por el Prof. Gustavo Acuña, cátedra Geodesia Geométrica de LUZ (GGenLUZ; https://ggenluz.blogspot.com), durante los meses transcuridos de 2022. El modelo es una versión mejorada de VVM17 (Acuña, 2020) producto del reprocesamiento de REGVEN2015 con GipsyX (https://gipsy-oasis.jpl.nasa.gov/), y de la inclusión de nuevas y más precisas velocidades reales estimadas en REMOS.

La Tabla 1 muestra un ejemplo de la calidad que puede obtenerse cuando se estiman velocidades con VVM20 en estaciones de control geodésico en Venezuela.
 
----------------------------------------------------------------------
       VELOCIDAD SIRGAS-REGVEN            VELOCIDAD VVM20 
ESTACIÓN   LAT(mm/a) LON(mm/a) VT(mm/a)   LAT(mm/a) LON(mm/a) VT(mm/a)
----------------------------------------------------------------------
  MARA       13.50      8.90    16.17       13.48      8.92    16.16 
 AMUA       11.60     11.90    16.62       13.82     10.53    17.37
 JUNQ       11.30      1.40    11.39       13.27      4.00    13.86
 AGUA       11.00     -3.90    11.67       11.48     -4.17    12.21
 CANO       12.10     -3.50    12.60       13.42     -2.87    13.72
----------------------------------------------------------------------
Tabla 1. Velocidades en estaciones SIRGAS-REGVEN de Venezuela
 obtenidas por mediciones GNSS y empleando VVM20.

La tabla presenta las velocidades "reales(medidas)" obtenidas por GNSS a través de la comparación de los resultados de las campañas REGVEN 1995, 2000 y 2015 en 5 estaciones SIRGAS-REGVEN (orden A) de Venezuela, a saber, MARACAIBO, AMUAY, JUNQUITO, AGUA LINDA y LA CANOA.

Las velocidades se muestran en mm/año para sus componentes en latitud y longitud, y velocidad total 2D. Además, la tabla presenta valores respectivos obtenidos del modelo VVM20.

En términos generales, las velocidades reales y aquellas obtenidas con VVM20 coinciden bastante bién en magnitud y dirección, presentando diferencias con RMS de sólo ±1.28 mm/a. Para esta selección de estaciones, las velocidades VVM20 se muestran levemente "más rápidas" que aquellas provenientes de las mediciones GNSS, un efecto del proceso de modelaje con LSC. A un nivel de precisión de ±1 mm/año (1s), VVM20 reproduce eficazmente las variaciones temporales de la posición en vértices geodésicos de Venezuela.

Tales variaciones también pueden  interpretarse como "deformaciones" que evolucionan con el tiempo en las principales redes geodésicas de control del país REGVEN y REMOS, en principio debidas a la téctónica regional de placas.

Una versión DEMO de VVM20 está disponible para los usuarios nacionales de GGenLUZ. Para descargarla contacte a su autor, Prof. Gustavo Acuña, a través del email gacuna@fing.luz.edu.ve.

Más detalles sobre VVM20 pueden consultarse en la siguiente publicación:

Acuña, G. (2022): VVM20, el más reciente modelo de velocidad para estaciones de control geodésico en VenezuelaTópicos de Geodesia Geométrica. Junio 12-18, 2022. Cátedra Geodesia Geométrica (GGenLUZ). Laboratorio de Geodesia Física y Satelital. Dpto. de Geodesia Superior. Facultad de Ingeniería. Universidad del Zulia. Maracaibo, Venezuela.