The influence of different turbulence schemes on modelling primary production in a 1D coupled physical–biological model

A one-dimensional (1D) coupled physical–microbiological model has been applied to a site in the central North Sea. The impact of the choice of the turbulence closure scheme on the modelling the primary production has been investigated.The model was run with four different parameterisations of vertical mixing of heat, momentum and dissolved and suspended matters, using M2 tidal forcing and the hourly mean meteorological forcing of 1989 to reproduce the annual thermal structure and primary production. The four mixing parameterisations are: Level 2 turbulence closure scheme [Mellor, G.L., Yamada, T., 1974. A hierarchy of turbulence closure models for planetary boundary layers. J. Atmos. Sci. 31, 1791–1806; Mellor, G.L., Yamada, T., 1982. Development of a turbulence closure model for geophysical Fluid problems. Rev. Geophys. Space Phys. 20 (4) 851–875] using an explicit numerical scheme [Sharples, J., Tett, P., 1994. Modelling the effect of physical variability on the midwater chlorophyll maximum. J. Mar. Res. 52, 219–238]; a version of the Level 2.5 turbulence closure scheme [Galperin, B., Kantha, L.H., Hassid, S., Rosati, A., 1988. A quasi-equilibrium turbulent energy model for geophysical flows. J. Atmos. Sci. 45, 55–62; Ruddick, K.G., Deleersnijder, E., Luyten, P.J., Ozer, J., 1995. Haline stratification in the rhine/meuse freshwater plume: a 3D model sensitivity analysis. Cont. Shelf Res. 15 (13) 1597–1630] simplified to use an algebraic mixing length by Sharples and Simpson [Sharples, J., Simpson, J.H., 1995. Semidiurnal and longer period stability cycles in the Liverpool Bay region of freshwater influence. Cont. Shelf Res. 15, 295–313], also solved explicitly; the same simplified L2.5 scheme with an implicit numerical solution and modified vertical discretisation scheme [Annan, J.D., 1999. Numerical methods for the solution of the turbulence energy equations in the shelf seas. Int. J. Numer. Methods Fluids 29, 193–206]; and another version of the same scheme (but using a different algebraic mixing length) as described by Xing and Davies [Xing, J., Davies, A.M., 1996a. Application of turbulence energy models to the computation of tidal currents and mixing intensities in the shelf edge regions. J. Phys. Oceanogr. 26, 417–447; Xing, J., Davies, A.M., 1996b. Application of a range of turbulence models to the computation of tidal currents and mixing intensities in shelf edge regions. Cont. Shelf. Res. 16, 517–547; Xing, J., Davies, A.M., 1998. Application of a range of turbulence energy models to the computation of the internal tide. Int. J. Numer. Methods Fluids 26, 1055–1084]. Various model outputs at the sea surface and in depth profiles have been compared with data collected in 1989 as part of the North Sea Project [Huthnance, J.M., 1990. Progress on North Sea Project. NERC News, vol. 12, pp. 25–29, UK]. It is shown that the biological results are extremely sensitive to the small changes in the physical conditions, which arise due to the different turbulence schemes tested. The timing of the spring bloom and the maintenance of the midwater chlorophyll maximum all differ greatly between model runs, and the gross primary production varies by a factor of two from the highest to lowest results. The simplified Level 2.5 scheme, implemented using the numerical methods of Annan [Annan, J.D., 1999. Numerical methods for the solution of the turbulence energy equations in the shelf seas. Int. J. Numer. Methods Fluids 29, 193–206], produces results, which give the best agreement with the available data.     

Late summer stratification, internal waves, and turbulence in the Yellow Sea

Microstructure profiling measurements at two locations in the Yellow Sea (a deeper central basin and a local shelf break) were analyzed focusing on tidal and internal-wave induced turbulence near the bottom and in the pycnocline. A classical three-layer density structure consisting of weakly stratified surface and bottom boundary layers and a narrow sharp pycnocline is developed by the end of warm season. Turbulence in the surface layer was not influenced by the tidal forcing but by the diurnal cycle of buoyancy flux and wind forcing at the sea surface. The enhanced dissipation and diffusivity generated by the shear stress at the seafloor was found in the water interior at heights 10–15 m above the bottom with a phase shift of ~ 5–6 m/h. No internal waves, turbulence, or mixing were detected in the pycnocline in the central basin, in contrast to the pycnocline near the local shelf break wherein internal waves of various frequencies were observed all the time. The thickness of the surface layer near the local shelf break slightly exceeded that of the bottom layer (20 vs. 18 m). A 5–6 m high vertical displacement of the pycnocline, which emerged during the low tide, was arguably caused by the passage of an internal soliton of elevation. During this episode, the gradient Richardson number decreased below 0.25 due to enhanced vertical shear, leading to local generation of turbulence with dissipation rates exceeding the background level by an order of magnitude.     

