Turbulent mixing and internal tides in Gaoping (Kaoping) Submarine Canyon, Taiwan

Turbulent overturning on scales greater than 10 m is observed near the bottom and in mid-depth layers within the Gaoping (formerly spelled Kaoping) Submarine Canyon (KPSC) in southern Taiwan. Bursts of strong turbulence coexist with bursts of strong sediment concentrations in mid-depth layers. The turbulence kinetic energy dissipation rate in some turbulence bursts exceeds 10^− 4 W kg^− 1, and the eddy diffusivity exceeds 10^− 1 m² s^− 1. Within the canyon, the depth averaged turbulence kinetic energy dissipation rate is ~ 7 × 10^− 6 W kg^− 1, and the depth averaged eddy diffusivity is ~ 10^− 2 m² s^− 1. These are more than two orders of magnitude greater than typical values in the open ocean, and are much larger than those found in the Monterey Canyon where the strong turbulent mixing has also been. The interaction of tidal currents with the complex topography in Gaoping Submarine Canyon is presumably responsible for the observed turbulent overturning via shear instability and the breaking of internal tides and internal waves at critical frequencies. Strong 1st-mode internal tides exist in KPSC. The depth averaged internal tidal energy near the canyon mouth is ~ 0.17 m² s^− 2. The depth integrated internal tidal energy flux at the mouth of the canyon is ~ 14 kW m^− 1, propagating along the axis of the canyon toward the canyon head. The internal tidal energy flux in the canyon is 3–7 times greater than that found in Monterey Canyon, presumably due to the more than 10 times larger barotropic tide in the canyon. Simple energy budget calculations conclude that internal tides alone may provide energy sufficient to explain the turbulent mixing estimated within the canyon. Further experiments are needed in order to quantify the seasonal and geographical distributions of internal tides in Gaoping Submarine Canyon and their effects on the sediment flux in the canyon.     

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.     

Turbulence in coastal fronts near the mouths of Block Island and Long Island Sounds

Measurements of turbulence were performed in four frontal locations near the mouths of Block Island Sound (BIS) and Long Island Sound (LIS). These measurements extend from the offshore front associated with BIS and Mid-Atlantic Bight Shelf water, to the onshore fronts near the Montauk Point (MK) headland, and the Connecticut River plume front. The latter feature is closely associated with the major fresh water input to LIS. Turbulent kinetic energy (TKE) dissipation rate, ε, was obtained using shear probes mounted on an autonomous underwater vehicle. Offshore, the BIS estuarine outflow front showed, during spring season and ebb tide, maximum TKE dissipation rate, ε, estimates of order 10^− 5 W/kg, with background values of order 10^− 6 to 10^− 9 W/kg. Edwards et al. [Edwards, C.A., Fake, T.A., and Bogden, P.S., 2004a. Spring–summer frontogenesis at the mouth of Block Island Sound: 1. A numerical investigation into tidal and buoyancy-forced motion. Journal of Geophysical Research 109 (C12021), doi:10.1029/2003JC002132.] model this front as the boundary of a tidally driven, baroclinically adjusted BIS flow around the MK headland eddy. At the entrance to BIS, near MK, two additional fronts are observed, one of which was over sand waves. For the headland site front east of MK, without sand waves, during ebb tide, ε estimates of 10^− 5 to 10^− 6 W/kg were observed. The model shows that this front is at the northern end of an anti-cyclonic headland eddy, and within a region of strong tidal mixing. For the headland site front further northeast over sand waves, maximum ε estimates were of order 10^− 4 W/kg within a background of order 10^− 7–10^− 6 W/kg. From the model, this front is at the northeastern edge of the anti-cyclonic headland eddy and within the tidal mixing zone. For the Connecticut River plume front, a surface trapped plume, during ebb tide, maximum ε estimates of 10^− 5 W/kg were obtained, within a background of 10^− 6 to 10^− 8 W/kg. Of all four fronts, the river plume front has the largest finescale mean-square shear, S² ~ 0.15 s^− 2. All of the frontal locations had local values of the buoyancy Reynolds number indicating strong isotropic turbulence at the dissipation scales. Local values of the Froude number indicated shear instability in all of the fronts.     

