Journal of Oceanology and Limnology   2022, Vol. 40 issue(5): 1881-1888     PDF       
http://dx.doi.org/10.1007/s00343-020-0044-3
Institute of Oceanology, Chinese Academy of Sciences
0

Article Information

GAO Le, SUN Hanwei, QI Jifeng, JIANG Qiufu
Effect of random phase error and baseline roll angle error on eddy identification by interferometric imaging altimeter
Journal of Oceanology and Limnology, 40(5): 1881-1888
http://dx.doi.org/10.1007/s00343-020-0044-3

Article History

Received Jan. 22, 2020
accepted in principle Mar. 22, 2020
accepted for publication Apr. 26, 2020
Effect of random phase error and baseline roll angle error on eddy identification by interferometric imaging altimeter
Le GAO1,3, Hanwei SUN2, Jifeng QI1, Qiufu JIANG1     
1 CAS Key Laboratory of Ocean Circulation and Waves, Institute of Oceanology, Chinese Academy of Sciences and Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China;
2 Beijing Radio Measurement Institute, Beijing 100854, China;
3 Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China
Abstract: To achieve better observation for sea surface, a new generation of wide-swath interferometric altimeter satellites is proposed. Before satellite launch, it is particularly important to study the data processing methods and carry out the detailed error analysis of ocean satellites, because it is directly related to the ultimate ability of satellites to capture ocean information. For this purpose, ocean eddies are considered a specific case of ocean signals, and it can cause significant changes in sea surface elevation. It is suitable for theoretical simulation of the sea surface and systematic simulation of the altimeter. We analyzed the impacts of random error and baseline error on the sea surface and ocean signals and proposed a combined strategy of low-pass filtering, empirical orthogonal function (EOF) decomposition, and linear fitting to remove the errors. Through this strategy, sea surface anomalies caused by errors were considerably improved, and the capability of satellite for capturing ocean information was enhanced. Notably, we found that the baseline error in sea surface height data was likely to cause inaccuracy in eddy boundary detection, as well as false eddy detection. These abnormalities could be prevented for "clean" sea surface height after the errors removal.
Keywords: interferometric altimetry    baseline error    random error    ocean eddy    theoretical and systematic simulation    error removal strategy    
1 INTRODUCTION

Oceanic eddies are important ubiquitous components of the ocean circulation system and play an essential role in the transportation and distribution of marine materials, energy, heat, and freshwater (Chelton et al., 2011; Chen et al., 2011; Nan et al., 2011; Wang et al., 2012; Sun et al., 2019). Ocean eddies can cause significant elevation differences among sea surface. According to size, oceanic eddies can be classified as either mesoscale or submesoscale eddies (Qiu et al., 2014; Penna and Gaube, 2019; Zhang et al., 2019). Existing high-resolution altimetry satellites have the capability to capture mesoscale oceanic eddies but not smaller eddies in the submesoscale. Therefore, better sampling (broader coverage, higher accuracy, higher spatiotemporal resolution) wide-swath interferometric imaging altimeter has been proposed, e.g., the Surface Water and Ocean Topography (SWOT) (Fu et al., 2012; Fu and Ubelmann, 2014; Xu et al., 2017) and Guanlan (Chen et al., 2019) satellite missions, which plan to launch in 2021 and 2022, respectively. The presence of instrument noise will affect the accuracy of ocean information extraction; therefore, the error analysis of the swath-wide interferometric altimeter plays an essential role before the launch. NASA SWOT workshop proposed an error analysis method based on wavenumber spectrum analysis (Fu et al., 2012; Fu and Ubelmann, 2014), and developed the SWOT Simulator software to directly evaluate the observation error from the perspective of wavenumber spectrum (Gaultier et al., 2016; Qiu et al., 2016). However, there is no complete simulation process of sea surface scattering modeling, radar altimeter echo modeling, imaging processing, and interference processing in the software. To facilitate data processing after the launch of the interferometric altimeter satellite, we shall analyze the ocean errors from the perspective of complete data processing. In this study, we carried out sea surface height extraction and a baseline and random error analysis based on theoretical and systematic simulation of interferometric altimetry. Further, we analyzed the influence of the error on the eddy extraction precision. Due to the limitation of the analysis depth of a single error source and article length, this article only discusses the random errors and baseline errors than other error sources, although they are equally important.

2 INTERFEROMETRIC ALTIMETRY PRINCIPLE

