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.     

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.     

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.     

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.     

The shoaling of nutrient-enriched subsurface waters as a mechanism to sustain primary productivity off Central Baja California during El Niño winters

Using CalCOFI data for coastal shallow stations (above 100 m depth), higher than expected nitrate concentrations were detected in near-surface high-temperature waters off of Central Baja California during some El Niño winters. Though recent data are not available for Central Baja California, past El Niño data, though limited, showed nitrate concentrations above 16 μM at temperatures above 16 °C, and nitrate concentrations between 1 and 2 μM at 19 °C, while the previously established relationship of temperature and nitrate for California Current waters predicts nitrate depletion above 14 or 15 °C. The anomalous, high temperature–high nitrate enrichment events documented in Central Baja California were detected as shallow as 9 m and as deep as 73 m, were associated with low-oxygen (between 2 and 4 ml/l) and high-salinity (between 33.8 and 34.3 psu) waters, and occurred during the winter months of an El Niño year. Using recent data for San Diego, CA, similar but weaker enrichment events were detected for the El Niño winter of 1997–1998. The periodic shoaling of a subsurface subtropical water mass of high temperature, high salinity, low oxygen and high nutrients during some El Niño winters is proposed to cause periodic enrichment and to maintain productivity during warming events in this area. Enrichment events were not detected off Ensenada, in Northern Baja California, possibly due to the amplification of the onshore flow during El Niño there, or due to the Ensenada front. The proposed mechanism of periodic enrichment of nutrient-depleted surface waters during some El Niño winters by subsurface waters from the California Undercurrent may explain the following: (1) survival of giant kelp forests at their southern limit in Central Baja California documented during past El Niño events in warm waters, (2) the rapid recovery and high carrying capacity of giant kelp documented after the mass disappearance during El Niño 1997–1998, and (3) the increase in the extent of mesotrophic chlorophyll detected in the area during the 1997–1998 and 1982–1983 El Niño events.     

Decadal-scale variations of trophic levels at high trophic levels in the Yellow Sea and the Bohai Sea ecosystem

A total of 2759 stomachs collected from a bottom trawl survey carried out by R/V “Bei Dou” in the Yellow Sea between 32°00 and 36°30N in autumn 2000 and spring 2001 were examined. The trophic levels (TL) of eight dominant fish species were calculated based on stomach contents, and trophic levels of 17 dominant species in the Yellow Sea and the Bohai Sea reported in later 1950s and mid-1980s were estimated so as to be comparable. The results indicated that the mean trophic level at high trophic levels declined from 4.06 in 1959–1960 to 3.41 in 1998–1999, or 0.16–0.19·decade^− 1 (mean 0.17·decade^− 1) in the Bohai Sea, and from 3.61 in 1985–1986 to 3.40 in 2000–2001, or 0.14·decade^− 1 in the Yellow Sea; all higher than global trend. The dominant species composition in the Yellow Sea and the Bohai Sea changed, with the percentage of planktivorous species increases and piscivorous or omnivorous species decreases, and this was one of the main reasons for the decline in mean trophic level at high tropic levels. Another main reason was intraspecific changes in TL. Similarly, many factors caused decline of trophic levels in the dominant fish species in the Yellow Sea and the Bohai Sea. Firstly, TL of the same prey got lower, and anchovy (Engraulis japonicus) as prey was most representative. Secondly, TLs of diet composition getting lower resulted in not only decline of trophic levels but also changed feeding habits of some species, such as spotted velvetfish (Erisphex pottii) and Trichiurus muticus in the Yellow Sea. Thirdly, species size getting smaller also resulted in not only decline of trophic levels but also changed feeding habits of some species, such as Bambay duck (Harpodon nehereus) and largehead hairtail (Trichiurus haumela). Furthermore, fishing pressure and climate change may be interfering to cause fishing down the food web in the China coastal ocean.     

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.     

Diatom stratigraphy and long-term dissolved silica concentrations in the Baltic Sea

