phosphate and iron concentrations in an ocean carbon cycle ... · phosphate and iron concentrations...

34
Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR), Gokhan Danabasoglu (NCAR), Scott Doney (WHOI), Jun Nishioka (Hokkaido Univ.), Takeshi Yoshimura (CRIEPI), Frank Bryan (NCAR), and Nakashiki Norikazu (CRIEPI)

Upload: others

Post on 17-Jan-2020

0 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Phosphate and iron concentrations in an ocean carbon cycle model

Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR), Gokhan Danabasoglu (NCAR), Scott Doney (WHOI),

Jun Nishioka (Hokkaido Univ.), Takeshi Yoshimura (CRIEPI), Frank Bryan (NCAR), and Nakashiki Norikazu (CRIEPI)

Page 2: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Objects

• To represent realistic global phosphate and iron distribution in comparison with observational database in an general circulation model

• To understand the reasonable parameters in a model for the current distribution of phosphate and iron

Page 3: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Physical model• CCSM3 POP (Keihl and Gent, 2004)• Horizontal; about 3 degree, Vertical; 25 layers• 3rd order upwind advection scheme• KPP scheme (Large et al., 1994)• Momentum; anisotropic GM scheme (Smith and

McWilliams, 2003)• Tracer; GM scheme (Gent and McWilliams,

1990)• Normal Year Forcing (Large and Yeager, 2005)• Lax and Wendroff advection scheme for tracer• Near-boundary eddy flux parameterization

(Danabasoglu, et al., 2006)

Page 4: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

CO2, O2 Air Sea Gas exchange Fe deposit ion

PO4, Fe

POCCaCO3Fe (Jprod)

DOC DOP, DOFe

Fe scavenging

75 m

DIC, O2, ALK(Ratio to PO4)

Remineralization

DOP PO4

Jprod

τ⎭⎬⎫

⎩⎨⎧

⋅⋅⋅⋅ cINT

zz

BFFF,1max( ) prodJσ−1

DOP

DOPτ

PO4 cycle

prodJσ

Remineralization

PO4DOP

DICALKO2

FeDOFe

Sinking

Productionzone

Consumptionzone

Schematic of ocean carbon cycle model (OCMIP’)

Calculation byUptake ratio x : PO4

Prognostic PO4 cycle

Page 5: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Parameter for PO4• Equation for PO4 cycle

ddt

PO4 − ∇D∇PO4 = −Jprod +DOPτDOP

+∂∂z

FPOP

Jprod = FT ⋅ FN ⋅ FI ⋅ B ⋅max 1, zml

zc

⎧ ⎨ ⎩

⎫ ⎬ ⎭

τ

FPOP =zzc

⎝ ⎜

⎠ ⎟

β

σ Jprod−zc

0∫ x,y,z,t( )dz

Production time scaleτ= 15, 5,1

FT ; Temperature limitation functionFN ; Nutrient limitation function, Michaelis-Menten limiting term, PO4 or FeFI ; Light (Irradiance) limitation functionB ; Proxy for biomass , PO4 or FeZml ; Mixed layer depthZc ; Compensation depth, 75m

Page 6: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Parameter for PO4

5.85 x 10-4Uptake ratio rFe:PO4

Iron concentration

Control

Variable uptake ratio

7.00 x 10-4

3.00 x 10-4

2 x Half Sat. const.6/7 x Half Sat. const.

• Uptake ratio Fe:PO4

Page 7: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Variable uptake ratio Fe : P

71m42m

21m4m

Uptake ratio change with iron concentrationsConstant uptake ratio is about 6.0e-4

Page 8: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Observation

Variable Uptake ratio

τ=5

τ=1

τ=15

τ=5

Constant Uptake ratio

PO4 surface concentrations

Page 9: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Variable Uptake ratio

τ=5

τ=1

τ=15

τ=5

Constant Uptake ratio

POP export fluxat 100m depth

Page 10: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Variable Uptake ratioτ=5τ=5

Constant Uptake ratio

Iron surface concentrations

Page 11: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Iron cycle (Doney et al., 2006)

DOFe

RFe:P ⋅ FT ⋅ FN ⋅ FI ⋅ B ⋅max 1, z

zc

⎧ ⎨ ⎩

⎫ ⎬ ⎭

τ

( ) prodJσ−1

DOFeτDOFe

prodJσJprod

RemineralizationFe

FeL Fe’

Complexation

POFe

Scavenging

Bottom

Atmospheric deposition

Page 12: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Iron scavenging (Archer and Johnson,

2000)

• Fefree ; Fe that is not bound to ligands• L ; concentration of ligand• KL ; Strength of the binding reaction• L = 1.0 nmol/L, KL = 300 L/nmol

