Modelling and forecasting long-term dynamics of Western Baltic macrobenthic fauna in relation to climate signals and environmental change

Long-term macrobenthos data from Kiel Bight in the Western Baltic collected between 1968 and 2000 have been correlated with the winter NAO index (North Atlantic Oscillation Index) and other environmental data such as temperature, salinity and oxygen content in the bottom water in order to detect systematic patterns related to so far unexplained abiotic signals in the dynamics of zoobenthic species assemblages. The benthos data come from a cluster of five stations (Süderfahrt/ Millionenviertel) in Kiel Bay. Our investigations concentrated on the macrobenthic dynamics with a focus on the number of species m^− 2 (species richness). Using logarithms and the time series analysis approach of Box/Jenkins (ARIMA modelling, transfer function modelling) it was shown that species richness was strongly influenced by the winter NAO (adjusted for a linear time trend within the 1968-2000 period) and salinity (with a shift/lag of four years). Bootstrapping experiments (i.e. sampling from the error process) and analysis of prediction power (by means of the one- or more-years leaving-out method) showed that the parameter estimates behaved in a stable way, leading to a relatively robust model.     

Developing a Marine Information System by Integrating Existing Ocean Models Using Object-Oriented Technology

Information systems developed for different applications within the environmental domain have common characteristics, which can potentially be abstracted for sharing and reuse of design and software modules. This article presents an approach to designing for reuse by abstracting commonalities in the design of a Marine Information System (MIS), facilitating data management in a prototype operational monitoring, forecasting, and management system for the North Atlantic and the Nordic Seas. A detailed study of the requirements and data analysis was carried out, and Object-Oriented Technology (OOT) is employed to encapsulate abstractions and to promote reuse of code and design. This article identifies the Object-Oriented Frameworks (OOFW) required to build the MIS. It also provides guidelines to environmental scientists for restructuring legacy software and employing modern programming techniques.     

A new form of generalized Boussinesq equations for varying water depth

A new set of equations of motion for wave propagation in water with varying depth is derived in this study. The equations expressed by the velocity potentials and the wave surface elevations include first-order non-linearity of waves and have the same dispersion characteristic to the extended Boussinesq equations. Compared to the extended Boussinesq equations, the equations have only two unknown scalars and do not contain spatial derivatives with an order higher than 2. The wave equations are solved by a finite element method. Fourth-order predictor–corrector method is applied in the time integration and a damping layer is applied at the open boundary for absorbing the outgoing waves. The model is applied to several examples of wave propagation in variable water depth. The computational results are compared with experimental data and other numerical results available in literature. The comparison demonstrates that the new form of the equations is capable of calculating wave transformation from relative deep water to shallow water.     

有风浪太湖水体的光照传输模式

根据辐射传输方程的近似解方法和光照传输的累加原理,建立了包括水面风浪效应、水中悬浮粒多次散射和底泥反射的垂直分层水体的光照传输数值模式。利用该模式和太湖悬浮粒实测资料,计算和讨论了太湖水体中光照随水深的变化以及悬浮粒浓度、风浪和太阳高度的影响。结果表明,悬浮粒的多次散射所形成的漫射光是湖水中光照不可忽视的重要部分。在水体的下层,漫射光甚至是光照的主要构成部分。风浪对水中光照的影响则主要在天顶角大于70°时才明显地表现出来,此时风浪水体对太阳直接辐射和全辐射的反射率比平静水面都有明显的减小。     

Evaluation of the Risk of Marine Slope Instability: A Pseudo-3D Approach for Application to Large Areas

This article presents a methodology developed to evaluate the instability of submarine slopes that extend over a large area. Special attention was paid to (1) the complex geometry (bathymetry) and the expanse of the slope, (2) the heterogeneity of the sediment, and (3) the distribution of the pore pressure. The safety factor was considered as a spatially varying quantity. The General Formulation (GLE, Fredlund and Krahn 1977), which fully satisfies equilibrium conditions, was used for evaluating the stability of the marine slope. The submarine slope failure, which occurred on 16 October 1979 during the construction of the new Nice airport, was studied in order to test the developed model. Geotechnical parameters were taken from experimental tests carried out by IFREMER on 19 cores extracted at different depths (from 27 m to 1300 m) (Cochonat, Bourillet, and Savoye, 1993; Mulder et al., 1994). Many scenarios were proposed in order to explain the cause of the Nice slope failure (Habib, 1994). In this article, two of those scenarios were tested. Simulation results are presented and discussed.     

