The biogeochemical cycling of methane in Ria de Vigo, NW Spain: Sediment processing and sea–air exchange

Methane (CH₄) concentrations were measured in the water column, in sediment porewaters, and in atmospheric air, in the Ría de Vigo, NW Spain, during both the onset (April 2003) and at the end of (September 2004) seasonal upwelling. In addition, CH₄ concentration and stable isotopic signatures (δ¹³CH₄) were measured in porewaters, and sediment methanogenesis and aerobic oxidation of CH₄ were determined in sediment incubations. Surface water column CH₄ (2 m depth) was in the range 3–180 nmol l^− 1 (110–8500% saturation) and followed a generally landward increase but with localised maxima in both the inner and middle Ría. These maxima were consistent with CH₄ inputs from underlying porewaters in which CH₄ concentrations were up to 3 orders of magnitude higher (maximum 350 μmol l^− 1). Surface water CH₄ concentrations were approximately three times higher in September than in April, consistent with a significant benthic CH₄ flux driven by enhanced sediment methanogenesis following the summer productivity maximum. CH₄ and δ¹³CH₄ in sediment porewaters and in incubated sediment slurries (20 °C) revealed significant sediment CH₄ oxidation, with an apparent isotopic fractionation factor (r_c) of  1.004. Using turbulent diffusion models of air–sea exchange we estimate an annual emission of atmospheric CH₄ from the Ría de Vigo of 18–44 × 10⁶ g (1.1–2.7 × 10⁶ mol). This estimate is approximately 1–2 orders of magnitude lower than a previous estimate based on a bubble transport model.     

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.     

Suspended sediment and erosion dynamics in Kugmallit Bay and Beaufort Sea during ice-free conditions

The Mackenzie River is the largest river on the North American side of the Arctic and its huge freshwater and sediment load impacts the Canadian Beaufort Shelf. Huge quantities of sediment and associated organic carbon are transported in the Mackenzie plume into the interior of the Arctic Ocean mainly during the freshet (May to September). Changing climate scenarios portend increased coastal erosion and resuspension that lead to altered river-shelf-slope particle budgets. We measured sedimentation rates, suspended particulate matter (SPM), particle size and settling rates during ice-free conditions in Kugmallit Bay (3–5 m depth). Additionally, measurements of erosion rate, critical shear stress, particle size distribution and resuspension threshold of bottom sediments were examined at four regionally contrasting sites (33–523 m depth) on the Canadian Beaufort Shelf using a new method for assessing sediment erosion. Wind induced resuspension was evidenced by a strong relationship between SPM and wind speed in Kugmallit Bay. Deployment of sediment traps showed decreasing sedimentation rates at sites along an inshore–offshore transect ranging from 5400 to 3700 g m^− 2 day^− 1. Particle settling rates and size distributions measured using a Perspex settling chamber showed strong relationships between equivalent spherical diameter (ESD) and particle settling rates (r² = 0.91). Mean settling rates were 0.72 cm s^− 1 with corresponding ESD values of 0.9 mm. Undisturbed sediment cores were exposed to shear stress in an attempt to compare differences in sediment stability across the shelf during September to October 2003. Shear was generated by vertically oscillating a perforated disc at controlled frequencies corresponding to calibrated shear velocity using a piston grid erosion device. Critical (Type I) erosion thresholds (u) varied between 1.1 and 1.3 cm s^− 1 with no obvious differences in location. Sediments at the deepest site Amundsen Gulf displayed the highest erosion rates (22–54 g m^− 2 min^− 1) with resuspended particle sizes ranging from 100 to 930 µm for all sites. There was no indication of biotic influence on sediment stability, although our cores did not display a fluff layer of unconsolidated sediment. Concurrent studies in the delta and shelf region suggest the importance of a nepheloid layer which transports suspended particles to the slope. Continuous cycles of resuspension, deposition, and horizontal advection may intensify with reduction of sea ice in this region. Our measurements coupled with studies of circulation and cross-shelf exchange allow parameterization and modeling of particle dynamics and carbon fluxes under various climate change scenarios.     

