CONEXÃO DE SISTEMAS VERTICAIS DE REFERÊNCIA LOCAIS AO SISTEMA GEODÉSICO BRASILEIRO COM BASE EM UM SISTEMA VERTICAL DE REFERÊNCIA GLOBAL

Connection of Local Vertical Reference Systems to Brazilian Geodetic System based on a Global Vertical Reference System

Resumos

Os Modelos Globais do Geopotencial (MGGs), em especial os oriundos da missão GOCE, associados a Modelos Digitais de Altitude (MDAs) com resoluções cada vez melhores, têm possibilitado novas alternativas para a obtenção de altitudes com significado físico. Neste trabalho, propõe-se a realização de um sistema local de altitudes, em uma área onde o Sistema Geodésico Brasileiro (SGB) é deficiente, e vinculá-lo a um Sistema Vertical de Referência Global (SVRG) e ao SGB. Os MGGs foram validados utilizando Referências de Nível associadas com posicionamento GPS. Para modelar o erro de omissão dos MGGs, foi utilizada a técnica Residual Terrain Model (RTM). Os números geopotenciais foram obtidos com a resolução simplificada do Problema do Valor de Contorno da Geodésia (PVCG) na forma fixada. A área de estudos é o Estado do Amapá, que utiliza o Datum Vertical Brasileiro de Santana (DVB-S), e para validação dos MGGs foram utilizadas estações do Estado do Pará, vinculadas ao Datum Vertical Brasileiro de Imbituba (DVB-I). Como resultados, obteve-se números geopotenciais e altitudes físicas consistentes. Na validação dos modelos, o RTM não teve influência significativa. Utilizando o MGG do GOCE, GO_CONS_SPW_R4, obteve-se que o segmento da rede vinculado ao DVB-S está 1,30m acima do segmento vinculado ao DVB-I.

Palavras-chave:
Modelos Globais do Geopotencial; Sistema Vertical de Referência Global; Problema do Valor de Contorno da Geodésia; RTM

Global Geopotential Models (GGMs), specially those provided by GOCE mission, associated with Digital Elevation Models (DEMs) with increasingly higher resolutions, has enable new alternatives for obtaining physical heights. The purpose of this work is the realization of a local height system, in a region with poor data from the Brazilian Geodetic System (BGS), and link it to a Global Vertical Reference System (GVRS). The GGMs were validated using GPS/leveling stations. The omission error of GGMs was modeled by using the Residual Terrain Model (RTM) technique. The geopotential numbers were obtained through a simplified solution of the Geodetic Boundary Value Problem (GBVP) in a fixed approach. The study area is the Amapá state, whose altitudes are linked with the Brazilian Vertical Datum of Santana (BVD-S), and for validation of GGMs some stations linked to Brazilian Vertical Datum of Imbituba (BVD-I) in the state of Pará were also used. As a result, geopotential numbers and heights with physical meaning were obtained. In the validation of the models, RTM technique had no significant influence. With the GGM from GOCE, GO_CONS_SPW_R4, it was obtained that the network linked to BVD-S is 1,30m above the network linked to BVD-I.

Keywords:
Global Geopotential Models; Global Vertical Reference System; Geodetic Boundary Value Problem; RTM

1. Introdução

As altitudes dotadas de significado físico, necessárias às aplicações em obras de engenharia, por exemplo, podem ser determinadas de maneira convencional por técnicas de levantamento como o nivelamento geométrico; estas altitudes são vinculadas a um Datum Vertical (DV), normalmente um marégrafo que fornece o Nível Médio do Mar (NMM) local referido a um determinado período de observações. Porém, os sistemas de altitudes dos países da América do Sul não levaram em consideração a influência da gravidade, mas introduziram, como no caso do Brasil, uma correção teórica da gravidade (correção do não paralelismo das superfícies equipotenciais utilizando o modelo da Terra Normal). Desta forma, o efeito sistemático do Datum e a não consideração da gravidade nos levantamentos gera inconsistências nos sistemas de altitudes, principalmente quando se trata de conexão entre países (De Freitas e Blitzkow, 1999; Sánchez, 2002).

Uma abordagem moderna do problema é possibilitada pela determinação da gravidade por meio dos Modelos Globais do Geopotencial (MGGs) e o conhecimento da topografia local tal como o expresso em Modelos Digitais de Altitude (MDAs), oriundos de plataformas orbitais. Neste contexto, os chamados Sistemas Verticais de Referência Locais (SVRLs) ganharam grande possibilidade de realização com base em redes com melhor adequação às redes nacionais e globais. No caso brasileiro, o vínculo de Sistemas e Redes Verticais Locais ao Sistema Geodésico Brasileiro (SGB) e a um Sistema Vertical de Referência Global (SVRG). Isto decorre justamente das novas possibilidades advindas a partir das ferramentas mais atuais da Geodésia, em particular dos MGGs derivados da missão GOCE (Gravity field and steady-state Ocean Circulation Explorer) em 2012 e 2013, em associação com posicionamento GNSS e o uso de MDAs (Gerlach e Rummel, 2012; Rummel, 2013).

Na presente pesquisa, a proposta é definir um Sistema Vertical de Referência Local, com sustentação em modelo global consistente conforme preconizado por (Rummel, 2013 e Gerlach e Rummel, 2012), em uma região com deficiência de informação altimétrica da rede oficial do país, ou seja, onde o Sistema Geodésico Brasileiro apresenta uma lacuna. Para isto, foi constituída uma base de dados que possibilitou a realização de testes de aderência na área de estudos, observando a consistência dos referenciais e sistemas de maré permanente. Os modelos globais do geopotencial foram avaliados em termos de sua compatibilidade com uma base geodésica conhecida de forma a validar suas aplicações, e o preenchimento dos vazios espectrais (erro de omissão) dos modelos do geopotencial foi efetuado com o emprego da técnica RTM - Residual Terrain Model. Os números geopotenciais foram obtidos a partir da solução do tipo Brovar para o Problema do Valor de Contorno da Geodésia (PVCG) na opção do problema vetorial fixo, de maneira simplificada. Foi também estimado o afastamento entre o Datum Vertical Brasileiro de Imbituba (DVB-I) e o Datum Vertical Brasileiro de Santana (DVB-S).

A solução do PVCG neste trabalho será o problema vetorial fixo simples, no qual a aproximação é linear, e as observações são referenciadas à superfície física. São conhecidas as coordenadas geodésicas φ,λ h e a gravidade g, e a incógnita é o geopotencial W. O problema é resolvido de acordo com a abordagem de Brovar (Brovar, 1964), usando superfície generalizada do potencial anômalo T. O diferencial neste caso é a independência das linhas de nivelamento; necessita-se apenas de dados gravimétricos com densidade adequada (Heck, 2011). Esta densificação precisa conduzir a uma resolução espacial determinada pela taxa de amostragem, que pode ser estabelecida em função da frequência de Nyquist - a melhor resolução espacial é igual à metade da frequência de amostragem.

Em relação aos Modelos Globais do Geopotencial, foi utilizado na solução do problema um MGG oriundo da missão GOCE, pois os modelos obtidos a partir dos dados gravimétricos do GOCE foram propostos por (Gerlach e Rummel, 2012) para servirem de base ao chamado Sistema Vertical de Referência Global (SVRG).

A realização do SVRG é possibilitada pela combinação de observações maregráficas, altimetria por radar a bordo de satélites, medidas terrestres, marinhas e aéreas de gravidadealém de técnicas espaciais, coordenadas ITRF e nivelamento geométrico (Ihde e Sánchez, 2005). Na atualidade, com o desenvolvimento das missões gravimétricas por satélite, os MGGs têm sido ferramentas essenciais na materialização dos SVRGs.

A modelagem das altas frequências do sinal gravitacional, o que na atualidade pode ser traduzido como a parcela do sinal gravitacional não contemplada pelos MGGs, é efetuada utilizando técnicas de decomposição espectral. Para a obtenção do sinal relativo aos curtos comprimentos de onda, duas técnicas são de interesse relevante: a clássica Decomposição Espectral ou Remove-Restore (RR) e a modelagem residual do terreno- Residual Terrain Model (RTM).

A técnica RTM consiste em calcular os curtos comprimentos de onda (altas frequências do campo da gravidade) não representados pelo MGG, utilizando uma superfície de referência (MDA) associada ao modelo do geopotencial; as massas topográficas acima da superfície são removidas e abaixo são preenchidas (Forsberg e Tscherning, 1981; Forsberg, 1984).

2. Materiais e Métodos

A área de estudos abrange parte dos Estados brasileiros do Amapá e do Pará, conforme mostra a Figura 1. Nesta figura também são apresentadas as estações que possuem dados de nivelamento e posicionamento GNSS utilizadas para validação dos MGGs, denominadas de estações GPS/NIV. As estações localizadas no Estado do Amapá têm origem no DVB-S, enquanto as estações do Estado do Pará são vinculadas ao DVB-I.

Figura 1
Área de estudos e estações GPS/NIV

Como dados para o trabalho, foram utilizados dois Modelos Globais do Geopotencial somente satélite do GOCE, GO_CONS_SPW_R4 e GO_CONS_TIM_R5, e os modelos combinados EIGEN-6C4 e EGM2008, disponibilizados pelo ICGEM (International Centre for Global Earth Models). As informações referentes aos modelos utilizados podem ser vistas na Tabela 1, na qual nos dados da terceira coluna (Fonte), S indica Satélite, G gravimetria e A altimetria.

Tabela 1:
Modelos Globais do Geopotencial utilizados

• FONTE: ICGEM (2015) RE - Resolução Espacial

Como forma de validar os modelos do geopotencial, foram utilizadas 17 estações GPS/NIV no Estado do Amapá (Tabela 2) e 15 estações no Estado do Pará (Tabela 3). Nestas tabelas são apresentadas as coordenadas, as altitudes elipsoidais e ortométricas (a rigor normais-ortométricas) conhecidas, com seus respectivos desvio-padrão. Para algumas estações, não constam no banco de dados do IBGE o desvio-padrão das altitudes ortométricas.

Tabela 2:
Estações GPS/NIV no Estado do Amapá (DVB-S)

• FONTE: IBGE,2014

Tabela 3:
Estações GPS/Niv no Estado do Pará (DVB-I)

• FONTE: IBGE(acesso em 2014)

2.1 Validação dos MGGs

Uma forma bastante usual de realizar a validação dos modelos do geopotencial é a comparação direta com dados levantados em campo, normalmente estações que possuem nivelamento e posicionamento GNSS associados.

Das estações GPS/NIV obteve-se a anomalia de altitude (ζ), que foi subtraída da anomalia de altitude obtida dos MGGs em seus graus máximos do desenvolvimento em séries de funções harmônicas esféricas, resultando nas diferenças das quais é possível inferir o erro absoluto (ε _absoluto).

Para obter a estimativa do erro relativo, foi calculada a diferença entre ζ_GPS/NIV-ζ_MGG e ζ_GPS/NIV -ζ_MGG+RTM e as distâncias entre estações; a distância foi considerada em concordância com a resolução espacial dos MGGs. Por exemplo, para o MGG GO_CONS_TIM_R5, que tem resolução de 71km, foram utilizadas para o cálculo apenas estações cujas distâncias entre si fossem iguais ou maiores que a resolução do modelo.

2.2 Compatibilização dos sistemas de maré permanente

As coordenadas ITRFxx e o posicionamento obtido com o sistema GPS estão no sistema livre de maré (tide-free). Em relação ao sistema de altitudes brasileiro, não foram aplicadas correções de maré ao nivelamento, portanto as altitudes normais-ortométricas estão no sistema médio de maré (mean-tide).

É necessário, portanto, realizar a compatibilização dos sistemas de maré aos dados utilizados. As funcionais dos MGGs foram obtidas no sistema livre de maré.. Para converter as altitudes normais-ortométricas das estações GPS/NIV do sistema médio de maré para o sistema livre de maré utilizou-se a formulação apresentada por (Tenzer et al., 2010):

k e h são os números de Love relativos à maré: k = 0.3 e h = 0.6

φ é a latitude geocêntrica

2.3 Abordagem do PVCG na forma fixada

Como solução do problema da obtenção de altitudes com significado físico sem dependência de rede de nivelamento, será abordada a solução do tipo Brovar para o PVCG na forma fixada, com a determinação de números geopotenciais baseados na modelagem do potencial anômalo.

O número geopotencial C_P é dado pela diferença entre o geopotencial no geoide e o geopotencial no ponto P considerado. A integral é independente do caminhamento utilizado (Torge, 2001 p. 80):

