Lagrangian analysis of satellite-derived currents: Application to the North Western Mediterranean coastal dynamics

Optimal interpolation methods for improving the reconstruction of coastal dynamics from along-track satellite altimetry measurements have been recently developed over the North Western Mediterranean Sea. Maps of satellite-derived geostrophic current anomalies are generated using these methods, and added to different mean circulation fields in order to obtained absolute geostrophic currents. The resulting fields are then compared to standard AVISO products, and their accuracies are assessed with Lagrangian diagnostics. The trajectories of virtual particle clusters are simulated with a Lagrangian code either with new current fields or with the AVISO ones. The simulated trajectories are then compared to 16 in situ drifter trajectories to evaluate the performance of the different velocity fields. The comparisons show that the new current fields lead to better results than the AVISO one, especially over the shallow continental shelf of the Gulf of Lion. However, despite the use of innovative strategies, some altimetry limitations still persist in the coastal domain, where small scale processes remain sub-sampled by conventional altimetry coverage but will benefit from technological development in the near future. Some of the limitations of the Lagrangian diagnostics presently used are also analyzed, but dedicated studies will be required for future further investigations.     

Estimates of vertical land motion along the southwestern coasts of Turkey from coastal altimetry and tide gauge data

The differences between coastal altimetry and sea level time series of tide gauges in between March 1993 and December 2009 are used to estimate the rates of vertical land motion at three tide gauge locations along the southwestern coasts of Turkey. The CTOH/LEGOS along-track coastal altimetry retrieves altimetric sea level anomalies closer to the coast than the standard along-track altimetry products. However, the use of altimetry very close to the coast is not found to improve the results. On the contrary, the gridded and interpolated AVISO merged product exhibits the best agreement with tide gauge data as it provides the smoothest variability both in space and time compared with along track altimetry data. The Antalya gauge to the south (in the Mediterranean Sea) and the Mentes/Izmir gauge to the west (in the Aegean Sea) both show subsidence while the Bodrum tide gauge to the south (in the Aegean Sea) shows no significant vertical land motion. The results are compared and assessed with three independent geophysical vertical land motion estimates like from GPS. The GIA effect in the region is negligible. The VLM estimates from altimetry and tide gauge data are in good agreement both with GPS derived vertical velocity estimates and those inferred from geological and archaeological investigations.     

Preliminary gravity recovery from CryoSat-2 data in the Baffin Bay

A number of geophysical phenomenons in the open ocean are still unresolved by conventional altimetry, but could be resolved through the potential improvements offered by Synthetic Aperture Radar (SAR), also called Delay-Doppler, altimetry. The SAR altimeter offers the following benefits with respect to conventional satellite altimetry: factor of 20 improvements in the along-track resolution, the along-track footprint length which does not vary with wave height (sea state), and improved precision in sea surface height measurements or sea surface slope measurements.     

Exploring the potential of Sentinel-3 delay Doppler altimetry for enhanced detection of coastal currents along the Northwest Atlantic shelf

In this study, we evaluate Sentinel-3A satellite synthetic aperture radar (SAR) altimeter observations along the Northwest Atlantic coast, spanning the Nova Scotian Shelf, Gulf of Maine, and Mid-Atlantic Bight. Comparisons are made of altimeter sea surface height (SSH) measurements from three different altimeter data processing approaches: fully-focused synthetic aperture radar (FFSAR), un-focused SAR (UFSAR), and conventional low-resolution mode (LRM). Results show that fully-focused SAR data always outperform LRM data and are comparable or slightly better than the nominal un-focused SAR product. SSH measurement noise in both SAR-mode datasets is significantly reduced compared to LRM. FFSAR SSH 20-Hz noise levels, derived from 80-Hz FFSAR data, are lower than 20-Hz UFSAR SSH with 25% noise reduction offshore of 5 km, and 55–70% within 5 km of the coast. The offshore noise improvement is most likely due to the higher native along-track data posting rate (80 Hz for FFSAR, and 20 Hz for UFSAR), while the large coastal improvement indicates an apparent FFSAR data processing advantage approaching the coastlines. FFSAR-derived geostrophic ocean current estimates exhibit the lowest bias and noise when compared to in situ buoy-measured currents. Assessment at short spatial scales of 5–20 km reveals that Sentinel-3A SAR data provide sharper and more realistic measurement of small-scale sea surface slopes associated with expected nearshore coastal currents and small-scale gyre features that are much less well resolved in conventional altimetric LRM data.     

