The importance of oceanographic fronts to marine birds and mammals of the southern oceans

During the last 30 years, at-sea studies of seabirds and marine mammals in the oceans south of the Subtropical Front have described an association with major frontal areas. More recently, the advancement in microtechnology has allowed the tracking of individuals and investigations into how these marine predators actually use the frontal zones. In this review, we examine 1) the relative importance to apex predators of the different frontal zones in terms of spatial distribution and carbon flux; 2) the processes that determine their preferential use; and 3) how the mesoscale dynamics of frontal structures drive at-sea foraging strategies of these predators. We review published results from southern waters and place them in a broader context with respect to what has been learned about the importance of fronts in oceans farther north.Some fronts constitute important boundaries for seabird communities in southern waters. At a mesoscale the maximum values of seabird diversity and abundance correspond to the location of the main fronts. At-sea surveys show a strong curvilinear correlation between seabird abundance and sea surface temperatures. High mean species richness and diversity for whales and seabirds are consistently associated with the southern water mass boundary of the Antarctic Circumpolar Current, the Subtropical Front and the Subantarctic Front; in the case of the Polar Front mean seabird densities are more variable. At small-scales, variation in seabird occurrence has been directly related to the processes at fronts in a limited number of cases. A significant positive relation was found between some plankton feeding species and frontal temperature gradient–phytoplankton variables.Telemetric studies have revealed that several apex predators (penguins, albatrosses, seals) perform long, directed foraging trips either to the Subtropical front or Polar Front, depending on locality. Seabirds with low flight costs, such as albatrosses, are able to reach fronts at long distances from colonies, showing variable foraging strategies as a function of the distances involved. Diving birds such as King penguins, that travel at a higher cost and lower speed, rely on the predictable spatial distribution of mesopelagic fish found close to the Polar Front. They may use the currents associated with eddies as oceanographic cues in the active search for frontal zones. Once in these areas they dive preferentially in and below the depth of the thermocline where catches per unit effort are high. Elephant seals concentrate foraging activity principally inside or at the boundary of cyclonic eddies. These mesoscale features appear to offer exceptional productivity favourable for foraging by various diving top predators.The connection between biophysical parameters at fronts and predators is likely to be made through biological enhancement. Top predators appear to forage at locations where prey are advected by physical processes and others where prey are produced locally. Long-term research on at-sea distributions and demographic parameters of top predators are essential to assess the consequences of potential shift in front distributions in relation to global warming. Such environmental changes would add to the impact of fish extraction by the industrial fisheries on the southern food webs.     

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.     

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.     

High-resolution surveys of the biogeochemistry of the New England shelfbreak front during Summer, 2002

We present the observations from a pair of field experiments at the New England shelfbreak front in June and August of 2002, each consisting of 14 cross-frontal surveys using the Lamont Pumping SeaSoar. Measurements of the front's physical, chemical, and bio-optical characteristics were made at high spatial and temporal resolution. The front, based on water-column hydrographic distributions, was found within a few km of the 200 m isobath during both cruises. We present here composite sections, based on averages of individual sections shifted in space to a common frontal location, of the cross-frontal distributions of these properties as a measure of the mean state of the front in both June and August. The observations show the familiar temperature, salinity, and density distributions of the summertime front, dominated by surface thermal heating. Nutrient and bio-optical distributions show the combined effects of water-mass exchange and biological processes. T, S, silicate, and phosphate distributions are suggestive of cross-frontal exchange of slope- and shelf waters, although transport mechanisms and pathways are not apparent. These properties, along with nitrate and optical measures of the suspended particle distributions, show vertical displacements of isopleths as the front is approached: property contours slope upwards toward the front from the shoreward side, and downwards toward the front from the seaward side. Again, actual water-movement pathways are not constrained by these suggestive patterns. Bio-optical distributions show elevated indicators of photosynthetic efficiency both seaward and shoreward of the front, but the front itself is a minimum in biomass. Accumulation of photosynthetic biomass appears to be controlled primarily by nitrate scarcity in waters within and above the pycnocline. At the base of the pycnocline, light limitation appears to be the controlling factor, although the base of the euphotic zone is deeper than the biomass maxima and the base of the pycnocline. Mechanisms explaining this phenomenon are unclear, but tenuous evidence suggests low stratification at the depth of the 1% light level may not allow phytoplankton to optimize for the low-light, high-nutrient conditions at depth. Cross-frontal differences in nutrient and bio-optical parameters, particularly in August, suggest distinct phytoplankton assemblages, and the presence of calcite-forming or nitrogen fixing groups to the community structure, especially in very shallow waters across the front and in pycnocline waters seaward of the front.     

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.     

