Model output for MAREMIP from REcoM2 (REgulated ecosystem model with two phytoplankton classes) coupled to MITgcm. --------------------------------------------------- Model equations, parameters and references: --------------------------------------------------- Model equations including list of model parameters will be uploaded as well or can be found on the web as an attachment to Hauck et al. (2013): Seasonally different carbon flux changes in the Southern Ocean in response to the southern annular mode. GBC, 27, 1236-1245, doi:10.1002/2013GB004600 . Hauck et al (2013) describes the latest model version. Please also cite Schartau et al. (2007): Modelling carbon overconsumption and the formation of extracellular particulate organic carbon, Biogeosciences, 4, 433-454, doi:10.5194/bg-4-433-2007 and Hohn (2009): Coupling and decoupling of biogeochemical cycles in marine ecosystems, PhD Thesis, Universitaet Bremen, http://elib.suub.uni-bremen.de/diss/docs/00011278pdf which I will upload as well. -------------------------------------------------- Atmospheric forcing of the ocean-only model: --------------------------------------------------- The ocean model was forced with MIROC5 atmospheric forcing fields from CMIP5 and was run from 2010 to 2100. -------------------------------------------------- Number of PFTs: ------------------------------------------------- PFTs: 1 zooplankton, 2 phytoplankton: small phytoplankton including a certain percentage of calcifiers (coccolithophores) and diatoms. Variable C:N:Si:Chl ratios -> all pools explicitly modelled for phy (no Si) and dia. C:Fe is fixed. -------------------------------------------------------------- Variables and original grid resolution, regridding -------------------------------------------------------------- Variables were renamed according to the MAREMIP convention (see MAREMIP table), units converted accordingly and variables were regridded onto a 1x1 degree grid (from a resolution of 2deg in longitude, varying between 0.38 to 2deg in latitude. (Higher latitudinal resolution at equator and in Southern Ocean.) Our model configuration is without the Arctic Ocean, so the grid is from 80S to 79N, so not 180 cells, but 159 in latitudinal direction. Name of files are: REcoM_monthly_"varname"_"year".nc with year = 2010 - 2100, varname see list below: Processed variables and their units: 2D fields: MMXL: m PCO2surf: muatm fgCO2: molC/m2/s ICEC: [] ICET: m INTNPPS: molC/m2/s integrated netpp MIX INTNPPD: molC/m2/s integrated netpp DIA 3D fields: SALT: psu THETA: degC UVELMASS: m/s VVELMASS: m/s HPLUS: molH+/m3 DIC: mol/m3 ALK: mol/m3 DIN: mol/m3 SiO3: mol/m3 POC: mol/m3 PON: mol/m3 DOC: mol/m3 DON: mol/m3 FER: mol/m3 BSI: mol/m3 (in detritus) BCaCO3: mol/m3 (in detritus) ZOO: mol/m3 PHYmix: mol/m3 PHYdia: mol/m3 PHYCHLmix: mg/m3 PHYCHLdia: mg/m3 PHYmixN: mol/m3 PHYdiaN: mol/m3 PHYdiaSi: mol/m3 GPPMIX: molC/m3/s GPPDIA: molC/m3/s NPPMIX: molC/m3/s NPPDIA: molC/m3/s REMIN: molC/m3/s gzZOOtot: molC/m3/s gzZOOMIX: molC/m3/s gzZOODIA: molC/m3/s PCACO3: molC/m3/s RDOCZOO: molC/m3/s RPOCZOO: molC/m3/s EXPPOC: molC/m2/s EXPCACO3: molC/m2/s EXPSI: molSi/m2/s EXPPON: molN/m2/s CDIAPOC: molC/m3/s CDIADOC: molC/m3/s CMIXPOC: molC/m3/s CMIXDOC: molC/m3/s DIALIMN: [] DIALIML: [] DIALIMFE: [] DIALIMSI: [] MIXLIMN: [] MIXLIML: [] MIXLIMFE: [] CRESPHYMIX: molC/m3/s CRESPHYDIA: molC/m3/s CRESPHYZOO: molC/m3/s CACO3DISS: molC/m3/s ETOT: W/m2 CRES: molC/m3/s ---------------------------------- NOTES ---------------------------------- NOTE: Export fluxes are not defined at box center, but across upper box limit. Upper box limits are given in these files as z-levels, but uppermost layer is filled with zeros, as this would be the export across the atmosphere-ocean-interface. To get export across 100 m, simply pick level 7 from the 3D export field. This applies for: EXPPOC, EXPCACO3, EXPSI, EXPPON NOTE, definitions on grazing diagnostics: * grazing flux on diatoms by zooplankton (gzZOODIA) grazing_flux_dia * grazing_efficiency grazing on MIX equally, (gzZOOMIX) * total zooplankton grazing (gzZOOtot), grazing efficiency not considered here: grazing_flux_dia + grazing_flux_mix NOTE: I don't think it makes sense to interpolate grid variables, such as grid cell area, volume etc. They would better be calculated from longitude, latitude and depth. If you're interested in the grid variables from our original grid, just let me know. ----------------------------------------------------------------------------- Further derived variables ----------------------------------------------------------------------------- Some variables were not processed, but can be converted from processed variables. This is described in the list below. Please go ahead and do so, if that's what you need. If there is any problem, do not hesitate to contact me and ask for help or for the data. * please extract 2D surface fields from 3D fields yourself. If this is a problem because the files are too big to download on a local computer, please contact me. (Use depth level 7 for EXPORT at 100 m, EXPORT is defined across the given depth levels.) * dpCO2 can be calculated from PCO2surf and the atmospheric CO2 concentration from RCP8.5 * potential density PRHO, calculate from depth, potential temperature and salinity * TOTPHY: add PHYmix + PHYdia * TOTCHL: add PHYCHLmix + PHYCHLdia * ZOO1: we only have one zooplankton, so that's the same as total zooplankton * PHYRCCHL: divide chlorophyll by biomass (convert biomass to mg/m3) * PHYiFE: multiply fixed qFe with PHYmix or PHYdia, iron ratio is fixed. * EXPPOF: multiply fixed qFe with EXPPOC * RCACO3 = PCACO3 * NPP: add NPPDIA + NPPMIX * GPP: add GPPDIA + GPPMIX * pH: convert H+ according to the pH scale that you want to use -------------------------------------------------------------------- Contact details --------------------------------------------------------------------- Judith Hauck, Christoph Voelker, Alfred Wegener Institute - Helmholtz Centre for Polar and Marine Research Am Handelshafen 12 27570 Bremerhaven, Germany judith.hauck@awi.de christoph.voelker@awi.de