|Numerical modeling of the central Black Sea ecosystem functioning during the eutrophication phase|Grégoire, M.; Raick, C.; Soetaert, K. (2008). Numerical modeling of the central Black Sea ecosystem functioning during the eutrophication phase. Prog. Oceanogr. 76(3): 286-333. dx.doi.org/10.1016/j.pocean.2008.01.002
In: Progress in Oceanography. Pergamon: Oxford,New York,. ISSN 0079-6611; e-ISSN 1873-4472, meer
ecosystem-hydrodynamic interactions; parameter estimation; mathematical
|Auteurs|| || Top |
- Grégoire, M.
- Raick, C.
- Soetaert, K.
A one-dimensional coupled physical-biogeochemical model has been developed to simulate the ecosystem of the central Black Sea at the end of the 1980s when eutrophication and invasion by gelatinous organisms seriously affected the stability and dynamics of the system. The physical model is the General Ocean Turbulence Model (GOTM) and the biogeochemical model describes the foodweb, from bacteria to gelatinous carnivores through 24 state variables including three groups of phytoplankton: diatoms, small phototrophic flagellates and dinoflagellates, two zooplankton groups: micro- and mesozooplankton, two groups of gelatinous zooplankton: the omnivorous and carnivorous forms, an explicit representation of the bacterial loop: bacteria, labile and semi-labile dissolved organic matter, particulate organic matter. The model simulates oxygen, nitrogen, silicate and carbon cycling. In addition, an innovation of this model is that it explicitly represents processes in the anoxic layer. Biogeochemical processes in anaerobic conditions have been represented using an approach similar to that used in the modeling of diagenetic processes in the sediments lumping together all the reduced substances in one state variable [Soetaert, K., Herman, P., 1996. A model of early diagenetic processes from the shelf to abyssal depths. Geo-chimica et Cosmochimica Acta 60 (6) 1019-1040]. In this way, processes in the upper oxygenated layer are fully coupled with anaerobic processes in the deep waters, allowing to perform longterm simulations.
The mathematical modeling of phytoplankton and zooplankton dynamics, detritus and the microbial loop is based on the model developed by Van den Meersche et a]. [Van den Meersche, K., Middelburg, J., Soetaert, K., van Rijswijk P.H.B., Heip, C., 2004. Carbon-nitrogen coupling and algal-bacterial interactions during an experimental bloom: Modeling a 13c tracer experiment. Limnology and Oceanography 49 (3), 862-878] and tested in the modeling of mesocosm experiments and of the Ligurian sea ecosystem [Raick, C., Dethez, E., Soetaert, K., Gregoire, M., 2005. Study of the seasonal cycle of the biogeochemical processes in the Ligurian sea using an ID interdisciplinary model. Journal of Marine Systems 55 (3-4) 177-203]. This model has been extended to simulate the development of top predators, the aggregation of detritus as well as the degradation and chemical processes in suboxic/anoxic conditions (e.g. denitrification, anoxic remineralization, redox reactions).
The coupled model extends down to the sediments (~= to 2000 m depth) and is forced at the air-sea interface by the 6 hourly ERA-40 reanalysis of ECMWF data. The model has been calibrated and validated using a large set of data available in the Black Sea TU Ocean Base. The biogeochemical model involves some hundred parameters which are first calibrated by hand using published values. Then, an identifiability analysis has been performed in order to determine a subset of identifiable parameters (i.e. ensemble of parameters that can be together estimated from the amount of data we have at our disposal, see later in the text). Also a subset of 10 identifiable parameters was isolated and an automatic calibration subroutine (Levenberg Marquart) has been used to fine tune these parameters. Additionally, in order to assess the sensitivity of model results to the parameterization of the two gelatinous groups, Monte Carlo simulations were performed perturbing all the parameters governing their dynamics.
In order to calibrate the particle dynamics and export, the chemical model was run off-line with the particle and microbial loop model in order to check its capacity of simulating anoxic waters. After a 104 year run, the model simulated NH4 and H2S profiles similar to observations but steady state was not reached suggesting that the Black Sea deep waters are not at steady state. The fully coupled model was then used to simulate the perio