The structure of dissipation in the western Irish Sea front

We report on an intensive campaign in the summer of 2006 to observe turbulent energy dissipation in the vicinity of a tidal mixing front which separates well mixed and seasonally stratified regimes in the western Irish Sea. The rate of turbulent dissipation ε was observed on a section across the front by a combination of vertical profiles with the FLY dissipation profiler and horizontal profiles by shear sensors mounted on an AUV (Autosub). Mean flow conditions and stratification were obtained from a bed mounted ADCP and a vertical chain of thermistors on a mooring. During an Autosub mission of 60 h, the vehicle, moving at a speed of ~ 1.2 m s^− 1, completed 10 useable frontal crossings between end points which were allowed to move with the mean flow. The results were combined with parallel measurements of the vertical profile of ε which were made using FLY for periods of up to 13 h at positions along the Autosub track. The two data sets, which show a satisfactory degree of consistency, were combined to elucidate the space–time variation of dissipation in the frontal zone. Using harmonic analysis, the spatial structure of dissipation was separated from the strong time dependent signal at the M₄ tidal frequency to yield a picture of the cross-frontal distribution of energy dissipation. A complementary picture of the frontal velocity field was obtained from a moored ADCP and estimates of the mean velocity derived from the thermal wind using the observed density distribution. which indicated the presence of a strong (0.2 m s^− 1) jet-like flow in the high gradient region of the front. Under neap tidal conditions, mean dissipation varied across the section by 3 orders of magnitude exceeding 10^− 2 W m^− 3 near the seabed in the mixed regime and decreasing to 10^− 5 W m^− 3. in the strongly stratified interior regime. The spatial pattern of dissipation is consistent in general form with the predictions of models of tidal mixing and does not reflect any strong influence by the frontal jet.     

Using high sampling-rate ADCP for observing vigorous processes above sloping [deep] ocean bottoms

Above sloping bottoms in the ocean mixing processes are not predominantly generated by shear-induced turbulence via bottom friction. Instead, the restratifying buoyancy forces and internal waves create a highly non-linearly varying environment including ‘stratified turbulence’. Most of the resulting vigorous mixing processes that dominate sediment resuspension occur during the passage of frontal bores or solitary boluses, ‘solibores’. Here, the observed evolution of different forms of highly non-linear strictly upslope moving ‘waves’, bores or boluses are reviewed from various NIOZ projects at deep sloping bottom sites ranging from 500 to 3000 m.Such fronts pass a fixed site within a few minutes, extending some 60 ± 30 m above the bottom and occurring over much larger periods at once per subinertial or meso-scale period or approximately, but not exactly, once per tidal harmonic period. In order to observe the details of such solibore one needs specific, high-sampling rate equipment. A suitable piece of equipment is a bottom-mounted 4-beam 300 kHz acoustic Doppler current profiler (ADCP), provided it samples at a rate of about once per second over a period of at least several weeks. Not just the three components of current velocity [u,v,w] are monitored over a range of some 80 m at 1 m intervals, but also the relative ‘echo intensity’ dI, which is a measure for suspended matter and stratified turbulence. Such ADCP-observations are combined and compared with high-resolution temperature measurements. Fine details show a turbulent inner core with more or less laminar streamlines outside it. Whether a front or a bolus, the bore is never observed as a completely closed contour, as swept up turbulent material is sucked into the core at the rear end.     

Mixing over the steep side of the Cycladic Plateau in the Aegean Sea