The challenges in long-term altimetry calibration for addressing the problem of global sea level change

Long-term change of the global sea level resulting from climate change has become an issue of great societal interest. The advent of the technology of satellite altimetry has modernized the study of sea level on both global and regional scales. In combination with in situ observations of the ocean density and space observations of Earth’s gravity variations, satellite altimetry has become an essential component of a global observing system for monitoring and understanding sea level change. The challenge of making sea level measurements with sufficient accuracy to discern long-term trends and allow the patterns of natural variability to be distinguished from those linked to anthropogenic forcing rests largely on the long-term efforts of altimeter calibration and validation. The issues of long-term calibration for the various components of the altimeter measurement system are reviewed in the paper. The topics include radar altimetry, the effects of tropospheric water vapor, orbit determination, gravity field, tide gauges, and the terrestrial reference frame. The necessity for maintaining a complete calibration effort and the challenges of sustaining it into the future are discussed.     

The X-TRACK/ALES multi-mission processing system: New advances in altimetry towards the coast

In the context of the ESA Climate Change Initiative project, a new coastal sea level altimetry product has been developed in order to support advances in coastal sea level variability studies. Measurements from Jason-1,2&3 missions have been retracked with the Adaptive Leading Edge Subwaveform (ALES) Retracker and then ingested in the X-TRACK software with the best possible set of altimetry corrections. These two coastal altimetry processing approaches, previously successfully validated and applied to coastal sea level research, are combined here for the first time in order to derive a 16-year-long (June 2002 to May 2018), high-resolution (20-Hz), along-track sea level dataset in six regions: Northeast Atlantic, Mediterranean Sea, West Africa, North Indian Ocean, Southeast Asia and Australia. The study demonstrates that this new coastal sea level product called X-TRACK/ALES is able to extend the spatial coverage of sea level altimetry data ~3.5 km in the land direction, when compared to the X-TRACK 1-Hz dataset. We also observe a large improvement in coastal sea level data availability from Jason-1 to Jason-3, with data at 3.6 km, 1.9 km and 0.9 km to the coast on average, for Jason-1, Jason-2 and Jason-3, respectively. When combining measurements from Jason-1 to Jason-3, we reach a distance of 1.2–4 km to the coast. When compared to tide gauge data, the accuracy of the new altimetry near-shore sea level estimations also improves. In terms of correlations with a large set of independent tide gauge observations selected in the six regions, we obtain an average value of 0.77. We also show that it is now possible to derive from the X-TRACK/ALES product an estimation of the ocean current variability up to 5 km to the coast. This new altimetry dataset, freely available, will provide a valuable contribution of altimetry in coastal marine research community.     

Analysis of coastal altimetry in the Mexican Caribbean

A reprocessing of sea-level anomalies (SLA) resulting from X-TRACK coastal altimetry was carried out for the ENVISAT (2002–2010) and TOPEX/POSEIDON-Jason (1992–2019) satellite missions in the coastal area of the Mexican Caribbean. This consisted of applying a tidal correction to coastal altimetry sea level observations. Harmonic analysis of five coastal tide gauge records was performed to estimate the most important tidal components of the area, resulting on M2, N2, O1, S2, K1, MF, and MM. The tidal signal was reconstructed with the seven tidal components using the TPXO9 model. The SLA signals corrected with the seven tidal components were validated with in situ data from coastal tide gauges. The validation showed that the TPXO9 tidal barotropic model (1/30° grid) used to reconstruct the tidal signal with the seven representative tidal components performed better than the FES2012 global model (1/16° grid) that uses 33 tidal components. The reprocessed SLAs showed clear seasonality with significant signals at 4, 6, and 12 months, with the annual signal being the dominant one. In the Mexican Caribbean coastal zone, oceanographic processes with different scales (from coastal to mesoscale) converge, showing their complexity in the different SLA signals observed. The aim of this work is to contribute to the analysis of coastal altimetry data and understanding the sea level variations in the Mexican Caribbean. This work is the first step in the implementation of methodologies that take advantage of coastal satellite altimetry in the Caribbean Sea.     

