Retrieval of Antarctic sea ice freeboard and thickness from HY-2B satellite altimeter data

Antarctic sea ice is an important part of the Earth’s atmospheric system,and satellite remote sensing is an important technology for observing Antarctic sea ice.Whether Chinese Haiyang-2B(HY-2B)satellite altimeter data could be used to estimate sea ice freeboard and provide alternative Antarctic sea ice thickness information with a high precision and long time series,as other radar altimetry satellites can,needs further investigation.This paper proposed an algorithm to discriminate leads and then retrieve sea ice freeboard and thickness from HY-2B radar altimeter data.We first collected the Moderate-resolution Imaging Spectroradiometer ice surface temperature(IST)product from the National Aeronautics and Space Administration to extract leads from the Antarctic waters and verified their accuracy through Sentinel-1 Synthetic Aperture Radar images.Second,a surface classification decision tree was generated for HY-2B satellite altimeter measurements of the Antarctic waters to extract leads and calculate local sea surface heights.We then estimated the Antarctic sea ice freeboard and thickness based on local sea surface heights and the static equilibrium equation.Finally,the retrieved HY-2B Antarctic sea ice thickness was compared with the CryoSat-2 sea ice thickness and the Antarctic Sea Ice Processes and Climate(ASPeCt)ship-based observed sea ice thickness.The results indicate that our classification decision tree constructed for HY-2B satellite altimeter measurements was reasonable,and the root mean square error of the obtained sea ice thickness compared to the ship measurements was 0.62 m.The proposed sea ice thickness algorithm for the HY-2B radar satellite fills a gap in this application domain for the HY-series satellites and can be a complement to existing Antarctic sea ice thickness products;this algorithm could provide long-time-series and large-scale sea ice thickness data that contribute to research on global climate change.

Validation of Chinese HY-2 satellite radar altimeter significant wave height

Chinese Haiyang-2（HY-2） satellite is the first Chinese marine dynamic environment satellite. The dual-frequency （Ku and C band） radar altimeter onboard HY-2 has been working effective to provide operational significant wave height （SWH） for more than three years （October 1, 2011 to present）.We validated along-track Ku-band SWH data of HY-2 satellite against National Data Buoy Center （NDBC） in-situ measurements over a time period of three years from October 1, 2011 to September 30, 2014, the root mean square error （RMSE） and mean bias of HY-2 SWH is 0.38 m and （-0.13±0.35） m, respectively. We also did cross validation against Jason-2 altimeter SWH data, the RMSE and the mean bias is 0.36m and （-0.22±0.28） m, respectively. In order to compare the statistical results between HY-2 and Jason-2 satellite SWH data, we validated the Jason-2 satellite radar altimeter along-track Ku-band SWH data against NDBC measurements using the same method. The results demonstrate the validation method in this study is scientific and the RMSE and mean bias of Jason-2 SWH data is 0.26 m and （0.00±0.26） m, respectively. We also validated both HY-2 and Jason-2 SWH data every month, the mean bias of Jason-2 SWH data almost equaled to zero all the time, while the mean bias of HY-2 SWH data was no less than -0.31m before April 2013 and dropped to zero after that time. These results indicate that the statistical results for HY-2 altimeter SWH are reliable and HY-2 altimeter along-track SWH data were steady and of high quality in the last three years. The results also indicate that HY-2 SWH data have greatly been improved and have the same accuracy with Jason-2 SWH data after April, 2013. SWH data provided by HY-2 satellite radar altimeter are useful and acceptable for ocean operational applications.

Absolute calibration of HY-2A and Jason-2 altimeters for sea surface height using GPS buoy in Qinglan, China

GPS buoy methodology is one of the main calibration methodologies for altimeter sea surface height calibration. This study introduces the results of the Qinglan calibration campaign for the HY-2A and Jason-2 altimeters. It took place in two time slices;one was from August to September 2014, and the other was in July 2015. One GPS buoy and two GPS reference stations were used in this campaign. The GPS data were processed using the real-time kinematic (RTK) technique. The fi nal error budget estimate when measuring the sea surface height (SSH) with a GPS buoy was better than 3.5 cm. Using the GPS buoy, the altimeter bias estimate was about -2.3 cm for the Jason-2 Geophysical Data Record (GDR) Version ‘D' and from -53.5 cm to -75.6 cm for the HY-2A Interim Geophysical Data Record (IGDR). The bias estimates for Jason-2 GDR-D are similar to the estimates from dedicated calibration sites such as the Harvest Platform, the Crete Site and the Bass Strait site. The bias estimates for HY-2A IGDR agree well with the results from the Crete calibration site. The results for the HY-2A altimeter bias estimated by the GPS buoy were verifi ed by cross-calibration, and they agreed well with the results from the global analysis method.

Preliminary calibration results of the HY-2B altimeter’s SSH at China’s Wanshan calibration site

Satellite altimeter needs to be calibrated to evaluate the accuracy of sea surface height data.The dedicated altimeter calibration field needs to establish a special calibration strategy and needs to evaluate its calibration ability.This paper describes absolute calibration of HY-2 B altimeter SSH using the GPS calibration method at the newly Wanshan calibration site,located in the Wanshan Islands,China.There are two HY-2 B altimeter passes across the Wanshan calibration site.Pass No.362 is descending and the ground track passes the east of Dan’gan Island.Pass No.375 is ascending and crosses the Zhiwan Island.The GPS data processing strategy of Wanshan calibration site was established and the accuracy of GPS calibration method of Wanshan calibration site was evaluated.Meanwhile,the processing strategies of the HY-2 B altimeter for the Wanshan calibration site were established,and a dedicated geoid model data were used to benefit the calibration accuracy.The time-averaged HY-2 B altimeter bias was approximately 2.12 cm with a standard deviation of 2.08 cm.The performance of the HY-2 B correction microwave radiometer was also evaluated in terms of the wet troposphere path delay and showed a mean difference-0.2 cm with a 1.4 cm standard deviation with respect to the in situ GPS radiosonde.

Measurement analyses and evaluations of sea-level heights using the HY-2A satellite’s radar altimeter

The HY-2A satellite is China’s first independent oceanic dynamic environmental satellite,and has been operating continuously for more than six years.The satellite’s radar altimeter,which is one of the main loads on the satellite,has the ability to realize all-weather and all-day observations of global sea-surface heights,as well as significant wave heights and sea-surface wind speeds.These observed data have been widely used in marine disaster prevention and reduction,along with resource development,maritime security and other fields.In order to achieve a comprehensive understanding of the multi-year overall observational performances of the HY-2A satellite’s radar altimeter,all of the observational data of the IGDR product from October 26,2012 to August 27,2017 were selected in this study for a comprehensive evaluation.The height measurement capability of the HY-2A satellite’s radar altimeter was evaluated using self-crossover and Jason-2 crossover methods.The height discrepancies at the self-crossover point of the HY-2A satellite’s ascending and descending orbits were also calculated.It was found that for the HY-2A satellite’s radar altimeter in global waters under the restriction conditions of ascending and descending orbits,the height anomaly differences were within a range of less than 30 cm.The absolute mean error was determined to be 5.81 cm,and the height anomaly standard deviation was 7.76 cm.Under the conditions of the observational areas being limited within a scope of 60°from the Equator,it was determined that the sea-level height anomaly differences were less than 10 cm at the junction of the ascending and descending orbits,the absolute mean error was 3.95 cm.In addition,the sea-level height anomaly standard deviation was observed to be 4.76cm.Using a mutual cross method with the Jason-2 satellite,it was found that under the conditions of the observational area being within the scope of 66°from the equator,the height anomaly differences at the junction were less than 30cm,and the absolute mean error of HY-2A and Jason-2 sea level height anomaly was 5.86 cm,with a standard deviation of 7.52 cm.It was observed that,if within the sea area the sea level height anomaly difference was limited to within 10cm,then the absolute mean error and standard deviation could reach 4.19cm and 4.98cm,respectively.It was confirmed that the HY-2A satellite’s radar altimeter had successfully reached the height measurement level of similar international altimeters.Therefore,it had the ability to meet the needs of marine scientific research and ocean circulation inversions.

Preliminary marine gravity field from HY-2A/GM altimeter data

HY-2 A(Haiyang-2 A),launched in 2011,is the first ocean dynamic environment satellite of China and is equipped with a radar altimeter as one of the primary payloads.HY-2 A shifted the drift orbit in March 2016 and has been accumulating geodetic mission(GM)data for more than three years with 168-day cycle.In this paper,we present the preliminary gravity field inverted by the HY-2 A/GM data from March 2016 to December 2017 near Taiwan(21°–26°N,119°–123°E).The gravity anomaly is computed by Inverse Vening Meinesz(IVM)formula with a onedimensional FFT method during remove-restore procedure with the EGM2008 gravity model as the reference field.For comparison,CryoSat-2 altimeter data are used to inverse the gravity field near Taiwan Island by the same method.Comparing with the gravity field derived from CryoSat-2,a good agreement between the two data sets is found.The global ocean gravity models and National Geophysical Data Center(NGDC)shipboard gravity data also are used to assess the performance of HY-2 A/GM data.The evaluations show that HY-2 A and CryoSat-2 are at the same level in terms of gravity field recovery and the HY-2 A/GM altimeter-derived gravity field has an accuracy of 2.922 mGal.Therefore,we can believe that HY-2 A will be a new reliable data source for marine gravity field inversion and has the potentiality to improve the accuracy and resolution of the global marine gravity field.

Assessment of reprocessed sea surface height measurements derived from HY-2A radar altimeter and its application to the observation of 2015–2016 El Nino

Haiyang-2A（HY-2A） is China＇s first ocean dynamic environment satellite and the radar altimeter is one of its main payloads. One of the main purposes of the radar altimeter is to measure the sea surface height（SSH）. The SSH determined from the altimeter range measurements includes some range and geophysical corrections. These corrections largely affect the accuracy of the SSH measurements. The range and the geophysical corrections are reprocessed and the altimeter waveforms in HY-2A sensor interim geophysical data set records（S-IGDR） are retracked from June 1, 2014 to June 14, 2014, and the accuracy of the reprocessed SSH measurements is evaluated.The methods of the range and geophysical corrections used to reprocess HY-2A altimeter data are validated by using these methods to reprocess the Jason-2 range and geophysical corrections and comparing the results with the range and geophysical corrections in Jason-2 geophysical dataset records（GDR） product. A crossover analysis is used to evaluate the accuracy of the reprocessed HY-2A SSH measurements. The standard deviation（STD） of the crossover SSH differences for HY-2A is around 4.53 cm while the STD of the SSH differences between HY-2A and Jason-2 is around 5.22 cm. The performance of the reprocessed HY-2A SSH measurements is significantly improved with respect to the SSH measurements derived from HY-2A interim geophysical dataset records（IGDR）product. The 2015–2016 El Ni？o has been the strongest El Ni？o event since 1997–1998. The range and the geophysical corrections in HY-2A IGDR are reprocessed and sea level anomalies are used to monitor the2015–2016 El Ni？o. The results show that the HY-2A altimeter can well observe the 2015–2016 El Ni？o.

Validation of the “HY-2” altimeter wet tropospheric path delay correction based on radiosonde data

Wet tropospheric path delay （PD） is a highly variable term for the altimeter measurement of a sea surface height, caused by the refraction effect of atmospheric water vapor and cloud liquid water. In order to esti- mate PD values, the ＂HY-2＂ system includes a calibration microwave radiometer （CMR） operating at 18.7, 23.8 and 37 GHz. The PD data of the CMR were compared and validated by coincident radiosonde profiles from ten globally distributed radiosonde stations during October 2011 to August 2012. The temporal interval was 1 h. In order to avoid land contamination, different spatial intervals between these two data sets were tested. The empirical fit function of PD uncertainty and spatial interval was found and extrapolated to the ideal situation that the data of CMR and radiosonde were totally coincident. The stability of the brightness temperature of the CMR and its impact on the PD correction was also studied. Consequently, the uncertainty of the PD algorithm of the CMR was estimated to be 2.1 cm.

Preliminary results of the global ocean tide derived from HY-2A radar altimeter data