A numerical model of the long term flow along the Malin-Hebrides shelf

A three dimensional hydrodynamic model of the Malin-Hebrides shelf region is used to investigate the spatial variability of the wind and tidally induced residual flow in the region and the influence of flow from the Irish Sea and along the shelf edge. By this means it is possible to understand the spatial variability in the long term observed flow fields in the region and the range of driving forces producing this flow. The model uses a sigma coordinate grid in the vertical with a finer grid in the near surface and near bed shear layers. The vertical diffusion of momentum in the model is parameterised using an eddy viscosity coefficient which is derived from turbulence energy closure models. Two different turbulence models are used to compute the eddy viscosity, namely a two-equation (itq²−q²ℓ) model which has prognostic equations for both turbulence energy and mixing length and a simpler model in which the mixing length is a specified algebraic function of the water depth.The wind induced response to spatially and temporally constant orthogonal wind stresses, namely westerly and southerly winds of 1 N m⁻², are derived from the model. By using orthogonal winds and assuming linearity, then to first order the response to any wind direction can be derived. Computed flows show a uniform wind driven surface layer of magnitude about 3% of the wind speed and direction 15 ° to the right of the wind, in deep water. Currents at depth particularly in the shelf edge and near coastal region show significant spatial variability which is related to variations in bottom topography and the coastline.Calculations show that tidal residual flows are only significant in the near coastal regions where the tidal current is strong and exhibits spatial variability. Flow into the region from the Irish Sea through the North Channel although having its greatest influence in the near coastal region, does affect currents near the shelf edge region. Again the spatial variability of the flow is influenced by topographic effects.A detailed examination of wind induced current profiles together with turbulence, mixing length and viscosity, at a number of locations in the model from deep ocean to shallow near coastal, shows that both turbulence models yield comparable results, with the mixing length in the two equation model showing a similar dependence to that specified in the simpler turbulence model.Calculations clearly show that flow along the shelf edge area to the west of Ireland and from the Irish Sea entering the region, together with local wind forcing can have a major effect upon currents along the Malin-Hebrides shelf. The flow fields show significant spatial variability in the region, comparable to those deduced from long term tracer measurements. The spatial variability found in the calculations suggests that a very intense measurement programme together with inflow measurements into the area is required to understand the circulation in the region, and provide data sets suitable for a rigorous model validation.     

Large eddy simulations for quasi-2D turbulence in shallow flows: A comparison between different subgrid scale models

In this study, the performance of the horizontal large eddy simulation module, developed at the University of Leuven (HLES-KULeuven module) is assessed. A comparison between different subgrid scale models has been carried out. The study is concerned with the non-rotating and unstratified flows. The results of the simulation for an oscillatory backward facing (BFS) flow are presented in case of an expanding flume based on a one-length scale approach and a two-length scale approach. Three subgrid scale (SGS) models have been tested: Smagorinsky SGS model (Smagorinsky, J., (1963). General circulation experiments with the primitive equations, I. the basic experiments. Monthly Weather Review, 91(3), 99–164), Uittenbogaard SGS model (Uittenbogaard, R.E., and van Vossen, B., (2004). Subgrid-scale model for quasi-2D turbulence in shallow water. Shallow Flows. Jirka and Uijttewaal (Eds.), Taylor & Francis Group, London, ISBN 90 5809 700 5) and a proposed two-length scale approach. The first two models are considered to be a one-length scale models. A simulation without a subgrid scale model for the horizontal mixing has also been conducted. In all simulations, a quadratic friction model parameterizes the dissipation produced by the 3D-subdepth scale turbulence. The two-length scale concept uses a newly mixing length formulation for the quasi-2D turbulence and doesn't depend on the filter width in contrast to the one-length scale approach, in which the mixing length is function of the filter width. The outputs of the HLES-KULeuven module have been compared with the experimental data taken from Stelling, G.S., and Wang, L.X., (1984). Experiments and computations on separating flow in an expanding flume. Dept. Civil Engineering, Delft University of Technology, Report 2–84.). The two-length scale approach has been validated with experimental data from SERC Flood Channel Facility at HR Wallingford. In general, there is a qualitative agreement with the experimental data. It has also been found that the two-length scale approach produces more elongated and less isotropic vortex than the one-length scale models.     

