Observation of a Fast Continental Shelf Wave Generated by a Storm Impacting Newfoundland Using Wavelet and Cross-Wavelet Analyses

Severin Thiebaut Department of Marine Science, University of Otago, Dunedin, New Zealand

Search for other papers by Severin Thiebaut in
Current site
Google Scholar
PubMed
Close
and
Ross Vennell Department of Marine Science, University of Otago, Dunedin, New Zealand

Search for other papers by Ross Vennell in
Current site
Google Scholar
PubMed
Close
Full access

We are aware of a technical issue preventing figures and tables from showing in some newly published articles in the full-text HTML view.
While we are resolving the problem, please use the online PDF version of these articles to view figures and tables.

Abstract

Wavelet and cross-wavelet power spectra of sea level records from tide gauges along the Atlantic coast of Canada showed a low-frequency barotropic response after Hurricane Florence crossed the Newfoundland shelf in September 2006. In comparison with two other storms, the results showed that Florence was the only one that excited a propagating sea level disturbance with a period range similar to the passage time of the storm over the shelf (26–30 h) and phase shifts consistent with a barotropic continental shelf wave (CSW). The high amplitude of the oscillations generated by Florence along the shore diminished from approximately 45 to 12 cm as the CSW propagated from the south coast of Newfoundland to the southern Nova Scotia seaboard. This paper presents the first direct measurement of a remarkably high alongshore group speed (11.4 ± 5.9 m s−1), in the manner of free-barotropic CSW, by examination of sea level wavelet power spectra at different locations. Furthermore, using cross-wavelet analysis of pairs of stations, an exceptional phase speed of 16.0 ± 5.1 m s−1 has been found, greater than had been previously observed for a free CSW. The results were consistent with dispersion curves for the first-mode barotropic CSW.

Corresponding author address: Severin Thiebaut, Dept. of Marine Science, University of Otago, Dunedin 9012, New Zealand. Email: thise135@student.otago.ac.nz

Abstract

Wavelet and cross-wavelet power spectra of sea level records from tide gauges along the Atlantic coast of Canada showed a low-frequency barotropic response after Hurricane Florence crossed the Newfoundland shelf in September 2006. In comparison with two other storms, the results showed that Florence was the only one that excited a propagating sea level disturbance with a period range similar to the passage time of the storm over the shelf (26–30 h) and phase shifts consistent with a barotropic continental shelf wave (CSW). The high amplitude of the oscillations generated by Florence along the shore diminished from approximately 45 to 12 cm as the CSW propagated from the south coast of Newfoundland to the southern Nova Scotia seaboard. This paper presents the first direct measurement of a remarkably high alongshore group speed (11.4 ± 5.9 m s−1), in the manner of free-barotropic CSW, by examination of sea level wavelet power spectra at different locations. Furthermore, using cross-wavelet analysis of pairs of stations, an exceptional phase speed of 16.0 ± 5.1 m s−1 has been found, greater than had been previously observed for a free CSW. The results were consistent with dispersion curves for the first-mode barotropic CSW.

Corresponding author address: Severin Thiebaut, Dept. of Marine Science, University of Otago, Dunedin 9012, New Zealand. Email: thise135@student.otago.ac.nz

1. Introduction

A continental shelf wave (CSW), first observed along the Australian coast by Hamon (1962, 1966), is a type of topographic planetary wave trapped over the continental margin with a frequency less than the local inertial frequency and a wavelength much greater than the depth (LeBlond and Mysak 1978). It is generally accepted that these waves are generated by large-scale weather systems moving across or along the shelf (Robinson 1964; Mysak 1967a,b; Adams and Buchwald 1969; Gill and Schumann 1974; LeBlond and Mysak 1978). In this study, we focus on a storm traveling across the shelf, which might produce a standing wave resulting from resonant motion excited by the forcing, and a CSW freely propagating with the coastline on its right in the North Hemisphere (LeBlond and Mysak 1978; Slørdal et al. 1994; Tang et al. 1998). Wavelet analysis is shown to be able to estimate group speeds by calculating the distance traveled by a wave envelope, which contains waves of a selected frequency interval. For calculating phase speed, cross-wavelet analysis is particularly adapted for studies of nonstationary signals.

Barotropic shelf waves have been observed in many locations; typically, the propagating speeds were within the range of 1–6 m s−1 (Hamon 1966; Isozaki 1968; Mysak and Hamon 1969; Cutchin and Smith 1973; Brooks and Mooers 1977; Enfield and Allen 1983; Camayo and Campos 2006). Higher speeds of barotropic CSW have been found, but were identified as forced waves, whose speed is determined by the speed of the forcing agent (Enfield and Allen 1983; Cuevas et al. 1986; Yankovsky and Garvine 1998). Several numerical simulations of the response of shelf seas to weather forcing have been done in the past and gave results of free-CSW phase and/or group speed (Martinsen et al. 1979; Enfield and Allen 1983; Gordon and Huthnance 1987; Gjevik 1991; Sheng et al. 2006). Nevertheless, to our knowledge, none of them found fast group and phase speeds for CSW by examining the sea level oscillations after a strong storm event, using a wavelet approach.

In this paper we focus on Hurricane Florence, which hit Newfoundland on 13 September 2006 and generated a significant sea level disturbance propagating freely along the eastern Canadian coast. The group and phase speeds of the free-barotropic CSW are estimated using wavelet and cross-wavelet techniques. The speeds were higher than previous observations in many other parts of the world. Wavelet methods are necessarily technical and require detailed explanation. However, the main focuses of this work are the observed high-CSW group and phase speeds. Results are compared with dispersion curves for the first-mode barotropic CSW, using two simplified shelf profiles and with two other storms with different durations (the passage times of the storms over a particular point on the coast) and intensities. Following a summary of the data used (section 2), wavelet and cross-wavelet analyses are presented in section 3, and the results and the conditions for propagation of the CSW generated by the events are discussed in section 4.

2. Data

Three events were studied: Hurricane Florence (12–15 September 2006), Hurricane Isaac (1–3 October 2006), and Tropical Storm Chantal (31 July–2 August 2007). Note that the general term “storm” is used for the three events hereafter. One-minute (min) sea level data records from the Integrated Science Data Management (ISDM), which is a branch of Canada’s Federal Department of Fisheries and Oceans (DFO), were used in this study at six different locations (Fig. 1). Storm tracks and surface atmospheric pressure under the storms (Table 1) were provided by the Canadian Hurricane Centre (CHC), which is a division of the Meteorological Service of Canada. For calculating the mean storm speed across the shelf, only data corresponding to the storm tracks located between latitudes 44° and 49°30′ are considered.

The sea level records were detided using version 1.2 of T_TIDE (Pawlowicz et al. 2002). The long-wave records were extracted by low-pass filtering, with a cutoff frequency of 1.41 cpd. We used an orthogonal quadrature mirror filter type “Symmlet,” with eight vanishing moments (Mallat 2001), to decompose the signal into wavelet coefficients and reconstruct a low-pass-filtered signal by removing the coefficients corresponding to frequencies higher than 1.41 cpd. The result of filtering the detided sea level records is illustrated by Fig. 2 for station Argentia (AR) during the Florence event.

For each station except Lawrence (LW), the 1-h surface atmospheric pressure data used to calculate the expected hydrostatic response of the ocean and the wind speed (Figs. 3 –5) were both provided by the Canadian Daily Climate Data (CDCD), which is operated and maintained by Environment Canada. The times when the storm centers hit the coast were estimated from the storm reports provided by the CHC. The times when the storm’s edges hit the seaboard are graphically approximated from the rise and fall of the hydrostatic ocean response curves (Figs. 3 –5) in AR (the closest station to the storm center), taking into account the longer-term trends. The storm conditions were characterized as the sections above 30% of the maximum hydrostatic sea level elevation under the storm relative to the 2–3-day general trend. This gives an approximation of the storm durations (Table 1), with an arbitrary absolute error of 2 h. The storm width was calculated from the storm duration and mean speed (Table 1).

The three storms had roughly the same path. The intensity of Florence was greater than both Isaac and Chantal (Figs. 3 –5 and Table 1). The great amplitude of the sea level disturbance within the range of 12–45 cm (Table 2), observed at the coast after the passage of Florence (Fig. 3), suggested that the response was barotropic. Furthermore, Clarke and Brink (1985) showed that the response of such a wide, nonequatorial shelf to fluctuating wind forcing should be barotropic.

3. Analysis

a. Wavelet analysis

The classical windowed Fourier transform (Fourier transform performed on a sliding segment of length T from a time series) represents an inaccurate and inefficient method of time–frequency localization because it imposes a scale or a response interval T into the analysis (Torrence and Compo 1998, hereafter TC98). Because the time series of sea level data is nonstationary, with a large range of frequencies, especially during storm events, a wavelet analysis was employed. The latter is a useful tool for extracting both time and frequency information from a time series. Measuring the distance between each station, the barotropic CSW group speed is evaluated by approximating the time corresponding to the maximum wavelet energy at each station. Such evaluations have been used previously for studies of dispersive waves in solid systems (Önsay and Haddow 1994; Kishimoto et al. 1995; Park and Kim 2001) and Yanai waves in equatorial regions of the Atlantic Ocean (Meyers et al. 1993). A mother wavelet (the Morlet wavelet) is used, which is given by
i1520-0485-40-2-417-eq1
where w0 denotes the dimensionless frequency of the mother wavelet, here, taken to be six to satisfy the admissibility condition (Farge 1992; TC98), and η is a nondimensional “time” parameter.
The continuous wavelet transform (WT) of a discrete sequence xn (n = 0, … , N − 1; N is the number of points in the time series) with equal time spacing δt is defined as the convolution of xn, with a scaled, translated, and normalized version of ψ(η) (e.g., TC98):
i1520-0485-40-2-417-eq2
where the asterisk denotes the complex conjugate. The wavelet is stretched in time (t) by varying its scale (s), so that η = t/s and normalized to have unit energy. Although it is possible to calculate the wavelet transform using this convolution, it is considerably faster to do the calculations in Fourier space (TC98). The Fourier transform of the Morlet wavelet is
i1520-0485-40-2-417-eq3
where H(w) is a Heaviside step function, H(w) = 1 if w > 0, H(w) = 0 otherwise; s is the scale.
After finding the WT for all scales and assuming a first-order Markov process, a chi-squared distribution is used to calculate the significance levels at 95% for the corresponding red-noise Fourier power spectrum (Gilman et al. 1963):
i1520-0485-40-2-417-eq4
where M is the maximum lag, k is the frequency, and α is the lag-1 autocorrelation parameter, which is assumed to be equivalent to the autocorrelation function ρ(d) = γ(d)/γ(0), computed from the autocovariance
i1520-0485-40-2-417-eq5
where d is the lag (d = 1 for a first-order Markov process), x(n) is the sea level value of the nth point, and x is the mean of the time series xn. Values for α were estimated prior to low-pass filtering and fell within the range of 0.970–0.999 for the Florence event (Table 2), indicating a high degree of autocorrelation for the detided records.

All the calculations for the normalizations, the wavelet power spectrum (WPS), and the scale-averaged wavelet power (SAWP) followed the methodology described by TC98. To reduce ringing, longer time series (23 days) centered in the event were used for wavelet analysis. Figure 6 shows only a few days before and after the passage of Florence, well within the cone of influence.

The WPSs of Fig. 6a show an equatorward propagation of a 26–30-h period wave along the coast from LW to Yarmouth (YA), after Florence moved onto the shelf. The energy decay, resulting from dissipative effects, is shown by Fig. 6b. Note that the energy was greater in Port aux Basques (PB) than North Sydney (NS), and this might be because station PB is more sheltered (Fig. 1). The maximum energy reached AR a few hours after LW, showing a wave propagation in the opposite direction. The two stations are located at less than half of the storm width (Table 1) from the storm center and may be affected by the forced wave under the storm.

For the other two storms, local responses were observed in the WPSs in AR and LW (not shown) over shorter period bands, including the storm durations (20–24 and 17–21 h for storms Isaac and Chantal, respectively). Nonetheless, there was no significant propagation of sea surface disturbance along the coast.