Characteristics of bubble plumes, bubble-plume bubbles and waves from wind-steepened wave breaking

Observations of breaking waves, associated bubble plumes and bubble-plume size distributions were used to explore the coupled evolution of wave-breaking, wave properties and bubble-plume characteristics. Experiments were made in a large, freshwater, wind-wave channel with mechanical wind-steepened waves and a wind speed of 13 m s^− 1. Bubble plumes exhibited a wide range of bubble distributions, physical extent and dynamics. A classification scheme was developed based on plume extent and “optical density” which is the ability of a plume to optically obscure the image of the background until maximum penetration of the plume. Plumes were classified as either dense (obscure) or diffuse (no-obscure). For each class, the plume bubble population size distribution, Φ(r,t), where r is the bubble radius and t the time, was determined. Dense plumes have a large radius peak in Φ and thus are enhanced in large bubbles. Diffuse plumes are well-described by a weakly size decreasing Φ(r,t) for r < 1000 μm and a more strongly size decreasing Φ(r,t) for r > 1000 μm.The bubble-plume formation rate, P, for each class, wave-breaking rate and wave characteristics were measured with respect to fetch. Wave-breaking rate and intensity are strongly fetch-dependent. In general, the trends in P and wave breaking are similar, reaching a maximum at the fetch of maximum wave breaking. The ratio of P for dense to diffuse plumes is even more sensitive to the occurrence of the most intense wave breaking, where dense plume formation is the greatest.Using P and the bubble size population distributions for each plume class, the global bubble-plume, injection size distribution, Ψ_i(r), was calculated. The volume injection rate for the study area was 640 cm³ s^− 1 divided approximately equally between bubbles smaller and larger than r  1700 μm.     

Carbon and nutrient dynamics in a hypertrophic lagoon in southwestern Taiwan

Dissolved and particulate phases of carbon (DIC, DOC, POC) and nutrients (DIN, DIP, DSi, DON, DOP, PN) were investigated bimonthly from August 1999 to August 2000 to study biogeochemical dynamics of carbon and nutrients in Tapong Bay, a small semi-enclosed and hypertrophic lagoon in southwestern Taiwan. The lagoon has only a tidal inlet for exchanging water between Tapong Bay and Taiwan Strait, which may result in low water exchange rates and various oxygen-deficient conditions in bottom water of the inner bay during warm seasons. The water exchange time of Tapong Bay ranges from 7 days (summer) to 13 days (winter) with a mean of 10 days. Nutrient dynamics were largely ascribed to allochthonous inputs, biological and exported removals in the lagoon. Diffusion fluxes from sediments to overlying water accounted for only about 7.6% of annual DIN inputs and 1.0% of annual DIP inputs. High primary productivity (89 mol C m⁻² year⁻¹) supported by abundant nutrients primarily drove the lagoon into a hypertrophic condition as particulate organic matter was derived mainly from biological production. Excess of DIP appeared to occur throughout the study period in the lagoon. Temperature, solar radiation and turbidity, rather than nutrients, perhaps controlled seasonal variations of primary productivity. The net ecosystem production (NEP) derived from daily changes of DOC and POC inventories was about 6.3 mmol C m⁻² day⁻¹ that was close to 6.7 mmol C m⁻² day⁻¹ simulated from the biogeochemical modeling. Therefore, the net ecosystem production (NEP) rate of organic carbon estimated from the biogeochemical model was reliable, and the NEP was temporally variable with an annual mean of 5.8 mol C m⁻² year⁻¹, implying that Tapong Bay was an autotrophic system. Although calcification proceeded pronouncedly in warm seasons, an invasion of CO₂ was significant in this system. In terms of nitrogen budget, the annual nitrogen fixation exceeded the annual denitrification with a magnitude of 1.30 mol N m⁻² year⁻¹, which may be supported by the abundance of nitrogen fixation microplanktons in the lagoon.     

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.     

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.     

Siliceous plankton dominate primary and new productivity during the onset of El Niño conditions in the Santa Barbara Basin, California

The potential for carbon export and the role of siliceous plankton in the cycling of C and N was assessed in natural plankton assemblages in the Santa Barbara Basin, California, by examining uptake rates of inorganic carbon, nitrate and silicic acid. In April–August 1997, the concentrations of chlorophyll a, particulate organic carbon, particulate organic nitrogen and biogenic silica were measured twice monthly, and results revealed the occurrence of at least three blooms, the largest in June. Particulate elemental ratios of C, N and Si were similar to ratios of nutrient-replete diatoms, suggesting that they dominated this bloom. Mean integrated rates of carbon, nitrate and silicon uptake during the 4-month study period are similar to other productive coastal and upwelling regions (103, 8.3 and 13 mmol m⁻² day⁻¹, respectively). New production rates were twice as high as previously reported in this region and indicate that high rates of new production along eastern boundary currents are not confined to the major coastal upwelling regions. C/NO₃⁻, Si/NO₃⁻ and Si/C uptake ratios varied widely, and mean integrated ratios were 14±5.4, 1.6±1.0 and 0.12±0.07 (S.D.), respectively. That mean C/NO₃⁻ uptake ratio corresponds to an f-ratio of about 0.5 indicating a large potential for particulate export. Based on the average Si/NO₃⁻ and Si/C uptake ratios, diatoms could perform all of the primary production and nitrate uptake that occurred during the study; these rates also suggest that export is controlled by diatoms in this system. The mean Si/C biomass ratio was lower than the mean Si/C uptake ratio, consistent with the preferential export of Si relative to C observed in sediment traps in the basin. The study took place during a period of surface-water warming, with nitrate and silicic acid concentrations decreasing throughout the onset of the 1997–1998 El Niño conditions. Although diatoms contributed less to particulate biomass during the low nutrient conditions, high f-ratios (0.33–0.66) were maintained.     

Influence of the Mackenzie River plume on the sinking export of particulate material on the shelf

We examined the influence of the Mackenzie River plume on sinking fluxes of particulate organic and inorganic material on the Mackenzie Shelf, Canadian Arctic. Short-term particle interceptor traps were deployed under the halocline at 3 stations across the shelf during fall 2002 and at 3 stations along the shelf edge during summer 2004. During the two sampling periods, the horizontal patterns in sinking fluxes of particulate organic carbon (POC) and chlorophyll a (chl a) paralleled those in chl a biomass within the plume. Highest sinking fluxes of particulate organic material occurred at stations strongly influenced by the river plume (maximum POC sinking fluxes at 25 m of 98 mg C m^− 2 d^− 1 and 197 mg C m^− 2 d^− 1 in 2002 and 2004, respectively). The biogeochemical composition of the sinking material varied seasonally with phytoplankton and fecal pellets contributing considerably to the sinking flux in summer, while amorphous detritus dominated in the fall. Also, the sinking phytoplankton assemblage showed a seasonal succession from a dominance of diatoms in summer to flagellates and dinoflagellates in the fall. The presence of the freshwater diatom Eunotia sp. in the sinking assemblage directly underneath the river plume indicates the contribution of a phytoplankton community carried by the plume to the sinking export of organic material. Yet, increasing chl a and BioSi sinking fluxes with depth indicated an export of phytoplankton from the water column below the river plume during summer and fall. Grazing activity, mostly by copepods, and to a lesser extent by appendicularians, appeared to occur in a well-defined stratum underneath the river plume, particularly during summer. These results show that the Mackenzie River influences the magnitude and composition of the sinking material on the shelf in summer and fall, but does not constitute the only source of material sinking to depth at stations influenced by the river plume.     

Spatial variability of the barotropic M2 constituent tidal current over the Rías Baixas Galician shelf (NW Spain)

