![Page 1: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/1.jpg)
MAMAMalta meeting, 27-30 January 2004
Expert Meeting
Towards Operational ecological models in coastal areas [email protected]
![Page 2: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/2.jpg)
The pelagic physical-biological interactions in the ocean
B
C
D
Stratification Mixing
lightlimitatio
n
New production
Regenerated production
Oceanic Ecosystems
Coastal
Ecosystems
Flagellatesand bacteria
Large phytoplankton
Microbial food web
Herbivorous food web
3
2
5
1
4
1 2
3
45
A
F
E
Legendre and Rassoulzadegan, 1995
Nutrient limitatio
n
![Page 3: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/3.jpg)
Ocean ecosystem dynamics strongly coupled with Ocean dynamics
Factors limiting predictability:
DataPredictability of the atmospheric forcing (coastal areas).Predictability of external inputs (River runoff and nutrient load)
ModelOpen boundary condition (Limited area nested models)Definition of initial conditions for forecast simulationsInitial adjustment problem for nested models.
To overcome (or reduce) such problems, the forecasting System must encompass both the open and the coastalOcean scales……
![Page 4: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/4.jpg)
The components of an interdisciplinary forecasting system
![Page 5: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/5.jpg)
Buoy stations
Adricosm“in situ”ObservingSystemCurrentlyRunning
![Page 6: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/6.jpg)
Adricosm remoteObserving System
SeaWifs
AVHRR
TOPEX
ERS-2
![Page 7: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/7.jpg)
The coupled physica-ecological modelling system
Need - Water column and sediment prognostic equations for
Physical state variablesMacro-scale: T, S, ρ, p, u, v, w (equation of motion equation of state equations for scalar properties conservation)
Sub-grid scale: Kv, KH, Iz (turbulence closure equations radiative transfer equations)
Air-sea fluxes:τw, Q, (E-P) (bulk formulae)
Water sediment interactions: τb, (bulk formulae)
![Page 8: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/8.jpg)
The coupled physical-ecological modelling system
Need - Water column and sediment prognostic equations for
chemical state variablesC, N, P, Si, (equation for non conservative scalar properties)
biological state variables (Functional groups): Phytoplankton, bacteria, zooplankton etcEach organism can be described by a 4D vectorVj=[VC, VN,VP,VSi] Where the subscripts C, N, P, Si are the “chemical currencies” or concentrations of chemical consituents in each organism
Basic Assumption: The dynamicsof the marine ecosystem can be expressed by the dynamics of the j-th element in each functional group V (biomass based model):
![Page 9: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/9.jpg)
Organism (C:N:P)organism
CO2
Basal activityStress respiration
Food components(C:N:P)food
Uptake Predation
Predators(C:N:P)food
Detritus fractions
MortalityExcretion Defaecation
Nutr.Nutrientexcretion
The “Standard Organism” (Functional group approach)
![Page 10: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/10.jpg)
Thus, the fundamental structure ofthe marine ecosystemModel Is:
1. Physical environment description (macro and micro-scales)2. Chemical currencies3. Functional groups (Different species in a single group)4. Closure hypothesis(or individual based modelling) for Higher trophic levels.
All components interacting in a deterministicway with bulk parameterizations
![Page 11: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/11.jpg)
The pelagic component of the MFSTEP Biogeochemical Fluxes Model
![Page 12: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/12.jpg)
The benthic component of the MFSTEP Biogeochemical Fluxes Model
![Page 13: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/13.jpg)
Mathematical formulation
biolphys
j
tV
tV
t
V
jHHHj
Vj
SjHH
phys
j VKz
VK
zz
VwwVU
t
V
jk
N
1kbiol
j Ft
V
Where N are the number of theBiogeochemical interactions forEach functional group
![Page 14: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/14.jpg)
EcologyPelagic Model
EcologyBenthic Model
CirculationModel
T (x, y, z, t)
S (x, y, z, t)
KH (x, y, z, t)
A (x, y, z, t)
u, v, w (x, y, z, t)
Nutrient inputParticulate
Inorganic Matter Qs Qb+Qe+Qhw (E-P-R)
PAR
Sedimentary andWater-Sediment
diffusive processes
THE GENERAL STRUCTURE OF THE MODELS FORCING AND COUPLING
TransportModel
Cp (x, y, z, t)
NumericalDriver
(Time Integration)
bio
p
t
C
phys
p
t
C
t
Cp
t
Cb
![Page 15: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/15.jpg)
Implementation towards operational use of ecological models
MFS strategy:
• Implementation of 1D models in data rich areas to validate/calibrate models and check the physical/ biological coupling (MFSPP task accomplished)
• Extend the implementation to 3D with climatological forcing and nesting approach (MFSTEP task underway) • Explore the use of data assimilation schemes for biogechemical state variables (MFSTEP task underway)
![Page 16: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/16.jpg)
1D implementations: Validation
Observed Seasonal Inorganic Suspended Matter Profiles(forcing functions in the light attenuation processes)
Chlorophyll Phosphate
![Page 17: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/17.jpg)
1D implementations: Validation under high frequency forcingBacterial biomass: 48 h simulation with 6hr atmospheric forcing
Observations Model
![Page 18: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/18.jpg)
S1
AA1 S3
Critical Depth
ML Depth
Chl-a(Cd ave.)
1DImplementations physical/ecological interactions:the Sverdrup-like mechanism
![Page 19: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/19.jpg)
O Data + stdev Standard model Improved model
Comparison with observedBacterial Carbon Production (BCP) rates
BCP = -b*f(T)*B +(1-BGE)*U(substrate)
BGE = 0.3 (standard)
BGE = c – a*T(Rivkin and Legendre, 2001)
1DImplementationimproving biological processes
![Page 20: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/20.jpg)
3D implementations: Nested approach based on MFSPPCirculation modelling
OGCMCoupled Model
RegionalCoupled Models
The MFSTEP Coupled Models Domain
The MFSTEP Coupled Models Domain
![Page 21: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/21.jpg)
The Adriatic modelling systemBased on the Princeton Ocean Model (POM)
And the Modular Ocean Model (MOM)
AIM (Adriatic Intermediate Model)POMWhole Adriatic Sea.5 km horizontal resolution, 21 sigma layersNested with the Mediterranean SeaGeneral Circulation Model
NASM (Northern Adriatic Shelf Model)POMNorthern Adriatic only1.5 km horizontal resolution11 Sigma layersNested with AIM
Mediterranean Sea OGCM (MOM)1/8° Horizontal resolution31 levels
![Page 22: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/22.jpg)
Preliminary results forthe Adriatic
Chlorophyll-a
![Page 23: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/23.jpg)
Surface DOC distribution
mgC/m3
winter
10 days 20 days
30 days
![Page 24: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/24.jpg)
Testing data assimilation schemes:The Singular evolutive Kalman Filter (Triantafyllou et al.2003)
![Page 25: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/25.jpg)
Testing data assimilation schemes:The Singular evolutive Kalman Filter (Triantafyllou et al.2003)
![Page 26: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/26.jpg)
Testing data assimilation schemes:The Singular evolutive Kalman Filter (Triantafyllou et al.2003)
![Page 27: MAMA Malta meeting, 27-30 January 2004 Expert Meeting](https://reader035.vdocument.in/reader035/viewer/2022081603/568138a0550346895da05cc9/html5/thumbnails/27.jpg)
CONCLUSIONS
• Operational ecological modelling lags (naturally) behind
operational circulation modelling• The nested modelling approach can potentially face
the problem of capturing and describing the many spatial and temporal scales manifested in marine ecosystem dynamics
• Potential for predictions is apparent• Data assimilation schemes can be successfully used