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

Changes to the partitioning of biomass to coarse woody debris classes #1008

Merged
merged 14 commits into from
Jun 6, 2023

Conversation

adamhb
Copy link
Contributor

@adamhb adamhb commented Mar 29, 2023

Description:

This pull request is associated with discussion 962 which was continued in discussion 989.

The changes proposed here update the partitioning of biomass (struct + sapw) to coarse woody debris classes (cwd; aka cwd fuel classes). Fates currently sends some biomass to all fuel classes regardless of cohort dbh. For example it sends stem biomass to the 1,000 hr fuel class (logs > 7.6 cm in diameter) even when the cohort dbh is much less than 7.6 cm.

The changes proposed here update the cwd partitioning for small cohorts to be based on cohort dbh such that stem biomass is sent to the appropriate cwd class based on the Rothermel fuel class diameter thresholds (see comment from @jkshuman in discussion 962. The proportions of stem biomass sent to each flammable cwd class (twig, small branch, and large branch) is kept constant, where possible, to minimize changes in fuel properties. Partitioning for cohorts larger than 7.6 cm in dbh is not changed.

Changes are applied to all subroutines where stem biomass is partitioned into cwd classes: CWDInput (ED physiology), mortality_litter_fluxes (patch dynamics), fire_litter_fluxes (patch dynamics), logging_litter_fluxes (logging), and SendCohortToLitter (cohort dynamics).

Collaborators:

@lmkueppers @ZacharyRobbins @ckoven @rgknox @jkshuman @jenniferholm

Expectation of Answer Changes:

Changes should produce more flammable coarse woody debris (twig, small branch, and large branch) and less non-flammable cwd (1,000 hr fuel). These changes should be most apparent a) during early forest development when cohorts are < 7.6 cm in diameter and, b) for pfts that generally stay small (i.e. shrubs). Fuel properties such as SAV and bulk density are calculated as weighted averages of the fuel classes where the weights are the fraction of total fuel belonging to each class. Therefore, when the proportion of these classes change on the landscape, the aggregate fuel properties will also change, thereby changing fire behavior.

Checklist:

  • My change requires a change to the documentation.
  • I have updated the in-code documentation .AND. (the technical note .OR. the wiki) accordingly.
  • 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:

I have tested how these changes affect the accumulation of cwd, fuel properties, and fire behavior for a lower montane mixed conifer forest site in California's Sierra Nevada.

CTSM (or) E3SM (specify which) test hash-tag:
ctsm5.1.dev115 (commit#17e2acb6a)

CTSM (or) E3SM (specify which) baseline hash-tag:
ctsm5.1.dev115 (commit#17e2acb6a)

FATES baseline hash-tag:
sci.1.61.0_api.25.0.0 (commit#e663a6e6)

Test Output:

These changes result in a larger accumulation of flammable coarse woody debris classes (twig, small branch, and large branch) and less accumulation of non-flammable cwd (1,000 hr fuel) during early forest development (0-20 simulation yrs). Differences attenuate as the simulation progresses when more of the cwd is coming from larger cohorts (which are not affected by these changes). Although there is more potentially burnable fuel, the landscape is less flammable on any given day for the same FDI because SAV ratio and fuel bulk density move towards less flammable values due to a higher large branch to leaf fuel ratio on the landscape (fuel properties are weighted averages of the classes). These changes to fuel properties result in less burned area on average, and changes to fire line intensity and vary through time. Fire line intensity was the same in some decades and increased in some decades, likely due to more build up of flammable fuels (see slide 8 of results).

Results can be found here.

@glemieux glemieux requested a review from rosiealice April 10, 2023 18:31
Copy link
Contributor

@ZacharyRobbins ZacharyRobbins left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Feel free to email me with any questions!

biogeochem/FatesLitterMod.F90 Outdated Show resolved Hide resolved
biogeochem/FatesLitterMod.F90 Show resolved Hide resolved
biogeochem/FatesLitterMod.F90 Outdated Show resolved Hide resolved
@ZacharyRobbins
Copy link
Contributor

Adam, if it is not too much work, could you cluster the results on slide 6 and 7 by decade, and plot them on a similar scale? Then if people adopt this change they will be able to understand by how much they could anticipate the change being.
Thank you,
-Zachary

@adamhb
Copy link
Contributor Author

adamhb commented May 2, 2023

Thanks Zachary for suggesting the decadal-scale figures which can be found on slide 8 of the results here

A new insight from these figures is that mean fire line intensity increases in some decades with the proposed changes and in other decades its very similar. Differences are likely due to having more combustible fuel on the landscape. (I've updated the "Test results" section on the PR above to reflect this insight).

@glemieux glemieux self-assigned this Jun 5, 2023
Copy link
Contributor

@rgknox rgknox left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

looks good

@glemieux
Copy link
Contributor

glemieux commented Jun 6, 2023

Regression testing on cheyenne is complete. All expected tests pass with multiple DIFFs. The set of differing variables appear to make sense given the changes made here, i.e. no unexpected variables not related to fire and litter. Note that some testmods do not have litter or fire variables reported and are reporting as passing b4b.

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

@glemieux glemieux merged commit fce15a9 into NGEET:main Jun 6, 2023
@adamhb adamhb deleted the cwd_partitioning branch July 7, 2023 18:09
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.

4 participants