Things needed in CTSM6 for Emissions driven simulations

[ekluzek]

For cesm3.0 we need to have emissions driven simulations run. This requires for the emission driven capabilities to be: in place, tested and verified to work correctly.

Emissions driven simulations have been run based off of the cesm2.1.5 release. Here's a repository for that:

https://github.com/lawrencepj1/CESM2.1.5-CCIS-Ens

for the tag: baselines_tag_06

The version of CTSM used is based off of release-clm5.0.37 and is tagged as active-bvocs_tag_01
The fork for CTSM is: https://github.com/lawrencepj1/CESM2.1.5-CLM5-New-Land-Use

There are some small namelist differences, bringing in a FUN fix (that's on CTSM master) and some changes to fire-emissions, that look like they need to come to master.

A git diff shows:

git diff --stat release-clm5.0.0 
 bld/CLMBuildNamelist.pm                       |  1 -
 bld/namelist_files/use_cases/1850_control.xml |  2 +-
 src/biogeochem/CNFUNMod.F90                   |  7 ++++---
 src/biogeochem/CNFireEmissionsMod.F90         | 27 +++++++++++++++++++--------
 src/biogeochem/CNPhenologyMod.F90             | 46 +++++++++++++++++++++++++++++++++++-----------
 5 files changed, 59 insertions(+), 24 deletions(-)

For the fully coupled case, the compset used is: BHIST_BPRPcmip6, which has the longname:

HIST_CAM60_CLM50%BGC-CROP-CMIP6DECK_CICE%CMIP6_POP2%ECO%ABIO-DIC_MOSART_CISM2%NOEVOLVE_WW3_BGC%BPRP

(So the CMIP6 user-mods are used to get the standard CMIP6 output)

Changes in user_nl_clm, are mostly changes to history output, but also turning on carbon isotopes. As well as this important change for irrigation:

! Increase Irrigation Target SMP to better match other crops besides wheat

irrig_target_smp = -1500.

Anb the following settings for fire-emissions and megan:

! Fire emission file and factor settings

fire_emis_factors_file = '/glade/campaign/cgd/tss/lawrence/ccisensinput/fire_emission_factors/fire_emission_factors_78PFTs_ACTIVEFIRE_c20240529.nc'
fire_emis_specifier = 'bc_a4 = BC', 'pom_a4 = 1.4 * OC', 'num_a4 = num_a4_bc + num_a4_pom', 'SO2 = SO2', 'DMS = DMS',
       'SOAG = 0.45 * BIGALK + 0.35 * BIGENE + 1.725 * TOLUENE + 1.4625 * BENZENE + 1.9875 * XYLENE + 0.34 * ISOP + 4.25 * C10H16'
fire_emis_elevated = .false.

! BVOC emission file and factor settings

megan_factors_file = '/glade/campaign/cgd/tss/lawrence/ccisensinput/megan_emission_factors/megan21_emis_factors_78pft_ACTIVEBVOC_c20240529.nc'
megan_mapped_emisfctrs = .false.
megan_specifier = 'SOAG = 0.25 * isoprene + 2.15 * (carene_3 + pinene_a + thujene_a + bornene +',
      'terpineol_4 + terpineol_a + terpinyl_ACT_a + myrtenal + sabinene + pinene_b + camphene +',
      'fenchene_a + limonene + phellandrene_a + terpinene_a + terpinene_g + terpinolene +',
      'phellandrene_b + linalool + ionone_b + geranyl_acetone + neryl_acetone + jasmone +',
      'verbenene + ipsenol + myrcene + ocimene_t_b + ocimene_al + ocimene_c_b + 2met_nonatriene) + ',
      '5.25 * (farnescene_a + caryophyllene_b + acoradiene + aromadendrene + bergamotene_a +',
      'bergamotene_b + bisabolene_a + bisabolene_b + bourbonene_b + cadinene_d + cadinene_g +',
      'cedrene_a + copaene_a + cubebene_a + cubebene_b + elemene_b + farnescene_b +',
      'germacrene_B + germacrene_D + gurjunene_b + humulene_a + humulene_g + isolongifolene +',
      'longifolene + longipinene + muurolene_a + muurolene_g + selinene_b + selinene_d +',
      'nerolidol_c + nerolidol_t)'

XML settings that are important are:

	CLM_CO2_TYPE: prognostic
	CCSM_BGC: CO2C

[ekluzek]

@lawrencepj1 is the irrigation change something that we'll want in the cesm3_0 release for ctsm6?

[ekluzek]

There are also changes in: POP, CICE, and cime. There are also a bunch of settings for CAM and new files that will need to be handled. For these purposes I don't think the POP or CICE changes matter. The cime changes have to do with the expressions for MEGAN and have already come into the cesm3_0 code base.

So I think the things that have to be handled for CTSM are:

• MEGAN namelist settings
• fire-emission namelist settings
• fire-emission code changes
• Running simulations to tune fire-emissions
• Bring in the new datasets from the tuning work

[tilmes]

I think, that we need to change SOAG -> SOAE and the factors are different from what we are using: 'SOAE = 0.5954isoprene + 5.1004(carene_3 + pinene_a + thujene_a + bornene +',
' terpineol_4 + terpineol_a + terpinyl_ACT_a + myrtenal + sabinene + pinene_b + camphene +',
' fenchene_a + limonene + phellandrene_a + terpinene_a + terpinene_g + terpinolene +',
' phellandrene_b + linalool + ionone_b + geranyl_acetone + neryl_acetone + jasmone +',
' verbenene + ipsenol + myrcene + ocimene_t_b + ocimene_al + ocimene_c_b + 2met_nonatriene) + ',
' 12.3942*(farnescene_a + caryophyllene_b + acoradiene + aromadendrene + bergamotene_a +',
' bergamotene_b + bisabolene_a + bisabolene_b + bourbonene_b + cadinene_d + cadinene_g +',
' cedrene_a + copaene_a + cubebene_a + cubebene_b + elemene_b + farnescene_b +',
' germacrene_B + germacrene_D + gurjunene_b + humulene_a + humulene_g + isolongifolene +',
' longifolene + longipinene + muurolene_a + muurolene_g + selinene_b + selinene_d +',
' nerolidol_c + nerolidol_t)'

[tilmes]

I am also not sure what changed in the new emissions file compared to the file we are using: /glade/campaign/cesm/cesmdata/inputdata/atm/cam/chem/trop_mozart/emis/megan21_emis_factors_78pft_c20161108.nc

@Emmons should we use that new emis factor file for our standard runs?

[ekluzek]

@tilmes thanks for chiming in here. Last week, @lawrencepj1 confirmed that the files he's using for CESM2.1.5 will need to be tuned based on running in cesm3_0 and using CMIP7 emissions. So the CESM2.1.5 HOPE simulation files won't be moved over to cesm3.0 without that work done first. I updated the top section to reflect that. So I think that means you probably do NOT want to use those files in your latest cesm3_0 work.

[lkemmons]

@lawrencepj1 Could you please describe how you decided to "tune" the megan emissions?

[lawrencepj1]

Hi Louisa, Erik and Simone,

Sorry for the slow response here. It has been a very busy week. Here is the basis for the tuning of the emission factors in both the fire model and the MEGAN model. I have attached a word doc with figures and the final tuned values.

Tuning of CESM 2.1.5 with MAM4 for Active Fire and BVOC

1. Initial Fire Emission and BVOC emission factors were provided by Simone Tilmes from simulations performed by ACOM.
2. Aerosol burdens were then tuned to match the LENS2 CMIP6 specified emissions maintaining the relative relationships between constituent emission factors.
3. New fire_emis_factors_file has been scaled by 0.68 across all components "Comp_EF" variables to reduce total fire aerosol load.
4. New megan emission factors file has been scaled by 0.5 across all components "betaT" variable values to reduce the temperature response of non isoprene emissions.
5. New megan_specifier values multiplied but 0.42 and mapped to SOAG rather than SOAE to account for the difference in pathway for SOAE than SOAG in MAM4.
6. No code changes were made to either the fire or megan models.

Happy to talk more
Best
Peter

firemegantuningemissionfactors.docx