Com esta abordagem, (Heck, 2011) afirma que na solução fixada do PVCG o número geopotencial C_P não depende de nivelamento e gravimetria ao longo das linhas de nivelamento, mas é derivado da solução do PVCG fixo para qualquer ponto na superfície terrestre.

Supondo que na região de interesse houvesse informação gravimétrica densa o suficiente (por exemplo, um ponto a cada 100km², o que corresponderia a uma resolução espacial de 20km, considerada a frequência de Nyquist, o que implica em grau/ordem aprox. 1000), a solução completa da técnica Remove-Restore, para a abrangência dos longos e curtos comprimentos de onda do campo da gravidade, poderia ser dada da forma apresentada na Tabela 4.

Tabela 4:
Exemplo da solução completa da técnica Remove-Restore

Entretanto, na área de estudos existem somente 8 estações que contêm dados de posicionamento GPS associado com gravimetria, e não estão bem distribuídas, concentrando-se apenas em uma parte do Estado do Amapá. Desta forma, não é possível inserir estas observações na solução, uma vez que distorceriam os resultados. Assim sendo, para a obtenção do número geopotencial em uma região isolada, a solução se resumiu na utilização de MGG do GOCE + RTM, pois não foi possível obter o potencial anômalo T observado, ou seja, baseado em distúrbios observados da gravidade δg.

Na aplicação da técnica Remove-Restore, portanto, as altas frequências do campo da gravidade são abrangidas pelo RTM. Tomando como exemplo a Tabela 4, os graus acima da resolução do MGG GOCE são preenchidos pelo RTM. Excluindo-se a parte residual (T_RES ), a formulação da técnica Remove-Restore, considerando o potencial anômalo T, tem a seguinte forma:

O subíndice P indica que os valores são obtidos para um determinado ponto na superfície física. O potencial anômalo T representa as irregularidades regionais e locais de W, correspondendo a aproximadamente 0,1% do sinal gravitacional real , enquanto o potencial normal U representa cerca de 99,9% do sinal gravitacional real W. Desta forma, o potencial anômalo é dado por (Torge, 2001 p.214):

Da fórmula de Bruns N = T γ /(Hofmann-Wellenhof e Moritz, 2005 p. 93), obtém-se a fórmula análoga para o cálculo das anomalias de altitude ζ.

Como se conhece o valor da anomalia de altitude obtida do cálculo do RTM e a gravidade normal γ - Equações 11 e 12 - pode-se calcular o valor de T_RTM da seguinte forma:

A gravidade normal γ para uma dada latitude pode ser calculada pela fórmula de Somigliana - Equação 9 (Torge, 2001 p. 106).

Sendo γ_e a gravidade normal no equador, γ_P a gravidade normal no polo, a o semi-eixo maior do elipsoide de referência e b o semi-eixo menor.

O valor de T_MGG é calculado também pela fórmula de Bruns, aplicando-se os valores de anomalia de altitude obtidos diretamente do MGG.

O potencial normal U_P no ponto da superfície física é dado por (Hofmann-Wellenhof e Moritz, 2005 p. 93):

como

U₀ é o potencial normal, h é a altitude geodésica e γ a gravidade normal.

O potencial normal no elipsoide de referência é uma constante física. Para o SGR80, U₀ = 62.636.860,850 m²s^-2 (Hofmann-Wellenhof e Moritz, 2005 p. 86). Entretanto, por convenção, faz-se W₀ = U₀. Será utilizado o valor de W₀ = 62.636.856,0 ± 0,5 m²s^-2, divulgado por (Burša et al., 2004) e bastante aceito pela comunidade científica.

2.3.1 Comparação entre o geopotencial global e local

Para indicar o geopotencial de acordo com a superfície de referência considerada, são adotadas neste trabalho as seguintes nomenclaturas: W₀ = geopotencial global (geoide); W_i = geopotencial no Datum local (neste caso DVB-I) e W_P = geopotencial no ponto P da superfície física.

O geopotencial no Datum é dado pela equação:

onde γ_DVB é a gravidade normal no Datum e δH₀ é a discrepância vertical entre o geopotencial do DVB-I e o geopotencial global. (Ferreira et al., 2010) realizaram a determinação da diferença de potencial e chegaram a δW = -0,31m ± 0,01m, ou diferença de potencial δW = 3,04 ± 0,08m²s^-2, ou seja, o geopotencial no DVB-I está abaixo do geoide global do EGM2008 em 31cm (δH₀ = 31cm). Obteve-se assim um valor de geopotencial W_i no DVB-I igual a 62.636.859,04 m²s^-2.

