Optimised gravity anomaly fields from along-track multi-mission satellite altimeter over Malaysian seas

Marine gravity anomalies are crucial parameters and elements for determining coastal and ocean geoid, tectonics and crustal structures, as well as offshore studies. This study aims to derive and develop a marine gravity anomaly model over Malaysian seas from multi-mission altimetry data. Universiti Teknologi Malaysia 2020 Mean Sea Surface Model is computed based on along-track data from nine satellite missions, incorporating TOPEX, Jason-1, Jason-2, ERS-2, Geosat Follow on (GFO), Envisat-1, CryoSat-2, SARAL/AltiKa, and Sentinel-3A. The data exploited are from 1993 to 2019 (27 years). Residual gravity anomaly is computed using Gravity Software, and two-dimensional planar Fast Fourier Transformation method is applied. The evaluation, selection, blunder detection, combination, and re-gridding of the altimetry-derived gravity anomalies and Global Geopotential Model data are demonstrated. Cross-validation procedure is employed for data cleaning and quality control using the Kriging interpolation method. Then, cross-validation procedure is applied to the tapering window width 200, which adopting the GECO model denotes the optimum gravity anomaly with root mean square errors in the range of ± 4.2472 mGal to ± 6.0202 mGal. The findings suggest that the estimated marine gravity anomaly is acceptable to be implemented in the marine geoid determination and bathymetry estimation over Malaysian seas. In addition, the results of this study are valuable for geodetic and geophysical applications in marine areas. Along-track altimetry data are used for mean sea surface derivation. Mean sea surface model is utilised in the estimation of marine gravity anomalies. Global Geopotential Model is crucial in the marine gravity estimation of a region. Along-track altimetry data are used for mean sea surface derivation. Mean sea surface model is utilised in the estimation of marine gravity anomalies. Global Geopotential Model is crucial in the marine gravity estimation of a region.


Introduction
Coastal areas are home to nearly half of the world's population, and the impacts of sea level rise (accompanied by tides, storm surges, erosion, and other effects) in the area is a significant concern in the near future (Urban et al. 2018). One of the imperative parameters required in the study of these phenomena is marine gravity anomaly. The parameter is crucial to develop the gravity model of

Open Access
Terrestrial, Atmospheric and Oceanic Sciences Page 2 of 16 Yazid et al. Terrestrial, Atmospheric and Oceanic Sciences (2022) 33:1 the Earth and for research related to global tectonics and continental margin structure. Conventionally, in marine areas, the airborne and shipborne were utilised to measure gravity. However, the information of gravity data are inadequate for the specified study area due to time constraints and the high cost necessitated to conduct the survey. Therefore, a better option is to use multi-mission satellite altimeters for comprehensive data acquisition, particularly for marine gravity estimations. Satellite altimetry is acknowledged as an invaluable source for providing homogeneous and economical data. Moreover, altimetry data also enable continuous, high-accuracy, high resolution, and broad coverage ocean study, making it significant for marine geodesy applications.
According to Liu et al. (2016), satellite altimetry data are utilised to compute and determine marine gravity anomaly, and it is significant in providing abundant marine gravity and seabed geophysical information. Besides, it is crucial in developing high-precision global gravity model for high-precison global gravity model for climate study and marine resources information. Moreover, Tanaka et al. (2019) implemented satellite altimetry measurements and GRACE to examine seismic gravity changes and sea-level changes associated with geoid height variability. Gravity models can be derived from ocean surface height or slope measured by satellite altimetry in the space domain (Rapp 1979;Haxby et al. 1983;Sandwell and Smith 1997;Andersen and Knudsen 1998;Hwang 1998;Fan et al. 2020).
Marine gravity field modelling depends on the accuracy and resolution of the assembled multi-mission satellite altimetry data. It is based on the following factors: (1) the precision of altimeter range; (2) the density of spatial track and along-track sampling rate; (3) the diversity of track orientation; (4) the accuracy of the modelled ocean tide corrections, particularly over coastal areas; and (5) the low-pass filters applied to the profile data .
As a result, regional marine gravity anomaly over Malaysian seas is computed based on the combination of altimetry-derived gravity anomaly and Global Geopotential Model (GGM)-derived gravity anomaly data. Thus, to determine precise gravity anomaly and estimate marine geoid model, several considerations have been taken into account, including the assessment of tapered window width. Computation of residual gravity anomaly has been conducted using the remove-compute-restore technique, planar estimation of Stokes' function, and Fast Fourier Transformation (FFT) method. Subsequently, the residual gravity anomaly is combined with GGMderived gravity anomaly from satellite only (GO_CONS_ GCF_2_DIR_R5) and combined solutions model (GECO) to obtain full spectrum gravity anomaly data. Altimetryderived gravity anomaly has been evaluated and refined based on cross-validation approach to detect and remove outliers from the data. Then, Kriging spatial interpolation has been used to perform the 0.06° integral of the altimetry-derived gravity anomaly data.
Concisely, this study demonstrates the computation and modelling of marine gravity using multi-mission satellite altimeters over Malaysian seas. The findings of this study are significant for marine geoid determination and bathymetry estimation over Malaysian seas. In addition, modelling of geodynamic phenomena like polar motion, Earth rotation, and crustal deformation can be predicted by implementing gravity data into the computations (Bogusz et al. 2015).

