Application of altimetry data assimilation on mesoscale eddies simulation

Mesoscale eddy plays an important role in the ocean circulation. In order to improve the simulation accuracy of the mesoscale eddies, a three-dimensional variation (3DVAR) data assimilation system called Ocean Variational Analysis System (OVALS) is coupled with a POM model to simulate the mesoscale eddies in the Northwest Pacific Ocean. In this system, the sea surface height anomaly (SSHA) data by satellite altimeters are assimilated and translated into pseudo temperature and salinity (T-S) profile data. Then, these profile data are taken as observation data to be assimilated again and produce the three-dimensional analysis T-S field. According to the characteristics of mesoscale eddy, the most appropriate assimilation parameters are set up and testified in this system. A ten years mesoscale eddies simulation and comparison experiment is made, which includes two schemes: assimilation and non-assimilation. The results of comparison between two schemes and the observation show that the simulation accuracy of the assimilation scheme is much better than that of non-assimilation, which verified that the altimetry data assimilation method can improve the simulation accuracy of the mesoscale dramatically and indicates that it is possible to use this system on the forecast of mesoscale eddies in the future.

Characterising errors in airborne laser altimetry data to extract soil roughness

Airborne laser altimetry has the potential to make frequent detailed observations that are important for many aspects of studying land surface processes. However, the uncertainties inherent in airborne laser altimetry data have rarely been well measured. Uncertainty is often specified as generally as 20cm in elevation, and 40cm planimetric. To better constrain these uncertainties, we present an analysis of several datasets acquired specifically to study the temporal consistency of laser altimetry data, and thus assess its operational value. The error budget has three main components, each with a time regime. For measurements acquired less than 50ms apart, elevations have a local standard deviation in height of 3.5cm, enabling the local measurement of surface roughness of the order of 5cm. Points acquired seconds apart acquire an additional random error due to Differential Geographic Positioning System (DGPS) fluctuation. Measurements made up to an hour apart show an elevation drift of 7cm over a half hour. Over months, this drift gives rise to a random elevation offset between swathes, with an average of 6.4cm. The RMS planimetric error in point location was derived as 37.4cm. We conclude by considering the consequences of these uncertainties on the principle application of laser altimetry in the UK, intertidal zone monitoring.

Corrected ICESat altimetry data, surface mass balance, and firn elevation change on Antarctic ice shelves

Accurate prediction of global sea-level rise requires that we understand the cause of recent, widespread and intensifying glacier acceleration along Antarctic ice-sheet coastal margins. Floating ice shelves buttress the flow of grounded tributary glaciers and their thickness and extent are particularly susceptible to changes in both climate and ocean forcing. Recent ice-shelf collapse led to retreat and acceleration of several glaciers on the Antarctic Peninsula. However, the extent and magnitude of ice-shelf thickness change, its causes and its link to glacier flow rate are so poorly understood that its influence on the future of the ice sheets cannot yet be predicted. Here we use satellite laser altimetry and modelling of the surface firn layer to reveal for the first time the circum-Antarctic pattern of ice-shelf thinning through increased basal melt. We deduce that this increased melt is the primary driver of Antarctic ice-sheet loss, through a reduction in buttressing of the adjacent ice sheet that has led to accelerated glacier flow. The highest thinning rates (~7 m/a) occur where warm water at depth can access thick ice shelves via submarine troughs crossing the continental shelf. Wind forcing could explain the dominant patterns of both basal melting and the surface melting and collapse of Antarctic ice shelves, through ocean upwelling in the Amundsen and Bellingshausen Seas and atmospheric warming on the Antarctic Peninsula. This implies that climate forcing through changing winds influences Antarctic Ice Sheet mass balance, and hence global sea-level, on annual to decadal timescales.

Combination of altimetry data from different satellite missions

Substantial altimetry datasets collected by different satellites have only become available during the past five years, but the future will bring a variety of new altimetry missions, both parallel and consecutive in time. The characteristics of each produced dataset vary with the different orbital heights and inclinations of the spacecraft, as well as with the technical properties of the radar instrument. An integral analysis of datasets with different properties offers advantages both in terms of data quantity and data quality. This thesis is concerned with the development of the means for such integral analysis, in particular for dynamic solutions in which precise orbits for the satellites are computed simultaneously. The first half of the thesis discusses the theory and numerical implementation of dynamic multi-satellite altimetry analysis. The most important aspect of this analysis is the application of dual satellite altimetry crossover points as a bi-directional tracking data type in simultaneous orbit solutions. The central problem is that the spatial and temporal distributions of the crossovers are in conflict with the time-organised nature of traditional solution methods. Their application to the adjustment of the orbits of both satellites involved in a dual crossover therefore requires several fundamental changes of the classical least-squares prediction/correction methods. The second part of the thesis applies the developed numerical techniques to the problems of precise orbit computation and gravity field adjustment, using the altimetry datasets of ERS-1 and TOPEX/Poseidon. Although the two datasets can be considered less compatible that those of planned future satellite missions, the obtained results adequately illustrate the merits of a simultaneous solution technique. In particular, the geographically correlated orbit error is partially observable from a dataset consisting of crossover differences between two sufficiently different altimetry datasets, while being unobservable from the analysis of altimetry data of both satellites individually. This error signal, which has a substantial gravity-induced component, can be employed advantageously in simultaneous solutions for the two satellites in which also the harmonic coefficients of the gravity field model are estimated.

Thermal front variability along the North Atlantic Current observed using satellite microwave and infrared data

Thermal fronts detected using multiple satellite sensors have been integrated to provide new information on the spatial and seasonal distribution of oceanic fronts in the North Atlantic. The branching of the North Atlantic Current (NAC) as it encounters the Mid-Atlantic Ridge (MAR) is reflected in surface thermal fronts, which preferentially occur at the Charlie Gibbs Fracture Zone (CGFZ) and several smaller fracture zones. North of the CGFZ there are few thermal fronts, contrasting with the region to the south, where there are frequent surface thermal fronts that are persistent seasonally and interannually. The alignment of the fronts confirms that the shallower Reykjanes Ridge north of the CGFZ is more of a barrier to water movements than the ridge to the south. Comparison of front distributions with satellite altimetry data indicates that the MAR influence on deep ocean currents is also frequently exhibited in surface temperature. The improved spatial and temporal resolution of the front analysis has revealed consistent seasonality in the branching patterns. These results contribute to our understanding of the variability of the NAC, and the techniques for visualising oceanic fronts can be applied in other regions to reveal details of surface currents that cannot be resolved using satellite altimetry or in situ measurements.

The rupture process of recent tsunamigenic earthquakes by geophysical data inversion

Subduction zones are the favorite places to generate tsunamigenic earthquakes, where friction between oceanic and continental plates causes the occurrence of a strong seismicity. The topics and the methodologies discussed in this thesis are focussed to the understanding of the rupture process of the seismic sources of great earthquakes that generate tsunamis. The tsunamigenesis is controlled by several kinematical characteristic of the parent earthquake, as the focal mechanism, the depth of the rupture, the slip distribution along the fault area and by the mechanical properties of the source zone. Each of these factors plays a fundamental role in the tsunami generation. Therefore, inferring the source parameters of tsunamigenic earthquakes is crucial to understand the generation of the consequent tsunami and so to mitigate the risk along the coasts. The typical way to proceed when we want to gather information regarding the source process is to have recourse to the inversion of geophysical data that are available. Tsunami data, moreover, are useful to constrain the portion of the fault area that extends offshore, generally close to the trench that, on the contrary, other kinds of data are not able to constrain. In this thesis I have discussed the rupture process of some recent tsunamigenic events, as inferred by means of an inverse method. I have presented the 2003 Tokachi-Oki (Japan) earthquake (Mw 8.1). In this study the slip distribution on the fault has been inferred by inverting tsunami waveform, GPS, and bottom-pressure data. The joint inversion of tsunami and geodetic data has revealed a much better constrain for the slip distribution on the fault rather than the separate inversions of single datasets. Then we have studied the earthquake occurred on 2007 in southern Sumatra (Mw 8.4). By inverting several tsunami waveforms, both in the near and in the far field, we have determined the slip distribution and the mean rupture velocity along the causative fault. Since the largest patch of slip was concentrated on the deepest part of the fault, this is the likely reason for the small tsunami waves that followed the earthquake, pointing out how much the depth of the rupture plays a crucial role in controlling the tsunamigenesis. Finally, we have presented a new rupture model for the great 2004 Sumatra earthquake (Mw 9.2). We have performed the joint inversion of tsunami waveform, GPS and satellite altimetry data, to infer the slip distribution, the slip direction, and the rupture velocity on the fault. Furthermore, in this work we have presented a novel method to estimate, in a self-consistent way, the average rigidity of the source zone. The estimation of the source zone rigidity is important since it may play a significant role in the tsunami generation and, particularly for slow earthquakes, a low rigidity value is sometimes necessary to explain how a relatively low seismic moment earthquake may generate significant tsunamis; this latter point may be relevant for explaining the mechanics of the tsunami earthquakes, one of the open issues in present day seismology. The investigation of these tsunamigenic earthquakes has underlined the importance to use a joint inversion of different geophysical data to determine the rupture characteristics. The results shown here have important implications for the implementation of new tsunami warning systems – particularly in the near-field – the improvement of the current ones, and furthermore for the planning of the inundation maps for tsunami-hazard assessment along the coastal area.

