Biogeochemical modelling at BSH

lameubiquityMécanique

21 févr. 2014 (il y a 3 années et 1 mois)

155 vue(s)

Biogeochemical modelling at BSH

Frank Janssen, Larissa Müller & Stephan Dick

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

2

What is it good for?

Daily forecasts:

focus on features of
public interest


HABs


oxygen depletion
areas


impact of river
floods


... or whatever might
be of public interest

...

Hindcast/scenario
simulations:

provide information for


HELCOM


OSPAR


Governmental
agencies/decision
makers, e.g. for
Marine Strategy
Framework
Directive
,...


...


25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

3

We are going for an

operational

biogeochemical

model

for the
North Sea/Baltic Sea
region

What are we aiming at?

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

4

How do we (try to) do it?

Stepwise Development:


Circulation model/Dispersion model


SPM module


biogeochemical module


data assimilation


physics (T/S)


biochemistry (Chl
-
a, light, nutrients)

wind, air pressure, air temperature,

cloud coverage, specific humidity


wave data

Circulation Model

(BSHcmod)

for North Sea and Baltic Sea


( 3dim., 5km + 900m grid)


Circulation + Wave Model
(BSHcmod.w)

for North Sea and Baltic Sea


( 3dim., 5km)


Surge Model

(BSHsmod)


for North Sea


(2dim., 5km, barotropic)

Model data archive:
currents, water levels, eddy coefficients,


salinity, temperature, ice data, wave data, meteo. data

Lagrangian Drift and Dispersion

Model
(BSHdmod.L)

for oil, drifting objects and
conservative substances

Eulerian Dispersion

(BSHdmod.E)


for conservative substances,
suspended matter and

biogeochemistry

Local Circulation Models

for estuaries

(Elbe, Weser, Ems)

Other forcing data:

tidal predictions, external surges
(BSHsmod.na), river input (BfG,
SMHI)

Meteorological Models
(GME + COSMO
-
EU)

forecasts up to 84 hours

Wave Model

(WAM)

forecasts up to 84 hrs

surge data

wave data

external surge, tidal constituents,

river runoff

Operational model system at BSH

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

6

Grid nesting:



-

10 km grid


-

5 km grid

-

900 m grid



Operational BSH Model, Version 4

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

7

Surface salinity

900 m grid





eddy formation

06
-
09. April 2008

BSHcmod V4



Cyanobacteria
2006/08/19

Image courtesy of MODIS
Rapid Response Project at
NASA/GSFC

SPM model

(Pleskachevski, Gayer)

u
*
sed
u
*
res
u
*
ero

equilibrium


erosion


sedimentation


resuspension


sedimentation and erosion caused by friction velocity


3 SPM
-
fractions with different sinking
velocities


4 sediment layers (20 cm)

friction velocity u*

depending on :


currents


waves

Transport

Sedimentation

Erosion

Resuspension

Bioturbation

Diffusion

Sand

SPM

Mixing

Sinken






u
*



h
ero

0

< u* <

< u* <

< u* <

< u*

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

9

Berechnete Schwebstoffverteilung
am 03.01.2005

seit 2003

Prä
-
operationeller

Modellbetrieb

BSH Eulerian transport model

with
SPM module of GKSS*

NOAA N
-
16 K1
-
K2


Hindcast (5
km grid)
:

Computed SPM
concentration

at surface on 13.04.2006


* Gayer et al., 2006

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

10

Biogeochemical module


Choose a general interface:
GOTM
-
BIO


Start with „
North Sea model



coupling of BSHdmod.E +
ECOHAM4 (IfM
-
HH)


PhD thesis Larissa Müller, finished in 2007


Implementation of „
Baltic Sea“ model


coupling of BSHdmod.E +
ERGOM (IOW)


Testing both models (and other models, e.g. BFM)


Combined solution for both regions

GOTM:

General Ocean Turbulence Model

www.gotm.net

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

11

Results from BSH biogeochemical model:
North Sea

Calculation based on:


BSHdmod.e


GOTM
-
BIO (ECOHAM4
modul)


two nested grids (3nm,
0.5 nm)

pre
-
operational since
01/2009

http://behemoth.nerc
-
essc.ac.uk/ncWMS/ecoop.html

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

12

Simulated Chl
-
a from ECOHAM4

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

13

Results from BSH biogeochemical model:
Baltic Sea

Calculation based on:


BSHdmod.e


GOTM
-
BIO (ERGOM
modul)


two nested grids (3nm,
0.5 nm)

pre
-
operational since
01/2009

http://behemoth.nerc
-
essc.ac.uk/ncWMS/ecoop.html

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

14

Simulated Chl
-
a from ERGOM

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

15

Data assimilation


In co
-
operation with AWI, DeMarine project


PDAF*) filter framework developed at AWI


Implementation of (local) SEIK filter


SST data from NOAA satellite processed at
BSH

*) Parallel Data Assimilation Framework

http://www.awi.de/de/forschung/wissenschaftliches_rechnen/forschungsthemen_wisse
nschaftliches_rechnen/sequential_data_assimilation/pdaf_parallel_data_assimilation_f
ramework
/

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

16

Stddev of Baltic Sea SST
differences

Application example: Assimilation of Satellite Chlorophyll
data into a Global Ocean
-
Biogeochemical model

mg/m
3

mg/m
3



Model: NASA Ocean


Biogeochemical Model (NOBM)



Observations: Daily chlorophyll


data from SeaWiFS 1997
-
2004



Assimilation method:


Local SEIK filter

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

18

What we have achieved so far ...

ultimate goal
:
fully operational integrated
coupled
physical
-
biogeochemical model system for the North Sea/Baltic Sea
region

steps
: [status from (
-

-

-
) to (





) ]


A) coupling of circulation and wave model (

)


B) coupling of A) with SPM
-
model (



)


C) coupling of an ecosystem model (low resolution, no)
with B) (



)


D) nesting of high resolution model (ku) in C).


physical (





), biogeochemical(





)


E) assimilation of ocean colour data in C, D (
-
)

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

19


Further development of ecosystem models will be done
in the MyOcean project (EU FP7 project , 2009
-
2011)



MyOcean is the implementation project of the GMES
Marine Core Service, aiming at deploying an integrated
pan
-
European capacity for Ocean Monitoring and
Forecasting.


*
GMES (Global Monitoring for Environment and Security) is a European initiative for the implementation


of information services dealing with environment and security.



Future development: MyOcean

MyOcean

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

21

MyOcean
-

Baltic MFC

Co
-
operation of 4 production units

Partner:
DMI, SMHI, BSH and FMI


Consortium approach:


new physical
-
biogechemical model 'HBM'


common model code with central code maintenance


common interface for ecosystem modules


ERGOM, BFM implemented by GOTM interface


different set
-
ups and forcing in different countries


‘ensemble forecasting’

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

22

