Queries on simulating SOM decomposition in aerobic and anaerobic conditions

19 views
Skip to first unread message

Bisesh Joshi

unread,
Jun 12, 2025, 11:47:40 AM6/12/25
to pflotra...@googlegroups.com

Hello all,
I'm trying to simulate the anaerobic and aerobic ammonium mineralization in Pflotran. Anaerobic mineralization should be from  decomposition of soil organic matter (SOM) coupled with Fe3+ reduction to Fe2+ in absence of O2, while aerobic ammonium mineralization , SOM decomposes in presence of O2. 
Please let me know if anybody has any suggestions to directly simulate this in pflotran.
Thanks in advance!

Best regards,
Bisesh 
--
Bisesh Joshi
PhD student
Water Science and Policy
Graduate College
University of Delaware

Hammond, Glenn E

unread,
Jun 13, 2025, 11:26:28 AM6/13/25
to pflotra...@googlegroups.com

Bisesh,

 

These reactions are possible. Would you be willing to share your mass action equations (A + B -> C + D) and the desired rate expressions. If proprietary, you are welcome to send them directly to me.

 

Glenn

 

From: pflotra...@googlegroups.com <pflotra...@googlegroups.com> on behalf of Bisesh Joshi <bjo...@udel.edu>
Date: Thursday, June 12, 2025 at 9:47
AM
To: pflotra...@googlegroups.com <pflotra...@googlegroups.com>
Subject: [pflotran-users: 8430] Queries on simulating SOM decomposition in aerobic and anaerobic conditions

Check twice before you click! This email originated from outside PNNL.

 

--
You received this message because you are subscribed to the Google Groups "pflotran-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to pflotran-user...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/pflotran-users/CAEvi1yQx1qs-9K6ho3X%2BUBb5VweeaDSZNdf4%2BP9VbKRqgP0VPA%40mail.gmail.com.

Bisesh Joshi

unread,
Jun 13, 2025, 1:27:41 PM6/13/25
to pflotra...@googlegroups.com
Hi Glenn,
Thanks for the response.
This is what I'm thinking: I want to use  SOM (C6H13NO5) as a mineral, producing the following products in aerobic and anaerobic conditions.
SOM + 6O2 + H2O→NH4+ + 6HCO3 + 5H+ (aerobic conditions)

SOM + 24Fe3+ + 13H2O→NH4+ + 6HCO3 + 24Fe2+ + 29H+ (anaerobic conditions, O2 inhibiting this reaction)

Additionally,  Fe3+ here also comes from the mineral Goethite. 

On a different note, I have got this mineral reaction in one database: 'SOM1' 100. 3 12.0 'DOM1' -12.0 O2(aq) 1.0 NH4+ 100. 100. 100. 100. 100. 100. 100. 100. 180.1566 'TAO-glucose’. Doesn't this reaction imply that SOM1 is stable and not producing DOM1 here?

Thanks again for your help.


Best regards

Bisesh


Hammond, Glenn E

unread,
Jun 17, 2025, 7:13:15 PM6/17/25
to pflotra...@googlegroups.com

From: pflotra...@googlegroups.com <pflotra...@googlegroups.com> on behalf of Bisesh Joshi <bjo...@udel.edu>
Date: Friday, June 13, 2025 at 10:27
AM
To: pflotra...@googlegroups.com <pflotra...@googlegroups.com>
Subject: Re: [pflotran-users: 8434] Queries on simulating SOM decomposition in aerobic and anaerobic conditions

Hi Glenn,

Thanks for the response.

This is what I'm thinking: I want to use  SOM (C6H13NO5) as a mineral, producing the following products in aerobic and anaerobic conditions.

SOM + 6O2 + H2O→NH4+ + 6HCO3− + 5H+ (aerobic conditions)

SOM + 24Fe3+ + 13H2O→NH4+ + 6HCO3− + 24Fe2+ + 29H+ (anaerobic conditions, O2 inhibiting this reaction)

Additionally,  Fe3+ here also comes from the mineral Goethite. 

 

These reactions can be modeled once SOM -> DOM using the MICROBIAL_REACTIONs with multiplicative Monod and inhibition factors. Otherwise, you will need a reaction sandbox to implement the rate expressions using SOM as the donor.

 

On a different note, I have got this mineral reaction in one database: 'SOM1' 100. 3 12.0 'DOM1' -12.0 O2(aq) 1.0 NH4+ 100. 100. 100. 100. 100. 100. 100. 100. 180.1566 'TAO-glucose’. Doesn't this reaction imply that SOM1 is stable and not producing DOM1 here?

Thanks again for your help.

 

The mass action equation for this reaction is

 

SOM + 12 O2(aq) -> 12 DOM + NH4+

 

This mineral reaction produces 12 moles of DOM for each mole of SOM consumed. If you use this mineral reaction, you need to specify DISSOLUTION_RATE_CONSTANT and PRECIPITATION_RATE_CONSTANT and set the precipitation constant to zero to avoid re-precipitating the SOM.

 

Glenn

 

Best regards

Bisesh

Bisesh Joshi

unread,
Jun 18, 2025, 9:53:17 AM6/18/25
to pflotra...@googlegroups.com
Thanks Glenn. 
Can you please share any similar reaction sandbox file if you have one?

Best regards,
Bisesh


Reply all
Reply to author
Forward
0 new messages