High-performance altimeter Doppler processing for measuring sea level height under varying sea state conditions

This paper describes an innovative method for processing nadir altimeter data acquired in Synthetic Aperture Radar (SAR) mode, enhancing the system performances over open ocean. Similarly to the current SAR data processing scheme, the so-called LR-RMC (Low Resolution with Range Migration Correction) method, originally designed by Phalippou and Demeester (2011), includes Doppler beam forming, Doppler shift correction and range correction. In LR-RMC, however, an alternative and less complex averaging (stacking) operation is used so that all the Doppler beams produced in a radar cycle (4 bursts of 64 beams for the open-burst Sentinel-3-mode altimeter) are incoherently combined to form a multi-beam echo. In that manner, contrarily to the narrow-band SAR technique, the LR-RMC processing enlarges the effective footprint to average out the effects of surface waves and particularly those from small sub-mesoscale structures (<1 km) that are known to impact SAR-mode performances. On the other hand, the number of averaged beams is as high as in current SAR-mode processing, thus providing a noise reduction at least equally good. The LR-RMC method has the added benefit of reducing the incoherent integration time with respect to the SAR-mode processing (50 ms compared to 2.5 s) limiting possible surface movement effects. By processing one year of Sentinel-3A SRAL SAR-mode data using the LR-RMC method, it is shown that the swell impact on the SAR altimeter performances is totally removed and that an improvement of 10–50% is obtained in the measurement noise of the sea surface height and significant wave height with respect to SAR mode. Additionally, observational capabilities over the middle scales are enhanced potentially allowing the ocean mesoscale features to be retrieved and observations assimilated more usefully in ocean models.     

Impact of rain on wave retrieval from Sentinel-1 synthetic aperture radar images in tropical cyclones

The main objective of our work was to investigate the impact of rain on wave observations from C-band (~5.3 GHz) synthetic aperture radar (SAR) in tropical cyclones. In this study, 10 Sentinel-1 SAR images were available from the Satellite Hurricane Observation Campaign, which were taken under cyclonic conditions during the 2016 hurricane season. The third-generation wave model, known as Simulating WAves Nearshore (SWAN) (version 41.31), was used to simulate the wave fields corresponding to these Sentinel-1 SAR images. In addition, rainfall data from the Tropical Rainfall Measuring Mission satellite passing over the spatial coverage of the Sentinel-1 SAR images were collected. The simulated results were validated against significant wave heights (SWHs) from the Jason-2 altimeter and European Centre for Medium-Range Weather Forecasts data, revealing a root mean square error (RMSE) of ~0.5 m with a 0.25 scatter index. Winds retrieved from the VH-polarized Sentinel-1 SAR images using the Sentinel-1 Extra Wide-swath Mode Wind Speed Retrieval Model after Noise Removal were taken as prior information for wave retrieval. It was discovered that rain did indeed affect the SAR wave retrieval, as evidenced by the 3.21-m RMSE of SWHs between the SAR images and the SWAN model, which was obtained for the ~1000 match-ups with raindrops. The raindrops dampened the wave retrieval when the rain rate was < ~5 mm/hr; however, they enhanced wave retrieval for higher rain rates. It was also found that the portion of the rain-induced ring wave with a wave number > 0.05 rad/m (~125 m wavelength) was clearly observed in the SAR-derived wave spectra.     

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.     

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.     

Exploiting the Sentinel-3 tandem phase dataset and azimuth oversampling to better characterize the sensitivity of SAR altimeter sea surface height to long ocean waves