For Florence, the maximum of each SAWP (Fig. 6b) gives a good approximation of the maximum amplitude of the barotropic CSW envelope (typically two or three oscillations). They are used to estimate the times when the CSW envelope reached each station. Thus, it is easy to calculate the group speed by estimating the lags for the pairs of stations and using the distance between them. However, it must be kept in mind that WT loses time resolution in the lower-frequency bands. The time and frequency spread of the WPSs are proportional to s and 1/s, respectively (Mallat 2001). Furthermore, the detided records contain red noise and residual tide, which increase the time error. In particular, prior to the low-pass filtering, the detided records showed a constant peak energy within the 11–13-h period band, indicating a residual of semidiurnal tide components not removed by detiding. For estimating the time uncertainty of the maximum amplitude of the CSW envelope, a Monte Carlo simulation is presumed to be fairly accurate. We created 1000 artificial 1-min-interval random red-noise signals frn of 215 min (∼546 h) with a zero mean for five different cases, with the α parameters previously calculated from the detided records related to the Florence event for stations LW, PB, NS, Halifax (HA), and YA (Table 2). To estimate the noise variance parameter, the noise was extracted (prior to filtering) by shrinkage of the wavelet coefficients, based on the Lorentz thresholding criterion (Katul and Vidakovic 1996). The advantage of this method is that the autocorrelation structure of the time series is taken into account. All the coefficients under the threshold are recomposed, and the variance of the noise was estimated (Table 2). To take into account the residual tide effects, a cosine wave fsd, from a period randomly taken within the range of 11–14 h, was added. Its amplitude was crudely estimated from the detided observation records by averaging the wave amplitudes of similar periods during Florence, using zero-crossing analysis. A Gaussian-modulated cosine wave (Y) of a period T0 randomly chosen between 26 and 30 h was added to roughly simulate the CSW from t = 240 h to t = 336 h (Fig. 7a). The signal had the form
i1520-0485-40-2-417-eq6
where A is the wave amplitude, t is the time variable, tm = 288 h is the time of the center of the peak induced by the Gaussian, c is a constant controlling the width of the Gaussian (here, c is arbitrarily chosen as 1440 by comparison with the detided time series), and θ is the wave phase (randomly taken within the range of 0 − π for each simulation). Each simulated time series was low-pass filtered with the same method used for the observation records. The WT, WPS, and SAWP (averaged between 26 and 30 h) of each residual were computed. The time uncertainty for each simulation was calculated as the difference between the time when the SAWP reaches its maximum amplitude and the time tm when the theoretical maximum of the cosine wave group occurs (Fig. 7b). The time uncertainty for each case/station was determinated as the 95th percentile of the simulations. The results are shown in Table 2 for different values of A, corresponding to the CSW amplitude observed at each station.

Table 3 shows the values used to calculate the barotropic CSW group speed along the coast after the passage of Florence. Lag uncertainties among station pairs were estimated by summing the time uncertainty at each station and were used for calculating group-velocity uncertainties (Table 3). Because of high uncertainties, results for station pairs HA–YA have not been used in the calculation. The mean CSW group speed was 11.4 ± 5.9 m s−1.

b. Cross-wavelet analysis

The cross-wavelet transform (XWT) of two signals exposes the regions in time–frequency space with high common power and further reveals information about the phase relationship (Grinsted et al. 2004). The XWTs of the pairs of stations are used to estimate the CSW phase speed. The XWT of two time series Xn and Yn is defined as WXY = WXWY*, where WX is the continuous WT of the time series Xn and the asterisk denotes the complex conjugate. The cross-wavelet power is defined as |WXY|. The complex argument of WXY describes the phase relationship between X and Y in a time–frequency space. As compared to traditional lagged cross-correlation methods, this technique offers the advantage that the time interval is reduced to the duration when the XWT of station pairs has significant power. Thus, the relative phase uncertainty is optimal for each pair and does not include additional errors for the studied event because of noise or other oscillations with similar frequencies, remaining out of the time boundary.

In this study, we calculated the XWT of pairs of the sea level time series in different locations to investigate the phase speed of the free CSW generated by Florence. As an example, results of these computations for station pair NS–YA are shown in Fig. 8. For calculating the mean relative phase between two time series, we used the circular mean of the phase over regions with higher than 5% statistical significance and within the 26–30-h period range. The circular mean of a set of angles (ai, i = 1, … , n) is given by Zar (1999)
i1520-0485-40-2-417-eq7

To estimate the confidence interval of the phase difference, we used the circular standard deviation defined as S = −2 ln(R/n), where R = X2 + Y2 (Grinsted et al. 2004).

The relative phase relationships can be considered as a measure of the time delay of two periodic signals expressed as a fraction of the wave period. Table 4 shows the values for different station pairs used to calculate the CSW phase velocities along the coast after the passage of Florence. The lags were calculated from the cross-wavelet phase angles (am) as lag = amTm, where Tm is the mean period. The lag errors were estimated as Δlag = lag(S/am + kσ̂T /Tm), where σ̂T is the standard deviation of the period and k is a constant that depends on the confidence level (here, k = 1.96 to satisfy the 95% confidence level).

Results of the XWT (Table 4) between LW and PB show no phase propagation of the CSW. Between LW and NS, the magnitude of the CSW phase speed is not realistic, and the large variability might result from direct influences of the storm. Indeed, it must be kept in mind that records at AR and LW were directly affected by the forced wave under the storm (Fig. 3). We conclude that all station pairs, which include LW, cannot be used in the calculations. The mean phase speed of the CSW calculated from Table 4 is about 16.0 m s−1 with a mean error of 5.1 m s−1.

4. Discussion

The wave observed after the passage of Florence across Newfoundland is expected to be a free CSW (Robinson 1969; Hamon 1966; LeBlond and Mysak 1978). Florence, the slowest storm, was the only one that excited a free CSW within a period band matching the passage of time of the event (28 ± 2 h) over the shelf (Figs. 3 –5). Sheng et al. (2006) observed from subtidal sea level records an equatorward propagation of shelf waves with similar periods (25–28 h) during the first few days after Hurricane Juan moved onto the Scotian shelf. The surface atmospheric pressure at the center of Juan when crossing the shelf was about 976 mb, slightly higher than Florence (964 mb). In terms of wind stress and barometric response of the sea surface, Florence, despite its relative slow speed, was much stronger than the two other storms, Isaac and Chantal. Consequently, one can expect a weaker response to the two faster but less intense storms. Tang et al. (1998)’s model showed that long barotropic shelf waves can easily be excited by a storm with barometric pressures similar to those under Isaac or Chantal. But surprisingly, no significant barotropic sea level disturbance propagation is observed for both Isaac and Chantal within the period ranges fitting the storm durations (Figs. 4 –5). Note that the storms might generate baroclinic disturbances whose surface manifestations are too weak to be observed with this dataset.

The dispersion relation of shelf waves can be determined from a given topographic shelf profile. The 2007 version of the computer program written by Brink and Chapman (1987) was used to compute the dispersion curves of the barotropic shelf wave (Fig. 10a) using the two cross-shelf sections A and B (Fig. 1a). Because of irregularities of the topography, simplified topographies with a cutoff water depth of 4000 m for section A and 3000 m for section B were used in the calculations (Fig. 9).

Low-frequency shelf waves have high group velocities and propagate faster than the maximum possible frequency, which has a zero group velocity. The significant regions on the WPSs (Fig. 6a) become narrower in period as one goes farther from LW, where the response includes periods from 16 to 64 h. The high-frequency waves (16–25 h) have a high wavenumber and a group velocity close or equal to zero (Fig. 10b). Their energy persists for a few hours on the southern–eastern Newfoundland shelf (AR and LW), with no propagation, and dissipates rapidly after the passage of Florence (Fig. 6a). The lower-frequency waves with high group velocities and smaller wavenumbers are weakly dispersive but lose energy resulting from dissipative effects while propagating from LW to YA. Although slowly dissipated, the six WPSs show a persistent period band (26–30 h), including Florence’s duration (28 h or 0.86 cpd), lower than the shelf wave maximum frequency that corresponds to the zero group-velocity regime. The group velocity (cg = dw/dk, where w is the frequency and k is the wavenumber) of a long shelf wave of frequency 0.86 ± 0.08 cpd, estimated from the dispersion curves (Fig. 10a), is within the range of 6.1–15.1 m s−1 and the phase velocity (w/k) is between 12.0 and 18.9 m s−1 (Fig. 10b) for shelf profiles at sections A and B. The mean speeds estimated in this study using wavelet and cross-wavelet analyses (11.4 ± 5.9 and 16 ± 5.1 m s−1 for group and phase speed, respectively) are consistent with these estimations.

Compared to previous observations, the speeds found were high for free-propagating shelf waves, but they are consistent with the computed dispersion curves and numerical models. Gordon and Huthnance (1987) and Tang et al. (1998) used numerical models and obtained dispersion relationships for the first mode of the shelf wave, with group velocities in the nondispersive regime of 13 m s−1 along the Scottish continental shelf and 11.6 m s−1 along the Labrador/Newfoundland shelf, respectively. Sheng et al. (2006) estimated from a nested grid model, a shelf wave phase speed of 15 m s−1 along the Scotian shelf. To our knowledge, speeds of this order of magnitude for a free-barotropic CSW with similar periods have not been detected in the past, using both the wavelet and cross-wavelet techniques.

Gordon and Huthnance (1987) and Tang et al. (1998) pointed out that the ocean response depends strongly on the duration of the atmospheric forcing. After Florence passed over LW and AR, the predominant frequency was slightly lower than the resonant nonpropagating regime for the five stations (Fig. 6a) and similar to the passage time for the storm over the shelf, which is consistent with the continental shelf wave theory. The duration of strong storms like Florence is the key factor that determines the propagation and the frequency band preferentially excited for the free CSW.

For the three storms, the nonpropagating regime is excited only at stations close to the storm track (not shown for storms Isaac and Chantal) in AR and LW (Fig. 6a). It lingers about one day after the storm passes. Section B is steeper than section A (Fig. 9). The former permits higher-frequency shelf waves with a maximum of 1.27 cpd (Fig. 10a), corresponding to the duration of Isaac (1.26 ± 0.08 cpd). Chantal’s duration (1.09 ± 0.08 cpd) was under this maximum but matched the nonpropagating regime in section A. Thus, the durations of Isaac and Chantal were too short to excite any propagating shelf waves.

The study points out the importance of topography and storm duration on the propagation of free CSWs generated by a moving low pressure system (Hamon 1966; Mooers and Smith 1968; Gordon and Huthnance 1987; Tang et al. 1998). The strong amplitude response of the coastal sea level (within the range of 12–45 cm) as a free-barotropic CSW after Florence impacted the Newfoundland shelf was clearly detected in time and frequency using wavelet and cross-wavelet analyses. The CSW period was similar to the storm duration, and the measured fast group and phase speeds were consistent with the dispersion curves for the first-mode CSW computed for the shelf profiles of Newfoundland and Nova Scotia.

Acknowledgments

