An Analysis of the 8.85- and 4.42-Year Cycles in the Gulf of Maine

: In the background of global warming and climate change, nuisance ﬂooding is only caused by astronomical tides, which could be modulated by the nodal cycle. Therefore, much attention should be paid to the variation in the amplitude of the nodal cycle. In this paper, we utilize the enhanced harmonic analysis method and the independent point scheme to obtain the time-dependent amplitudes of the 8.85-year cycle of N 2 tide and the 4.42-year cycle of 2N 2 tide based on water level records of four tide gauges in the Gulf of Maine. Results indicate that the long-term trends of N 2 and 2N 2 tides vary spatially, which may be affected by the sea-level rise, coastal defenses, and other possible climate-related mechanisms. The comparison between Halifax and Eastport reveals that the topography greatly inﬂuences the amplitudes of those cycles. Moreover, a quasi 20-year oscillation is obvious in the 8.85-year cycle of N 2 tide. This oscillation probably relates to a 20-year mode in the North Atlantic Ocean.


Introduction
The number of extreme sea-level events is expected to grow due to the influence of global warming and climate change. This poses a more damaging threat to coastal areas, especially for the low elevation coastal zones, where up to 310 million people reside [1] and are accompanied by many infrastructures. If adaptation measures are not implemented, 0.2% of the global population and 0.3% of the global domestic product are going to be flooded and lost annually in 2100 at the least [2]. An extreme sea-level indicator for the contiguous United States coastline has been recommended, which is comprised of separate indicators for mean sea level and storm surge climatology [3].
Extreme sea-level events are considered as a compounding effect of three factors: mean sea level (MSL), storm surge, and astronomical tides, if the effect of surface gravity wave is ignored [4]. Much focus has been paid to the change of MSL, and it is noted that MSLs have been rising during the 20th century [5] with multi-decadal fluctuations [6]. Recently, Frederikse et al. [7] reached an agreement between the rise of global mean sea level in the 20th century and its underlying causes, and offered a great explanation to the multidecadal variability of global MSL during the 1940s and 1970s. Storm surges are induced by tropical cyclones or winter storms in coastal zones, which are the direct response to wind stress and perturbation of atmospheric pressures on the sea surface [8].
What should be noted is that storm surges are often modulated through the nonlinear interaction between astronomical tides [9,10]. Astronomical tides are the rise and fall of sea levels influenced by the mutual position of the Moon, the Sun, and the rotating Earth, into yearly windows, and the independent tidal analysis was performed by the T_TIDE for each window [28]. In this process, the calculated amplitude was set to NaN if the missing data exceeded 25% in one year. The total numbers of abandoned months at the four stations are listed in Table 1. Ignoring the long-period and shallow water constituents, 7 main tidal constituents (M 2 , S 2 , K 1 , O 1 , N 2 , P 1 , 2N 2 ) were achieved by the T_TIDE. The average amplitudes of the main tidal constituents for the timespan at the four tide gauges are listed in Table 2. It can be found that the M 2 , N 2 , and S 2 tides are the three major constituents in the Gulf of Maine. Amplitudes at Eastport are larger than the other three gauges for any main tidal constituent. Eastport and Halifax are located at similar latitudes, but they show a remarkable difference in the amplitude of the tidal constituents. This shows the effect of the semi-enclosed topography on the tide, and the amplitudes of M 2 , N 2 , and 2N 2 tides increase by 4 times from Halifax to Eastport. In this work, the spatial and temporal variations in N 2 and 2N 2 tidal amplitudes were analyzed. N2, and 2N2 tides increase by 4 times from Halifax to Eastport. In this work, the spatial and temporal variations in N2 and 2N2 tidal amplitudes were analyzed.

