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.

Sign up for GitHub

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

Disable hydraulic failure mortality for leafless plants #1130

Merged
merged 6 commits into from
Jan 30, 2024
Merged

Disable hydraulic failure mortality for leafless plants #1130

Changes from 2 commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
da1b798
Prevent hydraulic failure mortality to occur in deciduous plants when…
mpaiao Dec 6, 2023
8756dc1
Fix typo.
mpaiao Dec 6, 2023
53b88e2
Edited hydraulic failure mortality so deciduousness prevents mortalit…
mpaiao Jan 24, 2024
6120cad
Merge tag 'sci.1.71.2_api.33.0.0' into mpaiao-pr-hfm+phenology
glemieux Jan 29, 2024
bee1769
implement @mpaiao suggested fix for deconflict
glemieux Jan 29, 2024
1777702
remove hmort setting missed in previous deconflict
glemieux Jan 29, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
67 changes: 45 additions & 22 deletions biogeochem/EDMortalityFunctionsMod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -13,6 +13,9 @@ module EDMortalityFunctionsMod
use EDTypesMod , only : ed_site_type
use EDParamsMod, only : maxpft
use FatesConstantsMod , only : itrue,ifalse
use FatesConstantsMod , only : ihard_stress_decid
use FatesConstantsMod , only : isemi_stress_decid
use FatesConstantsMod , only : leaves_off
use FatesAllometryMod , only : bleaf
use FatesAllometryMod , only : storage_fraction_of_target
use FatesInterfaceTypesMod , only : bc_in_type
Expand All @@ -25,6 +28,7 @@ module EDMortalityFunctionsMod

use PRTGenericMod, only : carbon12_element
use PRTGenericMod, only : store_organ
use PRTParametersMod , only : prt_params
use shr_log_mod , only : errMsg => shr_log_errMsg

implicit none
Expand Down Expand Up @@ -89,14 +93,27 @@ subroutine mortality_rates( cohort_in,bc_in,btran_ft, mean_temp, &
real(r8) :: min_fmc_ar ! minimum fraction of maximum conductivity for absorbing root
real(r8) :: min_fmc ! minimum fraction of maximum conductivity for whole plant
real(r8) :: flc ! fractional loss of conductivity
logical :: is_decid_dormant ! Flag to signal that the cohort is deciduous and dormant


real(r8), parameter :: frost_mort_buffer = 5.0_r8 ! 5deg buffer for freezing mortality
logical, parameter :: test_zero_mortality = .false. ! Developer test which
! may help to debug carbon imbalances
! and the like

! Check if the PFT is deciduous and leaves are completely abscised. If this is the case,
! we prevent hydraulic failure mortality to occur as plants are leafless. For now,
! semi-deciduous plants with senescing leaves may still die of hydraulic failure, but in
! the future we could accelerate senescence to avoid mortality. Note that both drought
! deciduous and cold deciduous are considered here to be consistent with the idea that
! plants without leaves cannot die of hydraulic failure.
is_decid_dormant = & !
( prt_params%stress_decid(cohort_in%pft) == ihard_stress_decid .or. & ! Drought deciduous
prt_params%stress_decid(cohort_in%pft) == isemi_stress_decid .or. & ! Semi-deciduous
prt_params%season_decid(cohort_in%pft) == itrue ) .and. & ! Cold deciduous
( cohort_in%status_coh == leaves_off ) ! ! Fully abscised
Copy link
Contributor

Choose a reason for hiding this comment

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

I think this looks good. I started to muse about what guidance we could offer in the future if we ever enable a long window of time for leaf drop, but that would be a full discussion and could not really be encapsulated in a couple of comments here.



! Size Dependent Senescence
! Size Dependent Senescence
! rate (r) and inflection point (ip) define the increase in mortality rate with dbh
mort_r_size_senescence = EDPftvarcon_inst%mort_r_size_senescence(cohort_in%pft)
mort_ip_size_senescence = EDPftvarcon_inst%mort_ip_size_senescence(cohort_in%pft)
Expand Down Expand Up @@ -139,29 +156,35 @@ subroutine mortality_rates( cohort_in,bc_in,btran_ft, mean_temp, &
bmort = EDPftvarcon_inst%bmort(cohort_in%pft)

! Proxy for hydraulic failure induced mortality.
hf_sm_threshold = EDPftvarcon_inst%hf_sm_threshold(cohort_in%pft)
hf_flc_threshold = EDPftvarcon_inst%hf_flc_threshold(cohort_in%pft)
if(hlm_use_planthydro.eq.itrue)then
!note the flc is set as the fraction of max conductivity in hydro
min_fmc_ag = minval(cohort_in%co_hydr%ftc_ag(:))
min_fmc_tr = cohort_in%co_hydr%ftc_troot
min_fmc_ar = minval(cohort_in%co_hydr%ftc_aroot(:))
min_fmc = min(min_fmc_ag, min_fmc_tr)
min_fmc = min(min_fmc, min_fmc_ar)
flc = 1.0_r8-min_fmc
if(flc >= hf_flc_threshold .and. hf_flc_threshold < 1.0_r8 )then
hmort = (flc-hf_flc_threshold)/(1.0_r8-hf_flc_threshold) * &
EDPftvarcon_inst%mort_scalar_hydrfailure(cohort_in%pft)
else
hmort = 0.0_r8
endif
! This is only applied on plants that have leaves. Deciduous plants in the
! "leaves off" season cannot die of hydraulic failure.
decid_check: if (is_decid_dormant) then
hmort = 0.0_r8
else
if(btran_ft(cohort_in%pft) <= hf_sm_threshold)then
hmort = EDPftvarcon_inst%mort_scalar_hydrfailure(cohort_in%pft)
hf_sm_threshold = EDPftvarcon_inst%hf_sm_threshold(cohort_in%pft)
hf_flc_threshold = EDPftvarcon_inst%hf_flc_threshold(cohort_in%pft)
if(hlm_use_planthydro.eq.itrue)then
!note the flc is set as the fraction of max conductivity in hydro
min_fmc_ag = minval(cohort_in%co_hydr%ftc_ag(:))
min_fmc_tr = cohort_in%co_hydr%ftc_troot
min_fmc_ar = minval(cohort_in%co_hydr%ftc_aroot(:))
min_fmc = min(min_fmc_ag, min_fmc_tr)
min_fmc = min(min_fmc, min_fmc_ar)
flc = 1.0_r8-min_fmc
if(flc >= hf_flc_threshold .and. hf_flc_threshold < 1.0_r8 )then
hmort = (flc-hf_flc_threshold)/(1.0_r8-hf_flc_threshold) * &
EDPftvarcon_inst%mort_scalar_hydrfailure(cohort_in%pft)
else
hmort = 0.0_r8
endif
else
hmort = 0.0_r8
if(btran_ft(cohort_in%pft) <= hf_sm_threshold)then
hmort = EDPftvarcon_inst%mort_scalar_hydrfailure(cohort_in%pft)
else
hmort = 0.0_r8
endif
endif
endif
end if decid_check

! Carbon Starvation induced mortality.
if ( cohort_in%dbh > 0._r8 ) then
Expand Down