-
Notifications
You must be signed in to change notification settings - Fork 1
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Create skinny topsoil layer for calculation of temperature and matric potentail at depth of microbial activity #535
Comments
One thing to add, the current setup that I'm working towards for the matric potential averaging will interpolate all of the soil layers apart from the top one. For the top one it just uses the value at the layer depth. This is obviously a bad assumption for a 50cm thick layer buts probably fine for a very thin layer. Main reason to mention this is that I think this means that there won't need to be any changes on the litter/soil model side once this skinny layer is implemented |
I don't completely get this. At the moment, we can configure any number of layers we want in |
no you are correct, this is not a difficult thing to add to the default ( it just needs to be done, hence the issue). |
Sounds good. I'm wondering if there is any mileage in defining the soil layers purely in terms of depths of particular processes:
That gives us soil layer depths in metres of [0.05, 0.1, 0.5, 1.5]. Is there any need to allow users to subdivide soil layers beyond those process driven depths? We'd generate indices to access the appropriate set of layers for each process. The order of them could in theory swap too. |
I'll have a think about this, at the moment it has mostly to do with the soil hydrology. For that I don't see why we should not define the core set by functionality. |
Is your feature request related to a problem? Please describe.
The depth of microbial activity does not necessarily align with the depth of the soil layers in the shared
LayerStructure
, see #348 . Thesoil
andlitter
model therefore need to get a weighted average between soil layers, see #532 . If the microbial activity falls above the depth of the topsoil layer, there needs to be an interpolation between the surface and the depth of the first deeper layer.The litter model currently implements a version of that which puts a constant value as the upper boundary condition; this issue asks for the roll out to the whole model. This links to updating the infiltration process #494 but could be done in two steps.
Describe the solution you'd like
Add a skinny soil layer, maybe 2-5 cm below ground, to default layer structure and update all models. For the infiltration process, the intermediate solution is to take the depth of the first two layers.
Describe alternatives you've considered
We also considered an above ground layer as an upper boundary (as this is also used in the litter model), but this became a bit tricky with the matric potential. However, it would be interesting in the future to think about litter moisture content or relative humidity in litter layer when it comes to flammability and simulations or fires.
The text was updated successfully, but these errors were encountered: