Organic matter transfer

[AlexisImbert]

Hi everyone,

I have concerns about the transfer of organic matter within the lake I want to model. I added in my inflow non-zero values of OGM_doc, OGM_don and OGM_dop but when I plot the evolution of nitrates, phosphates or chlorophyll the chronicles are identical, with or without inflow. The only parameter that changes (and this is logical) is the three components of the DOM. However in the equations of NO3 or PO4 the DOM appears so I don't understand why the transformation is not done.

Thanks for your help,

Alexis

[matthipsey]

Hi Alexis Are you setting Rdom_minerl in &aed_organic_matter? From: AlexisImbert ***@***.***> Date: Wednesday, 31 August 2022 at 4:14 pm To: AquaticEcoDynamics/libaed2 ***@***.***> Cc: Subscribed ***@***.***> Subject: [AquaticEcoDynamics/libaed2] Organic matter transfer (Issue #20) Hi everyone, I have concerns about the transfer of organic matter within the lake I want to model. I added in my inflow non-zero values of OGM_doc, OGM_don and OGM_dop but when I plot the evolution of nitrates, phosphates or chlorophyll the chronicles are identical, with or without inflow. The only parameter that changes (and this is logical) is the three components of the DOM. However in the equations of NO3 or PO4 the DOM appears so I don't understand why the transformation is not done. Thanks for your help, Alexis — Reply to this email directly, view it on GitHub<https://aus01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2FAquaticEcoDynamics%2Flibaed2%2Fissues%2F20&data=05%7C01%7Cmatt.hipsey%40uwa.edu.au%7C03df3eae59eb4f7d940b08da8b28dd90%7C05894af0cb2846d8871674cdb46e2226%7C0%7C0%7C637975304877738845%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=HVheT1EFnO3XPam0LXe401GUEH5VeWXxnI6xAUn0Bzc%3D&reserved=0>, or unsubscribe<https://aus01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fnotifications%2Funsubscribe-auth%2FAB2QIDM7JO7E6TP6SYZMWHLV34H7HANCNFSM6AAAAAAQBEFDHM&data=05%7C01%7Cmatt.hipsey%40uwa.edu.au%7C03df3eae59eb4f7d940b08da8b28dd90%7C05894af0cb2846d8871674cdb46e2226%7C0%7C0%7C637975304877895091%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=yLzk1mbF9kVedHCYH0ciaEZlFOawv8Vqgi9tozgJqlQ%3D&reserved=0>. You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>

[AlexisImbert]

Thanks for your reply,

yes i tried many values but it doesn't seems to work.

Here my latest configuration :

!-- Breakdown and mineralisation (basic pool)
Rpoc_hydrol = 0.5 ! 0.05
Rdoc_minerl = 0.5 ! 0.001
Rpon_hydrol = 0.5 ! 0.05
Rdon_minerl = 0.5 ! 0.005
Rpop_hydrol = 0.5 ! 0.05
Rdop_minerl = 0.5 !0.001
theta_hydrol = 1.08
theta_minerl = 1.08
Kpom_hydrol = 31.25
Kdom_minerl = 31.25
simDenitrification = 1

[matthipsey]

OK, some time back note there was a change to use Rdom_minerl, which has a common rate for C, N and P, in order to try to reduce parameters. I think some older examples are still in circulation where this wasn't updated properly,

libaed2/src/aed2_organic_matter.F90

Line 719 in 45e9d95

dom_mineralisation = R_dom(data%Rdom_minerl,data%Kdom_minerl,data%f_an,data%theta_minerl,oxy,temp)

If you set this in the nml, then you will get mineralisation. You can check the exact amount by plotting the diagnostic variables.

Please also note that libaed2 has been replaced by libaed-water (and others) a couple of years ago, and now users are recommended to use this new AED (instead of AED2). This code is in GLM v3.2+, and better documented here.

[AlexisImbert]

Thanks for your reply. I'm now working with the new AED and all is working well, thanks again !

Last question, do you know how could we use AED for predicting Secchi Depth ?