Implementing Super-Resolution of Non-Stationary Tides with Wavelets: An Introduction to CWT_Multi

Matthew Lobo aAtmospheric and Oceanic Sciences Program, Princeton University, Princeton, New Jersey
bDepartment of Civil and Environmental Engineering, Portland State University, Portland, Oregon

Search for other papers by Matthew Lobo in
Current site
Google Scholar
PubMed
Close
,
David A. Jay bDepartment of Civil and Environmental Engineering, Portland State University, Portland, Oregon

Search for other papers by David A. Jay in
Current site
Google Scholar
PubMed
Close
,
Silvia Innocenti cEnvironment and Climate Change Canada, Meteorological Research Division, Québec City, Québec, Canada

Search for other papers by Silvia Innocenti in
Current site
Google Scholar
PubMed
Close
,
Stefan A. Talke dCivil and Environmental Engineering, California Polytechnic State University, San Luis Obispo, California

Search for other papers by Stefan A. Talke in
Current site
Google Scholar
PubMed
Close
,
Steven L. Dykstra eSchool of the Earth, Ocean and Environment, University of South Carolina, Columbia, South Carolina
fCollege of Fisheries and Ocean Science, University of Alaska Fairbanks, Fairbanks, Alaska

Search for other papers by Steven L. Dykstra in
Current site
Google Scholar
PubMed
Close
, and
Pascal Matte cEnvironment and Climate Change Canada, Meteorological Research Division, Québec City, Québec, Canada

Search for other papers by Pascal Matte in
Current site
Google Scholar
PubMed
Close
Open access

Abstract

Tides are often non-stationary due to non-astronomical influences. Investigating variable tidal properties implies a tradeoff between separating adjacent frequencies (using long analysis windows) and resolving their time variations (short windows). Previous continuous wavelet transform (CWT) tidal methods resolved tidal species. Here, we present CWT_Multi, a Matlab code that: a) uses CWT linearity (via the “Response Coefficient Method”) to implement super-resolution (Munk and Hasselman 1964); b) provides a Munk-Hasselman constituent-selection criterion; and c) introduces an objective, time-variable form of inference (“dynamic inference”) based on time-varying data properties. CWT_Multi resolves tidal species on time-scales of days and multiple constituents per species with fortnightly filters. It outputs astronomical phase-lags and admittances, analyzes multiple records, and provides power spectra of the signal(s), residual(s) and reconstruction(s), confidence limits, and signal-to-noise ratios. Artificial data and water-levels from the Lower Columbia River Estuary (LCRE) and San Francisco Bay Delta (SFBD) are used to test CWT_Multi and compare it to harmonic analysis programs NS_Tide and UTide. CWT_Multi provides superior reconstruction, detiding, dynamic analysis utility, and time-resolution of constituents (but with broader confidence limits). Dynamic inference resolves closely spaced constituents (like K1, S1, and P1) on fortnightly time scales, quantifying impacts of diel power-peaking (with a 24-hour period, like S1) on water levels in the LCRE. CWT_Multi also helps quantify impacts of high flows and a salt-barrier closing on tidal properties in the SFBD. On the other hand, CWT_Multi does not excel at prediction, and results depend on analysis details, as for any method applied to non-stationary data.

© 2024 American Meteorological Society. This is an Author Accepted Manuscript distributed under the terms of the default AMS reuse license. For information regarding reuse and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: David A. Jay, djay@pdx.edu

Abstract

Tides are often non-stationary due to non-astronomical influences. Investigating variable tidal properties implies a tradeoff between separating adjacent frequencies (using long analysis windows) and resolving their time variations (short windows). Previous continuous wavelet transform (CWT) tidal methods resolved tidal species. Here, we present CWT_Multi, a Matlab code that: a) uses CWT linearity (via the “Response Coefficient Method”) to implement super-resolution (Munk and Hasselman 1964); b) provides a Munk-Hasselman constituent-selection criterion; and c) introduces an objective, time-variable form of inference (“dynamic inference”) based on time-varying data properties. CWT_Multi resolves tidal species on time-scales of days and multiple constituents per species with fortnightly filters. It outputs astronomical phase-lags and admittances, analyzes multiple records, and provides power spectra of the signal(s), residual(s) and reconstruction(s), confidence limits, and signal-to-noise ratios. Artificial data and water-levels from the Lower Columbia River Estuary (LCRE) and San Francisco Bay Delta (SFBD) are used to test CWT_Multi and compare it to harmonic analysis programs NS_Tide and UTide. CWT_Multi provides superior reconstruction, detiding, dynamic analysis utility, and time-resolution of constituents (but with broader confidence limits). Dynamic inference resolves closely spaced constituents (like K1, S1, and P1) on fortnightly time scales, quantifying impacts of diel power-peaking (with a 24-hour period, like S1) on water levels in the LCRE. CWT_Multi also helps quantify impacts of high flows and a salt-barrier closing on tidal properties in the SFBD. On the other hand, CWT_Multi does not excel at prediction, and results depend on analysis details, as for any method applied to non-stationary data.

© 2024 American Meteorological Society. This is an Author Accepted Manuscript distributed under the terms of the default AMS reuse license. For information regarding reuse and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: David A. Jay, djay@pdx.edu
Save