Interannual variability of the shelf-slope front position between 75° and 50° W

The shelf-slope front (SSF) is a continuous shelf-break front running from the Tail of the Grand Banks to Cape Hatteras, North Carolina, separating colder and less-saline continental shelf waters from warmer and more saline slope waters. Time series containing mean monthly SSF positions were produced along each of 26 longitude lines between 75° and 50°W by workers located at Bedford Institute of Oceanography by digitizing individual frontal charts and computing mean monthly latitudinal positions over a 29-year (1973–2001) period. After removing seasonal variability at each longitude, interannual variability (IAV) of the SSF position at each longitude was computed as the annual mean of all monthly SSF position anomalies for each year over the 29-year period. Despite some missing data, a longitude-time plot reveals alternating bands of offshore (late-1970s, late-1980s, late-1990s) and onshore (early-1980s, early-1990s, early-2000s) annual mean SSF anomaly values, exhibiting a period of approximately 10 years. Annual mean SSF anomaly amplitudes are largest in the east, with maxima of O (± 100 km) located east of 60° W for years when data are available. Empirical orthogonal function (EOF) modes 1–4 (accounting for > 90% of the variance) form a set of basis functions that describe the SSF anomaly data and allow reconstruction of the entire data set since missing data are relatively few (14%). A complex empirical orthogonal function (CEOF) analysis using the “reconstructed” data reveals a wavelength scale of approximately 20° of longitude, a distance nearly equal to the entire study domain, along with steady, westward phase propagation of SSF anomalies over approximately the same distance. Speed calculations for the westward-propagating features yield a value of approximately 1.2 to 2.4 cm s^− 1 (1 to 2 km d^− 1), with annual mean SSF anomalies thus requiring about 4 years to propagate from the Tail of the Grand Banks in the east to Cape Hatteras, North Carolina, in the west. This propagation speed and the timing of the SSF positional anomalies at the Tail of the Grand Banks for the 29-year study period are in agreement with speeds computed for the propagation of quasi-decadal salinity anomalies through the Labrador Sea and the time of their arrival at the Tail of the Grand Banks. The small westward SSF anomaly propagation speed is an order of magnitude smaller than the associated currents, in agreement with a highly damped flow-through system originating from both Davis Strait and the West Greenland Current as discussed by other workers. Observations from both southern and northern portions of the study domain, within both continental shelf and slope waters, show that interannual changes in the volume of shelf water along with shelf water bulk properties exhibit a strong relationship with IAV of the SSF position over long time periods.     

An algorithm for oceanic front detection in chlorophyll and SST satellite imagery

An algorithm is described for oceanic front detection in chlorophyll (Chl) and sea surface temperature (SST) satellite imagery. The algorithm is based on a gradient approach: the main novelty is a shape-preserving, scale-sensitive, contextual median filter applied selectively and iteratively until convergence. This filter has been developed specifically for Chl since these fields have spatial patterns such as chlorophyll enhancement at thermohaline fronts and small- and meso-scale chlorophyll blooms that are not present in SST fields. Linear Chl enhancements and localized (point-wise) blooms are modeled as ridges and peaks respectively, whereas conventional fronts in Chl and SST fields are modeled as steps or ramps. Examples are presented of the algorithm performance using modeled (synthetic) images as well as synoptic Chl and SST imagery. After testing, the algorithm was used on > 6000 synoptic images, 1999–2007, to produce climatologies of Chl and SST fronts off the U.S. Northeast.     

Surface patterns of zooplankton spatial variability detected by high frequency sampling in the NW Mediterranean. Role of density fronts

