ctsm5.4.004_bgc_f09_131_ [AD, SASU, pSASU, HIST]

[wwieder]

Description:
CLM 5.4 release tag + changes in CESM alpha08a tag on ne30 grid

• Crop model parameter cals, GDD cals, & GDD baselines from @samsrabin, which reduces wheat LAI (among others)
• Corrected CMIP7 CO2 isotope streams @ekluzek
• Boreal fire tunings suggested by @lifang0209
• Minor parameter changes (leaf angle, soil C) from @linniahawkins
• RTTMGP LAI bug fix from @olyson
• Revert to use Sturm over glaciers as requested by the Land Ice team

LDF plots vs. 124

---

Case directory:
Locally:
/glade/u/home/wwieder/CTSM (ctsm5.4.004*)

---

Sandbox:
Locally:
/glade/work/wwieder/ctsm5.4.0/

---

Case directory:
Locally:

• /glade/work/wwieder/ctsm5.4.0/ctsm5.4.004_bgc_f09_131_AD
• /glade/work/wwieder/ctsm5.4.0/ctsm5.4.004_bgc_f09_131_SASU
• /glade/work/wwieder/ctsm5.4.0/ctsm5.4.004_bgc_f09_131_pSASU
• /glade/work/wwieder/ctsm5.4.0/ctsm5.4.004_bgc_f09_131_HIST

./create_newcase --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.004_bgc_f09_131_AD --compset I1850Clm60BgcCropCrujra --res f09_t232 --run-unsupported

./create_clone --clone /glade/work/wwieder/ctsm5.4.0/ctsm5.4.004_bgc_f09_131_AD --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.004_bgc_f09_131_SASU --keepexe

./create_clone --clone /glade/work/wwieder/ctsm5.4.0/ctsm5.4.004_bgc_f09_131_AD --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.004_bgc_f09_131_pSASU --keepexe

./create_newcase --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.004_bgc_f09_131_HIST --compset IHISTClm60BgcCropCrujra --res f09_t232 --run-unsupported

---

1. AD

xml changes

./xmlchange MOSART_MODE=NULL
./xmlchange STOP_N=80,STOP_OPTION=nyears
./xmlchange CLM_ACCELERATED_SPINUP=on,CLM_FORCE_COLDSTART=on
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=3

2. SASU

./xmlchange MOSART_MODE=NULL
./xmlchange RUN_TYPE=hybrid
./xmlchange STOP_N=60,STOP_OPTION=nyears
./xmlchange CLM_ACCELERATED_SPINUP=sasu,CLM_FORCE_COLDSTART=off
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=5
./xmlchange GET_REFCASE=TRUE
./xmlchange RUN_REFCASE=ctsm5.4.004_bgc_f09_131_AD
./xmlchange RUN_REFDATE=0361-01-01
./xmlchange DRV_RESTART_POINTER=rpointer.cpl.0361-01-01-00000
./xmlchange RUN_REFDIR=/glade/derecho/scratch/wwieder/archive/ctsm5.4.004_bgc_f09_131_AD/rest/0361-01-01-00000

3. pSASU

./xmlchange MOSART_MODE=NULL
./xmlchange STOP_N=60,STOP_OPTION=nyears
./xmlchange CLM_ACCELERATED_SPINUP=off,CLM_FORCE_COLDSTART=off
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=2
./xmlchange RUN_TYPE=hybrid
./xmlchange GET_REFCASE=TRUE
./xmlchange RUN_REFDATE=0381-01-01
./xmlchange DRV_RESTART_POINTER=rpointer.cpl.0381-01-01-00000
./xmlchange RUN_REFCASE=ctsm5.4.004_bgc_f09_131_SASU
./xmlchange RUN_REFDIR=/glade/derecho/scratch/wwieder/archive/ctsm5.4.004_bgc_f09_131_SASU/rest/0381-01-01-00000

4. HIST
   user_nl_datm_streams

!! manually update datm_streams to get through 2023
co2tseries.20tr:year_last = 2024
co2tseries.20tr:datafiles = /glade/campaign/cesm/cesmdata/inputdata/atm/datm7/CO2/fco2_datm_globalSSP3-7.0_simyr_1750-2501_CMIP6_c201101.nc

! No CMIP6 aero files? 
! We'll just use default CMIP6 HIST files and cycle after 2014 
!presaero.hist:datafiles = /glade/campaign/cesm/cesmdata/inputdata/atm/cam/chem/trop_mozart_aero/aero/aerosoldep_rcp4.5_monthly_1849-2104_0.9x1.25_c100407.nc
!presaero.hist:year_last = 2024

user_nl_clm

finidat = '/glade/derecho/scratch/wwieder/archive/ctsm5.4.004_bgc_f09_131_pSASU/rest/0161-01-01-00000/ctsm5.4.004_bgc_f09_131_pSASU.clm2.r.0161-01-01-00000.nc'

use_init_interp = .true.
init_interp_fill_missing_urban_with_HD = .true.
ndep_taxmode = 'extend'

! "Survivability" output
hist_dov2xy(1) = .true.
hist_dov2xy(2) = .false.
hist_fields_list_file = .false.
hist_fincl2 = 'TLAI', 'TSA', 'TREFMNAV', 'TREFMXAV', 'BTRANMN', 'Vcmx25Z', 'FSH', 'VEGWP', 'FCTR', 'FCEV', 'FGEV',
        'FIRE', 'FSR', 'FIRA', 'FSA', 'GSSUNLN', 'GSSHALN', 'TSKIN', 'GPP', 'NPP', 'AGNPP', 'TOTVEGC',
        'NPP_NUPTAKE', 'AR', 'HR', 'HTOP'