Na obtenção dos números geopotenciais, C_P é o número geopotencial com caráter global, pois relaciona o geoide global W₀ e o geopotencial W_P , que tem referência à superfície global por estar diretamente relacionado a um MGG (C_P = W₀ - W_P ). Para calcular o número geopotencial local C_i , foi utilizado o geopotencial W_i referente ao Datum de Imbituba(C_{i IMB} =W_i - W_{P IMB} ). Para que o geopotencial W_P no ponto P tenha referência ao geoide local, subtrai-se a diferença de geopotencial δW, obtendo-se o valor de W_{P IMB} .

Para obter os números geopotenciais com caráter local das estações que têm origem no DVB-S (C_{i STN} ) é preciso considerar a discrepância existente entre os dois Data (δW_IMB-STN ), determinado neste trabalho e igual a 1,30m para o MGG GO_CONS_SPW_R4.

Também foi possível, aplicando a Equação 11, obter um valor de geopotencial W₀ no DVB-S; constatou-se que o geoide local no DVB-S está 0,99m acima do geoide global adotado nesta pesquisa. O valor obtido para W_i em Santana foi de 62.636.846,29m²s^-2.

A discrepância entre o Datum Vertical local (geopotencial W_i ), e o Datum Vertical global (geopotencial W₀ ) é dada em termos dos números geopotenciais na forma:

Os resultados das discrepâncias em termos dos números geopotenciais são apresentados na Tabela 3, tanto para o DVB-S quanto para o DVB-I.

3. Resultados

Antes de proceder à validação dos MGGs, efetuou-se a compatibilização dos sistemas de maré permanente. A conversão do sistema de maré médio de maré para sistema livre de maré resultou em 6,7cm, ou seja, as altitudes normais-ortométricas foram elevadas em média 6,7cm; estas altitudes corrigidas foram utilizadas na validação dos modelos.

3.1 Erro absoluto

A média das diferenças (ζ_GPS/NIV -ζ_MGG) são mostradas na Figura 2. Verifica-se que a maioria dos modelos tem comportamento semelhante, com diferenças médias de +0,65m para as estações do DVB-S e -0,68m para as estações do DVB-I. Fazendo a diferença das médias entre as estações com vínculo no DVB-S e as estações com vínculo no DVB-I, obtém-se a estimativa de afastamento entre os dois Data. Considerando o MGG GO_CONS_SPW_R4 até o grau 280, a diferença é de 1,30m. Incorporando a solução RTM, o afastamento é de 1,32m. O resultado obtido confirma a diferença existente entre os dois Data Verticais, conforme detectado por (Montecino, 2011).

É importante ressaltar que nas diferenças encontradas estão embutidas incertezas inerentes à precisão da própria rede de nivelamento, à limitação da precisão do posicionamento GNSS, à resolução alcançada pelo MGG, entre outros; outro fator relevante é que no Estado do Amapá os dados de nivelamento apresentam inconsistências, pois a rede de nivelamento não forma circuitos e portanto as altitudes não são ajustadas.

Figura 2:
ζ_GPS/NIV-ζ_MGG - Análise da Média

A qualidade dos modelos é avaliada com o erro médio quadrático (RMS - Root Mean Square). Para que o RMS reflita a qualidade real dos modelos, deve-se retirar o efeito sistemático, traduzido pela média das diferenças ζ_GPS/NIV- ζ_MGG. Desta forma, tem-se que a precisão dos modelos (p) é igual ao efeito sistemático (ē), calculado com a diferença de cada valor em relação à média, subtraído do erro absoluto (e):

O RMS para cada um dos Data, isento do efeito sistemático, é apresentado na Figura 3.

Figura 3:
ζ_GPS/NIV-ζ_MGG - Análise do RMS

3.2 Erro relativo

Na Tabela 5 são apresentados os erros relativos, considerados os graus máximos de cada MGG, para cada um dos Data. Foram consideradas estações aleatórias, com distâncias entre si superiores à resolução espacial dos modelos globais.

As mesmas observações descritas na seção 4.1.1 a respeito das incertezas existentes nos levantamentos e nos MGGs devem ser consideradas na análise do erro relativo.