Do observations adequately resolve the natural variability of oceanic turbulence?

Especially in high Reynolds number, naturally-occurring flows, turbulence is a highly variable process. It is challenging to measure yet it is vital that we do so in order to quantify the internal transports of mass, nutrients, energy and momentum. Isolated turbulence profiles are difficult to interpret; systematic sampling and subsequent averaging are necessary. Confidence in our ability to properly sample turbulence arises from intergroup comparisons, comparisons with other methods to assess mixing coefficients and, most fundamentally, the constraints imposed by the governing fluid dynamics on both energy losses via viscous dissipation caused by turbulence and on the mixing that results from turbulence. Several examples in which fluid processes have been isolated from the full range of oceanic motions are reviewed in this light. These examples show how observationally-derived estimates of turbulence dissipation or mixing are consistent with larger scale constraints. The larger oceanographic problem of defining the full geographic variability of mixing remains.     

Improving the physics of a coupled physical–biogeochemical model of the North Atlantic through data assimilation: Impact on the ecosystem

Several studies on coupled physical–biogeochemical models have shown that major deficiencies in the biogeochemical fields arise from the deficiencies in the physical flow fields. This paper examines the improvement of the physics through data assimilation, and the subsequent impact on the ecosystem response in a coupled model of the North Atlantic. Sea surface temperature and sea surface height data are assimilated with a sequential method based on the SEEK filter adapted to the coupling needs. The model domain covers the Atlantic from 20°S to 70°N at eddy-permitting resolution. The biogeochemical model is a NPZD-DOM model based on the P3ZD formulation. The results of an annual assimilated simulation are compared with an annual free simulation.With assimilation, the representation of the mixed layer depth is significantly improved in mid latitudes, even though the mixed layer depth is generally overestimated compared to the observations. The representation of the mean and variance of the currents is also significantly improved.The nutrient input in the euphotic zone is used to assess the data assimilation impact on the ecosystem. Data assimilation results in a 50% reduction of the input due to vertical mixing in mid-latitudes, and in a four- to six-fold increase of the advective fluxes in mid-latitudes and subtropics. Averaged zonally, the net impact is a threefold increase for the subtropical gyre, and a moderate (20–30%) decrease at mid and high latitudes.Surface chlorophyll concentration increases along the subtropical gyre borders, but little changes are detected at mid and high latitudes. An increase of the primary production appears along the Gulf Stream path, but it represents only 12% on average for mid and high latitudes. In the subtropical gyre centre, primary production is augmented but stays underestimated (20% of observations). These experiments show the benefits of physical data assimilation in coupled physical–biogeochemical applications.     

Simulations of biological production in the Rhodes and Ionian basins of the eastern Mediterranean

The biological production characteristics of the Rhodes and western Ionian basins of the eastern Mediterranean are studied by a one-dimensional, coupled physical–biological model. The biological model involves single aggregated compartments of phytoplankton, zooplankton, detritus as well as ammonium and nitrate forms of the inorganic nitrogen. It interacts with the physical model through the vertical eddy diffusivity which is calculated using the Mellor–Yamada level 2.5 turbulence parameterization. The model simulations demonstrate the importance of the contrasting physical oceanographic characteristics of these two basins on affecting their yearly planktonic structures. The annual primary production in the Rhodes basin is estimated as 97 g C m² yr⁻¹ which is comparable with the northwestern Mediterranean. The western Ionian basin, on the contrary, possesses only 10% of the Rhodes' productivity and therefore represent a most oligotrophic site in the eastern Mediterranean. The Rhodes basin reveals a strong bloom in early spring, typically in March, a weaker bloom in early winter, typically in January, and a subsurface production below the seasonal thermocline during summer. This structure is slightly modified in the western Ionian basin, and the early winter and early spring blooms are merged to cover the entire winter. These results are supported favorably by the available observations both in their magnitudes and timing.     

Description of a flexible and extendable physical–biogeochemical model system for the water column

