Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

b.e30_beta04.BLT1850.ne30_t232_wgx3.116 #19

Open
cecilehannay opened this issue Nov 12, 2024 · 18 comments
Open

b.e30_beta04.BLT1850.ne30_t232_wgx3.116 #19

cecilehannay opened this issue Nov 12, 2024 · 18 comments

Comments

@cecilehannay
Copy link
Collaborator

Purpose:
Add the new dust + latest dust tuning, adjust sea-salt and reduce c8


Description:
Same as #17 except:

  • add to user_nl_cam
seasalt_emis_scale = 1.0D0
clubb_c8=  4.4
dust_emis_method = 'Leung_2023'
dust_emis_fact = 2.3

Case directory:


Sandbox:


Diagnostics:


Output:

  • Raw output on derecho (if still available): /glade/derecho/scratch/hannay/archive/b.e30_beta04.BLT1850.ne30_t232_wgx3.116/
  • CAM Climos: /glade/campaign/cgd/amp/amwg/climo/b.e30_beta04.BLT1850.ne30_t232_wgx3.116

Contacts:
@cecilehannay @dlawrenncar, @JulioTBacmeister, @adamrher, @PeterHjortLauritzen, @wwieder


Extra details:

@wwieder
Copy link

wwieder commented Nov 12, 2024

I think we want this run to use the new land initial conditions files that are consistent with the updated surface dataset and parameter file that were developed for the beta04 tag. @gustavo-marques can this run be stopped with the following namelist changes for CLM, or are you wanting this simulation to be consistent with #18 and #17, @cecilehannay?

finidat = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/initdata_esmf/ctsm5.3/ctsm53n04ctsm52028_ne30pg3t232_BgcCrop_exice_pSASU.clm60.r.0121-01-01.nc'

reseed_dead_plants = .true.

! turn off shifting cultivation
do_grossunrep = .false.

! do this for isotopes
use_c13 = .true.
use_c14 = .true.
use_c13_timeseries = .true.
use_c14_bombspike = .true.

@NCAR NCAR deleted a comment from PeterHjortLauritzen Nov 12, 2024
@cecilehannay
Copy link
Collaborator Author

@wwieder: I can do another run if I you want (if you think it is worth it).

@wwieder
Copy link

wwieder commented Nov 12, 2024

At some point it think it would be helpful to use land initial conditions that are reflective of the parameters and surface dataset we're running with, but maybe it's not that critical since the climate the land is seeing is likely also different in these coupled runs.

@cecilehannay
Copy link
Collaborator Author

@adamrher: This run is still short but from a quick look, it looks like RESTOM is going to be out of balance. Should we bring c8 to where it was in 114? I am not even sure it will be enough. Thoughts?

@adamrher
Copy link
Collaborator

I expected RESTOM would be not be in balance -- in particular I am guessing it's positive. That's OK, let's just keep running at least 20 years so we can understand where we are at.

@dlawrenncar
Copy link
Collaborator

dlawrenncar commented Nov 14, 2024 via email

@adamrher
Copy link
Collaborator

The RESTOM in #116 starts at +2 W/m2, and quickly descends down to +1 W/m2. So it's definitely problematic. However I'd suggest that we continue this run for 60 years or so because we can test the hypothesis that the cool SSTs in #112 degrades ENSO variability

I prefer to evaluate biases in PD, and so I'll plan to have for you Tuesday some bias plots of the radiative fluxes based on an F-case equivalent of 116. Are you separately able to compare the current land albedo to obs?

I'm also going to do another F-case like 116, but cranking up the seasalt_emis_fact=1.5. This degrades the AODVIS skill, but what I've discovered today is our clear-sky net shortwave at the TOA is about +3 W/m2 too high compared to ceres-ed4.1, even though the AOD is "correct". I'd argue it's more important to have the correct radiative fluxes rather than the correct AOD, as how the model uses a given AOD to compute the radiative fluxes is highly uncertain.

