Improving Surveying Accuracy and Efficiency in Connecticut: An Accuracy Assessment of GEOID03 and GEOID09

Comparing published NAVD 88 Helmert orthometric heights of First-Order bench marks against GPS-determined orthometric heights showed that GEOID03 and GEOID09 perform at their reported accuracy in Connecticut. GPS-determined orthometric heights were determined by subtracting geoid undulations from ellipsoid heights obtained from a network least-squares adjustment of GPS occupations in 2007 and 2008. A total of 73 markers were occupied in these stability classes: 25 class A, 11 class B, 12 class C, 2 class D bench marks, and 23 temporary marks with transferred elevations. Adjusted ellipsoid heights were compared against OPUS as a check. We found that: the GPS-determined orthometric heights of stability class A markers and the transfers are statistically lower than their published values but just barely; stability class B, C and D markers are also statistically lower in a manner consistent with subsidence or settling; GEOID09 does not exhibit a statistically significant residual trend across Connecticut; and GEOID09 out-performed GEOID03. A "correction surface" is not recommended in spite of the geoid models being statistically different than the NAVD 88 heights because the uncertainties involved dominate the discrepancies. Instead, it is recommended that the vertical control network be re-observed.

What Does Height Really Mean? Part III: Height Systems

This is the third paper in a four-part series considering the fundamental question, “what does the word “height” really mean?” The first paper reviewed reference ellipsoids and mean sea level datums. The second paper reviewed the physics of heights culminating in a simple development of the geoid and explained why mean sea level stations are not all at the same orthometric height. This third paper develops the principle notions of height, namely measured, differentially deduced changes in elevation, orthometric heights, Helmert orthometric heights, normal orthometric heights, dynamic heights, and geopotential numbers. We conclude with a more in-depth discussion of current thoughts regarding the geoid.

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

Sediment volume variations inferred from DGPS measurements within Bossons glacier proglacial area (Mont-Blanc massif, France)

This dataset presents Differential Global Positioning System data (DGPS) acquired within the Bossons glacier proglacial area. Bossons glacier is a rapidly retreating glacier and its proglacial area is deglaciated for ~30 years. Bossons stream is one of the outlets of the subglacial drainage system. It starts as a 800 m steep cascade reach, then flows through an area with gentler slope : the Plan des Eaux (PdE). PdE is a 300 m long, 50 m wide proglacial alluvial plain with an increasing channel mobility in the downstream direction but decreasing slope gradient and incision. As it may act a sediment trap, studying periglacial and proglacial erosion processes in the Bossons catchment requires to quantify PdE sediment volume evolution. A several meter-sized block located within Bossons proglacial area was set up as GPS base : its location was measured by one antenna (Topcon Hyper Pro) by performing 600 consecutive measurements throughout one day. A second antenna (Topcon Hyper Pro) was then used to measure XYZ location of points in the proglacial area with a ~2 m grid. Radio communication between the two antennas allowed differential calculations to be automatically carried out on field using the Topcon FC-250 hand controller. This methodology yields 3 cm XY and 1.5 cm Z uncertainties. DGPS data have been acquired through 10 campaigns from 2004 to 2014; campaigns from 2004 to 2008 cover a smaller area than those from 2010 to 2014. Digital Elevation Model (DEM) have been interpolated from DGPS data and difference between two DEMs yields deposited and eroded volume within PdE. Maps of PdE volume variation between two campaigns show that incision mainly occurs in the upper and lower sections where as deposition dominates in the middle section. Deposition, denudation and net rate (deposition rate - denudation rate) are calculated by normalizing volumes by DEM areas. Deposition dominates results with a mean net rate of 29 mm/yr. However, strong inter-annual variability exists and some years are dominated by denudation : -36 mm/yr and -100 mm/yr for 2006 and 2011, respectively. Nonetheless, oldest campaigns (2004 to 2008) were carried out on the lower part part of the alluvial plain and ruling them out to keep only complete DEM (2010 to 2014) yields a mean net rate of ~15 mm/yr. This results is coherent with field observations of both strong deposition (e.g. flood deposits) and strong erosion (e.g. 30 cm incision) evidences. Bossons glacier proglacial area is thus dynamic with year-to-year geormorphological changes but may leans toward increasing its mean elevation through a deposition dominated system.

