Optimal interpolation of sea surface temperature for the North Sea and Baltic Sea

Sea surface temperature fields of the North Sea and Baltic Sea have been constructed for the year 2001 using a multiplatform Optimal Interpolation scheme. The analyzed fields are constructed every 12 h on a 10 km spatial grid. The product is based upon observations from the three NOAA satellites 12, 14 and 16 together with a large amount of in situ observations. Space dependent covariance functions are estimated from the satellite observations and account for spatial and temporal lags. Several independent methods have been used to assess the error on the sea surface temperature product. Compared against independent in situ observations, the mean RMS difference for the year 2001 is 0.78 °C. The spatial distribution of the errors reveals that the Baltic Sea in general show higher errors than the North Sea. The error statistics throughout the year show a temporal variation of the errors with maximum during summer and winter. Tests with a varying number of satellite observations show that the accuracy of the satellite observations is the most important parameter in terms of reducing the errors on the interpolated sea surface temperature product.     

Validation of a hybrid optimal interpolation and Kalman filter scheme for sea surface temperature assimilation

A hybrid data assimilation scheme designed for operational assimilation of satellite sea surface temperatures (SST) into an ocean model has been developed and validated against in-situ observations. The scheme consists of an optimal interpolation (OI) part and a greatly simplified Kalman filter (KF) part.The OI is performed only in the longitudinal and latitudinal directions. A climatological field is used as a background field for the interpolation. It is constructed by fitting daily averages of satellite SST to the annual mean, annual, and semiannual harmonics in a 20 km by 20 km grid. The background error covariance is approximated by a spatially varying two-dimensional exponential covariance model. The parameters of the covariance model are fitted to the deviations of the satellite data from the background field using data from a full year.The simplified KF uses ocean model forecasts as a background field. It is based on the assumption that it is possible to neglect horizontal SST covariances in the filter and that the typical time scale for vertical mixing in the mixed layer is much shorter than the average time between observations. We therefore assume that the error variance in a column of water is evenly spread out throughout the mixed layer. The result of these simplifications is a computationally very efficient KF.A one year validation of the scheme is performed for year 2001 using an operational eddy resolving ocean model covering the North Sea and the Baltic Sea. It is found that assimilation of sea surface temperature data reduces the model root mean square error from 1.13 °C to 0.70 °C. The hybrid scheme is found to reduce the root mean square error slightly more than the simplified KF without OI to 0.66 °C. The inclusion of spatially varying satellite error variances does not improve the performance of the scheme significantly.     

1/32° real-time global ocean prediction and value-added over 1/16° resolution

A 1/32° global ocean nowcast/forecast system has been developed by the Naval Research Laboratory at the Stennis Space Center. It started running at the Naval Oceanographic Office in near real-time on 1 Nov. 2003 and has been running daily in real-time since 1 Mar. 2005. It became an operational system on 6 March 2006, replacing the existing 1/16° system which ceased operation on 12 March 2006. Both systems use the NRL Layered Ocean Model (NLOM) with assimilation of sea surface height from satellite altimeters and sea surface temperature from multi-channel satellite infrared radiometers. Real-time and archived results are available online at http://www.ocean.nrlssc.navy.mil/global_nlom. The 1/32° system has improvements over the earlier system that can be grouped into two categories: (1) better resolution and representation of dynamical processes and (2) design modifications. The design modifications are the result of accrued knowledge since the development of the earlier 1/16° system. The improved horizontal resolution of the 1/32° system has significant dynamical benefits which increase the ability of the model to accurately nowcast and skillfully forecast. At the finer resolution, current pathways and their transports become more accurate, the sea surface height (SSH) variability increases and becomes more realistic and even the global ocean circulation experiences some changes (including inter-basin exchange). These improvements make the 1/32° system a better dynamical interpolator of assimilated satellite altimeter track data, using a one-day model forecast as the first guess. The result is quantitatively more accurate nowcasts, as is illustrated by several model-data comparisons. Based on comparisons with ocean color imagery in the northwestern Arabian Sea and the Gulf of Oman, the 1/32° system has even demonstrated the ability to map small eddies, 25–75 km in diameter, with 70% reliability and a median eddy center location error of 22.5 km, a surprising and unanticipated result from assimilation of altimeter track data. For all of the eddies (50% small eddies), the reliability was 80% and the median eddy center location error was 29 km. The 1/32° system also exhibits improved forecast skill in relation to the 1/16° system. This is due to (a) a more accurate initial condition for the forecast and (b) better resolution and representation of critical dynamical processes (such as upper ocean – topographic coupling via mesoscale flow instabilities) which allow the model to more accurately evolve these features in time while running in forecast mode (forecast atmospheric forcing for the first 5 days, then gradually reverting toward climatology for the remainder of the 30-day forecast period). At 1/32° resolution, forecast SSH generally compares better with unassimilated observations and the anomaly correlation of the forecast SSH exceeds that from persistence by a larger amount than found in the 1/16° system.     