The HY-2A satellite,which is equipped with a radar altimeter and was launched on August 16,2011,is the first Chinese marine dynamic environmental monitoring satellite.Extracting ocean tides is one of the important applications of the radar altimeter data.The radar altimeter data of the HY-2A satellite from November 1,2011 to August 16,2014 are used herein to extract global ocean tides.The constants representing the tidal constituents are extracted by HY-2A RA data with harmonic analysis based on the least squares method.Considering tide aliasing issues,the analysis of the alias periods and alias synodic periods of different tidal constituents shows that only the tidal constituents M_(2),N_(2),and K_(2)are retrieved precisely by the HY-2A RA data.The derived tidal constants of the tidal constituents M_(2),N_(2)and K_(2)are compared to those of tidal gauge data and the TPXO tide model results.The comparison between the derived results and the tidal gauge data shows that the RMSEs of the tidal amplitude and phase lag are 9.6 cm and 13.34°,2.4 cm and 10.47°,and 8.1 cm and 14.19°for tidal constituents M_(2),N_(2),and K_(2),respectively.The comparisons of the semidiurnal tides with the TPXO model results show that tidal constituents have good consistency with the TPXO model results.These findings confirm the good performance of HY-2A RA for retrieving semidiurnal tides in the global ocean.

HY-2C卫星散射计海面风场产品的检验评估

基于中国近海浮标的逐小时海面风矢量实测数据、ASCAT散射计风场数据以及欧洲中期天气预报中心的第五代全球再分析数据(ERA5)的逐小时风场数据,使用统计分析方法对HY-2C卫星散射计海面风场数据产品进行检验评估,检验使用的是2~24 m/s风速区间内的数据。结果表明:使用11个浮标共计3351个样本数据作为真值进行检验,风速和风向的均方根误差分别为1.01 m/s和20.7°;使用8304963个ASCAT样本数据进行检验,风速和风向的均方根误差分别为0.66 m/s和11.5°;使用ERA5风场数据进行全球区域检验,风速和风向的均方根误差分别为1.05 m/s和12.0°;使用ERA5风场数据进行西北太平洋区域检验,风速和风向的均方根误差分别为1.19 m/s和15.5°。由此可见,HY-2C卫星散射计海面风场产品质量的可信度较高,能够较好地满足业务化应用的精度要求。

HY-2D卫星散射计海面风场产品在西北太平洋区域的适用性分析

基于中国近海浮标数据、欧洲中期天气预报中心第五代全球再分析资料ERA5和MetOp-B卫星散射计ASCAT海面风场数据,对HY-2D卫星散射计反演的风场数据在西北太平洋区域的质量进行了检验评估,用以分析其在该区域的适用性。统计结果表明:HY-2D风场数据与中国近海浮标数据、ERA5数据和ASCAT数据均有较好的一致性,3组检验结果的风速均方根误差分别为0.98 m/s、1.23 m/s和0.58 m/s,误差都很好地控制在2 m/s以内;风向均方根误差分别为18.8°、16.8°和12.8°,误差均控制在20°以内;HY-2D风场数据与ASCAT数据有着最好的一致性。对比HY-2D风场数据和中国近海浮标站点风矢量数据,未发现离岸距离和误差之间的明显规律。在HY-2D风场数据与ERA5和ASCAT的对比结果中,近9成数据的风速绝对误差≤2 m/s、风向绝对误差≤20°。由此可见,HY-2D卫星散射计反演的海面风场数据在西北太平洋区域质量较高,能够较好地满足业务化应用的精度要求。

Performance evaluation of HY-2 series satellites in marine gravity field recovery

China has successfully launched four Haiyang-2(HY-2)series altimetry satellites.HY-2A has attracted significant attention in gravity field recovery,but the performance of other HY-2 series satellites,including HY-2B/C/D,is seldom discussed.This study evaluated the performance of all the HY-2 series satellites in recovering marine gravity field.First,the crossover discrepancies in sea surface height of the four satellites,HY-2A,HY-2B,HY-2C,and HY-2D,were analyzed to assess their altimetry stability.It was found that HY-2B had the best altimetry quality,followed by HY-2D.Subsequently,different combina-tions of altimetry data were used to calculate vertical deflections and gravity anomalies in the South China Sea(112°E-119°E,12°N-20°N).The results showed that combining data from HY-2B,HY-2C,and HY-2D improved the inversion accuracy of gravity anomalies by 0.3 mGal compared to using HY-2A data alone.HY-2C and HY-2D contributed to enhancing the accuracy of the east component of vertical deflections.

