Difficulties in separating hurricane induced effects from natural benthic succession: Hurricane Isabel,a case study from Eastern Virginia,USA

Hurricane Isabel reached the Eastern seaboard of North America on 18 September 2003 causing estimated damage >3 billion US dollars and the death of ∼50 people. Isabel is considered to be one of the most significant tropical cyclones to affect Virginia, since the Chesapeake Potomac Hurricane of 1933 and Hurricane Hazel in 1954. A study of the temporal changes in the benthic fauna pre- and post-hurricane was conducted on an intertidal sandflat within the dynamic barrier island system near Wachapreague, Eastern Virginia. Replicate sediment cores were collected 3 weeks before Isabel made landfall and further samples were collected on 5 occasions over the following 20 months. An immediate effect of Isabel was a doubling in the number of species, a significant increase in invertebrate species diversity (H′) and a rise in opportunistic species and deposit feeders, but a non-significant increase in the total number of organisms. Changes in infauna occurred such that by the end of the study there were significantly increased numbers of species, faunal abundances and community diversity measures, as compared with pre-hurricane samples, suggesting a potentially positive medium-term effect of this hurricane perturbation. The most notable direct effects of the hurricane were on the relative abundances of feeding guilds with a reduction in interface feeders from 87% pre-hurricane to 64% post-hurricane, and an increase in surface deposit feeders from 7% pre-hurricane to 20% post-hurricane. The study highlights potential problems in interpreting post-perturbation data when insufficient pre-perturbation data exist.     

Fundamental study for morphodynamic modelling: Sand mounds in oscillatory flows

Experiments on sand mounds in oscillatory flow, undertaken in controlled, large-scale laboratory conditions, have produced well-defined data sets for model comparison. Three bathymetries with different levels of submergence, including a surface-piercing case, were tested. The maximum slope was about 1:5.5. Sediment transport is due to bed load with ripple formation. The principal time-dependent bulk parameters are the vertical distance of the centre of gravity above the base and the volume of the mound. A semi-implicit finite-volume depth-averaged hydrodynamic model is used to drive morphodynamics, using van Rijn's sediment flux model generalized to take account of bed slope, and some justification is given for depth-averaged modeling in these conditions. Starting the model runs with the conditions at the end of the first cycle avoided initial atypical physical behaviour. In general good predictions were obtained with an angle of repose reduced from the standard value of about 30° for stationary beds to 15°. For these situations, morphodynamics was largely unaffected by a hydrodynamic roughness height in the range 2.5D₅₀ to 51D₅₀, with larger values accounting for ripple roughness. The reduced angle of repose may be physically expected with mobile beds but this specific value is only expected to be suited to this form of bed motion. In one case an exaggerated ripple formed near the top of the mound reducing agreement with experiment. For the submerged case with normal ripple structure excellent predictions were obtained. For the initially surface-piercing mound, the time of submergence was better predicted with a 30° angle of repose, presumably due to the prominent influence of the near stationary bed near the wet/dry interface, although long term predictions were better predicted with 15°. The occurrence of vortex shedding in the first cycle modeled was in agreement with experimental observation.     

Nitrification in Kochi backwaters

Nitrification rates, as oxidation of ¹⁵N-labelled ammonium and loss of nitrite from N-Serve treated samples, were measured in Kochi backwaters during three seasons. Nitrification rates ranged from undetectable to 166 nmol N L⁻¹ h⁻¹ in the water column and up to 17 nmol N (g wet wt)⁻¹ h⁻¹ in sediments. Nitrification rates were higher in intermediate salinities than in either freshwater or seawater end. Within this salinity range, nitrification rates could be related to ammonium concentrations. As shown by the relation between ammonification and nitrification rates, it is also likely that nitrification is more regulated by renewal rates, rather than by in situ concentrations, of substrate. Among other environmental parameters, temperature and pH may have an influence on nitrification. Potential nitrification rates calculated from loss of nitrite from N-Serve treated, nitrite-enriched samples were about 800 nmol N L⁻¹ h⁻¹ in the water column and 40 nmol N (g wet wt)⁻¹ h⁻¹ in sediments. While these rates are in balance with those of biological ammonium production they may be inadequate to mitigate ammonium pollution in this estuary.     

Developing a building policy for the erosion zone: Solutions to some key (Dutch) questions

