Annual Sea Level Amplitude Analysis over the North Pacific Ocean Coast by Ensemble Empirical Mode Decomposition Method

: Understanding spatial and temporal changes of seasonal sea level cycles is important because of direct inﬂuence on coastal systems. The annual sea level cycle is substantially larger than semi-annual cycle in most parts of the ocean. Ensemble empirical mode decomposition (EEMD) method has been widely used to study tidal component, long-term sea level rise, and decadal sea level variation. In this work, EEMD is used to analyze the observed monthly sea level anomalies and detect annual cycle characteristics. Considering that the variations of the annual sea level variation in the Northeast Paciﬁc Ocean are poorly studied, the trend and characteristics of annual sea level amplitudes and related mechanisms in the North Paciﬁc Ocean are investigated using long-term tide gauge records covering 1950–2016. The average annual amplitude of coastal sea level exhibits interannual-to-decadal variability within the range of 14–220 mm. The largest value of ~174 mm is observed in the west coast of South China Sea. In the other coastal regions of North Paciﬁc Ocean, the mean annual amplitude is relatively low between 77 and 124 mm for the western coast and 84 and 87 mm for the eastern coast. The estimated trend values for annual sea level amplitudes in the western coastal areas of South China Sea and Northeast Paciﬁc Ocean have statistically decreased over 1952–2014 with a range of − 0.77 mm · yr − 1 to − 0.11 mm · yr − 1 . Our results suggested that the decreasing annual amplitude in the west coast of South China Sea is in good agreement with the annual mean wind stress associated with the Paciﬁc Decadal Oscillation (PDO). This wind phenomenon also explains the temporal variations of annual sea level amplitude in Northeast Paciﬁc Ocean, especially the high correlations since 1980 (R = 0.61 − 0.72).


Introduction
Seasonal cycle is a common signal in sea level change and fluctuates due to various physical forces, including atmospheric pressure, wind stress, precipitation, river runoff, ice melting, ocean current, and steric component [1]. Thus, the temporal sea level change in seasonality (annual and semi-annual) is not constant with time. Understanding seasonal sea level cycles and its spatial and temporal changes is important because its temporal changes (amplitude or phase) may directly affect the coastal systems [2]. An increase in the amplitude of seasonal sea level cycle during storm surge season also increases the risk of coastal flooding and coastal erosion. Thus, accurately estimating the seasonal amplitude and phase in sea level cycle and understanding the physical mechanisms of the seasonal sea level cycle are helpful in assessing the impact of future climate change on coastal ocean environments. Early studies of seasonal sea level cycle are based on tide gauge records [3], which are limited in coastal areas. Since the early 1990s, satellite altimetry has provided near-global and high-quality observations of sea surface heights. Accuracy of satellite altimetric measurements is higher than 3 cm in open ocean, but rapidly decreases in marginal seas and when approaching coastal area [4,5]. This technique examines seasonal sea level cycle with high spatial coverage and resolution in coastal and open oceans. Chen et al. [6] used satellite altimetry and steric component to analyze seasonal sea level cycle on a global scale. Their results indicated that thermal effect has a dominant impact on observed seasonal variability, especially for zonal averages over regions. Vinogradov et al. [7] investigated seasonal sea level cycle in a near-global ocean by using satellite altimetry and an ocean circulation model. Their results indicated that annual cycle is substantially larger than semi-annual cycle in most parts of the ocean, except in tropical regions. They suggested that surface wind stress and heat flux are the main driving factors for seasonal sea level variations in tropical and mid-latitude regions, respectively, and these two mechanisms play a role at high-latitude regions. Vinogradov and Ponte [8] found that the spatial variability of annual sea level cycle varies from the shoreline to the adjacent shallow and deep ocean areas. The annual amplitude from tide gauge records is critically greater than that of nearby altimetry possibly due to the observations being conducted at different locations, inadequacy of some corrections, or a smoothing of the gridded altimetric data [9]. Wahl et al. [2] reported that the annual amplitude of sea level cycle has gradually increased along the United States Gulf coast since 1990s. They suggested that sea surface temperature and sea level pressure are the reason for the amplitude increase. Feng et al. [10] indicated that wind forcing and steric component are dominant impacts on the marginal seas and deep ocean (or open ocean) of the northwest Pacific Ocean, respectively. Amiruddin et al. [11] also reported similar results in the South China Sea. The temporal variation of seasonal sea level cycle has been investigated at regional scales, such as by Wahl et al. [2] for the United States Gulf coast, Feng et al. [10] and Amiruddin et al. [11] for the northwest Pacific coast, Torres and Tsimplis [12] for the Caribbean Sea, and Passaro et al. [13] and Cheng et al. [14] for the Baltic Sea and North Sea, respectively. Only a few studies have been conducted on seasonal sea level cycle in the Northeast Pacific Ocean. Hence, the driving factors behind temporal and spatial variations remain unrevealed. The amplitude of seasonal sea level cycle exhibits decadal variability, such as in the Caribbean Sea [12], the United States Gulf coast [2], the South China Sea [11], and the northwest Pacific [10]. Related studies are based on traditional linear regression analysis (TLRA) method with 5-year overlapping periods (shifts month-by-month or year-by-year). This technique assumes that the amplitude and frequency of annual and semi-annual sea level cycles are constant within each time segment; hence, these parameters cannot reflect instantaneous changes and may be contaminated by other frequency signals when the window size is not appropriate [10]. The present study aims to detect changes in annual sea level cycle by using a developed ensemble empirical mode decomposition (EEMD) method [15] that could supplement the TLRA. Wu et al. [15] proposed a definition of a modulated annual cycle that allows the annual cycle to change from year to year by using a developed EEMD method [16]. This technique has been used for analyzing annual cycle in other geophysical data, such as surface air temperature and sea surface temperature [15,17].
Tide gauge stations are sparsely distributed in the early 20th century but rapidly increased in the 1950s [18]. In order to compute more accurate sea level trends and minimize the impact of low-frequency variability, all gauges needed at least 50 years of data [19,20]. EEMD is introduced to extract annual cycle data from observed sea level anomalies to analyze the trends and characteristics of annual sea level amplitude derived from long-term tide gauge records along the North Pacific Ocean coast covering 1950-2016. The dominant force is regionally identified to understand physical mechanisms responsible for sea level variability at the annual frequency. The rest of the paper is organized as follows. Data sets and methodology are described in Sections 2 and 3, respectively. Results and discussion are reported in Sections 4 and 5, and Section 6 offers conclusions.

