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.     

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.     

Application of an Ensemble Kalman filter to a 1-D coupled hydrodynamic-ecosystem model of the Ligurian Sea

The Ensemble Kalman filter (EnKF) has been applied to a 1-D complex ecosystem model coupled with a hydrodynamic model of the Ligurian Sea. In order to improve the performance of the EnKF, an ensemble subsampling strategy has been used to better represent the covariance matrices and a pre-analysis step for correcting the non-normality of the members distribution has been implemented. Twin experiments have been realized to assess the performance of the developed tool and a real data assimilation experiment has been conducted to hindcast the ecosystem at the Dyfamed site during the year 2000. Finally the performance of the EnKF has been compared with a Singular Evolutive Extended Kalman (SEEK) filter with a fixed basis. We conclude that, on one hand, there is a benefit in using the subsampling strategy and the lognormal transformation with the EnKF, and on the other hand, this filter presents better performance than the fixed basis version of the SEEK filter. However, it also incurs a large computational cost.     

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

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

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.     

Hydrodynamic comparison of a semi-submersible, TLP, and Spar: Numerical study in the South China Sea environment

The South China Sea contains tremendous oil and gas resources in deepwater areas. However, one of the keys for deepwater exploration, the investigation of deepwater floating platforms, is very inadequate. In this paper, the authors studied and compared the hydrodynamics and global motion behaviors of typical deepwater platforms in the South China Sea environment. The hydrodynamic models of three main types of floating platforms, e.g. the Semi-submersible, tension leg platform (TLP), and Truss Spar, which could potentially be utilized in the South China Sea, were established by using the 3-D potential theory. Additionally, some important considerations which significantly influence the hydrodynamics were given. The RAOs in frequency domains as well as global motions in time domains under time-varying wind, random waves, and current in 100-y, 10-y, and 1-y return period environment conditions were predicted, compared, and analyzed. The results indicate that the heave and especially the pitch motion of the TLP are favorable. The heave response of the Truss Spar is perfect and comparable with that of the TLP when the peak period of random waves is low. However, the pitch motion of Truss Spar is extraordinarily larger than that of Semi-submersible and TLP.     

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

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

Importance of ocean circulation in ecological modeling: An example from the North Sea

There is an increasing number of ecological models for the North Sea around. Skogen and Moll (2000) [Skogen, M.D., Moll, A. 2000. Interannual variability of the North Sea primary production: comparison from two model studies. Continental Shelf Research 20 (2), 129–151] compared the interannual variability of the North Sea primary production using two state-of-the-art ecological models, NORWECOM and ECOHAM1. Their conclusion was that the two models agreed on an annual mean primary production, its variability and the timing and size of the peak production. On the other hand, there was a low (even negative dependent of area) correlation in the production in different years between the two models.In the present work, these conclusions are brought further. To try to better understand the observed differences between the two models, the two ecological models are run in an identical physical setting. With such a set-up also the interannual variability between the two models is in agreement, and it is concluded that the single most important factor for a reliable modeling of phytoplankton and nutrient distributions and transports within the North Sea is a proper physical model.     

A model study with light-dependent mortality rates of copepod stages

This paper is based on an advanced ecosystem model of the Baltic Sea (ERGOM [J. Mar. Sys. 25 (3–4) (2005) 405]), but with an increased resolution of the zooplankton stage variable [J. Plankton Res. 23 (2001) 1217; ICES Marine Science 219 (2003) 208]. The model copepods are represented by five stages: eggs, an aggregated variable of nauplii, two aggregated groups of copepodites and adults. The transfer among the stages, i.e., hatching, molting and reproduction, is controlled by food availability and temperature.As usual, the model food web is truncated at the level of zooplankton. The study explores the effects of different parametrization of zooplankton mortality and looks in particular on light-dependent rates. The light climate may serve a proxy for the effects of visual feeding of fish larvae and fish. Different choices of the mortality parameters can result in remarkable differences in abundances and biomass of the model zooplankton and in the timing of its development.It is found that the different choices of mortality affect the development of populations in several ways: Relative small initial differences of abundances at the beginning of the spring bloom are important for the development of the model populations. Higher mortality rates are less important at food rich conditions than at scarce resources. At low phytoplankton levels, the individual development of the copepods through the stages can be faster for elevated mortality rates because then less animals have to share the available food.     

