Seasonal and interannual variability of primary and export production in the South China Sea: a three-dimensional physical-biogeochemical model study

To investigate the seasonal and interannual variations in biological productivity in the South China Sea (SCS), a Pacific basin-wide physical - biogeochemical model has been developed and used to estimate the biological productivity and export flux in the SCS. The Pacific circulation model, based on the Regional Ocean Model Systems (ROMS), is forced with daily air-sea fluxes derived from the NCEP (National Centers for Environmental Prediction) reanalysis between 1990 and 2004. The biogeochemical processes are simulated with a carbon, Si(OH)(4), and nitrogen ecosystem (CoSiNE) model consisting of silicate, nitrate, ammonium, two phytoplankton groups (small phytoplankton and large phytoplankton), two zooplankton grazers (small micrograzers and large mesozooplankton), and two detritus pools. The ROMS-CoSiNE model favourably reproduces many of the observed features, such as ChI a, nutrients, and primary production (PP) in the SCS. The modelled depth-integrated PP over the euphotic zone (0-125 m) varies seasonally, with the highest value of 386 mg C m (-2) d (-1) during winter and the lowest value of 156 mg C m (-2) d (-1) during early summer. The annual mean value is 196 mg C m (-2) d (-1). The model-integrated annual mean new production (uptake of nitrate), in carbon units, is 64.4 mg C m (-2) d (-1) which yields an f-ratio of 0.33 for the entire SCS. The modelled export ratio (e-ratio: the ratio of export to PP) is 0.24 for the basin-wide SCS. The year-to-year variation of biological productivity in the SCS is weaker than the seasonal variation. The large phytoplankton group tends to dominate over the smaller phytoplankton group, and likely plays an important role in determining the interannual variability of primary and new production.

Organic matter deposition/resuspension in a one-dimensional physical-biogeochemical model. A modelling study

The shallow water configuration of the gulf of Trieste allows the propagation of the stress due to wind and waves along the whole water column down to the bottom. When the stress overcomes a particular threshold it produces resuspension processes of the benthic detritus. The benthic sediments in the North Adriatic are rich of organic matter, transported here by many rivers. This biological active particulate, when remaining in the water, can be transported in all the Adriatic basin by the basin-wide circulation. In this work is presented a first implementation of a resuspension/deposition submodel in the oceanographic coupled physical-biogeochemical 1-dimensional numerical model POM-BFM. At first has been considered the only climatological wind stress forcing, next has been introduced, on the surface, an annual cycle of wave motion and finally have been imposed some exceptional wave event in different periods of the year. The results show a strong relationship between the efficiency of the resuspension process and the stratification of the water column. During summer the strong stratification can contained a great quantity of suspended matter near to the bottom, while during winter even a low concentration of particulate can reach the surface and remains into the water for several months without settling and influencing the biogeochemical system. Looking at the biologic effects, the organic particulate, injected in the water column, allow a sudden growth of the pelagic bacteria which competes with the phytoplankton for nutrients strongly inhibiting its growth. This happen especially during summer when the suspended benthic detritus concentration is greater.

Modelling coupled physical-biogeochemical processes in ice-covered oceans

The last decades have seen a large effort of the scientific community to study and understand the physics of sea ice. We currently have a wide - even though still not exhaustive - knowledge of the sea ice dynamics and thermodynamics and of their temporal and spatial variability. Sea ice biogeochemistry is instead largely unknown. Sea ice algae production may account for up to 25% of overall primary production in ice-covered waters of the Southern Ocean. However, the influence of physical factors, such as the location of ice formation, the role of snow cover and light availability on sea ice primary production is poorly understood. There are only sparse localized observations and little knowledge of the functioning of sea ice biogeochemistry at larger scales. Modelling becomes then an auxiliary tool to help qualifying and quantifying the role of sea ice biogeochemistry in the ocean dynamics. In this thesis, a novel approach is used for the modelling and coupling of sea ice biogeochemistry - and in particular its primary production - to sea ice physics. Previous attempts were based on the coupling of rather complex sea ice physical models to empirical or relatively simple biological or biogeochemical models. The focus is moved here to a more biologically-oriented point of view. A simple, however comprehensive, physical model of the sea ice thermodynamics (ESIM) was developed and coupled to a novel sea ice implementation (BFM-SI) of the Biogeochemical Flux Model (BFM). The BFM is a comprehensive model, largely used and validated in the open ocean environment and in regional seas. The physical model has been developed having in mind the biogeochemical properties of sea ice and the physical inputs required to model sea ice biogeochemistry. The central concept of the coupling is the modelling of the Biologically-Active-Layer (BAL), which is the time-varying fraction of sea ice that is continuously connected to the ocean via brines pockets and channels and it acts as rich habitat for many microorganisms. The physical model provides the key physical properties of the BAL (e.g., brines volume, temperature and salinity), and the BFM-SI simulates the physiological and ecological response of the biological community to the physical enviroment. The new biogeochemical model is also coupled to the pelagic BFM through the exchange of organic and inorganic matter at the boundaries between the two systems . This is done by computing the entrapment of matter and gases when sea ice grows and release to the ocean when sea ice melts to ensure mass conservation. The model was tested in different ice-covered regions of the world ocean to test the generality of the parameterizations. The focus was particularly on the regions of landfast ice, where primary production is generally large. The implementation of the BFM in sea ice and the coupling structure in General Circulation Models will add a new component to the latters (and in general to Earth System Models), which will be able to provide adequate estimate of the role and importance of sea ice biogeochemistry in the global carbon cycle.

