Swimming for survival: A role of phytoplankton motility in a stratified turbulent environment

We investigate a role for vertical migration in stratified coastal water, where the swimming speed is generally significantly less than the typical turbulent fluctuations in a tidally-mixed bottom layer. In our modelling approach we use a k- turbulence model to describe the physical forcing, a Lagrangian random walk model to describe the vertical displacement of individual cells in response to turbulence and due to cell motility, and a phytoplankton growth model to direct the swimming behaviour of the phytoplankton according to their light and nutrient requirements. The model results show how the cells form a stable subsurface chlorophyll maximum (SCM) at the base of the thermocline where episodic tidal turbulence causes erosion of part of the SCM biomass into the bottom mixed layer (BML). We then focus on the question of whether an ability to swim (weakly, compared to typical bottom layer turbulent intensities) provides any advantage by allowing return to the SCM. Our results show that tidal turbulence in the BML helps both motile and neutrally-buoyant cells by periodically pushing them into the base of the thermocline. Motile cells then have the advantage that they can swim further into the thermocline towards higher light which also reduces the likelihood of being re-mixed back into the BML.     

A numerical investigation of the impact of turbulence on the feeding rates of Oithona davisae

Individual based numerical simulations of the copepod, Oithona davisae, feeding on motile prey, Oxyrrhis marina, under variable turbulent conditions are performed. These simulations correspond to laboratory observations conducted by Saiz et al. [Saiz, E., Calbet, A., and Broglio, E., 2003. Effects of small-scale turbulence on copepods: the case of Oithona Davisae. Limnol. Oceanogr., 48:1304–1311.].The flow field in the simulation is reconstructed by a kinematic simulation whose characteristic scales are derived from the grid mesh and the dissipation rates of the laboratory experiments. The kinematic simulation provides a simplified model, which while not fully realistic, captures the basic relevant feature of turbulence. A hop and sink swimming behaviour is prescribed for O. davisae, while O. marina moves along helical paths with random changes of directions.Three possible effects are tested: the existence of a time threshold in the duration of the contacts between predator and prey, a progressive reduction of the perceptive distance with increasing turbulence level and an abrupt reduction in feeding of O. davisae when the flow speed, in relation to the copepod position, is higher than a prescribed threshold. This last approach introduces an intermittency in the feeding which depends on the variations of velocity both in space and time within the numerical box.The introduction of the time threshold causes a dome-shaped relationship between the simulated enhancement factor and the dissipation rate, while with the other two effects, a monotonic decrease in the enhancement factor is observed, with values reasonably close to the ones observed in the laboratory experiment. In all the cases, the use of realistic values of biological parameters (e.g. swimming behaviour) reproduces response curves in the range of the observations.     

Quantifying copepod kinematics in a laboratory turbulence apparatus

We describe application of a new apparatus that permits simultaneous detailed observations of plankton behavior and turbulent velocities. We are able to acquire 3D trajectories amenable to statistical analyses for comparisons of copepod responses to well-quantified turbulence intensities that match those found in the coastal ocean environment. The turbulence characteristics consist of nearly isotropic and homogeneous velocity fluctuation statistics in the observation region. In the apparatus, three species of copepods, Acartia hudsonica, Temora longicornis, and Calanus finmarchicus were exposed separately to stagnant water plus four sequentially increasing levels of turbulence intensity. Copepod kinematics were quantified via several measures, including transport speed, motility number, net-to-gross displacement ratio, number of escape events, and number of animals phototactically aggregating per minute. The results suggest that these copepods could control their position and movements at low turbulence intensity. At higher turbulence intensity, the copepods movement was dominated by the water motion, although species-specific modifications due to size and swimming mode of the copepod influenced the results. Several trends support a dome-shaped variation of copepod kinematics with increasing turbulence. These species-specific trends and threshold quantities provide a data set for future comparative analyses of copepod responses to turbulence of varying duration as well as intensity.     

Lagrangian modelling of plankton motion: From deceptively simple random walks to Fokker–Planck and back again

