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.     

The Gulf Stream, rings and North Atlantic eddy structures from remote sensing (Altimeter and SeaWiFS)

Seasonal SeaWiFS chlorophyll a concentrations cycles and annual changes of altimeter Sea Level Anomaly are derived for the subtropical North Atlantic near  35°N and along a Gulf Stream axis. Spatial structure of SeaWiFS, is defined in terms of deviations from a local seasonal cycle and examined in relation to altimeter eddy structure. In the subtropical region near 35°N, SeaWiFS structure is evident during the spring bloom period with a scale of  430 km, or about twice the eddy scale. A Gulf Stream axis has been selected as a region where the Sea Level Anomaly variance is a maximum. Eddy propagation speeds and scales are examined. Cold-core (cyclonic) rings correspond to areas of high SeaWiFS chlorophyll a. Warm-core (anticyclonic) rings relate to areas of low chlorophyll concentration. Both SeaWiFS structure and eddy structure have a spatial scale of  450 km or twice the ring scale along the Gulf Stream axis. SeaWiFS chlorophyll anomalies and Altimeter Sea Level Anomaly structure have an overall negative correlation coefficient of r = − 0.34. Swirl currents between eddies redistribute surface chlorophyll concentrations and can spatially bias maximum and minimum concentration levels off eddy centre.     

A trace metal (Pb, Cd, Zn, Cu) balance for surface waters in the eastern Gotland Basin, Baltic Sea

Trace metal inventories and transboundary fluxes were estimated for the surface waters of the eastern Gotland Basin, and a mass balance for this water body was derived.The study area was delineated by a box which bordered vertically by the sea surface and the halocline and horizontally by the coast of Gotland and the coast of Latvia/Lithuania. For the calculation of the trace metal inventories in the box, a high-resolution sampling was carried out during different seasons. Additional sediment trap studies were performed to calculate the vertical particulate trace metal fluxes through the halocline.The following transboundary fluxes have been estimated: (a) atmospheric input of trace metals; (b) lateral transport from adjacent bodies of surface water; (c) diapycnal fluxes of dissolved trace metals at the halocline; (d) sinking of trace metals, associated with particles, through the halocline, (e) vertical advection through the halocline.The lateral transport into and out of the surface box is important for the metals (Cd, Cu, Zn) with “nutrient-like” behaviour. For particle reactive elements like lead, vertical sedimentation and lateral transport out of the box as much as atmospheric input and lateral input into the box are in the same order of magnitude. Turbulent diapycnal mixing plays a minor role compared to the advective fluxes.The estimated range for the residence times for the region under investigation are  0.5–1.3 years for Pb;  2–11 years for Zn,  7–36 years for Cd and 42–89 years for Cu. This demonstrates that the system reacts very fast for particle reactive elements like Pb, while for copper sedimentation processes are not the preferential sink and can be neglected.     

Air–sea gas exchange at extreme wind speeds measured by autonomous oceanographic floats

