H2 effects on CH4 lifetime #758
Open
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
…rs in the gcam coupling
…ead in from teh ini file
Differences in Hector outputsHello, this is The current pull request's outputs do not differ from 3.1.1 (d931a00). |
|
Some specific questions for @ssmithClimate & @bpbond
|
kdorheim
commented
Nov 26, 2024
DESCRIPTION
Outdated
There was a problem hiding this comment.
@bpbond just a heads up I did increase this for my own sanity purposes related to the gcam-hector coupling...
kdorheim
commented
Nov 26, 2024
| @@ -165,7 +165,7 @@ void CH4Component::run(const double runToDate) { | |||
| << std::endl; | |||
|
|
|||
| // Permafrost thaw produces CH4 emissions | |||
| #define PG_C_TO_TG_CH4 (1000.0 * 16.04 / 12.01) | |||
| const double PG_C_TO_TG_CH4 (1000.0 * 16.04 / 12.01); | |||
There was a problem hiding this comment.
@bpbond is this the correct usage of the const double?
kdorheim
commented
Nov 26, 2024
| H2_emissions.get(H2_emissions.firstdate()).value(U_TG_H2)); | ||
|
|
||
|
|
||
| toh = a + b + c + d + e; |
There was a problem hiding this comment.
@bpbond do we have strong feelings on the use of e as a variable name?
Differences in Hector outputsHello, this is The current pull request's outputs do not differ from 3.1.1 (d931a00). |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
NOTE: this PR is targeting the dev-h2 branch
This branch implements the H2 effects on CH4, it will close #749. We have added equation (5) into Hector to account for how$H_2$ interacts with OH, decreasing the amount of OH that can interact with $[CH_4]$ . We derived the value for the $c_{H_2}$ from Bertagni et al. 2022.
Since we don't have$H_2$ emissions for the historical or future periods, we shouldn't see any changes in Hector results (Leeyabot or the old new tests). In order to get a feel for how this development affects Hector results I've included some results from emission impulse runs.
This first figure shows results from three Hector runs, a constant emission-driven run, and then two impulse runs varying in magnitude. It is intended to show how this model development cascades throughout the model. The pulse of$H_2$ emissions causes $\tau_{OH}$ (methane's lifetime with respect to OH) to change, which causes $[CH_4]$ to change, which is what we were expecting. Since in Hector the trop. $O_3$ burden and strat. $H_2O$ vapor radiative forcing computations are proportional to $[CH_4]$ the $H_2$ emissions pulse affects these variables as well.
This second figure demonstrates the sensitivity of impulse response to our new parameter$c_{H_2}$ . As we would hope, increasing and decreasing this parameter changes the system response.
In this third figure, we compare the RF response from unit emission of$CO_2$ and $H_2$ 100 years after the emissions impulse, which can be used to calculate the GWP100. Hector's GWP100 for H2 ~ 5.17 is consistent with the Sand et al. results included in Supplementary Table 6. Sand et al. found the GWP100 of hydrogen due to methane ranges from 4.9 to 5.7 with a model mean of 5.1, which is not too far off what we are getting from Hector at the moment!
References
Bertagni, M. B., Pacala, S. W., Paulot, F., & Porporato, A. (2022). Risk of the hydrogen economy for atmospheric methane. Nature Communications, 13(1), 7706.
Sand, M., Skeie, R. B., Sandstad, M., Krishnan, S., Myhre, G., Bryant, H., Derwent, R., Hauglustaine, D., Paulot, F., Prather, M., & Stevenson, D. (2023). A multi-model assessment of the Global Warming Potential of hydrogen. Communications Earth & Environment, 4(1), 1–12.
Materials Used
PR-758.zip