Skip to content
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

Add mortality and productivity aboveground carbon fluxes as history variables #926

Merged
merged 9 commits into from
Feb 2, 2023

Conversation

JessicaNeedham
Copy link
Contributor

@JessicaNeedham JessicaNeedham commented Oct 21, 2022

Description:

This PR adds two new history variables, FATES_ABOVEGROUND_MORT_SZPF and FATES_ABOVEGROUND_PROD_SZPF - aboveground mortality and productivity carbon fluxes. These history variables will allow direct comparisons with the results from Piponiot et al. 2022, and potentially other forest plot data.

https://nph.onlinelibrary.wiley.com/doi/10.1111/nph.17995

Piponiot et al. output AWP and AWM - which they describe as “aboveground woody productivity (AWP), the flux in AGB associated with tree growth and recruitment) and aboveground woody mortality (AWM, the flux from AGB to necromass due to mortality)”.

I think woody here refers to woody plants, rather than woody parts. AWP and AWM are described as AGB fluxes, and AGB is calculated using equations from Chave et al. 2014, and others, which I believe includes leaves (although I’m not totally sure after reading it again). Leaves are probably a very small part of biomass anyway (less than 5% according to Chave et al. 2014). To be consistent with AGB calculations in FATES, I included leaves but not reproductive tissues. Although I'm curious why reproductive tissue isn't included in AGB....

I wanted these variables to be FATES_ABOVEGROUND_MORTALITY_CFLUX_SZPF and FATES_ABOVEGROUND_PRODUCTIVITY_CFLUX_SZPF to be consistent with other history variables but the names were too long and the HLM failed to recognise them. I’m open to suggestions on alternatives.

Results here are based on a 400 year run at BCI. New variables are compared with existing flux variables and with results from Piponiot et al.

Collaborators:

@ckoven @jenniferholm

Expectation of Answer Changes:

No changes expected.

Checklist:

  • My change requires a change to the documentation.
  • I have updated the in-code documentation .AND. (the technical note .OR. the wiki) accordingly.
  • [x ] I have read the CONTRIBUTING document.
  • FATES PASS/FAIL regression tests were run
  • If answers were expected to change, evaluation was performed and provided

Test Results:

CTSM (or) E3SM (specify which) test hash-tag:

CTSM (or) E3SM (specify which) baseline hash-tag:

FATES baseline hash-tag:

Test Output:
Fluxes_by_size_v_BCI

Fluxes_through_time

Jessica F Needham and others added 6 commits October 14, 2022 17:19
aboveground woody productivity for comparison of FATES with
forest plot data.
[ Change aboveground to agw (aboveground woody) in history variable
names ]

Fixes: [NGT-ED Github issue #]

User interface changes?: [Yes (describe what changes), No]

Code review: [Names]

Test suite: [suite name, machine, compilers]
Test baseline:
Test namelist changes:
Test answer changes: [bit for bit, roundoff, climate changing]
Test summary:no testing
[ Since AWM and AWP are fluxes of AGB, and AGB includes leaves,
add leaves to AWP and AWM calculations. Tidy up names as a result
of this change.  ]

Fixes: [NGT-ED Github issue #]

User interface changes?: [Yes (describe what changes), No]

Code review: [Names]

Test suite: [suite name, machine, compilers]
Test baseline:
Test namelist changes:
Test answer changes: [bit for bit, roundoff, climate changing]
Test summary: no testing
[ ]

Fixes: [NGT-ED Github issue #]

User interface changes?: [Yes (describe what changes), No]

Code review: [Names]

Test suite: [suite name, machine, compilers]
Test baseline:
Test namelist changes:
Test answer changes: [bit for bit, roundoff, climate changing]
Test summary:no testing
[  ]

Fixes: [NGT-ED Github issue #]

User interface changes?: [Yes (describe what changes), No]

Code review: [Names]

Test suite: [suite name, machine, compilers]
Test baseline:
Test namelist changes:
Test answer changes: [bit for bit, roundoff, climate changing]
Test summary: no testing
[ To be consistent with AGB reproductive carbon isn't included
so no need to add it here.  ]

Fixes: [NGT-ED Github issue #]

User interface changes?: [Yes (describe what changes), No]

Code review: [Names]

Test suite: [suite name, machine, compilers]
Test baseline:
Test namelist changes:
Test answer changes: [bit for bit, roundoff, climate changing]
Test summary: no testing
Conflicts:
	main/FatesHistoryInterfaceMod.F90
	main/FatesRestartInterfaceMod.F90
@glemieux glemieux self-assigned this Jan 30, 2023
@glemieux
Copy link
Contributor

glemieux commented Feb 2, 2023

Regression testing on cheyenne is complete. I tested this and #971 at the same time by merging them together. All expected tests pass b4b against fates-sci.1.63.0_api.25.1.0-ctsm5.1.dev116 with the exception of field list differences.

Folder location: /glade/u/home/glemieux/scratch/ctsm-tests/tests_pr971_926-dev116

@glemieux glemieux merged commit bcddaf6 into NGEET:main Feb 2, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants