wiki:DevelopmentActivities/ORCHIDEE-CNP

Version 128 (modified by dgoll, 7 years ago) (diff)

--

ORCHIDEE-CN-P (former branch ORCHIDEE-CNP)

This page describes the phosphorus cycle in ORCHIDEE-CN-P. It is based on ORCHIDEE-CN, which was extended and corrected as described here: https://forge.ipsl.jussieu.fr/orchidee/wiki/Branches/MergeOCN/Goll as well as additional non-documented bugfixes to avoid negative pools due to machine precisions as well as bugs.

The phosphorus cycle is an adaptation of the the model described by http://www.biogeosciences.net/9/3547/2012/bg-9-3547-2012.html . Nonetheless, the complexity was substantially increased due to the more detailed representation of the C and N cycle in ORCHIDEE compared to JSBACH.

 Schematic representation of the key processes represented in ORCHIDEE-CNP Schematic representation of the key processes represented in ORCHIDEE-CNP

3. Conceptual modifications to the nitrogen cycle

3.1 soil mineral N concentration in soil solution

Following Smith et al (2014), I introduced the use of the maximum water holding capacity of soils (max_var_eau) to approximate pore space which to derive the average soil mineral N concentration in solution. The use of the actual water volume can not be recommended as we this would lead to high N concentration in soil water when soil water is very low. As we do not account for the inhibition of replenishment of mineral N in the soil solution around roots when soil water is scarce.

3.2 soil mineral N dynamics

ORCHIDEE-CN performs poorly in respect to soil mineral N dynamics. This is due to two reason. First, ORCHIDEE-CN is based on Zaehle & Friend (2010) and does not incorporate the modification to the DNDC formulations which avoid instabilities when run on global scale (Zaehle et al., 2011). Second, most of the modifications in ORCHIDEE-CN I assume are not tested, as most of the environmental scaling function "explode" with the parameter set used, or are not sound imo.

I step by step replaced all modification by using the formulation in Zaehle et al. (2011) (http://://www.nature.com/ngeo/journal/v4/n9/extref/ngeo1207-s1.pdf), testing for the size of the mineral N stocks as well as ensuring all scaling function work well in an global environmental space.

Final mineral N dynamics The mineral N dynamics in ORCHIDEE-CN-P are now exactly as described in Zaehle et al. (2011), except that I did not take over the rain drainage formulation as we use an improved hydrological scheme compared to Zaehle et al. (2011) which is able to simulate the rain drainage as a result of more soil layers (11 layers).

The mineral N stocks of the test sites are all in realistic ranges. Pattern correlation test with N15 derived fraction of denitrification show poor correlation due to high NOx emissions. This needs to be corrected in the near future.

4. Re-calibration of the carbon cycle

The calibration of ORCHIDEE-CAN of the carbon cycle in combination with the Vcmax being derived from leaf nitrogen content, leads to unrealistic C cycle (CUE, biomass stocks, biomass turnover, etc.). Using as a guide rough ranges of observation-based estimates of aboveground biomass, LAI, CUE, and biomass turnover, we manipulated parameters of the pipe model (TAU_SAP, KLATOSA_MIN, K_LATOSA_MAX) as well as the inversion of the rate constant for tree turnover RESIDENCE_TIME and the factor relating maintenance respiration to tissue N, COEFF_MAINT.

The main issue with the previous calibration (Naudts et al GMD) is, that FLUXNET GPP & energy fluxes as well as LAI estimates are poor constraints on biomass and biomass turnover when calibrating the ORCHIDEE-CAN parameter. The clear cuts in the management simulations in the GMD hide the potential biomass stocks. Even the biomass stocks estimates which exclude managed forest are smalled than the ones in ORCHIDEE.

The range of K_LATOSA prescribed makes things tricky as total biomass is very sensitive to the actual values of K_LATOSA. The dynamics of K_LATOSA were several times revised and it is not clear what is going on.

I derived a set of parameters which result in kind of okay LAI, total biomass, biomass turnover and CUE. Still some PFTs behave not really realistic, like DBF PFTs.

ISSUE with my calibration: To facilitate the calibration I derive maintenance respiration only from leaf N and not from sapwood N too. This allows to manipulate total biomass via K_LATOSA without affecting the CUE. It would be better to have it coupled as this adds a constraint to K_LATOSA. Best would be to include even more constraints on K_LATOSA for example hydrological constraints (like the stuff Emilie develops).

5. New input files

You can find information on the new input files needed for the P cycle here: https://forge.ipsl.jussieu.fr/orchidee/wiki/DevelopmentActivities/ORCHIDEE-CNP/newInput

6. Notes regarding Vcmax and leaf nutrients

https://forge.ipsl.jussieu.fr/orchidee/wiki/DevelopmentActivities/ORCHIDEE-CNP/Vcmax

7. unresolved issues (ideas for future developments)

There are still yet unresolved issue with the model. They are listed here: https://forge.ipsl.jussieu.fr/orchidee/wiki/DevelopmentActivities/ORCHIDEE-CNP/issues

8. Howto install, compile & run the model

You can find information on how to setup the simulations here: https://forge.ipsl.jussieu.fr/orchidee/wiki/DevelopmentActivities/ORCHIDEE-CNP/howtoUse

9. Prescribing tissue stoichiometry: WARNING

At the moment prescribing tissue CN and NP ratio have only a purpose when spinning up the model, but cannot be used for the analysis of nutrient limitation. The effect of nutrient limitation on growth is not handled consistently when prescribing leaf CN ratios which are not optimal ratios (CNmin and NPmin). However, the effect of nutrient limitation on growth is substantial when CN and/or NP deviate from the optimal ratios, often they are larger than the effect on GPP (as maintenance respiration is coupled to the leaf nitrogen concentration, a lower N concentration means lower Vcmax and lower respiration). As growth limitation is not handled correctly, it is not possible to derive any conclusion from simulations prescribing non-optimal ratios regarding nutrient limitation on vegetation. (apparently you just modified the value of Vcmax via leaf nitrogen concentration).

At the moment, I have no idea how to fix this in a feasible way.

Attachments (11)

Download all attachments as: .zip