Methodology
Eliot pointed out that conventional harmonic and extreme analysis always obscures 18.61-, 8.85-, and 4.42-year cycles, and dedicated techniques are required for their identification [12]. Therefore, the enhanced harmonic analysis (EHA) method is applied in this paper for analyzing the water level time series of N 2 and 2N 2 tides. This method is proposed and used to investigate the temporal variations in internal tides in the South China Sea [29]. Then, it was conducted by the S_TIDE MATLAB toolbox [30], which is developed from the well-known T_TIDE. Besides, this method was also used to explore the tidal-fluvial interaction in the Columbia River Estuary [30] and the 18.61-year cycle of M 2 tide in the Gulf of Maine [27]. Here, we applied it to analyze the variations in the 8.85-and the 4.42-year cycles in the Gulf of Maine and the Bay of Fundy.
Traditionally, the cycle with the fixed period is estimated by a least-squares method [23], and the equation can be expressed as: where P(t) is the estimated value of the tidal amplitude or phase at the time t (unit: year), and A 0 and A 1 are a constant and the linear trend, respectively. a and b are the amplitudes of the cosine and sine functions of the cycle. c 0 represents the period of the cycle in years. Subsequently, the amplitude (H) and phase (G) of the cycle can be obtained as: However, it is assumed that the amplitude (H) and phase (G) are time-varying and can be estimated by the EHA. Therefore, Equation (1) is modified for the N 2 tide as: Similarly, the equation for the 2N 2 tide is: The a c 0 (t) and b c 0 (t) (c 0 = 18.61, 8.85, and 4.42, respectively) can be obtained by: where a c 0 ,i and b c 0 ,i are values at the i-th independent point for the c 0 cycle, m is the number of the independent points (IPs), and w c 0 ,t,i is the weighted coefficient for the i-th independent point at the time t for the c 0 cycle. Values of a c 0 and b c 0 at other points can be obtained by interpolation between the IPs. Equations (6) and (7) represent the independent point scheme, which has been examined successfully in many previous works [27,31,32]. Equations (2) and (3) are modified as: To guarantee stabilization, convergence, and smoothness in the interpolation process, the cubic spline interpolation is implemented in this work [30,31]. The calculation method of the weighted coefficient refers to Equations (20) and (21) in Zong et al. [32].
Finally, A 0 , A 1 , a c 0 (t), and b c 0 (t) can be obtained from many years' observations. The detailed matrix form of the solution refers to Equations (6)-(8) in Pan et al. [27].

