An eddy-permitting, coupled ecosystem-circulation model of the Arabian Sea: comparison with observations

A nitrogen-based, pelagic ecosystem model has been coupled with an eddy-permitting ocean general circulation model of the Arabian Sea, and the results are compared with observations. The seasonal variability simulated by the model is in good agreement with observations: during the southwest monsoon season, phytoplankton increases in the western Arabian Sea due to upwelling along the coast; during the northeast monsoon season, phytoplankton abundance is large in the northern Arabian Sea because of the enhanced nitrate entrained by relatively deep vertical mixing. Two major differences are, however, found in the basin-wide comparison between model results and observations: an unrealistic nitrate maximum in the subsurface layer of the northern Arabian Sea and too low primary production in oligotrophic regimes. The former may be attributed to the lack of denitrification in the model. Possible causes for the latter include the present model's underestimation of fast nutrient recycling, the neglect of carbon fixation decoupled from nitrogen uptake and of nitrogen fixation, and inadequate nitrate entrainment by mixed layer deepening. The rate at which simulated nitrate increases in the northern Arabian Sea is 11–24 TgN/year, and should correspond to the denitrification rate integrated over the northern Arabian Sea assuming that the loss of nitrogen through denitrification is balanced by advective input. The model does not reproduce the observed phytoplankton bloom in the late southwest monsoon season. Possible causes are that the mixed layer may be too shallow in summer and that the horizontal transport of nitrate from the coast of Oman may be too weak. Sensitivity experiments demonstrate a strong dependence of the simulated primary productivity on the vertical mixing scheme and on the inclusion of a fast recycling loop in the ecosystem model.     

Modelling the spatial and temporal variability of the Cretan Sea ecosystem

The ecosystem function of the oligotrophic Cretan Sea is explored through the development and application of a 3D ecological model. The simulation system comprises of two on-line coupled submodels: the 3D Princeton Ocean Model (POM) and the 1D European Regional Seas Ecosystem Model (ERSEM) adapted to the Cretan Sea. For the tuning and initialisation of the ecosystem parameters, the 1D version of the biogeochemical model is used.After a model spin up period of 10 years to reach a quasi-steady state, the results from an annual simulation are presented. A cost function is used as validation method for the comparison of model results with field data. The estimated annual primary and bacteria production are found to be in the range of the reported values. Simulation results are in good agreement with in situ data illustrating the role of the physical processes in determining the evolution and variability of the ecosystem.     

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).     

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.     

Is DOC the main source of organic matter remineralization in the ocean water column?

Recent interpretations of carbon flux data and deep-sea processes have led to a reconsideration of the role of dissolved organic carbon (DOC) in supporting water column remineralization and other mid-water biogeochemical transformations (Suzuki et al., 1985; Cho and Azam, 1988; Karl et al., 1988; Christensen et al., 1989; Naqvi and Shailaja, 1993). To date, there have been no direct comparisons of particulate carbon flux data with water column metabolic rates. Here, for the first time, particulate carbon flux and respiratory electron transport activity (from which metabolic CO₂ production is derived), have been monitored simultaneously for one year in the same area of the northwestern Mediterranean Sea. In the aphotic layer (200–1000 m), particulate organic carbon (POC) can support only 20% of the overall organic matter remineralization. Remineralization rates are consistent with recent calculations of DOC exported from the euphotic layer in this area, confirming the vital importance of DOC in maintaining deep-water metabolism. This finding would apply to other regions of mesotrophic and oligotrophic production and thus affect our understanding of carbon recycling in the water column, new production and O₂ utilization.     

Applications of a new ecosystem model to study the dynamics of phytoplankton and nutrients in the Ariake Sea,west coast of Kyushu,Japan

In this study, we present the development and application of a new ecosystem model coupled with a hydrodynamic model to describe the important physical, chemical and biological processes of an ecosystem in the marine environment, the Ariake Sea in the west coast of Kyushu, Japan. The model was calibrated and validated using in-situ field measurements from various monitoring stations in the sea. The presented results covered the period from January 1991 to December 2000. The results showed that chlorophyll-a, nutrients and dissolved oxygen levels varied seasonally in response to weather and boundary condition. Through this study, the model was shown to be able to handle the flooding and drying processes that usually exist and play an essential role over the estuarine-tidal flats of the sea.     

Annual uptake of atmospheric CO2 by the Weddell Sea derived from a surface layer balance, including estimations of entrainment and new production

Data from two cruises, one in April/May 1996 and one in December/January 1993, covering the same wide area in the offshore Weddell Sea, were used to derive the annual extent of entrainment and the capacity of the biological pump. The former property was obtained with the help of dissolved oxygen data, whereas the latter was approximated with nutrients. Especially the data from April/May, representing the initial state of the winter surface layer, were crucial to assess the annual extent of these processes. The results were applied to our carbon dioxide data. The annual increase of the Total CO₂ (TCO₂) concentration in the surface layer due to vertical transport amounts to 16.3 μmol kg⁻¹. An entrainment rate of deep water in the surface layer amounting to 35±10 m yr⁻¹ was deduced. The compensating, biologically mediated TCO₂ reduction was calculated to be larger than the TCO₂ increase due to vertical transport. Since the balance of these two processes determines whether the Weddell Sea is a source or a sink of CO₂, this indicates that the Weddell Sea, albeit upwelling area, is definitely a sink for atmospheric CO₂ on an annual basis. This conclusion is further supported by contemplations that the biological drawdown of CO₂ in the Weddell Sea as a whole is probably underestimated by our calculations. The new production for the Weddell Sea on a per unit area basis was found to be much higher than that for the Antarctic Ocean, when the latter value is being obtained by traditional biological methods. On the other hand, the CO₂ uptake by the Weddell Sea on a per unit area basis is somewhat smaller than the CO₂ uptake by the world ocean.     

Improving the physics of a coupled physical–biogeochemical model of the North Atlantic through data assimilation: Impact on the ecosystem

Several studies on coupled physical–biogeochemical models have shown that major deficiencies in the biogeochemical fields arise from the deficiencies in the physical flow fields. This paper examines the improvement of the physics through data assimilation, and the subsequent impact on the ecosystem response in a coupled model of the North Atlantic. Sea surface temperature and sea surface height data are assimilated with a sequential method based on the SEEK filter adapted to the coupling needs. The model domain covers the Atlantic from 20°S to 70°N at eddy-permitting resolution. The biogeochemical model is a NPZD-DOM model based on the P3ZD formulation. The results of an annual assimilated simulation are compared with an annual free simulation.With assimilation, the representation of the mixed layer depth is significantly improved in mid latitudes, even though the mixed layer depth is generally overestimated compared to the observations. The representation of the mean and variance of the currents is also significantly improved.The nutrient input in the euphotic zone is used to assess the data assimilation impact on the ecosystem. Data assimilation results in a 50% reduction of the input due to vertical mixing in mid-latitudes, and in a four- to six-fold increase of the advective fluxes in mid-latitudes and subtropics. Averaged zonally, the net impact is a threefold increase for the subtropical gyre, and a moderate (20–30%) decrease at mid and high latitudes.Surface chlorophyll concentration increases along the subtropical gyre borders, but little changes are detected at mid and high latitudes. An increase of the primary production appears along the Gulf Stream path, but it represents only 12% on average for mid and high latitudes. In the subtropical gyre centre, primary production is augmented but stays underestimated (20% of observations). These experiments show the benefits of physical data assimilation in coupled physical–biogeochemical applications.     

Enrichment, concentration and retention processes in relation to anchovy (Engraulis ringens) eggs and larvae distributions in the northern Humboldt upwelling ecosystem

A Lagrangian model is used to simulate and quantify in the northern Humboldt upwelling ecosystem the processes of enrichment, concentration and retention, identified by Bakun [Bakun, A., 1996. Patterns in the ocean. Ocean processes and marine population dynamics. University of California Sea Grant, California, USA, in cooperation with Centro de Investigaciones Biologicas de Noroeste, La Paz, Baja California Sur, Mexico, 323 pp.] as being important for the survival and recruitment of early life stages of pelagic fish. The method relies on tracking the positions of particles within water velocity fields generated by a three-dimensional hydrodynamic model. Simple criteria for considering particles as participating to enrichment, concentration or retention are used to derive indices of the three processes. We analyse the spatial distribution of and seasonal variability in these indices. The results are discussed in relation to anchovy (Engraulis ringens) eggs and larvae distributions off Peru, and to a comparable study conducted in the southern Benguela upwelling ecosystem.     

Oceanic CO2 uptake and biogeochemical variability during the formation of the Eastern North Atlantic Central water under two contrasting NAO scenarios

Long-term variability of the biogeochemical properties during the formation of central waters in the Eastern North Atlantic were analyzed between 42–47°N and 10–20°W from the dataset gathered during the Galicia VII (GVII) and C. Darwin 58/59 (CD58/59) cruises. These cruises that showed important changes in the thermohaline properties and the nutrient abundance of the upper layers were carried out under contrasting conditions of the North Atlantic Oscillation (NAO) index. The different climate forcing led a meridional shift of the transition zone between the formation regions of subpolar and subtropical Eastern North Atlantic Central Water (ENACWp and ENACWt, respectively). This displacement conditioned the presence of each ENACW in the study region and so the thermohaline and biogeochemical properties. The effect of the observed variability at decadal scale on the air–sea CO₂ gradient (ΔfCO₂) and exchange (F_CO2) was analyzed using 1D model approach throughout 11 weekly-steps that simulated the development of a spring bloom during the shoaling of the mixed layer. The outputs of the model showed an intensification of the ocean CO₂ uptake due to higher biological CO₂ drawdown, during positive NAO conditions and its weakening under negative NAO influence.     

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.     

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.     

Carbon and nutrient dynamics in a hypertrophic lagoon in southwestern Taiwan

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

The influence of upwelling and entrainment on the algal bloom in the Baltic Sea

Hydrodynamic processes control many geochemical and ecological processes in the sea. In this paper, the influence of up- and downwelling and entrainment on the ecosystem components are studied. The ecohydrodynamic model was initially used to simulate the whole Baltic Sea to get the boundary conditions for the Gulf of Riga. Then, to study the influence of different hydrodynamic conditions on the algal bloom, three simulations were made for the Gulf of Riga using different boundary and entrainment conditions. It appears that upwelling in the gulf was strongly dependent on open boundary conditions between the Baltic Proper and the gulf. The vertical transport in the Gulf of Riga was many times more intensive in the calculation system Baltic Proper and Gulf of Riga, than in the case where only the Gulf of Riga was simulated. The blue–green algal bloom was influenced by the vertical transport due to different nutrients' limitation mechanism.     

Carbon flux in ice–ocean–plankton systems of the Bellingshausen Sea during a period of ice retreat

Most analyses of marine microbial systems in the seasonally ice covered areas of the Southern Ocean have been based on data from the major embayment areas of the Ross and Weddell Seas. In this study data were collected at stations covering a range of regimes from full ice cover through to open water in the Bellingshausen Sea. A major feature of the production system was a rapid retreat of the ice-edge, which uncoupled the marginal ice zone from a phytoplankton bloom which remained associated with a frontal system. This bloom was maintained, and probably initiated, in an unusual environment generated by the interaction between the marginal ice zone and the front. Size-based analyses of the microbial system were derived for ice-covered, recently ice-covered and open water sites. Estimates of standing stocks and key rate processes were combined to produce a single food web network for each station. The under-ice system was one of low production and low recycling but apparently high retention. As the ice retreated the microbial systems to the north began to develop, but these were constrained by grazing pressure. The bloom in the area appeared to be sustained even though estimated losses were far higher than production, although the high sedimentation losses expected were not observed. The carbon flow networks are discussed in relation to the environmental changes and the interaction of the marginal ice zone and the frontal system appears crucial to the phytoplankton. Microzooplankton grazing is implicated as a major controlling factor. The local microbial dynamics are strongly influenced by material which was produced at an earlier time and somewhere else in the Southern Ocean.     

pH variability and CO2 induced acidification in the North Sea