A three-dimensional coupled physical-biological model study in the spring of 1993 in the Bohai Sea of China

A three-dimensional (3-D) coupled physical and biological model was used to investigate the physical processes and their influence on the ecosystem dynamics of the Bohai Sea of China. The physical processes include M-2 tide, time - varying wind forcing and river discharge. Wind records from I to 31 May in 1993 were selected to force the model. The biological model is based on a simple, nitrate and phosphate limited, lower trophic food web system. The simulated results showed that variation of residual currents forced by M, tide, river discharge and time-varying wind had great impact on the distribution of phytoplankton biomass in the Laizhou Bay. High phytoplankton biomass appeared in the upwelling region. Numerical experiments based on the barotropic model and baroclinic model with no wind and water discharge were also conducted. Differences in the results by the baroclinic model and the barotropic model were significant: more patches appeared in the baroclinic model comparing with the barotropic model. And in the baroclinic model, the subsurface maximum phytoplankton biomass patches formed in the stratified water.

Decadal reanalysis of biogeochemical indicators and fluxes in the North West European shelf-sea ecosystem

In this paper we present the first decadal reanalysis simulation of the biogeochemistry of the North West European shelf, along with a full evaluation of its skill and value. An error-characterized satellite product for chlorophyll was assimilated into a physical-biogeochemical model of the North East Atlantic, applying a localized Ensemble Kalman filter. The results showed that the reanalysis improved the model predictions of assimilated chlorophyll in 60% of the study region. Model validation metrics showed that the reanalysis had skill in matching a large dataset of in situ observations for ten ecosystem variables. Spearman rank correlations were significant and higher than 0.7 for physical-chemical variables (temperature, salinity, oxygen), ∼0.6 for chlorophyll and nutrients (phosphate, nitrate, silicate), and significant, though lower in value, for partial pressure of dissolved carbon dioxide (∼0.4). The reanalysis captured the magnitude of pH and ammonia observations, but not their variability. The value of the reanalysis for assessing environmental status and variability has been exemplified in two case studies. The first shows that between 340,000-380,000 km2 of shelf bottom waters were oxygen deficient potentially threatening bottom fishes and benthos. The second application confirmed that the shelf is a net sink of atmospheric carbon dioxide, but the total amount of uptake varies between 36-46 Tg C yr−1 at a 90% confidence level. These results indicate that the reanalysis output dataset can inform the management of the North West European shelf ecosystem, in relation to eutrophication, fishery, and variability of the carbon cycle.

Plant-wide predictive control for a thermal power plant based on a physical plant model

A constrained non-linear, physical model-based, predictive control (NPMPC) strategy is developed for improved plant-wide control of a thermal power plant. The strategy makes use of successive linearisation and recursive state estimation using extended Kalman filtering to obtain a linear state-space model. The linear model and a quadratic programming routine are used to design a constrained long-range predictive controller One special feature is the careful selection of a specific set of plant model parameters for online estimation, to account for time-varying system characteristics resulting from major system disturbances and ageing. These parameters act as nonstationary stochastic states and help to provide sufficient degrees-of-freedom to obtain unbiased estimates of controlled outputs. A 14th order non-linear plant model, simulating the dominant characteristics of a 200 MW oil-fired pou er plant has been used to test the NPMPC algorithm. The control strategy gives impressive simulation results, during large system disturbances and extremely high rate of load changes, right across the operating range. These results compare favourably to those obtained with the state-space GPC method designed under similar conditions.

Decadal Variability in the Northeast Pacific in a Physical-Ecosystem Model: Role of Mixed Layer Depth and Trophic Interactions

A basin-wide interdecadal change in both the physical state and the ecology of the North Pacific occurred near the end of 1976. Here we use a physical-ecosystem model to examine whether changes in the physical environment associated with the 1976-1977 transition influenced the lower trophic levels of the food web and if so by what means. The physical component is an ocean general circulation model, while the biological component contains 10 compartments: two phytoplankton, two zooplankton, two detritus pools, nitrate, ammonium, silicate, and carbon dioxide. The model is forced with observed atmospheric fields during 1960-1999. During spring, there is a similar to 40% reduction in plankton biomass in all four plankton groups during 1977-1988 relative to 1970-1976 in the central Gulf of Alaska (GOA). The epoch difference in plankton appears to be controlled by the mixed layer depth. Enhanced Ekman pumping after 1976 caused the halocline to shoal, and thus the mixed layer depth, which extends to the top of the halocline in late winter, did not penetrate as deep in the central GOA. As a result, more phytoplankton remained in the euphotic zone, and phytoplankton biomass began to increase earlier in the year after the 1976 transition. Zooplankton biomass also increased, but then grazing pressure led to a strong decrease in phytoplankton by April followed by a drop in zooplankton by May: Essentially, the mean seasonal cycle of plankton biomass was shifted earlier in the year. As the seasonal cycle progressed, the difference in plankton concentrations between epochs reversed sign again, leading to slightly greater zooplankton biomass during summer in the later epoch.

Physical and Biological Controls on the Latitudinal Asymmetry of Surface Nutrients and pCO(2) in the Central and Eastern Equatorial Pacific

Surface nutrients and dissolved inorganic carbon (DIC) in the central (CEP) and eastern equatorial Pacific (EEP) show much higher concentrations to the south than to the north of the equator. In this study, the physical and biological controls on this asymmetry are investigated using a coupled physical-biogeochemical model. Two numerical experiments are conducted to examine the effects of asymmetrical photosynthetic efficiency (a) due to asymmetrical iron supply about the equator. The experiment with asymmetrical photosynthesis produces improved results as compared with historical observations. A nitrate budget analysis suggests that in the EEP the divergence of upwelling waters controls the surface nitrate asymmetry with additional contribution from the South Equatorial Current (SEC) carrying nutrient-rich Peru upwelling water. The changes of a affect the surface nitrate distribution but not the overall asymmetry. The SEC further carries excess nitrate to the west and thus extends the asymmetry in the east to the CEP. In the CEP, however, stronger northward than southward transport tends to reduce the nitrate asymmetry, while the asymmetrical photosynthesis would help to maintain it. Similar processes also control the distributions of surface silicate and DIC in the equatorial Pacific, which is also affected by the air-sea CO(2) exchange. The asymmetrical photosynthesis influences the distribution of surface DIC, pCO(2), and the air-sea CO(2) flux, by redistributing about 20% CO(2) flux from the north to the south of the equator. Owing to the adjustment of air-sea CO(2) flux, however, the net surface DIC change is smaller than the direct change associated with primary production.

Potential consequences of climate change for primary production and fish production in large marine ecosystems

Existing methods to predict the effects of climate change on the biomass and production of marine communities are predicated on modelling the interactions and dynamics of individual species, a very challenging approach when interactions and distributions are changing and little is known about the ecological mechanisms driving the responses of many species. An informative parallel approach is to develop size-based methods. These capture the properties of food webs that describe energy flux and production at a particular size, independent of species' ecology. We couple a physical-biogeochemical model with a dynamic, size-based food web model to predict the future effects of climate change on fish biomass and production in 11 large regional shelf seas, with and without fishing effects. Changes in potential fish production are shown to most strongly mirror changes in phytoplankton production. We project declines of 30-60% in potential fish production across some important areas of tropical shelf and upwelling seas, most notably in the eastern Indo-Pacific, the northern Humboldt and the North Canary Current. Conversely, in some areas of the high latitude shelf seas, the production of pelagic predators was projected to increase by 28-89%.

Deep chlorophyll maximum and upper ocean structure interactions: Case of the Guinea Thermal Dome