The movement of plankton, either by turbulent mixing or their own inherent motility, can be simulated in a Lagrangian framework as a random walk. Validation of random walk simulations is essential. There is a continuum of mathematically valid stochastic integration schemes upon which random walk simulations depend, each of which lead to radically different macro-scale dynamics as expressed in their corresponding Fokker–Planck equations. In addition, diffusivity is not a unique parameter describing a random walk and its corresponding Fokker–Planck equation. Spatially varying translation speed and turn frequency have different effects on population distributions. Validation requires extra information either in the form of the well-mixed condition for physical diffusion, or in detailed information on the sensing ability, internal state modulation and swimming response for plankton motility.     

Effects of small-scale turbulence on the growth of two diatoms of different size in a phosphorus-limited medium

The effect of turbulence on the nutrient flux towards osmotrophic cells is predicted to be size dependent. This should translate into growth. We experimentally followed and modelled the growth of two marine diatoms of different size (Thalassiosira pseudonana, 6 μm in diameter and Coscinodiscus sp., ca. 109 μm in diameter) under still water and turbulent conditions, using a shaker table. Experiments were done with phosphorus-limited cultures and lasted for ca. 5 days. Turbulence enhanced the growth of Coscinodiscus sp. in agreement with theory but not the growth of T. pseudonana, which was actually slightly lower under turbulence. At the end of the experiments there were about 1.7 times as many Coscinodiscus sp. cells in the turbulent treatment than in the still treatment, while for T. pseudonana almost the same cell concentration was found in both conditions. In addition, the Coscinodiscus sp. cells growing under still conditions presented a higher specific alkaline phosphatase activity than those growing in turbulence which indicates a higher need for phosphorus in the still cultures. A simple dynamic model, based on Michaelis–Menten nutrient uptake kinetics, needed nearly no optimisation other than using observed initial conditions of phosphate and cell concentrations. The model showed how an increased nutrient flux towards the cells translates non-linearly into cell growth, most likely by affecting the half-saturation constant (K_M). However, since Coscinodiscus sp. experienced significant mortality and cells partially settled to the bottom of the containers, unequivocal support for the size-dependent effect of turbulence on nutrient uptake will require further experiments and more sophisticated modelling. The mechanisms to connect an increased nutrient flux towards cells with population growth and whether this process is size dependent are important in parameterizing the effects of turbulence on marine plankton in coupled physical–biological models.     

基于脆性熵理论的船舶电力系统脆性评价

将脆性理论引用到船舶电力系统研究中,通过对船舶电力系统脆性概率熵、脆性综合概率熵、脆性风险熵、脆性综合风险熵定义和研究来分析船舶电力系统的脆性.通过仿真表明以上定义能很好地评价船舶电力系统的脆性.     

基于元胞自动机的舰船电力系统脆性建模

为了分析复杂舰船电力系统的脆性,提高系统的可靠性,以及评价系统的脆性程度,基于熵理论,给出了包括系统自身熵变、脆性联系熵、负熵的舰船电力系统脆性熵的定义,以及舰船电力系统的脆性概率熵、脆性综合概率熵、脆性风险熵和脆性综合风险熵的定义。采用一维扩展元胞自动机(CA),对舰船电力系统脆性熵变和脆性评价分别制定了元胞自动机的演化规则,对舰船电力系统进行脆性建模,以元胞自动机构行序列形式给出了脆性建模的结果,模拟了复杂舰船电力系统在受到内外干扰后脆性演变过程。采用元胞自动机建模可避免求解复杂的微分方程组,使求解简单直观。     

On the parameterization of interleaving and turbulent mixing using CTD data from the Azores Frontal Zone

CTD-data obtained in the Azores Frontal Zone using a towed undulating vehicle are analyzed to study the relationship between characteristics of intrusions and mean parameters of the thermohaline field. A self-similar dependence between intrusion intensity and hydrological parameters is obtained. The most well-founded interpretation of the empirical dependence is as follows: (a) the main source supporting intrusive layering is the salt finger convection; (b) the abrupt decrease of intrusion intensity with the reduction of geostrophic Richardson number obtained from the analysis is explained by the beginning of turbulence when salt fingers do not work any longer, so the “driving force” for intrusive motion disappears. These results are consistent with the conclusions of the paper [Kuzmina N.P., Rodionov V.B., 1992. About the influence of baroclinicity upon generation of the thermohaline intrusions in the oceanic frontal zones. Izvestiya Akad. Nauk SSSR, Atmosperic and Oceanic Physics 28 (10–11), 1077–1086]. These conclusions imply that there are three main mechanisms of intrusive layering at oceanic fronts, namely the 2D baroclinic instability of geostrophic flow, the vertical shear instability and the thermohaline instability where the driving source of intrusive motion is double diffusive convection. The baroclinic and thermohaline instabilities can generate intrusions of large vertical scale, while vertical shear instability usually gives rise to thin turbulent layers. Turbulence in these thin layers can prevent salt finger convection and thus destroy the energy source of the intrusive motion conditioned by thermoclinicity. Therefore, the baroclinicity plays two parts in the processes of the intrusive layering: (1) it prevents double-diffusion interleaving by means of turbulence, and (2) it generates intrusions due to the 2D baroclinic instability of geostrophic current. Using features of thermohaline interleaving as a specific tracer of turbulent mixing, we have estimated turbulent mixing coefficient as k_tRi^−0.8 (Ri>1), where Ri is the geostrophic Richardson number. Application of the proposed approach to other frontal zones is discussed.     

