Energy input and dissipation in a temperate lake during the spring transition

ADCP and temperature chain measurements have been used to estimate the rate of energy input by wind stress to the water surface in the south basin of Windermere. The energy input from the atmosphere was found to increase markedly as the lake stratified in spring. The efficiency of energy transfer (Eff), defined as the ratio of the rate of working in near-surface waters (RW) to that above the lake surface (P10), increased from ∼0.0013 in vertically homogenous conditions to ∼0.0064 in the first 40 days of the stratified regime. A maximum value of Eff∼0.01 was observed when, with increasing stratification, the first mode internal seiche period decreased to match the diurnal wind period of 24 h. The increase in energy input, following the onset of stratification was reflected in enhancement of the mean depth-varying kinetic energy without a corresponding increase in wind forcing. Parallel estimates of energy dissipation in the bottom boundary layer, based on determination of the structure function show that it accounts for ∼15% of RW in stratified conditions. The evolution of stratification in the lake conforms to a heating stirring model which indicates that mixing accounts for ∼21% of RW. Taken together, these estimates of key energetic parameters point the way to the development of full energy budgets for lakes and shallow seas

Effect of coastal boundary resolution and mixing upon internal wave generation and propagation in coastal regions

A non-linear two-dimensional vertically stratified cross-sectional model of a constant depth basin without rotation is used to investigate the influence of vertical and horizontal diffusion upon the wind-driven circulation in the basin and the associated temperature field. The influence of horizontal grid resolution, in particular the application of an irregular grid with high resolution in the coastal boundary layer is examined. The calculations show that the initial response to a wind impulse is downwelling at the downwind end of the basin with upwelling and convective mixing at the opposite end. Results from a two-layer analytical model show that the initial response is the excitation of an infinite number of internal seiche modes in order to represent the initial response which is confined to a narrow near coastal region. As time progresses, at the downwind end of the basin a density front propagates away from the boundary, with the intensity of its horizontal gradient and associated vertical velocity determined by both horizontal and vertical viscosity values. Calculations demonstrate the importance of high horizontal grid resolution in resolving this density gradient together with an accurate density advection scheme. The application of an irregular grid in the horizontal with high grid resolution in the nearshore region enables the initial response to be accurately reproduced although physically unrealistic short waves appear as the frontal region propagates onto the coarser grid. Parameterization of horizontal viscosity using a Smagorinsky-type formulation acts as a selective grid size-dependent filter, and removes the short-wave problem although enhanced smoothing can occur if the scaling coefficient in the formulation is too large. Calculations clearly show the advantages of using an irregular grid but also the importance of using a grid size-dependent filter to avoid numerical problems

A model study of tidal distributions in the Celtic and Irish Sea regions determined with finite volume and finite element models

An unstructured mesh model of the west coast of Britain, covering the same domain and using topography and open boundary forcing that are identical to a previous validated uniform grid finite difference model of the region, is used to compare the performance of a finite volume (FV) and a finite element (FE) model of the area in determining tide–surge interaction in the region. Initial calculations show that although qualitatively both models give comparable tidal solutions in the region, comparison with observations shows that the FV model tends to under-estimate tidal amplitudes and hence background tidal friction in the eastern Irish Sea. Storm surge elevations in the eastern Irish Sea due to westerly, northerly and southerly uniform wind stresses computed with the FV model tend to be slightly higher than those computed with the FE model, due to differences in background tidal friction. However, both models showed comparable non-linear tide–surge interaction effects for all wind directions, suggesting that they can reproduce the extensive tide–surge interaction processes that occur in the eastern Irish Sea. Following on from this model comparison study, the physical processes contributing to surge generation and tide–surge interaction in the region are examined. Calculations are performed with uniform wind stresses from a range of directions, and the balance of various terms in the hydrodynamic equations is examined. A detailed comparison of the spatial variability of time series of non-linear bottom friction and non-linear momentum advection terms at six adjacent nodes at two locations in water depths of 20 and 6 m showed some spatial variability from one node to another. This suggests that even in the near coastal region, where water depths are of the order of 6 m and the mesh is fine (of order 0.5 km), there is significant spatial variability in the non-linear terms. In addition, distributions of maximum bed stress due to tides and wind forcing in nearshore regions show appreciable spatial variability. This suggests that intensive measurement campaigns and very high-resolution mesh models are required to validate and reproduce the non-linear processes that occur in these regions and to predict extreme bed stresses that can give rise to sediment movement. High-resolution meshes will also be required in pollution transport problem