Tabela 5:
Avaliação do Erro relativo DVB-I: ζ_GPS/NIV-ζ_MGG (cm/km)

Em geral, os MGGs apresentaram um bom desempenho relativo. Destaca-se que o EIGEN-6C4 teve o melhor resultado para a maioria das linhas de base consideradas, provavelmente porque este modelo possui os dados de gradiometria do GOCE, além dos dados de campo.

Para as estações vinculadas ao DVB-I, constata-se que para a menor linha de base (82,4km), o erro relativo é maior, chegando a mais que 1cm/km para os MGGs somente-satélite. Entretanto, em sua maioria, os MGGs também têm erros relativos aceitáveis.

3.3 Resultado do RTM

As elevações RTM foram construídas como a diferença entre as elevações do GMRT (Global Multi-Resolution Topography), que tem resolução espacial de 30m no terreno e 100m nas áreas oceânicas, e um MDA desenvolvido em séries de funções harmônicas esféricas, o ETOPO1, com resolução de 1' de arco e desenvolvido até o grau e ordem 2250. O RTM foi calculado truncando a topografia do ETOPO1 e os MGGs nos graus e ordem 180, 210, 250, 360 e graus máximos dos modelos globais considerados. Na Tabela 6 são mostrados os resultados para os diferentes graus e o efeito RTM médio para as estações em cada um dos Data. Observa-se que o efeito, em geral, é milimétrico, sendo que os maiores valores obtidos são para o grau e ordem 180 de desenvolvimento em séries de funções harmônicas esféricas.

Tabela 6:
Resultado do RTM - Média para as estações

Na análise da qualidade dos modelos, o efeito do RTM foi incorporado ao RMS para verificar qual a contribuição do RTM, retirado o efeito sistemático. De forma geral, a adição do efeito RTM à diferença ζ_GPS/NIV- ζ_MGG não melhorou a discrepância existente entre as duas técnicas, provavelmente devido à topografia plana da região de estudos, com pouca variação no relevo. Em regiões de relevo acentuado, há experimentos que comprovam a eficácia do RTM no preenchimento das feições relativas à topografia (Hirt et al., 2010).

Para algumas estações a aplicação do efeito do terreno residual não teve o resultado esperado, ou seja, houve piora na solução, provavelmente devido às imperfeições do modelo batimétrico disponível para as análises. Nas estações em que houve melhora, esta não foi significativa.

3.4 Obtenção dos números geopotenciais

Para a obtenção dos números geopotenciais foram aplicadas as formulações e a metodologia apresentados na Seção 3.4. O MGG utilizado foi o somente-satélite do GOCE, GO_CONS_SPW_R4, do ano de 2014. Os testes de validação dos MGGs não mostraram diferenças significativas entre os modelos, porém foi escolhido este MGG por ser o mais recente do GOCE e devido à metodologia de cálculo empregada, baseada em grades. Destaca-se que o erro de omissão deste modelo é muito próximo ao existente para os demais modelos do GOCE (30cm), uma vez que o erro de omissão depende principalmente do grau máximo dos harmônicos esféricos considerados (Reguzzoni, 2015 - comunicação pessoal). O MGG foi utilizado no grau e ordem 210 de desenvolvimento em séries de funções harmônicas esféricas, pois se expandido até o grau 280, não existe melhora significativa na modelagem do campo da gravidade pelo modelo.

Os resultados para todas as estações são mostrados na Tabela 7. A unidade de medida dos números geopotenciais é o kilogal x metro (kGal.m), que equivale a 10m²s^-2; para derivar altitudes a partir dos números geopotenciais, divide-se por um valor da gravidade que depende do tipo de altitude que se deseja determinar. Na referida tabela é apresentada também a diferença obtida nos números geopotenciais com vínculo no geopotencial global e local (C_P - C_i).

Tabela 7:
Números geopotenciais obtidos

A média calculada pelas diferenças entre os números geopotenciais referidos ao Datum Vertical local (C_i ) e os com caráter global (C_P ), derivados da solução do PVCG, reflete a discrepância existente entre a rede de altitudes oficial do país e o sistema proposto, materializado pelo GOCE (δW=W₀-W_i ou, analogamente, δW=C_P-C_i ), para ambos os Data Verticais de Santana e Imbituba.

