Seasonal and interannual variation of physical and biological processes during 1994–2001 in the Sea of Japan/East Sea: A three-dimensional physical–biogeochemical modeling study

A Pacific basin-wide physical–biogeochemical model has been used to investigate the seasonal and interannual variation of physical and biological fields with analyses focusing on the Sea of Japan/East Sea (JES). The physical model is based on the Regional Ocean Model System (ROMS), and the biogeochemical model is based on the Carbon, Si(OH)₄, Nitrogen Ecosystem (CoSiNE) model. The coupled ROMS–CoSiNE model is forced with the daily air–sea fluxes derived from the National Centers for Environmental Prediction (NCEP) and the National Center for Atmospheric Research (NCAR) reanalysis for the period of 1994 to 2001, and the model results are used to evaluate climate impact on nutrient transport in Mixed Layer Depth (MLD) and phytoplankton spring bloom dynamics in the JES.The model reproduces several key features of sea surface temperature (SST) and surface currents, which are consistent with the previous modeling and observational results in the JES. The calculated volume transports through the three major straits show that the Korea Strait (KS) dominates the inflow to the JES with 2.46 Sv annually, and the Tsugaru Strait (TS) and the Soya Strait (SS) are major outflows with 1.85 Sv and 0.64 Sv, respectively. Domain-averaged phytoplankton biomass in the JES reaches its spring peak 1.8 mmol N m^− 3 in May and shows a relatively weak autumn increase in November. Strong summer stratification and intense consumption of nitrate by phytoplankton during the spring result in very low nitrate concentration at the upper layer, which limits phytoplankton growth in the JES during the summer. On the other hand, the higher grazer abundance likely contributes to the strong suppression of phytoplankton biomass after the spring bloom in the JES. The model results show strong interannual variability of SST, nutrients, and phytoplankton biomass with sudden changes in 1998, which correspond to large-scale changes of the Pacific Decadal Oscillation (PDO). Regional comparisons of interannual variations in springtime were made for the southern and northern JES. Variations of nutrients and phytoplankton biomass related to the PDO warm/cold phase changes were detected in both the southern and northern JES, and there were regional differences with respect to the mechanisms and timing. During the warm PDO, the nutrients integrated in the MLD increased in the south and decreased in the north in winter. Conversely, during the cold PDO, the nutrients integrated in the MLD decreased in the south and increased in the north. Wind divergence/convergence likely drives the differences in the southern and northern regions when northerly and northwesterly monsoon dominates in winter in the JES. Subjected to the nutrient change, the growth of phytoplankton biomass appears to be limited neither by nutrient nor by light consistently both in the southern and northern regions. Namely, the JES is at the transition zone of the lower trophic-level ecosystem between light-limited and nutrient-limited zones.     

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.     

Modelling Pseudocalanus elongatus stage-structured population dynamics embedded in a water column ecosystem model for the northern North Sea

This paper outlines an approach to couple a structured zooplankton population model with state variables for eggs, nauplii, two copepodites stages and adults adapted to Pseudocalanus elongatus into the complex marine ecosystem model ECOHAM2 with 13 state variables resolving the carbon and nitrogen cycle. Different temperature and food scenarios derived from laboratory culture studies were examined to improve the process parameterisation for copepod stage dependent development processes. To study annual cycles under realistic weather and hydrographic conditions, the coupled ecosystem–zooplankton model is applied to a water column in the northern North Sea. The main ecosystem state variables were validated against observed monthly mean values. Then vertical profiles of selected state variables were compared to the physical forcing to study differences between zooplankton as one biomass state variable or partitioned into five population state variables. Simulated generation times are more affected by temperature than food conditions except during the spring phytoplankton bloom. Up to six generations within the annual cycle can be discerned in the simulation.     

Numerical analysis of water circulation and thermohaline structures in the Caspian Sea