TanDEM-X ice thickness dataset for Wilkins Ice Shelf, Antarctic Peninsula, with link to GeoTIFF

The dataset shows the ice thickness over Wilkins Ice Shelf, Antarctic Peninsula derived from TanDEM-X Interferometry. The data has been acquired between June and August 2012. The TanDEM-X heights have been linked to CryoSAT-2 heights (V. Helm) from the respective time stamp. Elevations have been transformed from WGS84 ellipsoidal heights to the EGM2008 geoid. The ice shelf thickness was estimated assuming hydrostatic equilibrium and a mean ice density of 915 kg/m³.

Coastal Urban Flood Resislience versus Climate Change

Processes of founding and expanding cities in coastal areas have undergone great changes over time driven by environmental conditions. Coastal settlements looked for places above flood levels and away from swamps and other wetlands whenever possible. As populations grew, cities were extending trying to avoid low and wet lands. No city has been able to limit its growth. The risk of flooding can never be eliminated, but only reduced to the extent possible. Flooding of coastal areas is today dramatically attributed to eustasic sea level rise caused by global climate change. This can be inaccurate. Current climate change is generating an average sea level upward trend, but other regional and local factors result in this trend being accentuated in some places or attenuated, and even reversed, in others. Then, the intensity and frequency of coastal flooding around the planet, although not so much as a unique result of this general eustasic elevation, but rather of the superposition of marine and crustal dynamic elements, the former also climate-related, which give rise to a temporary raising in average sea level in the short term. Since the Little Ice Age the planet has been suffering a global warming change leading to sea level rise. The idea of being too obeying to anthropogenic factors may be attributed to Arrhenius (1896), though it is of much later highlight after the sixties of the last century. Never before, the human factor had been able of such an influence on climate. However, other types of changes in sea levels became apparent, resulting from vertical movements of the crust, modifications of sea basins due to continents fracturing, drifting and coming together, or to different types of climate patterns. Coastal zones are then doubly susceptible to floods. Precipitation immediately triggers pluvial flooding. If it continues upland or when snow and glaciers melt eventually fluvial flooding can occur. The urban development presence represents modifying factors. Additional interference is caused by river and waste water drainage systems. Climate also influences sea levels in coastal areas, where tides as well as the structure and dynamic of the geoid and its crust come into play. From the sea, waters can flood and break or push back berms and other coastline borders. The sea level, controlling the mouth of the main channel of the basin's drainage system, is ultimately what governs flood levels. A temporary rise in sea level acts as a dam at the mouth. Even in absence of that global change, so, floods are likely going to increase in many urban coastal areas. Some kind of innovative methodologies and practices should be needed to get more flood resilience cities

Comparación de la precisión de los Modelos Geopotenciales EGM96 y EGM08 en la zona del Caribe.

Recientemente fue publicado un nuevo modelo geopotencial global, el EGM08. Este modelo ha mostrado una notable mejoría en la calidad de sus tres fuentes de datos; las observaciones del movimiento perturbado de los satélites artificiales,la altimetría por satélite y la gravimetría terrestre, por lo que se ha conseguido mejorar su precisión. En Puerto Rico, nuestra área de estudio, encontramos que al comprar las diferencias de los incrementos de la ondulación del geoide geométrico (calculado con medidas de campo) con los valores de los incrementos de la ondulación del geoide obtenidos utilizando estos modelos geopotenciales globales, la precisión del EGM08 fue ± 0,029 metros mientras que la precisión del EGM96 fue ± 0,055 metros. Estos resultados demuestran que en nuestra región, el modelo EGM08 ha presentado una mejoría considerable sobre su predecesor el EGM96 al momento de determinar los valores de los incrementos de la ondulación del geoide. Abstract: Recently, the new global geopotential model, the EGM08 was published. This model has shown a marked improvement in the quality of its three sources of data; the observations of the disturbed motion of artificial satellites, satellite altimetry and terrestrial gravity, so it has improved its precision. In our study area, Puerto Rico, we found that when we compare the differences of the increments of the geometric geoid undulation (computed with field data) with the values of the increments of the geoid undulation obtained using these models, the EGM08 accuracy was ± 0,029 meters, while the EGM96 accuracy was ± 0,055 meters. These results confirm that in our region, the EGM08 model has presented a significant improvement over its predecessor the EGM96 when determining the values of the increments of the geoid undulation.