Remote sensing using GNSS signals: Current status and future directions

The refracted, reflected and scattered signals of global navigation satellite systems (GNSS) have been successfully used to remotely sense the Earth’s surface and atmosphere. It has demonstrated its potential to sense the atmosphere and ionosphere, ocean, land surfaces (including soil moisture) and the cryosphere. These new measurements, although in need of refinement and further validation in many cases, can be used to complement existing techniques and sensors, e.g., radiosonde, ionosonde, radar altimetry and synthetic aperture radar (SAR). This paper presents the current status and new developments of remote sensing using GNSS signals as well as its future directions and applications. Some notable emerging applications include monitoring sea ice, dangerous sea states, ocean eddy and storm surges. With the further improvement of the next generation multi-frequency GNSS systems and receivers and new space-based instruments utilizing GNSS reflections and refractions, new scientific applications of GNSS are expected in various environment remote sensing fields in the near future.     

Mapping mean lake surface from satellite altimetry and GPS kinematic surveys

Lake water height is a key variable in water cycle and climate change studies, which is achievable using satellite altimetry constellation. A method based on data processing of altimetry from several satellites has been developed to interpolate mean lake surface (MLS) over a set of 22 big lakes distributed on the Earth. It has been applied on nadir radar altimeters in Low Resolution Mode (LRM: Jason-3, Saral/AltiKa, CryoSat-2) in Synthetic Aperture Radar (SAR) mode (Sentinel-3A), and in SAR interferometric (SARin) mode (CryoSat-2), and on laser altimetry (ICESat). Validation of the method has been performed using a set of kinematic GPS height profiles from 18 field campaigns over the lake Issykkul, by comparison of altimetry’s height at crossover points for the other lakes and using the laser altimetry on ICESat-2 mission. The precision reached ranges from 3 to 7 cm RMS (Root Mean Square) depending on the lakes. Currently, lake water level inferred from satellite altimetry is provided with respect to an ellipsoid. Ellipsoidal heights are converted into orthométric heights using geoid models interpolated along the satellite tracks. These global geoid models were inferred from geodetic satellite missions coupled with absolute and regional anomaly gravity data sets spread over the Earth. However, the spatial resolution of the current geoid models does not allow capturing short wavelength undulations that may reach decimeters in mountaineering regions or for rift lakes (Baikal, Issykkul, Malawi, Tanganika). We interpolate in this work the geoid height anomalies with three recent geoid models, the EGM2008, XGM2016 and EIGEN-6C4d, and compare them with the Mean Surface of 22 lakes calculated using satellite altimetry. Assuming that MLS mimics the local undulations of the geoid, our study shows that over a large set of lakes (in East Africa, Andean mountain and Central Asia), short wavelength undulations of the geoid in poorly sampled areas can be derived using satellite altimetry. The models used in this study present very similar geographical patterns when compared to MLS. The precision of the models largely depends on the location of the lakes and is about 18 cm, in average over the Earth. MLS can serve as a validation dataset for any future geoid model. It will also be useful for validation of the future mission SWOT (Surface Water and Ocean Topography) which will measure and map water heights over the lakes with a high horizontal resolution of 250 by 250 m.     

An Integer Precise Point Positioning technique for sea surface observations using a GPS buoy

GPS data dedicated to sea surface observation are usually processed using differential techniques. Unfortunately, the precision of resulting kinematic positions is baseline-length dependent. So, high precision sea surface observations using differential GPS techniques are limited to coasts, lakes, and rivers. Recent improvements in GPS satellite products (orbits, clocks, and phase biases) make phase ambiguity fixing at the zero difference level achievable and opens up the observation of the sea surface without geographical constraints. This paper recalls the concept of the Integer Precise Point Positioning technique and discusses the precision of GPS buoy positioning. A sequential version of the GINS software has been implemented to achieve single epoch GPS positioning. We used 1 Hz data from a two week GPS campaign conducted in the Kerguelen Islands. A GPS buoy has been moored close to a radar gauge and 90 m away from a permanent GPS station. This infrastructure offers the opportunity to compare both kinematic Integer Precise Point Positioning and classical differential GPS positioning techniques to in situ radar gauge data. We found that Precise Point Positioning results are not significantly biased with respect to radar gauge data and that horizontal time series are consistent with differential processing at the sub-centimetre precision level. Nevertheless, standard deviations of height time series with respect to radar gauge data are typically [4–5] cm. The dominant driver for noise at this level is attributed to errors in tropospheric estimates which propagate into position solutions.     