Variability of mesozooplankton spatial distribution in the North Aegean Sea, as influenced by the Black Sea waters outflow

The North Aegean Sea constitutes an important region of the Mediterranean Sea since in its eastern part the mesotrophic, low salinity and relatively cold water from the Black Sea (outflowing from the Dardanelles strait) meets the oligotrophic, warm and very saline water of Levantine origin, thus forming a thermohaline front. Mesozooplankton samples were collected at discrete layers according to the hydrology of the upper 100 m, during May 1997 and September 1998. In May highest biomass and abundance values (up to 66.82 mg m^− 3 and 14,157 ind m^− 3) were detected in the 10–20 m layer (within the halocline) of the stations positioned close to the Dardanelles strait. The front moved slightly southwards in September, characterized by high biomass and abundance values within the halocline layer. The areas moderately or non influenced by Black Sea water revealed lower standing stock values than the frontal area in both cruises and maxima were detected in the uppermost low salinity layer. Samples collected at the stations and/or layers more influenced by Black Sea water were distinguished from those collected at layers and/or stations more affected by Levantine waters in both periods. In May the former samples were characterized by the copepods Acartia clausi, Centropages typicus, Paracalanus parvus. The abundance of the above species decreased gradually with increasing salinity, in the horizontal and/or in the vertical dimension, with a parallel increase of the copepods Oithona plumifera, Oithona copepodites, Oncaea media, Ctenocalanus vanus, Farranula rostrata. During September the frontal area as well as that covered by the modified Black Sea water, were highly dominated by the cladoceran Penilia avirostris and doliolids. For both seasons, MDS plots, issued from the combination of mesozooplankton and water-type data, revealed the gradual differentiation of zooplankton composition from the frontal area towards the area covered by Levantine water, following the spreading and mixing of the Black sea water. The observed temporal and spatial variability in the distribution pattern of mesozooplankton standing stock and species composition seems to depend considerably on the variability of circulation and frontal flows.     

Simulation of the oil spill processes in the Sea of Japan with regional ocean circulation model

A simulation of the movement of spilled oil after the incident of the Russian tanker Nakhodka in the Sea of Japan, in January 1997, was performed by a particle tracking model incorporating advection by currents, random diffusion, the buoyancy effect, the parameterization of oil evaporation, biodegradation, and beaching. The currents advecting spilled oil were defined by surface wind drift superposed on the three-dimensional ocean currents obtained by the Geophysical Fluid Dynamics Laboratory modular ocean model (GFDL MOM), which was forced by the climatological monthly mean meteorological data, or by the European Center for Medium Range Weather Forecasts (ECMWF) daily meteorological data, and assimilated sea surface topography detected by satellite altimeter. A number of experiments with different parameters and situations showed that the wide geographical spread of oil observed is not explained by wind drift alone, and that including the simulated climatological currents gives better results. The combination of surface wind drift and daily ocean currents shows the best agreement between the model and observations except in some coastal areas. The daily meteorological effect on the ocean circulation model results in a stronger variability of currents that closely simulates some features of the nonlinear large-scale horizontal turbulent diffusion of oil. The effect of different parameterizations for the size distribution of model oil particles is discussed. Received for publication on July 26, 1999; accepted on Nov. 17, 1999     

Error quantification of a high-resolution coupled hydrodynamic-ecosystem coastal-ocean model: Part 2. Chlorophyll-a, nutrients and SPM

Marine systems models are becoming increasingly complex and sophisticated, but far too little attention has been paid to model errors and the extent to which model outputs actually relate to ecosystem processes. Here we describe the application of summary error statistics to a complex 3D model (POLCOMS-ERSEM) run for the period 1988–1989 in the southern North Sea utilising information from the North Sea Project, which collected a wealth of observational data. We demonstrate that to understand model data misfit and the mechanisms creating errors, we need to use a hierarchy of techniques, including simple correlations, model bias, model efficiency, binary discriminator analysis and the distribution of model errors to assess model errors spatially and temporally. We also demonstrate that a linear cost function is an inappropriate measure of misfit. This analysis indicates that the model has some skill for all variables analysed. A summary plot of model performance indicates that model performance deteriorates as we move through the ecosystem from the physics, to the nutrients and plankton.     

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.     