Determine the Numbers of IPs
Different IP numbers represent oscillations at different time scales in the EHA method [30]. If the number of IPs is set equal to 1, the oscillation is a constant. If it is set equal to 0, the oscillation represents the long-term trend of the data. Moreover, the long-term trend consists of linear and nonlinear trends, which are decided by a parameter in the S_TIDE MATLAB toolbox. When the number of IPs is equal to 0, the result is the linear trend of data if the parameter is equal to 2. Additionally, if the parameter is equal to 3, the oscillation represents the nonlinear trend. It can be found that the S_TIDE is a powerful tool to analyze time-series data. Furthermore, it should be realized that the number of IPs plays a vital role in the results of the EHA. Either fewer or more IPs will affect the outcome of the EHA method. Small IPs will introduce large root-mean-square error (RMSE) and large IPs will cause the overfitting and fictitious conclusion. As a result, a primary and vital procedure is to develop a reasonable evaluation criterion on the number of IPs.
In this paper, three parameters are used to determine the number of IPs. The first is the RMSE between the hindcasts obtained by the EHA and the primary data. It can be expressed as follows: where P i represents the estimated value of the tidal amplitude or phase by the EHA and P i represents the real tidal amplitude or phase. N is the length of P i . The second is the 95% confidence interval. The shorter confidence interval can give rise to a higher accuracy result based on the same confidence level. Therefore, the length of the 95% confidence interval can lead to a restriction for the result of EHA.
The third is the signal-to-noise ratio (SNR). It can be calculated as follows: where H represents the time-varying amplitude H(t), and H int represents the 95% confidence interval for H(t). It is credible when the SNR is greater than 2 [33]. A series of sensitivity experiments are implemented to select the optimum number of IPs. The N 2 tide is affected by 18.61-and 8.85-year cycles while 2N 2 tides are mainly influenced by 8.85-and 4.42-year cycles. These two tides can be analyzed for these two cycles by the EHA method simultaneously. Then, the process of determining the numbers of IPs can be revealed distinctly by following analysis taking the 8.85-year cycle of N 2 tide at Boston as an example.
Tables 3 and 4 display the SNR and RMSE for the 8.85-year cycle of N 2 tide at Boston, respectively. The IP1 and IP2 represent the numbers of IPs for the 18.61-and the 8.85-year cycle of N 2 tide, respectively. The numbers of IP1 and IP2 are greater than 1 because the amplitudes of those cycles (18.61 and 8.85) change with time. In other words, the amplitudes are constant in the case of IPs equal to 1, which turns into the traditional method (Equation (1)) and are not the objective of this paper (Equation (4)). However, it is noteworthy that the number of IPs is set equal to 1 for the 8.85-year cycle when the EHA method is used to analyze the 2N 2 tide because the major influence is introduced by the 4.42-year cycle.
The values of SNR and RMSE both reduce gradually with the increase of IP2 when the number of IP1 is fixed. The traditional standard that the SNR is greater than 2 is deemed to be credible. However, a more stringent criterion is applied in this paper: the last number of IP2 can keep the reduced tendency of SNR and the corresponding SNR value must be greater than 2.
According to the above stringent criterion, the applicable numbers of IP1 and IP2 are selected in Table 3 (bold). They are unique for the number of IP1 equal to 2, 4, 5, 7, and 8, respectively. However, there are two applicable choices for IP2 when IP1 = 3 and IP1 = 6, respectively. We may apply the 2nd criterion in determining the number of IPs, which is the length of the 95% confidence interval. Taking IP1 = 3 as an example, the corresponding values of RMSE and the length of the 95% confidence interval in the case of IP2 equal to 10 and 16 are listed in Table 5. Taking IP2 = 2 as a reference, it shows that the RMSE reduces by 4.8%, and the length increases by 84.6% when IP2 raises from 2 to 10. The RMSE reduces by 22.6%, and the length increases by 506.7% when IP2 raises from 2 to 16. Thus, 10 is the best choice for IP2 when IP1 = 3. According to the above method, 9 is the best choice for IP2 when IP1 = 6. The IP2 could be determined when IP1 ranges from 2 to 8 by the same method. Then, all the RMSE values listed in Table 4 are compared to each other, and 4 and 16 are the best-fit numbers for IP1 and IP2, respectively.
All the most suitable numbers for IPs in analyzing 8.85-year cycles of N 2 tide and 4.42-year cycles of 2N 2 tide at 4 stations are listed in Table 6. Meanwhile, the comparisons between the observational and hindcast amplitudes by the EHA methods of the N 2 and 2N 2 tide at four stations are exhibited in Figures 2 and 3, respectively. These hindcasts are obtained with the best-fit numbers of IPs listed in Table 6. Their amplitudes agree well with the observations except for extreme values, which may be caused by weather factors and tidal resonance. The errors between hindcasts and observations are small and reasonable. However, hindcasts at Halifax (Figure 3d) show a large difference with observations. This phenomenon implies that the number of IP3 = 1 in analyzing the 2N 2 tide is inapplicable at Halifax. In other words, the amplitude of the 8.85-year cycle is not a constant and its variation affects the water level remarkably at Halifax. Taking the numbers of IP3 for the 8.85-year cycle and IP4 for the 4.42-year cycle to be equal to 8 and 35, respectively, the results would improve remarkably (Figure 4 vs. Figure 3d). It is noted that the long-term trend shown in Figure 5 is not affected by the number of IPs. The numbers of IP3 and IP4 listed in Table 6 at Halifax are applied in Figures 6b and 7b to keep the consistency of the method.       Table 6).  Table 6).

