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.     

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.     

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.     

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.     

基于QZSS L1信号的岸基GNSS-R码延迟海面测高

为了探究GNSS L1波段中信号不同调制方式的GNSS-R码相位延迟海面的测高性能，开展了双天线岸基GNSS-R海面测高实验，收集了采样率为40 MHz的原始中频数据。利用自主开发的GNSS-R测高软件接收机和反演软件对实验数据进行了处理分析，同时获得了基于QZSS L1C/A码和L1C码的GNSS-R海面测高结果，并分别与岸基同步观测的雷达高度计的测高值进行对比，以对反演结果进行精度评定。实验结果表明:基于QZSS L1C/A码和L1C码的海面测高精度最优可达0.63 m和0.4 m，L1C码的延迟GNSS-R测高精度明显高于L1C/A码。此外，GNSS-R测高精度会随着卫星高度角的增加而有所增加。      

Detection of Envisat RA2/ICE-1 retracked radar altimetry bias over the Amazon basin rivers using GPS

Altimetry is now routinely used to monitor stage variations over rivers, including in the Amazon basin. It is desirable for hydrologic studies to be able to combine altimetry from different satellite missions with other hydrogeodesy datasets such as leveled gauges and watershed topography. One requirement is to accurately determine altimetry bias, which could be different for river studies from the altimetry calibrated for deep ocean or lake applications. In this study, we estimate the bias in the Envisat ranges derived from the ICE-1 waveform retracking, which are nowadays widely used in hydrologic applications. As a reference, we use an extensive dataset of altitudes of gauge zeros measured by GPS collocated at the gauges. The thirty-nine gauges are spread along the major tributaries of the Amazon basin. The methodology consists in jointly modeling the vertical bias and spatial and temporal slope variations between altimetry series located upstream and downstream of each gauge. The resulting bias of the Envisat ICE-1 retracked altimetry over rivers is 1.044 ± 0.212 m, revealing a significant departure from other Envisat calibrations or from the Jason-2 ICE-1 calibration.     

Validity assessment of SAMOSA retracking for fully-focused SAR altimeter waveforms

We demonstrate in this work how we can take advantage of known unfocused SAR (UF-SAR) retracking methods (e.g. the physical SAMOSA model) for retracking of fully-focused SAR (FF-SAR) waveforms. Our insights are an important step towards consistent observations of sea surface height, significant wave height and backscatter coefficient (wind speed) with both UF-SAR and FF-SAR. This is of particular interest for SAR altimetry in the coastal zone, since coastal clutter may be filtered out more efficiently in the high-resolution FF-SAR waveform data, which has the potential to improve data quality. We implemented a multi-mission FF-SAR altimetry processor for Sentinel-3 (S3) and Sentinel-6 Michael Freilich (S6), using a back-projection algorithm, and analysed ocean waveform statistics compared to multilooked UF-SAR. We find for Sentinel-3 that the averaged power waveforms of UF-SAR and FF-SAR over ocean are virtually identical, while for Sentinel-6 the FF-SAR power waveforms better resemble the UF-SAR zero-Doppler beam. We can explain and model the similarities and differences in the data via theoretical considerations of the waveform integrals. These findings suggest to use the existing UF-SAR SAMOSA model for retracking S3 FF-SAR waveforms but the SAMOSA zero-Doppler beam model for S6 FF-SAR waveforms, instead. Testing the outlined approach over short track segments, we obtain range biases between UF-SAR and FF-SAR lower than 2 mm and significant wave height biases lower than 5 cm.     

基于DDMR辅助的GNSS-R载波相位差测高方法

提出了一种基于时延-多普勒映射接收机（DDMR,Delay-Doppler Map Receiver）辅助的载波相位差提取方法,给出了系统结构及信号处理方法.该方法将DDMR中所观测到的码相位差作为直射信号与反射信号的码相位延时量,将完成跟踪的直射信号扩频码进行对应的延时用于完成对反射信号的开环码跟踪.该方法省去了码相位延时搜索的过程,且可以准确地对反射信号扩频码进行同步.为了验证系统的可行性及实际性能,进行了针对水面高度测量的岸基试验并给出了试验结果.岸基试验证明采用该方法的GNSS-R（Global Navigation Satellite System Reflection）接收机可以稳定地对反射信号进行跟踪并提取直射与反射信号的载波相位差,测高精度约为2.5 cm,经过0.5 s的数据平均后精度可达0.6 cm.      

Evaluating the feasibility of GPS measurements of SSH on board a ship along the Portuguese West Coast

Aiming to validate coastal altimetry data, an exploratory experiment was carried out, using a survey ship with onboard GPS and motion compensation systems. The ship navigation plan followed a track as coincident as possible with the passage of two altimetry satellites, Envisat and Jason-1, along the Portuguese coast, in June 2007.     

Properties of transit trajectory in the restricted three and four-body problem