Semi-automatic determination of the Azores Current axis using satellite altimetry: Application to the study of the current variability during 1995–2006

Satellite altimetry has been widely used to study the variability of the ocean currents such as the Azores Current (AzC) in the North Atlantic. Most analyses are performed over the region that encloses the current, thus being somehow affected by other oceanographic signals, e.g., eddies. In this study, a new approach for extracting the axis of a zonal current solely based on satellite altimetry is presented. This is a semi-automatic procedure that searches for the maximum values of the gradient of absolute dynamic topography (ADT), using the geostrophic velocity as auxiliary information. The advantage of this approach is to allow the analyses to be performed over a buffer centered on the current axis instead of using a wider region. It is here applied to the AzC for the period June 1995–October 2006.     

基于多种观测手段的东亚低纬电离层不规则体事件分析

利用海南台站（19.5°N,109.1°E,dip:13.6°N）和磁赤道区的多种地基和天基观测数据,对2011年11月20日观测到的电离层不规则体事件进行了分析.海南台站VHF雷达、电离层闪烁和数字测高仪的综合观测结果表明,当天日落附近发生了强的电离层不规则体事件,主要表现为雷达羽和强闪烁的形态.结合磁赤道区GPS和C/NOFS卫星观测结果进行分析可知,海南台站日落附近出现的雷达羽和强闪烁与南海磁赤道区产生的主等离子体泡存在明显联系.      

A 2-D comparison of ionospheric convection derived from SuperDARN and DMSP measurements

Ion drift vectors measured by the DMSP satellites are compared with plasma convection vectors obtained by the SuperDARN HF radars through the standard Map Potential algorithm of Ruohoniemi and Baker [Ruohoniemi, M., Baker, K.B. Large-scale imaging of high-latitude convection with super dual auroral radar network HF radar observations. J. Geophys. Res. 103, 20797–20811, 1998]. Despite significant data spread, the agreement can be qualified as reasonable for a data set comprising of 149 satellite passes over the Northern Hemisphere at high latitudes. The slope of the best-fit line relating SuperDARN and DMSP velocity magnitudes is of the order of 0.3 with a tendency for the SuperDARN velocities to be smaller. The agreement between the azimuths of the ion drift and convection is better with the slope of the best-fit line being close to 1. It is shown that consistency between the radar and satellite measurements is much better if the SuperDARN line-of-sight velocities are compared with the DMSP cross-track ion drifts for events showing slow spatial and temporal variations of the convection. If areas of strong convection changes are included into comparison, the degree of agreement deteriorates drastically. This result implies that differences in the spatial and temporal resolutions of DMSP and SuperDARN measurements are crucial factors contributing to the observed discrepancies. In addition, some differences are introduced when the SuperDARN line-of-sight velocities are filtered and reprocessed into vectors with the application of a background convection model.     

A comparison of seasonal variations of sea level in the southern Baltic Sea from altimetry and tide gauge data

In this paper, seasonal sea level variations have been determined at five locations in the Baltic Sea from satellite altimetry for the period 1993–2015. The results were compared to tide gauge water level data. Annual and semi-annual amplitudes have been investigated for both sea level anomalies and tide gauge water level. It was found that the two independent observations of sea level variations along the Polish coast are in good agreement both in terms of their annual and semi-annual amplitudes and their annual and semi-annual phases. The annual cycles in the sea level variations measured by altimetry and tide gauge reach maximum values at approximately the same month (November/December).Moreover, this article shows the differences between the annual and semi-annual amplitudes and phases in the sea level anomalies and water level data within the same time frame. The difference in the annual amplitudes between the satellite altimetry and the tide gauge results is between 0.33?cm and 1.53?cm. The maximum differences in the annual cycle of the sea level changes were found at the Swinoujscie station. The correlations between the original series and the calculated curves were determined, and the relationship between the amplitudes and the phases were investigated. The correlation between the annual variations observed from the two independent observation techniques is 0.92.To analyse the dynamics of the change in sea level, the linear trend was estimated from the satellite altimetry and tide gauge time series both in the original time series of the data and in the time series in which seasonal variations were removed. In addition, we calculated the estimated errors of regression and how many years’ worth of data are needed to obtain an accuracy of 0.1?mm per year. The estimated errors of regression showed that to get an accuracy of 0.1?mm per year, we need 100?years of data.     

