Revising the paradigm of tidal analysis – the uses of non-stationary data

Surface elevation and current records contain non-tidal variance, often dismissed as noise. The processes responsible for the non-tidal component may also modulate the tidal signal, altering its strength and frequency structure. Because of their manner of generation and propagation, internal tides are inherently irregular. The non-stationary character of these and other tidal processes provides an integral and useful property of tidal records, because it provides an opportunity to obtain insights into tidal dynamics and the interaction of tidal and non-tidal processes. It is, moreover, productive to use multiple approaches in analyzing coastal and estuarine tidal processes so that both the time-varying and average frequency content are determined. Only by confronting the causes of non-stationary behaviour in this way can some of the remaining challenges in tidal analysis and prediction be overcome, e.g. shelf and estuarine currents, river tides, internal tides, tide-surge interactions and tidally influenced ecological processes. Several examples illustrate the utility of non-stationary tidal analysis methods.Responsible Editor: Jens Kappenberg     

A methodology to estimate renewal time scales in estuaries: the Tagus Estuary case

High-resolution hydrodynamic models are a common tool to simulate water dynamics in estuaries. Results from these models are, however, difficult to interpret without the aid of additional parameters to integrate the information. In this paper a methodology to understand the transport patterns in the Tagus Estuary is proposed. It is based on the computation of two renewal time scales: residence time and integrated water fraction. This last parameter is used to build a dependency matrix that gives the integrated influence of each region of the estuary at a selected point. The parameters are computed using a Lagrangian transport model coupled to the hydrodynamic model. Results show that Tagus Estuary has two different types of regions: the central part of the estuary, with low renewal efficiency, and three regions with higher renewal efficiency. Renewal mechanisms are, however, different for each region as shown by the dependency matrix. Comparison of renewal time scales with results from a water-quality model revealed that residence time is not a limiting parameter for primary production in the Tagus Estuary.Responsible Editor: Hans Burchard     

A theoretical investigation of bed-form shapes

The deformation of movable boundaries under the action of an applied turbulent shear stress is well known. The resulting bed forms often are highly organized and nearly two-dimensional, which makes them an intriguing focus of study considering that they are generated in both steady and oscillatory turbulent flows. Many past studies share a common approach in which an infinitesimal perturbation is prescribed and the resulting growth or decay patterns are examined. In this approach, the bed forms are usually sinusoidal and the perturbation analysis does not provide a theoretical prediction of equilibrium bed-form geometry. An alternative approach is suggested here in which the forcing terms (pressure and stress) are prescribed parametrically and the governing equations are solved for the flow velocity and the associated boundary deformation. Using a multilayered approach, in which the bottom boundary layer is divided into a discrete, yet, arbitrary number of finite layers, analytical solutions for the horizontal current and bed profile are derived. The derivations identify two nondimensional parameters, p₀/u₀² and ₀/kh₀u₀², which modulate the amplitude of the velocity fluctuations and boundary deformation. For the case of combined pressure and stress divergence anomalies, the magnitude of the front face and lee slopes exhibit an asymmetry that is consistent with observed bed forms in steady two-dimensional flows.Responsible Editor: Jens Kappenberg     

Migrating sand waves

A simple mathematical model is described, which reproduces the major features of sand waves' appearance and growth and in particular predicts their migration speed. The model is based on the linear stability analysis of the flat configuration of the sea bottom subject to tidal currents. Attention is focused on the prediction of the complex growth rate that bottom perturbations undergo because of both oscillatory fluid motions and residual currents. While the real part _r of controls the amplification or decay of the amplitude of the bedforms, the imaginary part _i is related to their migration speed. Previous works on the migration of the sand waves (Németh etal. 2002) consider a forcing tide made up by the M2 constituent (oscillatory period equal to 12 h) plus the residual current Z0 and predict always a downcurrent migration of the bedforms. However, field cases exist of upcurrent-migrating sand waves (downcurrent/upcurrent-migrating sand waves mean bedforms moving in the direction of the steady residual tidal current or in the opposite direction, respectively). The inclusion of a tide constituent characterized by a period of 6 h (M4) is the main novelty of the present work, which allows for the prediction of the migration of sand waves against the residual current Z0. Indeed, the M4 tide constituent, as does also the residual current Z0, breaks the symmetry of the problem forced only by the M2 tide constituent, and induces sand-wave migration. The model proposed by Besio etal. (2003a) forms the basis for the present analysis. Previous works on the subject (Gerkema 2000; Hulscher 1996a,b; Komarova and Hulscher 2000) are thus improved by using a new solution procedure (Besio etal. 2003a) which allows for a more accurate evaluation of the growth rate for arbitrary values of the parameter r, which is the ratio between the horizontal tidal excursion and the perturbation wavelength. Responsible Editor: Jens Kappenberg     

Implementation of configuration dependent stiffness proportional damping for the dynamics of rigid multi-block systems

The distinct element method (DEM) has been used successfully for the dynamic analysis of rigid block systems. One of many difficulties associated with DEM is modeling of damping. In this paper, new procedures are proposed for the damping modeling and its numerical implementation in distinct element analysis of rigid multi-block systems. The stiffness proportional damping is constructed for the prescribed damping ratio, based on the non-zero fundamental frequency effective during the time interval while the boundary conditions remain essentially constant. At this time interval, the fundamental frequency can be estimated without complete eigenvalue analysis. The damping coefficients will vary while the damping ratio remains the same throughout the entire analysis. A new numerical procedure is developed to prevent unnecessary energy loss that can occur during the separation phases. These procedures were implemented in the development of the distinct element method for the dynamic analyses of piled multi-block systems. The analysis results for the single-block and two-block systems were in a good agreement with the analytic predictions. Applications to the seismic analyses of piled fourblock systems revealed that the new procedures can make a significant difference and may lead to much-improved results.     

The Ensemble Kalman Filter: theoretical formulation and practical implementation

The purpose of this paper is to provide a comprehensive presentation and interpretation of the Ensemble Kalman Filter (EnKF) and its numerical implementation. The EnKF has a large user group, and numerous publications have discussed applications and theoretical aspects of it. This paper reviews the important results from these studies and also presents new ideas and alternative interpretations which further explain the success of the EnKF. In addition to providing the theoretical framework needed for using the EnKF, there is also a focus on the algorithmic formulation and optimal numerical implementation. A program listing is given for some of the key subroutines. The paper also touches upon specific issues such as the use of nonlinear measurements, in situ profiles of temperature and salinity, and data which are available with high frequency in time. An ensemble based optimal interpolation (EnOI) scheme is presented as a cost-effective approach which may serve as an alternative to the EnKF in some applications. A fairly extensive discussion is devoted to the use of time correlated model errors and the estimation of model bias.Responsible Editor: Jörg-Olaf Wolff     

A consistent derivation of the wave-energy equation from basic hydrodynamic principles

Based on a decomposition of the velocity into mean flow, turbulent and wave components, momentum and hereafter a wave-energy equation is derived. It contains a turbulent energy dissipation term which is closed by applying a wave-related mixing length model and linear wave theory solutions. This closure produces a non-linear turbulent wave-energy dissipation including the wave energy in a 5/2 power law. The theory is able to predict correctly the shape of deep-water wave spectra according to Phillips' similarity law.Responsible Editor: Hans Burchard