Skip to content

Conversation

@mingxuanwupnnl
Copy link
Contributor

@mingxuanwupnnl mingxuanwupnnl commented Oct 14, 2025

Add tunable aerosol diameter thresholds for renaming to strat coarse mode in MAM5

When using CMIP7 volcanic SO2 emission, the model will produce larger stratospheric
AOD and sulfate burden due to increased emission amount as well as changing injection
height. The stratospheric sulfate burden during Pinatubo eruption when using CMIP7
volcanic SO2 emission is larger than HIRS observations.

For v3.HR simulations, we identified two tunable aerosol mode diameter thresholds
for renaming accumulation mode (mode 1) sulfate to stratospheric coarse mode (mode 5).
By changing them from 0.4 um to 1 um, changing the dgnumhi in mode 1 physprop file,
and changing dgnumlo in mode 5 physprop file, the stratospheric AOD and sulfate
burden during Pinatubo are reduced, the RESTOM and SW flux at TOA are also reduced
to the level when using default CMIP6 volcanic emission.

In this PR, we only introduce two namelist variables for the two tunable thresholds
and keep values the same as v3.LR simulations.

[BFB]
[NML] Two new namelist variables introduced affecting all tests using EAMv3


Testing results.

v3 LR F20TR Pinatubo_CMIP6-CMIP7

@mingxuanwupnnl mingxuanwupnnl added the EAM Fortran-based E3SM Atmosphere Model label Oct 14, 2025
@mingxuanwupnnl mingxuanwupnnl changed the title Add tunable aerosol diameter thresholds for renaming to strat coarse mode in MAM5 [WIP] Add tunable aerosol diameter thresholds for renaming to strat coarse mode in MAM5 Oct 14, 2025

<entry id="dp_xferall_thresh_accum_rename" type="real" category="scavenging"
group="phys_ctl_nl" valid_values="" >
accum mode aerosol diameter threshold for renaming to strat coarse mode
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@mingxuanwupnnl . Can you make the descriptions for these two namelist variables distinguishable?

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@mingxuanwupnnl @wlin7
dp_cut
Cutoff (geometric-mean) diameter for a lognormal mode. When growth pushes a mode’s mean size above dp_cut, the excess tail is shifted to the next larger mode. Applied each timestep to keep each mode’s mean/width physically reasonable.

dp_xferall_thresh
“Full-transfer” trigger. If the updated mean diameter of a mode exceeds this threshold, the entire mode (mass/number) is transferred to the next larger mode. This is a safety valve for runaway growth that dp_cut could not contain over prior steps.
Example: If the new mean diameter is 0.8 µm and dp_xferall_thresh = 0.6 µm, the mode is fully transferred—preventing an unphysical, oversized mean in that mode.

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@mingxuanwupnnl Based on Ziming's explanation, I suggest using the following brief description for the two namelist variables:
dp_cut_accum_rename: Cutoff diameter for larger accum mode particles renamed to the strat coarse mode
dp_xferall_thresh_accum_rename: Threshold diameter for all accum mode particles transferred to the strat coarse mode

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@wlin7 have revised the description according to @keziming 's explanation in three places

Copy link
Contributor

@keziming keziming left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

looks good!

@rljacob
Copy link
Member

rljacob commented Oct 15, 2025

Remove the [WIP] from the title. If it still is in progress, convert it to a draft PR. (link is under Reviewers).

@mingxuanwupnnl mingxuanwupnnl changed the title [WIP] Add tunable aerosol diameter thresholds for renaming to strat coarse mode in MAM5 Add tunable aerosol diameter thresholds for renaming to strat coarse mode in MAM5 Oct 15, 2025
@mingxuanwupnnl
Copy link
Contributor Author

Remove the [WIP] from the title. If it still is in progress, convert it to a draft PR. (link is under Reviewers).

removed.

add_default($nl, 'mam_amicphys_optaa'); #for enabling amicphys code in eam
add_default($nl, 'n_so4_monolayers_pcage'); #number of so4(+nh4) monolayers needed to "age" a carbon particle
add_default($nl, 'dp_cut_accum_rename');
add_default($nl, 'dp_xferall_thresh_accum_rename');
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@mingxuanwupnnl please add a brief description of the two namelist variables here.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

added


<entry id="dp_xferall_thresh_accum_rename" type="real" category="scavenging"
group="phys_ctl_nl" valid_values="" >
accum mode aerosol diameter threshold for renaming to strat coarse mode
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@mingxuanwupnnl Based on Ziming's explanation, I suggest using the following brief description for the two namelist variables:
dp_cut_accum_rename: Cutoff diameter for larger accum mode particles renamed to the strat coarse mode
dp_xferall_thresh_accum_rename: Threshold diameter for all accum mode particles transferred to the strat coarse mode

@mingxuanwupnnl mingxuanwupnnl added BFB PR leaves answers BFB Atmosphere labels Oct 16, 2025
Copy link
Contributor

@wlin7 wlin7 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The changes look good and consistent. e3sm_atm_integration tests for default configuration passed BFB against master baselines. The new tuning also tested in a new coupled configuration.

@wlin7 wlin7 added the NML label Oct 21, 2025
wlin7 added a commit that referenced this pull request Oct 21, 2025
…7795)

Add tunable aerosol diameter thresholds for renaming to strat coarse mode in MAM5

When using CMIP7 volcanic SO2 emission, the model will produce larger stratospheric
AOD and sulfate burden due to increased emission amount as well as changing injection
height. The stratospheric sulfate burden during Pinatubo eruption when using CMIP7
volcanic SO2 emission is larger than HIRS observations.

For v3.HR simulations, we identified two tunable aerosol mode diameter thresholds
for renaming accumulation mode (mode 1) sulfate to stratospheric coarse mode (mode 5).
By changing them from 0.4 um to 1 um, changing the dgnumhi in mode 1 physprop file,
and changing dgnumlo in mode 5 physprop file, the stratospheric AOD and sulfate
burden during Pinatubo are reduced, the RESTOM and SW flux at TOA are also reduced
to the level when using default CMIP6 volcanic emission.

In this PR, we only introduce two namelist variables for the two tunable thresholds
and keep values the same as v3.LR simulations.

[BFB]
[NML] Two new namelist variables introduced affecting all tests using EAMv3
@wlin7
Copy link
Contributor

wlin7 commented Oct 21, 2025

Merged to next.

@wlin7 wlin7 merged commit ebe5d0c into master Oct 22, 2025
8 checks passed
@wlin7 wlin7 deleted the mingxuanwupnnl/atm/MAM5_tuning_forv3HR branch October 22, 2025 15:17
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

Atmosphere BFB PR leaves answers BFB EAM Fortran-based E3SM Atmosphere Model NML

Projects

None yet

Development

Successfully merging this pull request may close these issues.

6 participants