Being the very first SAR mode altimeter tandem phase, the Sentinel-3 A/B tandem phase has provided an unprecedented opportunity to better characterize the sensitivity of SAR altimetry retrievals to high-frequency processes, such as long ocean waves. In this paper, we show that for some sea-state conditions, that are still to be precisely characterized, long ocean waves are responsible for high-frequency (spatial and temporal) coherent Sea Level Anomaly (SLA) signals. It is found that the peak wavelength corresponds to the dominant swell wavelength. Furthermore, the short time lag between S3-A and S3-B acquisitions allows performing cross-spectral analyses that reveal phase shifts consistent with waves travelling according to the wave dispersion relation. It is also demonstrated that the classical 20 Hz sampling frequency is insufficient to properly sample most swell-induced SLA signals and that aliasing can generate errors over the entire frequency spectrum, including at long wavelengths. These results advocate for the use of azimuth oversampling (40 Hz or 80 Hz). Low-pass filtering should be applied prior to any down-sampling to 20 Hz, in order to prevent long-wavelength errors induced by spectral leakage.     

First quality assessment of the Cryosat-2 altimetric system over ocean

This paper presents the results of a Calval analysis performed for the Cryosat-2 mission over ocean. The data set used in this analysis consists of products generated by the Cryosat-2 Processor Prototype developed by CNES (Center National d’Etudes Spatiales). This data set has been analysed focusing on LRM (Low Resolution Mode) mode only. One major objective of this paper is to illustrate the potential of Cryosat-2 data over ocean, mainly for waves and Sea Level Anomaly applications. All the results indicate very good performances of the SIRAL (SAR/Interferometric Altimeter) altimeter over ocean. Crossover standard deviation is close to 6.5 cm over the analysed period (3 months) which is close to the Jason-2 and ENVISAT performance. All these results confirm that Cryosat’s altimeter can provide data almost as valuable as other flying altimetric missions, and that it has the potential to contribute to oceanography (e.g. multi-mission climate record, mesoscale monitoring in near real time) and to geodesy (e.g. mean sea surface, bathymetry).     

Spiral eddies in the Baltic,Black and Caspian seas as seen by satellite radar data

Despite spiral eddies were first seen on the sea surface more than 40 years ago, there is still a lot of uncertainty concerning these eddies. The present paper is aimed to provide the comprehensive results on the occurrence and statistics of small-scale eddies in the three inner seas (the Baltic, Black and Caspian) using satellite synthetic aperture radar (SAR) images. The dataset used includes over 2000 medium resolution Envisat ASAR and ERS-2 SAR images obtained in 2009–2010 in the different parts of the seas mentioned. As a result of the analysis performed ∼14,000 vortical structures were detected. 71% of them were visualized due to surfactant films (“black” eddies), while 29% due to wave/current interactions (“white” eddies). Practically all the eddies detected were cyclonically rotating. Their diameter was within 1–20 km. Characteristic size of the “black” eddies in all the basins was discovered to be less than that of the “white” eddies. Characteristic eddy size for the Baltic, Black and Caspian seas proved to be strictly proportional to the values of the baroclinic Rossby radius of deformation typical for these basins. The “black” eddies did not demonstrate a significant connection with the basin- and meso-scale surface circulation of the seas. Most of the “white” eddies detected were attributed to the zones with the most intense drift currents, i.e. those along the western boundaries and (in the Baltic Sea only) in the elongated parts of the basin.     

机载三通道SAR/GMTI性能分析及改进方法

在分析经典的三通道SAR(Synthetic Aperture Radar)地面动目标指示(GMTI,Ground Moving Target Indication)方法原理的基础上,提出一种较全面三通道SAR/GMTI性能评价方法.该方法同时考虑了载机存在速度误差、通道不匹配、杂波内部运动和系统噪声等主要误差因素,并分析了各种因素对系统性能的影响.经典的三通道SAR/GMTI方法存在的缺点是不能同时提高最大不模糊可检测速度、最小可检测速度和盲速指标的性能以满足系统设计的要求.给出了一种新的天线布置方案和动目标检测方法,改进的检测方法在不增加通道数的前提下,较好地解决了该问题,具有很好的实用性.最后通过计算机仿真验证了新方法的可行性.     

月地双基SAR成像的快速后向投影算法研究

