Developing a conceptual model of marine farming in New Zealand

Survey and Geographic Information System (GIS) data analysis describes the relative influence of biophysical and human variables on site choices made by marine farmers in New Zealand. Community conflicts have grown in importance in determining farm location and different government planning strategies leave distinct signature patterns. Recent legislation empowers local governments to choose among three strategies for future regional aquaculture development. This paper suggests each strategy could result in different spatial outcomes. Simulation modelling of the type described here can provide a better understanding of farmer responses to management approaches and the range of futures that could result from planning choices made today.     

A temporally second-order accurate Godunov-type scheme for solving the extended Boussinesq equations

A numerical scheme for solving the class of extended Boussinesq equations is presented. Unlike previous schemes, where the governing equations are integrated through time using a fourth-order method, a second-order Godunov-type scheme is used thus saving storage and computational resources. The spatial derivatives are discretised using a combination of finite-volume and finite-difference methods. A fourth-order MUSCL reconstruction technique is used to compute the values at the cell interfaces for use in the local Riemann problems, whilst the bed source and dispersion terms are discretised using centred finite-differences of up to fourth-order accuracy. Numerical results show that the class of extended Boussinesq equations can be accurately solved without the need for a fourth-order time discretisation, thus improving the computational speed of Boussinesq-type numerical models. The numerical scheme has been applied to model a number of standard test cases for the extended Boussinesq equations and comparisons made to physical wave flume experiments.     

Influence of discrete step size on wind setup component of storm surge

Open coast storm surge water levels consist of wind setup due to wind shear at the water surface; a wave setup component caused by wind induced waves transferring momentum to the water column; an atmospheric pressure head component due to the atmospheric pressure deficit over the spatial extent of the storm system; a Coriolis forced setup or setdown component due to the effects of the rotation of the earth acting on the wind driven alongshore current at the coast; a possible seiche component due to resonance effects initiated by moving wind system, and, if astronomical tides are present, an astronomical tide component (although the tide is typically considered to be a forced astronomical event and not really a direct part of the external wind-driven meteorological component of storm surge). Typically the most important component of a storm surge is the wind setup component, especially on the U.S. East Coast and the Gulf of Mexico shorelines. In many approaches to storm surge modeling, a constant depth approximation is invoked over a limited step size in the computational domain. The use of a constant depth approximation has received little attention in the literature although can be very important to the resulting magnitude of the computed storm surge. The importance of discrete step size to the wind setup storm surge component is considered herein with a simple case computation of the wind setup component on a linear slope offshore profile. The present study findings show that the constant depth approximation to wind setup storm surge estimation is biased on the low side (except in extremely shallow water depths) and can provide large errors if discrete step size is not sufficiently resolved. Guidance has been provided on the error that one might encounter for various step sizes on different slopes.     

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.     

Second-order resonance of sloshing in a tank

Sloshing in a two-dimensional rectangular tank in horizontal motion is analysed based on the velocity potential theory. It is found that even when the first-order excitation is away from all the natural frequencies of the tank, second-order resonance can still occur when the sum-frequency or the difference-frequency is equal to one of the natural frequencies corresponding to the even mode. However, such resonance is not excited when the sum or difference frequency is equal to the natural frequency of an odd mode.     

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.     

A B-spline based higher order panel method for analysis of steady flow around marine propellers

A higher order panel method based on B-spline representation for both the geometry and the solution is developed for the analysis of steady flow around marine propellers. The self-influence functions due to the normal dipole and the source are desingularized through the quadratic transformation, and then shown to be evaluated using conventional numerical quadrature. By selecting a proper order for numerical quadrature, the accuracy of the present method can be increased to the machine limit. The far- and near-field influences are shown to be evaluated based on the same far-field approximation, but the near-field solution requires subdividing the panels into smaller subpanels continuously, which can be effectively implemented due to the B-spline representation of the geometry. A null pressure jump Kutta condition at the trailing edge is found to be effective in stabilizing the solution process and in predicting the correct solution. Numerical experiments indicate that the present method is robust and predicts the pressure distribution on the blade surface, including very close to the tip and trailing edge regions, with far fewer panels than existing low-order panel methods.     

Spectrophotometric determination of pH in seawater off Taiaroa Head,Otago, New Zealand: Full-spectrum modelling and prediction of pCO2 levels

Accurate measurement of seawater pH has long been sought by marine chemists (for example: [Dickson, A.G. 1993a. The measurement of sea water pH. Marine Chemistry, 44, 131–142, Dickson, A.G. 1993b. pH buffers for sea water media based on the total hydrogen ion concentration scale. Deep-Sea Research, 40, 107–118; Zhang, 1996; Tapp, M., Hunter, K.A., Currie, K. and Macaskill, B. 2000. Apparatus of continuous-flow underway spectrophotometric measurement of surface water pH. Marine Chemistry 72(2–4), 193–202; Friis, K., Koetzinger, A., Wallace, D.W.R. 2004. Spectrophotometric pH measurement in the ocean: Requirements, design and testing of an autonomous charge-coupled device detector system. Limnology and Oceanography: Methods 2, 126–136]. Recently, such attempts have taken on greater significance as anthropogenic carbon dioxide emissions may create rapidly changing oceanic pH. Spectrophotometric techniques have been accepted generally as the best for determination of seawater pH. Here we report a new technique using thymol blue as the indicator dye and fitting the entire spectrum from 400 to 900 nm rather than measuring the absorbance values at only two or three points in the spectrum. This full-spectrum modelling enables a reduction in signal to noise over other techniques. In the laboratory, we find with seawater samples a pH precision increase of five-fold “within” a sample and seven-fold “between” samples when comparing the full spectrum to the three-point method of analysis [Zhang, H., Byrne, R.H. 1996. Spectrophotometric pH measurements of surface seawater at in-situ conditions: absorbance and protonation behaviour of thymol blue. Marine Chemistry 52, 17–25].     

Assessment of a Sr isotope vital effect (87Sr/86Sr) in marine taxa from Lee Stocking Island, Bahamas

 Data from Lee Stocking Island, Bahamas, confirms the hypothesis that there are no vital effects with the uptake of Sr isotopes (⁸⁷Sr/⁸⁶Sr) at the present mass spectrometer resolution [±2×10^-5 (2σ)]. Our data set contains analyses of 40 samples derived from 37 different calcareous taxa inhabiting a wide range of carbonate subenvironments (i.e., reefal, intertidal, supratidal, mangrove). The mean value of our analyses was 0.709179 with a standard deviation of 2.4×10^-5 (2σ) which was very close to the long-term uncertainty of the strontium isotope methodology [±2.0×10^-5 (2σ)] and to the widely reported ⁸⁷Sr/⁸⁶Sr value of seawater, which clusters around 0.709175. Received: 14 October 1997 / Revision received: 20 March 1998