Fe2free + L +

1KL

− Fe⎛

⎝ ⎜

⎠ ⎟ Fefree −

FeKL

= 0

Fwscav = FeFree ⋅ C0 ⋅ 1+ α exp −z

zscav

⎝ ⎜

⎠ ⎟

⎝ ⎜

⎠ ⎟

• C0=0.2 y-1, α=200, zscav = 250m• Scavenging Fe is attached to the sinking particles to form POFe• A fraction (40%) of the scavenged Fe is assumed to be insoluble

and is directly lost to the sediments.• The remaining 60% can be remineralized back to dissolved form

below zc

Page 13: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Iron supply from atmospheric dust

• Monthly data (Mahowald et al., 2003) • 3.5% Fe by weight with 2% of the Fe bioavailable (Fung

et al., 2000)。

Annual mean dust flux (kg/m2/s)

Page 14: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Surface distribution of iron

Page 15: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Iron observation5500 samples Based on Parekh et al., 2005 (Moore et al., in preparation)

Page 16: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Iron observationRemove coastal zoneWestern NP is higher than eastern NPSO is lowerIndian Ocean is higherAtlantic Ocean is not so higher

Page 17: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Vertical distribution of iron

Cal.Obs.

Western NP is higher than eastern NPSO is lowerIndian Ocean is higherAtlantic Ocean is not so higher

Modeling of scavenging with complexation is important for reasonable vartical profile.Not different between basins.

Page 18: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Surface distribution of iron

Parekh et al., 2005

Observation Calculation

Page 19: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Iron distribution at the depth of 1000m

Parekh et al., 2005

Observation Calculation

Page 20: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Parameter of iron dust deposition

• 3.5% Fe by weight with 2% of the Fe bioavailable (Fung et al., 2000)。

1% 2% 3%

Page 21: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Iron concentration dust deposition

1%

3%2%Parekh et al., 2005

Page 22: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Iron concentration dust deposition

1%

3%2%Parekh et al., 2005

Fe concentration at 1000m

Page 23: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

PO4 concentration dust deposition

1%

3%2%

Obs.

Page 24: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Scavenging parameters

Fwscav = FeFree ⋅ C0 ⋅ 1+ α exp −z

zscav

⎝ ⎜

⎠ ⎟

⎝ ⎜

⎠ ⎟

C0=1.0 y-1

α=9zscav = 500m

C0=0.5 y-1

α=100zscav = 250m

C0=0.2 y-1

α=200zscav = 250m

Page 25: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Scavenging parametersFe concentration at 1000m

C0=1.0 y-1 α=9 zscav = 500m

C0=0.5 y-1 α=100 zscav =250m

C0=0.2 y-1 α=200 zscav = 250m

Parekh et al., 2005

Page 26: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Conclusions• Variable uptake ratio provides the reasonable distribution

of phosphate and dose not affect the iron distribution.• Iron scavenging model with complexation provides the

reasonable vertical distribution of iron.• When only iron input from atmosphere is considered, the

difference of iron concentration between basins can not be represented in this model.

• Iron dust deposition rate affect the iron distribution but it is difficult to understand which is better rate due to still sparse observational database of iron concentrations.

• We acquired the reasonable phosphate and iron distribution for future numerical experiments to understand the effect of input from Amour river via Okhotck ocean to iron concentration in the western North Pacific ,and artificial iron fertilization experiments.

Page 27: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

end

Page 28: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Sinking flux to the sea bottom ZbArcher and Johnson (2000)

( ) ( )( )( )min6.0 rescavtb

tb

POFeFezzzPOFezPOFe−−+

=

POFeremin =POFe zt( )− POFe zb( )

zb − zt

= POFe zt( )1− (zb zt )−β

zb − zt

Page 29: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Vertical distribution of iron

Observation Calculation

Page 30: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Long term trend

Temperature Salinity

Page 31: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Iron concentration dust deposition

1%

3%2%

Page 32: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Scavenging parametersC0=1.0 y-1

α=9zscav = 500m

C0=0.5 y-1

α=100zscav = 250m

C0=0.2 y-1

α=200zscav = 250m

Fwscav = FeFree ⋅ C0 ⋅ 1+ α exp −z

zscav

⎝ ⎜

⎠ ⎟

⎝ ⎜

⎠ ⎟

Page 33: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

Scavenging parameters

Parekh et al., 2005

Page 34: Phosphate and iron concentrations in an ocean carbon cycle ... · Phosphate and iron concentrations in an ocean carbon cycle model Daisuke Tsumune (CRIEPI), Keith Lindsay (NCAR),

PO4 concentration scavenging paremeter

Obs.