Based on the Interferometric Synthetic Aperture Radar (In-SAR) technique (shown as Fig. 1a), an interferometric altimeter can obtain an observation swath extending tens to hundreds of kilometers with a small incidence angle (usually within 10°) near nadir. The angle measurement error of the baseline roll angle can cause significant errors in the height measurement of the sea surface elevation. Derived from the interferometric principle, the relationship between the baseline roll angle measurement error and the interferometric height error is (Jin et al., 2014; Kong et al., 2017):

    (1)
Fig.1 Schematic diagram of interferometric altimetry (Jin et al., 2014) (a); ocean surface height error induced from baseline roll angle error of 1 arcsec (b)

where, θ is the incidence angle, r is the distance between antenna phase center and the sea surface, is the angle measurement error of the baseline roll angle, ∆α is the angle measurement error of the baseline roll angle.

Figure 1b shows the elevation error caused by the baseline roll angle error increases with the increase of the incidence angle in the cross-track direction. The baseline error of 1 arcsec results in an average height measurement error of 45 cm in the swath, but it can even reach ~80 cm in the distal swath (incidence angle of 10°). The elevation error caused by the baseline roll angle error can be simply recorded as the baseline error.

The random error is the height error caused by the interferometric phase error, and the influencing factors include signal-to-noise ratio (SNR) decorrelation, geometric decorrelation, and angular decorrelation, etc. In this study, ocean signal decoherence is mainly caused by primary and auxiliary SAR imaging error or their registration error. The random error is also related to the incidence angle shown as Eq.2 (Kong et al., 2017).

    (2)

where, k is radar wave number, σphase is interferometric phase error, B is baseline length, θ is the incidence angle.

The random height error induced from the interferometric phase error is shown in Fig. 2. The error is small in the middle and large on both sides of the edge in the interferometric swath, and the root-mean-square error (RMSE) is 0.85 cm at the gird of 3 km×3 km.

Fig.2 Ocean surface height error induced from the random interferometric phase error
3 METHOD

Our method consists of two parts, theoretical simulation and systematic simulation (Fig. 3). The goal of the theoretical simulation is to obtain an arbitrary ideal sea surface with the eddy signals or significant elevation difference. This ideal sea surface can be obtained using a numerical ocean model. Here, we took Regional Ocean Modeling System (ROMS) (Shchepetkin and McWilliams, 2005) sea surface model data as the ideal ocean surface height. The system simulation is the whole process simulation from altimeter observation to data processing, and the simulation process is consistent with the altimeter remote sensing observation and data processing. In this study, the goal of the systematic simulation is to extract the elevation from the ideal sea surface of theoretical simulation and analyze the eddy signal based on the designed parameters of the interferometric altimetry system. The systematic simulation process mainly includes sea surface scattering modeling (Vandemark et al., 2016), radar altimeter echo modeling (Zeng et al., 2010), primary and auxiliary SAR imaging processing (e.g., Range-Doppler algorithm in Cumming and Wong (2005)), and interference processing (Kong et al., 2017). We obtain the baseline error according to the altimeter principle in Eq.1 and add it to the modeling echo in the radar altimeter. Besides, the random phase error is also induced in the interference processing because of the SAR imaging inaccuracy and the registration error of the primary and secondary SAR images. The systematic simulation accuracy of the altimetry system will be analyzed by considering the influence of the random and baseline errors. Finally, the ocean eddies are extracted (Isern-Fontanet et al., 2003; Liu et al., 2016) from theoretical and systematic ocean surfaces and make a comparison in detail.

Fig.3 Flow chart of theoretical and systematic simulation of interferometric system

To analyze the effect of baseline error on ocean signal, a combined-strategy of low-pass filtering, empirical orthogonal function (EOF) decomposition, and linear fitting, is proposed for error removal of the finally systematic sea surface. The strategy is as follows: compared with the sea elevation signal, the random noise is a type of high-frequency signal. Therefore, as far as possible, useful ocean signals can be retained through low-pass filtering. Then, the single sea surface height (SSH) image is divided into several sub-images in along-track direction to form a sub-image sequence. Then, EOF analysis is applied to the sequence to extract the signal mode. The first few valid modal components can be selected using EOF decomposition and modal signal validity analysis (Li et al., 2000; Wenzel and Schröter, 2014). From the retained modal signal, the baseline error can be estimated by linear fitting in the cross-track direction, because of the characteristic of the approximate linear increase of the baseline error. Then, the estimated baseline error can be removed from the retained EOF signal to obtain a "clean?sea surface elevation. The "clean?sea surface elevation can be compared with the originally theoretical simulated surface.