Intensive microstructure sampling over the southern slope of the Cycladic Plateau found very weak mixing in the pycnocline, centered on a thin minimum of diapycnal diffusivity with Kρ=1.5×¹⁰−6 m² s^− 1. Below the pycnocline, Kρ increased exponentially in the bottom 200 m, reaching 1 × 10^− 4 m² s^− 1 a few meters above the bottom. Near-bottom mixing was most intense where the bottom slope equaled the characteristic slope of the semi-diurnal internal tide. This suggests internal wave scattering and/or generation at the bottom, a conclusion supported by near-bottom dissipation rates increasing following rising winds and with intensifying internal waves. Several pinnacles on the slope were local mixing hotspots. Signatures included a vertical line of strong mixing in a pinnacle's wake, an hydraulic jump or lee wave over a downstream side of the summit, and a ‘beam’ sloping upward at the near-inertial characteristic slope. Because dissipation rate averages were dominated by strong turbulence, ?/νN² > 100, the effect on Kρ of alternate mixing efficiencies proposed for this range of turbulent intensity is explored.     

白令海潮汐能通量和底边界能耗散

The spatial distribution of the energy flux, bottom boundary layer (BBL) energy dissipation, surface elevation amplitude and current magnitude of the major semidiurnal tidal constituents in the Bering Sea are examined in detail. These distributions are obtained from the results of a three-dimensional numerical simulation model (POM). Compared with observation data from seven stations, the root mean square errors of tidal height are 2.6 cm and 1.2 cm for M₂ and N₂ respectively, and those of phase-lag are 21.8° and 15.8° respectively. The majority of the tidal energy flux off the deep basin is along the shelf edge, although some of this flux crosses the shelf edge, especially in the southeast of the shelf break. The total M₂ energy dissipation in the Bering Sea is 30.43 GW, which is about 10 times of that of N₂ and S₂. The semidiurnal tidal energy enters mainly to the Bering Sea by Samalga Pass, Amukta Pass and Seguam Pass, accounting more than 60% of the total energy entering the Being Sea from the Pacific.     

Seasonal upwelling underneath the Tsushima Warm Current along the Japanese shelf slope

Seasonal upwelling variations are examined in the eastward coastal boundary branches of the Tsushima Warm Current (TWC). The climatological pattern and the fundamental hydrographic structure of the seasonal appearance of cold water are revealed by analyzing the temperature profile data around the Japanese shelf area. Seasonal maps of temperature at the intermediate layer around 200 m depth show the rise of the main pycnocline along the Japanese coast due to seasonal subsurface cooling from May to September. The cold water areas appear around the strong curvature in the continental shelf break. These areas are confined to the south of the TWC thermal front, i.e., to the coastal boundary region. The seasonal appearance of the localized cooling areas implies that the seasonal upwelling is induced by horizontal variations in shelf topography and the intensifying TWC from May to September.     

Internal tidal currents in the Gaoping (Kaoping) Submarine Canyon

Data from five separate field experiments during 2000–2006 were used to study the internal tidal flow patterns in the Gaoping (formerly spelled Kaoping) Submarine Canyon. The internal tides are large with maximum interface displacements of about 200 m and maximum velocities of over 100cm/s. They are characterized by a first-mode velocity and density structure with zero crossing at about 100 m depth. In the lower layer, the currents increase with increasing depth. The density interface and the along-channel velocity are approximately 90° out-of-phase, suggesting a predominant standing wave pattern. However, partial reflection is indicated as there is a consistent phase advance between sea level and density interface along the canyon axis.     

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.     

A 3D, finite element model for baroclinic circulation on the Vancouver Island continental shelf