In those coastal communities where the most seaward strip of mainland consists of dunes, these dunes often serve as a flexible sea defence. In addition, this strip offers large potential for housing and commercial enterprises. Unfortunately, due to severe storm surges part of this strip (the erosion zone) is subject to erosion, and as a result of which any buildings or infrastructure located here, are destroyed. Therefore, as we will illustrate in this paper, a building policy for this zone should reflect a compromise between two opposite interests: exploitation of the existing potential and, prevention of an unacceptable high risk due to erosion. Accordingly, the authors have developed a framework for such a building policy on the basis of which the desirability of various different types of investments and the location within the erosion zone of such investments can be determined. The examples that are used to illustrate this framework in this paper are limited to experiences in The Netherlands as relevant data and experiences are available and relatively easy accessible here. Nevertheless, the approach as is described is generic and applicable worldwide suggesting that the discovered unused potential for exploitation is not just limited to The Netherlands.     

The effect of a topographic high on the morphological stability of a two-inlet bay system

The cross-sectional stability of two tidal inlets connecting the same back-barrier lagoon to the ocean is investigated. The condition for equilibrium is the cross-sectional area tidal prism relationship. In an earlier study [Van de Kreeke, J., 1990. Can multiple inlets be stable? Estuarine, Coastal and Shelf Science 30: 261–273.], using the same equilibrium condition, it was concluded that where two inlets connect the same basin to the ocean ultimately one inlet will close. One of the major assumptions in that study was that the water level in the basin fluctuated uniformly. In hindsight this assumption might be too restrictive. For example, in the Wadden Sea the back barrier lagoon consists of a series of basins, rather than one single basin, separated by topographic highs. These topographic highs limit but do not exclude the exchange of water between the sub-basins. For this reason in the present study, a topographic high in the form of a weir was added, separating the back-barrier lagoon in two sub-basins. The water level in the sub-basins, rather than in the back-barrier as a whole, is assumed to fluctuate uniformly. For this schematization the hydrodynamic equations are solved using a finite difference method. The results, together with the equilibrium condition, yield the equilibrium flow curve for each of the inlets. The intersections of the two equilibrium flow curves represent combinations of cross-sectional areas for which both inlets are in equilibrium. The stability of the equilibriums was investigated using a non-linear stability analysis resulting in a flow diagram. Calculations were carried out for different inlet and weir characteristics. Sinussoidal tides were the same for both inlets. The results show that for relatively large wetted cross-sectional areas over the topographic high, approaching the situation of a single basin, there are no combinations of inlet cross-sectional areas for which both inlets are in a stable equilibrium. This supports the conclusion in the earlier study. For relatively small wetted cross-sectional areas over the topographic high there is one set of stable equilibriums. In that case the two-inlet bay system approaches that of two single-inlet bay systems.     

Large-scale dune erosion tests to study the influence of wave periods

Large-scale physical model tests were performed to quantify the effects of the wave period on dune erosion. Attention was focussed on 2D cross-shore effects in a situation with sandy dunes and extreme water levels and wave conditions. Besides profile measurements, detailed measurements in time and space of water pressure, flow velocities and sediment concentrations were performed in the near near-shore area. It was concluded that a longer wave period leads to a larger dune erosion volume and to a larger landward retreat of the dune face. Tests with double-peaked wave spectra showed that the influence of the spectral shape on dune erosion was best represented by the T_m − 1,0 spectral mean wave period, better than the peak wave period, T_p. The effect of the wave period on dune erosion was implemented in a dune erosion prediction method that estimates erosion volumes during normative storm conditions for the Dutch coast. More details of the measurements and additional analyses of physical processes are described in an accompanying paper by Van Thiel de Vries et al. [Van Thiel de Vries, J.S.M., van Gent, M.R.A., Reniers, A.J.H.M. and Walstra, D.J.R., submitted for publication. Analysis of dune erosion processes in large scale flume experiments, In this volume of Coastal Engineering.].     

Assessment of boundary-element method for modelling a free-floating sloped wave energy device. Part 1: Numerical modelling

The boundary-element method has been widely used as a design tool in the offshore and ship building industry for more than 30 years. Its application to wave energy conversion is, however, more recent. This paper deals with the numerical modelling of a free-floating sloped wave energy device. The power take-off mechanism of the device consists of an immersed tube with a piston sliding inside. The modelling is done using the boundary-element method package WAMIT. The model is first worked out for the case where the axis of the tube is vertical. It is then derived for the tube inclined and successfully verified against numerical benchmark data. A companion paper presents results of a detailed comparison with a physical model study.     

