星载SAR图像方位角计算及其对风速反演影响评估

星载合成孔径雷达（Synthetic Aperture Radar, SAR）是海洋常规观测的重要技术手段,然而在其数据处理和关键物理参数反演中,SAR图像方位角$ {phi }_{mathrm{I}mathrm{M}mathrm{G}} $和卫星飞行方向角$ {phi }_{mathrm{S}mathrm{A}mathrm{T}} $相混淆的情况普遍存在,影响SAR海面风速的反演精度。对此,根据SAR图像方位角的定义,结合空间卫星坐标系转换关系推导了$ {phi }_{mathrm{I}mathrm{M}mathrm{G}} $,并进一步提出了基于SAR图像地面控制点的$ {phi }_{mathrm{I}mathrm{M}mathrm{G}} $计算方法。通过对比分析2016年全球近70万景哨兵一号卫星波模式SAR图像,量化了$ {phi }_{mathrm{I}mathrm{M}mathrm{G}} $与$ {phi }_{mathrm{S}mathrm{A}mathrm{T}} $的系统偏差,发现该偏差与雷达入射角和轨道方向紧密相关。选用五次多项式拟合该偏差随纬度的变化规律,同时发现将卫星飞行方向角近似等于SAR图像方位角会引起海面风场反演误差,该误差空间分布不均,且升轨降轨有所差异,风速误差最大可达 ±0.5 m·s–1。研究结果表明,$ {phi }_{mathrm{I}mathrm{M}mathrm{G}} $的正确计算和使用对于SAR地球科学研究具有重要意义,本文提出的$ {phi }_{mathrm{I}mathrm{M}mathrm{G}} $计算方法对其他SAR卫星系列也具有实际参考价值。     

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.     

Interaction of Rossby waves with the Gulf Stream and Kuroshio using altimetry in a framework of a vortex layer model

In this paper, we show the potential of satellite altimetry to study the interaction of Rossby waves with the shear flow. The Miles-Ribner approach, which was developed in gas dynamics in the 1960 s, is used to describe Rossby waves interacting with the Gulf Stream and Kuroshio areas. The region of interaction is approximated by a nonzonal vortex layer. We apply the main formulations of the problem of a nonzonal vortex layer on the β-plane in the formulation of Miles-Ribner to observations in the real ocean. Earlier, we showed that the interaction of waves with a nonzonal flow gives rise to a new class of solutions, which is absent in the case of a zonal flow. This new class of solutions can be interpreted as the pure emission of Rossby waves by the nonzonal flow. We apply this theoretical approach to the areas of the Gulf Stream and Kuroshio as well. We use for analysis altimetry data available at Copernicus Marine Environment Monitoring Service. The analysis of Hovmöller diagrams in the areas under consideration confirms the previously obtained theoretical conclusions of the problem of the interaction of planetary waves with a nonzonal flow on the β-plane in the formulation of Miles-Ribner. The incident waves, as well as refracted and reflected waves are distinguished. The speed of refracted and reflected waves exceeds the speed of incident waves, which confirms the conclusions about the existence of mechanisms for the amplification of planetary waves when they interact with a nonzonal flow.     

GNSS-R open-loop difference phase altimetry: Results from a bridge experiment

The paper explores a method to obtain accurate lake surface heights using measurements of the Global Navigation Satellite System (GNSS) carrier phase reflected from the lake surface. The method is referred to as Global Navigation Satellite System-Reflection (GNSS-R) open-loop difference phase altimetry method. It consists of two key technologies: one is the open-loop tracking method to track the GNSS-R signals, where the direct GNSS signal’s frequency is used as a reference frequency to obtain the carrier phases of the GNSS-R signals; the other key technology is time difference phase altimetry method to invert the lake surface heights using two or more carrier phases of GNSS-R signals received simultaneously. A validation experiment is carried out on the SANYING bridge over GUANTING lake using a GNSS-R receiver developed by the Center for Space Science and Applied Research (CSSAR), processing the data with GNSS-R open-loop difference phase altimetry method. The lake surface height results are consistent with the height results of GPS dual-frequency differential positioning altimetry. The results show that we can achieve centimeter level height in one minute average, by using 11 minutes carrier phase data of three GNSS-R signals received simultaneously.     

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.     

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.     

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).     

The multi-source data fusion global ionospheric modeling software—IonoGim

