# Large-scale changes of the semidiurnal tide along North Atlantic coasts from 1846 to 2018

Lucia Pineau-Guillou, Pascal Lazure, Guy Woppelmann

► To cite this version:

Lucia Pineau-Guillou, Pascal Lazure, Guy Woppelmann. Large-scale changes of the semidiurnal tide along North Atlantic coasts from 1846 to 2018. Ocean Science, 2021, 17 (1), pp.17-34. 10.5194/os-17-17-2021. hal-03126135

HAL Id: hal-03126135

https://hal.science/hal-03126135

Submitted on 12 Aug 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

# Large-scale changes of the semidiurnal tide along North Atlantic coasts from 1846 to 2018

Lucia Pineau-Guillou¹, Pascal Lazure¹, and Guy Wöppelmann²

¹IFREMER, CNRS, IRD, UBO, Laboratoire d'Océanographie Physique et Spatiale, UMR 6523, IUEM, Brest, France

²LIENSS, Université de la Rochelle-CNRS, La Rochelle, France

**Correspondence:** Lucia Pineau-Guillou (lucia.pineau.guillou@ifremer.fr)

Received: 29 May 2020 – Discussion started: 15 June 2020

Revised: 5 November 2020 – Accepted: 10 November 2020 – Published: 4 January 2021

**Abstract.** We investigated the long-term changes of the principal tidal component $M_2$ along North Atlantic coasts, from 1846 to 2018. We analysed 18 tide gauges with time series starting no later than 1940. The longest is Brest with 165 years of observations. We carefully processed the data, particularly to remove the 18.6-year nodal modulation. We found that $M_2$ variations are consistent at all the stations in the North-East Atlantic (Cuxhaven, Delfzijl, Hoek van Holland, Newlyn, Brest), whereas some discrepancies appear in the North-West Atlantic. The changes started long before the 20th century and are not linear. The secular trends in $M_2$ amplitude vary from one station to another; most of them are positive, up to 2.5 mm/yr at Wilmington since 1910. Since 1990, the trends switch from positive to negative values in the North-East Atlantic. Concerning the possible causes of the observed changes, the similarity between the North Atlantic Oscillation and $M_2$ variations in the North-East Atlantic suggests a possible influence of the large-scale atmospheric circulation on the tide. Our statistical analysis confirms large correlations at all the stations in the North-East Atlantic. We discuss a possible underlying mechanism. A different spatial distribution of mean sea level (corresponding to water depth) from one year to another, depending on the low-frequency sea-level pressure patterns, could impact the propagation of the tide in the North Atlantic basin. However, the hypothesis is at present unproven.

## 1 Introduction

Tides have been changing due to non-astronomical factors since the 19th century (Haigh et al., 2019; Talke and Jay, 2020). In the North Atlantic, secular variations have been observed at individual tide gauge stations, e.g. Brest (Cartwright, 1972; Wöppelmann et al., 2006; Pouvreau et al., 2006; Pouvreau, 2008), Newlyn (Araújo and Pugh, 2008; Bradshaw et al., 2016), New York (Talke et al., 2014), and Boston (Talke et al., 2018), but also at regional scale, e.g. Gulf of Maine (Doodson, 1924; Godin, 1995; Ray, 2006; Ray and Talke, 2019), at the North Atlantic basin scale (Müller, 2011), and at a quasi-global scale (Woodworth, 2010; Müller et al., 2011; Mawdsley et al., 2015). Long-term changes in tidal constituents are rather small at some coastal stations but tend to be statistically significant. The order of magnitude of these changes varies spatially and may reach a few centimetres per century for $M_2$ amplitude. For example, Ray and Talke (2019) found trends varying from –1 to 8 cm per century in the Gulf of Maine over the last century. Woodworth et al. (2010) and Müller et al. (2011) found trends of a few percent per century in the Atlantic. The changes can be larger in many estuaries and rivers (Talke and Jay, 2020).

The physical causes of these changes can be multiple and difficult to disentangle. In particular, the complexity comes from the possible interaction between local and large-scale causes. Changes may have a local-scale origin, such as changes in the nearby environment (e.g. harbour development, deepening of channels, dredging, siltation) or changes in the instrumentation (e.g. tide gauge technology, observational location, instrumental errors). For example, Familkhalili and Talke (2016) show that mean tidal range at Wilmington has doubled since the 1880s, due to channel deepening in
the Cape Fear River estuary. Changes may also have a large-scale origin, i.e. regional or global. Haigh et al. (2019) reported several possible large-scale mechanisms: (1) tectonics and continental drift, (2) water depth changes due to mean sea level rise or geological processes such as the Earth's surface glacial isostatic adjustment (Müller et al., 2011; Pickering et al., 2017; Schindelegger et al., 2018), (3) shoreline position, (4) extent of sea-ice cover (Müller et al., 2014), (5) sea-bed roughness, (6) ocean stratification which may modify the internal tides and bottom friction over continental shelves (Müller, 2012), (7) non-linear interactions, and (8) radiational forcing (Ray, 2009).

Several authors have explored mean sea level (MSL) rise as a potential mechanism to explain $M_2$ changes. For example, simulations by Pickering et al. (2012) show that a 2 m sea level rise could modify $M_2$ from -20 to 20 cm around the whole ocean. Idier et al. (2017) show that depending on the location, the changes can account for ±15% of the regional sea level rise. Schindelegger et al. (2018) find changes of about 1%–5% of the sea level rise. Beyond MSL rise, other mechanisms have been explored to explain $M_2$ changes. For example, Colosi and Munk (2006) attribute the changes of $M_2$ amplitude at Honolulu, Hawaii, to a 28° rotation of the internal tide vector in response to ocean warming. Ray and Talke (2019) suggest that long-term changes in stratification could play a role in the Gulf of Maine. Müller (2011) suggests a possible link between $M_2$ changes and atmospheric dynamics in the North Atlantic; he reported that the time series of the North Atlantic Oscillation (NAO) show similar characteristics to those of the tidal amplitudes and phases. In the Gulf of Maine, Pan et al. (2019) suggest that changes in the response of the nodal modulation of the $M_2$ tide from 1970s to 2013 may be linked with the NAO. In Southeast Asian waters, Devlin et al. (2018) show that the impact of atmospheric circulation (via the wind stress, through Ekman current) on the $M_2$ seasonal cycle may be significant and comparable to the effect of permanent (geostrophic) currents. In the North Sea, Huess and Andersen (2001) explain a large part of $M_2$ seasonal cycle by the role of atmospheric dynamics, whereas Müller et al. (2014) and Gräwe et al. (2014) suggest a major role of the thermal stratification. These examples show the diversity of mechanisms that play a role in tide changes. In the present paper, we focus on the role of MSL and atmospheric dynamics.

This paper has two main objectives. The first is to characterize the secular changes of the $M_2$ tide over the North Atlantic. We focus on the longest time series, i.e. starting no later than 1940. This approach is complementary to previous studies investigating $M_2$ changes focusing on smaller spatial scales, e.g. Brest (Pouvreau et al., 2006; Pouvreau, 2008), Gulf of Maine (Ray, 2006; Ray and Talke, 2019), or focusing on shorter temporal scales, i.e. recent decades (Woodworth, 2010; Müller, 2011). The second objective is to detect if there is any large-scale coherence in the observed changes in the North Atlantic and investigate the possible link with the atmospheric circulation, already mentioned by Müller et al. (2011), on the basis of qualitative criteria. Here, we further provide quantitative insights into the possible influence of the NAO and discuss a possible NAO-related climate mechanism that can partly explain the observed changes.

The paper is organized as follows. The first section below describes the data: the sea level data (i.e. tide gauges and their processing) and the atmospheric data (i.e. climate indices and sea level pressure data). The following section presents the results (i.e. $M_2$ variations and trends). We then discuss a possible link between the observed tidal changes and MSL, as well as climate indices.

## 2 Data

### 2.1 Sea level data

#### 2.1.1 Tide gauge selection

The tide gauge data were retrieved from the University of Hawaii Sea Level Center (UHSLC, website accessed April 2020). The dataset consists of 249 stations in the Atlantic Ocean, with hourly sea level observations. Two additional long-term stations – Delfzijl and Hoek van Holland – were provided by Rijkswaterstaat (RWS) in the Netherlands.

We selected the stations following three criteria: time series (1) starting before 1940, (2) with at least 80 years of data, and (3) with tidal amplitude significant enough to detect trends, i.e. $M_2$ amplitude larger than 10 cm. Note that we selected only years with at least 75% of data (see Sect. 2.1.2). Only 24 stations among the 249 from UHSLC fulfilled the two first criteria (Fig. 1). They are all located in the Northern Hemisphere. On the east side of the North Atlantic, Stockholm, Gedser, Hornbaek, Tregde and Marseille were discarded due to too small of an $M_2$ amplitude (i.e. lower than 10 cm). These stations are located in the Baltic Sea (Stockholm, Gedser), in the strait separating the Baltic and the North Sea (Hornbaek), in the North Sea (Tregde), and in the Mediterranean Sea (Marseille). On the west side of the North Atlantic, Galveston, Pensacola and Cristobal were also discarded due to too small of a tidal amplitude (i.e. lower than 10 cm). These stations are located in the Gulf of Mexico (Galveston, Pensacola) and the Caribbean Sea (Cristobal).