Description of the study area
The study area is concentrated over the Malaysian marine areas, including Peninsular Malaysia, Sabah, and Sarawak, which are bordered between latitude 0° to 14° and longitude 95° to 126° (see Fig. 1).

Data used 2.2.1 Along-track altimetry data
The present Universiti Teknologi Malaysia 2020 Mean Sea Surface, known as UTM20 MSS model, developed by Hamden et al. (2021) over Malaysian seas, is applied in this study. The UTM20 MSS model is computed based on data from nine (9) satellite missions comprising TOPEX, Jason-1, Jason-2, ERS-2, Geosat Follow On (GFO), Envisat-1, CryoSat-2, SARAL/AltiKa, and Sentinel-3A, from 1993 to 2019 (27 years). ERS-1 and Geosat-3 missions are not used for this model because they are old geodetic mission that have low precision range (Andersen et al. 2015). TOPEX-class (TOPEX, Jason-1 and Jason-2) are established as a reference to ESA-class (Envisat-1 and Cryosat-2). To improve the spatial resolution of the MSS model, data from Jason-1 Phase C GM and Cryosat-2 are implemented in this study. Then, the development of the MSS model is conducted. It involves the combination of Exact Repeat Mission (ERM) and Geodetic Mission (GM) data.
To provide the highest resolution of the model, three (3) types of mission are considered in the development of the UTM20 MSS model, including Exact Repeat Mission (ERM), Interleaved Mission (IM), and Geodetic Mission (GM). Figure 1 represents along-track Sea Surface Height (SSH) for each satellite mission and along-track SSH for the combination of multi-mission satellite altimetry involved in this study. All satellite missions are adjusted to the TOPEX reference throughout the data processing phase. The multi-mission satellite altimetry data implemented in determining UTM20 MSS are listed in Table 1.
Satellite altimetry data obtained in this study are presented by the Technical University of Delft (Netherlands). The data can be assessed through the Radar Altimeter Database System (RADS) server at Universiti Teknologi Malaysia (UTM), which provides the latest orbital information and geophysical corrections. Satellite altimetry data implemented in the computation have been preprocessed based on the optimal range and geophysical corrections for the Malaysian region. Most of the ranges and geophysical corrections implemented in this study are based on user manuals and gradual experiences from prior studies (Scharroo et al. 2012;Din et al. 2014Din et al. , 2019Hamid et al. 2018;Zulkifle et al. 2019).
In general, MSS model is determined by temporal average method, which depends on the following procedure: data selection and pre-processing, crossover adjustment, ERM mean track derivation, removal of seasonal variability, and data gridding. After geophysical correction and bias elimination are conducted in pre-processing section, crossover adjustment procedures are performed.  Yazid et al. Terrestrial, Atmospheric and Oceanic Sciences (2022) 33:1 Crossover adjustment for dual-satellite missions is performed to adjust the discrepancy between two satellite observations at similar location. This procedure is significant when combining different satellite altimetry data, including ERM and GM. According to Hamid et al. (2018), crossover adjustment is a practical method to reduce errors and enhance the accuracy of the multimission satellite altimetry measurement. This could also minimise the height differences at the crossover between ascending and descending and limit track errors. Moreover, 19-year moving average method is applied in the UTM20 MSS model as recommended by Yuan et al. (2020) to certify that the residual errors of tide models are more degraded on the MSS model. Moving average technique is significant to remove the annual and semiannual variations from altimetry data and the formula is expressed in Smith (2003). For further studies, the integration of gravity anomaly with other marine and land gravity data, marine gravity anomaly data of 15 km from the coastal are excluded area to avoid using low-quality data near the coast. Figure 2 presents the overview of the processing flows in developing the MSS model over the Malaysian seas and the map of the UTM20 MSS model is illustrated in Fig. 3.