Determinación de un modelo de geoide gravimétrico para Puerto Rico como sistema de referencia para las altitudes ortométricas

El geoide, definido como la superficie equipotencial que mejor se ajusta (en el sentido de los mínimos cuadrados) al nivel medio del mar en una determinada época, es la superficie que utilizamos como referencia para determinar las altitudes ortométricas. Si disponemos de una superficie equipotencial de referencia como dátum altimétrico preciso o geoide local, podemos entonces determinar las altitudes ortométricas de forma eficiente a partir de las altitudes elipsoidales proporcionadas por el Sistema Global de Navegación por Satélite (Global Navigation Satellite System, GNSS ). Como es sabido uno de los problemas no resueltos de la geodesia (quizás el más importante de los mismos en la actualidad) es la carencia de un dátum altimétrico global (Sjoberg, 2011) con las precisiones adecuadas. Al no existir un dátum altimétrico global que nos permita obtener los valores absolutos de la ondulación del geoide con la precisión requerida, es necesario emplear modelos geopotenciales como alternativa. Recientemente fue publicado el modelo EGM2008 en el que ha habido una notable mejoría de sus tres fuentes de datos, por lo que este modelo contiene coeficientes adicionales hasta el grado 2190 y orden 2159 y supone una sustancial mejora en la precisión (Pavlis et al., 2008). Cuando en una región determinada se dispone de valores de gravedad y Modelos Digitales del Terreno (MDT) de calidad, es posible obtener modelos de superficies geopotenciales más precisos y de mayor resolución que los modelos globales. Si bien es cierto que el Servicio Nacional Geodésico de los Estados Unidos de América (National Geodetic Survey, NGS) ha estado desarrollando modelos del geoide para la región de los Estados Unidos de América continentales y todos sus territorios desde la década de los noventa, también es cierto que las zonas de Puerto Rico y las Islas Vírgenes Estadounidenses han quedado un poco rezagadas al momento de poder aplicar y obtener resultados de mayor precisión con estos modelos regionales del geoide. En la actualidad, el modelo geopotencial regional vigente para la zona de Puerto Rico y las Islas Vírgenes Estadounidenses es el GEOID12A (Roman y Weston, 2012). Dada la necesidad y ante la incertidumbre de saber cuál sería el comportamiento de un modelo del geoide desarrollado única y exclusivamente con datos de gravedad locales, nos hemos dado a la tarea de desarrollar un modelo de geoide gravimétrico como sistema de referencia para las altitudes ortométricas. Para desarrollar un modelo del geoide gravimétrico en la isla de Puerto Rico, fue necesario implementar una metodología que nos permitiera analizar y validar los datos de gravedad terrestre existentes. Utilizando validación por altimetría con sistemas de información geográfica y validación matemática por colocación con el programa Gravsoft (Tscherning et al., 1994) en su modalidad en Python (Nielsen et al., 2012), fue posible validar 1673 datos de anomalías aire libre de un total de 1894 observaciones obtenidas de la base de datos del Bureau Gravimétrico Internacional (BGI). El aplicar estas metodologías nos permitió obtener una base de datos anomalías de la gravedad fiable la cual puede ser utilizada para una gran cantidad de aplicaciones en ciencia e ingeniería. Ante la poca densidad de datos de gravedad existentes, fue necesario emplear un método alternativo para densificar los valores de anomalías aire libre existentes. Empleando una metodología propuesta por Jekeli et al. (2009b) se procedió a determinar anomalías aire libre a partir de los datos de un MDT. Estas anomalías fueron ajustadas utilizando las anomalías aire libre validadas y tras aplicar un ajuste de mínimos cuadrados por zonas geográficas, fue posible obtener una malla de datos de anomalías aire libre uniforme a partir de un MDT. Tras realizar las correcciones topográficas, determinar el efecto indirecto de la topografía del terreno y la contribución del modelo geopotencial EGM2008, se obtuvo una malla de anomalías residuales. Estas anomalías residuales fueron utilizadas para determinar el geoide gravimétrico utilizando varias técnicas entre las que se encuentran la aproximación plana de la función de Stokes y las modificaciones al núcleo de Stokes, propuestas por Wong y Gore (1969), Vanicek y Kleusberg (1987) y Featherstone et al. (1998). Ya determinados los distintos modelos del geoide gravimétrico, fue necesario validar los mismos y para eso se utilizaron una serie de estaciones permanentes de la red de nivelación del Datum Vertical de Puerto Rico de 2002 (Puerto Rico Vertical Datum 2002, PRVD02 ), las cuales tenían publicados sus valores de altitud elipsoidal y elevación. Ante la ausencia de altitudes ortométricas en las estaciones permanentes de la red de nivelación, se utilizaron las elevaciones obtenidas a partir de nivelación de primer orden para determinar los valores de la ondulación del geoide geométrico (Roman et al., 2013). Tras establecer un total de 990 líneas base, se realizaron dos análisis para determinar la 'precisión' de los modelos del geoide. En el primer análisis, que consistió en analizar las diferencias entre los incrementos de la ondulación del geoide geométrico y los incrementos de la ondulación del geoide de los distintos modelos (modelos gravimétricos, EGM2008 y GEOID12A) en función de las distancias entre las estaciones de validación, se encontró que el modelo con la modificación del núcleo de Stokes propuesta por Wong y Gore presentó la mejor 'precisión' en un 91,1% de los tramos analizados. En un segundo análisis, en el que se consideraron las 990 líneas base, se determinaron las diferencias entre los incrementos de la ondulación del geoide geométrico y los incrementos de la ondulación del geoide de los distintos modelos (modelos gravimétricos, EGM2008 y GEOID12A), encontrando que el modelo que presenta la mayor 'precisión' también era el geoide con la modificación del núcleo de Stokes propuesta por Wong y Gore. En este análisis, el modelo del geoide gravimétrico de Wong y Gore presento una 'precisión' de 0,027 metros en comparación con la 'precisión' del modelo EGM2008 que fue de 0,031 metros mientras que la 'precisión' del modelo regional GEOID12A fue de 0,057 metros. Finalmente podemos decir que la metodología aquí presentada es una adecuada ya que fue posible obtener un modelo del geoide gravimétrico que presenta una mayor 'precisión' que los modelos geopotenciales disponibles, incluso superando la precisión del modelo geopotencial global EGM2008. ABSTRACT The geoid, defined as the equipotential surface that best fits (in the least squares sense) to the mean sea level at a particular time, is the surface used as a reference to determine the orthometric heights. If we have an equipotential reference surface or a precise local geoid, we can then determine the orthometric heights efficiently from the ellipsoidal heights, provided by the Global Navigation Satellite System (GNSS). One of the most common and important an unsolved problem in geodesy is the lack of a global altimetric datum (Sjoberg, 2011)) with the appropriate precision. In the absence of one which allows us to obtain the absolute values of the geoid undulation with the required precision, it is necessary to use alternative geopotential models. The EGM2008 was recently published, in which there has been a marked improvement of its three data sources, so this model contains additional coefficients of degree up to 2190 and order 2159, and there is a substantial improvement in accuracy (Pavlis et al., 2008). When a given region has gravity values and high quality digital terrain models (DTM), it is possible to obtain more accurate regional geopotential models, with a higher resolution and precision, than global geopotential models. It is true that the National Geodetic Survey of the United States of America (NGS) has been developing geoid models for the region of the continental United States of America and its territories from the nineties, but which is also true is that areas such as Puerto Rico and the U.S. Virgin Islands have lagged behind when to apply and get more accurate results with these regional geopotential models. Right now, the available geopotential model for Puerto Rico and the U.S. Virgin Islands is the GEOID12A (Roman y Weston, 2012). Given this need and given the uncertainty of knowing the behavior of a regional geoid model developed exclusively with data from local gravity, we have taken on the task of developing a gravimetric geoid model to use as a reference system for orthometric heights. To develop a gravimetric geoid model in the island of Puerto Rico, implementing a methodology that allows us to analyze and validate the existing terrestrial gravity data is a must. Using altimetry validation with GIS and mathematical validation by collocation with the Gravsoft suite programs (Tscherning et al., 1994) in its Python version (Nielsen et al., 2012), it was possible to validate 1673 observations with gravity anomalies values out of a total of 1894 observations obtained from the International Bureau Gravimetric (BGI ) database. Applying these methodologies allowed us to obtain a database of reliable gravity anomalies, which can be used for many applications in science and engineering. Given the low density of existing gravity data, it was necessary to employ an alternative method for densifying the existing gravity anomalies set. Employing the methodology proposed by Jekeli et al. (2009b) we proceeded to determine gravity anomaly data from a DTM. These anomalies were adjusted by using the validated free-air gravity anomalies and, after that, applying the best fit in the least-square sense by geographical area, it was possible to obtain a uniform grid of free-air anomalies obtained from a DTM. After applying the topographic corrections, determining the indirect effect of topography and the contribution of the global geopotential model EGM2008, a grid of residual anomalies was obtained. These residual anomalies were used to determine the gravimetric geoid by using various techniques, among which are the planar approximation of the Stokes function and the modifications of the Stokes kernel, proposed by Wong y Gore (1969), Vanicek y Kleusberg (1987) and Featherstone et al. (1998). After determining the different gravimetric geoid models, it was necessary to validate them by using a series of stations of the Puerto Rico Vertical Datum of 2002 (PRVD02) leveling network. These stations had published its values of ellipsoidal height and elevation, and in the absence of orthometric heights, we use the elevations obtained from first - order leveling to determine the geometric geoid undulation (Roman et al., 2013). After determine a total of 990 baselines, two analyzes were performed to determine the ' accuracy ' of the geoid models. The first analysis was to analyze the differences between the increments of the geometric geoid undulation with the increments of the geoid undulation of the different geoid models (gravimetric models, EGM2008 and GEOID12A) in function of the distance between the validation stations. Through this analysis, it was determined that the model with the modified Stokes kernel given by Wong and Gore had the best 'accuracy' in 91,1% for the analyzed baselines. In the second analysis, in which we considered the 990 baselines, we analyze the differences between the increments of the geometric geoid undulation with the increments of the geoid undulation of the different geoid models (gravimetric models, EGM2008 and GEOID12A) finding that the model with the highest 'accuracy' was also the model with modifying Stokes kernel given by Wong and Gore. In this analysis, the Wong and Gore gravimetric geoid model presented an 'accuracy' of 0,027 meters in comparison with the 'accuracy' of global geopotential model EGM2008, which gave us an 'accuracy' of 0,031 meters, while the 'accuracy ' of the GEOID12A regional model was 0,057 meters. Finally we can say that the methodology presented here is adequate as it was possible to obtain a gravimetric geoid model that has a greater 'accuracy' than the geopotential models available, even surpassing the accuracy of global geopotential model EGM2008.