Measurements of the air–sea fluxes of N₂ and O₂ were made in winds of 15–57 m s^− 1 beneath Hurricane Frances using two types of air-deployed neutrally buoyant and profiling underwater floats. Two “Lagrangian floats” measured O₂ and total gas tension (GT) in pre-storm and post-storm profiles and in the actively turbulent mixed layer during the storm. A single “EM-APEX float” profiled continuously from 30 to 200 m before, during and after the storm. All floats measured temperature and salinity. N₂ concentrations were computed from GT and O₂ after correcting for instrumental effects. Gas fluxes were computed by three methods. First, a one-dimensional mixed layer budget diagnosed the changes in mixed layer concentrations given the pre-storm profile and a time varying mixed layer depth. This model was calibrated using temperature and salinity data. The difference between the predicted mixed layer concentrations of O₂ and N₂ and those measured was attributed to air–sea gas fluxes F_BO and F_BN. Second, the covariance flux F_CO(z) = wO₂′(z) was computed, where w is the vertical motion of the water-following Lagrangian floats, O₂′ is a high-pass filtered O₂ concentration and (z) is an average over covariance pairs as a function of depth. The profile F_CO(z) was extrapolated to the surface to yield the surface O₂ flux F_CO(0). Third, a deficit of O₂ was found in the upper few meters of the ocean at the height of the storm. A flux F_SO, moving O₂ out of the ocean, was calculated by dividing this deficit by the residence time of the water in this layer, inferred from the Lagrangian floats. The three methods gave generally consistent results. At the highest winds, gas transfer is dominated by bubbles created by surface wave breaking, injected into the ocean by large-scale turbulent eddies and dissolving near 10-m depth. This conclusion is supported by observations of fluxes into the ocean despite its supersaturation; by the molar flux ratio F_BO/F_BN, which is closer to that of air rather than that appropriate for Schmidt number scaling; by O₂ increases at about 10-m depth along the water trajectories accompanied by a reduction in void fraction as measured by conductivity; and from the profile of F_CO(z), which peaks near 10 m instead of at the surface.At the highest winds O₂ and N₂ are injected into the ocean by bubbles dissolving at depth. This, plus entrainment of gas-rich water from below, supersaturates the mixed layer causing gas to flux out of the near-surface ocean. A net influx of gas results from the balance of these two competing processes. At lower speeds, the total gas fluxes, F_BO, F_BN and F_CO(0), are out of the ocean and downgradient.     

Estimation of TKE dissipation rates in dense bottom plumes using a Pulse Coherent Acoustic Doppler Profiler (PC-ADP) — Structure function approach

During a hydrographic survey in January 2006 the spreading of inflowing saline water was observed in the Arkona Basin (Western Baltic Sea). Two bottom mounted ‘pulse coherent’ acoustic Doppler profilers (PC-ADP) were used to measure the near-bottom current field of the dense plume with a high temporal (1 s) and spatial resolution (5 cm). In order to estimate the dissipation rate of turbulent kinetic energy () a structure function approach was applied to the beam velocity data. Simultaneous measurements with a microstructure shear profiler (MSS) and an acoustic Doppler velocimeter (ADV) supplied independent data for the verification of the structure function method. Additional measurements with standard CTD, near-bottom towed and vessel mounted acoustic Doppler current profilers (ADCP) completed the data set.The estimated dissipation rates from the structure function approach fit well with the values derived from the ADV and the MSS probe. It is shown that the structure function approach is a reliable and easily applicable method to derive estimates of TKE dissipation rates from PC-ADP beam velocities. The observed dissipation rates ranged between 5 · 10^− 6 and 1 · 10^− 8 W kg^− 1 depending on the hydrographic conditions. Inside the plume the dissipation rates exceeded that of the overlaying brackish water by two orders of magnitude. Since the noise level of velocity data in pulse coherent mode is considerably lower than in the Doppler mode the PC-ADP can also be used for estimates in marine environments with low turbulence level. Reynolds stresses estimated from the PC-ADP and the ADV agreed well at the same depth level. TKE production derived from PC-ADP measurements compared reasonably well with the dissipation rate of TKE in a varying environment.     

Relationship between phytoplankton bloom and wind stress in the sub-polar frontal area of the Japan/East Sea

The seasonal dynamics of phytoplankton blooms in the central Japan/East Sea (JES) show pronounced year-to-year variability based on Sea-viewing Wide Field-of-view Sensor (SeaWiFS; 19972003) and Moderate Resolution Imaging Spectroradiometer (MODIS)/Terra (20002003) observations. Wind seems to strongly influence this variability. To study the relationship between wind and bloom initiation, we analyzed daily, remotely sensed wind stress data (Active Microwave Instrument–wind [AMI–wind], NASA Scatterometer [NSCAT], and Quick Scatterometer [QuickSCAT]: 19972003) and daily chlorophyll concentrations based on ocean color data (SeaWiFS and MODIS). The results agreed well with the hypotheses; in spring, blooms began 615 days after wind stress weakened. Fall blooms started 39 days after wind strengthened. We also simulated seasonal changes using a simple light–nutrient model using two values for the respiration ratio: 10% and 20%. The use of 20% seemed to reproduce the timing of the spring bloom quite well but underestimated the absolute level of chlorophyll concentration. On the other hand, using 10% produced a better estimation of the chlorophyll concentration but failed to match the timing. Neither of the model runs reproduced the timing of the fall bloom well.     

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.     

