Dynamic Green Ocean Project

Description of some of the input files


Forcing of the biogeochemical model

Deposition of Fe and Si from dust
Dust deposition fields are usually given in g of dust/m2/month. The DGOM PlankTOM5 uses monthly iron deposition values (Tegen & Fung 1995) as input in kg Fe/m2/month (variable dust(i,j)). To convert we assume there are 0.035 g of iron per g of dust. The solubility of Fe in water is taken into account directly in the model.

In the model, we use the following values: 2% solubility for Fe, 7.5% solubility for Si, and 8.8 gSi/gFe. Variables irondep and sidep are in mol/L/time_step for Fe and Si, respectively. The molecular weights of Fe and Si are 55.85 and 28.01, respectively.
Download the dust deposition file on an orthogonal grid with unequal spacing in latitude

Input of nutrients and carbon from rivers
River input of NO3, PO4, SiO3, DIC, Fe (including Fe in DOM), DOC, POC, and POFe (download on the ORCA2 grid). This file is described in da Cunha et al. (2007).

Gas exchange with the atmosphere
Gas exchange of CO2 and O2 is calculated according to Wanninkhof et al (1992), correcting for 100% water vapor according to Sarmiento et al. JGR (1992).


PlankTOM5 Initialisation

The model is initialised with GLODAP observations of DIC, Alkalinity, World Ocean Atlas 2005 NO3, PO4, SiO3, O2.
For the "biological" state variables we use the output of the previous model version:
Download the restart file on an orthogonal grid from TOM5_ET_NCCL9 2005
Download the restart file on the ORCA2 grid
Download restart files on the ORCA2 grid for 1948 from our OpenDAP server.


Physical forcing files
We force the model with NCEP reanalysis. The forcing files have been produced from data downloaded from the NOAA Earth System Research Laboratory. The data have been interpolated to the required grid using the INTERP tool available from Laboratoire d'Océanographie Dynamique et de Climatologie.
  • for each YEAR: file /esdata/env/e031/Input/ncep_bulk_YEAR.nc contains daily values for
    model variable NetCDF variable description units
    tatm air mean Daily Air temperature at 2 m degrees C
    hatm humidity mean Daily Humidity (see below) -
    vatm wspd mean Daily Windspeed at 10 m m/s
    catm cloud daily mean Total cloud cover -
    watm rain mean Daily Precipitation Rate kg/m2/s
  • for each YEAR: file /esdata/env/e031/Input/taux_1d_YEAR.nc
    for each YEAR: file /esdata/env/e031/Input/tauy_1d_YEAR.nc
    contain daily values for
    tauxuflxZonal wind stressN/m2
    tauyvflxMeridional wind stressN/m2
  • averaged : file /esdata/env/e031/Input/runoff_1m_nomask.nc contains monthly values for
    sorunoffsorunofffreshwater input from rivers and icekg/m2/s
In NEMO, humidity is calculated offline:
The humidity values are calculated from the specific humidity at 2m (shum) using the following relationship:
humidity = (shum*pres)/(.622 +.378*shum/sat_p)
where:
pres is the air pressure and sat_p = 611*10(7.5(T-273.15)/(T-35.86))
In OPA8, the forcing files contained the specific humidity at 2m, the air temperature at 2 m (T), and the air pressure, and humidity was calculated online.

Download forcing files from 1948 onwards on the ORCA2 grid from our OpenDAP server.


Initialisation of the physical model
The NEMO GCM is initialised with temperature and salinity from WOA 2005.
The LIM ice model is initialised with ice cover fraction, ice thickness and snow thickness from GIOMAS.
Download restart files on the ORCA2 grid for 1948 from our OpenDAP server.