Finally, 18 stations followed the three criteria detailed above and were selected for this study (see stations in bold in Fig. 1, 16 stations are from UHSLC, and 2 from RWS). Among them, 5 are located on the North-East Atlantic coasts (Newlyn, Brest, Hoek van Holland, Delfzijl and Cuxhaven – note that Hoek van Holland, Delfzijl and Cuxhaven are located in the North Sea) and 13 are located on the North-West Atlantic coasts (Halifax, Eastport, Portland, Boston, Newport, New London, New York, Atlantic City, Lewes, Wilmington, Charleston, Fort Pulaski and Key West).
<figure><img src="image_4.png" /><figcaption>Figure 1. Tide gauges in the North Atlantic. Stations with time series starting before 1940 and longer than 80 years are labelled. Stations selected for this study are in bold.</figcaption></figure>

The main characteristics of the 18 selected stations are summarized in Table 1. Among them, only Brest, Hoek van Holland and Halifax started in the 19th century, in 1846, 1879 and 1896 respectively (Table 1, column 2). The number of years with data for each station varies between 81 and 165 years, Brest being the longest time series (Table 1, column 3).

## 2.1.2 Data processing

Harmonic analysis was performed in order to compute the $M_2$ amplitude. We used the MAS program (Simon, 2007, 2013), developed by the French Hydrographic Office (SHOM). This program gives results similar to the T_Tide harmonic analysis toolbox (Pawlowicz et al., 2002). For instance, Pouvreau et al. (2006) found no differences in the yearly amplitudes of $M_2$ at Brest over the period 1846 to 2005 using either T_Tide or MAS. Hourly time series were analysed yearly. Note that at Delfzijl and Hoek van Holland, data had to be interpolated every hour before 1970, as the temporal sampling was 3 h. (We checked with hourly time series from recent years (1971–2018) that 3-hourly sampling did not result in a significant reduction of $M_2$ amplitude in a tidal analysis compared to hourly sampling.) We processed only years with at least 75 % of data, to avoid seasonal modulation affecting the computed amplitudes. In the North Atlantic, $M_2$ is affected by a seasonal variation of a few percent (Pugh and Vassie, 1976; Huess and Andersen, 2001; Müller et al., 2014; Gräwe et al., 2014). Considering only years with at least 75 % of data resulted in excluding up to 15 years for a given station (Table 1, columns 3 and 4). We carefully

removed the nodal modulation of $M_2$ amplitude (Simon, 2007, 2013), as described briefly in Appendix A. Finally, 3 station-years were discarded due to problems in the record (1953 and 1962 at Delfzijl, 1953 at Hoek van Holland), and 2 more station-years due to doubtful $M_2$ values (1972 at Eastport, 1978 at Newport).

At all the stations, we computed the normalized $M_2$ amplitude, removing the average and dividing by the standard deviation over the period 1910–2010:

$$ \text{normalized } M_2(t) = \frac{M_2(t) - \overline{M_2[1910,2010]}}{\sigma_{M_2[1910,2010]}} \quad (1) $$

The average $\overline{M_2}$ and standard deviation $\sigma_{M_2}$ over the 1910–2010 period are given in Table 1 (column 5). The idea is to scale the data in order to compare all the stations together.

## 2.2 Atmospheric data

### 2.2.1 Climate indices

We investigated the correlation between secular changes in the tide and climate indices, such as the North Atlantic Oscillation (NAO) or the Arctic Oscillation (AO) – also called Northern Annular Mode (NAM) (Hurrell, 1995; Hurrell and Deser, 2009; Thompson and Wallace, 2000; Thompson et al., 2000). These climate indices are related to the distribution of atmospheric masses. They are based on the difference of average sea-level pressure between two centres of actions (i.e. stations) over long periods (e.g. monthly, seasonal, annual).