4 RESULT AND DISCUSSION

Through three scenarios, we analyzed the effects of the altimeter system's baseline error and random error in data processing for the ocean surface obtained from the theoretical simulation and systematic simulation, and ocean surface signal extraction accuracy. The designed parameters of the interferometric altimetry system are shown in Table 1.

Table 1 Numerical simulation parameters of interferometric altimetry

Scenario 1: Altimeter system does not contain baseline error, only random error is considered because of imaging error or registration error of primary and secondary images.

Figure 4 shows the intermediate result of the theoretical and systematic simulation of the altimeter system. It can be seen that the random error only makes the sea surface obtained from the systematic simulation rough compared with the theoretical simulation. However, the ocean eddy signal is still clear and distinct.

Fig.4 Random phase error analyses in data processing of altimeter system

Scenario 2: Altimeter system contains baseline error, and the random error in data processing is also taken into account at the same time.

Figure 5 shows the theoretical simulated surface (Fig. 5a), systematic simulated surface (Fig. 5b), and the surface after error removal from Fig. 5b using the proposed combined-strategy (Fig. 5c & d). To estimate the error in the sea surface obtained from the systematic simulation and evaluate the error removal effect of the combined-strategy, we subtracted Fig. 5bd from Fig. 5a to get the sea surface error and the residual error (Fig. 6). Figure 6a shows that the error distribution pattern is consistent with that interferometric principle in Section 2; therefore, the baseline error is the main error source in the sea surface obtained from the systematic simulation. However, it can be seen that Fig. 6b is almost the same as Fig. 6a, and it shows the effect of random error removal using low-pass filtering is relatively weak because of the low proportion of simulated high-frequency noise. Table 2 also confirms the poor remove effect, and residual error decreased by only 0.2 (0.1) cm on average (RMSE). However, the proposed EOF and linear fitting combined-strategy have a noticeable effect for baseline error elimination shown in Fig. 6c. Table 2 shows the final residual error decreases from 2.7 (3.1) cm to 0.1 (1.0) cm on average (RMSE), and the residual error is mainly distributed at image edges in the swath. And, more remarkable, the final residual error conforms to the distribution characteristics of random error in Fig. 2, and the residual error in the margin in the swath is significant, the middle is small. It indicates that the random error is the noise floor in the systematic processing, and it cannot be eliminated by a particularly effective method at present. However, the low-pass filtering method was tried in this study and played a weak error removal effect, and other methods need to be further considered.

Fig.5 Theoretically simulated the sea surface from ocean model (a); systematic simulated sea surface (b); sea surfaces after error removal from (b) using a combined-strategy of low-pass filtering (c), EOF decomposition, and linear fitting (d), individually
Fig.6 System error obtained by Fig. 5b minus Fig. 5a (a); residual error by Fig. 5c minus Fig. 5a (b); final residual error by Fig. 5d minus Fig. 5a (c) (a) and (b) are almost the same. The comparison between (a) and (b) show the baseline error is the main system error and poor removal effect of low-pass filter for random error because of low proportion; (c) shows the noticeable effect of EOF and linear fitting for baseline error, and the final residual error conforms to the distribution characteristics of random error, and it is proved to be the noise floor in the systematic processing.
Table 2 Statistics of baseline error and residual error after the error removal

Scenario 3: Eddy extraction and analysis based on Scenario 2

Ocean eddies can be extracted by various methods (Isern-Fontanet et al., 2003; Mason et al., 2014; Liu et al., 2016). We extracted ocean eddies from the theoretical simulated surface, the systematic simulated surface, and the surface after error removal in Fig. 5, separately. In the theoretical simulated surface, Eddy 1 and Eddy 2 were detected (Fig. 7a). A dark blue area of low sea surface height can be seen separating the two eddies. In Fig. 7b, the existence of the baseline error slightly enlarges the boundary of Eddy 1. In contrast, it dramatically reduces the boundary of Eddy 2, mainly because of the reduced area of low sea surface height values between the two eddies in the right swath. Figure 7c shows that the effect of error removal using low-pass filtering is weak, i.e., the detected eddies are almost the same as Fig. 7b. Figure 7d shows that the height characteristics of the theoretical simulated surface are entirely restored, and the range of the area of low sea surface height is also restored; consequently, the extracted eddies are almost identical to those of the ideal sea surface shown in Fig. 7a. The sea surface height after baseline error removal through the proposed combined-strategy is perfect.