全球海表流场多尺度结构观测卫星计划

全球海表流场多尺度结构观测卫星计划（Ocean Surface Current multiscale Observation Mission, OSCOM）首次提出海表流场、海面风场和海浪谱（简称 “流–风–浪”）一体化探测的多普勒散射计（Doppler Scatterometer, DOPS）测量原理和系统体制。OSCOM采用Ka-Ku双频多波束圆锥扫描体制的真实孔径雷达,将实现超过1000 km观测刈幅、公里级分辨率的“流–风–浪”一体化卫星直接观测。OSCOM将突破海洋亚中尺度非平衡态动力学、海洋多尺度相互作用、海气耦合的研究瓶颈,支撑实现海洋系统科学、气候变化等理论研究的重大突破。未来,应用OSCOM海表流速观测的模式改进,将奠定海洋非平衡态过程数值模拟、同化和预报的动力学基础,实现海洋和海气耦合模式的重大改进。通过与多源数据融合,OSCOM海流观测的应用将为海洋生物地球化学循环、碳收支研究和国家重大任务提供支撑。OSCOM科学卫星的实施对于我国地球系统科学和卫星对地观测重大应用的突破有至关重要的意义,有望带动我国应用卫星的发展从追赶、并行走向领跑。      

Synergetic use of SAR and Thermal Infrared data to study the physical properties of the lunar surface

The surface layer of the Moon preserves vital evidences of lunar impact and cratering processes due to the absence of any Aeolian and fluvial erosion processes acting on it. By examining these evidences, which are recorded throughout the evolutionary history of the Moon, several basic aspects of lunar science can be understood, and this has direct relevance to the surfaces of other airless bodies within the solar system. In this study, rock abundance data obtained from Thermal Infrared (TIR) observations and radar Circular Polarization Ratio (CPR) data sets obtained from polarimetric SAR observations were correlated at some sample sites on the lunar surface. Preliminary results yielded qualitative and quantitative estimates for surface rock abundances. Except at distal ejecta deposits of young, bright craters a general correlation was observed between the two datasets. Mixed results were observed from the impact melt flows where the situation is complex due to the possible subsurface-volume and volume-subsurface interactions of the radar waves. But the flow features were clearly separated from the interior and ejecta regions of their parent craters in terms of CPR and rock abundances. The extent and distributions of pyroclastic deposits and dark haloed regions could not be distinctly identified at the resolution of datasets utilized. Near Gerasimovich D crater, the Diviner Radiometer has provided the first TIR observations of a newly discovered impact melt flow which was not visible in the optical imagery. This facilitated the first ever quantitative comparisons of the radar CPR and rock abundance values near such a region. Also, significant differences in spatial patterns between the radar and rock concentration data sets were observed, owing to the differences in the sensitivity of the two observations.     

Atmospheric correction of satellite altimetry observations and sea-level variability in the NE Atlantic

Satellite altimetry provides continuous and spatially regular measurements of the height of the sea surface. Sea level responds to density changes of the water, to mass changes, due to addition or reduction of water mass, and to changes in the atmosphere above it. The present study examines the influence of atmospheric effects on sea-level variability in the North-East Atlantic. The association between the height of the sea surface and the North Atlantic Oscillation (NAO) is investigated by considering different sets of altimetry measurements for which the atmospheric effects have been handled differently. Altimetry data not corrected for atmospheric effects are strongly anti-correlated with the state of the NAO, reflecting the hydrostatic response of sea-level to the NAO pressure dipole. The application of an atmospheric correction to satellite altimetry observations in the NE Atlantic decreases variability of the height time series by more than 70% and reduces the amplitude of the seasonal cycle by ∼5 cm. Altimetry data for which atmospheric effects are removed via an inverse barometer correction show a non-negligible correlation with the NAO index at some locations suggesting further indirect non-hydrostatic influences of the state of the NAO on sea level variability.     