Tide Gauge Records
Monthly mean sea level records at 316 tide gauge stations in the North Pacific Ocean were obtained from the Permanent Service for Mean Sea Level database [21,22]. Monthly records for Revised Local Reference were selected as the primary data set. If the data do not have substantial and unknown datum shifts, then some metric data were included [23]. All records provided at least 50 years of data between 1950 and 2016. The following tide gauge stations were excluded to avoid the influence of data gaps and large nonlinear signals due to typhoons, earthquakes, or groundwater pumping on the EEMD decomposition results [24]: those with data completeness of less than 90% within their time span, and those with records including nonlinear vertical land motion compared with neighboring tide gauge and altimeter measurements. Tide gauge data can be influenced by inverted barometer (IB) effect and contaminated by anomalous jumps owing to external forces and environmental changes [25]. In this study, IB effect and datum shift in tide gauge records were corrected according to tide gauge processing, introduced by Lan et al. [25]. In addition, individual tide gauge record with strong local effects was excluded to avoid the influence of local effect on the final outcome.
Monthly mean values of sea level time series from 107 tide gauge stations along North Pacific coast regions fulfilled the above criteria and are used in this study. Figure 1 shows the locations and record lengths of 107 tide gauge stations in the North Pacific Ocean. These tide gauges are classified into eight subregions (Table 1), namely the west coast of the South China Sea (SCSW), the East China Sea (ECS), the Sea of Japan (SJ), the northeast coasts of Japan (JNE), the southeast coasts of Japan (JSE), the Gulf of Alaska (GA), the west coast of Canada (WC), and the west coast of America (WA) (Figure 1). for sea level variability at the annual frequency. The rest of the paper is organized as follows. Data sets and methodology are described in Sections 2 and 3, respectively. Results and discussion are reported in Sections 4 and 5, and Section 6 offers conclusions.

Tide Gauge Records
Monthly mean sea level records at 316 tide gauge stations in the North Pacific Ocean were obtained from the Permanent Service for Mean Sea Level database [21,22]. Monthly records for Revised Local Reference were selected as the primary data set. If the data do not have substantial and unknown datum shifts, then some metric data were included [23]. All records provided at least 50 years of data between 1950 and 2016. The following tide gauge stations were excluded to avoid the influence of data gaps and large nonlinear signals due to typhoons, earthquakes, or groundwater pumping on the EEMD decomposition results [24]: those with data completeness of less than 90% within their time span, and those with records including nonlinear vertical land motion compared with neighboring tide gauge and altimeter measurements. Tide gauge data can be influenced by inverted barometer (IB) effect and contaminated by anomalous jumps owing to external forces and environmental changes [25]. In this study, IB effect and datum shift in tide gauge records were corrected according to tide gauge processing, introduced by Lan et al. [25]. In addition, individual tide gauge record with strong local effects was excluded to avoid the influence of local effect on the final outcome.
Monthly mean values of sea level time series from 107 tide gauge stations along North Pacific coast regions fulfilled the above criteria and are used in this study. Figure 1 shows the locations and record lengths of 107 tide gauge stations in the North Pacific Ocean. These tide gauges are classified into eight subregions (Table 1)