Steric sea-level variations inferred from combined Topex/Poseidon altimetry and GRACE gravimetry

In this study, we propose to estimate the steric sea-level variations over a < 2-year period (April 2002 through December 2003) by combining global mean sea level (GMSL) based on Topex/ Poseidon (T/P) altimetry with time-variable geoid averaged over the oceans, as observed by the GRACE (Gravity Recovery and Climate Experiment) satellite. In effect, altimetry-derived GMSL changes results from two contributions: Steric (thermal plus salinity) effects due to sea water density change and ocean mass change due to water exchange with atmosphere and continents. On the other hand, GRACE data over the oceans provide the ocean mass change component only. The paper first discusses the corrections to apply to the GRACE data. Then the steric contribution to the GMSL is estimated using GRACE and T/P data. Comparison with available thermal expansion based on in situ hydrographic data is performed.

Resonance studies of artificial earth satellites

Orbit determination from artificial satellite observations is a key process in obtaining information about the Earth and its environment. A study of the perturbations experienced by these satellites enables knowledge to be gained of the upper atmosphere, the gravity field, ocean tides, solid-Earth tides and solar radiation. The gravity field is expressed as a double infinite series of associated Legendre functions (tesseral harmonics). In contemporary global gravity field models the overall geoid is well determined. An independent check on these gravity field harmonics of a particular order may be made by analysis of satellites that pass through resonance of that order. For such satellites the perturbations of the orbital elements close to resonance are analysed to derive lumped harmonic coefficients. The orbital parameters of 1984-106A have been determined at 43 epochs, during which time the satellite was close to 14th order resonance. Analysis of the inclination and eccentricity yielded 6 lumped harmonic coefficients of order 14 whilst analysis of the mean motion yielded additional pairs of lumped harmonics of orders 14, 28 and 42, with the 14th order harmonics superseding those obtained from analysis of the inclination. This thesis concentrates in detail on the theoretical changes of a near-circular satellite orbit perturbed by the Earth's gravity field under the influence of minimal air-drag whilst in resonance with the Earth. The satellite 1984-106A experienced the interesting property of being temporarily trapped with respect to a secondary resonance parameter due to the low air-drag in 1987. This prompted the theoretical investigation of such a phenomenon. Expressions obtained for the resonance parameter led to the determination of 8 lumped harmonic coefficients, coincidental to those already obtained. All the derived lumped harmonic values arc used to test the accuracy of contemporary gravity field models and the underlying theory in this thesis.

