the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 30 Apr 2024
Lambda-PFLOTRAN 1.0: Workflow for Incorporating Organic Matter Chemistry Informed by Ultra High Resolution Mass Spectrometry into Biogeochemical Modeling
Abstract. Organic matter (OM) composition plays a central role in microbial respiration of dissolved organic matter and subsequent biogeochemical reactions. Here, a direct connection of organic carbon chemistry and thermodynamics to reactive transport simulators has been achieved through the newly developed Lambda-PFLOTRAN workflow tool that succinctly incorporates carbon chemistry data generated from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) into reaction networks to simulate organic matter degradation and the resulting biogeochemistry. Lambda-PFLOTRAN is a python-based workflow, executed through a Jupyter Notebook interface, that digests raw FTICR-MS data, develops a representative reaction network based on substrate-explicit thermodynamic modeling (also termed lambda modeling due to its key thermodynamic parameter λ used therein), and completes a biogeochemical simulation with the open source, reactive flow and transport code PFLOTRAN. The workflow consists of the following five steps: configuration, thermodynamic (lambda) analysis, sensitivity analysis, parameter estimation, and simulation output and visualization. Two test cases are provided to demonstrate the functionality of the Lambda-PFLOTRAN workflow. The first test case uses laboratory incubation data of temporal oxygen depletion to fit lambda parameters (i.e., maximum utilization rate and microbial carrying capacity). A slightly more complex second test case fits multiple lambda formulation and soil organic matter release parameters to temporal greenhouse gas generation measured during a soil incubation. Overall, the Lambda-PFLOTRAN workflow facilitates upscaling by using molecular-scale characterization to inform biogeochemical processes occurring at larger scales.
- Preprint
(977 KB) - Metadata XML
-
Supplement
(573 KB) - BibTeX
- EndNote
Status: closed
-
RC1: 'Comment on gmd-2024-34', Anonymous Referee #1, 26 Jun 2024
General comments:
The presented work is a good contribution to improving the reactive transport modeling approach. It develops a representative reaction network based on substrate-explicit thermodynamic modeling and completes a biogeochemical simulation with a reactive transport model.
However, I think the current version cannot be considered for publication in GMD. There are some possible reasons:
(a) Poor model validation. A strict validation checks whether a new model approach reflects the actual or desired physical and biochemical behavior, considering using a reactive transport model. In this study, I noticed that there are very limited observations. For example, the work only involved limited observed O2 (Fig.4a), CO2 (Fig.5a), and total organic carbon datasets (Fig.5b). I also noticed that the model cannot capture observed total organic carbon (Fig.5b).
(b) Incomplete simulation experiments. The work just considered two test cases (oxygen depletion and respiration). Too small a test size may lack sufficient statistical power to detect new model capabilities.
(c) Poor analysis. The work shows limited biogeochemical analysis. An advantage of RTM is its ability to elucidate complex biogeochemical dynamics. For example, How do oxygen dynamics and total carbon influence the variations of pH and C, N, and P species over time (Fig.5a)? Although I understand this work is technical, the justification for these time-series variations should be carefully discussed.
Citation: https://doi.org/10.5194/gmd-2024-34-RC1 -
AC2: 'Reply on RC1', Katherine Muller, 21 Aug 2024
We understand the reviewer’s comment on the importance of model validation. We agree that more complete data for the test cases and additional datasets would be greatly beneficial, yet high-resolution organic matter characterization (e.g., FTICR-MS) is not commonly measured, which makes it not feasible to thoroughly validate the model given the lack of required data with both incubation and organic matter speciation. We hope that building this modeling framework will encourage others to measure organic matter speciation via FTICR-MS in the future. As more data becomes available, the Lambda-PFLOTRAN workflow can be further validated and refined.
We have clarified the scope of this paper on line 77, emphasizing that our work aims to detail a new capability and associated workflow for incorporation of organic matter characterization into reactive transport models. Model validation is important and new data collection could be designed to make it possible in future studies.
The illustrative model results presented here consider only aerobic respiration through the FTICR-MS informed reaction network. If desired, the PFLOTRAN input deck can be expanded and customized to include additional processes and full geochemistry.
We have added the following statement on lines 325-330 to capture this point:
“Three additional processes were added to describe the experimental conditions of Test Case 2 more accurately (i.e., release of carbon, nitrogen and sustained aerobic conditions); however, a PFLOTRAN input deck can be expanded and customized to include a host of additional processes and full geochemistry for a specific system of interest. For instance, aqueous complexation, mineral dissolution and precipitation, sorption, and redox reactions can be added, all of which can influence the resultant pH and carbon, nitrogen, and other nutrient dynamics.”
Citation: https://doi.org/10.5194/gmd-2024-34-AC2
-
AC2: 'Reply on RC1', Katherine Muller, 21 Aug 2024
-
RC2: 'Comment on gmd-2024-34', Anonymous Referee #2, 28 Jun 2024
This manuscript describes a new workflow for integrating detailed organic matter data from FTICR-MS into a reactive transport modeling framework for simulating biogeochemical interactions. The workflow includes sensitivity analysis and parameter estimation capabilities. I think this has the potential to be a highly valuable tool for connecting biogeochemical models with organic matter measurements at a level of chemical detail that has historically been a major challenge for biogeochemical model-data integration. The workflow covers many of the pain points for this type of model-data connection, including aggregating large numbers of compounds into a smaller, tractable set of representative compounds, assigning meaningful model parameters to those compounds, and converting the resulting reaction network into a reactive transport model simulation configuration automatically.
The manuscript describes the underlying assumptions and theoretical framework in clear terms and acknowledges important caveats such as the lack of quantitative information about the relative amounts of different compounds. The steps for using the workflow are also clearly described, which fits with the goal of making the workflow accessible and reproducible for other scientists. The two provided test cases are useful concrete examples of the workflow, important parameters, and useful data types for parameter estimation.
I did think that the results shown from the test cases did not fully make the case for how the FTICR-MS data add value to the resulting simulations. While the comparison of the full reaction network to simulations using multiple compounds does show a difference in overall predicted concentrations relative to a simplified reaction network with only one organic matter compound, it is difficult to see how the high resolution of chemical compounds, which is the hallmark of the FTICR-MS analysis, specifically contributes to the simulation results. It would be helpful to include a visualization of the distribution of compounds, how they translate to reaction networks with different properties, and how that ultimately affects the model results. For example, is the reaction network derived from the Test Case 1 FTICR-MS measurements meaningfully different from that of Test Case 2 in terms of its stoichiometry or distribution of reaction rates? Does the specific FTICR-MS data from each experiment yield better model results than FTICR-MS data from a different experiment? What if the null hypothesis was 10 compounds derived from a different set of measurements, rather than 1 generic compound? This comparison could more directly show the value of those measurements from a particular experiment.
Other comments:
Line 97: is yOM here the same as yOC in Equation 2?
Line 101: I think “1000” is missing from this sentence
Line 123: HS should be HS-
Line 192: Given the importance of specifying parameters to be estimated for the workflow, I think it would be helpful to list all the relevant parameters in one table, with explanations of what they mean and perhaps with suggestions for what kind of measurements could constrain them or what would be a good context for including them in the parameter estimation and what reasonable ranges of values might be.
Section 2.3.4: When discussing the parameter optimization, it would be helpful to start with some explanation of what kind of observations are useful for comparing with the model. While examples are provided in the test cases, I think some general explanation earlier on would be helpful as well.
Line 256: I’m curious how accurate the assumption of equal division of C mass across bins is. Do any data exist that could inform this? The approach might be quite sensitive to this assumption.
Figure 3: I was confused why k_deg and C_inhibit are shown in this figure, even though they aren’t part of the parameters being tested. And I did not understand the significance of the time dimension. Does this indicate whether model ensemble members are diverging or converging over the course of the simulation? It could use some more explanation of how to interpret that part.
Figure 5: The model generally fails to reproduce the observed total carbon concentrations. I might guess that the assumption of how much carbon gets solubilized from the solid state over time, or the characteristics of that soluble carbon, are not correct for this system. I would not necessarily expect a model to get every aspect of the system right, especially if it relates to a steady state process involving solid organic matter, which this framework is not designed to replicate. But I think this should be at least acknowledged and explained in the text and currently it feels like this part of the result is ignored.
Citation: https://doi.org/10.5194/gmd-2024-34-RC2 -
AC3: 'Reply on RC2', Katherine Muller, 21 Aug 2024
The following are responses and updates based on reviewer comments:
- We have revised Figure 3, the associated figure description, and companion text to help highlight the impact of OM speciation for Test Case 1 (lines 291-298).
- We have added a new section titled “Variability and Impact of Organic Matter Speciation” to address these comments. Section 4 explores the differences in OM chemistry between samples (Figure 5), the effects of OM chemistry on the lambda-derived reaction networks (Figure 6), and how lambda-derived reaction networks compare to an assumed generic CH2O bulk carbon species (Figure 7). Additionally, this section also examines the influence of OM chemistry through forward simulations, where only FTICR-MS input data is varied (Figure 8).
- We have updated the minor edits including: yOM, adding 1000, and fixing the subscript on HS-.
- Model parameters are defined in the Methods section (Section 2). While we understand the interest in gaining further insight into optimal model parameterization, we believe this is beyond the scope of the current manuscript. However, we agree that additional work focused on parameterization, experimental design, and data collection would be valuable and could be pursued in future research.
- We agree about the potential impact of the equal division of carbon mass across bins. Unfortunately, to our knowledge, there is currently no information available regarding the assumption of mass distribution among organic matter species. This presents an excellent opportunity for follow-up research. Any data-backed updates on mass distribution can be easily incorporated into the Lambda-PFLOTRAN framework.
- In response to the sensitivity output, the Lambda-PFLOTRAN sensitivity is now being shown for Test Case 2 in Figure 5. Additional details have been added on lines 344 -350 as well.
- We have added our hypothesis for why the model is unsuccessful at capturing the total organic carbon concentrations on lines 353 - 355.
Citation: https://doi.org/10.5194/gmd-2024-34-AC3
-
AC3: 'Reply on RC2', Katherine Muller, 21 Aug 2024
-
AC1: 'Comment on gmd-2024-34', Katherine Muller, 21 Aug 2024
We thank the reviewers for their thoughtful review and comments. In response, we have made substantial improvements to the manuscript. Most notably, we have added a new section exploring the influence and variability of OM chemistry (Section 4- Variability and Impact of Organic Matter Speciation), which we believe has significantly strengthened our manuscript. Additional responses to specific reviewer comments are detailed below.
Citation: https://doi.org/10.5194/gmd-2024-34-AC1
Status: closed
-
RC1: 'Comment on gmd-2024-34', Anonymous Referee #1, 26 Jun 2024
General comments:
The presented work is a good contribution to improving the reactive transport modeling approach. It develops a representative reaction network based on substrate-explicit thermodynamic modeling and completes a biogeochemical simulation with a reactive transport model.
However, I think the current version cannot be considered for publication in GMD. There are some possible reasons:
(a) Poor model validation. A strict validation checks whether a new model approach reflects the actual or desired physical and biochemical behavior, considering using a reactive transport model. In this study, I noticed that there are very limited observations. For example, the work only involved limited observed O2 (Fig.4a), CO2 (Fig.5a), and total organic carbon datasets (Fig.5b). I also noticed that the model cannot capture observed total organic carbon (Fig.5b).
(b) Incomplete simulation experiments. The work just considered two test cases (oxygen depletion and respiration). Too small a test size may lack sufficient statistical power to detect new model capabilities.
(c) Poor analysis. The work shows limited biogeochemical analysis. An advantage of RTM is its ability to elucidate complex biogeochemical dynamics. For example, How do oxygen dynamics and total carbon influence the variations of pH and C, N, and P species over time (Fig.5a)? Although I understand this work is technical, the justification for these time-series variations should be carefully discussed.
Citation: https://doi.org/10.5194/gmd-2024-34-RC1 -
AC2: 'Reply on RC1', Katherine Muller, 21 Aug 2024
We understand the reviewer’s comment on the importance of model validation. We agree that more complete data for the test cases and additional datasets would be greatly beneficial, yet high-resolution organic matter characterization (e.g., FTICR-MS) is not commonly measured, which makes it not feasible to thoroughly validate the model given the lack of required data with both incubation and organic matter speciation. We hope that building this modeling framework will encourage others to measure organic matter speciation via FTICR-MS in the future. As more data becomes available, the Lambda-PFLOTRAN workflow can be further validated and refined.
We have clarified the scope of this paper on line 77, emphasizing that our work aims to detail a new capability and associated workflow for incorporation of organic matter characterization into reactive transport models. Model validation is important and new data collection could be designed to make it possible in future studies.
The illustrative model results presented here consider only aerobic respiration through the FTICR-MS informed reaction network. If desired, the PFLOTRAN input deck can be expanded and customized to include additional processes and full geochemistry.
We have added the following statement on lines 325-330 to capture this point:
“Three additional processes were added to describe the experimental conditions of Test Case 2 more accurately (i.e., release of carbon, nitrogen and sustained aerobic conditions); however, a PFLOTRAN input deck can be expanded and customized to include a host of additional processes and full geochemistry for a specific system of interest. For instance, aqueous complexation, mineral dissolution and precipitation, sorption, and redox reactions can be added, all of which can influence the resultant pH and carbon, nitrogen, and other nutrient dynamics.”
Citation: https://doi.org/10.5194/gmd-2024-34-AC2
-
AC2: 'Reply on RC1', Katherine Muller, 21 Aug 2024
-
RC2: 'Comment on gmd-2024-34', Anonymous Referee #2, 28 Jun 2024
This manuscript describes a new workflow for integrating detailed organic matter data from FTICR-MS into a reactive transport modeling framework for simulating biogeochemical interactions. The workflow includes sensitivity analysis and parameter estimation capabilities. I think this has the potential to be a highly valuable tool for connecting biogeochemical models with organic matter measurements at a level of chemical detail that has historically been a major challenge for biogeochemical model-data integration. The workflow covers many of the pain points for this type of model-data connection, including aggregating large numbers of compounds into a smaller, tractable set of representative compounds, assigning meaningful model parameters to those compounds, and converting the resulting reaction network into a reactive transport model simulation configuration automatically.
The manuscript describes the underlying assumptions and theoretical framework in clear terms and acknowledges important caveats such as the lack of quantitative information about the relative amounts of different compounds. The steps for using the workflow are also clearly described, which fits with the goal of making the workflow accessible and reproducible for other scientists. The two provided test cases are useful concrete examples of the workflow, important parameters, and useful data types for parameter estimation.
I did think that the results shown from the test cases did not fully make the case for how the FTICR-MS data add value to the resulting simulations. While the comparison of the full reaction network to simulations using multiple compounds does show a difference in overall predicted concentrations relative to a simplified reaction network with only one organic matter compound, it is difficult to see how the high resolution of chemical compounds, which is the hallmark of the FTICR-MS analysis, specifically contributes to the simulation results. It would be helpful to include a visualization of the distribution of compounds, how they translate to reaction networks with different properties, and how that ultimately affects the model results. For example, is the reaction network derived from the Test Case 1 FTICR-MS measurements meaningfully different from that of Test Case 2 in terms of its stoichiometry or distribution of reaction rates? Does the specific FTICR-MS data from each experiment yield better model results than FTICR-MS data from a different experiment? What if the null hypothesis was 10 compounds derived from a different set of measurements, rather than 1 generic compound? This comparison could more directly show the value of those measurements from a particular experiment.
Other comments:
Line 97: is yOM here the same as yOC in Equation 2?
Line 101: I think “1000” is missing from this sentence
Line 123: HS should be HS-
Line 192: Given the importance of specifying parameters to be estimated for the workflow, I think it would be helpful to list all the relevant parameters in one table, with explanations of what they mean and perhaps with suggestions for what kind of measurements could constrain them or what would be a good context for including them in the parameter estimation and what reasonable ranges of values might be.
Section 2.3.4: When discussing the parameter optimization, it would be helpful to start with some explanation of what kind of observations are useful for comparing with the model. While examples are provided in the test cases, I think some general explanation earlier on would be helpful as well.
Line 256: I’m curious how accurate the assumption of equal division of C mass across bins is. Do any data exist that could inform this? The approach might be quite sensitive to this assumption.
Figure 3: I was confused why k_deg and C_inhibit are shown in this figure, even though they aren’t part of the parameters being tested. And I did not understand the significance of the time dimension. Does this indicate whether model ensemble members are diverging or converging over the course of the simulation? It could use some more explanation of how to interpret that part.
Figure 5: The model generally fails to reproduce the observed total carbon concentrations. I might guess that the assumption of how much carbon gets solubilized from the solid state over time, or the characteristics of that soluble carbon, are not correct for this system. I would not necessarily expect a model to get every aspect of the system right, especially if it relates to a steady state process involving solid organic matter, which this framework is not designed to replicate. But I think this should be at least acknowledged and explained in the text and currently it feels like this part of the result is ignored.
Citation: https://doi.org/10.5194/gmd-2024-34-RC2 -
AC3: 'Reply on RC2', Katherine Muller, 21 Aug 2024
The following are responses and updates based on reviewer comments:
- We have revised Figure 3, the associated figure description, and companion text to help highlight the impact of OM speciation for Test Case 1 (lines 291-298).
- We have added a new section titled “Variability and Impact of Organic Matter Speciation” to address these comments. Section 4 explores the differences in OM chemistry between samples (Figure 5), the effects of OM chemistry on the lambda-derived reaction networks (Figure 6), and how lambda-derived reaction networks compare to an assumed generic CH2O bulk carbon species (Figure 7). Additionally, this section also examines the influence of OM chemistry through forward simulations, where only FTICR-MS input data is varied (Figure 8).
- We have updated the minor edits including: yOM, adding 1000, and fixing the subscript on HS-.
- Model parameters are defined in the Methods section (Section 2). While we understand the interest in gaining further insight into optimal model parameterization, we believe this is beyond the scope of the current manuscript. However, we agree that additional work focused on parameterization, experimental design, and data collection would be valuable and could be pursued in future research.
- We agree about the potential impact of the equal division of carbon mass across bins. Unfortunately, to our knowledge, there is currently no information available regarding the assumption of mass distribution among organic matter species. This presents an excellent opportunity for follow-up research. Any data-backed updates on mass distribution can be easily incorporated into the Lambda-PFLOTRAN framework.
- In response to the sensitivity output, the Lambda-PFLOTRAN sensitivity is now being shown for Test Case 2 in Figure 5. Additional details have been added on lines 344 -350 as well.
- We have added our hypothesis for why the model is unsuccessful at capturing the total organic carbon concentrations on lines 353 - 355.
Citation: https://doi.org/10.5194/gmd-2024-34-AC3
-
AC3: 'Reply on RC2', Katherine Muller, 21 Aug 2024
-
AC1: 'Comment on gmd-2024-34', Katherine Muller, 21 Aug 2024
We thank the reviewers for their thoughtful review and comments. In response, we have made substantial improvements to the manuscript. Most notably, we have added a new section exploring the influence and variability of OM chemistry (Section 4- Variability and Impact of Organic Matter Speciation), which we believe has significantly strengthened our manuscript. Additional responses to specific reviewer comments are detailed below.
Citation: https://doi.org/10.5194/gmd-2024-34-AC1
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 399 | 338 | 37 | 774 | 34 | 22 | 22 |
- HTML: 399
- PDF: 338
- XML: 37
- Total: 774
- Supplement: 34
- BibTeX: 22
- EndNote: 22
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1