Deep Chlorophyll Maximum (DCM) modifies the upper ocean heat capture distribution and thus impacts water column temperature and stratification, as well as biogeochemical processes. This energetical role of the DCM is assessed using a 1 m-resolution 1D physical-biogeochemical model of the upper ocean, using climatological forcing conditions of the Guinea Dome (GD). This zone has been chosen among others because a strong and shallow DCM is present all year round. The results show that the DCM warms the seasonal thermocline by +2 degrees C in September/October and causes an increase of heat transfer from below into the mixed layer (ML) by vertical diffusion and entrainment, leading to a ML warming of about 0.3 degrees C in October. In the permanent thermocline, temperature decreases by up to 2 degrees C. The result is a stratification increase of the water column by 0.3 degrees C m(-1) which improves the thermocline realism when compared with observations. At the same time, the heating associated with the DCM is responsible for an increase of nitrate (+300%, 0.024 mu M), chlorophyll (+50%, 0.02 mu g l(-1)) and primary production (+45%: 10 mg C m(-2) day(-1)) in the ML during the entrainment period of October. The considered concentrations are small but this mechanism could be potentially important to give a better explanation of why there is a significant amount of nitrate in the ML. The mechanisms associated with the DCM presence, no matter which temperature or biogeochemical tracers are concerned, are likely to occur in a wide range of tropical or subpolar regions; in these zones a pronounced DCM is present at least episodically at shallow or moderate depths. These results can be generalized to other thermal dome regions where relatively similar physical and biogeochemical structures are encountered. After testing different vertical resolutions (10 m, 5 m, 2.5 m, 1 m and 0.5 m), we show that using at least a 1 to vertical resolution model is mandatory to assess the energetical importance of the DCM.

A Census of Eddy Activities in the South China Sea During 1993-2007

Numerous mesoscale eddies occur each year in the South China Sea (SCS), but their statistical characteristics are still not well documented. A Pacific basin-wide three dimensional physical-biogeochemical model has been developed and the result in the SCS subdomain is used to quantify the eddy activities during the period of 1993-2007. The modeled results are compared with a merged and gridded satellite product of sea level anomaly by using the same eddy identification and tracking method. On average, there are about 32.9 +/- 2.4 eddies predicted by the model and 32.8 +/- 3.4 eddies observed by satellite each year, and about 52% of them are cyclonic eddies. The radius of these eddies ranges from about 46.5 to 223.5 km, with a mean value of 87.4 km. More than 70% of the eddies have a radius smaller than 100 km. The mean area covered by these eddies each year is around 160,170 km(2), equivalent to 9.8% of the SCS area with water depths greater than 1000 m. Linear relationships are found between eddy lifetime and eddy magnitude and between eddy vertical extent and eddy magnitude, showing that strong eddies usually last longer and penetrate deeper than weak ones. Interannual variations in eddy numbers and the total eddy-occupied area indicate that eddy activities in the SCS do not directly correspond to the El Nino-Southern Oscillation events. The wind stress curls are thought to be an important but not the only mechanism of eddy genesis in the SCS.

Modeling Responses of Diatom Productivity and Biogenic Silica Export to Iron Enrichment in the Equatorial Pacific Ocean

Using a three-dimensional physical-biogeochemical model, we have investigated the modeled responses of diatom productivity and biogenic silica export to iron enrichment in the equatorial Pacific, and compared the model simulation with in situ (IronEx II) iron fertilization results. In the eastern equatorial Pacific, an area of 540,000 km(2) was enhanced with iron by changing the photosynthetic efficiency and silicate and nitrogen uptake kinetics of phytoplankton in the model for a period of 20 days. The vertically integrated Chl a and primary production increased by about threefold 5 days after the start of the experiment, similar to that observed in the IronEx II experiment. Diatoms contribute to the initial increase of the total phytoplankton biomass, but decrease sharply after 10 days because of mesozooplankton grazing. The modeled surface nutrients (silicate and nitrate) and TCO(2) anomaly fields, obtained from the difference between the "iron addition'' and "ambient'' (without iron) concentrations, also agreed well with the IronEx II observations. The enriched patch is tracked with an inert tracer similar to the SF6 used in the IronEx II. The modeled depth-time distribution of sinking biogenic silica (BSi) indicates that it would take more than 30 days after iron injection to detect any significant BSi export out of the euphotic zone. Sensitivity studies were performed to establish the importance of fertilized patch size, duration of fertilization, and the role of mesozooplankton grazing. A larger size of the iron patch tends to produce a broader extent and longer-lasting phytoplankton blooms. Longer duration prolongs phytoplankton growth, but higher zooplankton grazing pressure prevents significant phytoplankton biomass accumulation. With the same treatment of iron fertilization in the model, lowering mesozooplankton grazing rate generates much stronger diatom bloom, but it is terminated by Si(OH)(4) limitation after the initial rapid increase. Increasing mesozooplankton grazing rate, the diatom increase due to iron addition stays at minimum level, but small phytoplankton tend to increase. The numerical model experiments demonstrate the value of ecosystem modeling for evaluating the detailed interaction between biogeochemical cycle and iron fertilization in the equatorial Pacific.