This paper describes the development and application of a 3-dimensional model of the barotropic and baroclinic circulation on the continental shelf west of Vancouver Island, Canada. A previous study with a 2D barotropic model and field data revealed that several tidal constituents have a significant baroclinic component (the K1 in particular). Thus we embarked on another study with a 3D model to study the baroclinic effects on the residual and several selected tidal constituents.The 3D model uses a harmonic expansion in time and a finite element discretization in space. All nonlinear terms are retained, including quadratic bottom stress, advection and wave transport (continuity nonlinearity). The equations are solved as a global and a local problem, where the global problem is the solution of the wave equation formulation of the shallow water equations, and the local problem is the solution of the momentum equation for the vertical velocity profile. These equations are coupled to the advection-diffusion equation for density so that density gradient forcing is included in the momentum equations. However, the study presented here describes diagnostic calculations for the baroclinic residual circulation only.The model is sufficiently efficient that it encourages sensitivity testing with a large number of model runs. In this sense, the model is akin to an extension of analytical solutions to the domain of irregular geometry and bottom topography where this parameter space can be explored in some detail.In particular, the consequences of the sigma coordinate system used by the model are explored. Test cases using an idealized representation of the continental shelf, shelf break and shelf slope, lead to an estimation of the velocity errors caused by interpolation errors inherent in the sigma coordinate system. On the basis of these estimates, the computational grid used in the 2D model is found to have inadequate resolution. Thus a new grid is generated with increased accuracy in the region of the shelf break. However, even with increased resolution, spurious baroclinic circulation seaward of the shelf break and in the vicinity of Juan de Fuca canyon remained a significant problem when the pressure gradient terms were evaluated using the σ coordinate system and using a realistic density profile.With the new grid, diagnostic calculations of the barotropic and baroclinic residual circulation are performed using forcing from the observed σ_t (density) field and from the gradient of this field.     

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.     

Cross-shelf circulation and phytoplankton distribution at the summertime New England shelfbreak front

We investigate aspects of the secondary (cross-shelf) circulation at the Middle Atlantic Bight shelfbreak front using high-resolution data collected on the New England Shelf in August 2002. The alongshore shelfbreak jet coincides with the front at the seaward edge of the cold pool (remnant winter shelf water) and there is a suggestion of a cross-stream convergence centered at the jet core. Despite indications of convergence we found no evidence of a surface subduction on the seaward side of the front. At depth 70 m near the shelfbreak there was a patch of chlorophyll, located within a local temperature–salinity maximum which, though significantly below the euphotic zone, appeared to be photo-acclimated and viable. The chlorophyll feature could be the result of a local subduction by a larger scale eddy circulation seaward of the front.Dye tracer experiments directly observed the convergence at the foot of the shelfbreak front and subsequent upwelling of bottom boundary layer water along the shoreward side of the shelfbreak front. But, we found no evidence that this upwelling influenced productivity at the front. Further, since there was no cross-shelf maximum in subsurface chlorophyll at the front, we conclude that this productivity is in general, sustained by a broadly distributed local vertical nutrient flux from an underlying nutrient reservoir.     

Numerical simulations of seasonal and interannual variability of the Black Sea thermohaline circulation

The Black Sea general circulation is simulated by a primitive equation model with active free surface. The forcing is seasonally variable and is based on available climatic data. The model reproduces the main features of the Black Sea circulation, including the river discharge effects on the mean sea level and the Bosphorus outflow. Model results show that the simulated sea surface elevation increases in spring over the whole sea, reaching a maximum in the Danube delta area. In the same region, a minimum is observed in winter. The amplitude of the seasonal oscillations (about 8–12 cm over the whole basin) is of the same order of magnitude as the maximum horizontal variations (about 15–18 cm between the coastal areas and the basin interior). This strong signal formed mostly by river discharges, along with the seasonal variability in the other forcing functions and the local dynamics creates a well-pronounced interannual variability. The performance of the model in simulating the seasonal and interannual variability is critically analyzed, with a special attention on the cold intermediate water formation and the circulation in the upper 150 m. The simulations demonstrate that the source of intermediate waters is on the shelf, and that the water mass in the core of cold intermediate layer changes with time as a response to the periodic forcing at sea surface. This type of variability is characterized by pronounced interannual changes, proving that important differences could exist between water mass structure in different years, even when using identical atmospheric forcings each year.     