Na comparação entre C_P e C_i , obteve-se um afastamento médio (δW) de -0,56m para as estações com vínculo no DVB-I, e de +0,72m para as estações com vínculo no DVB-S. Este resultado é coerente com a solução ζ_GPS/NIV - ζ_MGG (erro absoluto, Seção 4.1), com a qual obteve-se um off-set de 1,30m entre os dois DVs. Com a solução dos números geopotenciais, resulta uma discrepância de 1,31m. Salienta-se que no afastamento δW obtido estão inseridos o erro de omissão do MGG usado na solução, além de efeitos locais provenientes do cálculo do W_i . Neste caso, o geopotencial local W_i está associado ao EGM2008, cuja acurácia não é uniforme globalmente, mas depende dos dados terrestres de altimetria e gravimetria disponíveis.

O afastamento maior obtido para as estações com vínculo no DVB-S, quando comparado ao obtido para as estações com vínculo no DVB-I, pode ser atribuído às inconsistências existentes nas estações do SGB no Estado do Amapá.

3.4.1 Altitudes obtidas a partir dos números geopotenciais

A partir dos números geopotenciais pode-se derivar diferentes tipos de altitudes físicas, dependendo do tipo de gravidade empregada (e.g gravidade real, gravidade normal, gravidade média). De acordo com os dados disponíveis, e a partir dos números geopotenciais calculados, obteve-se as altitudes dinâmicas e normais para todas as estações da área de estudo. Não é possível obter as altitudes ortométricas, pois estas requerem a formulação de hipóteses simplificativas sobre a distribuição das massas internas da terra. Desta forma, calculou-se as altitudes de Helmert para 8 estações com vínculo no DVB-S e 2 estações com vínculo no DVB-I, pois somente estas estações possuem valores de gravidade observada (Tabela 8). Utilizou-se a redução de Poincaré-Prey para obter o valor de gravidade média (g_m ).

Tabela 8:
Altitudes calculadas a partir dos números geopotenciais

4. Conclusões e recomendações

A validação dos modelos do geopotencial com base em estações GPS/NIV mostrou que, em geral, os resultados foram similares para todos os MGGs. Em termos da média, as diferenças ζ_GPS/NIV - ζ_MGG para as estações com vínculo no DVB-S foram de +0,65m, e para as estações relativas ao DVB-I estas diferenças foram de -0,68m. Desta forma, considerando o MGG somente satélite do GOCE, GO_CONS_SPW_R4, pode-se estimar que a rede vinculada ao DVB-S tem sua origem 1,30m acima do segmento vinculado ao DVB-I.

Em relação à aplicação da técnica RTM para o preenchimento dos comprimentos de onda não abrangidos pelos MGGs, em geral o RTM não teve contribuição significativa; este comportamento pode ser justificado pelo relevo da área de estudos, praticamente sem variações acentuadas e com altitudes próximas ao nível do mar.

A obtenção dos números geopotenciais no contexto do PVCG na forma fixa mostrou-se bastante adequada, apesar da simplificação da solução devido à falta de dados gravimétricos. A metodologia proposta pode ser aplicada com sucesso principalmente em regiões com escassez de dados de nivelamento. Como não foi possível desenvolver a solução completa do PVCG por falta de dados na área de estudo, recomenda-se a aplicação da solução completa utilizando a gravimetria que está sendo executada pelo Projeto Amapá - Revitalização das Redes Geodésicas do Estado do Amapá, quando do término deste.

O emprego dos números geopotenciais é vantajoso porque a partir destes é possível derivar diferentes tipos de altitudes, dependendo das necessidades e das convenções estabelecidas. No contexto do SIRGAS, por exemplo, o princípio básico estabelecido pelo Grupo de Trabalho III - Datum Vertical - é o cálculo dos números geopotenciais referidos à mesma superfície equipotencial W₀ , e a partir daí cada país deriva altitudes físicas de acordo com suas necessidades.

Os números geopotenciais calculados com caráter global estão referenciados a um MGG oriundo do GOCE, que tem sido apontado para servir de base a um SVRG e também na questão da unificação de sistemas de altitudes. Conclui-se assim a eficácia da utilização dos números geopotenciais no contexto do PVCG para definição e realização de um sistema de altitudes local, em local remoto sem rede de nivelamento existente, vinculado ao SGB e ao SVRG.

REFERÊNCIAS BIBLIOGRÁFICAS

• Brovar, Vsevolod V. "On the solutions of Molodensky's boundary value problem". Bulletin géodésique. v. 72, Issue 1, p 167-173, June 1964.
• Burša, M., Kenyon, S., Kouba, J., Šíma, Z., Vatrt, V., Vojtíšková, M. "A global vertical reference frame based on four regional vertical datums". Studia Geophysica et Geodaetica. v. 48, p. 493 - 502, Praga, 2004.
• De Freitas, Silvio. R. C., e Blitzkow, Denizar. "Altitudes e Geopotencial". IGeS Bulletin N°.9 - International Geoid Service, p. 47 - 62, Milan, June 1999.
• Ferreira, Vagner G., De Freitas, Silvio R.C. e Heck, Bernhard. "Determinação da diferença de geopotencial do Datum Vertical Brasileiro relativa ao Earth Gravity Model 2008". Revista Brasileira de Cartografia nº 62/03, 2010.
• Forsberg, Rene. "A study of terrain reductions, density anomalies and geophysical inversion methods in gravity field modelling". The Ohio State University, Scientific Report n° 5, April 1984.
• Forsberg, Rene. and Tscherning, Carl C. "The Use of Height Data in Gravity Field Approximation by Collocation". Journal of Geophysical Research, v. 86, nºB9, p. 7843-7854, September 10, 1981.
• Gerlach, Christian, and Rummel, Reiner. "Global height system unification with GOCE: a simulation study on the indirect bias term in the GBVP approach". Journal of Geodesy; Springer, 2012.
• Heck, Bernhard. "A Brovar-Type Solution Of The Fixed Geodetic Boundary Value Problem". Studia Geophysica et Geodaetica. v. 55, Issue 3, p. 441-454, July 2011.
• Hirt, Christian., Featherstone, W.E., and Marti, U. "Combining EGM2008 and SRTM/DTM2006.0 residual terrain model data to improve quasi-geoid computations in mountainous areas devoid of gravity data"..Journal of Geodesy 84:557-567. doi 10.1007/s00190-010-0395-1. Springer-Verlag, 2010.
  » https://doi.org/10.1007/s00190-010-0395-1
• Hofmann-Wellenhof, B., and Moritz, H. Physical Geodesy. Springer WienNewYork, 2005.
• IBGE 2014 - Instituto Brasileiro de Geografia e Estatística - Banco de Dados Geodésicos. Acesso em 2014. http://www.ibge.gov.br/home/geociencias/geodesia/sgb.shtm
  » http://www.ibge.gov.br/home/geociencias/geodesia/sgb.shtm
• ICGEM 2015 - International Centre for Global Earth Models - Calculation of Gravity Field Functionals on Ellipsoidal Grids Acesso em 2015. http://icgem.gfz-potsdam.de/ICGEM/
  » http://icgem.gfz-potsdam.de/ICGEM/
• Ihde, Johannes, and Sánchez, Laura. "A unified global height reference system as a basis for IGGOS". Journal of Geodinamics. v. 40,, Issues 4-5 p. 400-413, November-December 2005.
• Montecino, Henry. "Proposta de alternativas para conexão dos Data Verticais Brasileiros de Imbituba e Santana". Dissertação de mestrado, Universidade Federal do Paraná, Curitiba, 2011.
• Reguzzoni, Mirko. Comunicação pessoal, 2015.
• Rummel, Reiner. GIM International Interview, February, 2013.
• Sánchez, Laura. "Untersuchung Physikalisher Höhen In Kolumbien" - Determinación de Alturas Físicas En Colombia (Tradução). Institut Für Planetare Geodäsie Technische Universität Dresden. Dresden, November 2002.
• Tenzer, Robert, Vatrt, V., Abdalla, A. and Dayoub, N. "Assessment of the LVD offsets for the normal-orthometric heights and different permanent tide systems - a case study of New Zealand". Appl. Geomat. 3: 1 - 8. doi 10.1007/s12518-010-0038-5. Springer, 2010.
  » https://doi.org/10.1007/s12518-010-0038-5
• Torge, Wofgang. Geodesy. 3rd completely rev. and extended ed. Berlin; New York: Walter de Gruyter, 2001.

Datas de Publicação

• Publicação nesta coleção
  Jun 2016

Histórico

• Recebido
  Jun 2015
• Aceito
  Set 2015