wetlands_fire.n01.ctsm5.4.029_f09noCrop_9HH_160_ [AD, SASU, pSASU, HIST]

[wwieder]

Description:
Clone of #159 that uses 9 column hillslope file.

• f09, Hillslope, noCrop w/
• WIEMIP surface data and DATM
• Code mods 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_9HH_160_AD

---

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

 ./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_9HH_160_AD \
   --keepexe

---

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.

! 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=50,STOP_OPTION=nyears
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=4
./xmlchange JOB_PRIORITY=regular
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=12:00:00

---

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]

test2

Sanity check to look at initial results from the 9 column 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_9HH_160_test2 \
   --keepexe

10 year test pointing to year 150 of 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_9HH_160_test2 \
  --case /glade/work/wwieder/wiemip/wetlands_fire.n01.ctsm5.4.029_f09noCrop_9HH_160_test3 \
  --keepexe

[wwieder]

From tests listed above annual burned area is still to low (5.8 M km2), but we can tune annual CH4 fluxes to ~180 Tg CH4.

@linniahawkins you may be interested in looking at LAI from these hillslope runs? Looking higher compared to #131, especially in boreal forests. Below are results from Eurasia, but many regions have:

• Wetter soils;
• Higher LAI and ET fluxes; with
• Lower burned area and runoff.

See full diagnostics here

[linniahawkins]

Thanks - seems pretty consistent. I'm a little surprised how different the TSA between the two runs (1 degree in the amazon). Is WIEMIP using CRUJRA or ERA5?

[wwieder]

nice spot, @linniahawkins. I'm assuming you're talking about this?

WIEMIP is using a version of CRUJRA from TRENDY2025, whereas we're using the version from 2024. The datasets are more different than I'd expect. That said, in #158 I looked differences in forcing data only and didn't notice anything for TBOT over this region. I wonder if it has something to do with the years being used in each of these 10 year simulations (not the full 20 years we typically spin up over)?