High frequency sampling was performed in daylight hours along a 35 km transect in the Ligurian Sea to investigate the upper layer zooplankton distribution during the spring phytoplankton bloom. The results show detailed spatial structure and biomass of key zooplankton functional groups, copepods, salps and krill larvae, within the different water masses characterizing this region. Although observed values of total copepod biomass distribution were rather constant along the transect, species-specific patterns were observed in the copepod spatial distribution. The larger species Calanus helgolandicus, as well as Centropages typicus, Oithona spp., and Oncaea spp., were associated with the frontal zone. However, Acartia spp. had a scattered distribution, and Clausocalanus/Paracalanus did not have a clear pattern. In addition, krill larvae were concentrated in the frontal area and salps had a scattered pattern. The cross-shore zooplankton distribution appeared strongly influenced by both the Northern Ligurian current governing inshore waters, which acts as a major flushing forcing, and the Ligurian front, which governs offshore waters and may act as retention area for zooplankton.     

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.     

A generalized model of pelagic biogeochemistry for the global ocean ecosystem. Part II: Numerical simulations

This paper presents a global ocean implementation of a multi-component model of marine pelagic biogeochemistry coupled on-line with an ocean general circulation model forced with climatological surface fields (PELAgic biogeochemistry for Global Ocean Simulations, PELAGOS). The final objective is the inclusion of this model as a component in an Earth System model for climate studies. The pelagic model is based on a functional stoichiometric representation of marine biogeochemical cycles and allows simulating the dynamics of C, N, P, Si, O and Fe taking into account the variation of their elemental ratios in the functional groups. The model also includes a parameterization of variable chlorophyll/carbon ratio in phytoplankton, carrying chl as a prognostic variable. The first part of the paper analyzes the contribution of non-local advective–diffusive terms and local vertical processes to the simulated chl distributions. The comparison of the three experiments shows that the mean chl distribution at higher latitudes is largely determined by mixing processes, while vertical advection controls the distribution in the equatorial upwelling regions. Horizontal advective and diffusive processes are necessary mechanisms for the shape of chl distribution in the sub-tropical Pacific. In the second part, the results have been compared with existing datasets of satellite-derived chlorophyll, surface nutrients, estimates of phytoplankton community composition and primary production data. The agreement is reasonable both in terms of the spatial distribution of annual means and of the seasonal variability in different dynamical oceanographic regions. Results indicate that some of the model biases in chl and surface nutrients distributions can be related to deficiencies in the simulation of physical processes such as advection and mixing. Other discrepancies are attributed to inadequate parameterizations of phytoplankton functional groups. The model has skill in reproducing the overall distribution of large and small phytoplankton but tends to underestimate diatoms in the northern higher latitudes and overestimate nanophytoplankton with respect to picoautotrophs in oligotrophic regions. The performance of the model is discussed in the context of its use in climate studies and an approach for improving the parameterization of functional groups in deterministic models is outlined.     

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.     

Spatial and temporal variation of bacterioplankton in a sub-Antarctic coastal area (Kerguelen Archipelago)

Bacterial abundance and production were measured monthly for one year along cross-shore transects in 3 sub-Antarctic fjords of the Kerguelen Archipelago (seven stations each). Mean values of the 3 most coastal (inside) and most offshore (outside) stations were used to describe the relationship between temperature, phytoplankton biomass, bacterial abundance and bacterial production over a one year annual cycle. The entire sampling protocol was repeated twice during each cruise: once at noon and once at midnight. Over the whole sampling period, the temperature ranged from 2.1 to 7.4 °C, while chlorophyll a concentrations varied by a factor of 10, and bacterial abundance and production varied by factors of 12 and 30, respectively. Within one day, all of these parameters sometimes varied by a factor of 4 between noon and midnight. A clear seasonality was observed for all of the parameters. However, while variations of phytoplankton and bacterial production paralleled those of temperature, bacterial abundance was low in midsummer and maximum in autumn. While no general pattern could be observed from the total data set, spatial gradients could interfere strongly with temporal changes.     

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.     

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.     