Dissolved organic matter in shelf waters off the Ría de Vigo (NW Iberian upwelling system)

Net in situ production and export of dissolved organic carbon (DOC) and nitrogen (DON) have been studied in shelf waters off the Ría de Vigo (NW Spain), as part of a comprehensive hydrographic survey carried out from September 1994 to September 1995 with a fortnight periodicity. DOC and DON correlated well (r=+0.78), the slope of the regression line being 12.0±0.7 mol-C mol-N⁻¹, about twice the Redfieldian slope of particulate organic matter, 6.5±0.2 mol-C mol-N⁻¹ (r=+0.95). Labile DOC and DON accumulated in the upper 50 m during the upwelling season (March–September), mainly after prolonged periods of wind relaxation, when horizontal flows were reduced. This labile material represented 50% and 35% of the total (dissolved+particulate) organic carbon and nitrogen susceptible of microbial utilisation, which assert the key contribution of dissolved organic matter (DOM) to the export of new primary production in the NW Iberian upwelling system. This surface excess in shelf waters appeared to be formed into the highly productive Ría de Vigo (a large coastal indentation) at net rates of 4.4 μM-C d⁻¹ and 1.3 μM-C d⁻¹ in the inner and outer segments of the embayment respectively, and subsequently exported to the shelf. Once in the shelf, simple dilution with the inert DOM pool of recently upwelled Eastern North Atlantic Central Water (ENACW) occurred. Eventually, the DOM excess produced during the upwelling season is exported to the adjacent open ocean waters by the coastal circulation. Conversely, during the unproductive downwelling season (October–February), the lowest DOC and DON levels were recorded and export was prevented by the characteristic downwelling front associated to the seasonal poleward slope current.     

The interaction of tidal advection, diffusion and mussel filtration in a tidal channel

Time series measurements of flow and pigment concentrations (Chl) in the Menai Strait have revealed that the strong residual flow in a tidal channel ( 500 m³ s^− 1) transports phytoplankton from the open sea into the channel where much of it is consumed by suspension feeders, mainly in commercial beds of Mytilus edulis. The progressive depletion of phytoplankton along the channel results in a strong horizontal gradient of plankton and hence Chl. Tidal displacement of this gradient causes large (± 50% of mean) oscillations of Chl in the vicinity of the mussel beds. Vertical mixing by the strong tidal flows is sufficiently vigorous for most of the tidal cycle to ensure that downward diffusion can resupply the near-bed layer although there are indications of some transient depletion around slack water.This paradigm of the interaction of advection, diffusion and filtration determining the distribution of plankton and its supply to mussels has been encapsulated in a series of simple models forced only by boundary values. In the first, a 1-D model of tidal flow in the channel reproduces the principal features of the observed currents including the unusually large spatial change in phase of the currents and the variation of the residual transport with tidal range. The flow field from this physical model is used to drive a second model based on the advection diffusion equation for Chl with a source at the Irish Sea boundary and a sink over the mussel bed. This model illustrates the formation of a strong Chl gradient along the channel and simulates the amplitude and phase of the M₂ oscillations of Chl and the development of the M₄ variation apparent in the observations. This second model has been extended to 2-D over the mussel beds to allow investigation of the effects of water column mixing. The model indicates that only for a short period ( 30 min), close to slack water, is mixing sufficiently reduced to permit the development of a depletion boundary layer and then only within  1 m from the bottom, a result which is consistent with the observations.     