Modern accumulation rates and a budget of sediment off the Gaoping (Kaoping) River, SW Taiwan: A tidal and flood dominated depositional environment around a submarine canyon

Ninety-two box cores collected during 2004–2006 from an area of ~ 3000 km² off the Gaoping (formerly spelled Kaoping) River, SW Taiwan, were analyzed for fallout radionuclides (²¹⁰Pb, ¹³⁷Cs and ⁷Be) to elucidate sedimentation rates and processes, and for the calculation of a sediment budget. The study area is located at an active collision margin with a narrow shelf and a submarine canyon extending essentially into the river's mouth. The results indicate fairly constant hemipelagic sedimentation in much of the open margin and for most of the time except in the inner shelf and along the axis of the canyon where sediment transport is more dynamic and is controlled by tidal current and wave activities constantly, and by fluvial floods or gravity-driven flows episodically. Sedimentation rates in the study area derived from ²¹⁰Pb and constrained by ¹³⁷Cs vary from 0.04 to 1.5 cm/yr, with the highest rates (> 1 cm/yr) flanking the Gaoping canyon over the upper slope (200–600 m) and the lowest rates (< 0.1 cm/yr) in the distal basin beyond the continental slope. The depocenter delineated from ²¹⁰Pb-based sedimentation rates overlaps with the area covered by a flood layer resulting from super-typhoon Haitang in July 2005. Such correspondence supports the notion that the processes operating on event timescale have bearing on the formation of the sediment strata over centennial or longer timescales.From the distribution of sedimentation rates, sediment deposited in the study area annually is estimated to be 6.6 Mton/yr, accounting for less than 20% of Gaoping River's sediment load. The calculated budget, coupled with the presence of the short-lived ⁷Be and non-steady-state distribution of low levels of ²¹⁰Pb in sediments along the canyon floor, suggests rapid transport of sediment from Gaoping River's mountainous watershed (the source) via the Gaoping (Kaoping) Submarine Canyon and adjacent channels (as the conduit and temporary sink) to the abyssal plain and the Manila Trench in the South China Sea (the ultimate sink).     

On the distribution of dissolved lead in the Loire estuary and the North Biscay continental shelf, France

The dissolved lead was studied in the whole salinity gradient of the system composed of the Loire estuary and the North Biscay continental shelf. About 130 samples were collected in winter 2001 and spring 2002 during Nutrigas and Gasprod campaigns (Programme PNEC-Golfe de Gascogne, RV Thalassa) and metal measurements were conducted on board by Potentiometric Stripping Analysis. In the Loire estuary, levels of dissolved lead ranged from 0.15 to 0.24 nM and from 0.04 to 0.26 nM in winter and spring, respectively. Compared to the concentrations reported in 1987 and 1990 (0.4–1.7 nM; Boutier, B., Chiffoleau, J.F., Auger, D., Truquet, I., 1993. Influence of the Loire river on dissolved lead and cadmium concentrations in coastal waters of Brittany. Estuar. Coast. Shelf S., 36:133–143, Estuarine, Coastal and Shelf Science 36, 133–143) our study indicated much lower values. The fall in concentration in the estuary could be attributed to the stopping of activity of Octel, a big manufacturer of tetra alkyl lead. Discharge in dissolved metal to the continental shelf by the Loire river was assessed as 7.5 and 1.9 kg day^− 1 for winter and spring, respectively. On the continental shelf, levels of dissolved lead varied within 0.06 and 0.27 nM in winter (0.15 ± 0.06 nM, sd = 1.96, n = 49), whereas concentrations measured in spring were in the range 0.06–0.17 nM (0.09 ± 0.03 nM, sd = 1.96, n = 60). This difference in metal concentration was related to the amounts of rainfall that have fallen over the continental shelf: estimations of inputs by this way (74 and 32 kg day^− 1 in winter and spring, respectively) appeared to be significantly higher than inputs from the Loire river (7.5 and 1.9 kg day^− 1 in winter and spring, respectively). The distributions of dissolved metal in the surface waters highlighted the role of suspended particular matter (SPM) for a rapid “trapping” of lead near the mouth of the estuary. The vertical distributions showed, in the stratified area, a biological transfer of lead between winter and spring from surface waters to the halocline.     