Long-term sea surface temperature baselines—time series, spatial covariation and implications for biological processes

Coastal areas such as estuaries, bays and fjords usually have hydrographic characteristics (e.g., temperature, salinity) which differ from those at larger spatial scales and in offshore areas. The differences can arise if the areas are subject to different climatic forcing or if they are relatively isolated from each other due to topographic and ocean circulation features which inhibit advective inputs of water mass properties. Local differences in hydrographic conditions can therefore potentially limit the applicability of existing long time series of coastally monitored temperatures for addressing questions at large spatial scales, such as the response of species distributions and phenologies to climate change. In this study we investigate the spatial synchrony of long-term sea surface temperatures in the North Sea–Baltic Sea region as measured daily at four coastal sites (Marsdiep, Netherlands; Torungen, Norway; Skagens Reef, Denmark; and Christiansø, Denmark) and in several large offshore areas. All time series, including two series reconstructed and intercalibrated for this study (Skagens Reef and Christiansø, Denmark), began during 1861–1880 and continue until at least 2001. Temperatures at coastal sites co-varied strongly with each other and with opportunistically measured offshore temperatures despite separation distances between measuring locations of 20–1200 km. This covariance is probably due to the influence of large-scale atmospheric processes on regional temperatures and is consistent with the known correlation radius of atmospheric fluctuations (ca. 1000 km). Differences (e. g, long-term trends, amplitude of seasonal variations) between coastal temperatures and those measured in adjacent offshore areas varied nonrandomly over time and were often significantly autocorrelated up to 2 years. These differences suggest that spatial variations in physical oceanographic phenomena and sampling heterogeneities associated with opportunistic sampling could affect perceptions of biological responses to temperature fluctuations. The documentation that the coastally measured temperatures co-vary with those measured opportunistically in offshore areas suggests that the coastal data, which have been measured daily using standardized methods and instruments, contain much of the variability seen at larger spatial scales. We conclude that both types of time series can facilitate assessments of how species and ecosystems have responded to past temperature changes and how they may react to future temperature changes.     

Surface chlorophyll in the Black Sea over 1978–1986 derived from satellite and in situ data

Absolute values of chlorophyll a concentration and its spatial and seasonal variations in the Black Sea were assessed by using satellite CZCS and in situ data. Since the satellite CZCS had operated for the 1978–1986 period, CZCS data was used for assessing the past state of the Black Sea just before the onset of drastic changes observed in late 1980s. The approach used for the calculation of the absolute values of chlorophyll a concentration from CZCS data was based on the direct comparison of in situ chlorophyll a data and those of CZCS and by applying the algorithm developed for the transformation of CZCS data into chlorophyll a values. CZCS Level 2 data related with pigment concentration having a spatial resolution of 1 km at nadir were used. The daily Level 3 files were derived by binning Level 2 values into 4-km grid cells and the monthly and seasonal Level 3 files were created by averaging the daily Level 3 files over the corresponding period. In situ chlorophyll a data were obtained by spectrophotometric and fluorometric methods in 15 scientific cruises over the 1978–1986 period. Total number of ship-measured data used for the comparison with those CZCS values was 590.Chlorophyll a concentration (Chl) was derived from CZCS values (C) with regression equations Chl=kC; the coefficient of transformation k was calculated from six different data sets by taking into account distinctions between subregions and seasons. The reasons for difference in the k values have been analyzed.Statistical comparison of the chlorophyll a values measured in situ and those derived from CZCS data was based on log-transformed data and gave the following results: regression SLOPE=0.842, regression INTERCEPT=−0.081, coefficient of determination (R²)=0.806, root–mean–square ERROR=0.195. The mean monthly chlorophyll a distributions derived from CZCS data over 1978–1986 have been constructed and the mean seasonal chlorophyll a values in different regions have been calculated and analyzed. The significant difference in chlorophyll concentration between the western shelf regions and the open part of the Black Sea has been demonstrated, especially in warm season. At almost all seasons, the highest chlorophyll concentration is observed in the western interior shelf region which is under strong influence of Danube. The summer mean chlorophyll concentration in this region is 18 times higher than that in the open parts and about nine times higher than in the eastern shelf region. The greatest seasonal variations are observed in the open part of the Black Sea: chlorophyll concentration in cold season is four to six times higher than in summer and three to five times higher than in April and October. To the contrary, in the western interior shelf regions, the concentration is higher in May–October (about twice than that in November–March). Seasonal variations in the western outer shelf regions are smoothed out as compared with both the western interior shelf and the open regions.     

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.     

