Biogenic heating and soil decomposition

[caroduro]

Are there any linked Issues or Pull Requests?

UM issue #69, JULES issue #59

Brief description

Catherine Luke and Joe Clarke used biogenic heating of soil carbon decomposition for their PhDs. We would like to add it here. There is initial code for it here (its fairly simple):
(Initial branch: https://github.com/eleanorgb/jules_wiemip/tree/biogenic_heating)

JULES Developing branch: https://github.com/caroduro/jules_for_develops/tree/59_biogenic_heating
JULES Test branch: https://github.com/caroduro/jules_for_develops/tree/test_59_biogenic_heating

UM Developing branch: https://github.com/caroduro/um_for_develops/tree/59_biogenic_heating
UM Test branch: https://github.com/caroduro/um_for_develops/tree/test_59_biogenic_heating

lfric_apps Developing branch: https://github.com/caroduro/lfric_apps_for_develops/tree/59_biogenic_heating
lfric_apps test branch https://github.com/caroduro/lfric_apps_for_develops/tree/test_59_biogenic_heating

However, I would like to change a couple of things.

1. Please could the switch (l_bgc_heat) go in the soil_biogeochemistry namelist not the soils namelist.
2. the heat_of_respiration = 3.9e7 is set in the code, please can we also set it in the soil_biogeochemistry namelist and trigger it if l_bgc_heat is true.
3. for the metadata - we need a warning that if l_q10 is false we wont have a peak in the function that shows the dependence of respiration on temperature - this is required in case we have a compost bomb and helps stabilise the model.
4. I dont know the reference to Catherine Lukes thesis but she was at Exeter University.
5. For the documentation - we need to add some text:
   l_bgc_heat=True:: Include the the release of heat associated with soil carbon decomposition on soil temperature. if l_q10 is True we also include the a peak in the dependence of respiration on soil temperature (Equation 4.4 and Figore 4.1) in Catherine Lukes thesis. This will eliminate the problem of unrealistically high soil temperature occuring in the model and represents the reduction in microbial respiration when the temperatures are very high. Parameters required heat_of_respiration and default value as above and units (J/kgC).

Im not entirely sure it makes sense without l_layeredc=true (maybe someone else can confirm)

I cant access this paper, but think its probably the one to reference: Luke and Cox - https://bsssjournals.onlinelibrary.wiley.com/doi/10.1111/j.1365-2389.2010.01312.x

I have done some runs and will add some results here shortly.

Further details of the issue.

No response