Geoacoustic Inversion via Genetic Algorithm and Its Application to Manganese Sediment Identification

An acoustic inversion method using a wide-band signal and two near field receivers is proposed and applied to multiple layered seabed models including a manganese sediment. The inversion problem can be formulated into a probabilistic model comprised of signals, a forward model, and additive noise. The forward model simulates wide-band signals, such as chirp signals, and is chosen to be the source-waveletconvolution plane wave modeling method. The wavelet matching technique, using weighted least-squares fitting, estimates the sediment sound-speed and thickness on which determination of the possible numerical ranges for a priori uniform distribution is based. The genetic algorithm is applied to a global optimization problem to find a maximum a posteriori solution for determined a priori search space. Here the object function is defined by an L 2 norm of the difference between measured and modeled signals. Not only the marginal pdf but also its statistics are calculated by numerical evaluation of integrals using the samples selected during importance sampling process of the genetic algorithm.     

Numerical simulation of solitary wave scattering by a circular cylinder array

The interaction of a solitary wave with an array of surface-piercing vertical circular cylinders is investigated numerically. The wave motion is modeled by a set of generalized Boussinesq equations. The governing equations are discretized using a finite element method. The numerical model is validated against the experimental data of solitary wave reflection from a vertical wall and solitary wave scattering by a vertical circular cylinder respectively. The predicted wave surface elevation and the wave forces on the cylinder agree well with the experimental data. The numerical model is then employed to study solitary wave scattering by arrays of two circular cylinders and four circular cylinders respectively. The effect of wave direction on the wave forces and the wave runup on the cylinders is quantified.     

The buckling of a corrugated carbon fibre cylinder under external hydrostatic pressure

A theoretical and an experimental investigation was carried out, where a carbon fibre corrugated circular cylinder was tested to destuction under external hydrostatic pressure. The theoretical investigation was via the finite element method, where the structure was modelled with several orthotropic axisymmetric thin-walled shell elements. The experimental observations were aided with strategically placed strain gauges. Comparison between theory and experiment showed that the experimentally observed buckling pressure was a little lower than the theoretical prediction. This may have been due to the fact that the model had slight initial geometrical imperfections in the circumferenential direction.     

A new efficient 3D non-hydrostatic free-surface flow model for simulating water wave motions

A new three-dimensional, non-hydrostatic free surface flow model is presented. For simulating water wave motions over uneven bottoms, the model employs an explicit project method on a Cartesian the staggered gird system to solve the complete three-dimensional Navier–Stokes equations. A bi-conjugated gradient method with a pre-conditioning procedure is used to solve the resulting matrix system. The model is capable of resolving non-hydrostatic pressure by incorporating the integral method of the top-layer pressure treatment, and predicting wave propagation and interaction over irregular bottom by including a partial bottom-cell treatment. Four examples of surface wave propagation are used to demonstrate the capability of the model. Using a small of vertical layers (e.g. 2–3 layers), it is shown that the model could effectively and accurately resolve wave shoaling, non-linearity, dispersion, fission, refraction, and diffraction phenomena.     

Remote estimation of surficial seafloor properties through the application Angular Range Analysis to multibeam sonar data

The variation of the backscatter strength with the angle of incidence is an intrinsic property of the seafloor, which can be used in methods for acoustic seafloor characterization. Although multibeam sonars acquire backscatter over a wide range of incidence angles, the angular information is normally neglected during standard backscatter processing and mosaicking. An approach called Angular Range Analysis has been developed to preserve the backscatter angular information, and use it for remote estimation of seafloor properties. Angular Range Analysis starts with the beam-by-beam time-series of acoustic backscatter provided by the multibeam sonar and then corrects the backscatter for seafloor slope, beam pattern, time varying and angle varying gains, and area of insonification. Subsequently a series of parameters are calculated from the stacking of consecutive time series over a spatial scale that approximates half of the swath width. Based on these calculated parameters and the inversion of an acoustic backscatter model, we estimate the acoustic impedance and the roughness of the insonified area on the seafloor. In the process of this inversion, the behavior of the model parameters is constrained by established inter-property relationships. The approach has been tested using a 300 kHz Simrad EM3000 multibeam sonar in Little Bay, NH. Impedance estimates are compared to in situ measurements of sound speed. The comparison shows a very good correlation, indicating the potential of this approach for robust seafloor characterization.