Infaunal burrow ventilation and pore-water transport in muddy sediments

The ventilation of burrows by tube-dwelling benthos is understood to be important in determining rates of exchange of solutes between the sediment and overlying water. However, few models have attempted to link the burrow ventilation behavior of tube-dwelling organisms with their geochemical consequences. The classic cylinder model of bioirrigation in muddy sediments (Aller, R.C., 1980. Quantifying solute distributions in the bioturbated zone of marine sediments by defining an average microenvironment. Geochimica et Cosmochimica Acta 44, 1955–1965) links pore-water processes and burrow sizes and distributions in the sediment by assuming that burrows are fully flushed. The equivalence between the cylinder model and the more commonly used one-dimensional non-local exchange model depends upon this assumption. However, this assumption has seldom been tested in the field. We have extended the cylinder model of bioirrigation to include burrow ventilation activities of organisms. Burrow ventilation is modeled as a simple non-local exchange of burrow water with overlying water. Model simulations indicate that burrow ventilation has a large effect on vertical profiles and fluxes of solute tracers. We collected data on burrow geometry in the field by CT-scanning freshly collected sediment cores. At the same study sites, we measured activity profiles of ²²²Rn, a naturally occurring radionuclide tracer of pore-water transport. With model geometry independently constrained, we tested the model by comparing our model-predicted profiles with measured profiles. Our results demonstrate that burrows in the field are not fully flushed. Our estimated burrow ventilation rates compare favorably with previous laboratory measurements. The inclusion of realistic burrow ventilation in this pore-water transport model strongly affects modeled solute profiles and fluxes. We demonstrate how model parameters can be determined from field samples and present a model that more realistically simulates pore-water transport processes in muddy sediments.     

Nonlinear transfer function modelling of beach morphology at Duck,North Carolina

This paper presents a simple nonlinear data-based modelling approach for predicting the beach profile volume at Duck, North Carolina, USA. The state-dependent parameter form of the general transfer function (SDP TF) model is used to describe nonlinearity influencing these morphological data in two case examples. Case 1 investigates the nonlinearity associated with the dependency of wave forcing on the preceding beach volume. Case 2 investigates the ability to model the variables within the well-known diffusion equation for beach volume using this data-based approach. The results of this study show that the SDP TF approach can be used successfully to develop statistically robust models for describing nonlinearity in beach morphological systems. Furthermore, these models are shown to predict the beach volumes over both short (1 month ahead) and long (2 years ahead) time periods, and thus show great potential for practical applications in coastal zone management and engineering.     

Velocity potential formulations of highly accurate Boussinesq-type models

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.     

Tsunami generation,propagation, and run-up with a high-order Boussinesq model

In this work we extend a high-order Boussinesq-type (finite difference) model, capable of simulating waves out to wavenumber times depth kh < 25, to include a moving sea-bed, for the simulation of earthquake- and landslide-induced tsunamis. The extension is straight forward, requiring only an additional term within the kinematic bottom condition. As first test cases we simulate linear and nonlinear surface waves generated from both positive and negative impulsive bottom movements. The computed results compare well against earlier theoretical, numerical, and experimental values. Additionally, we show that the long-time (fully nonlinear) evolution of waves resulting from an upthrusted bottom can eventually result in true solitary waves, consistent with theoretical predictions. It is stressed, however, that the nonlinearity used far exceeds that typical of geophysical tsunamis in the open ocean. The Boussinesq-type model is then used to simulate numerous tsunami-type events generated from submerged landslides, in both one and two horizontal dimensions. The results again compare well against previous experiments and/or numerical simulations. The new extension compliments recently developed run-up capabilities within this approach, and as demonstrated, the model can therefore treat tsunami events from their initial generation, through their later propagation, and final run-up phases. The developed model is shown to maintain reasonable computational efficiency, and is therefore attractive for the simulation of such events, especially in cases where dispersion is important.     

Effect of wave amplitude on turbulent flow in a wavy channel by direct numerical simulation

The present study numerically investigates the characteristics of three-dimensional turbulent flow in a wavy channel. For the purpose of a careful observation of the effect of the wave amplitude on the turbulent flow, numerical simulations are performed at a various range of the wave amplitude to wavelength ratio (0.01?α/λ?0.05), where the wavelength is fixed with the same value of the mean channel height (H). The immersed boundary method is used to handle the wavy surface in a rectangular grid system, using the finite volume method. The Reynolds number (Re=U_bH/ν) based on the bulk velocity (U_b) is fixed at 6760. The present computational results for a wavy surface are well compared with those of references. When α/λ=0.02, the small recirculating flow occurs near the trough at the instant, but the mean reverse flow is not observed. In the mean flow field, the reverse flow appears from α/λ=0.03 among the wave amplitude considered in this study. The domain of the mean reverse flow defined by the locations of separation and reattachment depends strongly on the wave amplitude. The pressure drag coefficient augments with increasing the wave amplitude. The friction drag coefficient shows the increase and decrease behavior according to the wave amplitude. The quantitative information about the flow variables such as the distribution of pressure and shear stress on the wavy surface is highlighted.     

The tidal range energy potential of the West Coast of the United Kingdom

With concerns mounting over the UK’s energy future and the effects of climate change, it will soon become paramount that all viable sources of renewable energy are fully exploited. This study has examined the scope for reliable and fully predictable tidal electricity generation from the conjunctive operation of 5 major estuary barrages on the West Coast of the UK in an attempt to establish the potential scale of the extractable resources. Two levels of investigation have been undertaken: simple 0-D (‘two-tank’) modelling of barrage energy generation under different operational modes, using the hydraulic characteristics of turbine performance; and 2-D modelling of tidal hydrodynamics over a wide sea area in a computational grid incorporating the barrages with turbines and sluices. It has been demonstrated that more than 33TWh per year of electricity should be attainable, from 22GW of installed capacity, this representing close to 10% of present UK demand.     

Response of a climate model to tidal mixing parameterization under present day and last glacial maximum conditions

Experiments with a climate model were conducted under present day and last glacial maximum conditions in order to examine the model’s response to a vertical mixing scheme based on internal tide energy dissipation. The increase in internal tide energy flux caused by a 120 m reduction in sea level had the expected effect on diffusivity values, which were higher under lower sea level conditions. The impact of this vertical diffusivity change on the Atlantic meridional overturning is not straightforward and no clear relationship between diffusivity and overturning is found. There exists a weak positive correlation between overturning and changes to the power consumed by vertical mixing. Most of the climatic response generated by sea level change was not related to alterations in the internal tide energy flux but rather to the direct change in sea level itself.     

Analysis of seabed instability using element free Galerkin method

Wave-induced seabed instability, either momentary liquefaction or shear failure, is an important topic in ocean and coastal engineering. Many factors, such as seabed properties and wave parameters, affect the seabed instability. A non-dimensional parameter is proposed in this paper to evaluate the occurrence of momentary liquefaction. This parameter includes the properties of the soil and the wave. The determination of the wave-induced liquefaction depth is also suggested based on this non-dimensional parameter. As an example, a two-dimensional seabed with finite thickness is numerically treated with the EFGM meshless method developed early for wave-induced seabed responses. Parametric study is carried out to investigate the effect of wavelength, compressibility of pore fluid, permeability and stiffness of porous media, and variable stiffness with depth on the seabed response with three criteria for liquefaction. It is found that this non-dimensional parameter is a good index for identifying the momentary liquefaction qualitatively, and the criterion of liquefaction with seepage force can be used to predict the deepest liquefaction depth.