In the restricted three-body problem if the Jacobi constant is just below the value corresponding to Lagrangian point only a little neck exists around the equilibrium point and capture trajectories are indicated as low-energy. Capture properties depend on the dynamics around these critical points and qualitative results can be obtained using linearized systems. In this paper, to study transit trajectory properties in the restricted three and four-body problem, the Earth–Moon–Sun–Satellite system is considered as example and studied using different models. In the circular restricted three-body problem (Earth–Moon–Satellite), transit, non transit and asymptotic trajectories, are easily identified by using the principal reference frame. Dynamics around Lagrangian point are then studied introducing the Moon eccentricity into the elliptical restricted three-body model. A preferential region for transit orbit is individuated and studied as a function of eigenvalue properties. To introduce the Sun effect, the bi-circular four-body model is considered and dynamics around Lagrangian points studied as a function of angular distance between Earth–Sun and Earth–Moon line. Finally, results obtained in the elliptical three-body model and bi-circular four-body model, are compared with numerical simulations using real Sun–Moon–Earth ephemeris.     

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.     

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.     

基于半无码的P(Y)码自相关GNSS-R海面测高方法

针对传统GPS C/A码延迟测高方法测高精度低和互相关测高方法时延波形信噪比低的问题,研究了基于半无码的P（Y）码自相关海面测高方法。首先,通过理论测高精度模型,对比分析了C/A码和P（Y）码的测高精度;然后,描述了L1、L2双频反射信号基于半无码的P（Y）码自相关处理架构,并论述了测高模型和L1、L2双频电离层时延误差修正方法;最后,利用模拟的L1、L2反射信号对比分析了基于半无码的P（Y）码自相关、C/A码自相关以及信号互相关方法的海面测高精度。结果表明:相比于C/A码自相关和信号互相关测高精度,基于半无码的P（Y）码自相关的测高精度分别提高了3.97倍和1.47倍。      

DORIS system: The new age

The boarding of the first DGXX DORIS instrument on Jason-2 mission gives us the opportunity to present the improvements that have been implemented on the DORIS system. The goal of this paper is to present information about the new capacities of the DORIS system and to give the current status of its components. An overview of the DORIS system, the International DORIS Service and the Jason-2 satellite mission are first presented. Then the new characteristics of the on-board instrument are detailed. The capacity to track up to seven ground beacons simultaneously dramatically increases the number of measurements performed: a factor of three increase over Jason-1 is observed at the altitude of 1330 km. It also increases the diversity of directions of observation and allows low elevation measurements from 0°. The new phase measurements capability allows now phase processing. The instability of the Jason-1 USOs (Ultra-Stable Oven-controlled quartz oscillator) while crossing the South Atlantic Anomaly has been solved by decreasing the sensitivity to radiation by a factor of 10. New features of the on-board software enhance the coastal and inland water altimetry and increase the robustness of the data. The new software also improves the real time orbit accuracy for operational altimetry. The improvements introduced concurrently on the ground segment have also significantly enhanced capability. The new RINEX exchange formats provide simultaneous phase and pseudo-range measurements. The maintenance of the DORIS Beacons Network and the work done by the DORIS Signal Integrity monitoring team lead to an increased availability of the Network from 75% to 90% and so to a more homogenous orbit coverage.     

航天测绘发展现状与展望

航天测绘是全球基础测绘的主要手段,包括卫星地形测量、卫星重力测量、卫星磁力测量、卫星测高、卫星海洋测量等。航天测绘工程的发展经历了“返回型测绘卫星”工程、“传输型测绘卫星”工程和部分“集成型测绘卫星”工程。首先梳理了国际卫星测量的发展历史和现状,阐述了我国航天测绘各个阶段所发挥的作用,并分析了存在的问题,最后对航天测绘的未来发展提出了思路。     

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.     

最优气动力辅助空间拦截和交会

航天器从高轨道向位于低地球圆轨道的靶目标实施拦截或进行交会 ,并用间接求解的直接伴随方法 ,即用D型拉格朗日方法求解其优化问题。拦截条件是指在最终时刻两者的位置相同 ,而交会条件是指不仅位置相同且最终的速度矢量亦相等。文章讨论了最小时间与最小脱靶量 ,或最小时间与最小燃料消耗的复合性能指标 ,以及推力协同机动的优化解。推力协同下 ,对于特定情况 ,如在热流约束边界部分推力弧是奇异的 ,气动力控制是正常的。这种复杂性和奇异问题的求解在拦截和交会问题中都是值得慎重处理的。     

Shaping low-thrust trajectories with thrust-handling feature

Shape-based methods are becoming popular in low-thrust trajectory optimization due to their fast computation speeds. In existing shape-based methods constraints are treated at the acceleration level but not at the thrust level. These two constraint types are not equivalent since spacecraft mass decreases over time as fuel is expended. This paper develops a shape-based method based on a Fourier series approximation that is capable of representing trajectories defined in spherical coordinates and that enforces thrust constraints. An objective function can be incorporated to minimize overall mission cost, i.e., achieve minimum $\mathrm{\Delta }V$. A representative mission from Earth to Mars is studied. The proposed Fourier series technique is demonstrated capable of generating feasible and near-optimal trajectories. These attributes can facilitate future low-thrust mission designs where different trajectory alternatives must be rapidly constructed and evaluated.