Three-dimensional thermohaline structure estimation derived from HY-2 satellite data over the Maritime Silk Road and its applications

Estimated ocean subsurface fields derived from satellite observations provide potential data sources for operational marine environmental monitoring and prediction systems.This study employs a statistic regression reconstruction method,in combination with domestic autonomous sea surface height and sea surface temperature observations from the Haiyang-2(HY-2)satellite fusion data,to establish an operational quasi-realtime three-dimensional(3D)temperature and salinity products over the Maritime Silk Road.These products feature a daily temporal resolution and a spatial resolution of 0.25°×0.25°and exhibit stability and continuity.We have demonstrated the accuracy of the reconstructed thermohaline fields in capturing the 3D thermohaline variations through comprehensive statistical evaluations,after comparing them against Argo observations and ocean analysis data from 2022.The results illustrate that the reconstructed fields effectively represent seasonal variations in oceanic subsurface structures,along with structural changes resulting from mesoscale processes,and the upper ocean’s responses to tropical cyclones.Furthermore,the incorporation of HY-2 satellite observations notably enhances the accuracy of temperature and salinity reconstructions in the Northwest Pacific Ocean and marginally improves salinity reconstruction accuracy in the North Indian Ocean when compared to the World Ocean Atlas 2018 monthly climatology thermohaline fields.As a result,the reconstructed product holds promise for providing quasi-real-time 3D temperature and salinity field information to facilitate fast decisionmaking during emergencies,and also offers foundational thermohaline fields for operational ocean reanalysis and forecasting systems.These contributions enhance the safety and stability of ocean subsurface activities and navigation.

Estimates of global M_2 internal tide energy fluxes using TOPEX/POSEIDON altimeter data

TOPEX/POSEIDON altimeter data from October 1992 to June 2002 are used to calculate the global barotropic M2 tidal currents using long-term tidal harmonic analysis. The tides calculated agree well with ADCP data obtained from the South China Sea （SCS）. The maximum tide velocities along the semi-major axis and semi-minor axis can be computed from the tidal ellipse. The global distribution of M2 internal tide vertical energy flux from the sea bottom is calculated based on a linear internal wave generation model. The global vertical energy flux of M2 internal tide is 0.96 TW, with 0.36 TW in the Pacific, 0.31 TW in the Atlantic and 0.29 TW in the Indian Ocean, obtained in this study. The total horizontal energy flux of M2 internal tide radiating into the open ocean from the lateral boundaries is 0.13 TW, with 0.06 TW in the Pacific, 0.04TW in the Atlantic, and 0.03 TW in the Indian Ocean. The result shows that the principal lunar semi-diurnal tide Me provides enough energy to maintain the large-scale thermohaline circulation of the ocean.

HY-2雷达高度计和微波散射计

海洋动力环境卫星(HY-2)是一种获取海洋动力环境信息的专用对地遥感卫星,计划中的我国第一颗海洋动力环境卫星(HY-2A)以微波遥感器为主,包括一台雷达高度计、一台微波散射计和一台多通道微波辐射计。介绍HY-2雷达高度计和微波散射计的系统设计和主要技术指标。

HY-2卫星雷达高度计时标偏差估算

卫星雷达高度计是海洋遥感监测的重要传感器之一,测高系统和定轨系统是高度计重要的组成部分。若两系统使用不同的系统时钟,则获得的轨道高度和卫星测距值之间可能会存在一个时标偏差,该时标偏差会降低卫星雷达高度计的海面高度测量精度。针对HY-2卫星雷达高度计的时标偏差问题,本文分析了时标偏差对测高精度的影响,介绍了一种使用自交叉点数据估算时标偏差值的方法,并基于HY-2卫星雷达高度计第21个周期数据开展了时标偏差修正实验。时标偏差修正后HY-2自交叉点的海面高度差值(也称"不符值")分布收敛程度有了明显的提高,其RMS均方根值从24.7cm减小到了7.0cm,HY-2与Jason-2互交叉点的不符值的RMS也从16.6cm减小到了7.3cm。这表明本文介绍的时标偏差修正方法可有效地提高HY-2卫星雷达高度计的测高精度。