基于熵理论的船舶电力系统脆性分析

将复杂系统脆性引用到船舶电力系统研究中,通过对船舶电力系统脆性概率熵、脆性综合概率熵、脆性风险熵、脆性综合风险熵定义和研究来分析船舶电力系统的脆性。通过仿真表明以上定义能很好地评价船舶电力系统的脆性。     

基于免疫遗传算法的核动力船舶设备故障诊断

针对传统遗传算法收敛速度慢且容易未成熟收敛引起的核动力船舶设备故障诊断响应延迟、误诊、漏诊问题,提出一种基于信息熵的免疫遗传算法用于核动力船舶设备的故障诊断:利用已知船舶设备故障征兆集合,选用概率因果模型,引入信息熵免疫遗传算法,求解具有最大后验概率的故障集合.某船用核动力蒸汽发生器与液压泵故障仿真结果表明,基于信息熵免疫遗传算法优化的概率因果模型不受故障样本的限制,具有较好的通用性,且模型故障诊断精度较高、寻优速度快.本方法同样适用于其他领域的故障诊断问题.     

Structural safety assessment of a pontoon-type VLFS considering damage to the breakwater

 A structural safety assessment of a pontoon-type very large floating structure (VLFS) surrounded by a gravity-type breakwater was carried out for extreme wave conditions by considering the damage to the breakwater. Bending and shear collapses are considered to be a failure mode of the floating structure, while overturning damages the breakwater. The probability of the breakwater overturning, and the transmitted wave height before and after damage to the breakwater, are evaluated using design formulae for port and harbor facilities in Japan. The ultimate bending and shear strengths of the floating structure are calculated by the idealized structural unit method (ISUM) and FEM, respectively. The calculated failure probability for the floating structure is compared with the specified target safety level. It was found that the floating structure under consideration is most likely to fail by bending in transverse waves, and that the corresponding failure probability satisfies the target level. Received: September 12, 2002 / Accepted: October 4, 2002 Acknowledgment. The authors are grateful to Dr. Shigeo Ohmatsu, National Maritime Research Institute, Japan, for allowing us to use the program of hydroelastic response analysis. Address correspondence to: M. Fujikubo (e-mail: fujikubo@naoe.hiroshima-u.ac.jp) Updated from the Japanese original, which won the 2002 SNAJ prize (J Soc Arthit Jpn 2002;190:337–345)     

Using a random displacement model to simulate turbulent particle motion in a baroclinic frontal zone: A new implementation scheme and model performance tests

Numerical artifacts can limit accurate simulation of turbulent particle motion when Lagrangian particle-tracking models are implemented in hydrodynamic models with stratified conditions like fronts. Yet, modeling of individual particle motion in frontal regions is critical for understanding sediment dynamics as well as the transport and retention of planktonic organisms. The objective of this research was to develop a numerical technique to accurately simulate turbulent particle motions in a particle-tracking model embedded within a hydrodynamic model of a frontal zone. A new interpolation scheme, the ‘water column profile’ scheme, was developed and used to implement a random displacement model for turbulent particle motions. A new interpolation scheme was necessary because linear interpolation schemes caused artificial aggregation of particles where abrupt changes in vertical diffusivity occurred. The new ‘water column profile’ scheme was used to fit a continuous function (a tension spline) to a smoothed profile of vertical diffusivities at the x–y particle location. The new implementation scheme was checked for artifacts and compared with a standard random walk model using (1) Well Mixed Condition tests, and (2) dye-release experiments. The Well Mixed Condition tests confirmed that the use of the ‘water column profile’ interpolation scheme for implementing the random displacement model significantly reduced numerical artifacts. In dye-release experiments, high concentrations of Eulerian tracer and Lagrangian particles were released at the same location up-estuary of the salt front and tracked for 4 days. After small differences in initial dispersal rates, tracer and particle distributions remained highly correlated (r = 0.84 to 0.99) when a random displacement model was implemented in the particle-tracking model. In contrast, correlation coefficients were substantially lower (r = 0.07 to 0.58) when a random walk model was implemented. In general, model performance tests indicated that the ‘water column interpolation’ scheme was an effective technique for implementing a random displacement model within a hydrodynamic model, and both could be used to accurately simulate diffusion in a highly baroclinic frontal region. The new implementation scheme has the potential to be a useful tool for investigating the influence of hydrodynamic variability on the transport of sediment particles and planktonic organisms in frontal zones.     