Figure 2.
The primary amplitudes (red dots) and the hindcasts amplitudes by the EHA method (black lines) of N2 tide at 4 stations with the best-fit numbers of IPs for the 8.85-year cycle (listed in Table 6).

Long-Term Trend of N2 and 2N2 Tidal Amplitudes
As Ray and Foster [22] pointed out, the nuisance flooding, which is only triggered by high spring tides, has become more frequent since 2011 at Boston. Consequently, much attention should be paid to secular variation in water levels, as well as the long-term trend of tidal components. Figure 5 displays the long-term nonlinear trends of the N2 (blue dash line) and 2N2 (red solid line) tidal amplitudes at Boston, Eastport, Portland, and Halifax, respectively. As indicated in Section 3.1, the number of IPs is set equal to 0 and the parameter is equal to 3 to obtain the long-term nonlinear trend of N2 and 2N2 tidal amplitudes by the S_TIDE.
The N2 and 2N2 tidal amplitudes kept growing for almost one century at Boston; however, the trends of these two tides have reversed since the 1980s at Halifax. The trends are more complex at Eastport and Portland. The N2 tidal amplitude increases in the exponential profile at Eastport, while the amplitude of 2N2 tide started to decrease from the

Long-Term Trend of N 2 and 2N 2 Tidal Amplitudes
As Ray and Foster [22] pointed out, the nuisance flooding, which is only triggered by high spring tides, has become more frequent since 2011 at Boston. Consequently, much attention should be paid to secular variation in water levels, as well as the long-term trend of tidal components. Figure 5 displays the long-term nonlinear trends of the N 2 (blue dash line) and 2N 2 (red solid line) tidal amplitudes at Boston, Eastport, Portland, and Halifax, respectively. As indicated in Section 3.1, the number of IPs is set equal to 0 and the parameter is equal to 3 to obtain the long-term nonlinear trend of N 2 and 2N 2 tidal amplitudes by the S_TIDE.
The N 2 and 2N 2 tidal amplitudes kept growing for almost one century at Boston; however, the trends of these two tides have reversed since the 1980s at Halifax. The trends are more complex at Eastport and Portland. The N 2 tidal amplitude increases in the exponential profile at Eastport, while the amplitude of 2N 2 tide started to decrease from the 1980s. At Portland, the trend of the N 2 tidal amplitude keeps growing all the time, and the 2N 2 tidal amplitude increased until it reached a constant value in the 1980s.
In summary, for the four stations, the N 2 tidal amplitudes have increased in almost one century except for Halifax and the 2N 2 tidal amplitudes have decreased in the recent 30 or 40 years except for Boston. These two tides' amplitudes kept growing for one century only at Boston. This may account partially for the event that the nuisance flooding has become frequent at Boston by Ray and Foster [22].