Data assimilation into a Princeton Ocean Model of the Mediterranean Sea using advanced Kalman filters

This study investigates the effectiveness of the Singular Evolutive Extended Kalman filter (SEEK) and its variants (SEIK and SFEK filters) for data assimilation into a Princeton Ocean Model (POM) of the Mediterranean Sea. The SEEK filters are sub-optimal Kalman filters based on the approximation of the filter's error covariance matrices by singular low-rank matrices, reducing in this way extensive computational burden. At the initialization, the filters error covariance matrix is parameterized by a set of multivariate empirical orthogonal functions (EOFs) which describe the dominant modes of the system's variability. The Mediterranean model is implemented on a 1/4° × 1/4° horizontal grid with 25 sigma levels and is forced with 6-hour ECMWF re-analysis atmospheric data. Several twin experiments, in which pseudo-observations of altimetric data and/or data profiles were assimilated, were first performed to evaluate the filters performances and to study their sensitivities to different parameters and setups. The results of these experiments were very encouraging and helped in setting up an effective configuration for the assimilation of real data in near-real time situation. In the hindcast experiments, Topex/Poseidon and ERS weekly sea level anomaly data were first assimilated during 1993 and the filters solution was evaluated against independent Reynolds sea surface temperature (SST) analysis. The assimilation system was able to significantly enhance the consistency between the model and the assimilated data, although the improvement with respect to independent SST data was significantly less pronounced. The model SST was only improved after including SST data in the assimilation system.     

Comparison between a reanalyzed product by 3-dimensional variational assimilation technique and observations in the Ulleung Basin of the East/Japan Sea

Reanalyzed products from a MOM3-based East Sea Regional Ocean Model with a 3-dimentional variational data assimilation module (DA-ESROM), have been compared with the observed hydrographic and current datasets in the Ulleung Basin (UB) of the East/Japan Sea (EJS). Satellite-borne sea surface temperature and sea surface height data, and in-situ temperature profiles have been assimilated into the DA-ESROM. The performance of the DA-ESROM appears to be efficient enough to be used in an operational ocean forecast system.Comparing with the results from Mitchell et al. [Mitchell, D. A., Watts, D. R., Wimbush, M., Teague, W.J., Tracey, K. L., Book, J. W., Chang, K.-I., Suk, M.-S., Yoon, J.-H., 2005a. Upper circulation patterns in the Ulleung Basin. Deep-Sea Res. II, 52, 1617-1638.], the DA-ESROM fairly well simulates the high variability of the Ulleung Warm Eddy and Dok Cold Eddy as well as the branching of the Tsushima Warm Current in the UB. The overall root-mean-square error between 100 m temperature field reproduced by the DA-ESROM and the observed 100-dbar temperature field is 2.1 °C, and the spatially averaged grid-to-grid correlation between the two temperature fields is high with a mean value of 0.79 for the inter-comparison period.The DA-ESROM reproduces the development of strong southward North Korean Cold Current (NKCC) in summer consistent with the observational results, which is thought to be an improvement of the previous numerical models in the EJS. The reanalyzed products show that the NKCC is about 35 km wide, and flows southward along the Korean coast from spring to summer with maximum monthly mean volume transport of about 0.8 Sv in August–September.     