An analysis of the characteristics of chlorophyll in the Sulu Sea

The seasonal and spatial variations of chlorophyll concentrations, Sea Surface Temperature (SST), wind fields and wind-induced Ekman pumping in the Sulu Sea are investigated using a set of new remote sensing measurements from October 1997 to December 2004. The results show the seasonality of chlorophyll, wind fields and SST and reveal the phytoplankton blooming events in the Sulu Sea basin during the northeast monsoon season. In summer, chlorophyll concentrations were relatively low (< 0.2 mg/m³) and distributed uniformly throughout the basin with a narrow belt of high chlorophyll concentrations along the coastal waters, particularly the coasts of Borneo and of the Sulu Archipelago. In winter, chlorophyll concentrations increased (> 0.2 mg/m³) throughout the entire basin, and phytoplankton bloomed southward to the central basin, while chlorophyll concentrations reached high levels (1 mg/m³) in the center of the blooms. One peak was observed during the northeast monsoon season each year. SSTs have significant negative correlations with chlorophyll concentrations; i.e., high and uniformly distributed in summer but lower with an obvious tongue of cold waters southward to the central basin in winter. The seasonal variation of chlorophyll concentrations and SST distribution were associated with the seasonally reversing monsoon. The winter phytoplankton blooming and the tongue of the cold waters were correlated to the vertical upwelling cold and nutrient-rich waters drawn by the northeast wind, with the center of the blooms and the location of cold tongues coinciding with the maximum of the wind speeds and the Ekman pumping velocities.     

Validation of the 3D biogeochemical model MIRO&CO with field nutrient and phytoplankton data and MERIS-derived surface chlorophyll a images

This paper presents results obtained with MIRO&CO-3D, a biogeochemical model dedicated to the study of eutrophication and applied to the Channel and Southern Bight of the North Sea (48.5°N–52.5°N). The model results from coupling of the COHERENS-3D hydrodynamic model and the biogeochemical model MIRO, which was previously calibrated in a multi-box implementation. MIRO&CO-3D is run to simulate the annual cycle of inorganic and organic carbon and nutrients (nitrogen, phosphorus and silica), phytoplankton (diatoms, nanoflagellates and Phaeocystis), bacteria and zooplankton (microzooplankton and copepods) with realistic forcing (meteorological conditions and river loads) for the period 1991–2003. Model validation is first shown by comparing time series of model concentrations of nutrients, chlorophyll a, diatom and Phaeocystis with in situ data from station 330 (51°26.00′N, 2°48.50′E) located in the centre of the Belgian coastal zone. This comparison shows the model's ability to represent the seasonal dynamics of nutrients and phytoplankton in Belgian waters. However the model fails to simulate correctly the dissolved silica cycle, especially during the beginning of spring, due to the late onset (in the model) of the early spring diatom bloom. As a general trend the chlorophyll a spring maximum is underestimated in simulations. A comparison between the seasonal average of surface winter nutrients and spring chlorophyll a concentrations simulated with in situ data for different stations is used to assess the accuracy of the simulated spatial distribution. At a seasonal scale, the spatial distribution of surface winter nutrients is in general well reproduced by the model with nevertheless a small overestimation for a few stations close to the Rhine/Meuse mouth and a tendency to underestimation in the coastal zone from Belgium to France. PO₄ was simulated best; silica was simulated with less success. Spring chlorophyll a concentration is in general underestimated by the model. The accuracy of the simulated phytoplankton spatial distribution is further evaluated by comparing simulated surface chlorophyll a with that derived from the satellite sensor MERIS for the year 2003. Reasonable agreement is found between simulated and satellite-derived regions of high chlorophyll a with nevertheless discrepancies close to the boundaries.     

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

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