We thank the DFO, the CHC, and Environment Canada for providing sea level data, storm details, atmospheric pressure, and wind data. The T_TIDE software package was provided by R. Pawlowicz et al. (available online at http://www.eos.ubc.ca/~rich/). The wavelet software was originally written by C. Torrence and G. Compo (available online at http://paos.colorado.edu/research/wavelets/) and the cross-wavelet software was provided by A. Grinsted et al. (available online at http://www.pol.ac.uk/home/research/waveletcoherence/). We are grateful for the advice of Dr. Kenneth Brink, who made the code for calculating the shelf wave dispersion relation (available online at http://www.whoi.edu/cms/files/Fortran_30425.htm). Suggestions provided by Pascal Sirguey were very helpful.

REFERENCES

  • Adams, J. K., and V. T. Buchwald, 1969: The generation of continental shelf waves. J. Fluid Mech., 35 , 815826.

  • Brink, K. H., and D. C. Chapman, 1987: Programs for computing properties of coastal-trapped waves and wind-driven motions over the continental shelf and slope. WHOI Tech. Rep. WHOI-87-24, 119 pp. [Available from Woods Hole Oceanographic Institution, Woods Hole, MA 02543].

    • Search Google Scholar
    • Export Citation
  • Brooks, D. A., and C. N. K. Mooers, 1977: Wind-forced continental shelf waves in the Florida Current. J. Geophys. Res., 82 , 25692576.

    • Search Google Scholar
    • Export Citation
  • Camayo, R., and E. Campos, 2006: Application of wavelet transform in the study of coastal trapped waves off the west coast of South America. Geophys. Res. Lett., 33 , L22601. doi:10.1029/2006GL026395.

    • Search Google Scholar
    • Export Citation
  • Clarke, A., and K. Brink, 1985: The response of stratified, frictional flow of shelf and slope waters to fluctuating large-scale, low-frequency wind forcing. J. Phys. Oceanogr., 15 , 439453.

    • Search Google Scholar
    • Export Citation
  • Cuevas, B., G. Brundrit, and A. Shipley, 1986: Low-frequency sea level fluctuations along the coasts of Namibia and South Africa. Geophys. J. Int., 87 , 3342.

    • Search Google Scholar
    • Export Citation
  • Cutchin, D. L., and R. L. Smith, 1973: Continental shelf waves: Low-frequency variations in sea level and currents over the Oregon continental shelf. J. Phys. Oceanogr., 3 , 7382.

    • Search Google Scholar
    • Export Citation
  • Enfield, D. B., and J. S. Allen, 1983: The generation and propagation of sea level variability along the Pacific coast of Mexico. J. Phys. Oceanogr., 13 , 10121033.

    • Search Google Scholar
    • Export Citation
  • Farge, M., 1992: Wavelet transforms and their applications to turbulence. Annu. Rev. Fluid Mech., 24 , 395457.

  • Gill, A. E., and E. H. Schumann, 1974: The generation of long shelf waves by the wind. J. Phys. Oceanogr., 4 , 8390.

  • Gilman, D. L., F. J. Fuglister, and J. M. Mitchell Jr., 1963: On the power spectrum of “red noise”. J. Atmos. Sci., 20 , 182184.

  • Gjevik, B., 1991: Simulations of shelf sea response due to travelling storms. Cont. Shelf Res., 11 , 139166.

  • Gordon, R. L., and J. M. Huthnance, 1987: Storm-driven continental shelf waves over the Scottish continental shelf. Cont. Shelf Res., 7 , 10151048.

    • Search Google Scholar
    • Export Citation
  • Grinsted, A., J. C. Moore, and S. Jevrejeva, 2004: Application of the cross-wavelet transform and wavelet coherence to geophysical time series. Nonlinear Processes Geophys., 11 , 561566.

    • Search Google Scholar
    • Export Citation
  • Hamon, B., 1962: The spectrums of mean sea level at Sydney, Coff’s Harbour, and Lord Howe Island. J. Geophys. Res., 67 , 51475155.

  • Hamon, B., 1966: Continental shelf waves and the effects of atmospheric pressure and wind stress on sea level. J. Geophys. Res., 71 , 28832893.

    • Search Google Scholar
    • Export Citation
  • Isozaki, I., 1968: An investigation on the variations of sea level due to meteorological disturbances on the coast of the Japanese Islands (II): Storm surges on the coast of the Japan Sea. J. Oceanogr. Soc. Japan, 24 , 178190.

    • Search Google Scholar
    • Export Citation
  • Katul, G., and B. Vidakovic, 1996: The partitioning of attached and detached eddy motion in the atmospheric surface layer using Lorentz wavelet filtering. Bound.-Layer Meteor., 77 , 153172.

    • Search Google Scholar
    • Export Citation
  • Kishimoto, K., H. Inoue, M. Hamada, and T. Shibuya, 1995: Time frequency analysis of dispersive waves by means of wavelet transform. J. Appl. Mech., 62 , 841.

    • Search Google Scholar
    • Export Citation
  • LeBlond, P., and L. Mysak, 1978: Waves in the Ocean. Elsevier Oceanography Series, Vol. 20, Elsevier, 602 pp.

  • Mallat, S., 2001: A Wavelet Tour of Signal Processing. 2nd ed. Academic Press, 637 pp.

  • Martinsen, E. A., B. Gjevik, and L. P. Roed, 1979: A numerical model for long barotropic waves and storm surges along the western coast of Norway. J. Phys. Oceanogr., 9 , 11261138.

    • Search Google Scholar
    • Export Citation
  • Meyers, S., B. Kelly, and J. O’Brien, 1993: An introduction to wavelet analysis in oceanography and meteorology: With application to the dispersion of yanai waves. Mon. Wea. Rev., 121 , 28582866.

    • Search Google Scholar
    • Export Citation
  • Mooers, C. N. K., and R. L. Smith, 1968: Continental shelf waves off Oregon. J. Geophys. Res., 73 , 549557.

  • Mysak, L. A., 1967a: On the theory of continental shelf waves. J. Mar. Res., 25 , 205227.

  • Mysak, L. A., 1967b: On the very low-frequency spectrum of the sea level on a continental shelf. J. Geophys. Res., 72 , 30433047.

  • Mysak, L. A., and B. V. Hamon, 1969: Low-frequency sea level behavior and continental shelf waves off North Carolina. J. Geophys. Res., 74 , 13971405.

    • Search Google Scholar
    • Export Citation
  • Önsay, T., and A. Haddow, 1994: Wavelet transform analysis of transient wave propagation in a dispersive medium. J. Acoust. Soc. Amer., 95 , 1441.

    • Search Google Scholar
    • Export Citation
  • Park, H., and D. Kim, 2001: Evaluation of the dispersive phase and group velocities using harmonic wavelet transform. NDT Int., 34 , 457467.

    • Search Google Scholar
    • Export Citation
  • Pawlowicz, R., B. Beardsley, and S. Lentz, 2002: Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE. Comput. Geosci., 28 , 929937.

    • Search Google Scholar
    • Export Citation
  • Robinson, A., 1964: Continental shelf waves and the response of sea level to weather systems. J. Geophys. Res., 69 , 367368.

  • Sheng, J., X. Zhai, and R. J. Greatbatch, 2006: Numerical study of the storm-induced circulation on the Scotian shelf during Hurricane Juan using a nested-grid ocean model. Prog. Oceanogr., 70 , 233254.

    • Search Google Scholar
    • Export Citation
  • Slørdal, L., E. Martinsen, and A. Blumberg, 1994: Modeling the response of an idealized coastal ocean to a traveling storm and to flow over bottom topography. J. Phys. Oceanogr., 24 , 16891705.

    • Search Google Scholar
    • Export Citation
  • Tang, C. L., Q. Gui, and B. M. DeTracey, 1998: Barotropic response of the Labrador/Newfoundland shelf to a moving storm. J. Phys. Oceanogr., 28 , 11521172.

    • Search Google Scholar
    • Export Citation
  • Torrence, C., and G. P. Compo, 1998: A practical guide to wavelet analysis. Bull. Amer. Meteor. Soc., 79 , 6178.

  • Yankovsky, A., and R. Garvine, 1998: Subinertial dynamics on the inner New Jersey shelf during the upwelling season. J. Phys. Oceanogr., 28 , 24442458.

    • Search Google Scholar
    • Export Citation
  • Zar, J., 1999: Biostatistical Analysis. 4th ed. Prentice-Hall, 931 pp.

Fig. 1.
Fig. 1.

Locations of sea level stations: AR (47°18′0″N, 53°59′0″W), LW (46°55′0″N, 55°23′24″W), PB (47°34′0″N, 59°08′0″W), NS (46°13′00″N, 60°15′0″W), HA (44°41′0″N, 63°37′0″W), and YA (43°50′0″N, 66°07′0″W) for the Atlantic coast of Canada (Newfoundland and Nova Scotia). Also shown are the tracks of storm centers (circles, crosses, and triangles for Florence, Isaac, and Chantal, respectively), the bathymetry contours of 100, 300, 2000, 3000, 4000, and 5000 m, and the cross-shelf sections A and B.

Citation: Journal of Physical Oceanography 40, 2; 10.1175/2009JPO4204.1

Fig. 2.
Fig. 2.

Detided sea level record relative to the mean sea level at AR for Florence before (thin line) and after (thick line) low-pass filtering (with a cutoff frequency of 1.41 cpd).

Citation: Journal of Physical Oceanography 40, 2; 10.1175/2009JPO4204.1

Fig. 3.
Fig. 3.

Detided sea level time series (low-pass filtered with a cutoff frequency of 1.41 cpd) relative to the mean sea level for each station (thick solid lines), barometric response of the sea surface (dotted lines) to the atmospheric pressure fluctuations, and wind speed (dashed lines) at stations AR, LW, PB, NS, HA, and YA. Note that neither the barometric response of the sea surface nor the wind speed is shown for LW because of the lack of data. Also shown are the times when Hurricane Florence (front, center, and back) hit the coast of Newfoundland (left, middle, and right thin vertical solid lines, respectively).

Citation: Journal of Physical Oceanography 40, 2; 10.1175/2009JPO4204.1

Fig. 4.
Fig. 4.

As in Fig. 3, but for storm Isaac.

Citation: Journal of Physical Oceanography 40, 2; 10.1175/2009JPO4204.1

Fig. 5.
Fig. 5.

As in Fig. 3, but for storm Chantal at stations AR, LW, NS, HA and YA. Note that the sea level record for AR is not represented because of lack of data.

Citation: Journal of Physical Oceanography 40, 2; 10.1175/2009JPO4204.1

Fig. 6.
Fig. 6.

(a) Morlet wavelet power spectrum of the detided low-pass-filtered sea level records for the six named stations, normalized by 1/(95% confidence level for the corresponding red-noise spectrum) at each scale (power relative to the significant level). Thick vertical lines indicate the time when Florence (front, center, and back of the storm) hit the coast of Newfoundland. (b) The scale-averaged wavelet power time series over the 26–30-h band for AR (Xs), LW (dashed line), PB (dash–dot line), NS (dotted line), HA (solid line), and YA (Os).

Citation: Journal of Physical Oceanography 40, 2; 10.1175/2009JPO4204.1

Fig. 7.
Fig. 7.

(a) One-min-interval artificial time series for LW Monte Carlo simulation, containing red noise (α = 0.970) and residual of semidiurnal tide components over the whole time series and a Gaussian-modulated cosine wave of amplitude 19 cm (period randomly chosen within the range of 26–30 h and random phase between 0 − π and from t = 240 h to t = 336 h). (b) Scale-averaged wavelet power time series of the low-pass residual of (a) over the 26–30-h period band. The solid and dashed vertical lines are the times of the maximum-averaged wavelet power and tm = 288 h, respectively.

Citation: Journal of Physical Oceanography 40, 2; 10.1175/2009JPO4204.1

Fig. 8.
Fig. 8.

XWT of detided sea level time series at NS and YA after the passage of Florence, normalized by 1/(95% confidence level for the corresponding red-noise spectrum) at each scale (power relative to the significant level). Also shown (as arrows) is the relative phase relationship (with in-phase pointing right, antiphase pointing left, and YA leading NS by 90° pointing straight down).

Citation: Journal of Physical Oceanography 40, 2; 10.1175/2009JPO4204.1

Fig. 9.
Fig. 9.

Simplified topography of sections A (solid lines) and B (dotted lines).

Citation: Journal of Physical Oceanography 40, 2; 10.1175/2009JPO4204.1

Fig. 10.
Fig. 10.

(a) Dispersion relations for the first mode of the CSW for sections A (solid lines) and B (dotted lines). The thin lines are the durations (represented as frequencies in cpd) of each storm with the estimated error. (b) The phase velocity c (thick solid and dashed line for sections A and B, respectively) and cg (thin solid and dashed line for sections A and B, respectively) are derived from the dispersion curves. The shaded region is the theoretical wavenumber range of the CSW generated by Florence, considering its duration’s error and dispersion relations for both sections A and B.

Citation: Journal of Physical Oceanography 40, 2; 10.1175/2009JPO4204.1

Table 1.

Characteristics of the three storms studied. Values are provided/calculated from the CHC and the CDCD data.

Table 1.
Table 2.

Values of lag-1 autocorrelation parameter (α), noise variance , amplitude of the semidiurnal tide residual (Asd), and amplitude of the CSW (A) related to Florence for the six stations. The 95% percentile for time uncertainty (Δt), resulting from the Monte Carlo simulation, are shown for all stations except AR.

Table 2.
Table 3.

Values used for the calculation of the CSW group speed (cg) after the passage of Florence. A positive value of the lag indicates that the second station of the pair leads the first. A positive value of cg indicates that the energy propagates alongshore with the coastline on its right.

Table 3.
Table 4.

Results of the XWT detided sea level records for different Station pairs. A positive value of the lag indicates that the second station of the pair leads the first. A positive value of c indicates that the wave phase propagates alongshore with the coastline on its right.

Table 4.
Save
  • Adams, J. K., and V. T. Buchwald, 1969: The generation of continental shelf waves. J. Fluid Mech., 35 , 815826.

  • Brink, K. H., and D. C. Chapman, 1987: Programs for computing properties of coastal-trapped waves and wind-driven motions over the continental shelf and slope. WHOI Tech. Rep. WHOI-87-24, 119 pp. [Available from Woods Hole Oceanographic Institution, Woods Hole, MA 02543].

    • Search Google Scholar
    • Export Citation
  • Brooks, D. A., and C. N. K. Mooers, 1977: Wind-forced continental shelf waves in the Florida Current. J. Geophys. Res., 82 , 25692576.

    • Search Google Scholar
    • Export Citation
  • Camayo, R., and E. Campos, 2006: Application of wavelet transform in the study of coastal trapped waves off the west coast of South America. Geophys. Res. Lett., 33 , L22601. doi:10.1029/2006GL026395.

    • Search Google Scholar
    • Export Citation
  • Clarke, A., and K. Brink, 1985: The response of stratified, frictional flow of shelf and slope waters to fluctuating large-scale, low-frequency wind forcing. J. Phys. Oceanogr., 15 , 439453.

    • Search Google Scholar
    • Export Citation
  • Cuevas, B., G. Brundrit, and A. Shipley, 1986: Low-frequency sea level fluctuations along the coasts of Namibia and South Africa. Geophys. J. Int., 87 , 3342.

    • Search Google Scholar
    • Export Citation
  • Cutchin, D. L., and R. L. Smith, 1973: Continental shelf waves: Low-frequency variations in sea level and currents over the Oregon continental shelf. J. Phys. Oceanogr., 3 , 7382.

    • Search Google Scholar
    • Export Citation
  • Enfield, D. B., and J. S. Allen, 1983: The generation and propagation of sea level variability along the Pacific coast of Mexico. J. Phys. Oceanogr., 13 , 10121033.

    • Search Google Scholar
    • Export Citation
  • Farge, M., 1992: Wavelet transforms and their applications to turbulence. Annu. Rev. Fluid Mech., 24 , 395457.

  • Gill, A. E., and E. H. Schumann, 1974: The generation of long shelf waves by the wind. J. Phys. Oceanogr., 4 , 8390.

  • Gilman, D. L., F. J. Fuglister, and J. M. Mitchell Jr., 1963: On the power spectrum of “red noise”. J. Atmos. Sci., 20 , 182184.

  • Gjevik, B., 1991: Simulations of shelf sea response due to travelling storms. Cont. Shelf Res., 11 , 139166.

  • Gordon, R. L., and J. M. Huthnance, 1987: Storm-driven continental shelf waves over the Scottish continental shelf. Cont. Shelf Res., 7 , 10151048.

    • Search Google Scholar
    • Export Citation
  • Grinsted, A., J. C. Moore, and S. Jevrejeva, 2004: Application of the cross-wavelet transform and wavelet coherence to geophysical time series. Nonlinear Processes Geophys., 11 , 561566.

    • Search Google Scholar
    • Export Citation
  • Hamon, B., 1962: The spectrums of mean sea level at Sydney, Coff’s Harbour, and Lord Howe Island. J. Geophys. Res., 67 , 51475155.

  • Hamon, B., 1966: Continental shelf waves and the effects of atmospheric pressure and wind stress on sea level. J. Geophys. Res., 71 , 28832893.

    • Search Google Scholar
    • Export Citation
  • Isozaki, I., 1968: An investigation on the variations of sea level due to meteorological disturbances on the coast of the Japanese Islands (II): Storm surges on the coast of the Japan Sea. J. Oceanogr. Soc. Japan, 24 , 178190.

    • Search Google Scholar
    • Export Citation
  • Katul, G., and B. Vidakovic, 1996: The partitioning of attached and detached eddy motion in the atmospheric surface layer using Lorentz wavelet filtering. Bound.-Layer Meteor., 77 , 153172.

    • Search Google Scholar
    • Export Citation
  • Kishimoto, K., H. Inoue, M. Hamada, and T. Shibuya, 1995: Time frequency analysis of dispersive waves by means of wavelet transform. J. Appl. Mech., 62 , 841.

    • Search Google Scholar
    • Export Citation
  • LeBlond, P., and L. Mysak, 1978: Waves in the Ocean. Elsevier Oceanography Series, Vol. 20, Elsevier, 602 pp.

  • Mallat, S., 2001: A Wavelet Tour of Signal Processing. 2nd ed. Academic Press, 637 pp.

  • Martinsen, E. A., B. Gjevik, and L. P. Roed, 1979: A numerical model for long barotropic waves and storm surges along the western coast of Norway. J. Phys. Oceanogr., 9 , 11261138.

    • Search Google Scholar
    • Export Citation
  • Meyers, S., B. Kelly, and J. O’Brien, 1993: An introduction to wavelet analysis in oceanography and meteorology: With application to the dispersion of yanai waves. Mon. Wea. Rev., 121 , 28582866.

    • Search Google Scholar
    • Export Citation
  • Mooers, C. N. K., and R. L. Smith, 1968: Continental shelf waves off Oregon. J. Geophys. Res., 73 , 549557.

  • Mysak, L. A., 1967a: On the theory of continental shelf waves. J. Mar. Res., 25 , 205227.

  • Mysak, L. A., 1967b: On the very low-frequency spectrum of the sea level on a continental shelf. J. Geophys. Res., 72 , 30433047.

  • Mysak, L. A., and B. V. Hamon, 1969: Low-frequency sea level behavior and continental shelf waves off North Carolina. J. Geophys. Res., 74 , 13971405.

    • Search Google Scholar
    • Export Citation
  • Önsay, T., and A. Haddow, 1994: Wavelet transform analysis of transient wave propagation in a dispersive medium. J. Acoust. Soc. Amer., 95 , 1441.

    • Search Google Scholar
    • Export Citation
  • Park, H., and D. Kim, 2001: Evaluation of the dispersive phase and group velocities using harmonic wavelet transform. NDT Int., 34 , 457467.

    • Search Google Scholar
    • Export Citation
  • Pawlowicz, R., B. Beardsley, and S. Lentz, 2002: Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE. Comput. Geosci., 28 , 929937.

    • Search Google Scholar
    • Export Citation
  • Robinson, A., 1964: Continental shelf waves and the response of sea level to weather systems. J. Geophys. Res., 69 , 367368.

  • Sheng, J., X. Zhai, and R. J. Greatbatch, 2006: Numerical study of the storm-induced circulation on the Scotian shelf during Hurricane Juan using a nested-grid ocean model. Prog. Oceanogr., 70 , 233254.

    • Search Google Scholar
    • Export Citation
  • Slørdal, L., E. Martinsen, and A. Blumberg, 1994: Modeling the response of an idealized coastal ocean to a traveling storm and to flow over bottom topography. J. Phys. Oceanogr., 24 , 16891705.

    • Search Google Scholar
    • Export Citation
  • Tang, C. L., Q. Gui, and B. M. DeTracey, 1998: Barotropic response of the Labrador/Newfoundland shelf to a moving storm. J. Phys. Oceanogr., 28 , 11521172.

    • Search Google Scholar
    • Export Citation
  • Torrence, C., and G. P. Compo, 1998: A practical guide to wavelet analysis. Bull. Amer. Meteor. Soc., 79 , 6178.

  • Yankovsky, A., and R. Garvine, 1998: Subinertial dynamics on the inner New Jersey shelf during the upwelling season. J. Phys. Oceanogr., 28 , 24442458.

    • Search Google Scholar
    • Export Citation
  • Zar, J., 1999: Biostatistical Analysis. 4th ed. Prentice-Hall, 931 pp.

  • Fig. 1.

    Locations of sea level stations: AR (47°18′0″N, 53°59′0″W), LW (46°55′0″N, 55°23′24″W), PB (47°34′0″N, 59°08′0″W), NS (46°13′00″N, 60°15′0″W), HA (44°41′0″N, 63°37′0″W), and YA (43°50′0″N, 66°07′0″W) for the Atlantic coast of Canada (Newfoundland and Nova Scotia). Also shown are the tracks of storm centers (circles, crosses, and triangles for Florence, Isaac, and Chantal, respectively), the bathymetry contours of 100, 300, 2000, 3000, 4000, and 5000 m, and the cross-shelf sections A and B.

  • Fig. 2.

    Detided sea level record relative to the mean sea level at AR for Florence before (thin line) and after (thick line) low-pass filtering (with a cutoff frequency of 1.41 cpd).

  • Fig. 3.

    Detided sea level time series (low-pass filtered with a cutoff frequency of 1.41 cpd) relative to the mean sea level for each station (thick solid lines), barometric response of the sea surface (dotted lines) to the atmospheric pressure fluctuations, and wind speed (dashed lines) at stations AR, LW, PB, NS, HA, and YA. Note that neither the barometric response of the sea surface nor the wind speed is shown for LW because of the lack of data. Also shown are the times when Hurricane Florence (front, center, and back) hit the coast of Newfoundland (left, middle, and right thin vertical solid lines, respectively).

  • Fig. 4.

    As in Fig. 3, but for storm Isaac.

  • Fig. 5.

    As in Fig. 3, but for storm Chantal at stations AR, LW, NS, HA and YA. Note that the sea level record for AR is not represented because of lack of data.

  • Fig. 6.

    (a) Morlet wavelet power spectrum of the detided low-pass-filtered sea level records for the six named stations, normalized by 1/(95% confidence level for the corresponding red-noise spectrum) at each scale (power relative to the significant level). Thick vertical lines indicate the time when Florence (front, center, and back of the storm) hit the coast of Newfoundland. (b) The scale-averaged wavelet power time series over the 26–30-h band for AR (Xs), LW (dashed line), PB (dash–dot line), NS (dotted line), HA (solid line), and YA (Os).

  • Fig. 7.

    (a) One-min-interval artificial time series for LW Monte Carlo simulation, containing red noise (α = 0.970) and residual of semidiurnal tide components over the whole time series and a Gaussian-modulated cosine wave of amplitude 19 cm (period randomly chosen within the range of 26–30 h and random phase between 0 − π and from t = 240 h to t = 336 h). (b) Scale-averaged wavelet power time series of the low-pass residual of (a) over the 26–30-h period band. The solid and dashed vertical lines are the times of the maximum-averaged wavelet power and tm = 288 h, respectively.

  • Fig. 8.

    XWT of detided sea level time series at NS and YA after the passage of Florence, normalized by 1/(95% confidence level for the corresponding red-noise spectrum) at each scale (power relative to the significant level). Also shown (as arrows) is the relative phase relationship (with in-phase pointing right, antiphase pointing left, and YA leading NS by 90° pointing straight down).

  • Fig. 9.

    Simplified topography of sections A (solid lines) and B (dotted lines).

  • Fig. 10.

    (a) Dispersion relations for the first mode of the CSW for sections A (solid lines) and B (dotted lines). The thin lines are the durations (represented as frequencies in cpd) of each storm with the estimated error. (b) The phase velocity c (thick solid and dashed line for sections A and B, respectively) and cg (thin solid and dashed line for sections A and B, respectively) are derived from the dispersion curves. The shaded region is the theoretical wavenumber range of the CSW generated by Florence, considering its duration’s error and dispersion relations for both sections A and B.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 628 331 104
PDF Downloads 249 42 5