We introduce a new global ionospheric modeling software—IonoGim, using ground-based GNSS data, the altimetry satellite and LEO (Low Earth Orbit) occultation data to establish the global ionospheric model. The software is programmed by C++ with fast computing speed and highly automatic degree, it is especially suitable for automatic ionosphere modeling. The global ionospheric model and DCBs obtained from IonoGim were compared with the CODE (Center for Orbit Determination in Europe) to verify its accuracy and reliability. The results show that IonoGim and CODE have good agreement with small difference, indicating that IonoGim owns high accuracy and reliability, and can be fully applicable for high-precision ionospheric research. In addition, through comparison between only using ground-based GNSS observations and multi-source data model, it can be demonstrated that the space-based ionospheric data effectively improve the model precision in marine areas where the ground-based GNSS tracking station lacks.     

An improved,near real time water vapor path delay correction of Cryosat-2 sea surface height measurements

This paper presents improvements of a method (Stum et al., 2011) aimed at computing the water vapor path delay correction of altimeter sea surface height, using total precipitable water measurements from scanning microwave radiometers. The main interest of this improved method is for the Cryosat-2 mission over the ocean. Focus is made on the applicability of the method in near real time. An experiment to produce an operational path delay correction for Jason-2 and Cryosat-2 Interim Geophysical Data Records (IGDR) has been set up. Results confirm that the new correction, although less accurate than the one attainable with an embarked radiometer, improves the Cryosat-2 sea surface height accuracy.     

The impact on tsunami detection from space using GNSS-reflectometry when combining GPS with GLONASS and Galileo

The devastating Sumatra tsunami in 2004 demonstrated the need for a tsunami early warning system in the Indian Ocean. Such a system has been installed within the German-Indonesian Tsunami Early Warning System (GITEWS) project. Tsunamis are a global phenomenon and for global observations satellites are predestined. Within the GITEWS project a feasibility study on a future tsunami detection system from space has therefore been carried out. The Global Navigation Satellite System Reflectometry (GNSS-R) is an innovative way of using GNSS signals for remote sensing. It uses ocean reflected GNSS signals for sea surface altimetry. With a dedicated Low Earth Orbit (LEO) constellation of satellites equipped with GNSS-R receivers, densely spaced sea surface height measurements could be established to detect tsunamis. Some general considerations on the geometry between LEO and GNSS are made in this simulation study. It exemplary analyzes the detection performance of a GNSS-R constellation at 900 km altitude and 60° inclination angle when applied to the Sumatra tsunami as it occurred in 2004. GPS is assumed as signal source and the combination with GLONASS and Galileo signals is investigated. It can be demonstrated, that the combination of GPS and Galileo is advantageous for constellations with few satellites while the combination with GLONASS is preferable for constellations with many satellites. If all three GNSS are combined, the best detection performance can be expected for all scenarios considered. In this case an 18 satellite constellation will detect the Sumatra tsunami within 17 min with certainty, while it takes 53 min if only GPS is considered.     

Impact of ITRS 2014 realizations on altimeter satellite precise orbit determination

This paper evaluates orbit accuracy and systematic error for altimeter satellite precise orbit determination on TOPEX, Jason-1, Jason-2 and Jason-3 by comparing the use of four SLR/DORIS station complements from the International Terrestrial Reference System (ITRS) 2014 realizations with those based on ITRF2008. The new Terrestrial Reference Frame 2014 (TRF2014) station complements include ITRS realizations from the Institut National de l’Information Géographique et Forestière (IGN) ITRF2014, the Jet Propulsion Laboratory (JPL) JTRF2014, the Deutsche Geodätisches Forschungsinstitut (DGFI) DTRF2014, and the DORIS extension to ITRF2014 for Precise Orbit Determination, DPOD2014. The largest source of error stems from ITRF2008 station position extrapolation past the 2009 solution end time. The TRF2014 SLR/DORIS complement impact on the ITRF2008 orbit is only 1–2 mm RMS radial difference between 1992–2009, and increases after 2009, up to 5 mm RMS radial difference in 2016. Residual analysis shows that station position extrapolation error past the solution span becomes evident even after two years, and will contribute to about 3–4 mm radial orbit error after seven years. Crossover data show the DTRF2014 orbits are the most accurate for the TOPEX and Jason-2 test periods, and the JTRF2014 orbits for the Jason-1 period. However for the 2016 Jason-3 test period only the DPOD2014-based orbits show a strong and statistically significant margin of improvement. The positive results with DTRF2014 suggest the new approach to correct station positions or normal equations for non-tidal loading before combination is beneficial. We did not find any compelling POD advantage in using non-linear over linear station velocity models in our SLR & DORIS orbit tests on the Jason satellites. The JTRF2014 proof-of-concept ITRS realization demonstrates the need for improved SLR+DORIS orbit centering when compared to the Ries (2013) CM annual model. Orbit centering error is seen as an annual radial signal of 0.4 mm amplitude with the CM model. The unmodeled CM signals show roughly a 1.8 mm peak-to-peak annual variation in the orbit radial component. We find the TRF network stability pertinent to POD can be defined only by examination of the orbit-specific tracking network time series. Drift stability between the ITRF2008 and the other TRF2014-based orbits is very high, the relative mean radial drift error over water is no larger than 0.04 mm/year over 1993–2015. Analyses also show TRF induced orbit error meets current altimeter rate accuracy goals for global and regional sea level estimation.     

