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

Coupled Sturm_vegMods (f2000.ne30_t232.BGC_CROP_Sturm_vegMods) #25

Closed
wwieder opened this issue Oct 17, 2023 · 12 comments
Closed

Coupled Sturm_vegMods (f2000.ne30_t232.BGC_CROP_Sturm_vegMods) #25

wwieder opened this issue Oct 17, 2023 · 12 comments

Comments

@wwieder
Copy link
Collaborator

wwieder commented Oct 17, 2023

Description:
Let's try to do everything we can to keep vegetation alive: Combining Sturm snow, phenology triggers, FUNfracfixers, and frootleaf (#2, #4, #7, and #15) in an F-case (similar to #13.)


Case directory:
/glade/u/home/slevis/cases_LMWG_dev/f2000.ne30_t232.BGC_CROP_Sturm_vegMods


Sandbox:
/glade/work/hannay/cesm_tags/cesm2_3_alpha16b


user_nl_ changes:

hist_fincl1 = 'TOTECOSYSC', 'TOTECOSYSN', 'TOTSOMC', 'TOTSOMN', 'TOTVEGC', 'TOTVEGN', 'TLAI', 'GPP', 'CPOOL', 'NPP', 'TWS', 'H2OSNO','NDEP_TO_SMINN'
use_c13 = .true.
use_c14 = .true.
use_c13_timeseries = .true.
use_c14_bombspike = .true.
fsurdat = '/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_ne30np4.pg3_SSP5-8.5_78pfts_CMIP6_1850-2100_c230227.nc'
finidat = '/glade/scratch/oleson/archive/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_SturmSnowtk_1850pAD/rest/0401-01-01-00000/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_SturmSnowtk_1850pAD.clm2.r.0401-01-01-00000.nc'
paramfile = '/glade/p/cgd/tss/people/oleson/modify_param/ctsm51_params.c211112.nc.ABCD.nc'
use_init_interp = .true.
reseed_dead_plants = .true.

Also user_nl_cam:

fincl2 = 'TREFHT'
nhtfrq = 0,-24
mfilt  = 1,1

SourceMods:
src.clm/atm2lndType.F90: Add some additional history fields
src.clm/clm_varpar.F90: nlevurb=5 -> nlevurb=10
src.clm/surfrdMod.F90: Remove read of PFTDATA_MASK
src.clm/Wateratm2lndType.F90: Add some additional history fields
src.clm/SoilTemperatureMod.F90 Sturm Snow tk
src.clm/CNPhenologyMod.F90 decreasing temperature thresholds to -3C, reducing the minimum critical daylength by 25%, doubling the snow threshold (I think to 20 cm)


Diagnostics:
https://webext.cgd.ucar.edu/F2000climo/f2000.ne30_t232.BGC_CROP_Sturm_vegMods/lnd/f2000.ne30_t232.BGC_CROP_Sturm_vegMods.77_96-f2000.ne30_t232.BGC_CROP.77_96/setsIndex.html


Output:
(if still available):
/glade/scratch/slevis/archive/f2000.ne30_t232.BGC_CROP_Sturm_vegMods/


Contacts:
@slevis-lmwg


Extra details:
./create_clone --clone /glade/u/home/slevis/cases_LMWG_dev/f2000.ne30_t232.BGC_CROP_SturmSnowtk --case /glade/u/home/slevis/cases_LMWG_dev/f2000.ne30_t232.BGC_CROP_Sturm_vegMods

Recommended mods #6 (comment) if also generating a new control, so I assume not for now.

@olyson
Copy link
Collaborator

olyson commented Oct 17, 2023

Seems like getting initial conditions from #24 is a good idea. That simulation is underway.

@olyson
Copy link
Collaborator

olyson commented Oct 31, 2023

@slevis-lmwg , since you ran the other F-cases, would you be able to run this one?
I have a spunup (4.94% disequilibrium) initial file here:
/glade/p/cgd/tss/people/oleson/CLM5_restarts/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_ABCD_1850pAD.clm2.r.0801-01-01-00000.nc

@slevis-lmwg slevis-lmwg self-assigned this Oct 31, 2023
@slevis-lmwg slevis-lmwg changed the title Coupled Sturm_vegMods (enter the name of the run here) Coupled Sturm_vegMods (f2000.ne30_t232.BGC_CROP_Sturm_vegMods) Oct 31, 2023
@wwieder
Copy link
Collaborator Author

wwieder commented Nov 3, 2023

Would ILAMB results on the 3 F-cases (#11, #14 #25) also be helpful here? Is this possible w/o a transient run?

@slevis-lmwg
Copy link
Collaborator

I plan on discussing iLamb questions (from here and the sp cases) with Keith on Weds. morning.

@wwieder
Copy link
Collaborator Author

wwieder commented Nov 14, 2023

Maybe ILAMB isn't super critical, but what about diagnostics for this case, e.g. vs. Control F-case (#14) or vs. Sturm F-case (#11)? Or polar plots of summer ELAI from these three cases with information on % survival?

@olyson
Copy link
Collaborator

olyson commented Nov 14, 2023

@slevis-lmwg , is this simulation complete? If so, I can generate some diagnostics.

@slevis-lmwg
Copy link
Collaborator

Just posted.

@wwieder
Copy link
Collaborator Author

wwieder commented Nov 14, 2023

Thanks, Sam. Can one you you also look at the summer LAI plots and (dead area for the control, STURM and SturmABCD f-cases?).

@olyson
Copy link
Collaborator

olyson commented Nov 14, 2023

I'll do that.

@wwieder
Copy link
Collaborator Author

wwieder commented Nov 14, 2023

Vegetation feedbacks with accelerated snowmelt are stronger in the F case. You can flip between this plot for the I-case vs. this one for the F-case

Permafrost losses from Sturm don't seem any worse in the F-case than they were in the I-case, although it's kind of hard to compare...

@olyson
Copy link
Collaborator

olyson commented Nov 14, 2023

Polar plots of boreal summer TLAI for the three F-cases. Top is control. Middle is Sturm only. Bottom is equivalent to the land-only ABCD. Improvements similar to those seen in the land-only simulations although not as pronounced. The improvements in dead vegetation compared to the control are:
(0) Percent Decrease in Sturm dead_veg (km2): -25.895
(0) Percent Decrease in ABCD dead_veg (km2): -45.1348

f2000 ne30_t232 BGC_CROP JJA TLAI polar
f2000 ne30_t232 BGC_CROP_SturmSnowtk JJA TLAI polar
f2000 ne30_t232 BGC_CROP_Sturm_vegMods JJA TLAI polar

Annual cycle plots with the F-cases designated by "_F".:
annual_cycle_new_Siberia1_1850

annual_cycle_radiation_new_Siberia1_1850

@wwieder
Copy link
Collaborator Author

wwieder commented Nov 14, 2023

It's fun to click through the maps in different windows and watch the arctic come alive!
Thanks for posting these @olyson and doing the runs, @slevis-lmwg.

The biophysical changes for snowmelt seem pretty subtle compared to the LAI and survival differences in the Sturm vs. Sturm-ABCD runs (blue vs. yellow)

@olyson olyson closed this as completed Oct 9, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

3 participants