为了提高对全球变化的大尺度观测能力,月基对地观测的概念应运而生。由于月基单站合成孔径雷达（Synthetic Aperture Radar,SAR）的双程传播会造成回波信噪比衰减,为此针对月地双基SAR工作模式展开研究。考虑到月、地、星复杂的相对运动为成像带来困难,根据月、地、星的几何关系建立平面直角回波模型。同时考虑到月球和卫星运行轨迹的弯曲性,分析模型误差,推导了精确的回波二维频谱表达式。在此基础上提出了一种基于快速后向投影的SAR成像算法,补偿了方位空变的模型误差,利用频域子孔径成像和频谱融合实现成像算法的快速化。通过分析可知月地双基SAR的理论回波信噪比与地月距离的平方成反比,与月基单站SAR相比有显著提升,最后通过X波段SAR仿真数据验证了文章算法处理月地双基SAR成像问题的效性。     

高精度星载SAR多普勒参数估算方法

提出了一种基于空间坐标转换,利用卫星位置、速度参数精确估算星载SAR(Synthetic Aperture Radar)全观测带多普勒参数的方法.利用卫星速度、位置,通过星载SAR空间几何模型和坐标转换关系,建立SAR图像中斜距同卫星下视角之间的四次方程,解出下视角并进一步计算出该斜距处的多普勒参数值.仿真结果表明,该方法在无卫星位置、速度误差情况下估算精度达到0.02Hz(多普勒中心频率)和2×10-4Hz/s(多普勒调频率);存在卫星位置测量误差(300m)以及速度测量误差(0.3m/s)的情况下,估算精度达到0.8Hz(多普勒中心频率)和0.07Hz/s(多普勒调频率).该方法适用于单星SAR以及分布式SAR高精度多普勒参数的估算.      

Semi-annual variations of chlorophyll concentration on the Eastern Tropical Pacific coast of Mexico

An analysis of the main remotly sensed oceanographic variables was conducted to characterize the dominant drivers that modulate the spatial-temporal variability of coastal phytoplankton biomass in the northern limit of the Eastern Tropical Pacific. The phytoplankton biomass was analyzed using monthly average satellite chlorophyll (Chlo) concentration data from MODIS sensors for the period 2003–2017. Oceanographic conditions were analyzed by considering (i) Sea Surface Temperature (SST) high-resolution data from the GHRSST project, (ii) wind stress calculated with data from the zonal and meridional components of the CCMP product, (iii) climatological precipitation, and (iv) climatological river flow. Chlo time series and spatial variability were analyzed using Hovmöller diagrams and Empirical Orthogonal Functions (EOF), respectively. A strong semi-annual signal in Chlo concentration along the coast was observed: the first peak occurs in winter-spring (5.0 mg·m⁻³) and a second one in summer-autumn (6.5 mg·m⁻³). A high year-round average of Chlo concentration (3.0–15.0 mg·m⁻³) was maintained by a 10 km wide along-shore coast, with maximum values spatially associated with river mouths, of which Santiago river registered the highest Chlo values (20.0 mg·m⁻³). Surface oceanographic conditions showed a marked annual cycle with warmer conditions (30–31 °C) from July to October and colder ones (23–24 °C) from December to April; during the cold period, seasonal wind stress stimulated coastal upwelling. The EOF showed that 70% of Chlo concentration variability was controlled by the semi-annual pattern, which responded to coastal upwelling conditions during cold period, while rivers outflows influenced high Chlo concentration during the warm period. These results highlighted the importance for land-ocean interface to sustain coastal ecosystems' biological production and the major role of watersheds as sources of nutrients to maintain high biological production during warm periods in transitional tropical-subtropical zones.     

Altimetry,bathymetry and geoid variations at the Gavdos permanent Cal/Val facility

The aim of this work has been to examine the relationship of steep bathymetry in the coastal areas around the permanent Cal/Val facility of Gavdos, and their influence on the produced calibration values for the Jason-2 satellite altimeter. The paper describes how changes in seafloor topography (from 200 to 3500 m depth over a distance of 10 km) are reflected on the determined altimeter parameters using different reference surfaces for satellite calibration. Finally, it describes the relation between these parameter trends and the region’s local characteristics.     