Seasonal variability in sea surface oceanographic conditions in the Aegean Sea (Eastern Mediterranean): an overview

Seasonal variability and the spatial distribution of sea surface temperatures (SST) and salinities (SSS) are reviewed, in relation to the prevailing climatological conditions, heat fluxes, water budget and general water circulation patterns. Within this context, consideration is given to: sea surface temperatures; air temperatures; precipitation; evaporation; wind speeds and directions; freshwater (mainly riverine) discharges throughout the Aegean; and the exchange of water masses with the Black Sea and eastern Mediterranean Sea. The investigation of satellite images, covering a 6-yr period (1988–1994), has enabled a synthesis of the monthly sea surface thermal distribution to be established.The climate of the Aegean Sea is characterised by annual air temperatures of 16–19.5°C, precipitation of about 500 mm yr⁻¹ and evaporation of some 4 mm d⁻¹. The Aegean has a negative heat budget (approximately −25 W m⁻²) and positive water balance (+ 1.0 m yr⁻¹), when inflow from the Black Sea is considered. During the summer, the (northerly) Etesians are the dominant winds over the Sea.Mean monthly sea surface temperatures (SST) vary from 8°C in the north during winter, up to 26°C in the south during summer. SST depends mainly upon air temperature; there is a month's delay between the former and latter maxima. The sea surface salinity (SSS) varies also spatially and seasonally, ranging from less than 31 psu, in the north, to more than 39 psu, in the southeast; lower values (< 25 psu) occur adjacent to the river mouths. SSSs present their maximum differences during summer, whilst during winter and autumn the distribution of SSS is more uniform. The overall spatial SST and SSS distribution pattern is controlled by: distribution of the (colder) Black Sea Waters; advection of the (warmer) Levantine Waters, from the southeastern part of the Aegean; upwelling and downwelling; and, to a lesser extent, but locally important, freshwater riverine inflows.     

Remotely-sensed chlorophyll a observations of the northern Red Sea indicate seasonal variability and influence of coastal reefs

The biological dynamics of the open northern Red Sea (21.5°–27.5° N, 33.5°–40° E) have not been studied extensively, due in part to both the inaccessibility of this desert region and political considerations. Remotely-sensed chlorophyll a data therefore provide a framework to investigate the primary patterns of biological activity in this oceanic basin. Monthly chlorophyll a data from the 8-year Sea-viewing Wide Field-of-View sensor (SeaWiFS) mission, and data from the Moderate Resolution Imaging Spectroradiometer (MODIS), were analyzed with the Goddard Earth Sciences Data and Information Services Center (GES DISC) online data analysis system “Giovanni”. The data indicate that despite the normal low chlorophyll a concentrations (0.1–0.2 mg m^− 3) in these oligotrophic waters, there is a characteristic seasonal bloom in March–April in the northernmost open Red Sea (24° to 27.5° N) concurrent with minimum sea surface temperature. The location of the highest chlorophyll concentrations is consistent with a linear box model [Eshel, G., and Naik, N.H., 1997. Climatological coastal jet collision, intermediate water formation, and the general circulation of the Red Sea. J. Phys. Oceanogr. 27(7), 1233–1257.] of Red Sea circulation. Two years in the data set exhibited a different seasonal cycle consisting of a relatively weak northern spring bloom and elevated chlorophyll concentrations to the south (21.5° to 24° N).The data also indicate that large coral reef complexes may be sources of either nutrients or chlorophyll-rich detritus and sediment, enhancing chlorophyll a concentration in waters adjacent to the reefs.     

Changes in the shelf macrobenthic community over large temporal and spatial scales in the Bohai Sea, China