The Caspian Sea has a unique ecosystem that consists of endemic species. The deterioration of the unique ecosystem has become increasingly worrisome since a wide variety of pollutants have been released into the water. Water circulation plays a key role in advection and diffusion of these pollutants. In the present study, water circulation and thermohaline structures in the Caspian Sea were analyzed by means of a three dimensional numerical simulation. The effects of meteorological changes, river inflow, and an icing event were taken into account as boundary conditions. Numerical simulation was carried out for 20 years to achieve stable seasonal variations in model variables. As a result, the horizontal distributions of water temperature and salinity could be reproduced; the gradient of water temperature in the north–south direction, the decrease in water temperature along the east coast of the middle Caspian Sea due to coastal upwelling, and low salinity in the northern Caspian Sea. The icing event kept the water temperature in the northern Caspian Sea from decreasing to an unrealistic value. The observed cyclonic gyres were basically formed by the density-driven current due to thermohaline structure.     

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.     

Study of the seasonal cycle of the biogeochemical processes in the Ligurian Sea using a 1D interdisciplinary model

A one-dimensional coupled physical–biogeochemical model has been built to study the pelagic food web of the Ligurian Sea (NW Mediterranean Sea). The physical model is the turbulent closure model (version 1D) developed at the GeoHydrodynamics and Environmental Laboratory (GHER) of the University of Liège. The ecosystem model contains 19 state variables describing the carbon and nitrogen cycles of the pelagic food web. Phytoplankton and zooplankton are both divided in three size-based compartments and the model includes an explicit representation of the microbial loop including bacteria, dissolved organic matter, nano-, and microzooplankton. The internal carbon/nitrogen ratio is assumed variable for phytoplankton and detritus, and constant for zooplankton and bacteria. Silicate is considered as a potential limiting nutrient of phytoplankton's growth. The aggregation model described by Kriest and Evans in (Proc. Ind. Acad. Sci., Earth Planet. Sci. 109 (4) (2000) 453) is used to evaluate the sinking rate of particulate detritus. The model is forced at the air–sea interface by meteorological data coming from the “Côte d'Azur” Meteorological Buoy. The dynamics of atmospheric fluxes in the Mediterranean Sea (DYFAMED) time-series data obtained during the year 2000 are used to calibrate and validate the biological model. The comparison of model results within in situ DYFAMED data shows that although some processes are not represented by the model, such as horizontal and vertical advections, model results are overall in agreement with observations and differences observed can be explained with environmental conditions.     

An integrated water/sediment approach to study plankton (a case study in the southern Adriatic Sea)

In marine coastal areas many planktonic species produce resting stages (cysts) that sink to the bottom. Integrated sampling from both the water column (to collect active stages), and sediments (to collect cysts), could be useful to achieve more complete information about plankton composition.In the framework of the “INTERREG II Albania-Italy Project” an oceanographic survey was carried out aboard the r/v “Italica” from 20 to 31 October 2000. The survey interested the northern Albanian coast (Gulf of Drin) and the northern Apulian coast (Gulf of Manfredonia) on the opposite sides of the South Adriatic Sea. The plankton was collected from 14 stations. A total of 188 categories were recognized in plankton samples. Among those categories, 130 species were recognized (87 of phytoplankton, 43 of microzooplankton), and only 53 (40.8%) resulted common to both the Adriatic sides. A total of 69 cyst morphotypes were recognized in sediment samples; 38 of them were classified at the level of species. A statistical analysis of the micro-zooplankton species abundance showed a segregation of the two areas better than that obtained with the phytoplankton. Cyst distribution in the sediments showed a good gulf-segregation too. In addition, they allowed us to find complementary information, particularly for dinoflagellates. The most abundant species in the water column were not equally dominant as resting stages in the sediments. Sediment sampling allowed further information about the composition of the plankton communities, and suggested us to search for a new method to enhance the yield of less abundant cysts.     

Hydrography and biogeochemistry of the north western Bay of Bengal and the north eastern Arabian Sea during winter monsoon

