wetlands_fire.n01.ctsm5.4.029_f09noCrop_HH_159_ [AD, SASU, pSASU]

[wwieder]

Description:
Hillslope hydrology 16 column run with fire modifications from coldstart

• f09, Hillslope, noCrop w/
• WIEMIP surface data and DATM
• Code mods suggested for n02 tag as in /glade/u/home/fangli/will2026/hh-test/SourceMods/src.clm that are [included in this PR] (branch_tags/wetlands.n02.ctsm5.4.029: updates for calculating fire with HH ESCOMP/CTSM#3934)

---

Case directory:
Locally:
/glade/work/wwieder/wiemip/wetlands_fire.n01.ctsm5.4.029_f09noCrop_HH_159_*

---

Sandbox:
Locally:
/glade/u/home/wwieder/CTSM (wetlands_fire)
branch_tags/wetlands.n01.ctsm5.4.029-1-g32048e5cb

---

./create_newcase --case /glade/work/wwieder/wiemip/wetlands_fire.n01.ctsm5.4.029_f09noCrop_HH_159_AD --compset I1850Clm60BgcCrujra --res f09_t232 --run-unsupported --user-mods-dirs /glade/u/home/wwieder/CTSM/cime_config/usermods_dirs/clm/WIEMIP

user_nl_ changes:

fsurdat = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/surfdata_esmf/ctsm5.4.0/surfdata_0.9x1.25_hist_1850_78pfts_WIEMIP_c260408.nc'
paramfile = '/glade/work/wwieder/modified_inputs/ctsm60_params.genericCropFix.ch4cal_c.c260410.nc'
use_c13 = .false.
use_c14 = .false.

! 16 column hillslope and vegetation data
hillslope_file = '/glade/derecho/scratch/swensosc/WIEMIP/hillslope_geo_params_trapezoid_worldcover_0.9x1.25_Global_FAB.nc'

! 9 column hillslope file
!hillslope_file = '/glade/work/swensosc/HillslopeHydrology/InputData/Global/hillslope_geo_params_3_col_3_aspect_trapezoid_worldcover_0.9x1.25_Global_FAB.nc'

use_hillslope = .true.
use_hillslope_routing = .true.

! with bgc active
finundation_method = 'h2osfc'

hillslope_fsat_equals_zero = .true.
n_dom_pfts = 2

---

SourceMods:

---

1. AD

./xmlchange MOSART_MODE=NULL
./xmlchange DATM_YR_START=1850,DATM_YR_END=1869
./xmlchange CLM_ACCELERATED_SPINUP=on,CLM_FORCE_COLDSTART=on
./xmlchange CCSM_CO2_PPMV=280.
./xmlchange STOP_N=35,STOP_OPTION=nyears
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=5
./xmlchange JOB_PRIORITY=regular
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=12:00:00

2. SASU

Diagnostics:
Diags (if available)

https://webext.cgd.ucar.edu/I1850/$CASE/lnd/

---

Output:
Output (if still available):
/glade/scratch/$USER/archive/$CASE/

---

Contacts:
@wwieder

---

Extra details:

[wwieder]

Quick look at the timing files suggest this

• 16 column setup gets 73 years / day; while the
• 9 column setup gets 103 ypd.

[wwieder]

159_test

Evaluate fire fluxes with hillslope hydrology turned on
/glade/work/wwieder/wiemip/wetlands_fire.n01.ctsm5.4.029_f09noCrop_HH_159_test

This case had to use the out of the box surface dataset to match initial conditions from 147.

./create_clone --clone /glade/work/wwieder/wiemip/wetlands_fire.n01.ctsm5.4.029_f09noCrop_HH_159_AD \
   --case /glade/work/wwieder/wiemip/wetlands_fire.n01.ctsm5.4.029_f09noCrop_HH_159_test \
   --keepexe

finidat = '/glade/derecho/scratch/wwieder/archive/ctsm5.4.022_BGC_hillslope_f09_147_pSASU/rest/0031-01-01-00000/ctsm5.4.022_BGC_hillslope_f09_147_pSASU.clm2.r.0031-01-01-00000.nc'

./xmlchange CLM_ACCELERATED_SPINUP=off,CLM_FORCE_COLDSTART=off
./xmlchange CCSM_CO2_PPMV=280.
./xmlchange JOB_PRIORITY=premium,STOP_N=10
./xmlchange RUN_TYPE=startup
./case.setup
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=06:00:00

Compare this with results from #147, diagnostics posted here, which show:

• higher burned area (maybe too high), but also
• very high CH4 fluxes (> 700 Tg CH4/y).

---

---

test2

Sanity check with initial conditions from results after 100 years of AD spinup & pointing to WIEMIP surface data.

./create_clone --clone /glade/work/wwieder/wiemip/wetlands_fire.n01.ctsm5.4.029_f09noCrop_HH_159_test \
   --case /glade/work/wwieder/wiemip/wetlands_fire.n01.ctsm5.4.029_f09noCrop_HH_159_test2 \
   --keepexe

Here pointing to initial conditions generated after 100 years of 159_AD spinup

test3

As above with f_ch4 = 0.04 to reduce apparently high methane fluxes from HH configuration.

./create_clone --clone /glade/work/wwieder/wiemip/wetlands_fire.n01.ctsm5.4.029_f09noCrop_HH_159_test2 \
  --case /glade/work/wwieder/wiemip/wetlands_fire.n01.ctsm5.4.029_f09noCrop_HH_159_test3 \
  --keepexe

[wwieder]

Quick test of results suggest that fire results look better in the hillslope model with the changes suggested for the wetland.n02 tag. Annual burned area is 11.4 M km2, diagnostics posted here

This also reduces LAI, as expected

[wwieder]

Compared to #131, an older f09 flat case the hillslope results in a larger annual burned are are higher 11.4 vs. 9.0 M km2, see diags posted here.

Despite higher burned area, especially in savannah regions, initially it seems LAI may be higher in the hillslope case than a similar 'flat' simulation. Will check if this holds when using initial conditions that reflect the configuration of this case.