ACCESS-ESM1.6 scope and use for CMIP7 Fast Track

ACCESS-ESM1.6 scope

ACCESS-ESM1.6 needs to be defined and tested during the remainder of 2024 so that spin-up can begin in 2025. Any upgrades from ACCESS-ESM1.5 will need to be relatively simple to implement and test. Please use this topic to

  • suggest model improvements
  • report any bugs that should be fixed
  • note any additional diagnostics that would be useful (or were incorrect for CMIP6).

Already in scope are

  • updates to WOMBAT, led by @pearseb
  • updates to the simple land-use scheme to better account for wood harvest, led by @tammasloughran
  • potential use of ACCESS-CM2 CABLE settings and use of the CABLE3 codebase instead of CABLE2.4

CMIP7 Fast Track experiments
Please also use this topic to indicate if you are interested in helping to run and analyse any of the Fast Track experiments. Listing the experiments that are of most interest to you would be really helpful. More information about the experiments is available
here - AR7 Fast Track experiment viewer.

Hello Rachel, all,

Thank you for leading this. It is great to see improvements to CABLE and WOMBAT! I understand the Fast track timeline is challenging and the ultimate goal is the development of the ACCESS-ESM3. Since there is more focus on the emission scenarios and bgc for CMIP7 and fast track, maybe improvements on the bgc only are sufficient.
We are planning to participate in PMIP with the LIGabrupt, but without much excitement. The ACCESS-ESM1.5 performances for the LIG were relatively good, with maybe too much sea-ice in the Arctic and very low sea-ice in the SO. The UK group suggests that the inclusion of melt-ponds in the sea-ice model is crucial to simulate changes in the Arctic sea-ice.
In our other experiments (i.e. AMOC weakening, 49ka exp), we have noted a very slow (~600 yrs) response of the SO sea-ice and little sensitivity (a slight sea-ice increase instead of decrease for an AMOC weakening and relatively small sea-ice increase at 49ka).
This leads me to wonder whether there would be merit in replacing CICE4.1 with CICE5.1.2 in the ACCESS-ESM1.6. That would mean that the ACCESS-ESM1.6 would then have the same structure as the ACCESS-OM2 series… just floating an idea really, realizing that the coupling to the UM/CABLE and re-tuning involved might not make it straightforward.

In terms of helping with the fast track in 2025, we can do the LIG abrupt, could consider the abrupt0p5CO2 and/or other experiments of the scenarioMIP/C4MIP that could help.

1 Like

Use the inland basin and lake evaporation corrections to river runoff from CM2?

Is this the same setup as CM2 @MartinDix? In which case would it be pretty straightforward to use in an ESM1.6, or were the UM versions different enough that the modifications to CICE5 wouldn’t be compatible?

@RachelLaw it wasn’t named in your request, but the build and deployment infrastructure should also be included in this plan. And that relies a bit on what code upgrades are required. The closer the upgrades are to existing models (like OM2 and CM2) the easier it will be.

How different code is incorporated is also important, particularly anything that is shared between multiple models but is a sub-component. I’m thinking particularly of WOMBAT.

@Aidan - another aspect is making sure we can run all the experiments under payu rather than the script-based runs that some of us are still using for ACCESS-ESM1.5.

@sofarrell will reply about CICE - basically the UM version means we wouldn’t get the advantages of newer CICE.

@Aidan, a few of us (@RachelLaw, @tiloz, @matthew.chamberlain, @pearseb, @AndyHoggANU) met to chat about this last week. The gist is this:

  • The WOMBAT code as it is currently used in ACCESS-OM2 and ACCESS-ESM1.5 is written directly into MOM5 src as a tracer package. Let’s call this the MOM5-version of WOMBAT.
  • ACCESS-NRI has rewritten WOMBAT so that it can be used by both MOM5 and MOM6. This is being tested with MOM6 in ACCESS-OM3. Let’s call this the generic-version of WOMBAT.
  • At the same time, @pearseb has been doing his science developments into the MOM5-version of WOMBAT.
  • The intention is to use the generic-version of WOMBAT in ACCESS-ESM1.6, but some additional work is still needed to couple it. In the meeting, it was decided that ACCESS-NRI would prioritize this work. At the same time, @pearseb will continue his developments in the MOM5-version (I think the code is mostly frozen now and he’s testing?). The hope is that both efforts will finish in time to port @pearseb’s developments into the generic-version of WOMBAT and use that in ACCESS-ESM1.6. This will mean changes to the way things are built relative to current ACCESS-ESM1.5, so all this will need to be done in consultation with you and your team.

ADDED: to be clear (or not) my plan is also to update ACCESS-ESM1.5 and ACCESS-OM2 to use the generic-version of WOMBAT, so there will be new releases of those.

1 Like

Very much yes.

HI @LaurieM, @RachelLaw,

It’s a nice idea to upgrade our models to align with ACCESS-OM2, but, it’s the coupling to the UM7.3 that’s the issue. We can couple CICE5.1.2 to UM7.3 but we would only get the old zero layer thermodynamics we had in CICE4.1. as the coupling for multilayer ice came in later in GA7 Met Office version.

This zero layer ice scheme has just been removed entirely from the CICE6.5 code after being depreciated in the last few versions that’s how old and out of date it is.

In terms of melt pond schemes there is a simple allowance for melt ponds in the reduction of albedo when the temperature is close to melting in ACCESS1.5 (its in the atmospheric part of the UM7.3 code). In CICE5.1.2 there are two schemes, we use the UK one from CPOM in Reading there is an NCAR one as well.

The complication there is the melt pond information needs to be then passed back through the coupler back to the UM so the atmospheric model can do the albedo calculation (it duplicates code that sits in the ice model) and then passed to the UM radiation code, that set up is not available to us in UM7.3 it came in GA6 not sure which version of UM that was Martin (it was 2015 timeframe), but in our code base that we still maintain with that in was UM10.6, so it’s not going to work for what @LaurieM needs in the paleo runs.

So Unfortunately the answer is no there is nor advantage in upgrading to CICE5.1.2.

Siobhan

2 Likes

HI @RachelLaw The bug fix that we should fix that also impacts on diagnostics is the OASIS MCT set up which was a test one that went into ACCESS1-4, @dhb599 updated it and corrected the bug for ACCESS-CM2 but that sits inside the CICE5.1.2 model, as OASIS_MCT links across to the ice-ocean combined.
The advantage of this is we can then do the time derivative ice statistics, there may be some for the land and ocean and BGC time derivative statistics that were impacted in ESM1.5 by this coupler bug as well that we could use if its fixed.

Also we need to set up the “CMIP6” ice statistics internally in the CICE4.1 code, which will improve the quality of diagnostics over what we could obtain through the post-processing APP4 output. We will use CICE6.4/5 as our framework for this rather than the CM2 code as there were some bugs that came through from a “Beta” GSI8.1 of the CMIP6 output set up.

@sofarrell what is this bug? Is it and the fix documented somewhere?

Hi @dougiesquire I think we may have made a one line mention of it in Tilo’s paper on the model. It basically took 2 atmospheric time steps before it did a coupling time step to the ocean right at the start, Thats why it messes up the d/dt type statistics. (In Jan and July as we did 6 month chunks). No it hasn’t been documented. Even worse the comments in the code match the earlier ACCESS1-3 coupling set up. So don’t even go and look at you will get very very confused.

Dave totally re-wrote that section of code for ACCESS-CM2 and thats what we need to replace into this code for ESM1.6. We were too hard pressed doing two model versions back in “CMIP6” timeframe, to put into ACESSS ESM1.5 and only realised about the diagnostics too late.

It’s not a simple translation CM2 has code for CICE5.1.2 (our version is actually GSI8.1 from the Met office with local #ifdef changes) and ESM1.5 (ESM1.6) is CICE4.1 (The coupling code passes through CICE across for all ice/ocean grid points; unlike what is being done with the mediator and NUOPC).

1 Like

Many of the land improvements since ESM1.5 are automatically available as part of the CABLE3 codebase (notably for water and energy balance improvements). Land centric improvements that were developed as part of CABLE-CM2 that should (ideally) be back-ported (and would need code changes):

  • revisions to the meteorological forcing used in the ‘implicit’ section
  • lake water/inland water basins/river runoff fix (if technically possible)
  • CABLE specific modifications to MOSES routine im_sf_pt2 that addresses the energy balance of coastal grid cells.
  • Changes to the CABLE-MOSES interface code that facilitate a broader range of diagnostic outputs (e.g. sublimation, soil evaporation and transpiration partitioning) as requested by the community/CMIP.

Discussions around adopting the CM2 Antarctic/Greenland run off and iceberg scheme in ESM1.6 have also commenced (requiring code changes in CABLE, MOM5, CICE and the OASIS).

We should also review the sequence of the calling point to CABLE in the explicit and implicit sections (they weren’t optimal in CM2) and the MOSES science options invoked by CABLE (some of these are not compatible with CABLE science).

Land science that could be assessed as part of a ESM1.6 calibration exercise include:

  • updates to PFT parameters
  • use of litter layer and/or Or soil evaporation model
  • use of Medlyn stomatal function
  • use PFT-dependent root distributions
  • use of the updated (pitot) functions which link thermal and hydraulic parameter values to the soil texture and composition.
  • use of revised soil and other land ancillaries