HY-2A卫星雷达高度计数据的全球统计评价及质量分析

自HY-2A卫星发射以来,针对HY-2A卫星雷达高度计产品的交叉定标、真实性检验及质量评估工作一直在持续开展。本文主要以HY-2A卫星高度计第44周期的IGDR产品数据为例,通过使用全球分布图、二维直方图和每日均值统计的方法完成了与Jason-2IGDR产品的比对验证,同时对主要环境校正参数及地球物理产品的数据质量稳定性进行了分析,结果显示高度计产品数据质量较稳定,此外利用HY-2A卫星升降轨交叉点海面高度差、与Jason-2卫星交叉点海面高度差以及沿轨海平面异常数据分析的方法进行了HY-2A卫星高度计观测系统的性能评估,结果显示,HY-2A卫星海面高度精度约为7.48cm,精度接近Jason-2,能满足海洋应用与科学研究的需要。

HY-2卫星雷达高度计风速反演验证

本文给出了海洋二号(HY-2)卫星雷达高度计风速反演的双参数算法模型,并利用Jason-1产品和HY-2卫星微波散射计产品定性验证了双参数算法模型的正确性,同时利用星星交叉比对方法、美国国家浮标数据中心(National Data Buoy Center,NDBC)浮标数据验证方法、美国国家环境预报中心(National Centers for Environmental Prediction,NCEP)再分析数据验证方法三种不同方法对该模型进行定量化验证,结果表明HY-2卫星雷达高度计风速反演的双参数算法模型的反演精度达到2.0 m/s的精度,满足海洋业务应用和科学研究的精度要求,可为用户提供满足观测要求的海面风产品。

HY-2卫星散射计热带气旋自动识别算法

对基于HY-2卫星散射计风矢量产品的热带气旋自动识别算法进行了研究。算法分为粗搜索与精搜索两部分。粗搜索利用热带气旋风场的风速与风向分布直方图特征确定搜索的阈值,快速剔除比较容易识别的非热带气旋区域。在此基础上,精搜索利用热带气旋风向的螺旋状分布特征,通过搜索目标区域内是否存在螺旋状流线的方法,确定目标区域的风向是否存在螺旋状流线特征,从而实现对热带气旋的准确自动识别。作为示例,将该方法应用到对HY-2散射计观测到的2012年6号强热带风暴"杜苏芮"的自动识别,结果表明,本文提出的算法可以从HY-2散射计风场数据中准确有效的自动识别出热带气旋。

基于GNSS浮标和验潮资料的HY-2A卫星高度计绝对定标

为探测我国HY-2A卫星高度计海面高测量绝对偏差及其在轨运行状态,本文利用GNSS浮标星下点同步测量和验潮资料海面高传递方法在山东千里岩和珠海担杆岛海域开展定标研究。为验证GNSS浮标定标方法的准确性,还对国外卫星Jason-2和Saral进行了定标实验。实验表明GNSS浮标绝对海面高测量精度达2cm,对Jason-2和Saral高度计多个周期定标得到的海面高偏差均值分别为5.7cm和-2.3cm,与国际专门定标场的结果符合较好。2014年9月和2015年5月HY-2A卫星高度计浮标定标结果分别是-65cm和-91cm,因两次结果差异显著,故又利用千里岩验潮站资料对HY-2A卫星高度计第56至73周期进行了定标分析,结果证明HY-2A卫星海面高存在约-51cm/a的漂移,置信度为95%的回归分析表明浮标和验潮定标结果符合。本文研究结果表明在我国尚无专门定标场的情况下,可利用GNSS浮标对我国高度计实施灵活、精准的在轨绝对定标,在有高度计轨迹经过验潮站的情况下可使用验潮资料结合精密大地水准面模型进行绝对定标。