The north eastern Arabian Sea and the north western Bay of Bengal within the Indian exclusive economic zone were explored for their environmental characteristics during the winter monsoons of 2000 and 2001 respectively. The two regions were found to respond paradoxically to comparable intensities of the atmospheric forcing. There is an asymmetry in the net heat exchange of these two basins with atmosphere because of the varying thickness of barrier layer. During winter, the convective mixing in the Arabian Sea is driven by net heat loss from the ocean, whereas the Bay of Bengal does not contribute to such large heat loss to the atmosphere. It appears that the subduction of high saline Arabian Sea water mass is the mechanism behind the formation of a barrier layer in the northeast Arabian Sea; whereas that in the Bay of Bengal and the southeast Arabian Sea are already established as due to low saline water mass. The weak barrier layer in the Arabian Sea yields to the predominance of convective mixing to bring in nitrate-rich waters from the deeper layers to the surface, thereby supporting enhanced biological production. On the other hand, the river discharge into the Bay of Bengal during this period results in the formation of a thick and stable barrier layer, which insulates vertical mixing and provide oligotrophic condition in the Bay.     

Short time-scale analysis of the NW Mediterranean ecosystem during summer–autumn transition: A 1D modelling approach

Modelling was used as a tool to better understand the physical and biological processes observed during the multidisciplinary cruise DYNAPROC 2 (DYNAmic of rapid PROCesses in the water column), which took place in the Ligurian Sea in September–October 2004. The aim of the cruise was to study the short time-scale physical and biological processes that occur when the ecosystem switches from summer oligotrophy to autumnal mesotrophy. In this study, we have tested two 1D physical–biological coupled models. The first was a classical model in which surface layer dynamics were obtained using the turbulent kinetic energy model of Gaspar [Gaspar et al., 1990]. The simulated food-web took into account ten state variables: three nutrients, three classes of phytoplankton, two classes of zooplankton and two types of detritus. The second model (called IDA, Isopycnals Depth Adjustment) was based on the initial one but it took into account the measured variations of isopycnals depths. The results showed that the IDA model most efficiently reproduced the observed ecosystem dynamics. We have therefore used the IDA model to show that physical processes observed during the cruise had a major effect on biological compartment, mainly on nano- and picophytoplankton.     

Planktonic community structure and carbon cycling in the Arabian Sea as a result of monsoonal forcing: the application of a generic model

The Arabian Sea exhibits a complex pattern of biogeochemical and ecological dynamics, which vary both seasonally and spatially. These dynamics have been studied using a one-dimensional vertical hydrodynamic model coupled to a complex ecosystem model, simulating the annual cycle at three contrasting stations. These stations are characterised by seasonally upwelling, mixed-layer-deepening and a-seasonal oligotrophic conditions, respectively, and coincide with extensively measured stations on the two JGOFS ARABESQUE cruises in 1994. The model reproduces many spatial and temporal trends in production, biomass, physical and chemical properties, both qualitatively and quantitatively and so gives insight into the main mechanisms responsible for the biogeochemical and ecological complexity. Monsoonal systems are typified by classical food web dynamics, whilst intermonsoonal and oligotrophic systems are dominated by the microbial loop. The ecosystem model (ERSEM), developed for temperate regions, is found to be applicable to the Arabian Sea system with little reparameterisation. Differences in in-situ physical forcing are sufficient to recreate contrasting eutrophic and oligotrophic systems, although the lack of lateral terms are probably the greatest source of error in the model. Physics, nutrients, light and grazing are all shown to play a role in controlling production and community structure. Small-celled phytoplanktons are predicted to be dominant and sub-surface chlorophyll maxima are robust centers of production during intermonsoon periods. Analysis of carbon fluxes indicate that physically driven outgassing of CO₂ predominates in monsoonal upwelling systems but ecological activity may significantly moderate CO₂ outgassing in the Arabian Sea interior.     

Environmental factors which affect growth of Japanese common squid, Todarodes pacificus, analyzed by a bioenergetics model coupled with a lower trophic ecosystem model

Bioenergetics model is applied to Japanese common squid, Todarodes pacificus. The temporal change of wet weight of common squid, which migrates in the Sea of Japan, is simulated. The time dependent horizontal distribution of prey is calculated a priori by 3-D coupled physical–biological model. The biological model NEMURO (North Pacific Ecosystem Model for Understanding Regional Oceanography) is used to simulate the lower-trophic ecosystem including three kinds of zooplankton biomass two of which is used as prey of common squid. A bioenergetics model reproduced appropriate growth curve of common squid, migrating in the North Pacific and the Sea of Japan. The results show that the wet weight of common squid in the northern Sea of Japan is heavier than that migrating in the central Sea of Japan, because prey density of the northern Sea of Japan is higher than that of the central Sea of Japan. We also investigate the wet weight anomaly for a global warming scenario. In this case, wet weight of common squid decreases because water temperature exceeds the optimum temperature for common squid. This result indicates that migration route and spawning area of common squid might change with global warming.     

