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


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

--

Legend:

Unmodified
Added
Removed
Modified

  • DevelopmentActivities/ORCHIDEE-DOFOCO

    v172 v173  
    384384ORCHIDEE-CN-CAN strengthen the links between sechiba and stomate. As in previous versions, stomate makes use of variables calculated in sechiba but in ORCHIDEE-CAN and ORCHIDEE-CN-CAN, sechiba requires information from stomate to work properly. For the moment set '''stomate_ok_stomate''' to y (_AUTOBLOCKER_). For the future it seems possible to prescribe LAI and assume a canopy structure but this code still needs to be restored and tested. For the time being set '''lai_map''' = n. 
    385385 
    386 == Parameterization and Evaluation (chronological order) == 
     386== Parameterization and Evaluation == 
    387387 
    388388=== Workflow === 
    … …  
    391391The work will likely require several iteration containing all or just a coupled of the following: 1) running the model, 2) tuning key parameters, 3) re-running the model 4) spin-up and 5) evaluating to model, to discern parameter values that allow us to reproduce as closely as possible the global patterns and trends for several key processes. The tests will move across scales starting at pixel scale moving to longitudinal band, European scale and ending at global scale. Although the longitudinal bands are of little use in the evaluation itself they can be considered pre-tests before global tests are run. It is hoped that longitudinal bands would speed up the tests and reduce the computational cost. 
    392392 
    393 Table 0. Workflow of the parameterization and evaluation of ORCHIDEE-CN-CAN. 
    394 Phase   Work to be done 
    395 Prepare and check       Check whether the model runs and spin-up are apparently bug free (Table 1). Check whether the scripts, tools and data are available and still working (Table 2).  
    396 Initial evaluation      Run longitudinal spin-up and check order of magnitude in soil carbon pools. Use the same runs to check response gradients to temperature, N-deposition, precipitation and management (§4.4).    
    397 Parameterization        Use FLUXNET and some additional data sources to parameterize the model at the site level (Table 3). 
    398 Final evaluation        Run a spin-up + transient simulation over Europe compare the simulation against the spatially explicit data (Table 4). 
     393Table 1. Workflow of the parameterization and evaluation of ORCHIDEE-CN-CAN. 
     394|| '''Phase''' || '''Work to be done''' ||  
     395|| Prepare and check || Check whether the model runs, the key functionality, and the spin-up are bug free (Table 2). Check whether the scripts, tools and data are available and still working (Table 3). || 
     396|| Initial evaluation || Run longitudinal spin-up and check order of magnitude in soil carbon pools. Use the same runs to check response gradients to temperature, N-deposition, precipitation and management (§5.4). ||    
     397|| Parameterization || Use FLUXNET and some additional data sources to parameterize the model at the site level (Table 4). || 
     398|| Final evaluation || Run a spin-up + transient simulation over Europe compare the simulation against the spatially explicit data (Table 5). || 
    399399 
    400400 
    401401=== Apparent bug-free === 
    402 Before starting the work proposed in Table 0 it needs to be confirmed that the model is technically capable of the tasks presented in Table 1. 
    403  
    404 Table 1. Essential technical capabilities before evaluating ORCHIDEE-CN-CAN. 
    405 Description of the task Status 
    406 Is the model stable? Can it be run for 500 - 1000 years? Can it be run for all PFTs?    Problems with grasslands 
    407 Is the N-cycle working? Can we run simulations with IMPOSE_CN=n and IMPOSE_CN=y?        OK 
    408 Is land cover change working? Can we run simulations with (changing) land cover maps (impose_veg=no)?   OK 
    409 Is forest management working? Can we read forest management maps, are the temporal trends in biomass what is expected for unmanaged, high-stand and coppice management? Still need to run idealized set-ups 
    410 Is litter raking working? Can we read litter raking maps and is the litter pool in forest decreasing and the litter pool in croplands increasing?       Still needs to be checked 
    411 Does the model run in parallel?         OK 
    412 Is the analytical spin-up working?      OK 
    413 Do we have all the driver maps for the years 1600 to 2000 for European simulations? Climate, land cover change, forest management, litter raking, and N-deposition?     OK 
    414 Do we have all the driver maps for the years 1600 to 2000 for global simulations?       Still need litter raking and FM 
    415 Read N-deposition maps through COMP     ??? 
     402Before starting the work proposed in Table 1 it needs to be confirmed that the model is technically capable of the tasks presented in Table 2. 
     403 
     404Table 2. Essential technical capabilities before evaluating ORCHIDEE-CN-CAN. 
     405|| '''Description of the task'''' || '''Status''''|| 
     406|| Is the model stable? Can it be run for 500 - 1000 years? Can it be run for all PFTs? || Recent problems with grasslands || 
     407|| Is the N-cycle working? Can we run simulations with IMPOSE_CN=n and IMPOSE_CN=y?     || OK || 
     408|| Is land cover change working? Can we run simulations with (changing) land cover maps (impose_veg=no)? || OK || 
     409|| Is forest management working? Can we read forest management maps, are the temporal trends in biomass what is expected for unmanaged, high-stand and coppice management? || Still need to run idealized set-ups || 
     410|| Is litter raking working? Can we read litter raking maps and is the litter pool in forest decreasing and the litter pool in croplands increasing? || Still needs to be checked || 
     411|| Does the model run in parallel? || OK || 
     412|| Is the analytical spin-up working? || OK || 
     413|| Do we have all the driver maps for the years 1600 to 2000 for European simulations? Climate, land cover change, forest management, litter raking, and N-deposition? ||       OK || 
     414|| Do we have all the driver maps for the years 1600 to 2000 for global simulations? || Still need litter raking and FM || 
     415|| Read N-deposition maps through COMP  || In progress || 
    416416 
    417417 
    418418=== Data availability === 
    419419Observational data products for model-data comparison can be found at:  /home/satellites5/maignan/ORCHIDEE/EVALUATION. 
    420 •       phenology  
     420* Phenology  
    421421◦       GIMMS (LAI and FPAR3g) 
    422 •       Forest structure 
     422* Forest structure 
    423423◦       Remote sensing products of biomass (temperate and boreal maps, i.e., Turner) 
    424424◦       Biomass of EU forest from JRC (Europe) 
    … …  
    427427◦       Forest basal area (Europe) 
    428428◦       Canopy height (Global) 
    429 •       NPP 
     429* NPP 
    430430◦       Site-level NPP database Luyssaert et al 2007 
    431 •       NEP 
     431* NEP 
    432432◦       FLUXNET site-level data 
    433 •       TER 
     433* TER 
    434434◦       FLUXNET site-level data 
    435 •       GPP 
     435* GPP 
    436436◦       FLUXNET site-level data 
    437437◦       EC-based upscaled GPP, i.e., Jung 
    438 •       NPP/GPP 
     438* NPP/GPP 
    439439◦       site-level data and regional patterns, i.e., Campioli et al 2015 
    440 •       Soil hydrology 
     440* Soil hydrology 
    441441◦       ESA CCI ECV 
    442442◦       measurements from Brazil (ABRACOS product) 
    443443◦       River discharge records from selected gauges (Global coverage) 
    444 •       Albedo 
     444* Albedo 
    445445◦       MODIS or GlobAlbedo for albedo 
    446 •       Energy (sensible and latent heat) 
     446* Energy (sensible and latent heat) 
    447447◦       GLEAM for evapotranspiration 
    448448◦       EC-based upscaled evapotranspiration, i.e., Jung 
    449 •       Tree ring width (not on /home/satellites5/) 
     449* Tree ring width (not on /home/satellites5/) 
    450450◦       ITRDB 
    451 •       LCC 
     451* LCC 
    452452◦       Luyssaert et al 2014 – FLUXNET changes in Rn, LE, H, G, albedo 
    453453◦       Duveiller et al 2018 – Remote sensing changes in Rn, LE, H+G, albedo  
    454 •       NFI 
     454* NFI 
    455455◦       France, Spain, Germany and Sweden (/home/satellites5/) 
    456456◦       EU-wide data through the VERIFY project? 
    … …  
    458458There are several data tools (ATLAS, Mapper and Jerome’s) to help compare model outputs with observation, which we might be able to use. If we will not be using the available tools for comparison, we need to preprocess the observational data products to produce global means, time series, decadal averages of spatial patterns etc. The analyses presented in Naudts et al 2015 are a good starting point for the evaluation of ORCHIDEE-CN-CAN. The following scripts are available and have been tested and re-activated 
    459459 
    460 Table 2. Model evaluations scripts available per November 2018 
    461 Description     Status  Path 
    462 Extract species-level productivity from French NFI data and compare against simulated productivity – looking to replace the French data by EU-wide data through the VERIFY project      OK      dofoco/dofoco/SCRIPTS/XXX 
    463 Extract Jung GPP and compare against spatially explicit simulations for 8 regions in Europe     ???     dofoco/dofoco/ SCRIPTS/XXX 
    464 Extract Jung evapotranspiration and compare against spatially explicit simulations for 8 regions in Europe      ???     dofoco/dofoco/ SCRIPTS/XXX 
    465 Extract MODIS albedo (NIR & VIS) and compare against spatially explicit simulations for 8 regions in Europe     ???     dofoco/dofoco/ SCRIPTS/XXX 
    466 Extract height and compare against spatially explicit simulations for 8 regions in Europe       ???     dofoco/dofoco/ SCRIPTS/XXX 
    467 Extract effective LAI and compare against spatially explicit simulations for 8 regions in Europe        ???     dofoco/dofoco/ SCRIPTS/XXX 
    468 Extract BA and compare against spatially explicit simulations for 8 regions in Europe   ???     dofoco/dofoco/ SCRIPTS/XXX 
    469 Extract half-hourly EC observations and compare with half-hourly site level simulations for NEP, GPP, TER, and evapotranspiration       ???     ENSEMBLE/FLUXNET 
    470 Extract half-hourly EC observations and compare with half-hourly site level simulations for albedo      ???     ENSEMBLE/FLUXNET 
    471  
     460Table 3. Model evaluations scripts available per November 2018 
     461|| '''Description''' || '''Status''' || '''Path''' || 
     462|| Extract species-level productivity from French NFI data and compare against simulated productivity – looking to replace the French data by EU-wide data through the VERIFY project || OK || dofoco/dofoco/SCRIPTS/XXX || 
     463|| Extract Jung GPP and compare against spatially explicit simulations for 8 regions in Europe  || OK || dofoco/dofoco/ SCRIPTS/XXX || 
     464|| 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 || 
     466|| Extract effective LAI and compare against spatially explicit simulations for 8 regions in Europe || OK || dofoco/dofoco/ SCRIPTS/XXX || 
     467|| Extract BA and compare against spatially explicit simulations for 8 regions in Europe || OK || dofoco/dofoco/ SCRIPTS/XXX || 
     468|| Extract half-hourly EC observations and compare with half-hourly site level simulations for NEP, GPP, TER, and evapotranspiration || ??? || ENSEMBLE/FLUXNET || 
    472469 
    473470=== Parameterization === 
    … …  
    485482The following parameterization approach – making use of parameters that were already shown to be sensitive to tuning – is proposed: 
    486483 
    487 Table 3. 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. 
     484Table 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. 
    488485Process Parameter(s)    Tool 
    489486Onset of growing season and start of senescence thresholds for phenology        FLUXNET