Over the past 20 years, the Bohai Sea has been subjected to a considerable human impact through over-fishing and pollution. Together with the influence of the Yellow River cut-off, the ecosystem experienced a dramatic change. In order to integrate available information to detect any change in macrobenthic community structure and diversity over space and time, data collected during the 1980s and the 1990s from 3 regions of the Bohai Sea (Laizhou Bay, 16 stations, 37–38°N, 119–120.5°E; central Bohai Sea, 25 stations, 38–39°N, 119–121°E; eastern Bohai Bay, 12 stations, 38–39°N, 118.5–119°E) were reanalyzed in a comparative way by means of a variety of statistical techniques. A considerable change in community structure between the 1980s and the 1990s and over the geographical regions at both the species and family level were revealed. After 10 years, there was a considerable increase in abundance of small polychaetes, bivalves and crustaceans but decreased number of echinoderms. Once abundant in Laizhou Bay in the 1980s, a large echinoderm Echinocardium cordatum and a small mussel Musculista senhousia almost disappeared from the surveying area in the 1990s. Coupled with the increased abundance was the increased species richness in general whereas evenness was getting lower in central Bohai Sea and Bohai Bay but increased in Laizhou Bay. K-dominance plot showed the same trend as evenness J′. After 10 years, the macrobenthic diversity in the Bohai Sea as a whole was slightly reduced and a diversity ranking of central Bohai Sea > Laizhou Bay > eastern Bohai Bay over space was also suggested. Sediment granulometry and organic content were the two major agents behind the observed changes.     

Preoperational ocean forecasting in the southeastern Mediterranean Sea: Implementation and evaluation of the models and selection of the atmospheric forcing

Within the framework of several local and international programs, a quasi-operational ocean-forecasting system for the Southeastern Mediterranean Sea has been established and evaluated through a series of preoperational tests. The Princeton Ocean Model (POM) is used for simulating and predicting the hydrodynamics while the Wave Model (WAM) is used for predicting surface waves. Both models were set up to allow varying resolution and multiple nesting. In addition, POM was set up to be easily relocatable to allow rapid deployment of the model for any region of interest within the Mediterranean Sea. A common requirement for both models is the need for atmospheric forcing. Both models require time varying wind or wind stress. In addition, the hydrodynamic model requires initial conditions as well as time dependent surface heat fluxes, fresh water flux, and lateral boundary conditions at the open boundaries. Several sources of atmospheric forcing have been assessed based on their availability and their impact on the quality of the ocean models' forecasts. The various sources include operational forecast centers, other research centers, as well as running an in-house regional atmospheric model. For surface waves, higher spatial and temporal resolution of the winds plays a central role in improving the forecasts in terms of significant wave height and the timing of various high wave events. For the hydrodynamics, using the predicted wind stress and heat fluxes directly from an atmospheric model can potentially produce short range ocean forecasts that are nearly as good as hindcasts forced with gridded atmospheric analyses. Finally, a high-resolution, nested version of the model has shown to be stable under a variety of forcing conditions and time scales, thus indicating the robustness of the selected nesting strategy. For the southeastern corner of the Mediterranean, at forecast lead times of up to 4 days the high-resolution model shows improved skill over the coarser resolution driving model when compared to satellite derived sea surface temperatures. Most of the error appears to be due to the analysis error inherent in the initial conditions.     

Recovery of North-East Atlantic temperature fields from profiling floats: Determination of the optimal float number from sampling and instrumental error analysis

Argo is an international project that is deploying an array of temperature and salinity profiling floats over the global ocean. Here we use the error formulation derived from Optimal Statistical Interpolation to estimate statistical errors associated with the recovery of the temperature field in the North-East Atlantic ocean. Results indicate that with the present distribution of floats (119 in the considered domain), scales of wavelength larger than 500 km can be recovered with a relative uncertainty (rms error relative to the standard deviation of the field) of about 7% at 50 m, 8% at 200 m and 10% at 1000 m. This corresponds to mean absolute errors of 0.111 °C at 50 m, 0.104 °C at 200 m and 0.073 °C at 1000 m.The splitting of total errors into instrumental and sampling contributions reveals that, in the present scenario, errors are more due to the small number of floats than to instrumental errors, especially at upper levels. For scales larger than 500 km this will hold true until 200–250 floats are deployed (less than 200 for deep levels). In such a simulated scenario, the number of observations and the technology become approximately equally limiting factors for the accuracy of the temperature field mapping, with total relative errors of less than 2% at upper levels and about 3% at 1000 m.     