The 8.85-and 4.42-Year Cycles
The changing amplitudes of the 8.85-year cycle of N2 tide and the 4.42-year cycl 2N2 tide are displayed in Figure 6 for the four stations of Boston, Eastport, Portland, a Halifax. The variations in the amplitudes of the 8.85-and 4.42-year cycles are the larg at Eastport among the four stations, and the least at Halifax. As mentioned in Section the water level record at Halifax is selected for contrasting the influence of topography the 8.85-and 4.42-year cycles. Halifax and Eastport are located at an approximate latitu and the most striking difference is the topography, and these two cycles at these two tions are compared in Figure 7. An obvious same phase can be seen in the 8.85-year cy at Eastport and Halifax from the middle 1930s to nowadays, and this phenomenon pears in the 4.42-year cycle from the middle 1930s to the middle 1980s. It is demonstra that the topography does not affect the phase of these cycles. However, the amplitude these cycles are greatly influenced by the topography and resonance.
The amplitude of the 8.85-year cycle has an approximately 20-year periodic osci tion at Eastport and Halifax with a consistent phase (Figure 6a). The 20-year oscillat appears at Boston too, but it lags that at Eastport and Halifax. At Portland, the 20-y oscillation has emerged since the 1950s, and the phase was in step with those at Eastp and Halifax between the 1990s and the 2010s. The variation in the 4.42-year cycle am tude of 2N2 tide is more complex than that in the 8.85-year cycle of N2 tide. The pheno enon of 20-year oscillation also appeared from the 1940s to the 2010s at Eastport (Fig  6b). From the early 1950s to the 1970s, there was a 10-year oscillation at Boston. Then, phase adjusted to 20 years. At Portland, there was no obvious oscillation until the 197 Afterwards, the phases kept in step with that at Eastport and Boston. At Halifax, the qu 20-year oscillation also appeared from the 1940s to the 1980s, which kept in step with t at Eastport.
The 95% confidence intervals for the 8.85-and 4.42-year cycles are presented in Fig   Figure 5. The long-term nonlinear trend of the N 2 (blue dash line) and 2N 2 (red solid line) tidal amplitudes at four stations.

The 8.85-and 4.42-Year Cycles
The changing amplitudes of the 8.85-year cycle of N 2 tide and the 4.42-year cycle of 2N 2 tide are displayed in Figure 6 for the four stations of Boston, Eastport, Portland, and Halifax. The variations in the amplitudes of the 8.85-and 4.42-year cycles are the largest at Eastport among the four stations, and the least at Halifax. As mentioned in Section 2.1, the water level record at Halifax is selected for contrasting the influence of topography on the 8.85-and 4.42-year cycles. Halifax and Eastport are located at an approximate latitude and the most striking difference is the topography, and these two cycles at these two stations are compared in Figure 7. An obvious same phase can be seen in the 8.85-year cycle at Eastport and Halifax from the middle 1930s to nowadays, and this phenomenon appears in the 4.42-year cycle from the middle 1930s to the middle 1980s. It is demonstrated that the topography does not affect the phase of these cycles. However, the amplitudes of these cycles are greatly influenced by the topography and resonance.
The amplitude of the 8.85-year cycle has an approximately 20-year periodic oscillation at Eastport and Halifax with a consistent phase (Figure 6a). The 20-year oscillation appears at Boston too, but it lags that at Eastport and Halifax. At Portland, the 20-year oscillation has emerged since the 1950s, and the phase was in step with those at Eastport and Halifax between the 1990s and the 2010s. The variation in the 4.42-year cycle amplitude of 2N 2 tide is more complex than that in the 8.85-year cycle of N 2 tide. The phenomenon of 20-year oscillation also appeared from the 1940s to the 2010s at Eastport (Figure 6b). From the early 1950s to the 1970s, there was a 10-year oscillation at Boston. Then, the phase adjusted to 20 years. At Portland, there was no obvious oscillation until the 1970s. Afterwards, the phases kept in step with that at Eastport and Boston. At Halifax, the quasi 20-year oscillation also appeared from the 1940s to the 1980s, which kept in step with that at Eastport.
The 95% confidence intervals for the 8.85-and 4.42-year cycles are presented in Figure 8. The confidence intervals are relatively large at the beginning and end of the time series for both cycles, of which the 8.85-year cycle is larger than those of the 4.42-year cycle. Nevertheless, they are all reasonable and demonstrate the validness of variations in the amplitude by the EHA method.