hist_mfilt(1) = 1
hist_mfilt(2) = 1
hist_nhtfrq(1) = 0
hist_nhtfrq(2) = 0
hist_type1d_pertape(1) = ' '
hist_type1d_pertape(2) = ' '

! Crop output
hist_dov2xy(3) = .false.
hist_mfilt(3) = 1
hist_nhtfrq(3) = -8760
hist_type1d_pertape(3) = ' '
hist_fincl3 = 'SDATES:I', 'SDATES_PERHARV:I', 'SYEARS_PERHARV:I', 'HDATES:I', 'GRAINC_TO_FOOD_PERHARV:I', 'GRAINC_TO_FOOD_ANN:I', 'GRAINN_TO_FOOD_PERHARV:I', 'GRAINN_TO_FOOD_ANN:I', 'GRAINC_TO_SEED_PERHARV:I', 'GRAINC_TO_SEED_ANN:I', 'GRAINN_TO_SEED_PERHARV:I', 'GRAINN_TO_SEED_ANN:I', 'HDATES:I', 'GDDHARV_PERHARV:I', 'GDDACCUM_PERHARV:I', 'HUI_PERHARV:I', 'SOWING_REASON_PERHARV:I', 'HARVEST_REASON_PERHARV:I', 'SWINDOW_STARTS:I', 'SWINDOW_ENDS:I', 'GDD20_BASELINE:I', 'GDD20_SEASON_START:I', 'GDD20_SEASON_END:I'

1850-1871

./xmlchange RUN_TYPE=startup
./xmlchange RUN_STARTDATE=1850-01-01
./xmlchange STOP_OPTION=nyears
./xmlchange STOP_N=21
./xmlchange CONTINUE_RUN=FALSE
./xmlchange RESUBMIT=0
./xmlchange DATM_YR_ALIGN=1901
./xmlchange DATM_YR_START=1901
./xmlchange DATM_YR_END=1920

1871-1901

./xmlchange STOP_N=30
./xmlchange CONTINUE_RUN=TRUE

1901-2023

./xmlchange STOP_N=31
./xmlchange DATM_YR_END=2023
./xmlchange RESUBMIT=3
./xmlchange DATM_SKIP_RESTART_READ=TRUE

---

SourceMods:
None

---

Diagnostics:
Diags (if available)
LDF plots

ILAMB diagnostics

---

Output:
Output (if still available):
/glade/derecho/scratch/wwieder/archive/ctsm5.4.004_bgc_f09_131_HIST/lnd

---

Contacts:
@wwieder

---

Extra details:

[wwieder]

Survival plots here, no real changes for the native PFTs compared to #124, as intended.

[wwieder]

Haven't plotted Vcmax in a while... here's where things stand. I guess we should have a look at what results looked with with CLM5? Who can help comment on the reasonableness of these values?

[linniahawkins]

Vcmax for the trees look reasonable to me, here is the distribution of values I have from the TRY database:

I don't have numbers for grasses but the relative values are as I would expect (highest for C3 in the Arctic, lowest for C4).

[linniahawkins]

Survivability looks very similar to cal115 to my eye. No concerns there.

[wwieder]

@samsrabin low urgency here, but wondering if this change in global food production is expected based on what you've seen? h2 files have crop outputs that you requested a while back. (note plots using climo data are from the last 20 years of the HIST case, despite what labels say on the ADF pages)

[wwieder]

@lifang0209, please have a look at burned area in these simulations when you're able. Simulations look pretty similar to simulations in #124, which is likely by design. Color bar here isn't too useful, but raw data are available on scratch (see above).

[wwieder]

@linniahawkins seems like changes to soil C turnover gave us a small increase in Arctic LAI. Still looks good OK to me, but I thought it's worth noting. LAI changes over agricultural regions reflect crop model changes (mostly for wheat). Will get ILAMB run on these too, but please see if LDF looks OK when you're able.

[wwieder]

@olyson, no rush on this, but can you add this case to the ILAMB plots you'd made for #124?

[olyson]

@olyson, no rush on this, but can you add this case to the ILAMB plots you'd made for #124?

ILAMB results are here

[samsrabin]

@samsrabin low urgency here, but wondering if this change in global food production is expected based on what you've seen? h2 files have crop outputs that you requested a while back. (note plots using climo data are from the last 20 years of the HIST case, despite what labels say on the ADF pages)

Yes, that's consistent with what I would expect. Compare the dark red/brown line here (before my new parameters) to the blue line (with my new parameters). The reduced yields of wheat and to a lesser extent rice are likely driving the difference seen in your plot.

[linniahawkins]

@linniahawkins seems like changes to soil C turnover gave us a small increase in Arctic LAI. Still looks good OK to me, but I thought it's worth noting. LAI changes over agricultural regions reflect crop model changes (mostly for wheat). Will get ILAMB run on these too, but please see if LDF looks OK when you're able.

Looks as expected, the soil carbon stocks came down from ~1650 to 1450 PgC and the increase in arctic LAI is as we expected. I skimmed through the ILAMB results and things look good there too. Thanks!