Forecast verification of a 3D model of the Mediterranean Sea. The use of discrete wavelet transforms and EOFs in the skill assessment of spatial forecasts

The quality assessment of a nested model system of the Mediterranean Sea is realised. The model has two zooms in the Provençal Basin and in the Ligurian Sea, realised with a two-way nesting approach. The experiment lasts for nine weeks, and at each week sea surface temperature (SST) and sea level anomaly are assimilated. The quality assessment of the surface temperature is done in a spatio-temporal approach, to take into account the high complexity of the SST distribution. We focus on the multi-scale nature of oceanic processes using two powerful tools for spatio-temporal analysis, wavelets and Empirical Orthogonal Functions (EOFs). We apply two-dimensional wavelets to decompose the high-resolution model and observed SST into different spatial scales. The Ligurian Sea model results are compared to observations at each of those spatial scales, with special attention on how the assimilation affects the model behaviour. We also use EOFs to assess the similarities between the Mediterranean Sea model and the observed SST. The results show that the assimilation mainly affects the model large-scale features, whereas the small scales show little or no improvement and sometimes, even a decrease in their skill. The multiresolution analysis reveals the connection between large- and small-scale errors, and how the choice of the maximum correlation length of the assimilation scheme affects the distribution of the model error among the different spatial scales.     

Air–sea CO2 fluxes in the Caribbean Sea from 2002–2004

Air–sea fluxes in the Caribbean Sea are presented based on measurements of partial pressure of CO₂ in surface seawater, pCO_2sw, from an automated system onboard the cruise ship Explorer of the Seas for 2002 through 2004. The pCO_2sw values are used to develop algorithms of pCO_2sw based on sea surface temperature (SST) and position. The algorithms are applied to assimilated SST data and remotely sensed winds on a 1° by 1° grid to estimate the fluxes on weekly timescales in the region. The positive relationship between pCO_2sw and SST is lower than the isochemical trend suggesting counteracting effects from biological processes. The relationship varies systematically with location with a stronger dependence further south. Furthermore, the southern area shows significantly lower pCO_2sw in the fall compared to the spring at the same SST, which is attributed to differences in salinity. The annual algorithms for the entire region show a slight trend between 2002 and 2004 suggesting an increase of pCO_2sw over time. This is in accord with the increasing pCO_2sw due the invasion of anthropogenic CO₂. The annual fluxes of CO₂ yield a net invasion of CO₂ to the ocean that ranges from − 0.04 to − 1.2 mol m^− 2 year^− 1 over the 3 years. There is a seasonal reversal in the direction of the flux with CO₂ entering into the ocean during the winter and an evasion during the summer. Year-to-year differences in flux are primarily caused by temperature anomalies in the late winter and spring period resulting in changes in invasion during these seasons. An analysis of pCO_2sw before and after hurricane Frances (September 4–6, 2004), and wind records during the storm suggest a large local enhancement of the flux but minimal influence on annual fluxes in the region.     

Recent observations in the straits of the East/Japan Sea: A review of hydrography, currents and volume transports

Recent observations of hydrography, currents and volume transports in the straits of the East/Japan Sea are reviewed. It is newly found that bottom cold water in the Korea/Tsushima Strait originating from the northern region of the East/Japan Sea appears not only in summer and autumn but also in winter. Intensive observations in the Korea/Tsushima Strait revealed two distinct cores of northeastward currents in the upper layer of the western and eastern channels. Mean volume transport through the Korea/Tsushima Strait is calculated as 2.5 ± 0.5 Sv from four-year direct and indirect measurements. As continuous monitoring has started in the Tsugaru and Soya Straits, understanding of temporal variability of currents and volume transports through the straits is in progress. For the first time, simultaneous time series of volume transports are available in the Korea/Tsushima and Tsugaru Straits during the winter of 1999–2000. Ouflow through the Tsugaru Strait accounts for about 70% of inflow through the Korea/Tsushima Strait for this period.     

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.     