Fig.7 Ocean eddy extraction based on ocean surface height of scenario 2 Eddy 1 and Eddy 2 were detected in the theoretical simulated surface.A dark blue area of low sea surface height can be seen separating the two eddies. Ocean signals are effectively extracted from the sea surface height after baseline error removal.
5 CONCLUSION

A new generation of altimetry satellites provides the possibility for better observation of the sea surface, especially for the submesoscale ocean signals. Before satellite launch, data processing and the related error analysis based on interferometric altimetry are particularly important, and they directly affect the satellite's ability to capture ocean signals. Because ocean eddies can cause significant elevation differences on the sea surface, this is very convenient for the theoretical simulation and systematic simulation for wide-swath interferometric altimetry. Therefore, in this study, the sea surface height extraction, sea surface error analysis, and ocean signal extraction are realized using the theoretical simulation and the systematic simulation. In particular, by analysis of three scenarios of the random error in data processing, baseline error in the altimetry system, and sea surface eddy extraction, the proposed combined-strategy of low-pass filtering, EOF decomposition, and linear fitting achieves a superior effect in systematic baseline error removal, thus significantly improving the ability of satellite altimetry to capture ocean signal. Through simulation experiments, the final error decreases from 2.7 (3.1) cm to 0.1 (1.0) cm on average (RMSE). The residual errors have the same distribution characteristics as random errors, and it is proved the noise floor in the altimeter data processing. It has also been proved that, after error removal, the sea surface can be recovered, and the ocean eddies are accurately extracted.

6 DATA AVAILABILITY STATEMENT

The theoretical sea surface can be simulated by ROMS model or other ocean models. NASA SWOT workshop use ROMS model data for error analysis and developed a SWOT error simulator. A large number of ideal sea surfaces were provided in its doc folder in SWOT error simulator at the website: https://github.com/SWOTsimulator/swotsimulator/tree/master/doc/images.