Biologically induced modification of seawater viscosity in the Eastern English Channel during a Phaeocystis globosa spring bloom

To identify the potential relationship between Pheaocystis globosa bloom conditions and seawater properties, a hydrobiological survey was performed in the inshore waters of the Eastern English Channel over the course of the phytoplankton spring bloom. Chlorophyll concentration, auto- and hetero/mixotrophic composition of protists and standing stock, and seawater viscosity were measured weekly from March to June 2004. The decline of the bloom is characterized by a massive foam formation in the turbulent surf zone. Before foam formation, seawater viscosity significantly increased, showing a significant positive correlation with chlorophyll concentration. In contrast, after foam formation this correlation was negative, seawater viscosity kept increasing despite a sharp decrease in chlorophyll concentrations. No significant correlation has been found between seawater viscosity and the composition of the phytoplankton assemblages observed during the survey. However, significant positive correlations have been found between seawater viscosity and both the size and the abundance of P. globosa colonies. From the correlation patterns observed between chlorophyll concentration and seawater viscosity, we suggest that the rheological properties of seawater are mainly driven by extracellular materials associated with colony formation and maintenance rather than by cell composition and standing stock.     

Analysis of velocity fluctuations and their intermittency properties in the surf zone using empirical mode decomposition

We consider here surf zone turbulence measurements, recorded in the Eastern English Channel using a sonic anemometer. In order to characterize the intermittent properties of their fluctuations at many time scales, we analyze the experimental time series using the Empirical Mode Decomposition (EMD) method. The series is decomposed into a sum of modes, each one narrow-banded, and we show that some modes are associated with the energy containing wave breaking scales, and other modes are associated with small-scale intermittent fluctuations. We use the EMD approach in association with a newly developed method based on Hilbert spectral analysis, representing the probability density function in an amplitude–frequency space. We then characterize the fluctuations in a stochastic framework using a cumulant generating function for all scales, and compare the results obtained from direct and classical structure function analysis, to EMD–Hilbert spectral analysis results, showing that the former method saturates at large scales, whereas the latter method is more precise in its scale approach. These results show the strength of the new EMD–Hilbert spectral analysis method for data presenting a strong forcing such as found in shallow water, wave dominated situations.     

Quasi-two-dimensional flow on the polar β-plane: Laboratory experiments

Geophysical turbulence is strongly affected by the variation of the Coriolis parameter with latitude. This variation results in the so-called β-effect, which forces energy from small-scales to be transferred preferentially into zonal motions. This effect results in the formation of narrow jet-like zonal flows that dominate the dynamics and act as transport barriers. Here, laboratory experiments are used to reproduce this effect in decaying turbulent flows. An electromagnetic cell is used to generate an initial field of vorticity in a rotating tank. Under conditions of quasi-geostrophic flow, the β-effect is produced by depth variation of the flow instead of variation of the Coriolis parameter. The effects of changing the container geometry and the overall fluid depth on the production of jets are investigated. The results suggest that this laboratory configuration can be used to model jet formation in the oceans and that increasing fluid depth is a practical way to decrease viscous effects.     

Real-time assessment and prediction on maritime risk state on the Arctic Route

ABSTRACT