Precise orbit determination and analysis from satellite altimetry and laser ranging

For optimum utilization of satellite-borne instrumentation, it is necessary to know precisely the orbital position of the spacecraft. The aim of this thesis is therefore two-fold - firstly to derive precise orbits with particular emphasis placed on the altimetric satellite SEASAT and secondly, to utilize the precise orbits, to improve upon atmospheric density determinations for satellite drag modelling purposes. Part one of the thesis, on precise orbit determinations, is particularly concerned with the tracking data - satellite laser ranging, altimetry and crossover height differences - and how this data can be used to analyse errors in the orbit, the geoid and sea-surface topography. The outcome of this analysis is the determination of a low degree and order model for sea surface topography. Part two, on the other hand, mainly concentrates on using the laser data to analyse and improve upon current atmospheric density models. In particular, the modelling of density changes associated with geomagnetic disturbances comes under scrutiny in this section. By introducing persistence modelling of a geomagnetic event and solving for certain geomagnetic parameters, a new density model is derived which performs significantly better than the state-of-the-art models over periods of severe geomagnetic storms at SEASAT heights. This is independently verified by application of the derived model to STARLETTE orbit determinations.

Satellite borne radar altimetry : empirical orbit refinement and ocean signal recovery techniques

Measurements of the sea surface obtained by satellite borne radar altimetry are irregularly spaced and contaminated with various modelling and correction errors. The largest source of uncertainty for low Earth orbiting satellites such as ERS-1 and Geosat may be attributed to orbital modelling errors. The empirical correction of such errors is investigated by examination of single and dual satellite crossovers, with a view to identifying the extent of any signal aliasing: either by removal of long wavelength ocean signals or introduction of additional error signals. From these studies, it was concluded that sinusoidal approximation of the dominant one cycle per revolution orbit error over arc lengths of 11,500 km did not remove a significant mesoscale ocean signal. The use of TOPEX/Poseidon dual crossovers with ERS-1 was shown to substantially improve the radial accuracy of ERS-1, except for some absorption of small TOPEX/Poseidon errors. The extraction of marine geoid information is of great interest to the oceanographic community and was the subject of the second half of this thesis. Firstly through determination of regional mean sea surfaces using Geosat data, it was demonstrated that a dataset with 70cm orbit error contamination could produce a marine geoid map which compares to better than 12cm with an accurate regional high resolution gravimetric geoid. This study was then developed into Optimal Fourier Transform Interpolation, a technique capable of analysing complete altimeter datasets for the determination of consistent global high resolution geoid maps. This method exploits the regular nature of ascending and descending data subsets thus making possible the application of fast Fourier transform algorithms. Quantitative assessment of this method was limited by the lack of global ground truth gravity data, but qualitative results indicate good signal recovery from a single 35-day cycle.