Using current meters and systematically repeated vessel mounted ADCP data legs, we describe the role of the semidiurnal and diurnal constituents in the tidal currents and investigate the mesoscale variability of the M₂ tidal currents field along three cross-shelf legs over the Galician shelf. We found that the estimated values of the semi-major axes of the main semidiurnal tidal constituents are closely related to those obtained in previous research and that the tidal currents are predominantly semidiurnal. Amplitudes of the M₂ and S₂ semi-major axes constituents are comparable, and important fortnightly variability should thus be expected in the amplitude of the semidiurnal tidal currents. Vertical profiles of the semi-major axes of the main tidal constituents reveal that, in the absence of stratification, and far from the bottom frictional effects, the tidal currents are mostly barotropic. Over the slope, the tidal ellipses tilt from a north–south orientation toward the northeast–southwest direction and, as we go further into the shelf, their amplitude is increased, a feature that is reasonably explained by the transition from the dynamics of a Poincaré wave offshore to the dynamics of a northward alongshore propagating Kelvin wave over the shelf. The Ría de Vigo exerts a notable influence over the tidal currents, acting as an extension of the shelf and introducing a meridional variation in the orientation of the ellipses from the alongshore direction in the southern straight coast toward the northeast–southwest direction at the latitude of the Rías. Finally, the spatial fitting methods we tried have proved to be successful in VMADCP data detiding in this region.     

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.     

Estimates of vertical mixing during a Lagrangian experiment off the Galician coast

An adjoint 1-D model was used to determine vertical diffusivity coefficients from temperature profiles collected within a filament escaping from the Galician coast following an upwelling event. The optimisation scheme ended with relatively high diffusivity values within the thermocline (9×10⁻⁵ m² s⁻¹). Such high values are relevant for biogeochemical exchanges between surface and deep waters in stratified areas.The optimised values were several orders of magnitude higher than the bulk of diffusivity measurements recorded with a free-falling device; however, the optimisation solution was consistent with the arithmetic mean of the measurements in the thermocline (7.7×10⁻⁵ m² s⁻¹), giving more weight to the few largest values. Below the thermocline, the data assimilation method failed because of the three-dimensional nature of the advective field of the upwelling system. Ignoring this advective forcing in the model led to estimates that were two orders of magnitude too high.The results suggest that turbulent mixing is a random process where a few intense events determine the average mixing that drives the long-term evolution of the water column structure. This statistical property is very important when one wants to use instantaneous diffusivity measurements for modelling purposes.     

Seasonality of picophytoplankton chlorophyll a and biomass in the central Cantabrian Sea, southern Bay of Biscay

Seasonal changes in the abundance and biomass of cyanobacteria (Synechococcus and Prochlorococcus) and picoeukaryotes were studied by flow cytometry in the upper layers of the central Cantabrian Sea continental shelf, from April 2002 to April 2006. The study area displayed the typical hydrographic conditions of temperate coastal zones. A marked seasonality of the relative contribution of prokaryotes and eukaryotes was found. While cyanobacteria were generally more abundant for most of the year (up to 2.4 10⁵ cells mL^− 1), picoeukaryotes dominated the community (up to 10⁴ cells mL^− 1) from February to May. The disappearance of Prochlorococcus from spring through summer is likely related to shifts in the prevailing current regime. The maximum total abundance of picophytoplankton was consistently found in late summer–early autumn. Mean photic-layer picoplanktonic chlorophyll a ranged from 0.06 to 0.53 µg L^− 1 with a relatively high mean contribution to total values (33 ± 2% SE), showing maxima around autumn and minima in spring. Biomass (range 0.58–40.16 mg C m^− 3) was generally dominated by picoeukaryotes (mean ± SE, 4.28 ± 0.27 mg C m^− 3) with an average contribution of cyanobacteria of 30 ± 2%. Different seasonality of pigment and biomass values resulted in a clear temporal pattern of picophytoplanktonic carbon to chlorophyll a ratio, which ranged from 10 (winter) to 140 (summer). This study highlights the important contribution of picoplanktonic chlorophyll a and carbon biomass in this coastal ecosystem.     