Global marine gravity field from DNSC08, DTU10, DTU13, DTU15, DTU17, and Sandwell models
To evaluate the derived gravity anomaly, six (6) global marine gravity fields from DNSC08, DTU10, DTU13, DTU15, DTU17, and Sandwell models from along-track altimetry data with 92,934 gravity points are used in the validation and evaluation processes. Hence, the optimal and relevant derived gravity anomaly are cross-validated to detect and remove blunders.

Satellite-derived gravity data from satellite-only and combined solutions
The importance of GGM in geoid computation has long been recognised, and this has resulted in the continuous development of new GGM. Efforts in developing GGM using satellite data has started since 1970s with the GEM-1 model developed by NASA/Goddard Space Flight Centre. Up until now, there are more than hundreds of GGMs available with free access provided to the scientific community. The International Centre for Global Earth Models (ICGEM) collects all existing GGMs and makes them publicly available on their website (http:// icgem. gfz-potsd am. de/ ICGEM/) (Barthelmes and Köhler 2012), where any interested users can download. The methodology of evaluating the appropriate GGM-derived gravity-related fields to be used for geoid computation in any region (in our case, the geoid for Malaysian seas) is necessary as it is a standard procedure, particularly for geoid computation using remove-restore technique (Karpik et al. 2016). Based on the evaluation with airborne-derived data from the Department of Survey and Mapping Malaysia, the satellite-only and combined solutions from GO_ CONS_GCF_2_DIR_R5 and GECO model represent the optimal and appropriate Global Geopotential Models (GGMs). GO_CONS_GCF_2_DIR_R5 is a GOCE satellite-only model based on a full combination of GOCE-SGG (Satellite Gravity Gradiometer) and GOCE-SST (Satellite-to-Satellite Tracking) that also encompasses GRACE (Gravity Recovery and Climatic Experiment) and LAGEOS (Laser GEOdynamics Satellite) data (Bruinsma et al. 2013). Besides, this model was attained by the direct approach with maximum degree/order, 300 of the harmonic expansion. However, GECO model was developed by Polytechnic University of Milan in 2015 up to a maximum degree/ order 2190 of the harmonic expansion. For the highfrequency gravity signal, GOCE-TIM-5R was used in conjunction with EGM2008, which provides precision enhancement at low and middle frequencies (Gilardoni et al. 2016). Thus, these models are applied and serve as the main parameters in deriving gravity anomaly from multi-mission satellite altimeters.

Residual sea surface height derivation
This section describes the methodology utilised and the data processing strategies applied in this study to compute marine gravity anomaly over Malaysian seas using altimetry data. The overall flow chart for computing marine gravity anomaly using the 2D FFT method is shown in Fig. 4.
Residual sea surface height, h residual , can be considered as residual geoid height representing the medium wavelength of the geoid height signal, which is the subsidiary computed quantity from altimetry measurements. It is crucial in gravity anomaly computation. According to Salam (2005), residual sea surface height is computed from the difference of mean sea surface and GGMderived geoid height, where geoid heights, N, approximately corresponds to the mean sea surface heights derived from satellite altimeter. Fig. 3 The map of the UTM20 mean sea surface model (Hamden et al., 2021) Page 6 of 16 Yazid et al. Terrestrial, Atmospheric and Oceanic Sciences (2022)  The residual sea surface height is computed by removing geoid height quantities from the mean sea surface. The process of removing long-wavelength geoid height is performed using the global geoid height from satellite-only and combined GGM solutions; GO_CONS_ GCF_2_DIR_R5 and GECO model, involving the remove-compute-restore method. This method has also been applied by Nguyen et al. (2020) in their calculation and determination of marine gravity using altimetry data. Hence, the approach is applied and assessed in the derivation of residual sea surface height and can be expressed as follows:

South China Sea Malacca Straits
where, h residual is defined as the residual sea surface height, MSS is the mean surface height derived from along-track altimetry data, N GGM is described as the GGM-derived geoid height (from GO_CONS_GCF_2_ DIR_R5 and GECO model). This method is used to assess and evaluate the accuracy of residual sea surface height. Subsequently, the residual sea surface height data is utilised to derive marine gravity anomalies over Malaysian seas. (1) Removing geoid field based on GGM provides residual field that can be identified as reference field, which is statistically more harmonised and smoother than the entire field. The effect of removing the reference is that the gravity information outside the data area is indirectly accounted for, and it will result in low correlation distance of covariance functions (Andersen 2013).

Determination of marine gravity anomalies using the 2D FFT method
Based on previous studies described by several researchers, it is highly recommended to exploit an abundance of altimetry data to estimate gravity anomalies. Discussions pertinent to this topic are presented as follows.
For instance, Andritsanos (2000) and Vergos (2002) mentioned that employing an abundance of data in gravity field estimation indicates better performance in the frequency domain to expedite the computation process. The frequency domain, also known as spectral approach, presents a time-efficient approach of processing abundance datasets and yielding similar results for numerical integration (Knudsen 1993;Tziavos et al. 1997Tziavos et al. , 1998. It is now recognised that altimetry data can be utilised to estimate marine gravity anomaly. In this study, the altimetry data employed is the residual sea surface Fig. 4 Marine gravity computation using along-track altimetry data Yazid et al. Terrestrial, Atmospheric and Oceanic Sciences (2022) 33:1 height, h residual . The h residual derived using along-track altimetry data is used to estimate gravity anomaly, which is then used as input data in marine geoid computation utilising the remove-compute-restore technique. The residual gravity anomaly is derived using the GRAVSOFT Fourier domain programs in either planar or sphere implementing discrete Fast Fourier Transform algorithm to estimate the numerous integrals of physical geodesy. The continuous two-dimensional Fourier transformation can be explained as (Forsberg and Tscherning 2008): where g and G represent gravity anomaly data. However, g is in the space domain and G is in the spectral domain, respectively. Then, k denotes the wavenumbers. Basically, FFT data and parameter must be in periodic and represented on a finite grid interval, and the continuous transform integral procedure is estimated in each direction by the fundamental discrete Fourier transform, as follows: Continuous estimation with discrete Fourier Transforms provides an increase to the difficulty of numbers, such as periodic effects Conventionally, these difficulties can be prevented by creating a data window, such as implementing a cosine taper to the data, where data near to the edge of the grid is multiplied by a function w decaying from 1 to 0 as a cosines curve. Otherwise, the implementation of zero-padding normally good in physical geodesy. Hence, zero-padding is applied for the major Fourier programs like GEOFOUR programs (GEOFOUR program is implemented in the computation of residual gravity anomaly using altimetry data). Normally, the zero-padding is organised by specifying the number of points along the grid margin where the computed data should efficiently close to zero. Figure 5 illustrates the zero-padding of grid procedure.
The task to obtain gravity anomaly in this study has been accomplished using GRAVSOFT software (Forsberg and Tscherning 2008). Equations (2) until (5) are implemented in GRAVSOFT software based on 2D Fast Fourier Transformation, which is required for computing residual gravity anomalies. The assessment and evaluation of the computed residual sea surface height, h residual , are implemented in these equations. Previously, GRAVSOFT software was for regional and local gravity field modelling, for example, geoid determination, computations of the vertical deflections, and enhancement of gravity anomalies from satellite altimetry. This software was established in the 1970s at the Geodetic Institute, later on its descendant, Geodetic and Geophysical Institute, University of Copenhagen, together with the National Survey and Cadastre of Denmark. Moreover, the GRAVSOFT programs performs fundamental operations of physical geodesy, fundamental arithmetic, and data files operations in point or grids formats. One of the modules in GRAVSOFT named Geofour, which is a routine program that can be employed for the computation of gravity anomaly using the FFT technique. The contributions of GGM is to provide low frequencies gravity field spectrum, whereby the contribution of GGM to the gravity field are determined depending on the coefficients computed in the spherical harmonic expansion configuration of the Earth's disturbing potential (Vergos 2002). However, the altimetry-derived gravity anomaly contributes to the medium frequencies during the gravity anomaly computation.
The implementation of Geofour gravity anomaly computation requires tapered window width (TWW) to be applied accordingly. This approach is performed to obtain accurate residual gravity anomaly results. TWW is defined as the width of the cosine-tapered window zone shown in grid points. According to Pilz et al. (2012), the mathematical principles entail an infinitely long time series to accomplish Fourier transform; hence, such windowing will initiate the Fourier transform method to acquire non-zero values, particularly at lower frequencies (commonly called spectral leakage, i.e., some frequencies manage to leak to other frequencies). Multiplying the data windows by a 'taper' is a regular exercise before performing the Fourier transform method. The taper comprises an operation that efficiently decays the residual gravity data to zero, closing the ends of each window. It is designed to minimise the consequences of discontinuation between the beginning and the last part of the time series. Even though spectral leakage is inevitable, it can be reduced by transforming the shape of the taper function in a way to minimise strong discontinuation close to the edges of the window (Pilz et al. 2012).
Mathematically, the cosine taper is specified with respect to time, t, and taper ratio, a. The cosine window denotes an effort to set the data efficiently to zero at the borders while not significantly decreasing the level of the windows transformed. The tapering will indicate the next consequence of decreasing in spectral power leakage from the spectral peak to the frequencies by far (Pilz et al. 2012). Based on Table 2, the TWW setting is selected randomly based on the data distribution. Generally, it is complex to present the recommendations with tapers for use in any particular situation (Bingham et al. 1967). Thus, in this study, the tapering window width is selected and stopped until it reaches 200 as the taper starts to decay to zero value and closes the end of each window. Two parameters are changed in the computations: tapering window width and removing the mean value from the input anomalies. The tapering window is utilised to eliminate the periodicity effect, which is related to the discrete Fourier transform. According to Forsberg and Tscherning (2008), this can be avoided by zero-padding or windowing the data on the borders (edges).

Restoring the long wavelength gravity anomaly
The reference gravity model used in this study is the GGM model from satellite-only and combined solutions model from GO_CONS_GCF_2_DIR_R5 and GECO model. Hence, the full free-air gravity anomaly obtained is referred as Δg full . The residual gravity anomaly, Δg altimeter , computed by Geofour is combined with the GGM-derived gravity anomaly, Δg GGM , as shown below:

Cross-validation procedure
In statistical terms, the cross-validation approach can be interpreted as the separation of data samples into sub-samples, so that the analysis can be primarily accomplished on a single sub-sample. Although, further sub-samples are kept "blind" for further use in the verification process of preliminary analysis. The cross-validation theory was initiated by Geisser and Eddy (1979).
The cross-validation approach estimates the value of g at data point by applying the values from the neighbouring data points and neglecting the value at the point in question (Kiamehr 2010). All data in the database can be evaluated by employing this procedure. On the other hand, points that are inadequately estimated by neighbouring data may be a symptom of anomalous values. Therefore, the cross-validation technique is applied in this study sequentially to validate the computed gravity anomaly data before determining marine geoid height.
In this study, the derived gravity anomaly from altimetry data that representing the lowest root mean square error (RMSE) are selected to perform a cross-validation procedure. Then, the cross-validation procedure is conducted repetitively on estimated gravity anomalies (Δg reduced ) in the database. There are three (3) techniques that have been followed: i. Eliminate the known point form dataset (the predicted point) ii. The remaining datasets (surrounding points) are employed to estimate the value of the previously removed point via interpolation method. iii. The error is computed by comparing the predicted and observed values at a similar point. The difference between the estimated gravity anomaly (interpolated value), Δg pre , and the estimated gravity anomaly (known value), Δg estimated , provide the information of interpolated residual for gravity anomaly, Δg residual There are numerous inaugurated gridding techniques, for instance, Kriging (Krige 1951), least-squares collocation (Moritz 1972), Bjerhammar method (Bjerhammar (6) g full = g altimeter + g GGM Table 2 The tapering window width for gravity anomaly computations  Subsequent to the cross-validation and gridding process, one histogram is constructed to illustrate the absolute values of differences between the estimated and the original reduced gravity anomaly (Δg reduced − Δg estimated ). Where the abrupt change of slope is clearly illustrated in the residual values below; thus, the expected value can be specified. The example of the residual histogram results before the cross-validation process of residual altimetryderived gravity anomaly is illustrated in Fig. 6. Based on the histogram in Fig. 6, the residual of altimetry-derived gravity anomaly plummet at 5mGal. Thus, any residual values greater than ± 5mGal are considered as outliers and removed. The analysis and validation of the gravity data are performed based on a cross-validation scheme (Kiamehr 2010).

Residual marine gravity anomalies from ten (10) examined tapering window widths
Gravity data play an important role in the determination of local marine geoids. In this study, the residual gravity anomalies are computed by applying the residual sea surface height results from Eq. 1. The residual gravity anomalies are derived using the Geofour program from GRAVSOFT software. Geofour program are developed to compute and determine the gravity field with the Fast Fourier Transformation (FFT) planar approach. Ten (10) tapering window widths have been used to compute residual gravity anomaly, as illustrated in Fig. 7. Where Fig. 7a-j represent the distributions of residual gravity anomaly from tapering window widths of 5, 10, 20, 50, 70, 100, 125, 150, 170 and 200, respectively. Tapering window widths involve multiplying the data windows by a taper before operating the Fourier Transform. Taper performs as a function to decay each window to zero near the end. A small block of tapering window width presents the narrow rectangular window. This Fig. 6 The residual histogram before the cross-validation process of residual altimetry-derived gravity anomaly Page 10 of 16 Yazid et al. Terrestrial, Atmospheric and Oceanic Sciences (2022) 33:1 indicates that the Fourier transform becomes wider; hence, more leakage is accomplished. Referring to Fig. 7, as the tapering window width increases, the distributions of residual gravity anomalies approach zero value. The decay process of the derived-residual gravity anomaly starts at tapering window width with block 50. Hence, at the end of tapering window width with block 200, the residual gravity anomaly almost completely decayed to zero near the  16 Yazid et al. Terrestrial, Atmospheric and Oceanic Sciences (2022) 33:1 end of each window, and the rectangular window becomes wider. There have been other assessments for altimetryderived residual gravity anomaly at tapering window width 200 with the implementation of GGM from EGM96, EGM2008, GO_CONS_GCF_2_DIRR5, and GECO, as shown in Fig. 8. In Fig. 8, the distributions of the residual gravity anomaly map for EGM2008 and GECO are almost similar, while the distributions of the residual gravity anomaly map for GO_CONS_ GCF_2_DIR_R5 is almost identical to the residual gravity anomaly map for EGM96. The reason for this is that degree and order plays a significant role in modelling. EGM2008 model and GECO model published high-degree global geopotential models, 2190°, while GO_CONS_GCF_2_DIR_R5 published low-degree global geopotential model, 300°, and 360° degree for EGM96 model. EGM2008 and GECO have relatively higher spatial resolutions. These models are derived by combining satellite gravimetric data, terrestrial and marine gravity data depending on spherical harmonic functions with full expansion degree and order of 2190° or 2159° (Wu et al., 2021). Furthermore, the residual altimetry-derived gravity anomaly over coastal areas are widely different compared to the open sea, as seen in Fig. 8.

The evaluation of full spectrum gravity anomalies with global marine gravity anomalies models
The GGM-derived gravity anomaly from GO_CONS_ GCF_2_DIR_R5 and GECO model are combined with the residual altimetry-derived gravity anomaly (starting from tapering window width with block 5 to 200) to provide the full information of gravity anomaly. Then, these altimetry-derived gravity anomalies data are evaluated and compared with the global marine gravity anomalies models from DNSC08, DTU10, DTU13, DTU15, DTU17, and Sandwell V29.1.The statistical analysis of the altimetry-derived gravity anomaly evaluation is depicted in Table 3. As illustrated by statistical analysis in Table 3, it indicates that the RMSE and standard deviation values decrease with increasing tapering window width (starting tapering window width 5 until 200). This result is related to the results discussed in Sect. 4.1, which is the residual altimetry-derived gravity anomaly are decayed to zero near the end of each window with the increasing block.
Moreover, RMSE values based on the implementation of GECO model from tapering window width 200 presents the lowest values in the ranges of ± 4.3317 mGal to ± 6.0726 mGal after verified with global marine gravity anomaly from DNSC08, DTU10, DTU13, DTU15, DTU17, and Sandwell models. However, there is a significant difference in RMSE value from the GO_ CONS_GCF_2_DIR_R5 model for tapering window width 200 with respect to their evaluation and verification with global marine gravity anomaly model that yielding RMSE value in the range of ± ± 18.5430 mGal to ± 19.2715 mGal. Hence, based on this results, it can also prove the statement from Yazid et al. (2016), who mentioned that the combined GGMs solutions provide better fit to terrestrial gravity data than satellite-only models, is accurate. Besides, these combinations lead towards the best approximation of the Earth's gravitational field. By combining all the gravity data, some limitations on higher degree expansion can be diminished. However, errors in terrestrial data remain.

Cross-validation results for marine gravity anomalies from tapering window width 200
As discussed in Sect. 4.2, the altimetry-derived gravity anomaly using tapering window width 200 presents the lowest RMSE, indicating good accuracy among other derivations. Hence, these derivations are selected to perform the cross-validation procedure since it is significant to detect and remove blunders in the altimetry-derived gravity anomaly data. Table 4 illustrates the statistical analysis for altimetry-derived gravity anomaly before and after the cross-validation procedure. Based on the analysis, the average differences, standard deviation, and RMSE values are decreased by approximately ± 0.0524 mGal to 0.0845 mGal. After performing the cross-validation procedure. Due to this result, altimetry-derived gravity anomaly with tapering window width 200 is selected as the altimetry-derived gravity anomaly model over Malaysian sea with RMSE values of ± 4.2472 mGal to ± 6.0202 mGal with respect to the evaluation and verification with global marine gravity anomaly from DNSC08, DTU10, DTU13, DTU15, DTU17, and Sandwell.
Moreover, the gravity anomalies in this paper have marked differences with respect to the anomalies from Scripps Institution of Oceanography (SIO) V29.1 about ± 6 mGal. In reference to Sandwell et al. (2014), this model was developed from 70 months of free-air marine gravity anomalies calculated from CryoSat-2 data and Jason-1 GM data, which are augmented by older altimetry calculations from GeoSat and ERS-1 with a 31-month operating period. Hence, these altimetry data determine more short-wavelength gravity features. However, as compared to DNSC08, DTU10, DTU13, DTU15 and DTU17 model published by Denmark Technical University (DTU), they are based on multi-mission satellite altimetry and include up to 10 different satellites with an operating period of 12 to 20 years. Boergens et al. (2018) state that single-mission altimetry data are spatially and temporally limited. Additionally, only missions with a short-repeat orbit such as Envisat, Jason-2 or SARAL, provide time series of sea level variation directly. As a trade-off, long or non-repeat orbit missions such as CryoSat-2 provide a very dense spatial resolution, but their repeat time is insufficient for extracting time series. As a result, multi-mission altimetric data allows for improved spatial and temporal resolution. This is the reason gravity anomalies relative to SIO V29.1 provide high RMSE.The map of the final product is illustrated in Fig. 9.

Conclusion
Determination of marine gravity anomaly with high resolution and high accuracy is vital for various implementations. Hence, the main focus of this study is to determine and optimise the altimetry-derived gravity anomaly over Malaysian seas. The 2D FFT method has been implemented to estimate and compute marine gravity anomalies with 0.06° grid resolution for the Malaysian seas using along-track altimeter mean sea surface data. The satellite-only and combined solutions model GGM from GO_CONS_GCF_2_DIR_R5 and GECO are implemented to derive residual sea surface height. Then, the combined full-spectrum gravity anomaly are utilised using the remove-compute-restore technique.
Hence, after the evaluation and verification procedure with the global marine gravity anomaly models, the analyses found that the altimetry-derived gravity anomaly presents RMSE value of ± 4.3317 mGal to ± 6.0726 mGal.
The bold values indicate the lowest values of the RMSE between the altimetry-derived gravity anomalies compared with the global marine gravity anomalies models from DNSC08, DTU10, DTU13, DTU15, DTU17, and Sandwell V29.1  Page 14 of 16 Yazid et al. Terrestrial, Atmospheric and Oceanic Sciences (2022) 33:1 After performing the cross-validation procedure, the RMSE value decreased by approximately ± 0.0524 mGal to 0.0845 mGal. Thus, the final output provides the altimetry-derived gravity anomaly model over Malaysian seas with good accuracy based on RMSE values of ± 4.2472 mGal to ± 6.0202 mGal.
To compute and develop the altimetry-derived gravity anomaly model with good quality and accuracy, evaluation of the optimal and appropriate GGM model with terrestrial data, such as airborne-derived gravity data, is necessary to align them with regional terrestrial gravity data. Therefore, in the future, along-track altimetry data can be implemented directly in the process of estimating marine gravity anomaly using vertical deflection technique without concerning about the accuracy of the GGM in the region.