Introducción
Actualmente, con el posicionamiento GNSS (Global Navigation Satellite System) es posible determinar la ubicación de un objeto sobre la superficie terrestre, mediante un sistema de coordenadas tridimensionales que describe con exactitud la posición de interés. GPS (Global Positioning System) es el sistema más usado a nivel mundial, cuya superficie de referencia es un elipsoide de revolución y permite a través de coordenadas geodésicas: latitud, longitud y altura elipsoidal, determinar la posición de un objeto en el espacio (Seeber, 2003). En la componente horizontal (latitud y longitud), las coordenadas pueden ser usadas de forma directa; sin embargo, en la componente vertical debido a que, por sus características geométricas, no es factible su aplicación para trabajos de ingeniería que requieren contar con alturas referidas a la superficie topográfica real.
Alturas de carácter físico están referidas al geoide, el cual, en países como Ecuador, que cuentan con escasos datos de gravedad para modelar el comportamiento dinámico de esta superficie, resulta difícil de obtener. Para suplir la necesidad de contar con alturas sobre la superficie topográfica, el Instituto Geográfico Militar (IGM) del Ecuador llevó a cabo la Red de Control Básico Vertical (RCBV) referida al datum vertical oficial del país, ubicado en el mareógrafo de la Libertad - Santa Elena (Paredes, 1986).
En términos prácticos, para trabajos de ingeniería, se considera el nivel medio del mar como una superficie equivalente al geoide. Por ende, a la altura nivelada y la altura ortométrica respectivamente se las puede considerar de la misma forma (Odera & Fukuda, 2015). La determinación de alturas sobre el nivel medio del mar (niveladas) es llevada a cabo a través de campañas de nivelación geométrica que recorren las vías estatales del país y que constituyen la RCBV (Figura 1). Esta técnica de medición de alturas es también usada para la generación de modelos de ondulación geoidal, que combinados con datos GPS y nivelación, permiten aproximar un geoide local (Palacios et al., 2021). A pesar de ser una técnica muy precisa, la nivelación geométrica demanda de tiempo y recursos técnico/ económicos dispendiosos, por lo que es necesario contar con métodos que permitan reducir los costos de obtención de alturas niveladas con precisiones suficientes para trabajos de ingeniería y las geociencias en general (Palacios, 2019).
El cálculo de alturas niveladas, se puede realizar a partir de la resta de la altura elipsoidal con la ondulación geoidal en el punto de interés (Kenyeres, 2016) (Ec. 1). Sin embargo, esta expresión representa un comportamiento ideal de las variables inmersas, lo cual no es correcto, ya que sus valores son obtenidos mediante instrumentos ópticos (nivel) y electrónicos (GNSS) para la altura nivelada y elipsoidal respectivamente. De ahí es lógico inferir la presencia de un error aleatorio subyacente en la medición, con lo que el modelo determinista cambia a uno probabilístico como el de la Ec. 2.
Hn ≅ h - n (1) n = h - Hn + ɛ (2)
Donde n es la ondulación geoidal, h altura elipsoidal, Hn altura nivelada, ɛ representa una cantidad pequeña de variabilidad aleatoria (Torge, 2001). Este error estocástico puede ser eliminado al considerar la separación de la ondulación (∆N), en lugar de la separación absoluta del geoide (n), dado por (Banerjee et al., 1999):
∆N = (h2 - h1) - (H2 - H1) (3)
Donde h2, h2 son alturas elipsoidales del punto a determinar y conocido, H2, H1 representan las alturas ortométricas del punto a determinar y conocido. Pero, como se mencionó anteriormente, en países como Ecuador donde no existen suficientes datos de gravedad para adoptar un sistema de alturas físicas, se considera para fines ingenieriles la equivalencia entre las alturas ortométricas y niveladas, por tanto H1 0Hn 1 y H2 ≅ Hn 2 . De esta manera, Hn 2 se puede obtener de la Ec. 3, con la siguiente expresión:
Hn 2 = Hn 1 + (h2- h1) - (n2- n1) (4)
Donde n2, n1 corresponde a la ondulación geoidal para los puntos desconocido y conocido, y Hn 1, Hn 2 representa la altura nivelada en los puntos desconocido y conocido respectivamente. Los valores de ondulación geoidal son derivados del modelo gravimétrico EGM08, de cuya diferencia se eliminan los errores sistemáticos de este modelo y por tanto del método de nivelación GPS (Fotopoulos, Kotsakis, & Sideris, 2003).
El presente estudio tiene por objetivo general, demostrar la aplicabilidad de la nivelación GPS para la obtención de alturas niveladas precisas con fines de ingeniería. De este objetivo macro, se dividen dos específicos: a) determinar rangos máximos en que se puede usar nivelación GPS dentro del área urbana de la ciudad de Macas, y b) validar y comparar las alturas generadas mediante nivelación GPS con alturas niveladas calculadas con nivelación geométrica, a partir de los puntos de la red geodésica local del cantón Morona, en la región amazónica del Ecuador.
Metodología
Zona de Estudio
Morona es uno de los doce cantones de la provincia de Morona Santiago, perteneciente a la región amazónica del Ecuador. La capital cantonal de Morona es la ciudad de Macas, ubicada geográficamente a 2°17’54.94” de latitud Sur y 78°09’59.03” de longitud Oeste, asentada sobre una planicie al occidente del río Upano, a ≈40 Km de distancia del volcán Sangay (Palacios, 2018) (Palacios & Arellano, 2021) y rodeada por abundante vegetación propia de la región oriental (Palacios, Castro, & Rodríguez, 2019) (Palacios & Rodríguez, 2021). Morona limita al norte con los cantones Pablo Sexto y Huamboya, hacia el sur con Sucúa, mientras que al este y oeste colinda con Taisha, Guamote y Alausí, respectivamente (ver Figura 2) (Palacios & Toulkeridis, 2020). Demográficamente, el cantón posee una población proyectada al año 2021 de 61019 habitantes, y su capital Macas de 28035 habitantes, siendo la más poblada de la provincia, y un importante nodo de conexión entre las regiones Sierra y Amazonía (Palacios, 2020).
Vértices Geodésicos
En el año 2015 se realizó la red geodésica del cantón, la misma que está conformada por 20 vértices materializados con placas de aluminio colocadas sobre aceras o hitos Tipo B (IGM, 2006).
Cada punto de la red, cuenta con un control geodésico suplementario, conformado por: nombre del punto, latitud, longitud, altura elipsoidal, coordenadas planas proyectadas (este y norte), y altura nivelada. En la Figura 3, se muestra la distribución de los vértices pertinentes a la red geodésica local.
La red geodésica cantonal está enlazada a la REGME (Red GNSS de Monitoreo Continuo del Ecuador) con época de
referencia 2016.44, datum horizontal ITRF 2008, datum vertical al nivel medio del mar del mareógrafo de La Libertad. Las coordenadas de la componente horizontal fueron obtenidas mediante posicionamiento estático diferencial con equipos GNSS de dos frecuencias (L1 y L2), marca Trimble modelo R8S en distintas sesiones de tres horas de rastreo por cada punto. Por otro lado, la altura nivelada fue medida a través de nivelación geométrica con líneas de nivelación tipo radial partiendo de vértices de la RCBV. El ajuste de la red cumple con los parámetros del IGM para ser considerada de primer orden, tanto en horizontal (±5 cm) como vertical (±4 mm√K) (IGM, 2006).
Gradiente de ondulación geoidal
Un análisis de la variabilidad de la ondulación geoidal fue necesario para determinar las zonas en que se puede aplicar la nivelación GPS de manera más confiable. Para ello, se calculó el gradiente de la ondulación geoidal con los valores del modelo EGM08 con una resolución de 1 minuto de arco. De esta forma se pueden estimar los errores que existen en función de la distancia a la que se realiza la nivelación. Según Oropeza (2010), se pueden considerar cinco zonas para la variación respecto a la inclinación del terreno (gradiente), que se resumen en la Tabla 1.
Zona | Descripción | Rango (mm/Km) |
---|---|---|
1 | Baja | 0.0001 - 0.0379 |
2 | Ligeramente suave | 0.0380 - 0.0739 |
3 | Suave | 0.0740 - 0.1126 |
4 | Moderada | 0.1127 - 0.1608 |
5 | Fuerte | 0.1609 - 0.2412 |
Nota: Fuente propia de la investigación.
En la Figura 4, se muestra el gradiente de ondulación geoidal en ppm (mm de error en altura / distancia en Km desde el punto con altura nivelada conocida), así como las zonas de variación que están presentes en el área de estudio.
Líneas de nivelación consideradas
Como se mencionó en la introducción del estudio, uno de los objetivos es determinar rangos o límites dentro del área urbana de la ciudad de Macas (escala más pequeña de superficie de análisis). Para ello, se consideraron tres escenarios posibles a diferentes distancias y orientaciones. Estos se resumen a continuación en la Tabla 2:
Escenario | Punto conocido | Punto desconocido | Distancia (Km) | Orientación |
1 | Macas - MR | Macas 2 | 1.07 | Sur-Norte |
Macas 4 | Macas 2 | 2.23 | ||
2 | Macas | Macas 3 | 1.04 | Norte-Sur |
Macas 1 | Macas 3 | 2.27 | ||
3 | Macas 1 | Macas 2 | 0.19 | Norte-Sur |
Nota: Fuente propia de la investigación.
Por otro lado, para validar la nivelación GPS a nivel cantonal (escala más grande de superficie analizada), se determinaron dos líneas de nivelación en las dos zonas de variación del gradiente de ondulación geoidal en que se contaron con puntos de control, las cuales se resumen en la Tabla 3.
Escenario | Punto conocido | Punto desconocido | Distancia (Km) | Orientación |
---|---|---|---|---|
1 | Río Blanco | Sinaí | 30.2 | Sur - Norte |
Río Blanco | Santa Rosa MR | 19.0 | ||
Río Blanco | Macas 3 | 5.1 | ||
2 | 9 de Octubre | 9 de Octubre | 2.7 | Oeste - Este |
MR | ||||
9 de Octubre | San Isidro MR | 9.5 | ||
9 de Octubre | Gral. Proaño MR | 13.0 |
Nota: Fuente propia de la investigación.
Análisis y resultados
Según los valores del gradiente de ondulación geoidal, en el cantón Morona existe una variabilidad considerable, con errores de hasta 24 mm/Km, que implica la necesidad de definir zonas homogéneas y analizar la influencia de la distancia para aplicar correctamente la técnica de nivelación GPS. Los valores más elevados de gradiente se encuentran en las estribaciones de las cadenas montañosas presentes en el cantón, en la cordillera de los Andes y del Kutukú, al occidente y oriente respectivamente; mientras que los relieves menos escarpados como la llanura amazónica y parte del valle interandino que se encuentra en Morona, reflejan los más bajos valores de gradiente. Esto denota la influencia de la topografía en el comportamiento de la ondulación geoidal como variable física (Sánchez, 2003), y por tanto, en la obtención de alturas mediante nivelación GPS.
Un total de cinco zonas fueron definidas en función de los rangos de variabilidad del gradiente de ondulación geoidal, los
cuales propenden a mejorar la aplicabilidad del método planteado. Las zonas 3, 4 y 5 son las predominantes en el territorio cantonal, y en menor medida las zonas 1 y 2, las cuales definen un rango teórico de error en nivelación GPS entre 0.074 0.241 ppm. En el caso del área urbana de la ciudad de Macas (escala micro), es abarcada en su totalidad por la zona 4, lo que resulta en una ventaja para aplicar la medición de alturas sobre la topografía mediante esta técnica de nivelación.
En cuanto a los escenarios evaluados a nivel cantonal (escala macro), los vértices escogidos se encuentran en dos zonas de variabilidad. El punto Río Blanco fue escogido como base para realizar la nivelación GPS en la zona 4; mientras que el punto 9 de Octubre se consideró como base para la nivelación GPS en la zona de variación 5. Los resultados de ambas líneas de nivelación se presentan en la Gráfica 1.
Los dos escenarios testados a nivel cantonal cumplen con el criterio de unidireccionalidad de la línea de nivelación; es decir, que la nivelación GPS se realizó en una sola orientación (Sur-Norte para el escenario 1 y Oeste-Este para el escenario 2), de esta forma se evita una posible divergencia de valores en altura debido al comportamiento anisotrópico del geopotencial (Heiskanen & Moritz, 1985). Como se observa en la Gráfica 1, los errores que presenta la nivelación GPS en la zona 4 son menores que los hallados en la zona 5, lo cual concuerda con la variabilidad de la ondulación geoidal en cada una de ellas. Por un lado, en zona 4 se alcanzó un máximo de 12.19 cm a 30 Km (≈ 4.04 ppm) y un mínimo de 1.27 cm a 5.1 Km (≈ 2.49 ppm) de distancia del punto base; mientras que en la zona 5, el error máximo fue de 38.48 cm a 13 Km (≈ 29.60 ppm) y un mínimo de 8.33 cm a 2.7 Km (≈ 30.85 ppm) de distancia del vértice conocido.
En los escenarios analizados a nivel urbano de la ciudad de Macas, todos cumplen con el criterio de unidireccionalidad de la línea de nivelación, al igual que en el caso cantonal. En la Tabla 4, se resumen los resultados de la nivelación GPS a diferentes distancias para los escenarios 1, 2 y 3.
Escenario | Punto conocido | Punto desconocido | Distancia (Km) | Error (cm) |
1 | Macas MR | Macas 2 | 1.07 | 2.48 |
Macas 4 | Macas 2 | 2.23 | 5.17 | |
2 | Macas | Macas 3 | 1.04 | 1.51 |
Macas 1 | Macas 3 | 2.27 | 4.69 | |
3 | Macas 1 | Macas 2 | 0.19 | 0.14 |
Nota: Fuente propia de la investigación.
Si se comparan los valores del error de nivelación GPS entre el escenario 1 y 2 de la Tabla 4 que presentan similitud en su planteamiento (distancia del punto base), los resultados fueron parecidos. Por tanto, es factible promediar y determinar un error medio de nivelación a las distancias evaluadas de 1 y 2.2 Km aproximadamente, con una equivalencia de 1.99 cm y 4.93 cm respectivamente. La distancia mínima de separación entre el punto conocido y desconocido (cerca de 200 metros) corresponde al escenario 3 en que se obtuvo un error de 1.4 mm.
De los resultados hallados tanto a nivel cantonal (Gráfica 1) como urbano (Tabla 4), se observa que el error en la nivelación GPS depende principalmente de dos factores condicionantes que son la distancia y la variación del gradiente de ondulación geoidal. Según Blanco (2010), esta técnica de nivelación se puede aplicar hasta una distancia de 15 Km; sin embargo, como se evidenció en este estudio, es necesario el análisis del comportamiento de gradiente para garantizar una mejor obtención de alturas. Además, trabajos anteriores sugieren que el error en nivelación GPS aumenta trigonométricamente (Sánchez & Martínez, 1997), aunque esto último estará en función nuevamente de la variación del gradiente en el área que se aplique.
Con el fin de evidenciar la divergencia del valor en altura que existe al partir desde distintos puntos, se plantearon dos ejemplos de nivelación GPS con trayectorias contrarias desde los vértices usados para arrastrar el valor en la componente vertical. El primer caso se consideró en el área urbana de Macas, con la nivelación desde los vértices Macas 3 Macas - MR (Sur - Norte) versus Macas 2 Macas - MR (Norte - Sur), en la cual se alcanzó una diferencia del valor en altura de 4.1 mm en una distancia de 1 Km aproximadamente. Un segundo caso fue la nivelación desde los vértices Sevilla
Santa Rosa (Sur - Norte) versus Sinaí Santa Rosa (Norte - Sur), donde se presentó una diferencia de 6.79 cm en el valor de altura con cerca de 13 Km de distancia. Este cambio en el error de la nivelación que existe al partir desde distintos puntos con trayectorias opuestas de nivelación, se puede atribuir al comportamiento anisotrópico del geopotencial, ya que, al tratarse de una técnica de obtención de alturas geométricas, la nivelación GPS depende de la posición de partida para arrastrar la altura cuyo valor de geopotencial es específico para cada punto base, a diferencia de alturas de tipo físicas que son independientes de la trayectoria seguida y con un único valor de potencial (Hofmann & Moritz, 2005).
Por otro lado, las zonas de variación de gradiente de ondulación geoidal que se determinaron en el territorio cantonal, tienen el objetivo de evitar errores al momento de aplicar la nivelación GPS, ya que resultaría erróneo calcular el valor de altura de un punto ubicado en una zona distinta a la del vértice base debido a la influencia del cambio de gradiente. Esto último se comprobó con tres ejemplos que se resumen en la Tabla 5.
Como se puede observar, los errores hallados en los ejemplos de la Tabla 5 son mayores a los resultados de las líneas de nivelación que se aplicaron en una misma zona de variación (Gráfica 1yTabla 4); además, se observa que los errores de la nivelación sugieren una cierta independencia en la distancia desde el punto base al punto desconocido cuando se trabaja en zonas heterogéneas de variación, sino que está en función del gradiente de gravedad de la ondulación geoidal en cada zona.
Con estos antecedentes, es posible establecer rangos o zonas de error conocido, en las que se puede aplicar la nivelación GPS para obtener alturas niveladas. En el caso del área urbana de la ciudad de Macas, al estar ubicada completamente en la zona 4 de variación de gradiente, resulta de gran ayuda mapear estos rangos de distancia para facilitar la medición de alturas niveladas con precisiones auspiciosas en trabajos de ingeniería, obtenidas de forma rápida con técnicas de posicionamiento GNSS, como se muestra en la Figura 5.
Punto conocido | Zona del punto conocido | Punto desconocido | Zona del punto desconocido | Distancia (Km) | Error (cm) |
---|---|---|---|---|---|
Gral. Proaño | Zona 5 | Sinaí MR | Zona 4 | 20.9 | 37.60 |
9 de Octubre | Zona 5 | Santa Rosa MR | Zona 4 | 17.3 | 62.54 |
Macas | Zona 4 | 9 de Octubre MR | Zona 5 | 15.4 | 58.36 |
San Isidro MR | Zona 5 | Macas 1 | Zona 4 | 10.3 | 46.08 |
Macas 3 | Zona 4 | Gral. Proaño | Zona 5 | 5.3 | 26.75 |
Nota: Fuente propia de la investigación.
A manera de resumen, en la Tabla 6 se muestran las diferencias en mm/Km entre el valor de altura nivelada calculado mediante nivelación geométrica versus el valor generado a través de nivelación GPS de los puntos desconocidos planteados en este estudio, tanto a escala urbana y cantonal.
Conclusiones
Debido a la topografía irregular del cantón Morona, el gradiente de ondulación geoidal presenta variaciones desde ligeras hasta fuertes, que resultan en cinco zonas de variación homogéneas que fueron determinadas, donde las zonas 3 (suave), 4 (moderada) y 5 (fuerte) son las predominantes en el territorio, ubicadas en la llanura amazónica, estribaciones montañosas y las cadenas montañosas, respectivamente.
La obtención de alturas niveladas fue evaluada en dos escalas de trabajo, una cantonal que abarca una superficie mayor de territorio y otra urbana correspondiente a una extensión menor de superficie. En el primer caso, los errores de nivelación GPS estuvieron en orden de 1.27 cm - 12.19 cm a distancias entre 5.1 Km - 30.2 Km para la zona 4; y 8.33 cm - 38.48 cm a 2.7 Km - 13 Km de líneas base para la zona 5. Por otro lado, en el área urbana de la ciudad de Macas, se alcanzó errores de 0.14 cm, 1.99 cm y 4.93 cm en distancia de 0.2 Km, 1 Km y 2.2 Km respectivamente.
Mediante la correcta aplicación de la técnica de nivelación GPS, es factible obtener alturas niveladas oficiales en el Ecuador (referidas al datum vertical del país), con precisiones auspiciosas para trabajos de ingeniería y geociencias en general, que requieren alturas sobre la superficie topográfica.
La nivelación GPS está en función de dos variables principales: la distancia de línea base y la variación del gradiente, por lo que es necesario considerarlas para aplicar esta técnica y obtener alturas niveladas más precisas.
Declaración de la contribución de los autores
Todos los autores afirmamos que se leyó y aprobó la versión final de este artículo.
El porcentaje total de contribución para la conceptualización, preparación y corrección de este artículo fue el siguiente:
I.F.P.O. 60 % y K.D.A.P. 40 %.
Escala | Punto desconocido | Altura nivelación geométrica | Altura nivelación GPS | Error (mm/Km) |
Urbano | Macas 2 | 1053.8416 | 1053.8168 | 23.18 |
Macas 3 | 1016.8931 | 1016.9082 | 14.52 | |
Cantonal | Sinaí | 1192.6969 | 1192.8188 | 4.04 |
Santa Rosa MR | 1141.7527 | 1141.7088 | 2.31 | |
Macas 3 | 1016.8931 | 1016.9058 | 2.49 | |
9 de Octubre MR | 1856.2138 | 1856.2971 | 30.85 | |
San Isidro MR | 1150.8208 | 1150.9951 | 18.35 | |
Gral. Proaño MR | 1110.8803 | 1111.2651 | 29.60 |
Nota: Fuente propia de la investigación.