A rational approach to intact ship stability assessment

In this study, asymptotic and total stability of the non-linear free and forced pure rolling motions of a ship are investigated. A ship performing a rolling motion is taken as a dynamical system. Lyapunov's direct method is employed in the analysis. By generating a time-invariant Lyapunov function, conditions and the domain of asymptotic stability are obtained for free rolling motion. Results of the work on “boundedness” and “uniform boundedness” of the solutions of the equation of forced rolling motion, done by Özkan (1977), that is, conditions of total (practical) stability and its domain in the phase-plane are given and illustrated.     

Total (practical) stability of ships

In this paper the “total (practical) stability” concept is introduced to nonlinear forced rolling motion of a ship. This is achieved by employing “boundedness” and “Lyapunov's function” approach. In this respect two new theorems are proved and conditions and domain of Practical Stability are evaluated. The Paper also contains a critical review of the present status of international intact ship stability regulations. A qualitative discussion of oscillatory rolling motion and the capsizing phenomena is presented.     

An experimental method for determining the frequency-dependent added mass and added mass moment of inertia for a floating body in heave and pitch motions

Most of the existing relevant materials have been obtained from experiments, in which evaluating the added mass at the resonant frequency corresponding to the peak of a frequency-response curve obtained from the “forced” vibration analysis is the most popular technique. In this paper, a simple experimental method was presented where the “free” vibration responses instead of the “forced” ones were used to determine the values of m_ah and I_ap. The main part of the experimental system is composed of a floating body (model) and a spring–shaft shaker. The “free” vibration of this main part was induced by imposing on it an initial displacement (and/or an initial velocity), and from the time histories of displacements information such as the “damped” natural frequencies, damping ratios, sectional added mass coefficients (C_V and C_P) were obtained. Since the displacements of the spring–shaft shaker are “translational” and those of the floating body due to pitch motions are “angular”, a technique for the transformation between the associated parameters of the two components of the main part was presented.     

Planing and impacting plate forces at large trim angles

Added mass theory has been shown to give excellent agreement with experimental measurements on planing surfaces at normal planing angles [e.g. Payne, P.R. (1982, Ocean Engng9, 515–545; 1988, Design of High-Speed Boats, Volume 1: Planning. Fishergate. Inc., Annapolis, Maryland)] and to agree exactly with more complex conformal transformations where such a comparison is possible. But at large trim angles, it predicts non-transient pressures that are greater than the free-stream dynamic pressure and so cannot be correct. In this paper, I suggest that the reason is because, unlike a body or a wing in an infinite fluid, a planing plate only has fluid on one side—the “high pressure” side. So the fluid in contact with the plate travels more slowly as the plate trim angle (and therefore static pressure) increases. This results in lower added mass forces than Munk, M. (1924) The Aerodynamic Forces on Airship Hulls (NACA TR-184) and Jones, R. T. (1946) Properties of Low-Aspect-Ratio Pointed Wings at Speeds Below and Above the Speed of Sound (NACA TR-835) originally calculated for wings and other bodies in an infinite fluid.For simplicity of presentation, I have initially considered the example of a triangular (vertex forward) planning plate. This makes the integration of elemental force very simple and so the various points are made without much trouble. But the penalty is that there seem to be no experimental data for such a configuration; at least none that I have been able to discover. But at least the equations obtained in the limits of zero and infinite aspect ratio, small trim angles (τ) and τ = 90° all agree with established concepts and the variation of normal force with trim angle looks like what we would expect from our knowledge of how delta wings behave in air.I then employed the new equation to calculate the force on a rectangular planing surface at a trim angle τ, having a constant horizontal velocity u_o and a vertical impact velocity of ż. This happens to have been explored experimentally by Smiley, R. F. [(1951) An Experimental Study of Water Pressure-Distributions During Landings and Planing of a Heavily Loaded Rectangular Flat Plate Model (NACA TM 2453)] up to trim angles of τ = 45°, and so a comparison between theory and experiment is possible. The results of this comparison are encouraging, as is also a comparison with the large trim angle planing plate measurements of Shuford, C. L. [(1958) A Theoretical and Experimental Study of Planing Surfaces Including Effects of Cross Section and Plan Form (NACA Report)].As two practical applications, I first employed the new equations to calculate the “design pressures” needed to size the plating of a transom bow on a high-speed “Wavestrider” hull. The resulting pressures were significantly different to those obtained using semi-empirical design rules in the literature. Then I used the theory to critically review data obtained from tank tests of a SES bow section during water impact to identify how the “real world” of resilient deck plating diverged from the “model world” of extreme structural rigidity.     