A coupled carbonate system–marine ecosystem–hydrodynamic model is used to simulate the temporal and spatial variability in pH across the southern North Sea as it relates to the environmental and biological processes affecting CO₂, namely, photosynthesis and respiration, riverine boundary conditions and atmospheric CO₂ concentrations. Annual pH ranges are found to vary from < 0.2 in areas of low biological activity to > 1.0 in areas influenced by riverine signals, consistent with observations and previous studies. It is shown that benthic, as well as pelagic, activity is an important factor in this variability. The acidification of the region due to increased fluxes of atmospheric CO₂ into the marine system is calculated and shown to exceed, on average, 0.1 pH units over the next 50 years and result in a total acidification of 0.5 pH units below pre-industrial levels at atmospheric CO₂ concentrations of 1000 ppm. The potential for measurable changes in biogeochemistry are demonstrated by simulating the observed inhibition of pelagic nitrification with decreasing pH. However, we conclude that there is a lack of knowledge of how acidification might affect the complex interaction of processes that govern marine biogeochemical cycles and a consequent need for further research and observations.     

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.     

Bayesian calibration of mechanistic aquatic biogeochemical models and benefits for environmental management

Aquatic biogeochemical models have been an indispensable tool for addressing pressing environmental issues, e.g., understanding oceanic response to climate change, elucidation of the interplay between plankton dynamics and atmospheric CO₂ levels, and examination of alternative management schemes for eutrophication control. Their ability to form the scientific basis for environmental management decisions can be undermined by the underlying structural and parametric uncertainty. In this study, we outline how we can attain realistic predictive links between management actions and ecosystem response through a probabilistic framework that accommodates rigorous uncertainty analysis of a variety of error sources, i.e., measurement error, parameter uncertainty, discrepancy between model and natural system. Because model uncertainty analysis essentially aims to quantify the joint probability distribution of model parameters and to make inference about this distribution, we believe that the iterative nature of Bayes' Theorem is a logical means to incorporate existing knowledge and update the joint distribution as new information becomes available. The statistical methodology begins with the characterization of parameter uncertainty in the form of probability distributions, then water quality data are used to update the distributions, and yield posterior parameter estimates along with predictive uncertainty bounds. Our illustration is based on a six state variable (nitrate, ammonium, dissolved organic nitrogen, phytoplankton, zooplankton, and bacteria) ecological model developed for gaining insight into the mechanisms that drive plankton dynamics in a coastal embayment; the Gulf of Gera, Island of Lesvos, Greece. The lack of analytical expressions for the posterior parameter distributions was overcome using Markov chain Monte Carlo simulations; a convenient way to obtain representative samples of parameter values. The Bayesian calibration resulted in realistic reproduction of the key temporal patterns of the system, offered insights into the degree of information the data contain about model inputs, and also allowed the quantification of the dependence structure among the parameter estimates. Finally, our study uses two synthetic datasets to examine the ability of the updated model to provide estimates of predictive uncertainty for water quality variables of environmental management interest.     

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.     

Winter phytoplankton blooms in the shallow mixed layer of the South China Sea enhanced by upwelling

The South China Sea (SCS) is the largest marginal sea in the world. Previous studies, including recent intensive paleo-oceanographic studies, suggest that the SCS is sensitive to many types of physical forcing on the short-term (e.g., internal waves and tides, mesoscale eddies, typhoons, etc.), annual (e.g., monsoon), inter-annual (e.g., El Niño), and very long-term (e.g., climate change) time scales. To better understand how various types of physical forcing influence biogeochemical cycles in the water column, a time-series study was initiated. Bimonthly hydrographic surveys occupied stations in the subtropical–tropic SCS at 19°N, 118.5°E. Results suggest that the Southeast Asian monsoons, northeasterly from October to April and southwesterly from May to September, have important effects on biogeochemical cycles in the upper water column. Hydrographic data showed that the mixed layer depth was much shallower in winter than in other seasons. During the winter monsoon period, the nitricline became shallower and upwelling sustained an elevated phytoplankton standing stock. Mean chlorophyll concentrations (0.65 mg Chl m^− 3) in winter were 8 times higher than in summer, and the integrated primary productivity over the euphotic zone reached as high as ca. 684 mg C m^− 2 day^− 1 in winter. The upwelling is produced by convergence of currents in the cyclonic gyre near the Luzon Strait, where the Kuroshio intrudes. In summer the current reverses following the wind change. The nitricline is depressed as downwelling occurs off northwest Luzon, resulting in strong nutrient limitation and very low chlorophyll concentrations.