Dynamics of suprabenthos-zooplankton communities around the Balearic Islands (western Mediterranean): Influence of environmental variables and effects on the biological cycle of Aristeus antennatus

Dynamics of suprabenthos and zooplankton were analyzed in two areas located in the NW (off Sóller harbour) and S (off Cabrera Archipelago) of Mallorca (Balearic Islands, western Mediterranean) at depths ranging between 135–780 m. Four stations situated respectively at 150 m (shelf-slope break), and at bathyal depths of 350, 650 and 750 m were sampled at bi-monthly intervals during six cruises performed between August 2003 and June 2004. Suprabenthos showed maximum biomass in both areas from late spring to summer (April to August), while minimum biomass was found in autumn (September–November). Though variable, temporal dynamics of zooplankton showed peaks of biomass in late winter and summer (February and June), while minimals occurred in autumn (August–September) and, at bathyal depths, in April. Suprabenthos (abundance; MDS analyses) showed a sample aggregation as a function of depth (3 groups corresponding to the shelf-slope break, upper slope — over 350 m; and the middle, deeper part of the slope — over 650–750 m), without any separation of hauls by season. By contrast, zooplankton samples were separated by season and not by depth. There was evidence of three seasonal groups corresponding to summer (June 2004–August 2003), autumn–winter (September and November 2003, February 2004), and spring (April 2004), being especially well established off Sóller. In general, suprabenthos was significantly correlated with the sediment variables (e.g. total organic matter content (% OM), potential REDOX), whereas zooplankton was almost exclusively dependent on Chl a at the surface, which suggests two different food sources for suprabenthos and zooplankton. The increase of suprabenthos abundance in April–June was paralleled by a sharp increase (ca. 2.8 times) in the %OM on sediment during the same period, coupled ca. 1–2 months of delay with the peak of surface Chl a recorded in February–March (from satellite imagery data). Suprabenthos biomass was also correlated with salinity close to the bottom, suggesting a link between suprabenthos abundance and changes in the oceanographic condition of water masses close to the bottom. It is suggested that a higher suprabenthos biomass recorded off Sóller in comparison to that off Cabrera in June could, in turn, be related to a seasonal inflow of Levantine Intermediate Water (LIW) in April–June in this area at mid bathyal depths (350–650 m). This trend would be based on: 1) it was evident only at mid-slope depths between 350–750 m, coinciding with the LIW distribution, and 2) it was not recorded among zooplankton (collected throughout the water column). The possible effect of the fluctuations of suprabenthos and zooplankton on higher trophic levels has been explored studying the diet and food consumption rates of the red shrimp Aristeus antennatus, as indicator species by its dominance in bathyal communities. A. antennatus increased its food consumption from February to April–June 2004 off Sóller, which in the case of large (CL > 40 mm) specimens was found in both areas. In addition, there was a shift of diet from winter to spring–early summer. In this last period, A. antennatus preyed upon euphausiids and mesopelagic decapods and fish, while benthos (e.g. polychaetes and bivalves) decreased in the diet. This indicates an increase in the food consumption and probably in the caloric content of the diet in pre-spawning females in April–June 2004, which is synchronized with the period when gonad development begins in A. antennatus females (May–June). Anyway, macrozooplankton, and not suprabenthos, is crucial as a high energetic food source in the coupling between food intake and reproduction in the red shrimp.     

Asymmetry of Columbia River tidal plume fronts