Multi-vessel seakeeping computations with linear potential theory

The wave-induced motion of two vessels in close proximity is studied using traditional two-dimensional (2D) and three-dimensional (3D) boundary element methods. Added mass, damping, and the behavior of the free surface between the vessels are examined in some detail. As expected based on previous work, the response of the water between the vessels is found to have a profound effect on the hydrodynamic forces. At “critical frequencies” corresponding to standing waves between the hulls, the hydrodynamic forces undergo significant drastic changes. The 2D and 3D results are compared, and the effects of a skew angle between the vessels are examined. Some of the consequences of the behavior at the critical frequencies for simulations in the time domain are examined, the most significant of which is a very lightly damped impulse response.     

Towing dynamics of a liferaft and fast rescue craft in a surface wave

The equations of motion for the coupled dynamics of a small liferaft and fast rescue craft in a surface wave are formulated in two dimensions using the methods of Kane and Levinson [1985. Dynamics: Theory and Applications. McGraw-Hill Inc., New York]. It is assumed that the motion normal to the wave surface is small and can be neglected, i.e. the bodies move along the propagating wave profile. The bodies are small so that wave diffraction and reflection are negligible. A Stokes second order wave is used and the wave forces are applied using Morison's equation for a body in accelerated flow. Wind loads are similarly modelled using drag coefficients. The equations are solved numerically using the Runge–Kutta routine “ode45” of MATLAB^®. The numerical model provides guidelines for predicting the tow loads and motions of small craft in severe sea states.     

A new view of near-shore dynamics based on observations from HF Radar

Since 1984 the OSCR HF Radar system has been used in over 50 deployments to measure near-shore surface currents for both scientific and engineering applications. The enhanced scope, resolution and accuracy of these measurements have yielded new insights into the tidal, wind and density driven dynamics of the near-shore zone.Tidal current ellipses obtained from these radar measurements have been shown to be in good aggrement with values calculated by numerical models both for the predominant constituents and also for higher harmonics. Coherent patterns of wind-forced currents ahve been determined with strong evidence of a “slab-like” surface response. In one deployment, with offshore winds blowing over relatively deep water, this “slab” rotated clockwise at near-inertial frequency. Strong (up to 20cm s⁻¹), persistent surface residual currents are commonly observed, these are almost certainly generated by (small) horizontal density gradients. These observed surface residuals provide ideal data for rigorous testing of 3-D numerical models.With a threatened rise in sea level, HF Radar is well-suited for observing the expected changes in the dynamics of near-shore regions. Continuing development of these radar systems offers exciting prospects of remote sensing of both surface waves and currents. Future applications may extend beyond the near-shore region to measurements along the shelf-edge, in oceanic gyres and for “beach-processes”.     

On extracting current profiles from vertically integrated numerical models

A transformation method is presented by which current profiles (of tidal or wind-induced origin) can be extracted at any horizontal position and moment in time from a vertically integrated, two-dimensional, hydrodynamic numerical model. An arbitrary vertical variation of eddy viscosity can be included in the method, which can incorporate a no-slip bottom boundary condition. The technique assumes that the sea is homogeneous.The method is used to improve the representation of bottom stress within the two-dimensional model, whereby the bottom stress is no longer related simply to the depth-mean current as in the “conventional” two-dimensional, vertically integrated model.Idealized calculations for a range of eddy viscosity profiles, show that elevations, current profiles, and time series of current extracted from this “enhanced” two-dimensional numerical model are in good agreement with currents obtained from a full three-dimensional model.     