In many parts of the world coastal waters with anthropogenic eutrophication have experienced a gradual depletion of dissolved silica (DSi) stocks. This could put pressure on spring bloom diatom populations, e.g. by limiting the intensity of blooms or by causing shifts in species composition. In addition, eutrophication driven enhanced diatom growth is responsible for the redistribution of DSi from the water phase to the sediments, and changes in the growth conditions may be reflected in the sediment diatom stratigraphy.To test for changes in diatom communities we have analyzed four sediment cores from the Baltic Sea covering approximately the last 100 years. The sediment cores originate from the western Gulf of Finland, the Kattegat, the Baltic Proper and the Gulf of Riga. Three out of the four cores reveal only minor changes in composition of diatom assemblages, while the Gulf of Riga core contains major changes, occurring after the second World War. This area is set apart from the other Baltic Sea basins by a high frequency of low after spring bloom DSi concentrations (< 2 µmol L^− 1) during a relatively well defined time period from 1991–1998. In 1991 to 1993 a rapid decline of DSi spring concentrations and winter stocks (down to 5 µmol L^− 1) in the Gulf was preceded by exceptionally intense diatom spring blooms dominated by the heavily silicified species Thalassiosira baltica (1991–1992; up to 5.5 mg ww L^− 1). T. baltica has been the principal spring bloom diatom in the Gulf of Riga since records began in 1975. DSi consumption and biomass yield experiments with cultured T. baltica suggest that intense blooms can potentially exhaust the DSi stock of the water column and exceed the annual Si dissolution in the Gulf of Riga. The phytoplankton time series reveals another exceptional T. baltica bloom period in 1981–1983 (up to 8 mg L^− 1), which, however, took place before the regular DSi measurements. These periods may be reflected in the conspicuous accumulation of T. baltica frustules in the sediment core corresponding to ca. 1975–1985.     

Large-scale physical controls on phytoplankton growth in the Irminger Sea, Part II: Model study of the physical and meteorological preconditioning

The upper water column in the Irminger Sea is characterized by cold fresh arctic and subarctic waters and warm saline North Atlantic waters. In this study the local physical and meteorological preconditioning of the phytoplankton development over an annual cycle in the upper water column in four physical zones of the Irminger Sea is investigated. Data from four cruises of the UK's Marine Productivity programme are combined with results from a coupled biological–physical nitrogen–phytoplankton–zooplankton–detritus model run using realistic forcing. The observations and model predictions are compared and analyzed to identify the key parameters and processes which determine the observed heterogeneity in biological production in the Irminger Sea. The simulations show differences in the onset of the bloom, in the time of the occurrence of the maximum phytoplankton biomass and in the length of the bloom between the zones. The longest phytoplankton bloom of 90 days duration was predicted for the East Greenland Current of Atlantic origin zone. In contrast, for the Central Irminger Sea zone a phytoplankton bloom with a start at the beginning of May and the shortest duration of only 70 days was simulated. The latest onset of the phytoplankton bloom in mid May and the latest occurrence of the maximum biomass (end of July) were predicted for the Northern Irminger Current zone. Here the bloom lasted for 80 days. In contrast the phytoplankton bloom in the Southern Irminger Current zone started at the same time as in Central Irminger Sea, but peaked end of June and lasted for 80 days. For all four zones relatively low daily (0.3–0.5 g C m^− 2d^− 1) and annual primary production was simulated, ranging between 35.6 g C m^− 2y^− 1 in the East Greenland Current of Atlantic origin zone and 45.6 g C m^− 2y^− 1 in the Northern Irminger Current zone. The model successfully simulated the observed regional and spatial differences in terms of the maximum depth of winter mixing, the onset of stratification and the development of the seasonal thermocline, and the differences in biological characteristics between the zones. The initial properties of the water column and the seasonal cycle of physical and meteorological forcing in each of the zones are responsible for the observed differences during the Marine Productivity cruises. The timing of the transition from mixing to stratification regime, and the different prevailing light levels in each zone are identified as the crucial processes/parameters for the understanding of the dynamics of the pelagic ecosystem in the Irminger Sea.     

