140
7. Modeling Atmosphere-Ocean Interactions and Primary Productivity
ATMOSPHERE-OCEAN INTERACTION MODEL
Calculation of pC02 in Surface Waters
AlphaO = 1 /(1+(l .28825E-06 /(10 A-pH))+(9.77237El6 /((10 A-pH) A2)))
Alphal = ((pH /lO A-6.3)+l+(10 A-10.3 /pH)) A_1
Alpha2 = (((pH) A2/10 A-l6.6)+(pH /lO A-10.3)+1) A_l
Ct_molar = Surface_Ct /1e5
pC02_Surface = (Surface_Ct /1e5*AlphaO /3.16e-2) *1e6
pH = -1805 .4*Ct_molar+12 .1
C02 Module
Atmospheric_C02(t) = Atmospheric_C02 (t-dt) + (Air_Sea_Exchange) * dt
INIT Atmospheric_C02 = 10.9
OUTFLOWS:
Air_Sea_Exchange
Piston_Velocity*((Atmospheric_C02 /le6)((Surface_Ct/1e5)*AlphaO))
Bottom_Ct(t) = Bottom_Ct(t-dt) + (Downwelled_Carbon +
PartC
Upwelled_CarbonSediment_Rate_of_Carbon) * dt
INIT Bottom_Ct = 7854
INFLOWS:
Downwelled_Carbon = (Surface_Ct /100)*Upwelling_Rate
PartC = (Part_N03*C:N03_Ratio)
Sediment_Remineralization =
Delay(Sediment_Rate_of_Carbon,2,0)
CaC03 = 15
OUTFLOWS:
Upwelled_Carbon = (Bottom_Ct /3400)*Upwelling_Rate
Sediment_Rate_of_Carbon = PartC*.01
Sediment(t) = Sediment(t-dt) +
(Sediment_Rate_of_CarbonSediment_Remineralization) * dt
INIT Sediment = 0.44
INFLOWS:
Sediment_Rate of Carbon = PartC*.Ol