Solution of the propagation of the waves in open channels by the transfer matrix method

Many problems in mechanics can be solved by the use of the transfer matrix method. The use of this method in hydraulics engineering is not widespread and only limited studies are available. In this study, linearized St. Venant equations were used and the use of transfer matrix in ocean engineering was investigated for long waves in open channels, and numerical application was carried out. The results obtained through the transfer matrix method, which is quite easy to use, program and comprehend, showed similar results obtained from the characteristics method and finite differences method.     

Trophic modeling of the Northern Humboldt Current Ecosystem, Part II: Elucidating ecosystem dynamics from 1995 to 2004 with a focus on the impact of ENSO

The Northern Humboldt Current Ecosystem is one of the most productive in the world in terms of fish production. Its location near to the equator permits strong upwelling under relatively low winds, thus creating optimal conditions for the development of plankton communities. These communities ultimately support abundant populations of grazing fish such as the Peruvian anchoveta, Engraulis ringens. The ecosystem is also subject to strong inter-annual environmental variability associated with the El Niño Southern Oscillation (ENSO), which has major effects on nutrient structure, primary production, and higher trophic levels. Here our objective is to model the contributions of several external drivers (i.e. reconstructed phytoplankton changes, fish immigration, and fishing rate) and internal control mechanisms (i.e. predator-prey) to ecosystem dynamics over an ENSO cycle. Steady-state models and time-series data from the Instituto del Mar del Perú (IMARPE) from 1995 to 2004 provide the base data for simulations conducted with the program Ecopath with Ecosim. In simulations all three external drivers contribute to ecosystem dynamics. Changes in phytoplankton quantity and composition (i.e. contribution of diatoms and dino- and silicoflagellates), as affected by upwelling intensity, were important in dynamics of the El Niño of 1997–98 and the subsequent 3 years. The expansion and immigration of mesopelagic fish populations during El Niño was important for dynamics in following years. Fishing rate changes were the most important of the three external drivers tested, helping to explain observed dynamics throughout the modeled period, and particularly during the post-El Niño period. Internal control settings show a mix of predator–prey control settings; however a “wasp-waist” control of the ecosystem by small pelagic fish is not supported.     

Seasonal changes in submarine groundwater discharge to coastal salt ponds estimated using Ra and Ra as tracers

Submarine groundwater discharge (SGD) to coastal southern Rhode Island was estimated from measurements of the naturally-occurring radioisotopes ²²⁶Ra (t_1/2 = 1600 y) and ²²⁸Ra (t_1/2 = 5.75 y). Surface water and porewater samples were collected quarterly in Winnapaug, Quonochontaug, Ninigret, Green Hill, and Pt. Judith–Potter Ponds, as well as nearly monthly in the surface water of Rhode Island Sound, from January 2002 to August 2003; additional porewater samples were collected in August 2005. Surface water activities ranged from 12–83 dpm 100 L^− 1 (60 dpm = 1 Bq) and 21–256 dpm 100 L^− 1 for ²²⁶Ra and ²²⁸Ra, respectively. Porewater ²²⁶Ra activities ranged from 16–736 dpm 100 L^− 1 (2002–2003) and 95–815 dpm 100 L^− 1 (2005), while porewater ²²⁸Ra activities ranged from 23–1265 dpm 100 L^− 1. Combining these data with a simple box model provided average ²²⁶Ra-based submarine groundwater fluxes ranging from 11–159 L m^− 2 d^− 1 and average ²²⁸Ra-derived fluxes of 15–259 L m^− 2 d^− 1. Seasonal changes in Ra-derived SGD were apparent in all ponds as well as between ponds, with SGD values of 30–472 L m^− 2 d^− 1 (Winnapaug Pond), 6–20 L m^− 2 d^− 1 (Quonochontaug Pond), 36–273 L m^− 2 d^− 1 (Ninigret Pond), 29–76 L m^− 2 d^− 1 (Green Hill Pond), and 19–83 L m^− 2 d^− 1 (Pt. Judith–Potter Pond). These Ra-derived fluxes are up to two orders of magnitude higher than results predicted by a numerical model of groundwater flow, estimates of aquifer recharge for the study period, and values published in previous Ra-based SGD studies in Rhode Island. This disparity may result from differences in the type of flow (recirculated seawater versus fresh groundwater) determined using each technique, as well as variability in porewater Ra activity.