The NAO is the major pattern of weather and climate variability over the Northern Hemisphere (Hurrell,
<table><caption>**Table 1.** Main characteristics of tide gauge records selected for this study. Name of the station, time span, number of years with data, number of years analysed (i.e. with at least 75% of data), *M*₂ average amplitude and standard deviation over the period 1910–2010, *M*₂ nodal modulation, estimated trends in *M*₂ amplitude since 1910 and since 1990 up to 2018 in each case (standard errors are 1σ, considering the noise content in the time series; see text).</caption><thead><tr><th>Name</th><th>Time span</th><th>Number of years with data</th><th>Number of years analysed</th><th>M&#773;<sub>2</sub> (cm) [1910–2010]</th><th>M<sub>2</sub> nod. mod. f<sub>nod</sub></th><th>M<sub>2</sub> trends since 1910 (mm/yr)</th><th>M<sub>2</sub> trends since 1990 (mm/yr)</th></tr></thead><tbody><tr><td>Cuxhaven</td><td>1918–2018</td><td>102</td><td>101</td><td>135.05 ± 3.68</td><td>1.8 %</td><td>0.68 ± 0.56</td><td>-0.47 ± 0.78</td></tr><tr><td>Delfzijl</td><td>1879–2018</td><td>138</td><td>138</td><td>125.58 ± 6.96</td><td>1.7 %</td><td>2.02 ± 0.59</td><td>-0.09 ± 0.28</td></tr><tr><td>Hoek van Holland</td><td>1900–2018</td><td>88</td><td>82</td><td>76.95 ± 2.63</td><td>0.8 %</td><td>0.85 ± 0.32</td><td>-0.45 ± 0.17</td></tr><tr><td>Newlyn</td><td>1916–2016</td><td>102</td><td>98</td><td>170.66 ± 0.75</td><td>3.3 %</td><td>0.14 ± 0.09</td><td>-0.28 ± 0.49</td></tr><tr><td>Brest</td><td>1846–2018</td><td>165</td><td>158</td><td>204.54 ± 0.91</td><td>3.8 %</td><td>0.13 ± 0.11</td><td>-0.36 ± 0.18</td></tr><tr><td>Halifax</td><td>1896–2012</td><td>99</td><td>95</td><td>62.83 ± 0.64</td><td>3.7 %</td><td>-0.15 ± 0.05</td><td>0.32 ± 0.35</td></tr><tr><td>Eastport</td><td>1930–2018</td><td>90</td><td>82</td><td>263.51 ± 2.50</td><td>2.5 %</td><td>0.80 ± 0.21</td><td>1.01 ± 0.65</td></tr><tr><td>Portland</td><td>1910–2018</td><td>109</td><td>104</td><td>135.07 ± 1.84</td><td>2.8 %</td><td>0.56 ± 0.06</td><td>0.72 ± 0.23</td></tr><tr><td>Boston</td><td>1922–2018</td><td>98</td><td>96</td><td>136.57 ± 1.03</td><td>2.9 %</td><td>0.27 ± 0.06</td><td>0.42 ± 0.24</td></tr><tr><td>Newport</td><td>1931–2018</td><td>89</td><td>84</td><td>50.86 ± 0.41</td><td>4.1 %</td><td>-0.09 ± 0.03</td><td>-0.03 ± 0.12</td></tr><tr><td>New London</td><td>1939–2018</td><td>81</td><td>76</td><td>35.93 ± 0.25</td><td>3.5 %</td><td>0.06 ± 0.02</td><td>0.03 ± 0.09</td></tr><tr><td>New York</td><td>1921–2018</td><td>95</td><td>80</td><td>65.13 ± 0.83</td><td>3.7 %</td><td>0.33 ± 0.07</td><td>0.93 ± 0.15</td></tr><tr><td>Atlantic City</td><td>1912–2018</td><td>107</td><td>101</td><td>58.48 ± 0.31</td><td>3.8 %</td><td>0.00 ± 0.03</td><td>-0.18 ± 0.09</td></tr><tr><td>Lewes</td><td>1919–2018</td><td>85</td><td>72</td><td>59.91 ± 0.43</td><td>3.1 %</td><td>-0.06 ± 0.09</td><td>-0.33 ± 0.07</td></tr><tr><td>Wilmington</td><td>1936–2018</td><td>84</td><td>82</td><td>56.84 ± 6.16</td><td>1.7 %</td><td>2.51 ± 0.46</td><td>1.80 ± 0.44</td></tr><tr><td>Charleston</td><td>1901–2018</td><td>101</td><td>100</td><td>76.40 ± 1.33</td><td>3.0 %</td><td>0.32 ± 0.18</td><td>-0.02 ± 0.10</td></tr><tr><td>Fort Pulaski</td><td>1936–2018</td><td>84</td><td>78</td><td>100.60 ± 1.01</td><td>3.1 %</td><td>0.18 ± 0.14</td><td>-0.01 ± 0.21</td></tr><tr><td>Key West</td><td>1913–2018</td><td>106</td><td>104</td><td>17.50 ± 0.36</td><td>2.9 %</td><td>0.08 ± 0.04</td><td>0.13 ± 0.02</td></tr></tbody></table>

1995; Hurrell and Deser, 2009). Variations of NAO drive the climate variability over Europe and North America (Hurrell et al., 2003). We used the wintertime (December to March) Hurrell station-based NAO Index (retrieved from https://climatedataguide.ucar.edu/climate-data/hurrell-north-atlantic-oscillation-nao-index-station-based, last access: April 2020). It is based on the difference of normalized average winter sea-level pressure between Lisbon (Portugal) and Stykkishólmur/Reykjavik (Iceland). The normalization involves removing the mean (1864–1983) and dividing by the long-term standard deviation. The NAO index covers the period 1864–2019.

The Arctic Oscillation (AO) is another index which resembles the NAO index. It is defined as the first EOF of Northern Hemisphere winter sea-level pressure data (Thompson and Wallace, 1998, 2000; Thompson et al., 2000). The AO index is highly correlated with the NAO. We used the wintertime Hurrell AO index (retrieved from https://climatedataguide.ucar.edu/climate-data/hurrell-wintertime-slp-based-northern-annular-mode-nam-index, last access: April 2020). The AO index covers the period 1899–2019.

To remove the interannual variability and estimate low-frequency variations, climate indices were low-pass filtered with a 9-year mean filter.

### 2.2.2 Sea level pressure

We employed the Twentieth Century Reanalysis (20CR version 3 dataset) (Compo et al., 2011; Slivinski et al., 2019), a historic weather reconstruction from 1836 to 2015, with a 1° gridded global coverage. However, we made use only of data from 1850 to be more consistent with the temporal coverage of the tide gauge measurements. This will be discussed in Sect. 4.

## 3 Results

### 3.1 $M_2$ variations

For the North-East Atlantic, the variations of normalized $M_2$ amplitude are presented in Fig. 2a.

The first result is that since 1910, the variations show similar patterns at all the stations; $M_2$ amplitude decreases up to the 1960s, then increases, and decreases again since the 1990s. This suggests that these changes are probably due to large-scale processes, rather than local effects due to changes in the environment (e.g. harbour development, dredging, silting) or instrumentation errors. The similar patterns between Brest and Cuxhaven may be surprising, as Cuxhaven is located in the North Sea, and not in the open Atlantic Ocean, and far away from Brest, around 1300 km. This indicates that the spatial scale of the processes responsible for these changes must be at least as large as the North-East Atlantic.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Normalized annual $M_2$ amplitude (a) in the North-East Atlantic (b) in the North-West Atlantic, stations with positive trends (c) in the North-West Atlantic, stations with negative or no trend. The stars in (b) in the 1860s correspond to $M_2$ amplitude at Eastport and Portland from Ray and Talke (2019), and New York from Talke et al. (2014), after normalization (Eq. 1).</figcaption></figure>

Different authors have noticed the increase of tidal range from 1960 to 1990 in the southern North Sea. Hollebrandse (2005) found a gradual increase during the period 1955–1980 at all the stations of the Dutch coast (five stations including Hoek van Holland) and the German coast (seven stations). Mudersbach et al. (2013) found a significant increase in $M_2$ amplitude at Cuxhaven since around the mid-1950s. Note that Cuxhaven is located in the German Bight; shallow depths and the shape of the coastline may induce some amplification. Variations in $M_2$ at Cuxhaven could therefore

be sensitive to local effects, such as the migration of the underwater channels and the evolution of the tidal flats (Jacob et al., 2016). Moreover, Cuxhaven is located in the Elbe estuary, and some river engineering works, such as narrowing and deepening, may induce tidal amplification (Winterwerp and Wang, 2013; Winterwerp et al., 2013).
Before 1910, normalized $M_2$ values are higher at Brest than at Delfzijl. The construction of dykes that have gradually closed the harbour of Brest since the end of the 19th century may have altered the tide at Brest. The high values before 1910 may be due to local changes, in addition to large-scale changes. To go further, the potential role of these successive constructions needs to be investigated (Wikipedia contributors, 2020). Cartwright (1972) made a first attempt to evaluate the influence of reducing the width of access to the harbour but did not take into account a potential role of dredging, for which we have no information. This example underlines the complexity of interpretation of the variations when changes of local and large-scale origin occur at the same time. Note that in the following, we focus mainly on the 20th century, as most of the stations start after 1900 (15 out of 18 stations).

The second result is that there is no obvious linear trend in $M_2$ variations, but rather break or change points, $M_2$ increasing and then decreasing, depending on the periods considered. Overall, $M_2$ decreases from 1910 until 1960, increases again until 1980–1990, to finally decrease since 1990; note that the curve flattens between 1920 and 1940. Pouvreau et al. (2006) already noticed these variations at Brest and Newlyn and suggested a long-period oscillation of around 140 years, rather than a steady secular trend. A careful analysis of the harmonic development of the tidal potential showed that no tidal component could explain this oscillation. Similarly, no linear combination of tidal harmonic components could explain it (Pouvreau et al., 2006). This indicates that these variations are not due to an astronomical component. However, in contrast to Brest, $M_2$ at Delfzijl stays flat between 1880 and 1920. The decrease observed at Brest between 1880 and 1920 may be due to harbour development and/or dredging (see above). This underlines the importance of sea level data archaeology, for research studies related to long-term changes (Pouvreau, 2008; Woodworth et al., 2010; Marcos et al., 2011; Talke and Jay, 2013, 2017; Ray and Talke, 2019; Bradshaw et al., 2015, 2020).

The third result is that changes in $M_2$ have not the same order of magnitude at each station (see Fig. B1 in Appendix B for time series of $M_2$). Note that Fig. 2 represents normalized $M_2$, i.e. removing the average and dividing by the standard deviation. The order of magnitude of unnormalized $M_2$ changes is roughly the same at Brest and Newlyn (standard deviations of 0.9 and 0.8 cm respectively, Table 1, column 5), but more than three times larger at Cuxhaven (standard deviation of 3.7 cm), and even larger at Delfzijl (standard deviation of 7 cm). This suggests that the North Sea may be more sensitive to the processes responsible for these changes. Note also that the environmental setting of Cuxhaven and Delfzijl in the Elbe and Ems estuaries, respectively, could introduce some amplification (Winterwerp and Wang, 2013; Winterwerp et al., 2013).

For the North-West Atlantic, the variations of normalized $M_2$ amplitude are presented in Fig. 2b and c. The first feature

is that $M_2$ amplitude varies differently in the North-West and in the North-East Atlantic. The second is that there are discrepancies between stations, even when close to each other (e.g. Atlantic City and Lewes). We split the stations into two groups, in order to facilitate the detection of patterns, each being consistent in terms of trends: one with positive trend (group 1 in Fig. 2b), the other one with negative or no trend (group 2 in Fig. 2c).

The first group (with positive trends) consists of nine stations (Fig. 2b). Three outcomes can be highlighted. The first is that $M_2$ amplitude has increased overall since 1900. However, between 1980 and 1990, all the stations slightly decrease, and since 1990 they have increased again. The second outcome is that the rate of increase is very different from one station to another (keeping in mind that $M_2$ is normalized by standard deviation in Fig. 2). Portland is increasing 1.4 times faster than Charleston (standard deviations being respectively of 1.82 and 1.33 cm) and 28 times faster than Key West (standard deviation being only 0.36 cm at Key West). The large increase in Portland may be explained by some amplification in the Gulf of Maine. In many semienclosed basins, resonance leads to tidal amplification (Talke and Jay, 2020; Haigh et al., 2019). In the Gulf of Maine, Ray and Talke (2019) reported that the tides in the Gulf are in resonance, with a natural resonance frequency close to the $N_2$ tide (Garrett, 1972; Godin, 1993). Tides may be then very sensitive to any changes in the environment (e.g. basin configuration – shape, depth – but also external forcing). The third outcome, and probably the most interesting one, is related to the values of $M_2$ at Eastport, Portland and New York in the 1860s, estimated from Ray and Talke (2019) and Talke et al. (2014), and represented (after normalization) as stars in Fig. 2b. These values are not consistent with the positive linear trends observed since 1900, which provides some consistency with the hypothesis formulated from the analyses of the data prior to the 20th century in Fig. 2a: long-term variations introduce some breaks or change points, $M_2$ increasing and then decreasing, depending on the periods considered. The decrease observed between the 1870s and 1920s at the four stations (Brest, Eastport, Portland, New York) suggests a possible large-scale signal, in addition to local processes.

The second group (with negative or no trend) consists of four stations (Fig. 2c). Two points can be highlighted. The first is that $M_2$ decreases overall for Halifax, Newport and Lewes. This is less clear for Atlantic City, which is quite noisy and shows no significant trend. The second point is that at Halifax, $M_2$ values in 1896–1897 are higher than those after 1920. This suggests that the decrease may have started before the 20th century. However note that at Halifax, there is a long gap in the data recording (1898–1919), which raises the possibility of an instrumentation origin in the observed decrease of the $M_2$ amplitude.
## 3.2 Estimated trends

We estimated the trends for $M_2$ amplitude at each station, using linear regression. We computed the trends over two periods: 1910–2018, which corresponds roughly to the whole period of data (only five stations start before 1910), and 1990–2018, which corresponds to recent decades. Some tests showed that the later results were not very sensitive to the start date (moving 1990 to 1985 or 1995). The trend uncertainties were estimated considering the noise content in the time series using SARI software (Santamaría-Gómez, 2019). The noise was modelled as a white plus power-law noise, whose spectral index was found to be close to -1 (flicker noise). The results are summarized in Table 1 (columns 7 and 8) and Figs. 3 and 4.

The trends estimated since 1910 vary significantly from one station to another (Fig. 3). They are positive overall (up to 2.5 mm/yr at Wilmington), which is consistent with previous findings (Araújo and Pugh, 2008; Ray, 2009; Woodworth, 2010; Müller et al., 2011; Ray and Talke, 2019). They are slightly negative at three stations (Halifax, Newport, Lewes), and one station shows no trend (Atlantic City). The estimates are statistically consistent with those found previously by different authors (e.g. 0.14 ± 0.09 mm/yr at Newlyn compared to 0.19 ± 0.03 mm/yr in Araújo and Pugh (2008), 0.56 ± 0.06 mm/yr in Portland, compared to 0.59 ± 0.04 mm/yr in Ray and Talke, 2019). Note that our error bars are larger, because we considered the noise content in the time series as a white noise plus power law noise (we obtained the same error bars considering white noise only). In the North-East Atlantic, the trends are consistent with each other (in terms of sign), which is not surprising as the stations vary similarly (Fig. 2a).

The largest trends since 1910 are mainly observed in semi-closed basins: Wilmington in the Cape Fear River estuary, Delfzijl in Ems estuary, Cuxhaven in Elbe estuary, and Eastport and Portland in the Gulf of Maine. This suggests a possible amplification due to resonance effects (e.g. Gulf of Maine) and/or propagation in shallow waters (e.g. Cuxhaven), in addition to local effects. The stations located in estuaries or in a harbour with a channel may have been subject to dredging. Channel deepening increases the water depths, which reduces the effective drag and leads to tidal range amplification. This effect may be particularly large in estuaries (Ralston et al., 2019; Talke and Jay, 2020) and may explain the larger trends at Wilmington (Familkhalili and Talke, 2016) and Delfzijl. Finally, the shifting locations of amphidromic points could also play a role (Haigh et al., 2019). In the North Sea, different authors show a possible migration of the present-day amphidromes, under a 2 m sea-level rise scenario (Pickering et al., 2012; Idier et al., 2017).

The trends estimated since 1990 are quite different from those estimated since 1910 (Figs. 3 and 4), with more stations with negative trends: 9 stations out of 18 have post-1990 negative trends, whereas only 3 stations out of 18 have

post-1910 negative trends (Table 1, columns 7 and 8). In the North-East Atlantic, they all switch from positive to negative trends. This underlines (1) some spatially coherent changes in recent decades (Müller, 2011; Ray and Talke, 2019) and (2) the difficulty in estimating long-term trends from short records (i.e. less than 30 years), especially if the data are noisy (interannual variability) and the underlying processes non-linear (change points).

The trends have to be interpreted very carefully as the $M_2$ variations are not linear and may increase or decrease depending on the years; as a consequence, the estimated trends depend strongly on the period considered to estimate it. The interannual variability also plays an important role, and when substantial, trends can vary depending on the computational period. For example, at Cuxhaven, the large interannual variability leads to a large uncertainty on the trend computed since 1990 (−0.47 ± 0.78 mm/yr).

# 4 Discussion

## 4.1 Possible link with mean sea level rise

The MSL rise could partly explain $M_2$ changes. Simulations show that MSL rise can result in an change of $M_2$ up to ±10% of the rise (Pickering et al., 2017; Idier et al., 2017; Schindelegger et al., 2018). Schindelegger et al. (2018) show that the sign of the observed $M_2$ trend is correctly reproduced at 80% of the tide gauges on a global scale, but their simulated trends tend to differ from observations by a factor of 3 to 5; i.e. their simulations underestimate the $M_2$ response to MSL rise in terms of magnitude. Schindelegger et al. (2018) conclude that “magnitudes of observed and modeled $M_2$ trends are within a factor of 4 (or less) from each other in nearly 50% of the considered cases”. The large discrepancies between the simulations and the observations strongly suggest that MSL rise is not the only process that may explain $M_2$ changes – other large-scale processes, in addition to local processes, may also play a role.

Figure 5 shows the annual MSL, after removing the average over the period 1910–2010 and filtering with 9-year time windows. The correlations between $M_2$ and MSL indicate that $M_2$ varies strongly with MSL (see Sect. 4.2). However, $M_2$ variations show some variability in the North-East Atlantic (Fig. 2a), which may not be explained with MSL rise alone.

## 4.2 Possible link with MSL and climates indices

Processes other than MSL rise may impact the tide (see Sect. 1), such as the atmospheric circulation and the ocean stratification. Ocean and atmosphere are fully coupled, and air-sea fluxes are responsible for the exchange of momentum, water (evaporation and precipitation budget) and heat at their interface. Among the wide range of possible interactions, two mechanisms have been explored for their ability
<figure><img src="image_6.png" /><figcaption>Figure 3. Estimated trends in $M_2$ amplitude over the period 1910–2018.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Figure 4. Estimated trends in $M_2$ amplitude over the period 1990–2018.</figcaption></figure>

to modify the tide: (1) the momentum flux (wind stress) and the gradient of sea level pressure which act on the barotropic tide and (2) the water and heat fluxes which induce changes in both temperature and salinity distribution in the ocean. The latter effect acts on the stratification, which in turn could impact the tide in two different ways. The first way is the internal tide generation which transfers energy from barotropic and baroclinic motion and modifies surface tidal expression (Colosi and Munk, 2006). However, in the present study, most of the observations come from coastal stations sheltered by wide continental shelves which dampen internal waves. More important is the second way: the stratification acts on the eddy viscosity profile by modifying current profiles and bottom drag over continental shelves, which in turn modifies the $M_2$ surface expression (Kang et al., 2002; Müller, 2012; Katavouta et al., 2016).

Here, we focus on the effect of the atmospheric circulation on the tide. We used pressure indices (NAO and AO) that are relevant to represent atmospheric circulation. The NAO index represents the difference of normalized sea level pressure between the Azores high pressure system and the Iceland low pressure one (Hurrell, 1995). It indicates the redistribution of atmospheric masses between the subtropical Atlantic and the Arctic (Hurrell and Deser, 2009). In the North-East Atlantic, the similarity between the variations of the low-frequency winter NAO index and those of $M_2$ (Fig. 6) suggests a possible impact of large-scale atmospheric circulation on the tide.

The NAO index varies from positive to negative phases. Filtering the interannual variability, the NAO index tends overall to decrease between 1910 and 1970, then increase until 1990, and once again decrease. In the same way, $M_2$ amplitude tends to decrease up to 1960, then increase until 1990, and once again decrease. These similar patterns raise a possible connection between NAO and $M_2$ variation, already mentioned by Müller (2011) on the basis of qualitative criteria. In the following, we provide quantitative insights into the possible influence of NAO.

We computed the correlations ($r$ value) between normalized $M_2$ and climate indices, NAO and AO (Fig. 7). $M_2$, NAO and AO are filtered using the same time window (9 years). The correlations are computed since 1910, to have similar periods for all the stations. The correlations are considered as significant only if the $p$ value is lower than 0.05 (95% significance level). (Note that other statistics to measure the degree of association between the $M_2$ and NAO (AO) quantities would be worth exploring, for instance, nonlinear association using Spearman's correlation coefficient. In this respect, our study should be regarded as a first step that identifies sites worth considering in future investigations, especially investigating causal relationships with physics-based modelling.) The results are the following: (1) for NAO, 14 stations out of 18 show significant correlation. Note that at Brest, the correlation is significant since 1910, but not since 1864 (the NAO index used in this study starts only in 1864). This can be ex-
<figure><img src="image_8.png" /><figcaption>Figure 5. Annual mean sea levels (MSLs), after removing the average over the period 1910–2010. MSL values are filtered using 9-year windows.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Figure 6. Low-frequency winter NAO and AO indices, obtained with a 9-year mean filter. Normalized annual $M_2$ amplitudes in the North-East Atlantic (from Fig. 2a) are also plotted in grey.</figcaption></figure>

plained by the $M_2$ larger amplitude over all the 19th century, which decreases between 1890 and 1910 (Fig. 2a), possibly due to harbour development and construction of dykes (see Sect. 3.1). (2) In the North-East Atlantic, all the stations are positively correlated with NAO. (3) The strongest correlations (i.e. greater than 0.5) are in the northern part of the North Atlantic, with strong positive correlations at Cuxhaven and Hoek van Holland and strong negative correlation at Halifax (−0.55). (4) For AO, we found similar, but overall larger, r values. This is not surprising as these two indices are closely related.

To go further in the relative contribution of MSL and NAO in $M_2$ variability, we fitted two linear regression models on $M_2$ variations. In the following, $M_2$, MSL and NAO are filtered over 9-year time windows and normalized. At all the stations, we fitted $M_2$ variations with a MSL linear regression model (model 1) and a MSL and NAO multiple linear regression model (model 2). Models 1 and 2 may be expressed as

$$ \text{Model 1} = \alpha_1 \text{MSL} \qquad (2) $$

$$ \text{Model 2} = \alpha \text{MSL} + \beta \text{NAO}. \qquad (3) $$

The correlations between $M_2$ and model 1 (MSL) and model 2 (NAO and MSL) are presented in Fig. 8. We checked if there was correlation between NAO and MSL at the stations: there is no correlation at six stations, and r value is between 0.2 and 0.6 at eight stations; see Fig. 8 and discussion below. The results are the following: (1) $M_2$ varies at first order with MSL (Fig. 8). (2) The introduction of the NAO (model 2) allows increasing the predictive performance of the model, beyond the inherent effect of adding an additional regression parameter. Indeed, on average, the Akaike information criterion (AIC) is 99.9 for model 2, instead of 112.7 for model 1. On average, the $r^2$ value is 0.67 for model 2 instead of 0.61 for model 1. At some stations, the increase
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Correlation (*r* value) since 1910 between *M*₂ and (a) North Atlantic Oscillation and (b) Arctic Oscillation. Black dots are stations with no significant correlation. *M*₂, NAO and AO are filtered using the same time window (9 years).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Variance explained (*r*² value) since 1910 between *M*₂ and NAO, *M*₂ and MSL, *M*₂ and fitted model αMSL + βNAO (model 2), NAO and MSL. *M*₂, NAO and MSL are filtered using the same time window (9 years). Note that there is no orange bar for NAO-MSL when the correlation is not significant (*p* > 0.05).</figcaption></figure>

is quite large. For example at Cuxhaven, the *r*² value jumps from 0.42 to 0.64 between model 1 and 2. (3) The ratio β/(α+β) represents roughly the relative contribution of the NAO compared to the total effect of MSL and NAO (Fig. 9), as MSL and NAO are normalized. We found a significant contribution at some stations (e.g. more than 30% at Cuxhaven and Halifax), whereas it is negligible at others (e.g. only 5% at Portland). A total of 8 stations out of 18 show large NAO contribution (> 20%). The North-East Atlantic seems to be more sensitive to the NAO. Note that the interpretation of the results is tricky when MSL-NAO correlation is significant

(orange bars in Fig. 8). For example, at Hoek van Holland, the relative NAO contribution is very small, mainly because MSL and NAO are highly correlated (*r* = 0.59). Figure 10 shows *M*₂ variations along with the predictions from the two models, at all four stations where the NAO contribution is significant (β/(α+β) > 0.25), and the correlation between *M*₂ and model 2 is large enough (*r* > 0.3). At Cuxhaven, Halifax and Key West, model 2 (MSL- and NAO-dependent) naturally captures the *M*₂ variations better than model 1 (MSL-dependent); at Brest, the improvement is less significant. The trend switch observed since 1990 in the North-East Atlantic could be partly explained by the influence of the NAO on the tide.

These results suggest that a NAO-related mechanism may explain part of the variability of $M_2$. As mentioned by Müller (2011), “it is shown that sea-level, sea surface temperature and Arctic ice thickness are correlated with the NAO index. Thus, changes in the dynamics of the atmosphere could affect both $M_2$ and $S_2$ tides by processes discussed under (1), (2) and (3).” An underlying mechanism linked with (2) – sea surface temperature – could be changes in the ocean stratification. This is one of main possible hypotheses invoked in Ray and Talke (2019) to explain secular changes in $M_2$ amplitude in the Gulf of Maine; this is also the main hypothesis in Müller et al. (2014) and Gräwe et al. (2014) to explain seasonal modulation of $M_2$ in the North Sea. The relationship between the NAO index and stratification is complex and spatially variable across the North-East Atlantic (Fromentin and Planque, 1996). In the North Sea, the sea surface temperatures are positively correlated with NAO (Becker and Pauly, 1996), while subsurface temperatures show no significant correlation with NAO (Tian et al., 2016). Stratification
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Relative contribution of $\alpha$ compared to $\beta$ in the fitted model $\alpha\text{MSL} + \beta\text{NAO}$. Black dots are stations with no significant $M_2$-NAO correlation. The size of each large dot is proportional to the correlation between $M_2$ and the fitted model. Stations with no MSL-NAO correlations are labelled in bold.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Variations since 1910 of $M_2$, $\alpha\text{MSL}$ (model 1), $\alpha\text{MSL}+\beta\text{NAO}$ (model 2). $M_2$, NAO and MSL are filtered using the same time window (9 years).</figcaption></figure>

could therefore be (positively) correlated with NAO, but a dedicated study, outside the scope of this paper, would be necessary. Another underlying mechanism linked this time with (1) – sea level – could be the difference of spatial distribution of water level, due to different sea-level pressure and wind stress patterns. This is the hypothesis invoked in Huess and Andersen (2011) to explain the seasonal modulation of $M_2$ in the North Sea. They ran a barotropic model, forced with tides only and with both tides and meteorological fields; their results show that the $M_2$ seasonal modulation is better captured when the model is forced with both

tides and meteorological fields rather than with tides only. Figure 11a shows the average sea-level pressure during the period 1850–2015, derived from the Twentieth Century Reanalysis (20CR) (Compo et al., 2011; Slivinski et al., 2019). A positive NAO winter (e.g. 1989) corresponds to a situation with a stronger pressure gradient than average, between the two pressure systems of Azores and Iceland (Fig. 11c). By contrast, a negative NAO winter (e.g. 1969) corresponds to a weaker gradient pressure than usual (Fig. 11b). (We define winter here as December–February.) This way, from one year to another, the large-scale atmospheric masses are dis-
tributed differently, and as a consequence, the water volumes
are also distributed differently in the North Atlantic. In a sit-
uation of NAO$^{+}$, the surface waters are pushed onshore by
westerly winds, moving from Iceland to the European coasts
of France, Spain and Portugal. Figure 12 shows the redistrib-
bution of the sea-level pressure, between two years with high
and low NAO indices (here 1989 and 1969). Note that this is
an extreme situation, as these years have strong positive and
negative indices. Assuming an inverse barometer response
of sea level, the changes in terms of water level may vary
from more than 24 cm in the northwestern part of the area to
around –12 cm in the region that includes most of the North-
East Atlantic tide gauges considered in this study. This vari-
ation of a few tens of centimetres is probably negligible off-
shore but may have some impact on tide propagation along
the continental shelves and in shallow waters. It could also
shift slightly the amphidromic points. Assuming that these
changes have a similar impact (in terms of magnitude) on
M₂ as MSL changes, that is, ±10% in shallow waters ac-
cording to recent simulations (Pickering et al., 2017; Idier
et al., 2017), we find that they can yield centimetric changes
in M₂ amplitude. In other words, their order of magnitude is
roughly in agreement with the changes observed in M₂ (Ta-
ble 1). However, it is difficult to disentangle the effects of
stratification and meteorological forcing (sea-level pressure
and wind stress) in M₂ changes, and possibly both mecha-
nisms coexist. Dedicated simulations should be conducted to
assess the effects of atmospheric forcing on M₂ variability.

5 Conclusions

We investigated the long-term changes of the principal tidal
component *M*₂ over the North Atlantic coasts. We analysed
18 tide gauges with time series starting no later than 1940.
The longest is Brest with 165 years of data. We carefully pro-
cessed the data, particularly to remove the 18.6-year nodal
modulation.

We found that M₂ variations were consistent at all the sta-
tions in the North-East Atlantic (Cuxhaven, Delfzijl, Hoek
van Holland, Newlyn, Brest), whereas variations appear be-
tween stations in the North-West Atlantic. The changes
started long before the 20th century and are not linear. The
trends vary significantly from one station to another; they
are overall positive, up to 2.5 mm/yr, or slightly negative.
Since 1990, in many stations, the trends switch from posi-
tive to negative values. The significant differences between
the trends since 1910 and 1990 indicate caution when inter-
preting trends based on short records, i.e. less than 30 years,
especially if the data are noisy (interannual variability) and
the underlying processes non-linear (change points).

Concerning the causes of the observed changes, M₂ varies
primarily with the MSL, but MSL rise is not sufficient to ex-
plain the variations alone. The similarity between the North
Atlantic Oscillation and M₂ variations in the North-East At-

lantic suggests a possible influence of the large-scale atmo-
spheric circulation on the tide. Our statistical analysis con-
firms large correlations at all the stations in the North-East
Atlantic. The trend switch observed since 1990 could be
the signature of the large-scale atmospheric circulation on
the M2 tide. The underlying mechanism would be a differ-
ent spatial distribution of water level from one year to an-
other, depending on the low-frequency sea-level pressure pat-
terns, and impacting the propagation of the tide in the North
Atlantic basin. In the future, dedicated modelling studies
should be undertaken to confirm or discard this hypothesis.
These simulations should also allow estimating the effect of
the wind (through the Ekman current) and currents on M2
changes (Devlin et al., 2018).

In this study, we focused only on M₂ amplitude. A similar
analysis on the phase lag would draw a more complete pic-
ture of the M₂ variations (Müller, 2011; Woodworth, 2010;
Ray and Talke, 2019). Other constituents are also affected.
Results show that S₂ amplitude decreases at all the stations
located in the North-West Atlantic and, in contrast, tends to
increase in the North-East Atlantic (not shown). The large-
scale decrease of S₂ observed in the North-West Atlantic is
consistent with previous studies (e.g. Ray, 2006, in the Gulf
of Maine). Further investigations should be definitely con-
ducted to extend this work to more constituents.

The historic data show that the changes started long before
the 20th century. This conclusion would not have been pos-
sible without the huge work of data rescue undertaken over
the past decades (e.g. Pouvreau et al., 2006; Pouvreau, 2008;
Bradshaw et al., 2016). This underlines the great importance
of sea level data archaeology, which allows extending and
improving historical datasets (Pouvreau, 2008; Woodworth
et al., 2010; Marcos et al., 2011; Talke and Jay, 2013, 2017;
Ray and Talke, 2019; Bradshaw et al., 2015, 2020; Haigh
et al., 2019). This is essential for studies related to climate
change.

Finally, we should mention several additional limitations
and perspectives in this study. (1) We processed the time se-
ries considering that they were quality controlled. A fuller
analysis of the data quality before processing would prob-
ably be valuable. (2) We did not investigate the history of
each station. There are probably some local changes (e.g.
environment or instrumentation) that may explain a part of
the variability of M₂ amplitude and some discrepancies be-
tween stations. (3) The tide gauges are located mainly in
harbours. They are affected at the same time by local- and
regional/global-scale changes, which are difficult to sepa-
rate. Moreover, they may not be representative of changes
offshore. A similar study based on satellite altimetry data
would probably be of great interest, even if temporal scale
for satellite data is still rather short (i.e. < 30 years) com-
pared to climate-scale processes. (4) We focused mainly on
the UHSLC dataset, which consists of 249 stations in the At-
lantic Ocean. Other relevant stations (which are not in this
dataset) may be considered in future studies. (5) We did not
<figure><img src="image_14.png" /><figcaption>**Figure 11.** Winter sea-level pressure over the North-East Atlantic (a) average over 1850–2015 (b) anomaly in 1969 (NAO⁻) (c) anomaly in 1989 (NAO⁺). Contour intervals are every 2 hPa.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Figure 12.** Difference of winter sea-level pressure between 1989 (NAO⁺) and 1969 (NAO⁻) over the North-East Atlantic. Contour intervals are every 4 hPa.</figcaption></figure>

investigate the impact of storminess on the tide. Dedicated studies are necessary to estimate if changes in storminess could affect significantly tidal constituents. (6) We used only winter AO and NAO indices, which show more variability than annual indices. A similar analysis with annual indices shows similar results for the correlation with AO or NAO (positive correlation on the North-East Atlantic). With annual rather than monthly indices, the difference of pressure fields will decrease, and as a consequence, the magnitude of the sea-level response will also decrease. Further investigations should be conducted on this point.
## Appendix A: Nodal modulation

The $M_2$ component is subject to a 18.6-year modulation, separated from a neighbouring line in the tidal potential ($m_2$) whose Doodson number differs in its fifth frequency (255 555 and 255 545 for $M_2$ and $m_2$, respectively) (Doodson and Warburg, 1941; Pugh and Woodworth, 2014). This fifth frequency corresponds to $N'$, the negative of the mean longitude of the Moon ascending node – hence the “nodal” term – whose period is 18.6 years. Note that there is also another component close to $M_2$, whose Doodson number differs only from the fifth frequency (255 565), but it is negligible, its amplitude in the tidal potential being only 0.05% of $M_2$, whereas $m_2$ amplitude is 3.7% of $M_2$ (Simon, 2007, 2013). With one year of hourly data, the two components $M_2$ and $m_2$ cannot be separated by a yearly harmonic analysis (at least 18.6 years are necessary). As a consequence, $M_2$ amplitude is modulated by $m_2$. However, we can estimate this modulation and remove it. The harmonic formulation is expressed schematically as a sum of harmonic components:

$$h(t) = \sum_{i} a_{i} \cos(V_{i}(t) - \kappa_{i}), \qquad (A1)$$

where $h(t)$ is the sea level height at time $t$, $V_i(t)$ is the astronomical argument (computed from Doodson number) and $a_i$, $\kappa_i$ the amplitude and phase lag of each component. Considering that $M_2$ and $m_2$ are very close in terms of frequency, we can assume that their phase lags are similar ($\kappa_{M_2} \simeq \kappa_{m_2}$). As their difference of astronomical arguments is $V_{m_2} - V_{M_2} = N' + \pi$, the $M_2$ and $m_2$ contributions to the total water level may be expressed as

$$h_{M_2}(t) + h_{m_2}(t) = h_{M_2}(t)[1 + f_{\text{nod}} \cos(N' + \pi)], \qquad (A2)$$

where $f_{\text{nod}}$, the nodal modulation, is the ratio of the amplitude of $m_2$ and $M_2$. As $M_2$ and $m_2$ are very close in terms of frequency, $f_{\text{nod}}$ is generally considered as close to the ratio of their amplitude in the tidal potential, $A_{m_2}$ and $A_{M_2}$:

$$f_{\text{nod}} = \frac{a_{m_2}}{a_{M_2}} \approx \frac{A_{m_2}}{A_{M_2}} \approx 0.037. \qquad (A3)$$

<figure><img src="image_16.png" /><figcaption>Figure A1. (a) Estimation of the nodal modulation of $M_2$ amplitude (mean removed) at Newlyn. (b) Impact on $M_2$ amplitude of the nodal modulation correction at Newlyn. $M_2$ is detrended in (a) to better fit the nodal modulation.</figcaption></figure>

The negative of the mean longitude of the Moon ascending node is expressed simply as a function of time (p. 116 in Simon, 2007, p. 112 in Simon, 2013):

$$N' = -N = 234.555 + 1934.1363T + 0.0021T^2, \qquad (A4)$$

with $N'$ in degrees, and $T$ the time elapsed since 1 January 2000 at 12:00, expressed in Julian centuries (36 525 d).

The tidal program we used (MAS) corrected $M_2$ applying the usual 3.7% nodal modulation (Eq. A3). However, this value may vary significantly from one station to another; Ray (2006) reported values ranging from 2.3% to 3.6% in the Gulf of Maine. Here, we computed directly $f_{\text{nod}}$ from the observed data, proceeding as follows. (1) We added the default nodal correction $1+0.037\cos(N'+\pi)$ to the $M_2$ variations. (2) We detrended the obtained signal removing the last intrinsic mode function (IMF) of an empirical mode decomposition (EMD) (Huang et al., 1998); note that the EMD is an analysis tool which partitions a series into “modes” (i.e. IMFs), the last one being the trend of the signal. (3) We fitted a function $a_{m_2} \cos(N'+\pi)$ to this detrended signal to estimate $a_{m_2}$, $N'$ being expressed as in Eq. (A4). (4) We finally computed $f_{\text{nod}}$ as the ratio between $m_2$ and $M_2$ amplitudes (Eq. A3). Figure A1a shows an example of estimate of $M_2$ modulation at Newlyn: the fit leads to a nodal modulation of 3.3%. Note that this value is consistent with Woodworth (2010) (3.2%), whereas Woodworth et al. (1991) gave a slightly different value (2.8%). Figure A1b shows the impact of this value rather than the default one: oscillations of 18.6 years are clearly reduced. Note that in this study, the $m_2$ amplitude – and then the nodal correction – could have been computed from the full time series harmonic analysis, as records are longer than 18.6 years. However, the method presented here to compute the nodal correction can be applied even for time series shorter than 18.6 years.

The computed nodal modulations are summarized in Table 1 (column 6). They vary from 0.8% to 4.1%. Note that these values are consistent with those obtained by previous authors (Ray, 2006; Müller, 2011; Woodworth, 2010; Ray and Talke, 2019). Only the value at Charleston differs significantly: 3.0% in our study compared to 3.7% in Müller (2011).
<figure><img src="image_17.png" /><figcaption>**Figure B1.** Annual *M*₂ amplitude at the 18 selected tide gauges.</figcaption></figure>
*Data availability.* The tide gauge data are available at University of Hawaii Sea Level Center ftp://ftp.soest.hawaii.edu/uhslc/rqds atlantic (Caldwell et al., 2015). The NAO climate index is available at https://climatedataguide.ucar.edu/climate-data/hurrell-north-atlantic-oscillation-nao-index-station-based (Hurrell and the National Center for Atmospheric Research Staff, 2020). The AO climate index is available at https://climatedataguide.ucar.edu/climate-data/hurrell-wintertime-slp-based-northern-annular-mode-nam-index (National Center for Atmospheric Research Staff, 2020). The Twentieth Century Reanalysis Project version 3 dataset is available at https://www.psl.noaa.gov/data/gridded/data.20thC_ReanV3.monolevel.html#caveat (Slivinski et al., 2019).

*Author contributions.* LPG analysed the data and wrote the paper. PL contributed to the interpretation of the data and the writing of the paper. GW contributed to the analysis and interpretation of the data and the writing of the paper.

*Competing interests.* The authors declare that they have no conflict of interest.

*Special issue statement.* This article is part of the special issue “Developments in the science and history of tides (OS/ACP/HGSS/NPG/SE inter-journal SI)”. It is not associated with a conference.

*Acknowledgements.* The sea level observations were provided by the University of Hawaii Sea Level Center – retrieved from ftp://ftp.soest.hawaii.edu/uhslc/rqds (last access: April 2020). The sea level data at Delfzijl and Hoek van Holland were provided by Rijkswaterstaat (RWS) Service Desk, Netherlands. The climate indices (NAO and AO indices) were provided by the Climate Analysis Section, NCAR, Boulder, USA – retrieved from https://climatedataguide.ucar.edu/climate-data/ (last access: April 2020). The harmonic analysis program MAS was provided by the French Hydrographic Office (SHOM). Support for the Twentieth Century Reanalysis Project version 3 dataset was provided by the US Department of Energy, Office of Science Biological and Environmental Research (BER), by the National Oceanic and Atmospheric Administration Climate Program Office, and by the NOAA Physical Sciences Laboratory. The authors very warmly thank the two reviewers (Stefan Talke and an anonymous reviewer) and the editor (Philip Woodworth) for their careful reading and their many constructive comments, which greatly improved the paper.

*Financial support.* This research has been supported by the French Research Institute for Exploitation of the Sea (IFREMER) and by the research theme “Long-term observing systems for ocean knowledge” of the ISblue project “Interdisciplinary graduate school for the blue planet”, co-funded by a grant from the French government under the program “Investissements d’Avenir” (ANR-17-EURE-0015).

*Review statement.* This paper was edited by Philip Woodworth and reviewed by Stefan Talke and one anonymous referee.

## References

Araújo, I. B. and Pugh, D. T.: Sea levels at Newlyn 1915–2005: Analysis of trends for future flooding risks, J. Coast. Res., 24, 203–212, https://doi.org/10.2112/06-0785.1, 2008.

Becker, G. A. and Pauly, M.: Sea surface temperature changes in the North Sea and their causes, ICES J. Mar. Sci., 53, 887–898, https://doi.org/10.1006/jmsc.1996.0111, 1996.

Bradshaw, E., Rickards, L., and Aarup, T.: Sea level data archaeology and the Global Sea Level Observing System (GLOSS), Georesj., 6, 9–16, https://doi.org/10.1016/j.grj.2015.02.005, 2015.

Bradshaw, E., Woodworth, P., Hibbert, A., Bradley, L., Pugh, D., Fane, C., and Bingley, R.: A century of sea level measurements at Newlyn, Southwest England, Mar. Geodesy, 39, 115–140, https://doi.org/10.1080/01490419.2015.1121175, 2016.

Bradshaw, E., Ferret, Y., Pons, F., Testut, L., and Woodworth, P.: Workshop on sea level data archaeology, Technical Report, Workshop Report No. 287, Intergovernmental Oceanographic Commission, Paris, France, 47 pp., 2020.

Caldwell, P. C., Merrifield, M. A., and Thompson, P. R.: Sea level measured by tide gauges from global oceans – the Joint Archive for Sea Level holdings (NCEI Accession 0019568), Version 5.5, NOAA National Centers for Environmental Information, Dataset, ftp://ftp.soest.hawaii.edu/uhslc/rqds atlantic, 2015.

Cartwright, D. E.: Secular changes in the oceanic tides at Brest, 1711–1936, Rev. Geophys., 57, 433–449, https://doi.org/10.1111/j.1365-246X.1972.tb05826.x, 1972.

Colosi, J. A. and Munk, W.: Tales of the venerable Honolulu tide gauge, J. Phys. Oceanogr., 36, 967–996, https://doi.org/10.1175/JPO2876.1, 2006.

Compo, G. P., Whitaker, J. S., Sardeshmukh, P. D., Matsui, N., Alvan, R. J., Yin, X., Gleason, B. E., Vose, R. S., Rutledge, G., Bessemoulin, P., Brönnimann, S., Brunet, M., Crouthamel, R. I., Grant, A. N., Groisman, P. Y., Jones, P. D., Kruk, M., Kruger, A. C., Marshall, G. J., Maugeri, M., Mok, H. Y., Nordli, Ø., Ross, T. F., Trigo, R. M., Wang, X. L., Woodruff, S. D., and Worley, S. J.: The twentieth century reanalysis project, Q. J. Roy. Meteorol. Soc., 137, 1–28, https://doi.org/10.1002/qj.776, 2011.

Devlin, A. T., Zaron, E. D., Jay, D. A., Talke, S. A., and Pan, J.: Seasonality of tides in Southeast Asian waters, J. Phys. Oceanogr., 48, 1169–1190, https://doi.org/10.1175/JPO-D-17-0119.1, 2018.

Doodson, A. T.: Perturbations of harmonic tidal constants, P. R. Soc. A, 106, 513–526, https://doi.org/10.1098/rspa.1924.0085, 1924.

Doodson, A. T. and Warburg, H. D.: Admiralty manual of tides, HMSO, London, UK, 1941.

Familkhalili, R. and Talke, S. A.: The effect of channel deepening on tides and storm surge: A case study of Wilmington, NC, Geophys. Res. Lett., 43, 9138–9147, https://doi.org/10.1002/2016GL069494, 2016.

Fromentin, J.-M. and Planque, B.: Calanus and environment in the eastern North Atlantic. II. Influence of the North Atlantic Oscillation on C. finmarchicus and C. helgolandicus, Mar. Ecol. Prog. Ser., 134, 111–118, https://doi.org/10.3354/meps13411, 1996.

Garrett, C.: Tidal resonance in the Bay of Fundy and Gulf of Maine, Nature, 238, 441–443, https://doi.org/10.1038/238441a0, 1972.
Godin, G.: On tidal resonance, Cont. Shelf Res., 13, 89–107, https://doi.org/10.1016/0278-4343(93)90037-X, 1993.

Godin, G.: Rapid evolution of the tide in the Bay of Fundy, Cont. Shelf Res., 15, 369–372, https://doi.org/10.1016/0278-4343(93)E0005-S, 1995.

Gräwe, U., Burchard, H., Müller, M., and Schuttelaars, H. M.: Seasonal variability in M₂ and M₄ tidal constituents and its implications for the coastal residual sediment transport, Geophys. Res. Lett., 41, 5563–5570, https://doi.org/10.1002/2014GL060517, 2014.

Haigh, I. D., Pickering, M. D., Green, J. A. M., Arbic, B. K., Arns, A., Dangendorf, S., Hill, D. F., Horsburgh, K., Howard, T., Idier, D., Jay, D. A., Jänicke, L., Lee, S. B., Müller, M., Schindelegger, M., Talke, S. A., Wilmes, S., and Woodworth, P. L.: The tides they are a-changin': A comprehensive review of past and future nonastronomical changes in tides, their driving mechanisms and future implications, Rev. Geophys., 57, e2018RG000636, https://doi.org/10.1029/2018RG000636, 2019.

Hollebrandse, F.: Temporal development of the tidal range in the southern North Sea, Master thesis, Faculty of Civil Engineering and Geosciences, Delft University of Technology, Delft, The Netherlands, 95 pp., 2005.

Huang, N. E., Shen, Z., Long, S. R., Wu, M. C., Shih, H. H., Zheng, Q., Yen, N.-C., Tung, C. C., and Liu, H. H.: The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis, P. Roy. Soc. Lond. A, 454, 903–995, https://doi.org/10.1098/rspa.1998.0193, 1998.

Huess, V. and Andersen, O. B.: Seasonal variation in the main tidal constituent from altimetry, Geophys. Res. Lett., 28, 567–570, https://doi.org/10.1029/2000GL011921, 2001.

Hurrell, J.: Decadal trends in the North Atlantic Oscillation, regional temperatures and precipitation, Science, 269, 676–679, https://doi.org/10.1126/science.269.5224.676, 1995.

Hurrell, J. and the National Center for Atmospheric Research Staff (Eds): The Climate Data Guide: Hurrell North Atlantic Oscillation (NAO) Index (station-based), available at: https://climatedataguide.ucar.edu/climate-data/hurrell-north-atlantic-oscillation-nao-index-station-based, last modified at: 24 April 2020.

Hurrell, J., Kushnir, Y., Ottersen, G., and Visbeck, M.: An overview of the North Atlantic Oscillation, in: The North Atlantic Oscillation: Climatic Significance and Environmental Impact, edited by: Hurrell, J. W., Kushnir, Y., Ottersen, G., Visbeck, M., and Visbeck, M. H., American Geophysical Union, Washington D.C., USA, 1–36, https://doi.org/10.1029/134GM01, 2003.

Hurrell, J. W. and Deser, C.: North Atlantic climate variability: The role of the North Atlantic Oscillation, J. Mar. Syst., 78, 28–41, https://doi.org/10.1016/j.jmarsys.2008.11.026, 2009.

Idier, D., Paris, F., Cozannet, G. L., Boulahyaa, F., and Dumas, F.: Sea-level rise impacts on the tides of the European Shelf, Cont. Shelf Res., 137, 56–71, https://doi.org/10.1016/j.csr.2017.01.007, 2017.

Jacob, B., Stanev, E. V., and Zhang, Y. J.: Local and remote response of the North Sea dynamics to morphodynamic changes in the Wadden Sea, Ocean Dyn., 66, 671–690, https://doi.org/10.1007/s10236-016-0949-8, 2016.

Kang, S. K., Foreman, M. G. G., Lie, H. J., Lee, J. H., Cherniawsky, J., and Yum, K. D.: Two-layer tidal modeling of the Yellow and East China Seas with application to seasonal vari-

ability of the M₂ tide, J. Geophys. Res.-Oceans, 107, 3020, https://doi.org/10.1029/2001JC000838, 2002.

Katavouta, A., Thompson, K. R., Lu, Y., and Loder, J. W.: Interaction between the tidal and seasonal variability of the Gulf of Maine and Scotian shelf region, J. Phys. Oceanogr., 46, 3279–3298, https://doi.org/10.1175/JPO-D-15-0091.1, 2016.

Marcos, M., Puyol, B., Wöppelmann, G., Herrero, C., and García-Fernández, M. J.: The long sea level record at Cadiz (southern Spain) from 1880 to 2009, J. Geophys. Res.-Oceans, 116, C12003, https://doi.org/10.1029/2011JC007558, 2011.

Mawdsley, R. J., Haigh, I. D., and Wells, N. C.: Global secular changes in different tidal high water, low water and range levels, Earths Future, 3, 66–81, https://doi.org/10.1002/2014EF000282, 2015.

Mudersbach, C., Wahl, T., Haigh, I. D., and Jensen, J.: Trends in high sea levels of German North Sea gauges compared to regional mean sea level change, Cont. Shelf Res., 65, 111–120, https://doi.org/10.1016/j.csr.2013.06.016, 2013.

Müller, M.: Rapid change in semi-diurnal tides in the North Atlantic since 1980, Geophys. Res. Lett., 38, L11602, https://doi.org/10.1029/2011GL047312, 2011.

Müller, M.: The influence of changing stratification conditions on barotropic tidal transport and its implications for seasonal and secular changes of tides, Cont. Shelf Res., 47, 107–118, https://doi.org/10.1016/j.csr.2012.07.003, 2012.

Müller, M., Arbic, B. K., and Mitrovica, J. X.: Secular trends in ocean tides: Observations and model results, J. Geophys. Res.-Oceans, 116, C05013, https://doi.org/10.1029/2010JC006387, 2011.

Müller, M., Cherniawsky, J. Y., Foreman, M. G. G., and von Storch, J. S.: Seasonal variation of the M₂ tide, Ocean Dyn., 64, 159–177, https://doi.org/10.1007/s10236-013-0679-0, 2014.

National Center for Atmospheric Research Staff (Eds): The Climate Data Guide: Hurrell wintertime SLP-based Northern Annular Mode (NAM) Index, available at: https://climatedataguide.ucar.edu/climate-data/hurrell-wintertime-slp-based-northern-annular-mode-nam-index, last modified at: 22 May 2020.

Pan, H., Zheng, Q., and Lv, X.: Temporal changes in the response of the nodal modulation of the M₂ tide in the Gulf of Maine, Cont. Shelf Res., 186, 13–20, https://doi.org/10.1016/j.csr.2019.07.007, 2019.

Pawlowicz, R., Beardsley, B., and Lentz, S.: Classical tidal harmonic analysis including error estimates in MAT-LAB using T_TIDE, Comput. Geosci., 28, 929–937, https://doi.org/10.1016/S0098-3004(02)00013-4, 2002.

Pickering, M., Wells, N. C., Horsburgh, K. J., and Green, J. A. M.: The impact of future sea-level rise on the European Shelf tides, Cont. Shelf Res., 35, 1–15, https://doi.org/10.1016/j.csr.2011.11.011, 2012.

Pickering, M. D., Horsburgh, K. J., Blundell, J. R., Hirschi, J. J.-M., Nicholls, R. J., Verlaan, M., and Wells, N.: The impact of future sea-level rise on the global tides, Cont. Shelf Res., 142, 50–68, https://doi.org/10.1016/j.csr.2017.02.004, 2017.

Pouvreau, N.: Trois cents ans de mesures marégraphiques en France: outils, méthodes et tendances des composantes du niveau de la mer au port de Brest, PhD thesis, University of La Rochelle, La Rochelle, France, 475 pp., 2008.
Pouvreau, N., Miguez, B. M., Simon, B., and Wöppelmann, G.: Évolution de l'onde semi-diurne M₂ de la marée à Brest de 1846 à 2005, C. R. Geosci., 338, 802–808, https://doi.org/10.1016/j.crte.2006.07.003, 2006.

Pugh, D. and Woodworth, P.: Sea-Level Science: Understanding tides, surges, tsunamis and mean sea-level changes, Cambridge University Press, Cambridge, UK, 2014.

Pugh, D. T. and Vassie, J. M.: Tide and surge propagation offshore in the Dowsing region of the North Sea, Dtsch. Hydrogr. Zeitschrift, 29, 163–213, https://doi.org/10.1007/BF02226659, 1976.

Ralston, D. K., Talke, S., Geyer, W. R., Al-Zubaidi, H. A. M., and Sommerfield, C. K.: Bigger tides, less flooding: Effects of dredging on barotropic dynamics in a highly modified estuary, J. Geophys. Res.-Oceans, 124, 196–211, https://doi.org/10.1029/2018JC014313, 2019.

Ray, R. D.: Secular changes of the M₂ tide in the Gulf of Maine, Cont. Shelf Res., 26, 422–427, https://doi.org/10.1016/j.csr.2005.12.005, 2006.

Ray, R. D.: Secular changes in the solar semidiurnal tide of the western North Atlantic ocean, Geophys. Res. Lett., 36, L19601, https://doi.org/10.1029/2009GL040217, 2009.

Ray, R. D. and Talke, S. A.: Nineteenth-century tides in the Gulf of Maine and implications for secular trends, J. Geophys. Res.-Oceans, 124, 7046–7067, https://doi.org/10.1029/2019JC015277, 2019.

Santamaría-Gómez, A.: SARI: interactive GNSS position time series analysis software, GPS Solut., 23, 52, https://doi.org/10.1007/s10291-019-0846-y, 2019.

Schindelegger, M., Green, J. A. M., Wilmes, S., and Haigh, I. D.: Can we model the effect of observed sea level rise on tides?, J. Geophys. Res.-Oceans, 123, 4593–4609, https://doi.org/10.1029/2018JC013959, 2018.

Simon, B.: La marée océanique côtière, Institut Océanographique Ed., Paris, France, available at: https://iho.int/iho_pubs/CB/C-33/C-33_maree_simon_fr.pdf (last access: May 2020), 2007.

Simon, B.: Coastal tides, Institut Océanographique Ed., Paris, France, available at: https://iho.int/iho_pubs/CB/C-33/C-33_maree_simon_en.pdf (last access: May 2020), 2013.

Slivinski, L. C., Compo, G. P., Whitaker, J. S., Sardeshmukh, P. D., Giese, B. S., McColl, C., Allan, R., Yin, X., Vose, R., Titchner, H., Kennedy, J., Spencer, L. J., Ashcroft, L., Brönnimann, S., Brunet, M., Camuffo, D., Cornes, R., Cram, T. A., Crouthamel, R., Domínguez-Castro, F., Freeman, J. E., Gergis, J., Hawkins, E., Jones, P. D., Jourdain, S., Kaplan, A., Kubota, H., Blancq, F. L., Lee, T., Lorrey, A., Luterbacher, J., Maugeri, M., Mock, C. J., Moore, G. K., Przybylak, R., Pudmenzky, C., Reason, C., Slonosky, V. C., Smith, C., Tinz, B., Trewin, B., Valente, M. A., Wang, X. L., Wilkinson, C., Wood, K., and Wyszyński, P.: Towards a more reliable historical reanalysis: Improvements for version 3 of the Twentieth Century Reanalysis system, Q. J. Roy. Meteorol. Soc., 145, 2876–2908, https://doi.org/10.1002/qj.3598, 2019 (data available at: https://www.psl.noaa.gov/data/gridded/data.20thC_ReanV3.monolevel.html#caveat, last access: July 2020).

Talke, S. A. and Jay, D. A.: Nineteenth century North American and Pacific tidal data: lost or just forgotten?, J. Coast. Res., 29, 118–127, https://doi.org/10.2112/JCOASTRES-D-12-00181.1, 2013.

Talke, S. A. and Jay, D. A.: Archival water-level measurements: Recovering historical data to help design for the future, Technical Report, Civil and Environmental Engineering Faculty Publications and Presentations, Portland State University, Portland, OR, USA, 50 pp., 2017.

Talke, S. A. and Jay, D. A.: Changing tides: The role of natural and anthropogenic factors, Ann. Rev. Mar. Sci., 12, 121–151, https://doi.org/10.1146/annurev-marine-010419-010727, 2020.

Talke, S. A., Orton, P., and Jay, D. A.: Increasing storm tides in New York harbor, 1844–2013, Geophys. Res. Lett., 41, 3149–3155, https://doi.org/10.1002/2014GL059574, 2014.

Talke, S. A., Kemp, A., and Woodruff, J.: Relative sea level, tides, and extreme water levels in Boston (MA) from 1825 to 2018, J. Geophys. Res.-Oceans, 123, 3895–3914, https://doi.org/10.1029/2017JC013645, 2018.

Thompson, D. W. J. and Wallace, J. M.: The Arctic oscillation signature in wintertime geopotential height and temperature fields, Geophys. Res. Lett., 25, 1297–1300, https://doi.org/10.1029/98GL00950, 1998.

Thompson, D. W. J. and Wallace, J. M.: Annular modes in the extratropical circulation. Part I: Month-to-month variability, J. Climate, 13, 1000–1016, https://doi.org/10.1175/1520-0442(2000)013<1000:AMITEC>2.0.CO;2, 2000.

Thompson, D. W. J., Wallace, J. M., and Hegerl, G. C.: Annular modes in the extratropical circulation. Part II: Trends, J. Climate, 13, 1018–1036, https://doi.org/10.1175/1520-0442(2000)013<1018:AMITEC>2.0.CO;2, 2000.

Tian, T., Su, J., Boberg, F., Yang, S., and Schmith, T.: Estimating uncertainty caused by ocean heat transport to the North Sea: experiments downscaling EC-Earth, Clim. Dynam., 46, 99–110, https://doi.org/10.1007/s00382-015-2571-8, 2016.

Wikipedia: Brest Arsenal, available at: https://en.wikipedia.org/wiki/Brest_Arsenal (last access: 2 November 2020), 2020.

Winterwerp, J. C. and Wang, Z. B.: Man-induced regime shifts in small estuaries – I: theory, Ocean Dyn., 63, 1279–1292, https://doi.org/10.1007/s10236-013-0662-9, 2013.

Winterwerp, J. C., Wang, Z. B., van Braeckel, A., van Holland, G., and Kösters, F.: Man-induced regime shifts in small estuaries – II: a comparison of rivers, Ocean Dyn., 63, 1293–1306, https://doi.org/10.1007/s10236-013-0663-8, 2013.

Woodworth, P.: A survey of recent changes in the main components of the ocean tide, Cont. Shelf Res., 30, 1680–1691, https://doi.org/10.1016/j.csr.2010.07.002, 2010.

Woodworth, P. L., Shaw, S. M., and Blackman, D. L.: Secular trends in mean tidal range around the British Isles and along the adjacent European coastline, Geophys. J. Int., 104, 593–609, https://doi.org/10.1111/j.1365-246X.1991.tb05704.x, 1991.

Woodworth, P. L., Pouvreau, N., and Wöppelmann, G.: The gyre-scale circulation of the North Atlantic and sea level at Brest, Ocean Sci., 6, 185–190, https://doi.org/10.5194/os-6-185-2010, 2010.

Wöppelmann, G., Pouvreau, N., and Simon, B.: Brest sea level record: a time series construction back to the early eighteenth century, Ocean Dyn., 56, 487–497, https://doi.org/10.1007/s10236-005-0044-z, 2006.