Application of ERS-1 and Topex/Poseidon altimetry using precise orbits

The work described in this thesis concerns the application of radar altimetry, collected from the ERS-1 and TOPEX/POSEIDON missions, to precise satellite orbits computed at Aston University. The data is analysed in a long arc fashion to determine range biases, time tag biases, sea surface topographies and to assess the radial accuracy of the generated orbits through crossover analysis. A sea surface variability study is carried out for the North Sea using repeat altimeter profiles from ERS-1 and TOPEX/POSEIDON in order to verify two local U.K. models for ocean tide and storm surge effects. An on-side technique over the English Channel is performed to compute the ERS-1, TOPEX and POSEIDON altimeter range biases by using a combination of altimetry, precise orbits determined by short arc methods, tide gauge data, GPS measurements, geoid, ocean tide and storm surge models. The remaining part of the thesis presents some techniques for the short arc correction of long arc orbits. Validation of this model is achieved by way of comparison with actual SEASAT short arcs. Simulations are performed for the ERS-1 microwave tracking system, PRARE, using the range data to determine time dependent orbit corrections. Finally, a brief chapter is devoted to the recovery of errors in station coordinates by the use of multiple short arcs.

Espessura da zona de transição do manto sob a Província Borborema, NE Brasil, através de conversões Ps telessísmicas