Picophytoplankton and nanophytoplankton abundance and distribution in the southeastern Beaufort Sea (Mackenzie Shelf and Amundsen Gulf) during Fall 2002

The distribution of picophytoplankton (0.2–2 µm) and nanophytoplankton (2–20 µm) in the Beaufort Sea–Mackenzie Shelf and Amundsen Gulf regions during autumn, 2002 is examined relative to their ambient water mass properties (salinity, temperature and nutrients: nitrate + nitrite, phosphate, and silicate) and to the ratio of variable to maximum fluorescence, Fv/Fm. Total phytoplankton and cell abundances (< 20 µm) were mainly correlated with salinity. Significant differences in picophytoplankton cell numbers were found among waters near the mouth of the Mackenzie River, ice melt waters and the underlying halocline water masses of Pacific origin. Picophytoplankton was the most abundant phytoplankton fraction during the autumnal season, probably reflecting low nitrate concentrations (surface waters average ~ 0.65 µM). The ratio Fv/Fm averaged 0.44, indicating that cells were still physiologically active, even though their concentrations were low (max Chl a = 0.9 mg m^− 3). No significant differences in Fv/Fm were evident in the different water masses, indicating that rate limiting conditions for photosynthesis and growth were uniform across the whole system, which was in a pre-winter stage, and was probably already experiencing light limitation as a result of shortening day lengths.     

The Dead Sea hydrography from 1992 to 2000