Discussion
Spatial and temporal variations are obvious among the four stations for the amplitudes and phases of 8.85-and 4.42-year cycles. Eastport kept in step with Halifax in the phase, but the modulation range is significantly different. A similar situation emerged along the Western Australian coast [12]. The lunar nodal cycle and the lunar perigean subharmonic occupy different areas. Haigh et al. [17] indicated that the main tidal constituents and tidal characteristics influence the spatial variations in the range and phase of the tidal modulations. The phase at Eastport agrees well with that at Halifax, and it may be proposed that the high water levels are in phase at these two tide gauges [17].
In addition, two key points are easily found in Figure 5. Firstly, different tidal constituents have different long-term nonlinear trends. Although the variations in the trends of the N2 and 2N2 tidal amplitudes are the same at Boston and Halifax, their change rates are different. Moreover, the trends are obviously different for the N2 and 2N2 tidal amplitudes at Eastport and Portland. Secondly, the nonlinear trend is localized and spatially varied. Taking the trend of N2 as an example, the monotonously increasing trend can only be found at Boston. The trends increased before the 1980s and decreased thereafter at Eastport, Portland, and Halifax, but the changing rates were different. A similar situation occurred for 2N2 tide.
Why did this phenomenon appear? Was there any change in 100 years at these four stations? The most obvious changes are the sea-level rise (SLR) and coastal defense. Figure  9 shows the long-term linear trend (the IPs is equal to 0 and the parameter is equal to 2) of the total sea levels at the four stations, and they increased evidently. Some papers have noticed the effect of SLR on the tide, and the localized responses have been explored on the European shelf [34][35][36]. Besides, the coastal defense could exacerbate the increases in the tidal range induced by the SLR, and flooding over land induces damping of the tidal amplitude [37][38][39]. In addition to local change induced by the SLR and coastal defense, Jay [40] proposed two climate-related mechanisms: variations in stratification and mean vorticity of the upper ocean induced by large-scale changes in wind-driven circulation. The reanalysis model data are useful to verify the proposed mechanisms, such as temperature, salinity, sea surface height, and horizontal flow velocity.

Discussion
Spatial and temporal variations are obvious among the four stations for the amplitudes and phases of 8.85-and 4.42-year cycles. Eastport kept in step with Halifax in the phase, but the modulation range is significantly different. A similar situation emerged along the Western Australian coast [12]. The lunar nodal cycle and the lunar perigean subharmonic occupy different areas. Haigh et al. [17] indicated that the main tidal constituents and tidal characteristics influence the spatial variations in the range and phase of the tidal modulations. The phase at Eastport agrees well with that at Halifax, and it may be proposed that the high water levels are in phase at these two tide gauges [17].
In addition, two key points are easily found in Figure 5. Firstly, different tidal constituents have different long-term nonlinear trends. Although the variations in the trends of the N 2 and 2N 2 tidal amplitudes are the same at Boston and Halifax, their change rates are different. Moreover, the trends are obviously different for the N 2 and 2N 2 tidal amplitudes at Eastport and Portland. Secondly, the nonlinear trend is localized and spatially varied. Taking the trend of N 2 as an example, the monotonously increasing trend can only be found at Boston. The trends increased before the 1980s and decreased thereafter at Eastport, Portland, and Halifax, but the changing rates were different. A similar situation occurred for 2N 2 tide.
Why did this phenomenon appear? Was there any change in 100 years at these four stations? The most obvious changes are the sea-level rise (SLR) and coastal defense. Figure 9 shows the long-term linear trend (the IPs is equal to 0 and the parameter is equal to 2) of the total sea levels at the four stations, and they increased evidently. Some papers have noticed the effect of SLR on the tide, and the localized responses have been explored on the European shelf [34][35][36]. Besides, the coastal defense could exacerbate the increases in the tidal range induced by the SLR, and flooding over land induces damping of the tidal amplitude [37][38][39]. In addition to local change induced by the SLR and coastal defense, Jay [40] proposed two climate-related mechanisms: variations in stratification and mean vorticity of the upper ocean induced by large-scale changes in wind-driven circulation. The reanalysis model data are useful to verify the proposed mechanisms, such as temperature, salinity, sea surface height, and horizontal flow velocity.  The reasons for this oscillation are proposed as follows. A 20-to 30-year timescale oscillation was reported in the North Atlantic by analyzing observational and model datasets, and it is triggered by the westward propagation of subsurface temperature anomalies [41]. Afterwards, a 20-year mode coupled ocean-sea ice-atmosphere variability in the North Atlantic was proposed by Escudier et al. [42]. This mode is driven by the westpropagation of the near-surface temperature and salinity anomalies, which leads to anomalous sea ice melting. Then, the anomalous surface atmospheric temperature forces sealevel pressure anomalies. It is hard to verify the relationship between the amplitude oscillation of the 8.85-year cycle of N2 tide and the 4.42-year cycle of 2N2 tide and the 20-year mode [42]. However, the western North Atlantic Ocean is a complex land-ocean-atmosphere system [43], and any variation that appeared at one place may interact with and propagate to other places in the coupled system.