A modelling system for coupled physical–biogeochemical simulations in the water column is presented here. The physical model component allows for a number of different statistical turbulence closure schemes, ranging from simple algebraic closures to two-equation turbulence models with algebraic second-moment closures. The biogeochemical module consists of models which are based on a number of state variables represented by their ensemble averaged concentrations. Specific biogeochemical models may range from simple NPZ (nutrient–phytoplankton–zooplankton) to complex ecosystem models. Recently developed modified Patankar solvers for ordinary differential equations allow for stable discretisations of the production and destruction terms guaranteeing conservative and non-negative solutions. The increased stability of these new solvers over explicit solvers is demonstrated for a plankton spring bloom simulation. The model system is applied to marine ecosystem dynamics the Northern North Sea and the Central Gotland Sea. Two different biogeochemical models are applied, a conservative nitrogen-based model to the North Sea, and a more complex model including an oxygen equation to the Baltic Sea, allowing for the reproduction of chemical processes under anoxic conditions. For both applications, earlier model results obtained with slightly different model setups could be basically reproduced. It became however clear that the choice for ecosystem model parameters such as maximum phytoplankton growth rates does strongly depend on the physical model parameters (such as turbulence closure models or external forcing).     

Parameterization of air–sea gas fluxes at extreme wind speeds