Assimilation of sea surface temperature in the POL Coastal Ocean Modelling System

A one-dimensional scheme is used to assimilate satellite Sea Surface Temperature data into the Proudman Oceanographic Laboratory Coastal Ocean Modelling System, set up in the Irish Sea with a fine resolution ( 1.8 km). The capabilities of the assimilation scheme are investigated using two different sets of satellite data, of lower and similar resolution to that of the model respectively. Comparison of results with independent data show that assimilation improves the modelled Sea Surface Temperature, but does not address model representation of the temperature vertical structure. It is concluded that for the Irish Sea and at the scales resolved by the model, the assimilation problem cannot be approached in a one-dimensional framework. It is also pointed out that forecast error needs to account explicitly for errors in the representation of the vertical structure of the thermal field.Three-dimensional methods that are suited for coastal systems are then suggested.     

A dense water outflow from the Ross Sea, Antarctica: Mixing and the contribution of tides

We use hydrographic, current, and microstructure measurements, and tide-forced ocean models, to estimate benthic and interfacial mixing impacting the evolution of a bottom-trapped outflow of dense shelf water from the Drygalski Trough in the northwestern Ross Sea. During summer 2003 an energetic outflow was observed from the outer shelf ( 500 m isobath) to the  1600 m isobath on the continental slope. Outflow thickness was as great as  200 m, and mean speeds were  0.6 m s^− 1 relative to background currents exceeding  1 m s^− 1 that were primarily tidal in origin. No outflow was detected on the slope in winter 2004, although a thin layer of dense shelf water was present on the outer shelf. When the outflow was well-developed, the estimated benthic stress was of order one Pascal and the bulk Froude number over the upper slope exceeded one. Diapycnal scalar diffusivity (K_z) values in the transition region at the top of the outflow, estimated from Thorpe-scale analysis of potential density and measurements of microscale temperature gradient from sensors attached to the CTD rosette, were of order 10^− 3−10^− 2 m² s^− 1. For two cases where the upper outflow boundary was particularly sharply defined, entrainment rate w_e was estimated from K_z and bulk outflow parameters to be  10^− 3 m s^− 1 ( 100 m day^− 1). A tide-forced, three-dimensional primitive equation ocean model with Mellor-Yamada level 2.5 turbulence closure scheme for diapycnal mixing yields results consistent with a significant tidal role in mixing associated with benthic stress and shear within the stratified ocean interior.     

A novel method for estimating vertical eddy diffusivities using diurnal signals with application to western Long Island Sound

We present an approach that allows the estimation of vertical eddy diffusivity coefficients from buoy measurements made at two or more depths. By measuring the attenuation and phase lag of a scalar signal generated periodically at the surface as it propagates downwards, the vertical eddy diffusivity coefficients can be calculated as K_v = ωΔz²/2ln²(α₂/α₁), where α₂/α₁ is the ratio of the real amplitudes at frequency ω at the two depths separated by Δz = z₂ − z₁; as K_V = ωΔz²/2φ², where φ is the phase lag at the frequency ω; or as K_v = iωΔz²/ln²(U₂/U₁), where U₂/U₁ is the ratio of the complex signal amplitudes at the two depths. The method requires that horizontal fluxes be small at the ω frequency and that the signal-to-noise ratios at the two depths allow the determination of the amplitude and phase of ω.Application of this method to summertime 2004 western Long Island Sound oxygen and temperature buoy measurements at two depths provides a time-series of two-day average vertical eddy diffusivity estimates. Using these eddy diffusivities in conjunction with measured vertical concentration gradients, we obtain a time-series of vertical transport rates for oxygen and heat and estimate mean downward fluxes for June and July as 150–260 mMol m^− 2 day^− 1 and 100–400 W m^− 2 respectively. These estimates are of a similar magnitude to sub-pycnocline O₂ and heat demands of 240 ± 200 mMol m^− 2 day^− 1 and 180 ± 60 W m^− 2 that we infer from simple budgets, implying that vertical transport is significant to both budgets.The eddy coefficients obtained from the independent O₂ and temperature measurements have a 68% correlation, and the O₂ flux estimates show a correlation of 41% to measured rates of change in bottom dissolved oxygen levels. Our results indicate that extended time-series of eddy diffusivity coefficients can be obtained from in situ buoy measurements and the method shows promise as a way to constrain the vertical transport variability in budgets of dissolved materials in estuaries.     