Sea level trends in the Mediterranean from tide gauges and satellite altimetry

Sea level variation is one of the parameters directly related to climate change. Monitoring sea level rise is an important scientific issue since many populated areas of the world and megacities are located in low-lying regions. At present, sea level is measured by means of two techniques: the tide gauges and the satellite radar altimetry. Tide gauges measure sea-level relatively to a ground benchmark, hence, their measurements are directly affected by vertical ground motions. Satellite radar altimetry measures sea-level relative to a geocentric reference and are not affected by vertical land motions. In this study, the linear relative sea level trends of 35 tide gauge stations distributed across the Mediterranean Sea have been computed over the period 1993-2014. In order to extract the real sea-level variation, the vertical land motion has been estimated using the observations of available GPS stations and removed from the tide gauges records. These GPS-corrected trends have then been compared with satellite altimetry measurements over the same time interval (AVISO data set). A further comparison has been performed, over the period 1993-2013, using the CCI satellite altimetry data set which has been generated using an updated modeling. The absolute sea level trends obtained from satellite altimetry and GPS-corrected tide gauge data are mostly consistent, meaning that GPS data have provided reliable corrections for most of the sites. The trend values range between +2.5 and +4 mm/yr almost everywhere in the Mediterranean area, the largest trends were found in the Northern Adriatic Sea and in the Aegean. These results are in agreement with estimates of the global mean sea level rise over the last two decades. Where GPS data were not available, information on the vertical land motion deduced from the differences between absolute and relative trends are in agreement with the results of other studies.

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.

Global and Regional Sea level budget components from GRACE and radar altimetry (2002-2014)

This dataset contains the result of a joint least squares inversion of GRACE and altimetry data. The results are evaluated in terms of sea level change for the global ocean as well as dedicated areas. In addition, some auxiliary data is provided to enable reproducibility of the results in Rietbroek et al. 2016, and a google Earth kmz file is provided which visualizes the trends derived from the inversion results.

Cross-calibrating ALES Envisat and CryoSat-2 Delay-Doppler: A coastal altimetry study in the Indonesian Seas

A regional cross-calibration between the first Delay Doppler altimetry dataset from Cryosat-2 and a retracked Envisat dataset is here presented, in order to test the benefits of the Delay-Doppler processing and to expand the Envisat time series in the coastal ocean. The Indonesian Seas are chosen for the calibration, since the availability of altimetry data in this region is particularly beneficial due to the lack of in-situ measurements and its importance for global ocean circulation. The Envisat data in the region are retracked with the Adaptive Leading Edge Subwaveform (ALES) Retracker, which has been previously validated and applied successfully to coastal sea level research. The study demonstrates that CryoSat-2 is able to decrease the 1-Hz noise of sea level estimations by 0.3 cm within 50 km of the coast, when compared to the ALES-reprocessed Envisat dataset. It also shows that Envisat can be confidently used for detailed oceanographic research after the orbit change of October 2010. Cross-calibration at the crossover points indicates that in the region of study a sea state bias correction equal to 5% of the significant wave height is an acceptable approximation for Delay-Doppler altimetry. The analysis of the joint sea level time series reveals the geographic extent of the semiannual signal caused by Kelvin waves during the monsoon transitions, the larger amplitudes of the annual signal due to the Java Coastal Current and the impact of the strong La Nina event of 2010 on rising sea level trends.

Cross-calibrating ALES Envisat and CryoSat-2 Delay-Doppler: A coastal altimetry study in the Indonesian Seas

