= Coupled vegetation - atmospheric chemistry model = == Context == [[Image(Diapositive1.GIF, 40%)]], [[Image(Diapositive1.2.GIF, 40%)]] The terrestrial biosphere is a significant source of reactive chemical compounds, such as isoprene, monoterpenes..., or other volatile organic compounds emitted by plants and mainly leaves, or nitrogen oxides emitted by nitrification and denitrification soil processes. Those compounds play a key role in the atmospheric chemical composition and especially in the ozone cycle or aerosol formation. Moreover, distribution and growth of terrestrial ecosystems is a strong driver of the deposition, and therefore the loss, of atmospheric chemical compounds at the surface. On the other hand, strong and persistent atmospheric concentrations of compounds such as ozone, nitrogen oxides or sulfur doixide can lead to leaf necrosis, decrease in photosynthetic activity and alter plan growth. There are therefore strong interactions between the terrestrial biosphere and the atmospheric chemical composition, that have to be considered when investigating the evolution of the Earth System. This objective of the coupling between the vegetation model ORCHIDEE and the chemistry-transport model LMDzINCA is to take into account those interactions and add consistency in the IPSL Earth System Model when related to the description of vegetation growth and distribution: * Biogenic emissions of reactive compounds (VOCs and NOx) calculated by ORCHIDEE sent to INCA, instead of having an input emission files * Distribution of vegetation types prescribed or calculated by ORCHIDEE sent to INCA, corresponding to the time period of interest, instead of having one vegetation distribution file only, read as input file by INCA for any scenario * In the purpose of investigating the impact of atmospheric chemical composition on plant growth and distribution: atmospheric concentrations of ozone (in a first step) and other chemical compounds (longer-term) calculated by INCA to be sent to ORCHIDEE where the impact on plant functional types will be calculated (see the PhD work carried out by Thomas Verbeke, LSCE, see this [wiki:Branches/ORCHIDEE-OzoneImpact page]). [[PageOutline]] == Technical notes == === in orchidee === '''SVN VERSION''' : branches/ORCHIDEE-BCOV [[BR]] * Add a routine '''sechiba_GetFromOrchideeToInca''' in sechiba.f90. This routine is called from INCA. * output :: nvm_inca --> Number of vegetation types * output :: veget_max_inca --> Max. fraction of vegetation type (LAI -> infty) * input :: field_out_names --> Names for emission variables. Inca give the names to Orchidee * output :: fields_out --> fields for emissions variables : to be sent by Orchidee to Inca * To be done : an other routine '''sechiba_GetFromIncaToOrchidee'''. This routine will be called from ORCHIDEE. And will ask for variables (ozone, and other compounds, atmospheric concentrations) from Inca. * Modify declaration of flux in diffuco. === in inca === '''SVN VERSION ''' : branches/INCA4_BCOV * Add a routine '''surf_chem_atm'''. This routine call sechiba_GetFromOrchideeToInca * Add a module for the variables management '''surf_chem_mod'''. This module use getin to read the flux names in inca.def parameter file. * Add output file veget.nc writing all flux and veget_max. * (08 august 2012) add title and unit for flux in output file veget.nc {{{ flx_iso: Isoprene emissions from vegetation, kgC/m²/s flx_mono: Monoterpene emissions from vegetation, kgC/m²/s flx_ORVOC: Other Volatile Organic Compound emissions from vegetation, kgC/m²/s flx_MBO: 2-methyl-3-buten-2-ol emissions from vegetation (mainly pines in America), kgC/m²/s flx_methanol: Methanol emissions from vegetation, kgC/m²/s flx_acetone: Acetone emissions from vegetation, kgC/m²/s flx_formal: Formaldehyde emissions from vegetation, kgC/m²/s flx_acetal: Acetaldehyde emissions from vegetation, kgC/m²/s flx_acetic: Acetic acid emissions from vegetation, kgC/m²/s flx_formic: Formic acid emissions from vegetation, kgC/m²/s flx_no_soil: Nitrogen Oxide emissions from soil, before deposition on canopy, ngN/m²/s flx_no: Net nitrogen Oxide emissions from soil, after deposition on canopy, ngN/m²/s flx_fertil_no: Nitrogen Oxide emission related ONLY to the use of fertilisers, before deposition on canopy, ngN/m²/s veget_max : SurfType_frac, 1 }}} === in config === * Add a parameter file inca.def {{{ #nombre de flux transferes depuis orchidee vers inca nbFlux_FromOrch=13 #noms des flus a transferer # iso - mono - ORVOC - MBO - methanol # acetone - acetal - formal - acetic - formic # no_soil - no - fertil_no emi_FromOrch=iso mono ORVOC MBO methanol acetone acetal formal acetic formic no_soil no fertil_no }}} * Add ''DIFFUCO_OK_INCA=y'' in orchidee.def === Extract this new configuration === For this you need to add in mod.def {{{ #-H- LMDZORINCA_BCOV LMDZ4 with ORCHIDEE and INCA (coupled vegetation #- and chemistry): #-H- LMDZORINCA_BCOV ORCHIDEE branches BCOV (close to trunk 951) #-H- LMDZORINCA_BCOV INCA branches BCOV (close to trunk 264) #-H- LMDZORINCA_BCOV LMDZ5 LMDZ5/trunk rev 1628 #-H- LMDZORINCA_BCOV IOIPSL/src svn tags/v2_2_1 #-H- LMDZORINCA_BCOV libIGCM tag libIGCM_v2.0_beta4 #-M- LMDZORINCA_BCOV Anne.Cozic@lsce.ipsl.fr #-C- LMDZORINCA_BCOV IOIPSL/trunk/src HEAD 8 IOIPSL/src modeles #-C- LMDZORINCA_BCOV branches/ORCHIDEE-BCOV/ORCHIDEE HEAD 14 ORCHIDEE modeles #-C- LMDZORINCA_BCOV LMDZ5/trunk 1628 11 LMDZ modeles #-C- LMDZORINCA_BCOV branches/INCA4_BCOV HEAD 9 INCA3 modeles #-C- LMDZORINCA_BCOV tags/libIGCM_v2.0_beta4 HEAD 10 libIGCM . #-C- LMDZORINCA_BCOV CONFIG/LMDZORINCA/branches/LMDZORINCA_BCOV HEAD 8 LMDZORINCA_v5 config }}} And compile with the resolution ''NMHC_AERxLMD9695-L39'' === Simulations === * first test : one year climatic with this new configuration (LMDZ5 1628 - Orchidee-BCOV - INCA4_BCOV). /ccc/store/cont003/dsm/p86cozic/IGCM_OUT/LMDZORINCA/NMHC_AER/PROD/BCOV/L4OI439.4