Validation of non-linear split window algorithm for land surface temperature estimation using Sentinel-3 satellite imagery: Case study; Tehran Province,Iran

In recent years, land surface temperature (LST) has become critical in environmental studies and earth science. Remote sensing technology enables spatiotemporal monitoring of this parameter on large scales. This parameter can be estimated by satellite images with at least one thermal band. Sentinel-3 SLSTR data provide LST products with a spatial resolution of 1 km. In this research, direct and indirect validation procedures were employed to evaluate the Sentinel-3 SLSTR LST products over the study area in different seasons from 2018 to 2019. The validation method was based on the absolute (direct) evaluation of this product with field data and comparison (indirect) evaluation with the MODIS LST product and the estimated LST using the non-linear split-window (NSW) algorithm. Also, two emissivity estimation methods, (1) NDVI thresholding method (NDVI-THM) and (2) classification-based emissivity method (CBEM), were used to estimate the LST using the NSW method according to the two thermal bands of Sentinel-3 images. Then, the accuracy of these methods in estimating LST was evaluated using field data and temporal changes of vegetation, which the NDVI-THM method generated better results. For indirect evaluation between the Sentinel-3 LST product, MODIS LST product, and LST estimated using NSW, four filters based on spatial and temporal separates between pairs of pixels and pixel quality were used to ensure the accuracy and consistency of the compared pairs of a pixel. In general, the accuracy results of the LST products of MODIS and Sentinel-3, and LST estimated using NSW showed a similar trend for LST changes during the seasons. With respect to the two absolute and comparative validations for the Sentinel-3 LST products, summer with the highest values of bias (?1.24 K), standard deviation (StDv = 2.66 K), and RMSE (2.43 K), and winter with the lowest ones (bias of 0.14 K, StDv of 1.13 K, and RMSE of 1.12 K) provided the worst and best results for the seasons in the period of 2018–2019, respectively. According to both absolute and comparative evaluation results, the Sentinel-3 SLSTR LST products provided reliable results for all seasons on a large temporal and spatial scale over our studied area.     

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 improved altimeter-derived gravity anomaly from shipborne gravity based on the mean sea surface height constraint factors method

Coastal marine gravity modeling faces challenges due to the degradation of the quality and poor coverage of altimeter data in coastal regions. The effective fusion of shipborne gravity data and altimeter-derived marine gravity data can make shipborne gravity data more useful for the accurate estimation of altimeter-derived coastal marine gravity. A mean sea surface height constraint factor (MSSHCF) method based on the ordinary kriging method and the remove-restore technique is proposed to fuse altimeter-derived gravity model with shipborne gravity data. In this method, all data are standardized during the interpolation process to reduce the error and mean sea surface as a vertical variable is added to the semi-variance function in ordinary kriging to obtain the residual shipborne gravity as corrected data source. The coastal marine gravity models V2.1 and V3.1 which fused altimeter-derived gravity data with shipborne gravity data and V1.1 without shipborne gravity data at a spatial resolution of 1′×1′ can be obtained. Validation experiments show that the accuracy of the gravity model V3.1 obtained by the MSSHCF method more closely agrees with the validated gravity model DTU17 and SS V31 than the model V2.1 obtained by the ordinary kriging interpolation method and the V1.1 model. Our results were validated against shipborne gravity data; the accuracy of model V3.1 was 4.95 % higher than the model V1.1 in South China Sea area A and 2.48 % higher in South China Sea area B. Meanwhile, the accuracy of model V3.1 was 2.07 % higher than model V2.1 in South China Sea area A and 2.42 % higher in South China Sea area B. The effects of distance from the coast and sea depth on the marine gravity model were also evaluated. The results show that the gravity model V3.1 has higher accuracy with the change in ocean distance and depth than the V2.1 and V1.1 gravity models. Thus, our study shows that the MSSHCF method effectively refines coastal altimeter-derived gravity using shipborne gravity data.     

Physical assessment of coastal vulnerability under enhanced land subsidence in Semarang,Indonesia, using multi-sensor satellite data

Research has been conducted in Semarang, Indonesia, to assess coastal vulnerability under enhanced land subsidence using multi-sensor satellite data, including the Advanced Land Observing Satellite (ALOS) Phased Array type L-band SAR (PALSAR), Landsat TM, IKONOS, and TOPEX/Poseidon. A coastal vulnerability index (CVI) was constructed to estimate the level of vulnerability of a coastline approximately 48.68?km in length using seven physical variables, namely, land subsidence, relative sea level change, coastal geomorphology, coastal slope, shoreline change, mean tidal range, and significant wave height. A comparison was also performed between a CVI calculated using seven parameters and a CVI using six parameters, the latter of which excludes the land subsidence parameter, to determine the effects of land subsidence during the coastal vulnerability assessment. This study showed that the accuracy of coastal vulnerability was increased 40% by adding the land subsidence factor (i.e., CVI 6 parameters?=?53%, CVI 7 parameters?=?93%). Moreover, Kappa coefficient indicated very good agreement (0.90) for CVI 7 parameters and fair agreement (0.3) for CVI 6 parameters. The results indicate that the area of very high vulnerability increased by 7% when land subsidence was added. Hence, using the CVI calculation including land subsidence parameters, the very high vulnerability area is determined to be 20% of the total coastline or 9.7?km of the total 48.7?km of coastline. This study proved that land subsidence has significant influence on coastal vulnerability in Semarang.     

Comprehensive comparison of VHR 3D spatial data acquired from IKONOS and TerraSAR-X imagery

With the launch of very high resolution (VHR) optical and synthetic aperture radar (SAR) remote sensing satellites such as IKONOS and TerraSAR-X, a new era has begun in 3D spatial data acquirement. IKONOS provided the first VHR data is still being preferred for many remote sensing applications. TerraSAR-X is considered as a revolution in SAR imaging as a result of 1 m resolution imaging capability. The imaging principles of these satellites are quite different and include advantages and disadvantages that have considerable effects on the quality of acquired 3D spatial data.     

Deramp range migration processing for space-borne spotlight synthetic aperture radar

The deramp range migration algorithm (DRMA) for space-borne spotlight synthetic aperture radar (SAR) is presented, which is based on the two-step algorithm and range migration algorithm. The presented algorithm combines the advantages of SPECAN algorithm and range migration algorithm. Firstly, the presented algorithm performs a bulk azimuth raw data compression to achieve a pixel spacing no larger than the expected azimuth resolution of the fully focused image. Therefore, the azimuth spectral folding phenomenon is overcome. Secondly, the residual and precise focusing of the synthetic aperture radar data is achieved by applying the range migration algorithm. The space-borne SAR focusing parameters’ computation is also described in the paper. At the end of this paper, results are presented which confirm the validity of the presented algorithm.     

星载SAR图像斑点抑制与质量评估系统

结合星载合成孔径雷达(SAR)单视图像斑点噪声抑制技术与质量评估,并基于为科学数据可视化及数学分析提供完备计算环境的IDL(Interactive Data Language)语言,完成了在IDL环境下开发的SAR图像斑点抑制与质量评估系统,获得了高空间分辨率和高辐射分辨率的星载SAR图像.      

高分辨率多模式星载SAR系统设计方法

针对当前星载合成孔径雷达(SAR,Synthetic Aperture Radar)朝高分辨率、多模式方向快速发展的趋势,兼顾星载SAR系统多种工作模式的特点,总结星载SAR系统的星地一体化指标体系;在此基础上提出了高分辨率、多模式星载SAR系统星地一体化设计方法,给出了系统星地一体化设计流程,对雷达天线参数设计、雷达系统参数设计、雷达发射功率选择进行重点分析;以国外高分辨率、多模式星载SAR系统为例,以其总体技术指标及轨道参数为依据,设计了系统的中央电子设备参数、天线参数及波位参数,并分析了系统的性能,从而对所提出的设计方法进行了验证.该方法克服了现有方法的局限性,为高分辨率、多模式星载SAR系统研制提供了技术支撑.