The modern hydrological regime of the Dead Sea is strongly affected by anthropogenic activity. The natural fresh water budget has changed mainly due to the drastic reduction of runoff. Since 1977, the surface level of the Dead Sea has been lowered by an average rate of about 60 cm/year and for the period from 1998 to 2000, the lowering rate has reached about 100 cm/year. As a result of the runoff reduction, the upper layer salinity of the Dead Sea has increased and the gravitational stability of the water body was diminished. Eventually, during the winter of 1978–1979, the lake waters overturned, bringing to an end the long-term stable meromictic¹ hydrological regime. The lake entered a new phase in which its hydrological regime switches between holomictic and meromictic regimes, depending on the size of the runoff into the lake (i.e. the amount of precipitation in the lake's watershed). The first holomictic period, 1979–1980, lasted for 2 months only. It was succeeded by a 4-year meromictic period (1980–1983). The second holomictic period lasted for 9 years (1983–1991). The rainy winter of 1991–1992 resulted in an almost 2-m sea level rise. The upper layer with a relatively low salinity was restored and a new meromictic period persisted for 4 years, until winter 1995–1996. During the last meromictic period, the hydrological regime of the Dead Sea was characterized by following long-term trends: the depth of the summer thermocline increased from 12–15 to 25–30 m; the quasi-salinity of the upper layer, initially of about 164 kg/m³, increased rapidly at a rate of about 16–18 kg/m³/year; the quasi-salinity of the deep water, initially of about 235 kg/m³, decreased slowly at a rate of about 0.08–0.10 kg/m³/year (for the sake of comparison, a quasi salinity of 235 kg/m³ is the equivalent of 280‰ “usual” salinity); and the winter minimal temperature of the upper layer, initially of about 16 °C, increased rapidly at a rate of about 2 °C/year. In November 1995, the latest meromictic period of the Dead Sea came to an end. During the present holomictic period, 1996–2000, the hydrological regime of the Dead Sea is also characterized by long-term trends: the quasi-salinity of the entire Dead Sea increased at a rate of about 0.5 kg/m³/year, with practically no decrease during the winters; the temperature of the deep water mass increased with a rate of about 0.25 °C/year; and the period of vertical convection of the entire water column, initially about 3 months, increased at a rate of about 1 week/year. Moreover, we observed that the temperature and salinity of the bottom layer in the deepest part of the Dead Sea raised by about 0.5–0.6 °C and 0.15–0.25 kg/m³ during each holomictic summer.     

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

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

Dynamics of inorganic nutrient species in the Bohai seawaters

Within the frame of a Sino-German Joint Research Program, two cruises of “R/V Dong Fang Hong 2” were carried out in September–October 1998 and April–May 1999, respectively, to understand the dynamics of nutrients in the Bohai. Nutrient species (NO₃⁻, NO₂⁻, NH₄⁺, PO₄³⁻ and SiO₃²⁻) are determined colorimetrically on board for five anchor and 30 grid stations. In situ incubation experiments are performed to determine planktonic nutrient uptake and benthic exchange flux. Nutrient concentrations display short-term variability and seasonal change in the Bohai, with higher levels in shallow coastal waters than in the Central Bohai. The influence of riverine discharge on nutrient levels can be seen from salinity isopleths, nutrient distribution and species ratios. Near-bottom (nb) waters have similar nutrient concentrations as to the surface waters in the Central Bohai, whereas stratification takes place in the Bohai Strait and North Yellow Sea. In situ incubation experiments provide evidence that the uptake ratio (i.e. N, P) by phytoplankton is proportional to the ratios among nutrient species in ambient waters. Based on the data of this study and previously publications, a preliminary estimate of nutrient budgets via riverine input and atmospheric deposition is established. The results indicate that atmospheric deposition gains importance over rivers in delivering nutrients into the Bohai and sustain the new production, following recent decrease in riverine inflow caused by drought periods in North China and damming practices. A historical review of nutrient data indicates that concentrations of nitrogen increase and phosphorus and silica decrease in the Central Bohai over last 40 years. This potentially has an important influence on the health of ecosystem in Bohai (e.g. food web and community structure), though further study is needed to examine the scenario in more detail.     

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.     

The trophic state of the Baltic Sea a century ago: a model simulation study

We apply a 3-D circulation model with a biogeochemical module (ERGOM) for the simulation of trophic conditions in the Baltic Sea a century ago. One aim is to provide reference or background data for nitrogen, phosphorus and chlorophyll, which is required for the implementation of the European Water Framework Directive (WFD). We assume that the situation a century ago serves this purpose well. Model input for this long-term simulation study are the regionally differentiated riverine and atmospheric nutrient loads to the Baltic Sea, which were compiled and calculated for a situation 100 years ago on the basis of various literature sources. For the mixed surface layer of the open Baltic Sea, we suggest maximum winter concentrations for dissolved phosphorus (dissolved inorganic nitrogen) of 0.23–0.35 (2.7–3.7) mmol/m³. Maximum chlorophyll-a concentrations are between 1.8 and 2.4 mg/m³. The concentrations of all parameters for different coastal waters vary in a wide range, depending on exposure to nutrient sources. Our nutrient concentrations for the situation a century ago are close to early measured data (1950–1960) and suggest that this data is suitable as reference data, as well.     

Environmental changes and the responses of the ecosystems of the Yellow Sea during 1976–2000

The changes in the environmental features of the Yellow Sea during the last 25 years of the 20th century were studied using a set of seasonally monitored data along a transect (at 36°N) maintained by the State Oceanic Administration of China. The data included the ocean temperature (T), salinity (S) and biogenic elements, such as dissolved oxygen (DO), phosphorus (PO₄-P), silicon (SiO₃-Si) and dissolved inorganic nitrogen (DIN).The seasonal (summer and winter) values and the annual mean of these elements showed significant changes during the monitored period. Time series of T, S, DIN and N:P ratios exhibited positive trends, while those of DO, P and Si exhibited negative trends. During this period, the annual mean of T and DIN in the Yellow Sea increased by 1.7 °C and 2.95 μmol L⁻¹, respectively, while those of DO, P and Si decreased by 59.1, 0.1 and 3.93 μmol L⁻¹, respectively. In the 1980s, particularly in between 1985 and 1989, concentrations of P and Si dropped to near the ecological threshold for growth of diatoms. The N:P ratio increased from 4 in 1984 to over 16 in 2000. The climate trend coefficients, R_xt, for these time series are all above 0.43 with significance levels of 95%, except for salinity. The increases in T were consistent with the recent climate warming in northern China and the adjacent seas, i.e. the Bohai Sea and the East China Sea. The reduction of DO was probably attributable to the increase in T and decrease in primary production in these regions. The positive trend of DIN was mainly attributable to precipitation and partly to Changjiang River discharge. The negative concentration trends of P and Si were due to the decreases in their concentrations in seawater that flowed to the Yellow Sea from the Bohai Sea. As a result, N:P ratios greatly increased in the seawater of the Yellow Sea.Moreover, some important responses of the Yellow Sea ecosystems to the changes in physical variables and chemical biogenic elements were obviously displayed. These responses include strengthening nutrient limitation, decreasing chlorophyll a, primary production and phytoplankton abundance, succession of dominant phytoplankton species from diatoms to non-diatoms, changes in fish community structure and species diversity.     

Gulf of Gdansk spring bloom physical, bio-optical, biological modelling and contact data assimilation

In this study we propose a model of phytoplankton population dynamics in the marine ecosystem, which includes physical, biological and bio-optical parts. As an example we simulate the abnormal 1993 Gulf of Gdansk spring bloom, when extremely high chlorophyll concentrations were observed. For the one-dimensional model we use two different methods of contact chlorophyll observation assimilation to fit a model of “in situ” data. The results are compared with two-dimensional ecosystem modelling based on a barotropic model of wind-driven circulation without assimilation.     

Primary productivity of the Palmer Long Term Ecological Research Area and the Southern Ocean

A major objective of the Palmer Long Term Ecological Research (Palmer LTER) project is to obtain a comprehensive understanding of the various components of the Antarctic marine ecosystem. Phytoplankton production plays a key role in this so-called high nutrient, low chlorophyll environment, and factors that regulate production include those that control cell growth (light, temperature, and nutrients) and those that control cell accumulation rate and hence population growth (water column stability, grazing, and sinking). Sea ice mediates several of these factors and frequently conditions the water column for a spring bloom which is characterized by a pulse of production restricted in both time and space. This study models the spatial and temporal variability of primary production within the Palmer LTER area west of the Antarctic Peninsula and discusses this production in the context of historical data for the Southern Ocean. Primary production for the Southern Ocean and the Palmer LTER area have been computed using both light-pigment production models [Smith, R.C., Bidigare, R.R., Prézelin, B.B., Baker, K.S., Brooks, J.M., 1987. Optical characterization of primary productivity across a coastal front. Mar. Biol. (96), 575–591; Bidigare, R.R., Smith, R.C., Baker, K.S., Marra, J., 1987. Oceanic primary production estimates from measurements of spectral irradiance and pigment concentrations. Global Biogeochem. Cycles (1), 171–186; Morel, A., Berthon, J.F., 1989. Surface pigments, algal biomass profiles and potential production of the euphotic layer—relationships reinvestigated in view of remote-sensing applications. Limnol. Oceanogr. (34), 1545–1562] and an ice edge production model [Nelson, D.M., Smith, W.O., 1986. Phytoplankton bloom dynamics of the western Ross Sea ice edge: II. Mesoscale cycling of nitrogen and silicon. Deep-Sea Res. (33), 1389–1412; Wilson, D.L., Smith, W.O., Nelson, D.M., 1986. Phytoplankton bloom dynamics of the Western Ross Sea ice edge: I. primary productivity and species-specific production. Deep-Sea Res., 33, 1375–1387; Smith, W.O., Nelson, D.M., 1986. Importance of ice edge phytoplankton production in the Southern Ocean. BioScience (36), 251–257]. Chlorophyll concentrations, total photosynthetically available radiation (PAR) and sea ice concentrations were derived from satellite data. These same parameters, in addition to hydrodynamic conditions, have also been determined from shipboard and Palmer Station observations during the LTER program. Model results are compared, sensitivity studies evaluated, and productivity of the Palmer LTER region is discussed in terms of its space time distribution, seasonal and interannual variability, and overall contribution to the marine ecology of the Southern Ocean.     

Chlorophyll estimation in the Southeastern Mediterranean using CZCS images: adaptation of an algorithm and its validation

Region-specific algorithm to derive chlorophyll concentration from Coastal Zone Color Scanner (CZCS) images was adopted following field measurements in the Southeastern Mediterranean Sea, in June 1992. Reflectance in the range 400 to 700 nm, chlorophyll concentration and seston dry weight were measured simultaneously at 21 stations. The optical model used based on algorithm for Case I waters, using shipboard simulation of CZCS channels for adaptation of parameters. The resulting model was validated with an independent data sets from series of cruises in the Southeastern Mediterranean (1981–1984). Concurrent CZCS images were used for retrieval of the radiance ratio Channel 1/Channel 3, for each station. The resulting ratios were then used to calculate chlorophyll concentrations. A comparison between those calculated concentrations and the actual chlorophyll measurements showed that an estimation error of chlorophyll as low as 0.013 mg m⁻³ was achieved.     

The Mercator global ocean operational analysis system: Assessment and validation of an 11-year reanalysis

This paper presents Prototype Système 2 Global (PSY2G), the first Mercator global Ocean General Circulation Model (OGCM) to assimilate along-track sea level anomaly (SLA) satellite data. Based on a coarse resolution ocean model, this system was developed mainly for climatic purposes and will provide, for example, initial oceanic states for coupled ocean-atmosphere seasonal predictions. It has been operational since 3 September 2003 and produces an analysis and a two-week forecast for the global ocean every week. The PSY2G system uses an incremental assimilation scheme based on the Cooper and Haines [Cooper, M., Haines, K., 1996. Data assimilation with water property conservation. J. Geophys. Res., 101, 1059-1077.] lifting–lowering of isopycnals. The SLA increment is obtained using an optimal interpolation method then the correction is partitioned into baroclinic and barotropic contributions. The baroclinic ocean state correction consists of temperature, salinity and geostrophic velocity increments and the barotropic correction is a barotropic velocity increment. A reanalysis (1993–2003) was carried out that enabled the PSY2G system to perform its first operational cycle. All available SLA data sets (TOPEX/Poséïdon, ERS2, Geosat-Follow-On, Jason1 and Envisat) were assimilated for the 1993–2003 period. The major objective of this study is to assess the reanalysis from both an assimilation and a thermodynamic point of view in order to evaluate its realism, especially in the tropical band which is a key region for climatic studies. Although the system is also able to deliver forecasts, we have mainly focused on analysis. These results are useful because they give an a priori estimation of the qualities and capabilities of the operational ocean analysis system that has been implemented. In particular, the reanalysis identifies some regional biases in sea level variability such as near the Antarctic Circumpolar Current, in the eastern Equatorial Pacific and in the Norwegian Sea (generally less than 1 cm) with a small seasonal cycle. This is attributed to changes in mean circulation and vertical stratification caused by the assimilation methodology. But the model's low resolution, inaccurate physical parameterisations (especially for ocean–ice interactions) and surface atmospheric forcing also contribute to the occurrence of the SLA biases. A detailed analysis of the thermohaline structure of the ocean reveals that the isopycnal lifting–lowering tends to diffuse vertically the main thermocline. The impact on temperature is that the surface layer (0–200 m) becomes cooler whereas in deeper waters (from 500 to 1500 m), the ocean becomes slightly warmer. This is particularly true in the tropics, between 30°N and 30°S. However it can be demonstrated that the assimilation improves the variability in both surface currents and sub-surface temperature in the Equatorial Pacific Ocean.