A regional cross-calibration between the first Delay Doppler altimetry dataset from Cryosat-2 and a retracked Envisat dataset is here presented, in order to test the benefits of the Delay-Doppler processing and to expand the Envisat time series in the coastal ocean. The Indonesian Seas are chosen for the calibration, since the availability of altimetry data in this region is particularly beneficial due to the lack of in-situ measurements and its importance for global ocean circulation. The Envisat data in the region are retracked with the Adaptive Leading Edge Subwaveform (ALES) Retracker, which has been previously validated and applied successfully to coastal sea level research. The study demonstrates that CryoSat-2 is able to decrease the 1-Hz noise of sea level estimations by 0.3 cm within 50 km of the coast, when compared to the ALES-reprocessed Envisat dataset. It also shows that Envisat can be confidently used for detailed oceanographic research after the orbit change of October 2010. Cross-calibration at the crossover points indicates that in the region of study a sea state bias correction equal to 5% of the significant wave height is an acceptable approximation for Delay-Doppler altimetry. The analysis of the joint sea level time series reveals the geographic extent of the semiannual signal caused by Kelvin waves during the monsoon transitions, the larger amplitudes of the annual signal due to the Java Coastal Current and the impact of the strong La Nina event of 2010 on rising sea level trends.

Bifurcation of Pacific North Equatorial Current at the surface

The grid altimetry data between 1993 and 2006 near the Philippines were analyzed by the method of Empirical Orthogonal Function (EOF) to study the variation of bifurcation of the North Equatorial Current at the surface of the Pacific. The relatively short-term signals with periods of about 6 months, 4 months, 3 months and 2 months are found besides seasonal and interannual variations mentioned in previous studies. Local wind stress curl plays an important role in controlling variation of bifurcation latitude except in the interannual timescale. The bifurcation latitude is about 13.3A degrees N in annual mean state and it lies at the northernmost position (14.0A degrees N) in January, at the southernmost position (12.5A degrees N) in July. The amplitude of variation of bifurcation latitude in a year is 1.5A degrees, which can mainly be explained as the contributions of the signals with periods of about 1 year (1.2A degrees) and 0.5 year (0.3A degrees).

Eddy generation and evolution in the North Pacific Subtropical Countercurrent (NPSC) zone

The seasonal generation and evolution of eddies in the region of the North Pacific Subtropical Countercurrent remain poorly understood due to the scarcity of available data. We used TOPEX/POSEIDON altimetry data from 1992 to 2007 to study the eddy field in this zone. We found that velocity shear between this region and the neighboring North Equatorial Current contributes greatly to the eddy generation. Furthermore, the eddy kinetic energy level (EKE) shows an annual cycle, maximum in April/May and minimum in December/January. Analyses of the temporal and spatial distributions of the eddy field revealed clearly that the velocity shear closely related to baroclinic instability processes. The eddy field seems to be more zonal than meridional, and the energy containing length scale shows a surprising lag of 2-3 months in comparison with the 1-D and 2-D EKE level. A similar phenomenon is observed in individual eddies in this zone. The results show that in this eddy field band, the velocity shear may drive the EKE level change so that the eddy field takes another 2-3 months to grow and interact to reach a relatively stable state. This explains the seasonal evolution of identifiable eddies.

Remote sensing of intertidal morphological change in Morecambe Bay, U.K., between 1991 and 2007

Tidal Flats are important examples of extensive areas of natural environment that remain relatively unaffected by man. Monitoring of tidal flats is required for a variety of purposes. Remote sensing has become an established technique for the measurement of topography over tidal flats. A further requirement is to measure topographic changes in order to measure sediment budgets. To date there have been few attempts to make quantitative estimates of morphological change over tidal flat areas. This paper illustrates the use of remote sensing to measure quantitative and qualitative changes in the tidal flats of Morecambe Bay during the relatively long period 1991–2007. An understanding of the patterns of sediment transport within the Bay is of considerable interest for coastal management and defence purposes. Tidal asymmetry is considered to be the dominant cause of morphological change in the Bay, with the higher currents associated with the flood tide being the main agency moulding the channel system. Quantitative changes were measured by comparing a Digital Elevation Model (DEM) of the intertidal zone formed using the waterline technique applied to satellite Synthetic Aperture Radar (SAR) images from 1991–1994, to a second DEM constructed from airborne laser altimetry data acquired in 2005. Qualitative changes were studied using additional SAR images acquired since 2003. A significant movement of sediment from below Mean Sea Level (MSL) to above MSL was detected by comparing the two Digital Elevation Models, though the proportion of this change that could be ascribed to seasonal effects was not clear. Between 1991 and 2004 there was a migration of the Ulverston channel of the river Leven north-east by about 5 km, followed by the development of a straighter channel to the west, leaving the previous channel decoupled from the river. This is thought to be due to independent tidal and fluvial forcing mechanisms acting on the channel. The results demonstrate the effectiveness of remote sensing for measurement of long-term morphological change in tidal flat areas. An alternative use of waterlines as partial bathymetry for assimilation into a morphodynamic model of the coastal zone is also discussed.