Microscale variability in the distributions of large fluorescent particles observed in situ with a planar laser imaging fluorometer

It has been known for decades that particle-size and biomass spectra show regular patterns in the ocean, and that these patterns often show systematic variations with other properties such as total biomass, nutrient concentration, season, and distance (both vertical and horizontal). The recent finding of the ubiquitous nature of layers of phytoplankton < 1 m thick prompted us to explore the fine- and microscale vertical variations of size- and fluorescence-abundance spectra in the ocean. Using a two-dimensional planar laser imaging system mounted on a free-falling platform, we quantified the properties of large fluorescent particles ( 20 μm–2 cm) through the water column, obtaining images every 10–30 cm. These images showed systematic relationships of the spectral properties to total chlorophyll: increased proportions of the smallest particles at high chlorophyll concentrations, and a lengthening of the spectral size range at high total chlorophyll concentrations (more large particles at high chlorophyll concentrations). Further, we observed significant variations of the spectral properties over scales of 1 m and less, and recorded the frequent occurrence of unusual layers of large particles. Our new instrument, which is sensitive to thin layers of enhanced phytoplankton biomass, shows the planktonic community to be highly structured vertically on scales of 1–2 m, particularly within the DCM.     

Pollution history of a tropical estuary revealed by combined hydrodynamic modelling and sediment geochemistry

Hydrodynamic modelling of water movement in Hunts Bay, a protected part of Kingston Harbour, Jamaica, shows that depth averaged tidal flows are very low. In the northeast corner of Hunts Bay, water is essentially stagnant. Even under high flow conditions, much of the Bay bottom water is ‘bypassed' by buoyant, lower salinity surface flows. The muddy sediments of Hunts Bay reflect these sluggish to stagnant conditions; sediment cores from the northeast corner of the Bay contain progressively higher amounts of organic matter in their upper parts (last 15–20 years sedimentation). Combined C/N ratios and stable carbon isotope compositions of this organic matter imply a sewage origin. Both lead and chromium metal concentrations and enrichment factors relative to average crustal shales show geographically related patterns that reflect hydrodynamic circulation predicted by modelling. In particular, metal concentrations and enrichment factors are highest at the northern end of the bay, especially in the northeast corner. Modelling confirms that stagnant conditions would occur in the northeast part of the bay even without the presence of a major causeway. The causeway may contribute to low flow conditions, but is not the principal cause of organic contamination, which is simply an excessive input of sewage.     

Circulation, hydrographic and nutrient characteristics of the Cilician Basin, Northeastern Mediterranean Sea

The water mass, circulation and chemical properties of the Cilician Basin, the northeastern Levantine Sea, are described on the basis of three hydrographic cruises performed during May 1997 (spring), July 1998 (summer) and October 2003 (autumn). The hydrographic data reveal the presence of Levantine Surface Water (LSW) and Modified Atlantic Water (MAW) within the upper 90 m layer, Levantine Intermediate Water (LIW) between 90 and 250 m, and Transitional Mediterranean Water (TMW) further below. The temporal variability of the circulation system is manifested by a change in shape, size and intensity of eddies as well as the pathways of the Lattakia Basin coastal current system. The nutrient concentrations varied between nitrate + nitrite = 0.16–0.31 μM, phosphate = 0.02–0.03 μM and silicate = 0.95–1.2 μM for the surface layer during sampling periods. Dissolved nutrient concentrations in the Transitional Mediterranean Water were: 2.1–5.3 μM for NO₃ + NO₂, 0.10–0.21 μM for PO₄ and 5.7–10 μM for Si. The molar ratios of nitrate to phosphate in the water column range between 5 and 20 in the surface layer and reach up to a value of 45 at the top of the nutricline at the depths of 29.05 kg/m³ isopycnal surface for most of the year. Below the nutricline the N / P ratios retain the values around 24–28.     