Satellite Altimetric Data
Satellite altimetry has provided high-quality, global observations of sea surface height since early 1990s [4,5]. The altimeter data in this study are monthly gridded sea level anomaly maps at 0.25 • × 0.25 • spatial resolution with respect to a 20-year (1993-2012) mean sea level, which are provided by the Archiving, Validation, and Interpretation of Satellite Oceanographic (AVISO) data center (http://www.aviso.altimetry.fr/en/home.html, accessed on 19 January 2021), which are sea level anomaly products merged from several satellite altimetry missions of TOPEX/Poseidon, Jason-1/2, ERS-1/ERS-2, ENVISAT, GFO, Cryosat-2, Saral/AltiKa, and HY-2A. All standard corrections recommended by AVISO, including instrumental, medium (ionosphere, wet, and dry troposphere), and geophysical (solid earth, polar, ocean and load tides, sea state bias, and IB) corrections, were applied to the altimeter data. The data cover the period from January 1993 to December 2016. Altimetric observations at the closest grid points to the tide gauge location were compared with tide gauge data for the analysis of annual sea level cycle. The average distance between tide gauges and altimetric observations is approximately 10 km.

In Situ Hydrographical Data
Steric sea level variation was calculated from an updated version (v7.3) of in situ measurements of ocean temperature and salinity data processed by Ishii et al. [27] (called Ishii data hereinafter). This product has passed through several quality controls and provided monthly gridded fields with a spacing resolution of 1 • × 1 • and 28 vertical levels from the surface to 3000 m depth covering the period of 1955-2019. Steric sea level variation represents sea level change due to seawater density variation caused by fluctuations in ocean temperature and salinity. Its changes over water depth, h steric , can be computed from seawater density change [28] expressed as: where S 0 and T 0 are the mean salinity and temperature, respectively; S and T are the monthly salinity and temperature of sea water, respectively; ρ and ρ 0 are the density of water and the reference density of water, respectively; and z is the depth.

Wind Data
Surface wind stress was calculated using zonal and meridional (U and V) component winds obtained from the National Centers for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR) Reanalysis Monthly Means database (downloaded from https://www.psl.noaa.gov/, accessed on 19 January 2021) [29] by using bulk Remote Sens. 2021, 13, 730 7 of 18 formulas [30]. This parameter provides monthly gridded fields from 1948 to present time with a spacing resolution of 2.5 • × 2.5 • . The IB effect was also removed from tide gauge records by using mean sea level pressure from NCEP/NCAR Reanalysis data set [29].

Climate Indices
Two natural modes of interannual/decadal variability, namely El Niño-Southern Oscillation (ENSO) and Pacific Decadal Oscillation (PDO), were analyzed to investigate the temporal variation of annual sea level amplitude in the North Pacific Ocean and the effect of low-frequency climate variability. These phenomena are typically considered the interannual and decadal variability with 2-8 and 10-30 years, respectively [31,32]. Their climate variabilities can be represented by the multivariate ENSO index (MEI) and PDO indices defined as the first unrotated principal component of six main observed variables (including sea-level pressure, zonal and meridional components of the surface wind, sea surface temperature, surface air temperature, and total cloudiness fraction of the sky) over the tropical Pacific (data source at: https://www.esrl.noaa.gov/psd/data/correlation/ mei.data, accessed on 19 January 2021) and the leading principal component of North Pacific (poleward of 20 • N in 1900-1993) monthly sea surface temperature anomalies (data source at: http://research.jisao.washington.edu/pdo/PDO.latest, accessed on 19 January 2021) [33,34], respectively. The two indices covered different time spans, that is, MEI for 1950 to 2016 and PDO for 1948 to 2016. For a clear description of the interannual and decadal characteristics of these two indicators, the MEI index was used to remove the MEI index with 65-month low-pass filter from low-passed MEI with 5-month filter, and the PDO index was smoothed with a 65-month running mean [35].

Methods
The analysis method is based on the Empirical Mode Decomposition (EMD), which is part of the Hilbert Huang Transform [36]. This method analyzes nonlinear and nonstationary time series, investigates the spatial and temporal characteristics of various data timescales [37,38], and it has been widely used to examine annual cycle in other geophysical data, such as surface air temperature and sea surface temperature [15,17]. The goal of EMD is to decompose a time series, x(t), into a set of components called intrinsic mode function (IMF), C i (t), which can be expressed as: where n is the total number of IMFs, and n is close to but no more than log 2 N, where N is the length of the input data, X(t) [39]. r n is the residual. Each IMF may represent different physical processes of various time scales from the highest frequency (IMF 1 ) to the lowest frequency (IMF n ), and the sum of all modes and residual can be used to reconstruct original data [40]. The decomposition process of EMD may produce mode mixing (one IMF component includes two or more different frequencies), so the ensemble EMD (EEMD) is applied for the considerable reduction of the mode mixing disparity [41], which can be expressed as: where X E is the added noise data, X is original data, noise is the added white noise, andεis the amplitude of the added white noise. According to Wu and Huang [41], εis basically set to 0.2 standard deviation of the original data. However, when the data are dominated by high-or low-frequency signals, the noise amplitude may be small or large, respectively. These two steps were performed repeatedly by adding white noise to the original data in Equation (3) and decomposing the added noise data into IMFs in Equation (2). This procedure was repeated M times with a different added white noise. M times were averaged for the corresponding IMFs of the decompositions to obtain the (ensemble) means of IMFs as the final values. The number of ensemble members, M, is approximately 100 times [41].
For the monthly tide gauge data, the second and third components by EEMD method are always the annual cycle signals but still include other long timescale variation. A single EMD must be applied to the combination of the second and third components to extract a relatively narrow band annual cycle [15]. The first IMF mode of the EMD is designated as the modulated annual cycle of monthly sea level data. The first and last 2 years of the decomposition results were excluded because the results of EEMD may produce end effects [17]. This analysis method was applied for the sea level data from Tajiri tide gauge records in the SJ coast as shown in Figure 2. The annual cycle results show a good performance and reflect the instantaneous changes of modulated annual sea level cycle from year to year. Figure 2 shows other low-frequency components at interannual-todecadal (IMF 4 to 5) and multi-decadal (IMF 6 to 8) time scales. Although this study focuses on the annual cycle (i.e., the second and third EEMD modes), the analysis also provides a tool to extract the different physical processes of various time scales (e.g., References [37,42]). The annual amplitude of sea level cycle derived from altimetry or tidal gauges can be extracted by calculating the half distance between the upper and lower envelopes of the annual cycle. Furthermore, the criterion that the correlation coefficient of annual amplitude of sea level between the individual tide gauge record and the regional average of individual anomalies in each subregion must be larger than 0.4 (generally considered correlation above 0.4 to be moderate to relatively strong [43]) was imposed to exclude individual tide gauge records with strong local effect. In North Pacific coast regions, 107 tide gauges fulfilled this criterion (five stations are removed), 80% have correlations > 0.6, and the average correlation coefficient of all 107 tide stations was 0.71. Figure 1  For the monthly tide gauge data, the second and third components by EEMD method are always the annual cycle signals but still include other long timescale variation. A single EMD must be applied to the combination of the second and third components to extract a relatively narrow band annual cycle [15]. The first IMF mode of the EMD is designated as the modulated annual cycle of monthly sea level data. The first and last 2 years of the decomposition results were excluded because the results of EEMD may produce end effects [17]. This analysis method was applied for the sea level data from Tajiri tide gauge records in the SJ coast as shown in Figure 2. The annual cycle results show a good performance and reflect the instantaneous changes of modulated annual sea level cycle from year to year. Figure 2 shows other low-frequency components at interannual-to-decadal (IMF 4 to 5) and multi-decadal (IMF 6 to 8) time scales. Although this study focuses on the annual cycle (i.e., the second and third EEMD modes), the analysis also provides a tool to extract the different physical processes of various time scales (e.g., References [37,42]). The annual amplitude of sea level cycle derived from altimetry or tidal gauges can be extracted by calculating the half distance between the upper and lower envelopes of the annual cycle. Furthermore, the criterion that the correlation coefficient of annual amplitude of sea level between the individual tide gauge record and the regional average of individual anomalies in each subregion must be larger than 0.4 (generally considered correlation above 0.4 to be moderate to relatively strong [43]) was imposed to exclude individual tide gauge records with strong local effect. In North Pacific coast regions, 107 tide gauges fulfilled this criterion (five stations are removed), 80% have correlations > 0.6, and the average correlation coefficient of all 107 tide stations was 0.71. Figure 1 displays the locations and record lengths of 107 tide gauge stations along the North Pacific Ocean coast.

Temporal Variability of Mean Annual Sea Level Cycle
The temporal variability of the annual sea level cycle from observed data was determined by using the developed EEMD method. Figure 3 shows that regional mean annual Remote Sens. 2021, 13, 730 9 of 18 amplitude was obtained by averaging all annual amplitudes derived from tide gauges or satellite altimetric data at the closest grid points to each tide gauge location in each subregion. In the study, we aimed to calculate reliable estimates of regional averages, so the monthly regional averages were determined from at least three individual data in a month. The mean value estimated by EEMD method varies with time and location. The range between maximum and minimum annual amplitudes in eight subregions varies from 14 mm to 220 mm. The annual amplitude of sea level in eight subregions exhibits interannual-to-decadal variability (Figure 3).  TLRA was applied to estimate the annual cycle (amplitude and phase) from sea level data (red and cyan lines in Figure 3). The mean annual amplitude of sea level from tide gauges in each subregion also shows decadal variability, but the interannual variability is relatively insignificant. For example, the mean annual amplitude estimated by EEMD method in the ECS displays a peak between 1962 and 1963 and reaches the maximum of 170 mm, but the result of TLRA is not evident (Figure 3). Kagoshima II tide gauge data with the highest correlation coefficient of 0.9 with regional mean annual amplitudes in ECS are displayed in Figure 4. The annual amplitude in sea level anomaly is approximately 100 mm in the period 1960-1961, which has remarkably increased to 150-200 mm in 1962-1963 but clearly decreased after 1963 (black line in Figure 4). This finding is consistent with the regional mean annual amplitude of wind stress in the ECS (orange line in Figure 4).
For the eight subregions, the temporal variability of annual sea level amplitude extracted by EEMD covering period of 1952-2014 are consistent with that of sea level estimated by TLRA using a 3-year-long data segment with the correlation coefficient ranging from 0.77 to 0.85. Although the decadal variations of the annual sea level amplitudes by the EEMD are remarkably consistent with the TLRA results, the latter may not clearly reflect the details of the interannual variation of annual sea level amplitude, especially when estimating annual amplitudes with long window size (e.g., 5-year), as shown in Figure 3.   The average annual amplitudes in eight subregions obtained by the developed EEMD are shown in Table 2. The largest annual amplitude of 174 mm over 1952-2011 was observed in the SCSW. Ko Sichang tide station located in Gulf of Thailand had the largest mean annual amplitude of more than 200 mm. In other regions, the average annual amplitude was between 77 and 124 mm for those in the western coast and 84 and 87 mm for those in the eastern coast (Table 2). Along the South China Sea coast, seasonal sea level variations are strongly influenced by the local monsoon system [44], which is strong monsoon blowing from the northeast to southwest in winter (December-February), and vice-versa in summer (June-August). In addition, the South China Sea has a special topography and is the largest semi-enclosed marginal sea in the western Pacific (Figure 1) [11]. Thus, the coastal sea level height along SCSW could reach the maximum in winter and minimum in summer; when the topography factor is considered, the annual amplitude will increase. Vinogradov and Ponte [8] indicated that annual sea level cycle mostly has large amplitudes in shallow coastal areas. Steric component, another important factor for annual sea level variability [1], has a relatively small contribution (~42 mm) to the annual sea level cycle for SCSW coast with less than quarter of the mean amplitude of annual sea level cycle ( Table 2). Note that a similar result can also be drawn, if the EN4 [45] is used other than the Ishii, since the annual amplitude is approximately 30 mm covering the same period. In the eight subregions, the contribution of the annual amplitude of steric component in the SCSW to the total annual sea level amplitude is the smallest ( Table 2). This finding may be associated with the smoothing of gridded ocean temperature and salinity data, the shallow continental shelf areas along coast, or the poor spatial sampling of Argo floats in the South China Sea [46]. A good agreement was found between the EEMD and TLRA estimated results by tide gauges ( Table 2). The mean annual amplitude differences between the results of EEMD and TLRA with 3-year-long data segment is 0 ± 2.6 mm in the eight subregions. Moreover, the regional average amplitudes of annual sea level cycle derived from the tide gauge in the eight subregions are larger than those from altimeter data with an average annual amplitude difference of 21 mm. For the same period, the regional average amplitudes in the eight subregions are also larger than those from altimeter data with an average annual amplitude difference of 19 mm. This finding is similar to previous studies [8]. These differences in annual amplitude could partly be attributed to the smoothing of the gridded AVISO data [9]. The regional average annual amplitude from altimeter data along the east coast of North Pacific Ocean is substantially lower by at least 23 mm than that for west coast ( Table 2). Further analysis is needed to understand the mechanisms responsible for the larger annual amplitude in the continental shelf of Northwest Pacific Ocean than in the Northeast coast as obtained from the altimeter data.

Trend of Annual Amplitude of Sea Level Cycle
Trends of regional average amplitude of annual sea level in the eight subregions of the North Pacific Ocean derived from long-term tide gauge data in 1952-2014 are presented in Table 3. Statistically significant negative trends were observed in the most areas of marginal seas, with rates of −0.77 mm·yr −1 to −0.11 mm·yr −1 . For the rest of coastal areas (ECS, SJ, and JSE), the trends were insignificant (Table 3). For each station, the estimated long-term trends of annual sea level amplitude exhibited geographical non-uniformity (Table 1 and Figure 5). The East Asia tide gauge in the areas of ECS and Japan (SJ, JNE, and JSE) coast showed an inconsistent trend in annual sea level amplitude with a range of −5.3 mm·yr −1 to 6.5 mm·yr −1 ; the maximum trend of 6.5 mm·yr −1 was obtained from the Shirahama tide gauge data in the JSE coast ( Figure 5). In the areas of ECS and Japan coast, the annual sea level variability was strongly influenced by various factors, including atmospheric pressure, wind stress, ocean current, and steric component [10]. Therefore, annual amplitude of sea level cycle in these areas still requires continuous monitoring and analysis to understand the extent of climate change effects on these regions and for disaster prevention. Table 3. Trends of mean annual sea level amplitude derived from tide gauges in the eight subregions of the North Pacific Ocean. C.C. refers to correlation coefficient of annual amplitudes between tide gauges and satellite altimetry.   Figure 6 shows the annual cycle amplitude in SCSW that is derived from tide gauges and wind stress and overlaid with the low-pass filtered MEI and PDO. The mean annual sea level cycle amplitude in SCSW decreased significantly from 1952 to 1998, with the smallest amplitude of 115 mm found in 1986-1987. No trend change was observed after 1998. The annual amplitude variability between wind stress and sea level in the SCSW is consistent. The mean annual amplitude of wind stress decreased significantly from 1950 to 1998 but slightly increased after 1998, with a correlation coefficient of +0.62 (Table 4). Cheng et al. [44] indicated that annual sea level variability in the SCSW is dominated by wind forces. The PDO index has an opposite phase to annual amplitude variabilities of sea level and wind stress in the SCSW, which increased significantly from 1950 to 1997, showed the largest index in 1987 and exhibited a slightly decreasing trend after 1997. Comparison between the mean annual amplitudes of sea level and wind stress and PDO index implied that the temporal variations of annual sea level and wind stress amplitudes are negatively correlated with the PDO with the correlation coefficients of −0.48 and −0.46 (Table 4). These results suggested that the changes of annual wind stress in SCSW are associated with PDO [47]. In 1952-2011, the average values of annual sea level amplitude in SCSW decreased significantly. The increases in PDO index prompted the negative trend of the annual amplitude of local monsoon system. In the SCSW, the regional average trend of annual sea level amplitude was significantly decreased with a rate of −0.77±0.04 mm·yr −1 for 1952-2011 (Table 3). A significant decreasing trend was also observed at each station (Table 1 and Figure 5) with a range of −1.15 mm·yr −1 to −0.18 mm·yr −1 . In addition, the average mean rate of annual sea level cycle for three subregions of the Northeast Pacific Ocean showed statistically significant negative trends with a range of −0.34 mm·yr −1 to −0.11 mm·yr −1 ( Table 3). The changes of annual sea level amplitudes for most tide gauge stations in the Northeast Pacific Ocean showed decreasing trends, except for the positive rate at the Cordova tide gauge in the GA with rate of 0.34±0.07 mm·yr −1 (Table 1 and Figure 5). This discrepancy may be caused by the observations covering different time periods. The time period of Cordova station covers 1964-2016, but the developed EEMD method removes the first and last 2 years from the tide gauge records. Almost no observational data were obtained in 1968-1969, resulting in no observational data before 1970. In the GA, the maximum average annual amplitude (~92 mm) was found in 1952-1969 (Figure 3f), leading to the positive trend at Cordova station. In the following section, we focus on the impact of several factors including steric component, wind stress, and climate indices on the long-term temporal variations of annual sea level amplitude in the areas of SCSW, GA, WC, and WA.

West Coast of South China Sea
C.C. between tide gauge and altimetry is statistically significant (at 95% confidence level) in each subregion. Figure 6 shows the annual cycle amplitude in SCSW that is derived from tide gauges and wind stress and overlaid with the low-pass filtered MEI and PDO. The mean annual sea level cycle amplitude in SCSW decreased significantly from 1952 to 1998, with the smallest amplitude of 115 mm found in 1986-1987. No trend change was observed after 1998. The annual amplitude variability between wind stress and sea level in the SCSW is consistent. The mean annual amplitude of wind stress decreased significantly from 1950 to 1998 but slightly increased after 1998, with a correlation coefficient of +0.62 (Table 4). Cheng et al. [44] indicated that annual sea level variability in the SCSW is dominated by wind forces. The PDO index has an opposite phase to annual amplitude variabilities of sea level and wind stress in the SCSW, which increased significantly from 1950 to 1997, showed the largest index in 1987 and exhibited a slightly decreasing trend after 1997. Comparison between the mean annual amplitudes of sea level and wind stress and PDO index implied that the temporal variations of annual sea level and wind stress amplitudes are negatively correlated with the PDO with the correlation coefficients of −0.48 and −0.46 (Table 4). These results suggested that the changes of annual wind stress in SCSW are associated with PDO [47]. In 1952-2011, the average values of annual sea level amplitude in SCSW decreased significantly. The increases in PDO index prompted the negative trend of the annual amplitude of local monsoon system.

Eastern Boundary of the North Pacific Ocean
The annual sea level amplitudes of the Northeast Pacific coasts exhibited considerable interannual-to-decadal variability (Figure 7). Those derived from tide gauges along the Northeast Pacific coast (GA, WC, and WA) were correlated with its steric component derived from Ishii data [27] with the correlation coefficients of 0.25-0.51 (Table 4). However, all annual amplitudes of steric component were less than half of the mean annual sea level amplitude derived from tide gauges (Figure 7, Table 2). These differences are partly attributed to the smoothing of the gridded temperature and salinity product and the observations being conducted at different locations. The regional averages of the annual amplitudes of wind stress in the Northeast Pacific coast were highly correlated with the corresponding quantity for the annual sea level amplitude [48] with correlation coefficients of 0.58 for GA, 0.42 for WC, and 0.56 for WA. Since 1980, the two annual amplitudes in the Northeast Pacific coasts have become consistent (Figure 7) with correlation coefficients of 0.61-0.72 compared with those in the period of 1952-1980 (0.25-0.57). This finding indicates that the influence of wind stress field on the annual sea level amplitude varies at different time periods, and the wind field system has a dominant effect on the annual sea level amplitudes in the Northeast Pacific coasts since 1980. Similar results are found for the comparison with altimeter data, with high correlations in the most areas of north Pacific marginal seas [49]. In addition, the annual wind stress and sea level amplitude from 1980 to 2014 showed similar trends in the Northeast Pacific coasts (Figure 7). The weakened annual sea level amplitude of WC and WA for the past 30 years is speculated to be caused by the annual amplitude of wind force manifesting a negative trend.
Wang et al. [50] links the coastal sea level changes to the PDO in the Northeast Pacific Ocean. Comparison between the mean annual sea level amplitude and PDO index revealed that the temporal variations of annual amplitudes in the Northeast Pacific coasts are lowly correlated with the PDO with correlation coefficients ranging from 0 to 0.26 (Table 4). This finding can be attribute to the PDO impact being negligible for the entire period. Hence, the annual sea level amplitude in the Northeast Pacific coasts is not dominated by the PDO. In addition, the variability of annual wind stress amplitude is low correlated with the PDO with correlation coefficients of 0.15 for GA, 0.31 for WC, and 0.31 for WA. This result indicates that PDO did not have a dominant influence on annual amplitude variability during 1952-2014. However, in the Northeast Pacific Coasts, PDO has crucial effects on the changes in other sea level signals (such as sea level trend [51]). (Table 4). This finding can be attribute to the PDO impact being negligible for the entire period. Hence, the annual sea level amplitude in the Northeast Pacific coasts is not dominated by the PDO. In addition, the variability of annual wind stress amplitude is low correlated with the PDO with correlation coefficients of 0.15 for GA, 0.31 for WC, and 0.31 for WA. This result indicates that PDO did not have a dominant influence on annual amplitude variability during 1952-2014. However, in the Northeast Pacific Coasts, PDO has crucial effects on the changes in other sea level signals (such as sea level trend [51]).

Conclusions
The temporal variations of the annual sea level cycle along the North Pacific Ocean coast were investigated by the developed EEMD method using 107 tide gauge records covering the period of 1950-2016. The results showed that average annual sea level am-

Conclusions
The temporal variations of the annual sea level cycle along the North Pacific Ocean coast were investigated by the developed EEMD method using 107 tide gauge records covering the period of 1950-2016. The results showed that average annual sea level amplitudes exhibit interannual-to-decadal variability in the eight subregions and range between 14 and 220 mm. Comparison of the average temporal variations of the annual sea level amplitudes in each subregion for 1952-2014 shows good agreement (R = 0.77-0.85) between the values computed in this analysis and using the TLRA method (with 3-yearlong data segment), especially at decadal time scales. The largest mean annual amplitude of 174 mm is found in the SCSW. For the rest of coastal regions, the average annual amplitude is between 77 and 124 mm in the western coast and 84 and 87 mm in the eastern coast. These results agree with the annual amplitude derived by TLRA. The mean difference of annual amplitude between the EEMD and TLRA results is 0 ± 2.6 mm in the eight subregions. Therefore, the EEMD method used in this study can sensitively reflect the annual sea level cycle to obtain a good physical understanding of processes involved in the variation of annual sea level cycle.
The areas of ECS and Japan coast (SJ, JNE, and JSE) are characterized by the estimated long-term trends of annual sea level amplitudes in 1952-2014. An inconsistency in spatial patterns was observed with a range of −5.3 mm·yr −1 to 6.5 mm·yr −1 . In these regions, the annual sea level variability is strongly influenced by various factors, including atmospheric pressure, wind stress, ocean current, and steric component [10]. Therefore, their annual amplitudes of sea level cycles still require continuous observation and analysis to understand the extent of climate change effects on these regions.
In the SCSW, the regional average trend of annual sea level amplitude has significantly decreased (−0.77 ± 0.04 mm·yr −1 ) over 1952-2011. For each station located in the SCSW, the annual sea level amplitudes present statistically decreasing trends with a range of −1.15 mm·yr −1 to −0.18 mm·yr −1 . Steric component has relatively minor impacts on the annual sea level cycle in the SCSW and is estimated to contribute~42 mm to sea level. This value is less than quarter of the average annual sea level amplitude, which is the smallest contribution of annual amplitude of steric component to the total annual sea level amplitude in the coastal region of the North Pacific Ocean. This difference may be caused by the smoothing of gridded ocean temperature and salinity data, the shallow continental shelf areas, or the poor spatial sampling of Argo floats in the South China Sea [46]. The temporal variation of the annual sea level amplitude is mainly associated with wind stress [44,52], with a correlation coefficient of +0.62. In addition, the PDO index has opposite phases to annual sea level and wind stress amplitudes for the SCSW with the correlation coefficients of −0.48 and −0.46, respectively. We speculated that the weakened local wind system may be associated with the increase in PDO index [44]. The mean annual sea level amplitude from SCSW gradually decreased as a consequence of weakened monsoon system over the past 60 years. The last PDO shift from cold to warm phase occurred in 2014 ( Figure 6) [50], leading to high and low sea levels in the Northeast and Northwest Pacific, respectively, in the next few decades [25,53,54].
For the Northeast Pacific Ocean coast (GA, WC, and WA), the temporal variations of average annual sea level amplitude show statistically significant negative trends over 1952-2014 with a range of −0.34 mm·yr −1 to −0.11 mm·yr −1 . Most of the tide gauge stations in the Northeast Pacific Ocean exhibited statistically decreasing trends for their annual sea level amplitudes (Table 1 and Figure 5). The regional mean annual sea level amplitudes derived from tide gauge records are in good agreement with the corresponding wind stress (R = 0.42-0.58), especially with its high correlation with the annual amplitude curves since 1980 (Figure 7) with the correlation coefficients of 0.61-0.72. In addition, annual wind stress and sea level amplitude show similar trends in the Northeast Pacific coasts from 1980 to 2014 (Figure 7). The wind force affecting the annual sea level cycle varies at different time periods and has played a dominant role in the Northeast Pacific coasts since 1980. Furthermore, the annual sea level cycle in the Northeast Pacific Ocean associated with PDO was explored using the tide gauges over 1952-2014. The result indicated that the annual sea level amplitude in the Northeast Pacific coasts is not dominated by the PDO.