Conclusions
This paper focused on the 8.85-year cycle of N2 tide and the 4.42-year cycle of 2N2 tide surrounding the Bay of Fundy and the Gulf of Maine at four stations: Boston, Portland, Eastport, and Halifax. To obtain the temporal variation in the amplitudes of the cycles, the enhanced harmonic analysis method and the independent point scheme were implemented. The method and scheme were effective as the hindcast results compared well with primary data. Enhanced harmonic analysis results showed that the long-term trends of N2 and 2N2 tidal amplitudes both increased at Boston, while they increased before the 1980s and decreased thereafter at Halifax. At Eastport and Portland, they increased at different rates before the 1980s but divided thereafter. The long-term trends of  The reasons for this oscillation are proposed as follows. A 20-to 30-year timescale oscillation was reported in the North Atlantic by analyzing observational and model datasets, and it is triggered by the westward propagation of subsurface temperature anomalies [41]. Afterwards, a 20-year mode coupled ocean-sea ice-atmosphere variability in the North Atlantic was proposed by Escudier et al. [42]. This mode is driven by the west-propagation of the near-surface temperature and salinity anomalies, which leads to anomalous sea ice melting. Then, the anomalous surface atmospheric temperature forces sea-level pressure anomalies. It is hard to verify the relationship between the amplitude oscillation of the 8.85-year cycle of N 2 tide and the 4.42-year cycle of 2N 2 tide and the 20-year mode [42]. However, the western North Atlantic Ocean is a complex land-oceanatmosphere system [43], and any variation that appeared at one place may interact with and propagate to other places in the coupled system.

Conclusions
This paper focused on the 8.85-year cycle of N 2 tide and the 4.42-year cycle of 2N 2 tide surrounding the Bay of Fundy and the Gulf of Maine at four stations: Boston, Portland, Eastport, and Halifax. To obtain the temporal variation in the amplitudes of the cycles, the enhanced harmonic analysis method and the independent point scheme were implemented. The method and scheme were effective as the hindcast results compared well with primary data. Enhanced harmonic analysis results showed that the long-term trends of N 2 and 2N 2 tidal amplitudes both increased at Boston, while they increased before the 1980s and decreased thereafter at Halifax. At Eastport and Portland, they increased at different rates before the 1980s but divided thereafter. The long-term trends of the tidal amplitudes at the four stations were probably triggered by the SLR, coastal defense, and other possible climate-related mechanisms. Moreover, the amplitude variations in the 8.85-year cycle of N 2 tide and the 4.42-year cycle of 2N 2 were extracted by the EHA method, and they varied spatially. The comparison at Eastport and Halifax revealed that the topography influences the amplitudes rather than the phase of these cycles. A quasi 20-year oscillation existed in the amplitude variations of the two cycles. For the 8.85-year cycle, the oscillation was remarkable at the four stations but with different amplitudes. For the 4.42-year cycle, the oscillation appeared throughout the century at Eastport and from the 1970s to the 2010s at Boston and Portland. It was also observed from the 1940s to the 1980s at Halifax. The quasi 20-year oscillation may be associated with the 20-year mode in the North Atlantic Ocean [43]. More measurements are needed to verify this hypothesis.