The differences between a wing and a planing plate in two-dimensional flow

The normal force coefficient on a flat planing surface having arbitrary heave and pitch motion in two-dimensional flow is compared with the lift coefficient of a thin wing in an infinite fluid. Despite the totally different derivations, they are found to be identical (at large Froude numbers and low trim angles and allowing for the wing's interaction with twice as much fluid) at low reduced frequencies. For higher frequency motions, the wing's angle of attack induced lift and its pitch and heave damping are less than those of a planing surface, but the acceleration terms remain identical. The differences at the higher reduced frequencies are due to the fact that, in invisad irrotational flow, the planning plate cannot leave a vortex wake, whereas a wing does.It seems to follow that the “virtual mass” planing hull analysis can be applied to “quasi-static” problems involving wings and bodies in an infinite fluid without the slenderness restriction originally imposed by Jones (1946). Certainly, it is remarkable that the so called “quasi-steady” forces on a two-dimensional wing can be obtained in a few lines of elementary analysis. On the other hand, the method fails entirely when used to compute the pitching moment on a two-dimensional plate, even though it has been found to give good results for the three-dimensional case (Payne, 1981c).This work is offered as a very incomplete study of an intriguing relationship between two very different bodies of analysis. Much more work will need to be done before the relationship between the two approaches will be fully understood.     

The discrete methods for free vibration analyses of an immersed beam carrying an eccentric tip mass with rotary inertia

In this paper, a beam without contact with water is called the “dry” beam and the one in contact with water is called the “wet” beam. For a partially (or completely) immersed uniform beam carrying an eccentric tip mass possessing rotary inertia, the conventional analytical (closed-form) solution is achieved by considering the inertial forces and moments of the tip mass and rotary inertia as the boundary conditions at the tip end of the beam. However, it has been found that the approximate solution for the last problem may be achieved by two techniques: Method 1 and Method 2. In Method 1, the basic concept is the same as the conventional analytical method; but in Method 2, the tip end of the beam is considered as a free end, while the inertial forces and moments induced by the tip mass and rotary inertia are considered as the external loads applied at the tip end of the beam. The main differences between the formulation of Method 1 and that of Method 2 are: In Method 1, the “normal” shapes of the “dry” beam are functions of the frequency-dependent boundary conditions but the external loads at the tip end are equal to zero; On the contrary, in Method 2, the “normal” mode shapes of the “dry” beam are determined based on the zero boundary conditions at the tip end of the beam but the external loads at the tip end due to the inertial effects of the tip mass and rotary inertia must be taken into consideration for the free vibration analysis of the “wet” beam. Numerical results reveal that the approximate solution obtained from Method 2 are very close to that from Method 1 if the tip mass moment of inertia is negligible. Besides, the two approximate solutions are also very close to the associated analytical (closed-form) solution or the finite element solution. In general, it is hoped that there exist several methods for tackling the same problem so that one may have more choices to incorporate with the specified cases. It is believed that the two approximate methods presented in this paper will be significant from this point of view.     

Forced vibration analysis of an offshore tower carrying an eccentric tip mass with rotary inertia due to support excitation

Dynamic responses of structures due to earthquake excitation are the important problems in engineering, thus, the information concerned is plenty. However, most of the literature is relating to the discrete methods, particularly to the finite element method (FEM), and the one relating to the method combining both the “continuous” and “discrete” models is rare. The objective of this paper is to provide some information in this respect. First, the analytical solution for the natural frequencies and normal mode shapes of a “continuous” tower, without contacting water (or “dry” tower), carrying an eccentric tip mass possessing rotary inertia is determined. Next, the partial differential equation of motion for the forced vibration of the tower, contacting water (or “wet” tower), subjected to support excitation is transformed into a matrix equation by using the last natural frequencies and normal modes shape of the freely vibrating dry tower. Finally, the numerical integration method is used to solve the matrix equation to yield the seismic response of the wet tower. In theory, the mode superposition method is correct only if the total number of modes considered approaches infinity, however, numerical results of this paper reveal that superposition of only the lowest six modes will yield excellent results to be very close to the corresponding ones obtained from the conventional FEM. For this reason, the CPU time required by the presented approach is less than 5% of that required by the conventional FEM.     

An efficient algorithm for strapdown accelerometer-based motion measurement

The paper is concerned with the measurement of motion of free floating bodies using accelerometers strapped to a cross-arm mounted on the body. The measurement system has been studied with the help of a formulation involving a set of coupled nonlinear initial-value equations involving the angular acceleration components. A PC-based software using the Bulirsch-Stoer technique has been developed to solve the initial-value problem so as to deduce the body motions from the measured accelerations. Suitable filtering strategy has been employed at every stage of numerical integration. The reliability of the strapdown accelerometer system together with software developed has been validated using a “dry” test. Typical motion measurements have been done in all the six degrees of freedom of a tug model in a wave flume. The method is deemed to be an efficient and cost-effective technique suitable for free floating bodies and for large motions.     

The BCI criterion for the initiation of breaking process in Boussinesq-type equations wave models

The Breaking Celerity Index (BCI) is proposed as a new wave breaking criterion for Boussinesq-type equations wave propagation models (BTE).The BCI effectiveness in determining the breaking initiation location has been verified against data from different experimental investigations conducted with incident regular and irregular waves propagating along uniform slope [Utku, M. (1999). “The Relative Trough Froude Number. A New Criteria for Wave Breaking”. Ph.D. Dissertation, Dept. of Civil and Enviromental Engineering, Old Dominion University, Norfolk, VA; Gonsalves Veloso dos Reis, M.T.L. (1992). “Characteristics of waves in the surf zone”. MS Thesis, Department of Civil Engineering, University of Liverpool., Liverpool; Lara, J.L., Losada, I.J., and Liu, P.L.-F. (2006). “Breaking waves over a mild gravel slope: experimental and numerical analysis”. Journal of Geophysical Research, VOL 111, C11019] and barred beaches [Tomasicchio, G.R., and Sancho, F. (2002). “On wave induced undertow at a barred beach”. Proceedings of 28th International Conference on Coastal Engineering, ASCE, New York, 557–569]. The considered experiments were carried out in small-scale and large-scale facilities. In addition, one set of data has been obtained by the use of the COBRAS model based upon the Reynolds Averaged Navier Stokes (RANS) equations [Liu, P.L.-F., Lin, P., Hsu, T., Chang, K., Losada, I.J., Vidal, C., and Sakakiyama, T. (2000). “A Reynolds averaged Navier–Stokes equation model for nonlinear water wave and structure interactions”. Proceedings of Coastal Structures ‘99, Balkema, Rotterdam, 169–174; Losada, I.J., Lara, J.L., and Liu, P.L.-F. (2005). “Numerical simulation based on a RANS model of wave groups on an impermeable slope”. Proceedings of Fifth International Symposium WAVES 2005, Madrid].Numerical simulations have been performed with the 1D-FUNWAVE model [Kirby, J.T., Wei, G., Chen, Q., Kennedy, A.B., and Dalrymple, R.A. (1998). “FUNWAVE 1.0 Fully Nonlinear Boussinesq Wave Model Documentation and User's Manual”. Research Report No CACR-98-06, Center for Applied Coastal Research, University of Delaware, Newark]. With regard to the adopted experimental conditions, the breaking location has been calculated for different trigger mechanisms [Zelt, J.A. (1991). “The run-up of nonbreaking and breaking solitary waves”. Coastal Engineering, 15, 205–246; Kennedy, A.B., Chen, Q., Kirby, J.T., and Dalrymple, R.A. (2000). “Boussinesq modeling of wave transformation, breaking and run-up. I: 1D”. Journal of Waterway, Port, Coastal and Ocean Engineering, 126, 39–47; Utku, M., and Basco, D.R. (2002). “A new criteria for wave breaking based on the Relative Trough Froude Number”. Proceedings of 28th International Conference on Coastal Engineering, ASCE, New York, 258–268] including the proposed BCI.The calculations have shown that BCI gives a better agreement with the physical data with respect to the other trigger criteria, both for spilling and plunging breaking events, with a not negligible reduction of the calculation time.     

Tracking the last sea-level cycle: seafloor morphology and shallow stratigraphy of the latest Quaternary New Jersey middle continental shelf

Seafloor geomorphology and surficial stratigraphy of the New Jersey middle continental shelf provide a detailed record of sea-level change during the last advance and retreat of the Laurentide ice sheet (120 kyr B.P. to Present). A NW–SE-oriented corridor on the middle shelf between water depths of 40 m (the mid-shelf “paleo-shore”) and 100 m (the Franklin “paleo-shore”) encompasses 500 line-km of 2D Huntec boomer profiles (500–3500 Hz), an embedded 4.6 km² 3D volume, and a 490 km² swath bathymetry map. We use these data to develop a relative stratigraphy. Core samples from published studies also provide some chronological and sedimentological constraints on the upper <5 m of the stratigraphic succession.The following stratigraphic units and surfaces occur (from bottom to top): (1) “R”, a high-amplitude reflection that separates sediment >46.5 kyr old (by AMS ¹⁴C dating) from overlying sediment wedges; (2) the outer shelf wedge, a marine unit up to 50 m thick that onlaps “R”; (3) “Channels”, a reflection sub-parallel to the seafloor that incises “R”, and appears as a dendritic system of channels in map view; (4) “Channels” fill, the upper portion of which is sampled and known to represent deepening-upward marine sediments 12.3 kyr in age; (5) the “T” horizon, a seismically discontinuous surface that caps “Channels” fill; (6) oblique ridge deposits, coarse-grained shelly units comprised of km-scale, shallow shelf bedforms; and (7) ribbon-floored swales, bathymetric depressions parallel to modern shelf currents that truncate the oblique ridges and cut into surficial deposits.We interpret this succession of features in light of a global eustatic sea-level curve and the consequent migration of the coastline across the middle shelf during the last 120 kyr. The morphology of the New Jersey middle shelf shows a discrete sequence of stratigraphic elements, and reflects the pulsed episodicity of the last sea-level cycle. “R” is a complicated marine/non-marine erosional surface formed during the last regression, while the outer shelf wedge represents a shelf wedge emplaced during a minor glacial retreat before maximum Wisconsin lowstand (i.e., marine oxygen isotope stage 3.1). “Channels” is a widespread fluvial subarial erosion surface formed at the late Wisconsin glacial maximum 22 kyr B.P. The shoreline migrated back across the mid-shelf corridor non-uniformly during the period represented by “Channels” fill. Oblique ridges are relict features on the New Jersey middle shelf, while the ribbon-floored swales represent modern shelf erosion. There is no systematic relationship between modern seafloor morphology and the very shallowly buried stratigraphic succession.     

Evaluating the structural condition of synthetic and metallic cables

In view of the fact that mechanical cables are very popular and powerful structural members, methods that provide information concerning their “health” are anxiously needed. In general, and because of their extended length (a rope used in a gold mine in South Africa is 9.3 miles long), methods which inspect the cable as a continuous function of the length variable are very expensive and time-consuming. The present paper surveys the state-of-the-art and discusses limitations and advantages of the techniques available today, the most recent one being the “induced, transverse wave propagation method”.     

Hydrothermal plumes in the eastern Manus Basin, Bismarck Sea: CH4, Mn, Al and pH anomalies

We present evidence for strong hydrothermal activity in the eastern Manus Basin (depth: 1700–2100 m), the existence of large scale triple-layered buoyant plumes at depths of 1100 m (“shallow plume”), 1700 m (“deep plume”), and 1400 m (“middle plume” with less extent than the other two plumes) that were revealed from water column anomalies of CH₄, Mn, Al and pH observed in November to December 1990. Judging from the horizontal distribution of these parameters, the deep plume seems to originate from two distinct hydrothermal sites (eastern and western sites) in the research area, the eastern site being visually ascertained with deep-tow observations at the same time. The CH₄/Mn ratio (mol mol⁻¹) of the deep plume (0.02–0.05) is the lowest yet observed in hydrothermal plumes. The order of magnitude difference of CH₄/Mn ratios between the shallow plume and the deep plume suggests that different kinds of fluid-rock interaction occurred to make the hydrothermal end members for the deep and shallow plumes. The shallow plume, which had an areal extent of more than 50 km, may be an episodic “megaplume”, because it was not recognized in the previous CH₄ profiles in 1986, and because it has a similar CH₄/Mn ratio as the megaplume observed in the North Fiji Basin. We found that the eastern deep plume is characterized by enormously high aluminium concentrations (0.6– 1.5 μmol kg⁻¹), pH anomalies (0.1) and high Al/Mn ratios (10–17). The endmember fluid for the eastern deep plume may have an unusually low pH value to dissolve this much aluminum during fluid-rock interaction, or this plume may originate from an eruption-influenced fluid.     

Geometry and architectural associations of co-genetic debrite–turbidite beds in basin-margin strata, Carboniferous Ross Sandstone (Ireland): Applications to reservoirs located on the margins of structurally confined submarine fans

Co-genetic debrite–turbidite beds are most commonly found in distal basin-plain settings and basin margins. This study documents the geometry, architectural association and paleogeographic occurrence of co-genetic debrite–turbidite beds in the Carboniferous Ross Sandstone with the goal of reducing uncertainty in the interpretation of subsurface data in similarly shaped basins where oil and gas is produced.The Ross Sandstone of western Ireland was deposited in a structurally confined submarine basin. Two outcrops contain co-genetic debrite–turbidite beds: Ballybunnion and Inishcorker. Both of the exposures contain strata deposited on the margin of the basin. An integrated dataset was used to characterize the stratigraphy of the Ballybunnion exposure. The exposure is divided into lower, middle, and upper units. The lower unit contains laminated shale with phosphate nodules, structureless siltstone, convolute bedding/slumps, locally contorted shale, and siltstone turbidites. The middle unit contains co-genetic debrite–turbidite beds, siltstone turbidites, and structureless siltstone. Each co-genetic debrite–turbidite bed contains evidence that fluid turbulence and matrix strength operated alternately and possibly simultaneously during deposition by a single sediment-gravity-flow event. The upper unit contains thin-bedded sandy turbidites, amalgamated sandy turbidites, siltstone turbidites, structureless siltstone, and laminated shale. A similar vertical facies pattern is found at Inishcorker.Co-genetic debrite–turbidite beds are only found at the basin-margin. We interpret these distinct beds to have originated as sand-rich, fully turbulent flows that eroded muddy strata on the slope as well as interbedded sandstone and mudstone in axial positions of the basin floor forming channels and associated megaflute erosional surfaces. This erosion caused the axially dispersing flows to laterally evolve to silt- and clay-rich flows suspended by both fluid turbulence and matrix strength due to a relative increase in clay proportions and associated turbulence suppression. The flows were efficient enough to bypass the basin center/floor, physically disconnecting their deposits from coeval lobes, resulting in deposition of co-genetic debrite–turbidite beds on the basin margin. The record of these bypassing flows in axial positions of the basin is erosional surfaces draped by thin siltstone beds with organic debris.A detailed cross-section through the Ross Sandstone reveals a wedge of low net-to-gross, poor reservoir-quality strata that physically separates sandy, basin-floor strata from the basin margin. The wedge of strata is referred to as the transition zone. The transition zone is composed of co-genetic debrite–turbidite beds, structureless siltstone, slumps, locally contorted shale, and laminated shale. Using data from the Ross Sandstone, two equations are defined that predict the size and shape of the transition zone. The equations use three variables (thickness of basin-margin strata, thickness of coeval strata on the basin floor, and angle of the basin margin) to solve for width (w) and trajectory of the basinward side of the low net-to-gross wedge (β). Beta is not a time line, but a facies boundary that separates sandy basin floor strata from silty basin-margin strata. The transition zone is interpreted to exist on lateral and distal margins of the structurally confined basin.Seismic examples from Gulf of Mexico minibasins reveal a wedge of low continuity, low amplitude seismic facies adjacent to the basin margin. Strata in this wedge are interpreted as transition-zone sediments, similar to those in the Ross Sandstone. Besides defining the size and shape of the transition zone, the variables “w” and “β” define two important drilling parameters. The variable “w” corresponds to the minimum distance a well bore should be positioned from the lateral basin margin to intersect sandy strata, and “β” corresponds to the deviation (from horizontal) of the well bore to follow the interface between sandy and low net-to-gross strata. Calculations reveal that “w” and “β” are related to the relative amount of draping, condensed strata on the margin and the angle of the basin margin. Basins with shallowly dipping margins and relatively high proportions of draping, clay-rich strata have wider transition zones compared to basins with steeply dipping margins with little draping strata. These concepts can reduce uncertainty when interpreting subsurface data in other structurally confined basins including those in Gulf of Mexico, offshore West Africa, and Brunei.     

Hydrological heterogeneity, nutrient dynamics and water quality of a non-tidal lentic ecosystem (Lesina Lagoon, Italy)

The dynamics of the Lesina coastal lagoon (Italy) in terms of nutrients, phytoplankton and chemical–physical parameters were evaluated, together with their functional relationships with freshwater inputs, in order to identify ecosystem responses to changes in driving forces in a Mediterranean non-tidal lentic environment. Lesina Lagoon is a shallow coastal environment characterised by limited exchange with coastal waters, which favours enrichment of nutrients and organic matter and benthic fluxes within the system. Lagoon–sea exchanges are influenced by human management. There is a steep salinity gradient from East to West. High nitrogen and silica values were found close to freshwater inputs, indicating wastewater discharges and agricultural runoff, especially in winter. Dissolved oxygen was well below saturation (65%) near sewage and runoff inputs in the western part of the lagoon during summer. Classification in accordance with EEA (2001) guidelines suggests the system is of “poor” or “bad” quality in terms of nitrogen concentrations in the eastern zone during the winter rainy period. In terms of phosphate concentrations, the majority of the stations fall into the “good” category, with only two stations (close to the sewage and runoff inputs) classed as “bad”. In both cases, the raw nitrogen levels make the lagoon a P-limited system, especially in the eastern part. There was wide space–time variability in chlorophyll a concentrations, which ranged from 0.25 to 56 μg l⁻¹. No relationships between chlorophyll a and nutrients were found, suggesting that autotrophic biomass may be controlled by a large number of internal and external forcing factors driving eutrophication processes. Water quality for this type of environment depends heavily on pressure from human activities but also on the management of sewage treatment plants, agricultural practices and the channels connecting the lagoon with the sea.     

A metamodel-based towed system simulation

The objective of this exercise was the development of a computer simulation of a US Navy towed system (towed body and line array) capable of near real-time predictive capabilities. The purpose of the computer simulation was to be able to determine the lateral motion of the towed body in response to ship maneuvers. The theoretical approach was to develop a spatial-domain computer simulation of each of the towed system components treated as a system, and then develop a “metamodel” of the response by using simple algebraic expressions. At-sea trials were conducted aboard US Navy vessels to verify the simulation and fine-tune its performance.