Selecting a representative optical deep-water area is crucial for accurate satellite-derived bathymetry (SDB) based on semi-theoretical and semi-empirical models. This study proposed a deep-water area selection method where potential areas were identified by integrating remote sensing imagery with existing global bathymetric data. Specifically, the effects of sun glint correction for deep-water areas on SDB estimation were investigated. The results indicated that the computed SDB had significant instabilities when different optical deep-water areas without sun glint correction were used for model training. In comparison, when sun glint correction was applied, the SDB results from different deep-water areas had greater consistency. We generated bathymetric maps for the Langhua Reef in the South China Sea and Buck Island near the U.S. Virgin Islands using Sentinel-2 multispectral images and 70% of the Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) bathymetry data. Additionally, 30% of the ICESat-2 bathymetry data and NOAA NGS Topo-bathy Lidar data served as the validation data to evaluate the qualities of the computed SDB, respectively. The results showed that the average quality of the SDB significantly improved with sun glint correction application by a magnitude of 0.60 m in terms of the root mean square error (RMSE) for two study areas. Moreover, an evaluation of the SDB data computed from different deep-water areas showed more consistent results, with RMSEs of approximately 0.4 and 1.4 m over the Langhua Reef and Buck Island, respectively. These values were consistently below 9% of the maximum depth. In addition, the effects of the optical image selection on SDB inversion were investigated, and the SDB calculated from the images over different time periods demonstrated similar results after applying sun glint correction. The results showed that this approach for optical deep-water area selection and correction could be used for improving the SDB, particularly in challenging scenarios, thereby enhancing the accuracy and robustness of SDB.

This paper discusses the sea level determination using raw Intermediate Frequency (IF) data transmitted from the Global Navigation Satellite System (GNSS) and received by spaceborne GNSS-Reflectometry satellites, TechDemoSat-1. The reflected signals scattered from a sea ice surface and a rough sea surface are investigated. The altimetry method based on the bistatic group delay (code phase) from GNSS signals for sea level estimation are introduced. The two raw IF datasets recorded on January 18 and 27, 2015 for a duration of 40 seconds are analyzed to drain more information than Level 1 data. The results show a good consistence with mean sea surface (MSS) model. The orbit error of the GNSS-R satellite is corrected by a proposed method that combines the MSS and the least squares solution, which help evaluate the actual altimetry precision. The defect of fixed temporal resolution and fixed 4 onboard processing channels of Level 1 data can be improved by postprocessing using Software Define Receiver (SDR) to mine more information, so as to explore the potential. At the end, fake high signal-to-noise ratio (SNR) DDMs from the raw data are analyzed, which provides a reference for the altimetry of GNSS-R technique using raw data.

A multilayer approach is set up for local gravity field recovery within the framework of multi-resolution representation, where the gravity field is parameterized as the superposition of multiple layers of Poisson wavelets located at different depths beneath the Earth's surface. The layers are designed to recover gravity signals at different scales, where the shallow and deep layers mainly capture the short- and long-wavelength signals, respectively. The depths of these layers are linked to the locations of different anomaly sources beneath the Earth's surface, which are estimated by wavelet decomposition and power spectrum analysis. For testing the performance of this approach, a gravimetric quasi-geoid model over the North Sea, QGNSea V1.0, is modeled and validated against independent control data. The results show that the multilayer approach fits the gravity data better than the traditional single-layer approach, particularly in regions with topographical variation. An Akaike information criterion (AIC) test shows that the multilayer model obtains a smaller AIC value and achieves a better balance between the goodness of fit of data and the simplicity of the model. Further, an evaluation using independent GPS/leveling data tests the ability of regional models computed from different approaches towards realistic extrapolation, which shows that the accuracies of the QGNSea V1.0 derived from the multilayer approach are better by 0.4, 0.9, and 1.1 cm in the Netherlands, Belgium, and parts of Germany, respectively, than that using the single-layer approach. Further validation with existing models shows that QGNSea V1.0 is superior with respect to performance and may be beneficial for studying ocean circulation between the North Sea and its neighboring waters.

The determination of high‐quality mean dynamic topography (MDT) over island areas faces multiple challenges, due to the error in a global geopotential model (GGM) as well as the contaminated waveforms of satellite altimetry data and scarcity of surveyed gravimetric observations. We study the potential to enhance regional quasi‐geoid/MDT over island areas by using recently released satellite altimetric gravity data and airborne gravimetric observations. Numerical experiments over the Paracel Islands in the northwestern South China Sea show that the combination of airborne gravity data enhances the local gravity field, and reduces the errors in the associated quasi‐geoid, in comparison to the result directly derived from a GGM. The results demonstrate that the direct application of a GGM may cause significant oscillations in MDT determination over island areas, by a magnitude exceeding 6 cm. The incorporation of altimetric gravity data slightly reduces these errors; however, the bubble‐like errors are prominent, suggesting that altimetry data alone may be insufficient for high‐accuracy MDT recovery. In contrast, combining airborne gravity data dramatically reduces these oscillations, by a magnitude of several centimeters, indicating that airborne gravimetric surveys may be necessary for MDT determination with centimeter‐level accuracy over coastal/island areas. Further comparison with independent ocean data shows that the geostrophic velocities computed by merging a high‐degree GGM and local airborne gravity data can reconstruct real surface ocean circulation; and the GOCE‐based (Gravity Field and Steady‐State Ocean Circulation Explorer) GGMs may be preferable as the reference models in MDT determination.

Monthly gravitational field solutions as spherical harmonic coefficients produced by the GRACE satellite mission require post-processing to reduce the effects of shortwave-length noises and north–south stripe errors. However, the spatial smoothing and de-striping filter commonly used in the post-processing step will either reduce spatial resolution or remove short-wavelength features of geophysical signals, mainly at high latitudes. Here, by using prior covariance information that reflects the spatial and temporal features of the geophysical signals and the correlated errors derived from the synthetic model, together with the covariance matrix of the formal errors for the monthly gravity spherical harmonic coefficients, we apply the Kalman filter to separate the geophysical signal from GRACE Level-2 data and simultaneously to estimate the correlated errors. By increasing the number of observations, the iterative process is applied to update the state vector and covariance in the Kalman filter because the prior information is not accurate. Due to the inevitable truncation error, multiple gridded-gain factors method considering different temporal frequencies has been developed to recover the geophysical signal. The results show that the Kalman filter can reduce the high-frequency noises and correlated errors remarkably. When compared with the commonly used filter, no spatial filter (such as Gaussian filter) is used in the Kalman filter. Therefore, the estimated signal preserves its natural resolution, and more detailed information is retained. It shows good consistency when compared with mascon solutions in both secular trend and annual amplitude.