Columbia River tidal plume dynamics can be explained in terms of two asymmetries related to plume-front depth and internal wave generation. These asymmetries may be an important factor contributing to the observed greater primary productivity and phytoplankton standing crop on the Washington shelf. The tidal plume (the most recent ebb outflow from the estuary) is initially supercritical with respect to the frontal internal Froude number F_R on strong ebbs. It is separated from the rotating plume bulge by a front, whose properties are very different under upwelling vs. downwelling conditions. Under summer upwelling conditions, tidal plume fronts are sharp and narrow (< 20–50 m wide) on their upwind or northern side and mark a transition from supercritical to subcritical flow for up to 12 h after high water. Such sharp fronts are a source of turbulent mixing, despite the strong stratification. Because the tidal plume may overlie newly upwelled waters, these fronts can mix nutrients into the plume. Symmetry would suggest that there should be a sharp front south of the estuary mouth under summer downwelling conditions. Instead, the downwelling tidal plume front is usually diffuse on its upstream side. Mixing is weaker, and the water masses immediately below are low in nutrients. There is also an upwelling–downwelling asymmetry in internal wave generation. During upwelling and weak wind conditions, plume fronts often generate trains of non-linear internal waves as they transition from a supercritical to a subcritical state. Under downwelling conditions, internal wave release is less common and the waves are less energetic. Furthermore, regardless of wind conditions, solition formation almost always begins on the south side of the plume so that the front “unzips” from south to north. This distinction is important, because these internal waves contribute to vertical mixing in the plume bulge and transport low-salinity water across the tidal plume into the plume bulge.F_R and plume depth are key parameters in distinguishing the upwelling and downwelling situations, and these two asymmetries can be explained in terms of potential vorticity conservation. The divergence of the tidal outflow after it leaves the estuary embeds relative vorticity in the emerging tidal plume water mass. This vorticity controls the transition of the tidal plume front to a subcritical state and consequently the timing and location of internal wave generation by plume fronts.     

Analysis of the hydrophysical structure of the Sea of Azov in the period of the bottom anoxia development

The hydrophysical and hydrochemical structure of the Sea of Azov, with developed bottom anoxia, was studied during the RV “Akvanavt” cruise from July 31 to August 03, 2001. The anoxic zone with a thickness from 0.5 to 4 m above the bottom was found in all deep regions of the Sea. Concentrations of hydrochemical parameters were similar to the pronounced anoxic conditions (about 90 mmol m^− 3 of hydrogen sulfide, 17 mmol m^− 3 of ammonia, 6 mmol m^− 3 of phosphate, 7 mmol m^− 3 of total manganese). The hydrophysical structure was characterized by the uniform distribution of temperature in the upper 6–7 m mixed layer (UML). Below this a thin (0.4–0.8 m) thermocline layer was observed, just above the anoxic waters. Formation of this phenomenon was connected with that summer weather conditions. Intensive rains led to increased influx of river waters in June. That resulted in large input of allochtonous organic matter (OM) and inorganic nutrients; the latter were consumed on the additional autochthonous organic matter production. In July the weather was characterized by a significant rise in the daily averaged air temperature and large oscillations of temperature during the day. In this period a wind of constant direction was absent, but wind bursts were observed. The completed analyses showed that the formation of such a structure could be connected with the following factors: (i) positive growth trends of the daily averaged temperature and the daily oscillations of temperature, (ii) presence of wind bursts. The joint action of these factors resulted in the formation of the UML. The amplitude of wind bursts determined the depth of UML, and the value of trend determined the value of the temperature change in the thermocline. An initial presence of bottom halocline (caused by the Black Sea water influx to the bottom of the Sea of Azov) prevented the heating of the bottom layer and therefore led to an increase of vertical gradient of temperature in the thermocline. The spatial distribution of the turbulent exchange coefficient confirmed the existence of a “stagnation” area located above the anoxia zone, which is also, apparently, the reason for its occurrence.