Air–sea flux measurements of O₂ and N₂ obtained during Hurricane Frances in September 2004 [D'Asaro, E. A. and McNeil, C. L., 2006. Measurements of air–sea gas exchange at extreme wind speeds. Journal Marine Systems, this edition.] using air-deployed neutrally buoyant floats reveal the first evidence of a new regime of air–sea gas transfer occurring at wind speeds in excess of 35 m s^− 1. In this regime, plumes of bubbles 1 mm and smaller in size are transported down from near the surface of the ocean to greater depths by vertical turbulent currents with speeds up to 20−30 cm s^− 1. These bubble plumes mostly dissolve before reaching a depth of approximately 20 m as a result of hydrostatic compression. Injection of air into the ocean by this mechanism results in the invasion of gases in proportion to their tropospheric molar gas ratios, and further supersaturation of less soluble gases. A new formulation for air–sea fluxes of weakly soluble gases as a function of wind speed is proposed to extend existing formulations [Woolf, D.K, 1997. Bubbles and their role in gas exchange. In: Liss, P.S., and Duce, R.A., (Eds.), The Sea Surface and Global Change. Cambridge University Press, Cambridge, UK, pp. 173–205.] to span the entire natural range of wind speeds over the open ocean, which includes hurricanes. The new formulation has separate contributions to air–sea gas flux from: 1) non-supersaturating near-surface equilibration processes, which include direct transfer associated with the air–sea interface and ventilation associated with surface wave breaking; 2) partial dissolution of bubbles smaller than 1 mm that mix into the ocean via turbulence; and 3) complete dissolution of bubbles of up to 1 mm in size via subduction of bubble plumes. The model can be simplified by combining “surface equilibration” terms that allow exchange of gases into and out of the ocean, and “gas injection” terms that only allow gas to enter the ocean. The model was tested against the Hurricane Frances data set. Although all the model parameters cannot be determined uniquely, some features are clear. The fluxes due to the surface equilibration terms, estimated both from data and from model inversions, increase rapidly at high wind speed but are still far below those predicted using the cubic parameterization of Wanninkhof and McGillis [Wannikhof, R. and McGillis, W.R., 1999. A cubic relationship between air–sea CO₂ exchange and wind speed. Geophysical Research Letters, 26:1889–1892.] at high wind speed. The fluxes due to gas injection terms increase with wind speed even more rapidly, causing bubble injection to dominate at the highest wind speeds.     

一种河道紊流三维数值模拟离散格式的应用探讨

通过引用气体动力学中的Osher格式计算控制体交界面通量,建立了一种基于有限体积法的河道紊流三维数学模型。该格式是对准确的非线性黎曼问题的近似解法,具有很好的守恒性、逆风性、单调性和高分辨率捕捉间断的能力,解决了河床变化剧烈处水流计算易发散的难题,在紊流模型的选取上,考虑到紊动能量和紊动尺度在水平方向和垂直方向上的差异,平面上采用大涡模型,垂向上采用k-ε模型。通过桥墩和丁坝局部水流实测资料校核模型模拟局部复杂水流结构的能力,模型计算得到的马蹄涡、回流结构及流速分布与实测基本一致。     

Krill behaviour as recorded by acoustic doppler current profilers in the Gullmarsfjord

Vertical distribution of sound scattering layers were observed using bottom deployed acoustic doppler current profilers (ADCP) during early spring of 1996 and autumn of 1997 in the Gullmarsfjord on the Swedish west coast. Variations in relative backscatter were interpreted in relation to horizontal water velocities, oxygen saturation as well as differences in the light, salinity and temperature regimes. Net catches revealed that much of the backscatter below 20-m depth was associated with the presence of krill, principally Meganyctiphanes norvegica.Horizontal currents seemed to influence the migration and distribution of krill, which showed weak vertical migration patterns with low abundance during periods of strong intermediate in- and outflows, while during periods with weaker currents, a more regular diel migration occurred. Horizontal water velocities >5 cm s⁻¹ seemed to have the potential to decrease the peak in the backscatter profile. Mean vertical migration rates of krill was 1 cm s⁻¹, while maximum vertical migration rates were estimated to be 2.5–3 cm s⁻¹. The range of the vertical migration was different in 1997 due to severe oxygen deficiency in the bottom water, which prevented the krill from descending >80 m. The commencement of vertical migration correlated closely to the seasonal light conditions. The descent was immediately triggered by sunrise, while ascent occurred with a delay of about 1 h at sunset.     

Temporal variations of water flow between the Venetian lagoon and the open sea

Long-term measurements of the water flow at three Venice Lagoon inlets with the bottom-mounted ADCPs show that the main part of the variance (>90%) is associated with the tidal variability. Semi-diurnal constituents (mainly M2 and S2) are responsible for about 80% of the flow variance. Phase-lag between the axial current and sea-level is on the order of 2 h for M2 and 4 h for the K1, the maximum inflow leading the sea-level maximum. Phase-difference of tidal flows between inlets shows that Chioggia leads both Malamocco and Lido, suggesting that tidal signal progresses northward, thus in the opposite direction of both the semi-diurnal and diurnal tidal signals in the open Adriatic currents. The sea-level slope between the open sea and the lagoon interior controls the inlet flow, which is due to the time lag being constant for all tidal constituents. It was shown that tidal oscillations at Punta Salute (lagoon interior) lag those in Lido by about 45 min. The pressure gradient due to the sea-level slope generates the flow acceleration. Only for large current speeds (>0.5 m/s), the bottom friction term becomes equally important as the local acceleration and the horizontal pressure gradient terms. Wind effects manifest as a remote forcing through Adriatic seiches at semi-diurnal and diurnal scales, and as a local forcing at very long time scales on the order of a month. This latter mechanism is limited to a winter period (November–January). Seiches are present over the entire year, being however, more energetic and frequent during autumn and winter.     

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.     

Study of the seasonal cycle of the biogeochemical processes in the Ligurian Sea using a 1D interdisciplinary model

A one-dimensional coupled physical–biogeochemical model has been built to study the pelagic food web of the Ligurian Sea (NW Mediterranean Sea). The physical model is the turbulent closure model (version 1D) developed at the GeoHydrodynamics and Environmental Laboratory (GHER) of the University of Liège. The ecosystem model contains 19 state variables describing the carbon and nitrogen cycles of the pelagic food web. Phytoplankton and zooplankton are both divided in three size-based compartments and the model includes an explicit representation of the microbial loop including bacteria, dissolved organic matter, nano-, and microzooplankton. The internal carbon/nitrogen ratio is assumed variable for phytoplankton and detritus, and constant for zooplankton and bacteria. Silicate is considered as a potential limiting nutrient of phytoplankton's growth. The aggregation model described by Kriest and Evans in (Proc. Ind. Acad. Sci., Earth Planet. Sci. 109 (4) (2000) 453) is used to evaluate the sinking rate of particulate detritus. The model is forced at the air–sea interface by meteorological data coming from the “Côte d'Azur” Meteorological Buoy. The dynamics of atmospheric fluxes in the Mediterranean Sea (DYFAMED) time-series data obtained during the year 2000 are used to calibrate and validate the biological model. The comparison of model results within in situ DYFAMED data shows that although some processes are not represented by the model, such as horizontal and vertical advections, model results are overall in agreement with observations and differences observed can be explained with environmental conditions.     

Mixing in the surface boundary layer of a tropical freshwater reservoir

A combined observational-modeling study was conducted to investigate turbulence mixing, and the relation to surface forcing, in the surface boundary layer (SBL) of a tropical, high-altitude, freshwater reservoir. A suite of vertical profiles of temperature microstructure, collected at three different stations of one-day duration each, provided estimates of dissipation rates of turbulence kinetic energy, , and temperature variance, χ. Numerical simulations of and χ, using state-of-the-art, public domain, two-equation turbulence closure models, compared favorably with the observations and reproduced the dynamics of daytime wind mixing as well as the vertical and temporal turbulence structure during nighttime convective conditions.Two independent estimates of vertical eddy diffusivities in the stably stratified (daytime) SBL, computed from the microstructure measurements, agreed closely, and the near surface heat and buoyancy fluxes, computed from the diffusivities, were similar to those computed independently from surface meteorology. Model generated eddy diffusivities agreed closely with the observed values, except those generated by K profile parameterization (KPP) model simulations. The good agreement provides confidence that nutrient fluxes in the SBL may be accurately computed from the models when forced with regularly measured surface meteorological parameters. The consequences are important for estimation of daily primary productivity rates in the euphotic zone and the ability to predict algal blooms such as those observed in the present reservoir.     

One-dimensional modelling of the plankton ecosystem of the north-western Corsican coastal area in relation to meteorological constraints

In order to study the influence of wind mixing on the spring variability of the plankton production of the north western Corsican coastal area, a one-dimensional (1D), vertical, coupled hydrodynamic/biological model (ECOHYDROMV) is used. A hydrodynamic 1D model of the water column with a k–l turbulent closure is applied. The biological model comprises six state variables, representing the plankton ecosystem in the spring period: phytoplankton, copepods, nitrate, ammonium, particulate organic matter of phytoplanktonic origin and particulate organic matter of zooplanktonic origin. The system is influenced by turbulence (expressed by the vertical eddy diffusivity), temperature and irradiance. The model takes into account momentum and heat surface fluxes computed from meteorological data in order to simulate a typical spring atmospheric forcing for the considered area. Results show that primary production vertical structure is characterised by a subsurface maximum which deepens with time and is regulated by the opposite gradients of nitrate concentration and irradiance. Surface plankton productivity is mainly controlled by turbulent vertical transport of nutrients into the mixed layer. The short time scale variability of turbulent mixing generated by the wind appears to be responsible for the plurimodal shape of plankton blooms, observed in the considered area. Furthermore, the model is applied to the study of the spring evolution of the plankton communities off the bay of Calvi (Corsica) for the years 1986 and 1988. In order to initiate and validate the model, time series of hydrological, chemical and biological data have been used. The model reproduces accurately the spring evolution of the phytoplankton biomass measured in situ and illustrates that its strong variability in those years was in close relation to the variability of the wind intensity.     

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.     

复杂大气环境下的飞行研究

飞机在海面上空飞行时比在陆地上空更容易受到大气环境的影响,文中针对海洋上空的大气环境,建立了风和大气紊流的工程化模型,研究了飞机在变化风场中飞行时,受到风水平加速度和垂直加速度影响后的响应情况,同时分析了大气紊流对飞行的影响,提出了一些抑制阵风扰动的控制方法.     

Steady-state single cell model simulations of photoacclimation in a vertically mixed layer: implications for biological tracer studies and primary productivity

A numerical single-cell photoacclimation–diffusion model was constructed and used to develop criteria regarding the use of individual phytoplankton cells as tracers for vertical mixing and to illustrate how rates of vertical mixing might affect phytoplankton physiology. Both first-order and logistic representations of photoacclimation kinetics were used. Steady state was assumed for simplicity and to provide a starting point for further investigations. The modeled variance and higher moments (within a phytoplankton population) of a generic photoacclimative parameter all show trends, which are diagnostic of mixing rates and/or boundary effects. This allowed the establishment of criteria by which frequency distributions of phytoplankton physiological properties (e.g., cell fluorescence) might be used as indicators of vertical mixing. The same model can be used to predict the effects of vertical mixing on phytoplankton productivity and growth. Application of the model to both photosynthesis and carbon to chlorophyll ratios suggested that a combination of vertical mixing and hysteresis (as represented in the logistic model of photoacclimation) in acclimation kinetics can enhance specific growth rates of phytoplankton. This enhanced growth occurred as a result of mixing-induced variation in carbon to chlorophyll ratios and is in contrast to chlorophyll-specific productivity, which was maximal at low mixing rates. Differential rates of photoacclimation to upward vs. downward shifts in irradiance, may enable phytoplankton cells to better survive in a turbulent environment.