Radar altimetry measurements over antarctic ice sheet: A focus on antenna polarization and change in backscatter problems

In this paper, we investigate the impact of the error due to the penetration of the altimetric wave within the snowpack. The phenomenon has two different impacts. The first one, due to temporal change in snow characteristics, affects the ice sheet volume trend as derived from altimetric series. The second one, because of both the anisotropy of the ice sheet surface properties and of the linear antenna polarization, introduces a difference in measurements at crossover points. These two phenomena are the cause of what are probably the most critical limitations to the interpretation of long-term altimetric series in term of mass balance and to the comparison between or data fusion of different missions. Moreover, they will lead to the largest error when comparing data from EnviSat with data from CryoSat, because of the different orbits, or with data from AltiKa, because of the different radar frequencies.     

DORIS-based point mascons for the long term stability of precise orbit solutions

In recent years non-tidal Time Varying Gravity (TVG) has emerged as the most important contributor in the error budget of Precision Orbit Determination (POD) solutions for altimeter satellites’ orbits. The Gravity Recovery And Climate Experiment (GRACE) mission has provided POD analysts with static and time-varying gravity models that are very accurate over the 2002–2012 time interval, but whose linear rates cannot be safely extrapolated before and after the GRACE lifespan. One such model based on a combination of data from GRACE and Lageos from 2002–2010, is used in the dynamic POD solutions developed for the Geophysical Data Records (GDRs) of the Jason series of altimeter missions and the equivalent products from lower altitude missions such as Envisat, Cryosat-2, and HY-2A. In order to accommodate long-term time-variable gravity variations not included in the background geopotential model, we assess the feasibility of using DORIS data to observe local mass variations using point mascons. In particular, we show that the point-mascon approach can stabilize the geographically correlated orbit errors which are of fundamental interest for the analysis of regional Mean Sea Level trends based on altimeter data, and can therefore provide an interim solution in the event of GRACE data loss. The time series of point-mass solutions for Greenland and Antarctica show good agreement with independent series derived from GRACE data, indicating a mass loss at rate of 210 Gt/year and 110 Gt/year respectively.     

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.     

Introduction to global short message communication service of BeiDou-3 navigation satellite system

Short Message Communication (SMC) is a featured service of BeiDou Navigation Satellite System (BDS). After its successful deployment in 2003, Regional Short Message Communication (RSMC) service has been continuously serving China and its neighboring countries and regions, especially in life safety scenarios. In this paper, the architecture of the Global Short Message Communication (GSMC) system is proposed based on the medium earth orbit (MEO) constellation and the crosslinks of the global BeiDou navigation system (BDS-3). Three subtypes of GSMC service, i.e. positioning report, emergency search and rescue (SAR) and regular SMC are designed in accordance with the technical characteristic of integration of navigation and communication in BDS-3, which supports future wide applications of GSMC. The performance of the designed GSMC system is analyzed by numerical calculations. As BDS-3 was officially announced completion on July 31, 2020, GSMC has been providing initial service. First test results of the in orbit GSMC payloads are also presented in the paper to verify the designed capabilities. Preliminary results also show that the requirements of Global Maritime Distress and Safety System (GMDSS) can also be fulfilled.