Distributions of oxygen and carbon stable isotopes and CFC-12 in the water masses of the Southern Ocean at 30°E from South Africa to Antarctica: results of the CIVA1 cruise

This study presents oceanic distributions of stable isotopes (δ¹⁸O of water and δ¹³C of ΣCO₂) and CFC-12 from samples collected during the CIVA1 cruise (February/March 1993), across the Southern Ocean, along a meridian section at 30°E, from South Africa (44°S) to Antarctica (70°S). The isotopic measurements show important variations between the subantarctic surface waters with low δ¹⁸O–high δ¹³C values and the antarctic surface waters with very low δ¹⁸O–low δ¹³C values. The surface distributions of δ¹³C values follow the major frontal oceanic structures; the vertical distribution shows the progressive upwelling from the subantarctic zone to the antarctic divergence of ¹³C-depleted CO₂ derived from remineralization of organic matter. Along the Antarctic continental shelf, between 2500 and 4000 m, a core of water with δ¹⁸O values close to −0.1‰ is associated with a relative maximum in CFC-12 concentration, although this core is not detected by its temperature and salinity parameters. This water mass, which corresponds to recently formed deep water, may originate from the eastward extension of the Weddell gyre or from bottom waters coming from the East and formed near Prydz Bay.     

A study of the response of chlorophyll-a biomass to a winter upwelling event off Western Iberia using SeaWiFS and in situ data

Observations of a winter upwelling event off Western Iberia shelf/slope in the area of influence of the Western Iberia Buoyant Plume (WIBP) were conducted in February 2000. Spatial patterns and time evolution of the chlorophyll-a (chl-a) biomass are analysed, based on in situ and satellite data. SeaWiFS-derived chl-a concentration L2 products were used to track the chlorophyll front and estimate its westward migration velocity (maximum up to 29 km day⁻¹), as well as to characterize the frontal system and its evolution. A method associating the type of spectral signature of a pixel to the fraction of chlorophyll probed by SeaWiFS enabled the estimation of the chl-a biomass within error intervals. High chlorophyll concentrations (for wintertime) were observed over the shelf and slope, up to large distances to the coast. Due to the WIBP, a shallow Ekman layer developed, being nearly coincident with the stratified upper meters. The transport comprised westward advection and stretching of the plume, with little entrainment with the offshore deep mixed layer waters. The relative enlargement of the total area of the Inside-Front Zone (IFZ) during the upwelling event was roughly accompanied by the maintenance of the average biomass per unit of area, considering the water column up to depths of interest. This suggests that there was a net increase of chl-a biomass inside the water column associated with the IFZ, roughly proportional to the increase in the IFZ area. Retention of phytoplankton in the shallow stratified nutrient-rich waters of the WIBP was a key factor for this increase in chl-a biomass.     

Sequential data assimilation in an upwelling influenced estuary

We have evaluated the impact of assimilating chlorophyll, nitrate, phosphate, silicate and ammonium into a coupled 1D hydrodynamic ecosystem model (GOTM-ERSEM) in an upwelling influenced estuary. The assimilation method chosen is the Ensemble Kalman Filter (EnKF), which has been demonstrated to improve field estimates of key variables (chlorophyll, nutrients) for bulk algal bloom prediction. The 1D model has been set up for a central station inside the Ría de Vigo (Spain). Data from bi-weekly surveys are used to constrain the model. Temperature and salinity profiles are used to ensure the correct representation of the water structure through a relaxation scheme. Chlorophyll extracts and nutrients at three depths are assimilated sequentially during 1 year simulation (1991). The assimilation period includes episodes of active upwelling and downwelling. All five assimilated variables are successfully constrained and represent a large improvement on the reference simulation (without assimilation). Small divergences can be related to poorly resolved physical processes in the model. The assimilation was further evaluated by comparing observed biomass partitioning with model results. Diatoms accounted for the largest biomass update and the largest improvement in terms of percentage of variance explained (R²). This is particularly significant as they represent the 46% of the yearly integrated observed biomass of the planktonic autotrophs. Nonetheless the R² value was low for all phytoplankton groups. Bacteria and nanoflagellates showed an improvement with respect to their yearly Root Mean Square (RMS), while the other functional groups worsen or remained unaffected. Chlorophyll assimilation was responsible for most of the impact on the phytoplankton biomass with small contributions from the silicate. It had minor impact on the updates of nutrients which in turn corrected the state variables related to the detrital pool. In this current setting, combined assimilation of chlorophyll and nutrients is not sufficient to produce a skillful simulation of the phytoplankton succession.     