Distribution of void fraction in bubbly flow through a horizontal channel: Bubbly boundary layer flow, 2nd report

A method of enveloping the hull with a sheet of microbubbles is discussed. It forms part of a study on means of reducing the skin friction acting on a ship's hull. In this report, a bubble traveling through a horizontal channel is regarded as a diffusive particle. Based on this assumption, an equation based on flow flux balance is derived for determining the void fraction in approximation. The equation thus derived is used for calculation, and the calculation results are compared with reported experimental data. The equation is further manipulated to make it compatible with a mixing length model that takes into account the presence of bubbles in the liquid stream. Among the factors contained in the equation thus derived, those affected by the presence of bubbles are the change of mixing length and the difference in the ratio of skin friction between cases with and without bubbles. These factors can be calculated using the mean void fraction in the boundary layer determined by the rate of air supply into the flow field. It is suggested that the ratio between boundary layer thickness and bubble diameter could constitute a significant parameter to replace the scale effect in estimating values applicable to actual ships from corresponding data obtained in model experiments.List of symbols a ₁ proportionality constant indicating directionality of turbulence - B law-of-the-wall constant - C _f local skin-friction coefficient in the presence of bubbles - C _f0 local skin-friction coefficient in the absence of bubbles - d _b bubble diameter (m) - g acceleration of gravity (m/s²) - j _g flow flux of gas phase accountable to buoyancy (m/s) - j _t flow flux of gas phase accountable to turbulence (m/s) - k ₄ constant relating reduction of liquid shear stress by bubble presence to decrease of force imparted to bubble by its displacement due to turbulence - l _b mixing length of gas phase (m) - l _m mixing length of liquid phase (m) - l _mb diminution of liquid phase mixing length by bubble presence (m) - Q _G rate of air supply to liquid stream (l/min) - q _/g velocity of bubble rise (m/s) - 2R height of horizontal channel (m) - T _* integral time scale (s) - U _m mean stream velocity in channel (m/s) - U friction velocity in channel (m/s) - V volume of a bubble (m³) - u, ¯ v time-averaged stream velocities inx- andy-directions, respectively (m/s) - u, v turbulent velocity components inx- andy-directions, respectively (m/s) - v root mean square of turbulence component in they-direction (m/s) - root mean square of bubble displacement iny-direction with reference to turbulent liquid phase velocity (m) - y displacement from ceiling (m) - local void fraction - _m mean void fraction in boundary layer - _m constant relating local void fraction to law-of-the-wall constant - _t reduction of turbulent stress (N/m²) - law-of-the-wall constant in turbulent liquid region in absence of bubbles - ₁ law-of-the-wall constant in turbulent liquid region in presence of bubbles - ₂ law-of-the-wall constant in gas phase - _m constant indicating representative turbulence scale (m) - viscosity (Pa × s) - v kinematic viscosity (m²/s) - density (kg/m³) Suffixes G gas - L liquid - 0 absence of bubbles     

舰船气泡尾流数值模拟(英文)

Based on a volume of fluid two-phase model imbedded in the general computational fluid dynamics code FLUENT6.3.26, the viscous flow with free surface around a model-scaled KRISO container ship (KCS) was first numerically simulated. Then with a rigid-lid-free-surface method, the underwater flow field was computed based on the mixture multiphase model to simulate the bubbly wake around the KCS hull. The realizable k-ε two-equation turbulence model and Reynolds stress model were used to analyze the effects of turbulence model on the ship bubbly wake. The air entrainment model, which is relative to the normal velocity gradient of the free surface, and the solving method were verified by the qualitatively reasonable computed results.     