MFC data
code
Field description
Units
Global
MFC
Arctic
MFC
Baltic
MFC
Med.
MFC
Black Sea
MFC
IBI
MFC
NWS
MFC
M-B-1 or
M-S-1
Potential
Temperature field
Degree
Celsius
V0,V1,
V2
V0,V1,
V2
V0,V1,
V2
V0,V1,
V2
V0 (surface)
V1,V2 (full
depth)
V0,V1
, V2
V0,V1,
V2
M-B-2 or
M-S-2
Salinity field
PSU
V0,V1,
V2
V0,V1,
V2
V0,V1,
V2
V0,V1,
V2
V0 (surface)
V1,V2 (full
depth)
V0,V1
, V2
V0,V1,
V2
M-B-4 or
M-S-4
Zonal velocity
component field
m/sec
V0,V1,
V2
V0,V1,
V2
V0,V1,
V2
V0,V1,
V2
V0 (surface)
V1,V2 (full
depth)
V0,V1
, V2
V0,V1,
V2
M-B-5 or
M-S-5
Meridional
velocity
component field
m/sec
V0,V1,
V2
V0,V1,
V2
V0,V1,
V2
V0,V1,
V2
V0 (surface)
V1,V2 (full
depth)
V0,V1
, V2
V0,V1,
V2
M-B-7 or
M-S-7
Sea level field
m
V0,V1,
V2
V0,V1,
V2
V0,V1,
V2
V0,V1,
V2
V1,V2
V1,
V2
V0,V1,
V2
M-B-8 or
MS-8
Chl-a
mg.m
-3
V2
V2
V2
V0,V1,
V2
V2
V0D,V1
, V2
M-B-8.a or
M-S-8.a
Attenuation
coefficient
m
-1
V2
V2
V1,V2
V0D,V1
, V2
M-B-9 or
M-S-9
Nitrate
mMol N.m
-
3,
V2
V2
V1, V2
V1, V2
V2
V0D,V1
, V2
M-B-9.a or
M-S-9.a
Phosphate
mMol P.m
-
3
V2
V2
V1, V2
V1, V2
V2
V0D,V1
, V2
M-B-10 or
M-S-10
Phytoplancton
biomass
mg C m
-3
V2
V2
V2
V0D,V1
D, V2
M-B-11 or
M-S-11
Dissolved oxygen
mMol O.m
-
3
V2
V2
V0D,V1
D, V2
M-B-12 or
M-S-12
Primary
production
mMolC.m
-
3
d
-1
V2
V2
V0D,V1
D,V2
M-B-16 or
M-S-16
Sea ice
concentration
%
V0,V1,
V2
V0,V1,
V2
V0,V1,
V2
M-B-17 or
M-S-17
Sea ice velocity
field (zonal &
meridonal)
m/sec
V0,V1,
V2
V0,V1,
V2
M-S-18
Sea ice thickness
m
V0, V1,
V2
V0,V1,
V2
V0,V1,
V2
M-S-19
Sea Ice type (FY,
MY)
V1,V2
[j1]
V1,V2
MyOcean FTSS

(Fast Track Service Specification)




Overview of MFC
-
products:

(baseline and standard products of

Monitoring and Forecasting Centres)



V0
-
services: by April 2009 (t0)

V1
-

services : by end of 2010

V2
-

services : in 2011




25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

23

Thank you for your attention!

Estimated Chlorophyll back
-
compared to SeaWiFS



Assimilation strongly improves


surface chlorophyll estimate



Intended deviations (Arabian


Sea, Congo, Amazon)


due to errors in observations



Other deviations in high
-


Chlorophyll regions

mg/m
3

mg/m
3

mg/m
3

GOTM

NPZD

Fasham

ECOHAM

ERGOM

BFM
,

...

BSHdmod.E

passiv

SPM

BIO

forcing data

initialisation data

output data

observations

Post
-
processing

validation

products

archive

Biogeochemical model environment

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

28

ERGOM

NO
3

O
2

P

NH
4

Zoopl.

Detrit.

Fixation

Uptake

Mortality

Grazing

Recycling

Settling

Resuspen
-

sion

Respira
-

tion

Denitrifi
-

cation

Atmosph.

Input

Nitrifica
-

tion

N

2

Solar

Radiation

O

2

P. only

N

2

Flagellates

Cyanobacteria

Diatoms

Sediment

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

29

GOTM
-
BIO: ERGOM

Structure of the
Neumann et al. (2002)
model
with cyanobacteria (cya), diatoms (dia),
dinoflagellates (a), detritus (det),zooplankton
(zoo), ammonium (amm), nitrate (nit) detritus
sediment (sed), oxygen (oxy) and phosphorus
(pho) as the ten state variables.

The concentrations are in mmol N m−3, mmol
N m−2, mmol P m−3 and l O2m−3.
Conservative fluxes are denoted by thin green
arrows, nonconservative

fluxes by bold arrows.

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

30

25.03.2010

WGPBI Aberdeen 2010, Frank
Janssen

31