Changes between Version 173 and Version 174 of DevelopmentActivities/ORCHIDEE-DOFOCO


Ignore:
Timestamp:
2018-11-28T22:48:04+01:00 (5 years ago)
Author:
luyssaert
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • DevelopmentActivities/ORCHIDEE-DOFOCO

    v173 v174  
    403403 
    404404Table 2. Essential technical capabilities before evaluating ORCHIDEE-CN-CAN. 
    405 || '''Description of the task'''' || '''Status''''|| 
     405|| '''Description of the task''' || '''Status''' || 
    406406|| Is the model stable? Can it be run for 500 - 1000 years? Can it be run for all PFTs? || Recent problems with grasslands || 
    407407|| Is the N-cycle working? Can we run simulations with IMPOSE_CN=n and IMPOSE_CN=y?     || OK || 
     
    419419Observational data products for model-data comparison can be found at:  /home/satellites5/maignan/ORCHIDEE/EVALUATION. 
    420420* Phenology  
    421 ◦       GIMMS (LAI and FPAR3g) 
     421GIMMS (LAI and FPAR3g) 
    422422* Forest structure 
    423423◦       Remote sensing products of biomass (temperate and boreal maps, i.e., Turner) 
     
    463463|| Extract Jung GPP and compare against spatially explicit simulations for 8 regions in Europe  || OK || dofoco/dofoco/ SCRIPTS/XXX || 
    464464|| Extract Jung evapotranspiration and compare against spatially explicit simulations for 8 regions in Europe || OK || dofoco/dofoco/ SCRIPTS/XXX || 
    465 || Extract MODIS albedo (NIR & VIS) and compare against spatially explicit simulations for 8 regions in Europe || In progress || dofoco/dofoco/ SCRIPTS/XXX || Extract height and compare against spatially explicit simulations for 8 regions in Europe ||     OK || dofoco/dofoco/SCRIPTS/XXX || 
     465|| Extract MODIS albedo (NIR & VIS) and compare against spatially explicit simulations for 8 regions in Europe || In progress || dofoco/dofoco/ SCRIPTS/XXX ||  
     466|| Extract height and compare against spatially explicit simulations for 8 regions in Europe || OK || dofoco/dofoco/SCRIPTS/XXX || 
    466467|| Extract effective LAI and compare against spatially explicit simulations for 8 regions in Europe || OK || dofoco/dofoco/ SCRIPTS/XXX || 
    467468|| Extract BA and compare against spatially explicit simulations for 8 regions in Europe || OK || dofoco/dofoco/ SCRIPTS/XXX || 
     
    470471=== Parameterization === 
    471472As we will compare the simulations to observations, the simulations need to be our best shot to resemble reality. Hence, all site-level simulations made for the parameterization will have the following configuration: 
    472 •       start from a spin-up 
    473 •       impose_cn = no (i.e., accounting for N-deposition ) 
    474 •       Forest management = 2 (i.e., accounting for a thin and fell type of forest management) 
     473* Start from a spin-up 
     474* impose_cn = no (i.e., accounting for N-deposition ) 
     475* Forest management = 2 (i.e., accounting for a thin and fell type of forest management) 
    475476 
    476477European simulations will have the following configuration: 
    477 •       start from a spin-up 
    478 •       impose_cn = no (i.e., accounting for N-deposition ) 
    479 •       impose_veg = no (i.e., accounting for PFT distribution) 
    480 •       Forest management from map 
     478* Start from a spin-up 
     479* impose_cn = no (i.e., accounting for N-deposition ) 
     480* impose_veg = no (i.e., accounting for PFT distribution) 
     481* Forest management from map 
    481482 
    482483The following parameterization approach – making use of parameters that were already shown to be sensitive to tuning – is proposed: 
    483484 
    484485Table 4. Proposed order, parameters and tools to parameterize ORCHIDEE-CN-CAN. Once all steps has been performed, the NPP to GPP ratio should be re-evaluated, possibly resulting in another tuning cycle for some of the PFTs. 
    485 Process Parameter(s)    Tool 
    486 Onset of growing season and start of senescence thresholds for phenology        FLUXNET 
    487 NPP/GPP ratio   coeff for maintenance respiration       FLUXNET + Campioli 
    488 Magnitude of LAI        k_latosa_min and k_latosa_max   NFI + Luyssaert et al 2007 
    489 Magnitude of GPP        LL_alpha, Vcmax, J_max  FLUXNET 
    490 Magnitude of NPP        Implicit through NPP/GPP and GPP        Luyssaert et al 2007 
    491 Evapotranspiration      LAI_top, water stress   FLUXNET 
    492 Water stress    To be discovered        FLUXNET 
    493 Magnitude of Rh To be checked - Rh      Luyssaert et al 2007 
    494 Magnitude of TER        To be checked - Rh      FLUXNET 
    495 Diameter, height, and biomass   form factor     NFI + Luyssaert et al 2007 
    496 density and biomass     self-thinning parameters        NFI + Luyssaert et al 2007 
    497 Harvest and biomass     Self-thinning and RDI parameters        NFI + ??? 
    498 Albedo  Implicit through LAI and forest structure       FLUXNET? 
    499 Tree ring width Self-thinning, recruitment      ITRDB or VERIFY Fig 4 
     486|| '''Process''' || '''Parameter(s)''' || '''Tool''' || 
     487|| Onset of growing season and start of senescence || thresholds for phenology || FLUXNET || 
     488|| NPP/GPP ratio || coeff for maintenance respiration || FLUXNET + Campioli || 
     489|| Magnitude of LAI || k_latosa_min and k_latosa_max || NFI + Luyssaert et al 2007 || 
     490|| Magnitude of GPP || LL_alpha, Vcmax, J_max || FLUXNET || 
     491|| Magnitude of NPP || Implicit through NPP/GPP and GPP || Luyssaert et al 2007 || 
     492|| Evapotranspiration || LAI_top, water stress || FLUXNET || 
     493|| Water stress || To be discovered || FLUXNET || 
     494|| Magnitude of Rh || To be checked - Rh || Luyssaert et al 2007 || 
     495|| Magnitude of TER || To be checked - Rh || FLUXNET || 
     496|| Diameter, height, and biomass || form factor || NFI + Luyssaert et al 2007 || 
     497|| Density and biomass || self-thinning parameters || NFI + Luyssaert et al 2007 || 
     498|| Harvest and biomass || Self-thinning and RDI parameters || NFI + ??? || 
     499|| Albedo || Implicit through LAI and forest structure || FLUXNET? || 
     500|| Tree ring width || Self-thinning, recruitment || ITRDB or VERIFY Fig 4 || 
    500501 
    501502 
     
    517518 
    518519Longitudinal bands can be run largely independent from the 1-pixel tests. After concluding the 1-pixels tests it probably makes sense to run the longitudinal bands to check whether the latest parameters still produce the expected response gradients. 
    519  
    520520         
    521521==== Settings for European tests ==== 
     
    538538•       N-deposition, we have N-deposition files from approximately 1860 until present. 
    539539 
    540 Table 4. Observational data products to evaluate European control simulation 
    541 Observations    Tool 
    542 EU-wide NFI data (Fig. 3) through the VERIFY project    ??? 
    543 Jung GPP – relies on FLUXNET thus partly circular       XXX 
    544 Jung evapotranspiration – relies on FLUXNET thus partly circular        XXX 
    545 MODIS albedo NIR & VIS  XXX 
    546 Remote sensing tree height      XXX 
    547 Extract effective LAI from MODIS        XXX 
    548 BA from JRC – relies on NFI thus partly circular        XXX 
    549 Biomass from JRC – relies on NFI thus partly circular   ??? 
    550 Biomass from Turner – Remote sensing?   ??? 
    551 Harvest from Schelhaas - data through the VERIFY project        ??? 
    552 LCC – Luyssaert et al 2014      XXX 
    553 LCC – Duveiller et al 2018      ??? 
    554  
     540Table 5. Observational data products to evaluate European control simulation 
     541|| '''Observations''' || '''Tool''' || 
     542|| EU-wide NFI data (Fig. 3) through the VERIFY project || ??? || 
     543|| Jung GPP – relies on FLUXNET thus partly circular || XXX || 
     544|| Jung evapotranspiration – relies on FLUXNET thus partly circular || XXX || 
     545|| MODIS albedo NIR & VIS || XXX || 
     546|| Remote sensing tree height || XXX || 
     547|| Extract effective LAI from MODIS || XXX || 
     548|| BA from JRC – relies on NFI thus partly circular || XXX || 
     549|| Biomass from JRC – relies on NFI thus partly circular || ??? || 
     550|| Biomass from Turner – Remote sensing? || ??? || 
     551|| Harvest from Schelhaas - data through the VERIFY project || ??? || 
     552|| LCC – Luyssaert et al 2014 || XXX || 
     553|| LCC – Duveiller et al 2018 || ??? || 
    555554 
    556555==== Settings for spin-up and re-parameterization ==== 
    557556As a spin-up is costly in both time and computer resources, we need a strategy to avoid wasting these resources. Thus, the spin-up will like the parameterization, be done across scale moving from pixel to global scale. The spin-up will be done in parallel with the parameterization. Often the problems with the spin-up have a technical characters and show up for the pixels with extreme climate conditions. Before launching a longitudinal, regional or global spin-up, we should agree on the model version to use, because structural changes to the code will necessitate re-running the spin-up. The model version to use, will most likely be the version ready once the parameterization at the 1-pixel level is satisfying, and no more changes to the code need be added.  
    558557 
    559 •       Identify the variables that are targeted by the spin-up (such as NEP, heterotrophic respiration, decomposition etc). The spin-up will reveal whether parameters affecting these variables need to be tuned. The spin-up itself is also an interesting test case that could be loosely compared against data. 
    560 •       We could compare the spin-ups to maps of soil carbon stocks to check the order of magnitude. Soil carbon maps should only be formally compared with the control run (spin-up + transient simulation) for the year 2000 because that run includes the simulated effects of N-deposition, management, litter raking and land cover change. The more simple configurations of the spin-up do not account for these processes or do not account for the right sequence of processes  
     558* Identify the variables that are targeted by the spin-up (such as NEP, heterotrophic respiration, decomposition etc). The spin-up will reveal whether parameters affecting these variables need to be tuned. The spin-up itself is also an interesting test case that could be loosely compared against data. 
     559* We could compare the spin-ups to maps of soil carbon stocks to check the order of magnitude. Soil carbon maps should only be formally compared with the control run (spin-up + transient simulation) for the year 2000 because that run includes the simulated effects of N-deposition, management, litter raking and land cover change. The more simple configurations of the spin-up do not account for these processes or do not account for the right sequence of processes