DORIS and GPS monitoring of the Gavdos calibration site in Crete

Due to its specific geographical location as well as its geodetic equipment (DORIS, GNSS, microwave transponder and tide gauges), the Gavdos station in Crete, Greece is one of the very few sites around the world used for satellite altimetry calibration. To investigate the quality of the Gavdos geodetic coordinates and velocities, we analyzed and compared here DORIS and GPS-derived results obtained during several years of observations. The DORIS solution is the latest ignwd11 solution at IGN, expressed in ITRF2008, while the GPS solution was obtained using the GAMIT software package. Current results show that 1–2 mm/yr agreement can be obtained for 3-D velocity, showing a good agreement with current geophysical models. In particular, the agreement obtained for the vertical velocity is around 0.3–0.4 mm/yr, depending on the terrestrial reference frame. As a by-product of these geodetic GPS and DORIS results, Zenith Tropospheric Delays (ZTDs) estimations were also compared in 2010 between these two techniques, and compared to ECMWF values, showing a 6.6 mm agreement in dispersion without any significant difference between GPS and DORIS (with a 97.6% correlation), but with a 13–14 mm agreement in dispersion when comparing to ECMWF model (with only about 90% correlation for both techniques). These tropospheric delay estimations could also provide an external calibration of the tropospheric correction used for the geophysical data of satellite altimetry missions.     

Integrating Non-Tidal Sea Level data from altimetry and tide gauges for coastal sea level prediction

The main objective of this paper is to integrate Non-Tidal Sea Level (NSL) from the joint TOPEX, Jason-1 and Jason-2 satellite altimetry with tide gauge data at the west and north coast of the United Kingdom for coastal sea level prediction. The temporal correlation coefficient between altimetric NSLs and tide gauge data reaches a maximum higher than 90% for each gauge. The results show that the multivariate regression approach can efficiently integrate the two types of data in the coastal waters of the area. The Multivariate Regression Model is established by integrating the along-track NSL from the joint TOPEX/Jason-1/Jason-2 altimeters with that from eleven tide gauges. The model results give a maximum hindcast skill of 0.95, which means maximum 95% of NSL variance can be explained by the model. The minimum Root Mean Square Error (RMSe) between altimetric observations and model predictions is 4.99 cm in the area. The validation of the model using Envisat satellite altimetric data gives a maximum temporal correlation coefficient of 0.96 and a minimum RMSe of 4.39 cm between altimetric observations and model predictions, respectively. The model is furthermore used to predict high frequency NSL variation (i.e., every 15 min) during a storm surge event at an independent tide gauge station at the Northeast of the UK (Aberdeen).     

Kuafu and the studies of CME initiation

We first briefly review the current trend in the studies of coronal mass ejections (CMEs), then summarize some recent efforts in understanding the CME initiation. Emphasis has been put on the studies of Earth-directed CMEs whose associated surface activity and large scale magnetic source have been well identified. The data analysis by combining the MDI full disc magnetograms, vector magnetograms of active regions, EUV waves and dimmings, non-thermal radio sources, and the SOHO LASCO observations has shed new light in understanding the CME magnetism. However, the current studies seem to invoke new observations in a few aspects: (1) The observations which enable us to trace CMEs from the earliest associated surface activity to its initial acceleration and key development in the low corona in the height of 1–3 R; (2) The imaging spectroscopic observations which can be used to diagnose the early plasma outflow and the line-of-sight velocity in understanding the kinematics of CMEs; (3) The accurate timing from primary magnetic energy release, manifested by chromospheric activity, non-thermal radio bursts, and EUV, X-ray and γ-ray emissions, to the CME initiation, early acceleration and propagation, and the consequences in the interplanetary space and magnetosphere. The Kuafu Mission will meet the basic requirement for the new observations in CME initiation studies and serve as a monitor of space weather of the Sun–Earth system.