Interannual variation of the Polar Front in the Japan/East Sea from summertime hydrography and sea level data

The Polar Front in the Japan/East Sea separates the southern warm water region from the northern cold water region. A merged TOPEX/POSEIDON and ERS-1/2 altimeter dataset and upper water temperature data were used to determine the frontal location and to examine the structure of its interannual variability from 1993 to 2001. The identified frontal location, where sea surface height gradient has a maximum about 10–20 cm over the horizontal distance of 100 km, corresponds well to the maximum subsurface horizontal temperature gradient. The front migrates more widely (36°N–41°N) in the western part of the sea than in the eastern part. The interannual migration induces large variability in upper water temperatures and sea surface height in the western region. Responsible physical mechanisms were studied using a reduced-gravity model. Differences between inflow and outflow change the total volume of warm water, and total warm water volume change in the warm water region uniformly pushes the front in the meridional direction across its mean position in the model simulation. Interannual variation of wind stress causes relatively wide migration of the modeled front in the western part.     

A study of the Alboran sea mesoscale system by means of empirical orthogonal function decomposition of satellite data

This paper presents the results of a combined empirical orthogonal function (EOF) analysis of Advanced Very High Resolution Radiometer (AVHRR) sea surface temperature (SST) data and sea-viewing wide field-of-view sensor (SeaWiFS) chlorophyll concentration data over the Alboran Sea (Western Mediterranean), covering a period of 1 year (November 1997–October 1998). The aim of this study is to go beyond the limited temporal extent of available in situ measurements by inferring the temporal and spatial variability of the Alboran Gyre system from long temporal series of satellite observations, in order to gain insight on the interactions between the circulation and the biological activity in the system. In this context, EOF decomposition permits concise and synoptic representation of the effects of physical and biological phenomena traced by SST and chlorophyll concentration. Thus, it is possible to focus the analysis on the most significant phenomena and to understand better the complex interactions between physics and biology at the mesoscale. The results of the EOF analysis of AVHRR-SST and SeaWiFS-chlorophyll concentration data are presented and discussed in detail. These improve and complement the knowledge acquired during the in situ observational campaigns of the MAST-III Observations and Modelling of Eddy scale Geostrophic and Ageostrophic motion (OMEGA) Project.     

About the seasonal variability of the Alboran Sea circulation

Data from a mooring line deployed midway between the Alboran Island and Cape Tres Forcas are used to study the time variability of the Alboran Sea from May 1997 to May 1998. The upper layer salinity and zonal velocity present annual and semiannual cycles characterised by a minimum in spring and autumn and a maximum in summer and winter. Temperature has the opposite behaviour to that of salinity indicating changes in the presence of the Atlantic water within the Alboran Passage. A large set of SST images is used to study these cycles. The decrease of salinity and velocity in our mooring location in spring and autumn seems to be related to the eastward drifting of the Western Alboran Gyre (WAG). The increase of salinity and velocity is caused by the Atlantic current flowing south of the Alboran Island and its associated thermohaline front. Conductivity–temperature–depth (CTD) data from two cruises along the 3°W are coherent with current meters and SST interpretations.During the period analysed, summer months are characterised by the stability of the two-gyre system, while in winter, the circulation is characterised by a coastal jet flowing close to the African shore. We use sea level differences across the Strait of Gibraltar for studying the variability of the Atlantic inflow. We discuss the changes in the Alboran Sea circulation and its relation with the variability of the inertial radius of the Atlantic inflow. Though our results are speculative, we find a possible relation between the disappearance of the two-gyre system and a reversal of the circulation in Gibraltar. Longer time series are needed to conclude, but comparison with previous works makes us think that the seasonal cycle described from May 1997 to May 1998 could be the most likely one for the Alboran Sea upper layer.     

Multigrid state vector for data assimilation in a two-way nested model of the Ligurian Sea