A numerical transport model for predicting the distributions of Cd, Cu, Ni, Pb and Zn in the southern North Sea: the sensitivity of model results to the uncertainties in the magnitudes of metal inputs

A new transport model for metals (named NOSTRADAMUS) has been developed to predict concentrations and distributions of Cd, Cu, Ni, Pb and Zn in the southern North Sea. NOSTRADAMUS is comprised of components for water, inorganic and organic suspended particulate matter transport; a primary production module contributes to the latter component. Metal exchange between dissolved (water) and total suspended particulate matter (inorganic + organic) phases is driven by distribution coefficients. Transport is based on an existent 2-D vertically integrated model, incorporating a 35 × 35 km grid. NOSTRADAMUS is largely driven by data obtained during the Natural Environment Research Council North Sea Project (NERC NSP). The sensitivity of model predictions to uncertainties in the magnitudes of metal inputs has been tested. Results are reported for a winter period (January 1989) when plankton production was low. Simulated ranges in concentrations in regions influenced by the largest inflows, i.e. the NE English coast and the Southern Bight, are similar to the ranges in the errors of the concentrations estimated at the northern and southern open sea boundaries of the model. Inclusion of uncertainties with respect to atmospheric (up to ± 54%) and riverine (± 30%) inputs makes little difference to the calculated concentrations of both dissolved and particulate fractions within the southern North Sea. When all the errors associated with the inputs are included there is good agreement between computed and observed concentrations, and that for dissolved and particulate Cd, Cu and Zn, and dissolved Ni and Pb, many of the observations fall within, or are close to, the range of values generated by the model. For particulate Pb, model simulations predict concentrations of the right order, but do not reproduce the large scatter in actual concentrations, with simulated concentrations showing a bias towards lower values compared to those observed. A factor which could have contributed to observed concentrations, and which is not included in the model, is considered to be a substantial benthic input of dissolved lead during this winter period, coupled to a rapid and extensive scavenging of the dissolved lead to particles. Significant reductions in riverine and aeolian inputs of total Cd and Cu of 70% and 50%, respectively, consistent with aims of North Sea Conferences, are predicted to lead to minor decreases (~ 10%) in water column concentrations of dissolved and particulate Cd and Cu, except near river sources, where maximum reductions of ~ 30–40% may occur.     

Regional Ocean Governance in China: An Appraisal of the Clean Bohai Sea Program

In 2001 China embarked on an ambitious program to improve water quality through creation of the Clean Bohai Sea Program (CBSP). It is the first regional ocean governance program in China. The joint conference of central and local government units directed CBSP under the leadership of the State Environment Protection Administration. To meet the environmental objectives established by the joint conference local governments must pay for the program and operate it, but to date they have not been able to fully fund the effort. Although data indicate a decline of polluted discharges to the Bohai, one tabulation shows that the area of the sea that is polluted has increased each year from 2001 to 2005. Future enhancement of CBSP depends on revising incentives to better link local environmental programs and development, establishing funding mechanisms to equitably spread costs throughout watersheds, and creating comprehensive monitoring programs to better appraise program results.     

A Mean Dynamic Topography of the Mediterranean Sea computed from altimetric data, in-situ measurements and a general circulation model