Production of Acartia omorii (Copepoda: Calanoida) in Ilkwang Bay, southeastern coast of Korea

Production of the marine calanoid copepod Acartia omorii was measured from 2 October 1991 to 8 October 1992 at a station in Ilkwang Bay on the southeastern coast of Korea. A. omorii (nauplii + copepodites + adults) were present in the plankton throughout the year, with seasonal variation in abundance. Biomass of A. omorii was averaged at 0.44 mgC m^− 3, with peaks in February and July, and relatively low biomass in late summer and fall. Egg production rate ranged from 2.4 to 151.9 μgC m^− 3 day^− 1, which was equivalent to 95–6075 eggs m^− 3 day^− 1. Fecundity of an adult female was averaged at 38 eggs female^− 1 day^− 1. Instantaneous growth rates of copepodites were higher than those of nauplii stages. Annual production of A. omorii ranged from 33.5 mgC m^− 3 year^− 1 to 221 mgC m^− 2 year^− 1, showing a seasonal variation of daily production rate with peaks in February and July. The daily production rate of A. omorii was significantly correlated with chlorophyll a concentration. These results suggest that standing stocks and/or productivity of phytoplankton are the major influencing factors, rather than water temperature for the seasonal variation of production of A. omorii in Ilkwang Bay.     

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.     

Fluxes of dissolved inorganic carbon in three estuarine systems of the Cantabrian Sea (north of Spain)

The diffusive and in situ fluxes of dissolved inorganic carbon (DIC) and total alkalinity (TA) have been measured and an estimation has been made of the water–atmosphere fluxes of CO₂ in three estuarine systems of the Cantabrian Sea during the spring of 1998. Each of these systems undergoes a different anthropogenic influence. The diffusive fluxes of dissolved inorganic carbon and total alkalinity obtained present values ranging between 0.54–2.65 and 0.0–2.4 mmol m⁻² day⁻¹, respectively. These ranges are in agreement with those of other coastal systems. The in situ fluxes are high and extremely variable (35–284 mmol TA m⁻² day⁻¹, 43–554 mmol DIC m⁻² day⁻¹ and 22–261 mmol dissolved oxygen (DO) m⁻² day⁻¹), because the systems studied are very heterogeneous. The values of the ratio of the in situ fluxes of TA and DIC show on average that the rate of dissolution of CaCO₃ is 0.37 times that of organic carbon oxidation. Equally, the interval of variation of the relationship between the benthic fluxes of inorganic carbon and oxygen (F_DIC/F_DO) is very wide (0.3–13.9), which demonstrates the different contributions made by the processes of aerobic and anaerobic degradation of the organic matter, as well as by the dissolution–precipitation of CaCO₃. The water–atmosphere fluxes of CO₂ present a clear dependence on the salinity. The brackish water of these systems (salinity<20), where maximum fluxes of 989 mmol m⁻² day⁻¹ have been estimated, act as a source of CO₂ to the atmosphere. The more saline zones of the estuary (salinity>30) act as a sink of CO₂, with fluxes between −5 and −10 mmol m⁻² day⁻¹.