A system of two nested models composed by a coarse resolution model of the Mediterranean Sea, an intermediate resolution model of the Provençal Basin and a high resolution model of the Ligurian Sea is coupled with a Kalman-filter based assimilation method. The state vector for the data assimilation is composed by the temperature, salinity and elevation of the three models. The forecast error is estimated by an ensemble run of 200 members by perturbing initial condition and atmospheric forcings. The 50 dominant empirical orthogonal functions (EOF) are taken as the error covariance of the model forecast. This error covariance is assumed to be constant in time. Sea surface temperature (SST) and sea surface height (SSH) are assimilated in this system.     

The annual cycle of nanoflagellates in the Central Cantabrian Sea (Bay of Biscay)

The annual cycle of nanoflagellates (NF) including autotrophic (ANF), heterotrophic (HNF) and mixotropic (MNF) flagellates carried out in a temperate sea (Central Cantabrian Sea, southern Bay of Biscay) is presented. Three stations with characteristics ranging from coastal to oceanic conditions were analysed in order to compare NF response to this gradient. Samples were monthly collected at each station at three different depths between February 2002 and December 2002. CTD profiles were also taken at each station. NF were grouped according to their trophic status into ANF, HNF and MNF. Abundance and biomass were determined for each group. The annual cycle showed a general pattern consisting in a maximum in July with secondary maxima in March and October and minimum values in May. ANF were the most important fraction, making a major contribution (nearly 75%) to total NF biomass in all stations. HNF represented over 20% along the cycle, except for a peak in spring found in every station. MNF reached less than 5%, showing low seasonability. Small flagellates (2–5 µm) dominated throughout the cycle. Microplankton community was also analysed in terms of abundance and biomass. A significant positive correlation (r² = 0.49) was obtained between 2–5 µm NF and 10–20 µm HNF–MNF biomasses, suggesting a possible trophic relationship between these groups which should be cautiously taken. No significant relationships were found between microplankton and NF or between nutrients and ANF, indicating that the regulation of NF numbers is complex and probably implicates other groups. In addition to this, the unexpected 2002 Chl a concentration pattern and the misplacing of upwelling events render necessary to perform additional studies to fully understand the precise behaviour of NF in the Cantabrian Sea. To the best of our knowledge, this is the first study of a NF cycle in a temperate sea that considers all functional groups.     

Th-derived particulate organic carbon fluxes in the northern Barents Sea with comparison to drifting sediment trap fluxes

Large-volume sampling of ²³⁴Th was conducted to estimate particulate organic carbon (POC) export in conjunction with drifting sediment trap deployments in the northern Barents Sea in July 2003 and May 2005. ²³⁴Th-derived POC fluxes averaged 42.3 ± 39.7 mmol C m^− 2 d^− 1 in 2003 and 47.1 ± 30.6 mmol C m^− 2 d^− 1 in 2005. Sediment trap POC fluxes averaged 13.1 ± 8.2 mmol C m^− 2 d^− 1 in 2003 and 17.3 ± 11.4 mmol C m^− 2 d^− 1 in 2005, but better reflected the transient bloom conditions that were observed at each station within a season. Although ²³⁴Th fluxes agreed within a factor 2 at most stations and depths sampled, sediment trap POC fluxes were lower than large-volume POC flux estimates at almost every station. This may represent an under-collection of POC by the drifting sediment traps or, conversely, an over-collection of POC by the large-volume sampling of ²³⁴Th. It is hypothesized that the offset between the two methods is partly due to the presence of the prymnesiophyte Phaeocystis pouchetii, which potentially causes a large variation in > 53-μm POC/²³⁴Th ratios. Due to the large proportion of dissolved carbon or mucilage released by P. pouchetii, and because it is thought that P. pouchetii does not contribute significantly to the vertical export of biogenic matter in the Barents Sea, the application of large-volume sampling of ²³⁴Th may yield relatively high, and possibly inaccurate POC/²³⁴Th ratios. Hence, POC fluxes derived from ²³⁴Th sampling may be inappropriate and drifting sediment traps might be a more reliable method to measure the vertical export of biogenic matter in regions that have recurrent P. pouchetii blooms, such as the Barents Sea.