@dlawrenncar
Copy link
Collaborator

OK. @wwieder Some evaluation of land albedo would be good for the next CESM Project meeting, perhaps looking at run 116 and also comparing it to land-only runs to get a sense of how the albedo bias is affected by the coupled climate.

@olyson
Copy link

olyson commented Nov 15, 2024

Land diagnostics are here

@wwieder
Copy link

wwieder commented Nov 15, 2024

You're amazing. Thanks @olyson

@cecilehannay
Copy link
Collaborator Author

cecilehannay commented Nov 15, 2024

@olyson or @slevis-lmwg: would you mind adding the land diags for

/glade/derecho/scratch/hannay/archive/b.e30_beta04.BLT1850.ne30_t232_wgx3.116
versus
/glade/derecho/scratch/hannay/archive/f.e30_alpha04a.FMTHIST.ne30_L93_tuning.008

We noticed large differences in AODDUST between our B cases and our F cases. That said the setting is very different between B and F cases:

  • B case: 1850, low top and CLM%BGC
  • F case: 2000-2010, high top top and CLM%SP
    TLAI, FV, and H2OSOI etc will all look different, so differences in AODDUST is expected.

However, Danny would like to look at the land diagnostics to make sure everything looks ok. This is the reason we are requesting this comparison. Thanks a bunch.

@dmleung @adamrher @tilmes

@slevis-lmwg
Copy link

I will try to address this request over the weekend.

@olyson
Copy link

olyson commented Nov 16, 2024

Land diagnostics, compared to a 1850 pSASU ctsm5.3 run (land-only) are here.
Either @slevis-lmwg or I will address the other diagnostics request (comparison to F-case).

@adamrher
Copy link
Collaborator

adamrher commented Nov 18, 2024

When we get a sense of what RESTOM is, could we see a map and/or a
longitudinal average. The land group would like to think about whether or
not albedos could be adjust to help, but it would be good to know if this
seems like a tropical problem or a high latitude problem, if it can be
established.

Here are some plots of the net clear sky SW biases at the surface and the top of the model. I added the latest coupled run, but the biases in FMTHIST (011, 012) 2000-2009 runs are more appropriate for comparison to the 2000-2009 ceres data shown. Note the FMT compsets run with SP mode whereas BLT runs with BGC.

Clear-sky net SW bias at surface
Clear-sky net SW bias at TOA

The NE Canada biases of too much absorption seem the most prominent to me in these plots. The biases of clear-sky net sw at the surface for S.America are large, but they are more attenuated at the TOA. Let me know if you want to see plots of the other terms in the radiation budget.

@olyson
Copy link

olyson commented Nov 19, 2024

Diagnostics of b.e30_beta04.BLT1850.ne30_t232_wgx3.116 compared to f.e30_alpha04a.FMTHIST.ne30_L93_tuning.008 are here

@wwieder
Copy link

wwieder commented Nov 27, 2024

I'm not sure where to put this comment, but I wonder if land initial conditions may be somewhat to blame for regional biases over the Amazon that seem to have low precip , low clouds, low humidity, and high incoming shortwave we're seeing in this case. We also know that leaf area index is also very low, but it seems like this may be related to the land initial conditions?
See also discussion on LMWG_dev

Below is the difference in LAI in the first month of the simulation between a new F2000-control case that Sam L is running (that uses new CTSM5.3 initial conditions) from this B-116 case (that uses land initial conditions from previous CESM3 test simulations).

image

If we initialize with a healthier (leafier) forest can we increase ET, cloud and precipitation and reduce incoming solar radiation in ways that even help RESTOM?

@wwieder
Copy link

wwieder commented Nov 27, 2024

Note, an 1850 initial condition files is available here
/glade/campaign/cgd/tss/people/oleson/CLM5_restarts/ctsm53n04ctsm52028_ne30pg3t232_SASU.clm2.r.0281-01-01-00000.nc

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

7 participants