ORCHIDEE-TRUNK OPTIMIZATION
ROUND8 (11/03/2024)

Select: config PFT10 PFT11 PFT12 PFT13 PFT14 PFT15

Update (29/03/2024)

1. Simulations are reproduced with ORC3 (rev8120) at the same sites
2. Spinup time-series are added for ORC3 (variable: total soil carbon)
3. Prior daily time-series and seasonal cycles across 5Y are added for ORC3 and ORC4 (variables: fluxlat, fluxsens, gpp, lai, nee)

Configuration (11/03/2024)

1. Model version: ORCHIDEE_trunk revision 8464.
2. Sites: selection of 5 sites from the ensemble used in SIF project (selected by C.Bacour, 254 sites across all PFTs).
3. Forcing: global CRUJRA v2.4 (closest pixel to each site location).
4. PFT composition: 100%-composition of a corresponding PFT.
5. Spinup: 240Y recycling 1901-1910 forcing with fixed ATM_CO2 = 350 ppm.
6. Transient: 100Y across 1901-2000 with fixed ATM_CO2 = 350 ppm.
7. Optimization: 20Y (2001-2020) with yearly outputs, the last 10Y are used in optimization.
8. Control variables with ranges of variation:

   Variable	PFT10	PFT11	PFT12	PFT13	PFT14	PFT15
   GPP [gC/m2/y]	500 - 1000	600 - 1400	1000 - 1400	600 - 1600	600 - 1400	500 - 1200
   LAI_MEAN_GS	0.5 - 1.5	0.5 - 1.5	0.5 - 1.5	0.5 - 1.5	0.5 - 2	0.5 - 1.5
   ng_ratio	0.4 - 0.6	0.32 - 0.58	0.53 - 0.7	0.51 - 0.62	0.32 - 0.58	0.32 - 0.5
   lr_ratio	0.5 - 1.5
   ge_ratio	1 - 3	3 - 8.5	1.2 - 2.8	3 - 8.5	1 - 2	1 - 3.2
   bm_ratio	0 - 0.1

   ng_ratio = NPP / GPP
   lr_ratio = LEAF_M_c / ROOT_M_c
   bm_ratio = (LABILE_M_c + RESERVE_M_c) / (SAP_M_BE_c + SAP_M_AB_c + HEART_M_BE_c + HEART_M_AB_c)
   ge_ratio = GPP / EVAP


9. Optimized parameters with default values and ranges of variation
   (vmin = 0.5*vdef and vmax = 2*vdef for all parameters except for VMAX_UPTAKE_XXX vmin=0.25*vdef):

   Parameter	PFT10	PFT11	PFT12	PFT13	PFT14	PFT15
   K_LATOSA	7000 [3500 - 14000]	5500 [2750 - 11000]	9500 [4750 - 19000]	11000 [5500 - 22000]	8500 [4250 - 17000]	10000 [5000 - 20000]
   K_BELOWGROUND	4.0E-7 [2.0E-7 - 8.0E-7]	4.0E-7 [2.0E-7 - 8.0E-7]	4.0E-7 [2.0E-7 - 8.0E-7]	4.0E-7 [2.0E-7 - 8.0E-7]	4.0E-7 [2.0E-7 - 8.0E-7]	4.0E-7 [2.0E-7 - 8.0E-7]
   VMAX_UPTAKE_NH4	9 [2.25 - 18]	9 [2.25 - 18]	9 [2.25 - 18]	9 [2.25 - 18]	9 [2.25 - 18]	9 [2.25 - 18]
   VMAX_UPTAKE_NO3	9 [2.25 - 18]	9 [2.25 - 18]	9 [2.25 - 18]	9 [2.25 - 18]	9 [2.25 - 18]	9 [2.25 - 18]
   COEFF_MAINT_INIT	0.06 [0.03 - 0.12]	0.05 [0.025 - 0.1]	0.05 [0.025 - 0.1]	0.08 [0.04 - 0.16]	0.07 [0.035 - 0.14]	0.13 [0.065 - 0.26]
   SLAINIT	0.031 [0.0155 - 0.062]	0.031 [0.0155 - 0.062]	0.02 [0.01 - 0.04]	0.02 [0.01 - 0.04]	0.031 [0.0155 - 0.062]	0.031 [0.0155 - 0.062]
   NUE_OPT	45 [22.5 - 90]	45 [22.5 - 90]	60 [30 - 120]	60 [30 - 120]	45 [22.5 - 90]	45 [22.5 - 90]
   LONGEVITY_SAP	180 [90 - 360]	180 [90 - 360]	480 [240 - 960]	4801 [2400.5 - 9602]	180 [90 - 360]	240 [120 - 480]

   K_LATOSA_MAX = K_LATOSA_MIN = K_LATOSA
   


10. Observation error (R) is set as RMSD for simulations time-series vs the mean between min/max values of range of variation.
11. Parameter error (B) is set as 15% of the parameter range of variation.
12. Cost-function is calculated with the standard formula, but the obs-sim difference is calculated as deviation from the prescribed range of variation and set to 0 for the time-series part falling inside the range.
13. Minimization method: Genetic Algorithm (population = 30, number of iterations = 15).
14. Computing platform: LSCE obelix with 15 processors.