References
Chelton D B, Schlax M G, Samelson R M. 2011. Global observations of nonlinear mesoscale eddies. Progress in Oceanography, 91(2): 167-216. DOI:10.1016/j.pocean.2011.01.002
Chen G X, Hou Y J, Chu X Q. 2011. Mesoscale eddies in the South China Sea: mean properties, spatiotemporal variability, and impact on thermohaline structure. Journal of Geophysical Research: Oceans, 116(C6): C06018. DOI:10.1029/2010JC006716
Chen G, Tang J W, Zhao C F, Wu S H, Yu F J, Ma C Y, Xu Y S, Chen W B, Zhang Y H, Liu J, Wu L X. 2019. Concept design of the "Guanlan" science mission: China's novel contribution to space oceanography. Frontiers in Marine Science, 6: 194. DOI:10.3389/fmars.2019.00194
Cumming I G, Wong F H. 2005. Digital Processing of Synthetic Aperture Radar Data: algorithms and Implementation. Artech House, Boston. 625p.
Fu L L, Alsdorf D, Morrow R, Rodriguez E, Mognard N. 2012. SWOT: the Surface Water and Ocean Topography Mission: Wide-Swath Altimetric Elevation on Earth. Jet Propulsion Laboratory, National Aeronautics and Space Administration, Pasadena. http://hdl.handle.net/2014/41996.
Fu L L, Ubelmann C. 2014. On the transition from profile altimeter to swath altimeter for observing global ocean surface topography. Journal of Atmospheric and Oceanic Technology, 31(2): 560-568. DOI:10.1175/JTECH-D-13-00109.1
Gaultier L, Ubelmann C, Fu L L. 2016. The challenge of using future SWOT data for oceanic field reconstruction. Journal of Atmospheric and Oceanic Technology, 33(1): 119-126. DOI:10.1175/JTECH-D-15-0160.1
Isern-Fontanet J, García-Ladona E, Font J. 2003. Identification of marine eddies from altimetric maps. Journal of Atmospheric and Oceanic Technology, 20(5): 772-778. DOI:10.1175/1520-0426(2003)20<772:IOMEFA>2.0.CO;2
Jin G W, Xu Q, Zhang H M. 2014. Synthetic Aperture Radar Interferometry. National Defend Industry Press, Beijing. (in Chinese)
Kong W Y, Chong J S, Tan H. 2017. Performance analysis of ocean surface topography altimetry by Ku-Band near-nadir interferometric SAR. Remote Sensing, 9(9): 933. DOI:10.3390/rs9090933
Li X F, Pietrafesa L, Lan S F, Xie L A. 2000. Significance test for Empirical Orthogonal Function (EOF) analysis of meteorological and oceanic data. Chinese Journal of Oceanology and Limnology, 18(1): 10-17. DOI:10.1007/BF02842536
Liu Y J, Chen G, Sun M, Liu S, Tian F L. 2016. A parallel SLA-based algorithm for global mesoscale eddy identification. Journal of Atmospheric and Oceanic Technology, 33(12): 2743-2754. DOI:10.1175/JTECH-D-16-0033.1
Mason E, Pascual A, McWilliams J C. 2014. A new sea surface height-based code for oceanic mesoscale eddy tracking. Journal of Atmospheric and Oceanic Technology, 31(5): 1181-1188. DOI:10.1175/JTECH-D-14-00019.1
Nan F, Xue H J, Chai F, Shi L, Shi M C, Guo P F. 2011. Identification of different types of Kuroshio intrusion into the South China Sea. Ocean Dynamics, 61(9): 1291-1304. DOI:10.1007/s10236-011-0426-3
Penna A D, Gaube P. 2019. Overview of (sub) mesoscale ocean dynamics for the NAAMES field program. Frontiers in Marine Science, 6: 384. DOI:10.3389/fmars.2019.00384
Qiu B, Chen S M, Klein P, Sasaki H, Sasai Y. 2014. Seasonal mesoscale and submesoscale eddy variability along the north pacific subtropical countercurrent. Journal of Physical Oceanography, 44(12): 3079-3098. DOI:10.1175/JPO-D-14-0071.1
Qiu B, Chen S M, Klein P, Ubelmann C, Fu L L, Sasaki H. 2016. Reconstructability of three-dimensional upper-ocean circulation from SWOT sea surface height measurements. Journal of Physical Oceanography, 46(3): 947-963. DOI:10.1175/JPO-D-15-0188.1
Shchepetkin A F, McWilliams J C. 2005. The regional oceanic modeling system (ROMS): a split-explicit, free-surface, topography-following-coordinate oceanic model. Ocean Modelling, 9(4): 347-404. DOI:10.1016/j.ocemod.2004.08.002
Sun B W, Liu C Y, Wang F. 2019. Global meridional eddy heat transport inferred from Argo and altimetry observations. Scientific Reports, 9: 1345. DOI:10.1038/s41598-018-38069-2
Vandemark D, Chapron B, Feng H, Mouche A. 2016. Sea surface reflectivity variation with ocean temperature at Ka-band observed using near-nadir satellite radar data. IEEE Geoscience and Remote Sensing Letters, 13(4): 510-514. DOI:10.1109/LGRS.2016.2520823
Wang X D, Li W, Qi Y Q, Han G J. 2012. Heat, salt and volume transports by eddies in the vicinity of the Luzon Strait. Deep Sea Research Part Ⅰ: Oceanographic Research Papers, 61: 21-33. DOI:10.1016/j.dsr.2011.11.006
Wenzel M, Schröter J. 2014. Global and regional sea level change during the 20th century. Journal of Geophysical Research — Oceans, 119(11): 7493-7508. DOI:10.1002/2014JC009900
Xu Y S, Gao L, Zhang Y H. 2017. New generation altimetry satellite SWOT and its reference to China's swath altimetrysatellite. Remote Sensing Technology and Application, 32(1): 84-94. (in Chinese) DOI:10.11873/j.issn.1004-0323.2017.1.0084
Zeng D Z, Sun H W, Zeng T, Long T. 2010. A high accuracy method for interference fringes suppression in SAR distributed targets' raw data simulation. In: 2010 IEEE International Geoscience and Remote Sensing Symposium. IEEE, Honolulu. p. 4675–4678, https://doi.org/10.1109/IGARSS.2010.5651373.
Zhang Y J, Hu C M, Liu Y G, Weisberg R H, Kourafalou V H. 2019. Submesoscale and mesoscale eddies in the Florida Straits: observations from satellite ocean color measurements. Geophysical Research Letters, 46(22): 13262-13270. DOI:10.1029/2019GL083999