The mantle transition zone is defined by two seismic discontinuities, nominally at 410 and 660 km depth, which result from transformations in the mineral olivine. The topography of these discontinuities provides information about lateral temperature changes in the transition zone. In this work, P-to-S conversions from teleseismic events recorded at 32 broadband stations in the Borborema Province were used to determine the transition zone thickness beneath this region and to investigate whether there are lateral temperature changes within this depth range. For this analysis, stacking and migration of receiver functions was performed. In the Borborema Province, geophysical studies have revealed a geoid anomaly which could reflect the presence of a thermal anomaly related to the origin of intraplate volcanism and uplift that marked the evolution of the Province in the Cenozoic. Several models have been proposed to explain these phenomena, which include those invoking the presence of a deep-seated mantle plume and those invoking shallower sources, such as small-scale convection cells. The results of this work show that no thermal anomalies are present at transition zone depths, as significant variations in the transition zone thickness were not observed. However, regions of depressed topography for both discontinuities (410 and 660 km) that approximately overlap in space were identified, suggesting that lower-thanaverage, lateral variations in seismic velocity above 410 km depth may exist below the the Borborema Province. This is consistent with the presence of a thermally-induced, low-density body independently inferred from analysis of geoid anomalies. Therefore, the magma source responsible for the Cenozoic intraplate volcanism and related uplift in the Province, is likely to be confined above the upper mantle transition zone.