On the key role of nutrient data to constrain a coupled physical–biogeochemical assimilative model of the North Atlantic Ocean

A sequential assimilative system has been implemented into a coupled physical–biogeochemical model (CPBM) of the North Atlantic basin at eddy-permitting resolution (1/4°), with the long-term goal of estimating the basin scale patterns of the oceanic primary production and their seasonal variability. The assimilation system, which is based on the SEEK filter [Brasseur, P., Verron, J., 2006. The SEEK filter method for data assimilation in oceanography: a synthesis. Ocean Dynamics. doi: 10.1007/s10236-006-0080-3], has been adapted to this CPBM in order to control the physical and biogeochemical components of the coupled model separately or in combination. The assimilated data are the satellite Topex/Poseidon and ERS altimetric data, the AVHRR Sea Surface Temperature observations, and the Levitus climatology for salinity, temperature and nitrate.In the present study, different assimilation experiments are conducted to assess the relative usefulness of the assimilated data to improve the representation of the primary production by the CPBM. Consistently with the results obtained by Berline et al. [Berline, L., Brankart, J-M., Brasseur, P., Ourmières, Y., Verron, J., 2007. Improving the physics of a coupled physical–biogeochemical model of the North Atlantic through data assimilation: impact on the ecosystem. J. Mar. Syst. 64 (1–4), 153–172] with a comparable assimilative model, it is shown that the assimilation of physical data alone can improve the representation of the mixed layer depth, but the impact on the ecosystem is rather weak. In some situations, the physical data assimilation can even worsen the ecosystem response for areas where the prior nutrient distribution is significantly incorrect. However, these experiments also show that the combined assimilation of physical and nutrient data has a positive impact on the phytoplankton patterns by comparison with SeaWiFS ocean colour data, demonstrating the good complementarity between SST, altimetry and in situ nutrient data. These results suggest that more intensive in situ measurements of biogeochemical nutrients are urgently needed at basin scale to initiate a permanent monitoring of oceanic ecosystems.     

Nutrient status of the Northeast Water Polynya

The nutrient distribution in the Northeast Water Polynya (NEW) was investigated intensively between the end of May and the beginning of August 1993 during the R/V Polarstern cruise ARK IX. The major characteristics were low initial nitrate concentrations (ca. 4 μM) in the surface mixed layer of the East Greenland Shelf Water, accompanied by high silicate values (ca. 10–14 μM). These concentrations were not reduced by phytoplankton growth. Silicate was rather homogeneously distributed in the entire water column, whereas nitrate increased continuously with depth to about 13 μM. Phosphate concentrations were about 1.1 μM and had a similar distribution to that of silicate. During the course of the summer, nutrients became depleted, and nitrate was exhausted in large parts of the NEW. Silicate was reduced to values of less than 2 μM at some stations which implies that diatom growth continued despite nitrate depletion, ammonium serving as a nitrogen source. The polynya is fertilised by water with the initial nutrient concentrations downstream of the Norske Øer Ice Shelf. This process continuously supplies nutrients to the surface throughout the year and these are transported northward by the anticyclonic surface circulation following the topography of the trough system. The northern boundary of this tongue of relatively nutrient-rich water is controlled by the uptake of nutrients by phytoplankton in summer. Its extemsion is variable due to interactions between biological processes, circulation and ice cover. In the Ob Bank region the nutrient distribution can be altered by the inflow of Polar Water from the north when strong northerly winds prevail as happened during the first part of the study.     

Error quantification of a high resolution coupled hydrodynamic-ecosystem coastal-ocean model: Part3, validation with Continuous Plankton Recorder data

