Soil Carbon is too low.

[rosiealice]

https://ns9560k.web.sigma2.no/datalake/diagnostics/noresm/kjetisaa/NorESM_Key_Simulations/n1850.ne16pg3_tn14.noresm3_0_beta09.Run4.2026-01-06/OBS_comparison/ANN/n1850.ne16pg3_tn14.noresm3_0_beta09.Run4.2026-01-06_compare_TOTSOMC_1m_HWSD_ANN_0221-0230.png

[rosiealice]

It's not -that- bad!

https://www.ilamb.org/CMIP6/historical/EcosystemandCarbonCycle/SoilCarbon/HWSD/HWSD.html

[rosiealice]

Just before christmas I did a simulation aimed at trying to improve the soil carbon:

diff /cluster/shared/noresm/inputdata/lnd/clm2/paramdata/fates_params_sci.1.88.6_api.42.0.0_14pft_nor_sci1_api1_c251204.cdl fates_params_sci.1.88.6_api.42.0.0_14pft_nor_sci3_api1_c251218.cdl
<     0.03, 0.03, 0.03, 0.04, 0.04 ;
---
>     0.022, 0.022, 0.022, 0.04, 0.04 ;
1111c1111
<  fates_allom_l2fr = 1.5, 1.5, 1, 1, 1, 1, 1, 0.7, 0.7, 1, 1.7, 2, 1.5, 2.3 ;
---
>  fates_allom_l2fr = 1.5, 0.75, 1, 1, 1, 1, 1, 0.7, 0.7, 0.75, 1.4, 1.5, 1.5, 2.3 ;
1482c1482
<  fates_mort_ip_size_senescence = 200, 190, 175, 175, 175, 175, 175, _, _, _, 
---
>  fates_mort_ip_size_senescence = 190, 190, 175, 175, 175, 175, 175, _, _, _,

The results and logic are described in here:
https://docs.google.com/presentation/d/1cyT6C82dEoXenx4AvJsvHto76KcoVcQWv0xjEyDlNkI/edit?slide=id.g3b193e1c7b5_0_15#slide=id.g3b193e1c7b5_0_15

[rosiealice]

To Do:

Set off AD spinup with these modified parameters.
Compare with Jessie's beta09 run.
@JessicaNeedham can you send me the script/path to the script?

[rosiealice]

@kjetilaas the parameters to increase the soil carbon are here:

/cluster/shared/noresm/inputdata/lnd/clm2/paramdata/fates_params_sci.1.88.6_api.42.0.0_14pft_nor_sci1_api1_soilCinc.nc

The differences are

diff fates_params_sci.1.88.6_api.42.0.0_14pft_nor_sci1_api1_soilCinc.cdl fates_params_sci.1.88.6_api.42.0.0_14pft_nor_sci1_api1_c251204.cdl 
1c1
< netcdf fates_params_sci.1.88.6_api.42.0.0_14pft_nor_sci3_api1_c251218 {
---
> netcdf fates_params_sci.1.88.6_api.42.0.0_14pft_nor_sci1_api1_c251204 {
1061c1061
fates_allom_d2bl1
<     0.022, 0.022, 0.022, 0.04, 0.04 ;
---
>     0.03, 0.03, 0.03, 0.04, 0.04 ;
1111,1112c1111
<  fates_allom_l2fr = 1.5, 0.75, 1, 1, 1, 1, 1, 0.7, 0.7, 0.75, 1.4, 1.5, 1.5, 
<     2.3 ;
---
>  fates_allom_l2fr = 1.5, 1.5, 1, 1, 1, 1, 1, 0.7, 0.7, 1, 1.7, 2, 1.5, 2.3 ;