Recently, the Arctic Route has become busier with the continuous melting of Arctic ice. However, navigation on the Arctic Route would be much more complex than in normal water as harsh environmental conditions, such as ice-covered water and scarce costal ports that may cause more uncertainty. Nowadays, with the rapid development of sensors on board, more related data has become available. Thus, implementing comprehensive Arctic maritime risk assessment is urgent and necessary in practice. This study proposes an Arctic maritime risk state assessment method including real-time risk state assessment and risk prediction. Specifically, real-time observation samples’ numerical risk state would be firstly obtained with projection pursuit method from 10 risk indicators. Due to the fuzzy uncertainty of single observation set, information diffusion would be applied to provide diffusion estimation on risk probability distribution in order to depict risk state precisely. Also, the accumulated distribution can be regarded as the risk prediction for next time slot and risk entropy is introduced to depict risk tendency directly. Case study based on ‘Yongsheng’ is conducted to demonstrate and verify the effectiveness of the proposed method. The findings can be useful for the operators and management on board during the Arctic voyage.     

Plankton distribution and the impact of copepod grazing on primary production in Fram Strait, Greenland Sea

Two hydrobiological transects across the East Greenland Shelf and the open waters of Fram Strait in summer were chosen to illustrate the distribution and production of phyto- and zooplankton in relation to water masses and ice cover. The parameters used were temperature and salinity, inorganic nutrients, chlorophyll a, primary production, phytoplankton species composition, abundance of the dominant herbivorous copepods Calanus finmarchicus, C. glacialis, C. hyperboreus, Metridia longa and egg production of C. finmarchicus and C. glacialis. Grazing impact of copepodites and adults of these four species was modelled for each station by using egg production rates as an index of growth. Seasonal development of plankton communities was closely associated with the extent of the ice cover, hydrographic conditions and the water masses typical of the different hydrographic domains. Four regions were identified from their biological activities and physical environment: The Northeast Water polynya on the East Greenland Shelf, with a springbloom of diatoms and active reproduction of herbivorous copepods. The pack ice region, dominated by small flagellates and negligible grazing activities. The marginal ice zone, with high variability and strong gradients of autotroph production related to eddies and ice tongues, an active microbial loop and low egg production. The open water, with high station-to-station variability of most of the parameters, probably related to hydrographic mesoscale activities. Here, Phaeocystis pouchetii was a prominent species in the phytoplankton communities. Its presence may at least partly be responsible for the generally low egg production in the open waters. Grazing impact on primary production was always small, due to low zooplankton biomass in the polynya and due to low ingestion in the remaining regions.     

A pseudofield model approach to simulate compartment-fire phenomena for marine fire safety design

In order to prevent the spread of marine compartment fires, it is necessary to understand the governing factors or characteristics of fire-spread phenomena. We present a pseudofield model approach to this problem. We first described a field model of turbulent heat convection based on a standard k – turbulence model. Two-dimensional numerical simulations of a two-linked compartment fire were carried out in order to predict the turbulent convection flow induced by the heat released from the fire. Then a more complicated fire-spread problem of multilinked compartment fires was analyzed by means of a zone model, in which the amounts of oxygen consumption and gas generation were solved by a gas-balance equations system. The effect of threshold conditions on fire propagation and the effect of the thickness of the heat insulation were investigated with numerical simulations.     

Long-term prediction method for the green water load and volume for an assessment of the load line

 To develop a practical prediction method for the green water load and volume on the bow deck in irregular waves, model tests were conducted using a tanker and a cargo ship on a domestic Japanese voyage. The relation between green water load and relative water height at the stem was considered. Based on the finding that the maximum value of the green water load is proportional to the square of the maximum value of the water elevation over the bow top, the probability density functions of the green water load and volume in short-term predictions were proposed. It was verified that the proposed functions show good agreement with the measured distributions, especially in the tails, and were better than conventional functions. Using these functions, long-term predictions of the green water load were carried out. It was confirmed that the present method is more rational than the conventional one for estimating the long-term probability of the green water load. An assessment of the bow height of a domestic Japanese ship from the viewpoint of deck wetness was carried out using these prediction methods. This research was used as the technical background for the revision of domestic rules on load lines, which was enforced in October 2001. Received: July 19, 2002 / Accepted: October 30, 2002 Acknowledgment. Some of the present study was carried out as part of a cooperation project (RR45) with the Shipbuilding Research Association of Japan, supported by the Nippon Foundation. Address correspondence to: Y. Ogawa (e-mail: ogawa@nmri.go.jp)