Numerical simulation of the barotropic tides in the Tunisian Shelf and the Strait of Sicily

We have investigated the barotropic tides in the Tunisian shelf and the Strait of Sicily using the Regional Ocean Modelling System (ROMS) with very high-resolution. Model performance was evaluated with respect to tide gauge, satellite data, and current meter measurements. The model fields faithfully reproduced the major feature of the barotropic tidal currents and agreed well with existing tidal elevation and phase observations. General features for the various semidiurnal constituents are nearly similar to each other with maximum amplitude in the Gulf of Gabes. The larger tidal currents occur over the continental shelves. In the Adventure Bank, the current is essentially of diurnal type whereas in the Gulf of Gabes it is of semidiurnal type.Tidal energy lost, which is primarily due to bottom stress dissipation, is predominantly in the Gulf of Gabes ( 61%), the Strait of Sicily, and the Strait of Messina. The forcing function for internal tides shows for both M₂ and K₁ constituents, significant spatial variability in the Strait of Sicily. This suggests that some internal tides will be generated in these regions and could thus explain the observed strong diurnal internal waves in the Adventure Bank.     

The numerical friction line

This paper presents a study on the numerical calculation of the friction resistance coefficient of an infinitely thin plate as a function of the Reynolds number. Seven eddy-viscosity models have been selected: the one-equation turbulence models of Menter and Spalart–Allmaras; the k-ω two-equation model proposed by Wilcox and its TNT, BSL and SST variants and the two-equation model. The flow has been computed at 14 Reynolds numbers in sets of seven geometrically similar Cartesian grids to allow a reliable estimation of the numerical uncertainty. The effect of the computational domain size has been reduced to negligible levels (below the numerical uncertainty). And the same holds for the iterative and round-off errors. In the finest grids of each set, the numerical uncertainty of the friction resistance coefficient is always below 1%. Special attention has further been given to the solution behaviour in the laminar-to-turbulent transition region. Curve fits have been applied to the data obtained at the 14 Reynolds numbers and the numerical friction lines are compared with four proposals from the open literature: the 1957 ITTC line, the Schoenherr line and the lines suggested by Grigson and Katsui et al. The differences between the numerical friction lines obtained with the seven turbulence models are smaller than the differences between the four lines proposed in the open literature.     

An approach to estimation of near-surface turbulence and CO2 transfer velocity from remote sensing data

The air–sea CO₂ exchange is primarily determined by the boundary-layer processes in the near-surface layer of the ocean since it is a water-side limited gas. As a consequence, the interfacial component of the CO₂ transfer velocity can be linked to parameters of turbulence in the near-surface layer of the ocean. The development of remote sensing techniques provides a possibility to quantify the dissipation of the turbulent kinetic energy in the near-surface layer of the ocean and the air–sea CO₂ transfer velocity on a global scale. In this work, the dissipation rate of the turbulent kinetic energy in the near-surface layer of the ocean and its patchiness has been linked to the air–sea CO₂ transfer velocity with a boundary-layer type model. Field observations of upper ocean turbulence, laboratory studies, and the direct CO₂ flux measurements are used to validate the model. The model is then forced with the TOPEX POSEIDON wind speed and significant wave height to demonstrate its applicability for estimating the distribution of the near-surface turbulence dissipation rate and gas transfer velocity for an extended (decadal) time period. A future version of this remote sensing algorithm will incorporate directional wind/wave data being available from QUIKSCAT, a now-cast wave model, and satellite heat fluxes. The inclusion of microwave imagery from the Special Sensor Microwave Imager (SSM/I) and the Synthetic Aperture Radar (SAR) will provide additional information on the fractional whitecap coverage and sea surface turbulence patchiness.