Adaptation of Classical Tidal Harmonic Analysis to Nonstationary Tides, with Application to River Tides
Abstract One of the most challenging areas in tidal analysis is the study of nonstationary signals with a tidal component, as they confront both current analysis methods and dynamical understanding. A new analysis tool has been developed, NS_TIDE, adapted to the study of nonstationary signals, in this case, river tides. It builds the nonstationary forcing directly into the tidal basis functions. It is implemented by modification of T_TIDE; however, certain concepts, particularly the meaning of a constituent and the Rayleigh criterion, are redefined to account for the smearing effects on the tidal spectral lines by nontidal energy. An error estimation procedure is included that constructs a covariance matrix of the regression coefficients, based on either an uncorrelated or a correlated noise model. The output of NS_TIDE consists of time series of subtidal water levels [mean water level (MWL)] and tidal properties (amplitudes and phases), expressed in terms of external forcing functions. The method was tested using records from a station on the Columbia River, 172 km from the ocean entrance, where the tides are strongly altered by river flow. NS_TIDE hindcast explains 96.4% of the signal variance with a root-mean-square error of 0.165 m obtained from 288 parameters, far better than traditional harmonic analysis (38.5%, 0.604 m, and 127 parameters). While keeping the benefits of harmonic analysis, its advantages compared to existing tidal analysis methods include its capacity to distinguish frequencies within tidal bands without losing resolution in the time domain or data at the endpoints of the time series.