Data collected during the Continuous Plankton Recorder (CPR) survey has been used to validate a three-dimensional hydrodynamic ecosystem model simulation of the North-west European Shelf for the years 1988–89. The CPR time series is unique to the North Atlantic region as a validation tool. Data were extracted from the model to correspond with those collected by the CPR survey, and both the model and survey plankton data were standardised to allow the comparison of model biomass with survey counts. Simple linear regression and absolute error maps provide a qualitative evaluation of spatio-temporal model performance of simulated diatoms, flagellates, total phytoplankton and omnivorous mesozooplankton. Comparisons of z-scores indicate that the model reproduces the main pelagic seasonal features, and there is good correlation between magnitudes of these features with respect to standard deviations from a long-term mean. The model is replicating up to 62% of the mesozooplankton seasonality across the domain, with variable results for the phytoplankton. There are, however, differences in the timing of patterns in plankton seasonality. The validation exercise has highlighted that the spring diatom bloom in the model is too early, suggesting the need to reparameterise the response of phytoplankton to changing light levels in the model. Errors in the north and west of the domain imply that model turbulence and vertical density structure need to be improved to more accurately capture plankton dynamics.     

Description of a flexible and extendable physical–biogeochemical model system for the water column

A modelling system for coupled physical–biogeochemical simulations in the water column is presented here. The physical model component allows for a number of different statistical turbulence closure schemes, ranging from simple algebraic closures to two-equation turbulence models with algebraic second-moment closures. The biogeochemical module consists of models which are based on a number of state variables represented by their ensemble averaged concentrations. Specific biogeochemical models may range from simple NPZ (nutrient–phytoplankton–zooplankton) to complex ecosystem models. Recently developed modified Patankar solvers for ordinary differential equations allow for stable discretisations of the production and destruction terms guaranteeing conservative and non-negative solutions. The increased stability of these new solvers over explicit solvers is demonstrated for a plankton spring bloom simulation. The model system is applied to marine ecosystem dynamics the Northern North Sea and the Central Gotland Sea. Two different biogeochemical models are applied, a conservative nitrogen-based model to the North Sea, and a more complex model including an oxygen equation to the Baltic Sea, allowing for the reproduction of chemical processes under anoxic conditions. For both applications, earlier model results obtained with slightly different model setups could be basically reproduced. It became however clear that the choice for ecosystem model parameters such as maximum phytoplankton growth rates does strongly depend on the physical model parameters (such as turbulence closure models or external forcing).     

Application of a SEEK filter to a 1D biogeochemical model of the Ligurian Sea: Twin experiments and real in-situ data assimilation

The Singular Evolutive Extended Kalman (SEEK) filter has been implemented to assimilate in-situ data in a 1D coupled physical-ecosystem model of the Ligurian Sea. The biogeochemical model describes the partly decoupled nitrogen and carbon cycles of the pelagic food web. The GHER hydrodynamic model (1D version) is used to represent the physical forcings. The data assimilation scheme (SEEK filter) parameterizes the error statistics by means of a set of empirical orthogonal functions (EOFs). Twin experiments are first performed with the aim to choose the suitable experimental protocol (observation and estimation vectors, number of EOFs, frequency of the assimilation,…) and to assess the SEEK filter performances. This protocol is then applied to perform real data assimilation experiments using the DYFAMED data base. By assimilating phytoplankton observations, the method has allowed to improve not only the representation of the phytoplankton community, but also of other variables such as zooplankton and bacteria that evolve with model dynamics and that are not corrected by the data assimilation scheme. The validation of the assimilation method and the improvement of model results are studied by means of suitable error measurements.     

用生态动力学模型分析海带养殖对浮游生态系统的影响

以我国近海典型海带养殖区桑沟湾为例,建立了"浮游植物-海带-浮游动物-营养盐-有机碎屑"生态动力学模型,较好地模拟了不同养殖密度下浮游生态系统结构的变化,通过灵敏度分析,研究状态变量对模型初值和参数变化的响应,提出了参数优化的思路.     

Tropical cyclone genesis over the south China sea