In the Mediterranean Sea, where the mean circulation is largely unknown and characterized by smaller scales and less intensity than in the open ocean, the interpretation of altimetric Sea Level Anomalies (SLA) is rather difficult. In the context of operational systems such as MFS (Mediterranean Forecasting System) or MERCATOR, that assimilate the altimetric information, the estimation of a realistic Mean Dynamic Topography (MDT) consistent with altimetric SLA to be used to reconstruct absolute sea level is a crucial issue. A method is developed here to estimate the required MDT combining oceanic observations as altimetric and in-situ measurements and outputs from an ocean general circulation model (OGCM).In a first step, the average over the 1993–1999 period of dynamic topography outputs from MFS OGCM provides a first guess for the computation of the MDT. Then, in a second step, drifting buoy velocities and altimetric data are combined using a synthetic method to obtain local estimates of the mean geostrophic circulation which are then used to improve the first guess through an inverse technique and map the MDT field (hereafter the Synthetic Mean Dynamic Topography or SMDT) on a 1/8° resolution grid.Many interesting current patterns and cyclonic/anticyclonic structures are visible on the SMDT obtained. The main Mediterranean coastal currents are well marked (as the Algerian Current or the Liguro–Provenço–Catalan Current). East of the Sicily channel, the Atlantic Ionian Stream divides into several main branches crossing the Ionian Sea at various latitudes before joining at 19°E into a unique Mid-Mediterranean Jet. Also, strong signatures of the main Mediterranean eddies are obtained (as for instance the Alboran gyre, the Pelops, Ierapetra, Mersa-Matruh or Shikmona anticyclones and the Cretan, Rhodes or West Cyprius cyclones). Independent in-situ measurements from Sea Campaigns NORBAL in the North Balearic Sea and the North Tyrrhenian Sea and SYMPLEX in the Sicily channel are used to validate locally the SMDT: deduced absolute altimetric dynamic topography compares well with in-situ observations. Finally, the SMDT is used to compute absolute altimetric maps in the Alboran Sea and the Algerian Current. The use of absolute altimetric signal allows to accurately follow the formation and propagation of cyclonic and anticyclonic eddies in both areas.     

Seasonal carbon distribution of copepods in the eastern Weddell Sea, Antarctica

Copepods were sampled by a multiple opening-closing net in the eastern Weddell Sea during various seasons (late winter/early spring, summer, autumn). Total copepod biomass integrated over the upper 1000 m varied seasonally between 1.7 mg C m⁻³ in late winter/early spring and 3.7 mg C m⁻³ in autumn. After the dark season the copepods were rather evenly distributed vertically and highest biomass levels were found in the mid-water layers between about 200 m and 500 m. By contrast, especially in summer but also in autumn copepod biomass concentrated in the uppermost water layer. A total of 64 calanoid species were identified in the upper 1000 m with maximum species numbers in the deepest layer. The large calanoids Calanus propinquus, Calanoides acutus, Metridia gerlachei, Euchaeta antarctica and the small calanoid Microcalanus pygmaeus prevailed and accounted for 60–70% of total copepod biomass, while the small poecilostomatoid Oncaea and the cyclopoid Oithona species comprised about 20%. Hence, the distribution pattern of the entire copepod biomass is strongly influenced by the life cycles of a few dominant species.     

Inorganic carbon fluxes through the boundaries of the Greenland Sea Basin based on in situ observations and water transport estimates

A carbon budget for the exchange of total dissolved inorganic carbon C_T between the Greenland Sea and the surrounding seas has been constructed for winter and summer situations. An extensive data set of C_T collected over the years 1994–1997 within the European Sub-polar Ocean Programmes (ESOP1 and ESOP2) are used for the budget calculation. Based on these data, mean values of C_T in eight different boxes representing the inflow and outflow of water through the boundaries of the Greenland Sea Basin are estimated. The obtained values are then combined with simulated water transports taken from the ESOP2 version of the Miami Isopycnic Coordinate Ocean Model (MICOM). The fluxes of inorganic carbon are presented for three layers; a surface mixed layer, an intermediate layer and a deep layer, and the imbalance in the fluxes are attributed to air–sea exchange, biological fixation of inorganic carbon, and sedimentation. The main influx of carbon is found in the surface and the deep layers in the Fram Strait, and in the surface waters of direct Atlantic origin, whereas the main outflux is found in the surface layer over the Jan Mayen Fracture Zone and the Knipovich Ridge, transporting carbon into the Atlantic Ocean via the Denmark Strait and towards the Arctic Ocean via the Norwegian Sea, respectively. The flux calculation indicates that there is a net transport of carbon out of the Greenland Sea during wintertime. In the absence of biological activity, this imbalance is attributed to air sea exchange, and requires an oceanic uptake of CO₂ of 0.024±0.006 Gt C yr⁻¹. The flux calculations from the summer period are complicated by biological fixation of inorganic carbon, and show that data on organic carbon is required in order to estimate the air–sea exchange in the area.     

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

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

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