The South China Sea (SCS) is among areas in the Northwest Pacific most frequented by tropical cyclones (TCs) with intensity reaching a tropical storm or stronger. It is also an area of significant TC genesis. In this study, TC genesis in SCS and its monsoonal variability for 1948–2003 are analyzed. Altogether, in May–September (southwest monsoon period) 157 TC geneses have occurred north of 12°N in SCS, while in October–December (northeast monsoon period) 64 out of 65 TC geneses have happened south of 18°N. It is found that the monsoonal characteristics of the SCS basically determine the region of TC genesis in each monsoon season. Winter TC genesis in the SCS happens over the region where the marine environment satisfies the four criterions on, respectively, the sea surface temperature (SST), mid-troposphere relative humidity, vertical shear of the horizontal winds and low-level atmospheric vorticity. During the summer, as the two criterions on SST and the mid-troposphere relative humidity are satisfied for the whole SCS, TC genesis occurs in the region where both the low-level vorticity and the vertical shear satisfy the criterion. In addition, there is likely more TC genesis in the winter during the onset of La Nina, and more TC genesis in the summer following the onset of El Nino.     

One-dimensional modelling of the plankton ecosystem of the north-western Corsican coastal area in relation to meteorological constraints

In order to study the influence of wind mixing on the spring variability of the plankton production of the north western Corsican coastal area, a one-dimensional (1D), vertical, coupled hydrodynamic/biological model (ECOHYDROMV) is used. A hydrodynamic 1D model of the water column with a k–l turbulent closure is applied. The biological model comprises six state variables, representing the plankton ecosystem in the spring period: phytoplankton, copepods, nitrate, ammonium, particulate organic matter of phytoplanktonic origin and particulate organic matter of zooplanktonic origin. The system is influenced by turbulence (expressed by the vertical eddy diffusivity), temperature and irradiance. The model takes into account momentum and heat surface fluxes computed from meteorological data in order to simulate a typical spring atmospheric forcing for the considered area. Results show that primary production vertical structure is characterised by a subsurface maximum which deepens with time and is regulated by the opposite gradients of nitrate concentration and irradiance. Surface plankton productivity is mainly controlled by turbulent vertical transport of nutrients into the mixed layer. The short time scale variability of turbulent mixing generated by the wind appears to be responsible for the plurimodal shape of plankton blooms, observed in the considered area. Furthermore, the model is applied to the study of the spring evolution of the plankton communities off the bay of Calvi (Corsica) for the years 1986 and 1988. In order to initiate and validate the model, time series of hydrological, chemical and biological data have been used. The model reproduces accurately the spring evolution of the phytoplankton biomass measured in situ and illustrates that its strong variability in those years was in close relation to the variability of the wind intensity.     

Modelling the silica pump in the Permanently Open Ocean Zone of the Southern Ocean

A coupled 1D physical–biogeochemical model has been built to simulate the cycles of silicon and of nitrogen in the Indian sector of the Permanently Open Ocean Zone of the Southern Ocean. Based on a simplified trophic network, that includes two size classes of phytoplankton and of zooplankton, and a microbial loop, it has been calibrated by reference to surface physical, chemical and biological data sets collected at the KERFIX time-series station (50°40′S–68°25′E). The model correctly reproduces the high nutrient low chlorophyll features typical of the studied area. In a region where the spring–summer mixed layer depth is usually deeper than 60 m, the maximum of chlorophyll never exceeds 1.5 mg m⁻³, and the annual primary production is only 68 g C m⁻² year⁻¹. In the surface layer nitrate is never exhausted (range 27–23.5 mmoles m⁻³) while silicic acid shows strong seasonal variations (range 5–20 mmoles m⁻³). On an annual basis 71% of the primary production sustained by nanophytoplankton is grazed by microzooplankton. Compared to North Atlantic, siliceous microphytoplankton is mainly prevented from blooming because of an unfavourable spring–summer light-mixing regime. Silicic acid limitation (high half saturation constant for Si uptake: 8 mmoles m⁻³) also plays a major role on diatom growth. Mesozooplankton grazing pressure excerpts its influence especially in late spring. The model illustrates the efficiency of the silica pump in the Southern Ocean: up to 63% of the biogenic silica that has been synthetized in the photic layer is exported towards the deep ocean, while only 11% of the particulate organic nitrogen escapes recycling in the surface layer.