Difference between revisions of "ISMIP6-Projections-Greenland"

From CliC Wiki
Jump to: navigation, search
m (Initial state and experiment duration)
(Greenland Standalone Ice Sheet Modeling)
 
(120 intermediate revisions by 4 users not shown)
Line 1: Line 1:
 
== Overview ==
 
== Overview ==
'''NOTE: THIS PAGE IS UNDER CONSTRUCTION. WHEN IT IS FINISHED, THIS STATEMENT WILL BE REMOVED'''
+
'''NOTE: THIS PAGE IS UNDER CONSTRUCTION. WHEN IT IS FINISHED, THIS STATEMENT WILL BE REMOVED. LAST UPDATE: 19 September  2019'''
  
 
This page describes the experimental protocol for the ISMIP6 projections that target the upcoming IPCC AR6 assessment. Due to the delay in CMIP6 climate simulations, the initial set of ISMIP6 simulations are based on CMIP5 projections. As CMIP6 model output become available, ISMIP6 will include simulations based on these models.  
 
This page describes the experimental protocol for the ISMIP6 projections that target the upcoming IPCC AR6 assessment. Due to the delay in CMIP6 climate simulations, the initial set of ISMIP6 simulations are based on CMIP5 projections. As CMIP6 model output become available, ISMIP6 will include simulations based on these models.  
Line 8: Line 8:
 
• Sampling CMIP scenarios: main focus is on the high emission RCP8.5, but ice sheet evolution in response to low emission RCP2.6 is also investigated.   
 
• Sampling CMIP scenarios: main focus is on the high emission RCP8.5, but ice sheet evolution in response to low emission RCP2.6 is also investigated.   
  
• Sampling CMIP models: 6 AOGCMs have been selected from the CMIP5 model ensemble. The AOGCMs were identified based on the following steps: 1) present plausible climates near Greenland (evaluated by model biases over the historical period), 2) have the data at the temporal resolution needed for RCM downscaling, 3) sample a diversity of forcing (evaluated by differences in projections and code similarities).
+
• Sampling CMIP models: 6 AOGCMs have been selected from the CMIP5 model ensemble. The AOGCMs were identified based on the following steps: 1) present plausible climates near Greenland (evaluated by model biases over the historical period), 2) have the data at the temporal resolution needed for RCM downscaling, 3) sample a diversity of forcing (evaluated by differences in projections and code similarities) and 4) allow only models with RCP8.5 and ideally with RCP2.6. The sampling methodology for the CMIP5 models is described in Barthel et al. (under review). As CMIP6 models became available, ISMIP6 prepared dataset for CNRM-CM6 (ssp585 and ssp126), CESM2 (ssp585, in prep) and UKEM1-CM6 (ssp585). Unlike for the rigorous analysis for the CMIP5 models, the CMIP6 models were selected because of their availability.
  
• Sampling ice sheet model uncertainty: "standard" and "open" experiments. The "standard" experiments are based on parameterizations developed by the ocean and atmospheric focus groups, while "open" experiments utilize the parameterizations already in use by respective ice sheet models.
+
• Sampling ice sheet model uncertainty: "standard" and "open" experiments. The "standard" experiments are based on parameterizations developed by the ocean focus group, while "open" experiments utilize the parameterizations already in use by respective ice sheet models. The open experiments are an important contribution as the use of different parameterizations reflect the uncertainty in our current understanding of ice-ocean interactions.
  
• Sampling forcing uncertainty: the standard experiments include "high", "mid" and "low" parameters.  
+
• Sampling ocean forcing uncertainty: the standard experiments include "high", "mid" and "low" parameters.  
  
• Experiment ranking: This experimental framework results in 12 to 72 experiments, and a control run. Not every ice sheet model will be able to carry out the full set of experiments. The experiments are therefore ranked, groups are encouraged to work shown the list and complete as many experiments as possible. This approach is based on Shannon et al. (2013): it ensures that all groups do a subset of identical experiments, while it also allows faster models to explore the experiment space more fully.
+
• Experiment ranking: This experimental framework results in a series of projections (divided into core and targeted experiments), a historical run and a control run. Not every ice sheet model will be able to carry out the full set of experiments, but they are strongly encouraged to participate in the full suite of core experiments (Table 1, which uses both the "open" and "standard" experiments/parameterizations). Given that the standard experiments requires new implementations, it is OK to only participate with only the open experiments. Similarly, it is OK to only participate with the "standard" experiments.  Groups are encouraged to work through the lists presented below (see Table 1 for core experiments, and targeted experiments will follow soon), starting from the top, and complete as many experiments as possible. This approach of combining core and targeted experiments is based on Shannon et al. (2013): it ensures that all groups do a subset of identical experiments, while it also allows faster models to explore the targeted experiment space more fully.
  
 
[[File:Greenland_exp_design.png|720x405px]]
 
[[File:Greenland_exp_design.png|720x405px]]
  
 
Figure 1: Overview of the Greenland experimental framework
 
Figure 1: Overview of the Greenland experimental framework
 
  
== Experiment Ranking ==
+
== List of Projections ==
With the help of the atmosphere and ocean focus groups, six CMIP5 AOGCMs have been selected for ISMIP6 standalone ice sheet model projections. The table below lists the initial number of experiments based on the first three AOGCMs: '''<AOGCM1>''', NorESM, MIROC5. This table is the minimum contribution expected from ISMIP6 models and should be worked on in order. Modeling groups that can run many simulations will be encouraged to explore the ice sheet response using three additional CMIP5 AOGCMs '''[HG: Should we have another table in the appendix with a schematic extension for models 4-6?]'''. For groups that already have their own methods for implementing ocean and atmosphere forcing in place, the suite of "open" experiments (1-4) allows for a direct start. The basic requirement for the "open" experiments is that the forcing is based on the given CMIP5 GCM and uses the given duration for the  reference period and scenario (see below). All groups are expected to submit results for the "standard" experiments, using forcing methods suggested by ISMIP6. Depending on the results of experiments 3 and 7, which consider RCP2.6, additional AOGCMs may be suggested (exploring RCP2.6 in detail has low priority if the ISM response is very limited). Additional RCP2.6 experiments may be reserved for models that are able to do many simulations, but these would be of lower priority than completing the set with the 6 AOGCMs with RCP8.5.  
+
With the help of the atmosphere and ocean focus groups, a number of CMIP5 AOGCMs have been selected for ISMIP6 standalone ice sheet model projections. Table 1 lists the ''core experiments'' which are the minimum contribution expected from ISMIP6 models, in addition to the initialization experiments listed in Table 2. All groups are encouraged to contribute to the "standard" experiments. Groups that have their own methods for implementing ocean and atmosphere forcing, are encouraged to do the suite with "open" experiments (1-4), but these are not compulsory. Models that perform the "open" experiments can use the parameterization of their choice to simulate atmospheric and oceanic forcings, but these parameterizations must use the given CMIP5 AOGCM outputs.  
  
 +
Modeling groups that can run many simulations will be encouraged to further explore the ice sheet response using ''targeted experiments'' (See [https://docs.google.com/spreadsheets/d/1YEiPV3Uc0K8EqD57mz-ZqBIF4XnVb5UbyRDRmox0KeE/edit?usp=sharing Table]). These include three additional CMIP5 AOGCMs under RCP8.5, and experiments that explore the ocean forcing uncertainty. Depending on the results of experiments 3 and 7, which consider RCP2.6, additional AOGCMs may be suggested with RCP2.6 for models that are able to do many simulations, but these would be a lower priority than the completing the set of experiments with the 6 AOGCMs for the RCP8.5 scenario. As CMIP6 AOGCMs are becoming available, we are preparing these datasets. The spreadsheet will be updated as new dataset become available. Because there is value in both completing the 6 CMIP5 AOGCMs (to sample the uncertainty in CMIP) and simulations with CMIP6 models, we encourage groups do to as many experiments as possible.
 
    
 
    
''Note: As of Jan 21, MAR has completed MIROC5 and NorESM for both RCP8.5 and RCP2.6. Greenland model selection team is actively working on final selection''
+
''Note: As of August 6, all the datasets needed for the experiments in Table 1, as well as the CMIP5 and CMIP6 AOGCMs listed in the targeted experiments (See [https://docs.google.com/spreadsheets/d/1YEiPV3Uc0K8EqD57mz-ZqBIF4XnVb5UbyRDRmox0KeE/edit?usp=sharing Table]) are ready. ''
  
  
 
{| class="wikitable" style="margin: 1em auto 1em auto;"
 
{| class="wikitable" style="margin: 1em auto 1em auto;"
| align="center" style="background:#f0f0f0;"|'''Expt'''
+
! colspan="7"| Table 1: Core Experiments based on  MIROC5, NorESM1-M and HadGEM2-ES
 +
 
 +
|-
 +
| align="center" style="background:#f0f0f0;"|'''Exp'''
 
| align="center" style="background:#f0f0f0;"|'''RCP'''
 
| align="center" style="background:#f0f0f0;"|'''RCP'''
 
| align="center" style="background:#f0f0f0;"|'''AOGCM'''
 
| align="center" style="background:#f0f0f0;"|'''AOGCM'''
 
| align="center" style="background:#f0f0f0;"|'''Std/open'''
 
| align="center" style="background:#f0f0f0;"|'''Std/open'''
| align="center" style="background:#f0f0f0;"|'''Forcing Unc'''
+
| align="center" style="background:#f0f0f0;"|'''Ocean Forcing Unc'''
 
| align="center" style="background:#f0f0f0;"|'''Note'''
 
| align="center" style="background:#f0f0f0;"|'''Note'''
 
|- align="center"
 
||0
 
||N/A
 
||N/A
 
||Control
 
||N/A
 
||Model drift evaluation
 
  
 
|- align="center"
 
|- align="center"
 
||1
 
||1
 
||8.5
 
||8.5
|| <AOGCM1>
+
||MIROC5
 
||Open
 
||Open
 
||Medium
 
||Medium
||Expected high SLR
+
||Expected largest response to SMB, median ocean warming
  
 
|- align="center"
 
|- align="center"
 
||2
 
||2
 
||8.5
 
||8.5
|| NorESM
+
||NorESM
 
||Open
 
||Open
 
||Medium
 
||Medium
||Expected low SLR
+
||Low atmosphere change, low ocean warming
  
 
|- align="center"
 
|- align="center"
 
||3
 
||3
 
||2.6
 
||2.6
|| <AOGCM1>
+
||MIROC5
 
||Open
 
||Open
 
||Medium
 
||Medium
||Expected high SLR
+
||Expected largest response to SMB, median ocean warming
  
 
|- align="center"
 
|- align="center"
 
||4
 
||4
 
||8.5
 
||8.5
|| MIROC5
+
||HadGEM2-ES
 
||Open
 
||Open
 
||Medium
 
||Medium
||Expected mid SLR
+
||Expected median response to SMB, median ocean warming
  
 
|- align="center"
 
|- align="center"
 
||5
 
||5
 
||8.5
 
||8.5
|| <AOGCM1>
+
||MIROC5
 
||Standard
 
||Standard
 
||Medium
 
||Medium
||Expected high SLR
+
||Expected largest response to SMB, median ocean warming
  
 
|- align="center"
 
|- align="center"
 
||6
 
||6
 
||8.5
 
||8.5
|| NorESM
+
||NorESM
 
||Standard
 
||Standard
 
||Medium
 
||Medium
||Expected low SLR
+
||Low atmosphere changes, low ocean warming
  
 
|- align="center"
 
|- align="center"
 
||7
 
||7
 
||2.6
 
||2.6
|| <AOGCM1>
+
||MIROC5
 
||Standard
 
||Standard
 
||Medium
 
||Medium
||Expected high SLR
+
||Expected largest response to SMB, median ocean warming
  
 
|- align="center"
 
|- align="center"
 
||8
 
||8
 
||8.5
 
||8.5
|| MIROC5
+
||HadGEM2-ES
 
||Standard
 
||Standard
 
||Medium
 
||Medium
||Expected mid SLR
+
||Expected median response to SMB, median ocean warming
  
 
|- align="center"
 
|- align="center"
 
||9
 
||9
 
||8.5
 
||8.5
|| <AOGCM1>
+
||MIROC5
 
||Standard
 
||Standard
 
||High
 
||High
||Forcing Uncertainty
+
||Ocean Forcing Uncertainty
  
 
|- align="center"
 
|- align="center"
 
||10
 
||10
 
||8.5
 
||8.5
|| <AOGCM1>
+
||MIROC5
 
||Standard
 
||Standard
 
||Low
 
||Low
||Forcing Uncertainty
+
||Ocean Forcing Uncertainty
  
 
|}
 
|}
  
'''[HG: I am wondering about exp 9 and 10. Shouldn't these be done with a medium model, i.e. MIROC5? Otherwise there is a lot of weight on the high end model]'''
+
== Initial state, control run, historical run and projections set up ==
  
== Experiment duration and initial state ==
+
The core and targeted experiments all start on January 2015 and end in December 2100. The start date follows the CMIP6 protocol for projections, while the end date is constrained by the availability of forcing. In many cases, modelers will need to run a short historical run to bring their models from the “initialization date” to the “projection start date” of January 2015.  
To facilitate analysis of the sea-level projections resulting from the ISMIP6 suite of ice sheet model simulations, ISMIP6 uses 1995 -2014 (January 1995 to December 2014) as reference period. The projections start on January 2015 and end in December 2100. The cutoff date is constrained by the availability of forcing data.
+
  
Groups can reuse their initMIP initial state. However, if the assigned date of the initial state is earlier than 2014, modellers have to find a way to prepare an updated starting state for end 2014.
+
The “initialization date” (or initial state) is left to the modeling groups discretion and can be any time prior to January 2015. The “initialization date” corresponds to the date assigned to the initialization procedure. Groups can reuse their initMIP initialization configuration or generate a new initial state. In the later case, it is important to redo the initMIP schematic experiments ([http://www.climate-cryosphere.org/wiki/index.php?title=InitMIP-Greenland 'asmb'], Goelzer et al. 2018), as it will help understanding how a novel initial state contributes to the uncertainty in ice sheet evolution.  
  
If modellers chose to produce a new initial state, rerunning the [http://www.climate-cryosphere.org/wiki/index.php?title=InitMIP-Greenland initMIP-Greenland] experiments ('ctrl', 'asmb') is required to place results with this new model version into context.
+
Note that ISMIP6 has a new standard grid (horizontal resolution of 1 km on projection EPSG:3413), which is different than that used in initMIP-Greenland. The revised datasets for the initMIP 'asmb' forcing for the new grid can be obtained from:
 +
sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6/Projections/GrIS/initMIP/dSMB_epsg3413  We have changed to the new grid  because its projection is standardised and used by many observational datasets that are key input for ice sheet models.  
  
== Control Run and schematic forcing ==
+
 
The control run ('ctrl') is needed to evaluate model drift. As in the initMIP setup, the control run is obtained by running the model forward without any anomaly forcing, such that whatever surface mass balance (SMB) and ocean forcing was used in the initialization technique would continue unchanged. The second initMIP experiment ('asmb') that also needs to be repeated for each new model version is forcing with a schematic SMB anomaly.
+
A control run (‘ctrl’) is also needed to evaluate model drift. As for initMIP, the control run is obtained by running the model forward, keeping the surface mass balance and ocean forcing used in the initialization technique unchanged. The control run starts from the initial state (typically before 2015) and should last a minimum of 100 years, the same duration as the schematic initMIP experiments. The control run should also be sufficiently long to reach 2100. (See examples in Table 2).
 +
 
 +
Note that in the event that initMIP schematic experiments and control run are redone as part of the projection setup, then consistency with the projections protocol is more important than consistency with the original initMIP setup. For example, in initMIP bedrock was not allowed to evolve. However, an ISM planning to run the projection with evolving bedrock AND planning to redo the initial state, would also rerun initMIP (ctrl,  asmb) with bedrock change. Similarly, the original initMIP requested 2D output every 5 years, whereas the projections protocol request for yearly values. Therefore if rerunning initMIP, the schematic experiments should be saved yearly.
 +
 
 +
A '''single''' "historical run" is required for each ice sheet model, from which all the projections will branch off. The historical run starts at the initMIP state and ends in December 2014. Groups are free to choose how to run the “historical run”, using a reanalysis, a historical run from an RCM, a historical run from an AOGCM or combination of multiple datasets. However, we reject multiple historical runs for each individual AOGCM, because this would very much complicate the forcing strategy and  interpretation. As a consequence, in case AOGCM data is used, please decide for one. ISMIP6 provides a climatology for the SMB and surface temperature for each of the AOGCMs used to generate the projection dataset, as well as anomalies. Groups that would like to use an ISMIP6 Greenland dataset for initialization and historical run are recommended to use MIROC5. Note that the SMB and temperature climatology corresponds to January 1960-December 1989 because the ice sheet is assumed in steady state with the climate during this period. To give a concrete example of a historical run that seems to work well: The ice sheet is initialised to a steady state in 1960. RCM anomalies (here RACMO) are applied from 1960 - 2014 on top of the SMB assumed during initialisation (here RACMO average 1960-1990). Other SMB products forced with reanalysis data (e.g. MAR) are expected to work equally well.
 +
 
 +
The "projection control" (ctrl_proj) is an unforced simulation that starts at the same ice sheet state as the projections. It will run until 2100, and is implemented with zero anomalies. It is meant to capture how much drift arises from the historical run.
 +
 
 +
In most cases the same historical run (and potentially spin-up) can be used for both the standard and open experiments. However, if the implementation of the open and standard experiments requires changes to the ISM that are substantially different (in terms of physics in the ISM), then modelers are allowed to carry out an historical run for the open experiments and an historical run for the standard experiments.
 +
 
 +
{| class="wikitable" style="margin: 1em auto 1em auto;"
 +
! colspan="5"| Table 2: Initialization experiments and examples of different initialization start date
 +
 
 +
|-
 +
|align="center" style="background:#f0f0f0;"|'''Experiment'''
 +
| align="center" style="background:#f0f0f0;"|'''Note'''
 +
|align="center" style="background:#f0f0f0;"|'''start 1 (duration)'''
 +
|align="center" style="background:#f0f0f0;"|'''start 2 (duration)'''
 +
|align="center" style="background:#f0f0f0;"|'''start 3 (duration)'''
 +
 
 +
|- align="center"
 +
||ctrl
 +
||Unforced control run, needed for model drift evaluation
 +
||2015  (100 years)
 +
||2005 (100 years)
 +
||1980 (120 years)
 +
 
 +
|- align="center"
 +
||asmb*
 +
||initMIP prescribed surface mass balance anomaly
 +
||2015 (100 years)
 +
||2005 (100 years)
 +
||1980 (100 years)
 +
 
 +
|- align="center"
 +
||historical
 +
||needed to bring model from initial state to projection start date
 +
||N/A (0 years)
 +
||2005 (10 years)
 +
||1980 (35 years)
 +
 
 +
|- align="center"
 +
||ctrl_proj
 +
||Unforced control run, starting from January 2015, needed for model drift evaluation following historical
 +
||2015 (85 years)
 +
||2015 (85 years)
 +
||2015 (85 years)
 +
|-
 +
! colspan="5"| *only needed if initial state is different from initMIP
 +
 
 +
|}
  
 
== Atmospheric forcing: SMB and temperature anomalies ==
 
== Atmospheric forcing: SMB and temperature anomalies ==
 
=== Introduction ===  
 
=== Introduction ===  
'''General intro on forcing approach, why anomalies of SMB? 1. More robust forcing signal. 2. applicable to a large ensemble of models'''
+
ISMIP6 provides anomalies of SMB and temperature, along with associated climatology, as it allows for an experimental framework that is applicable to a diverse set of ice sheet models.
  
Before applying SMB anomalies, ISMs will need to be initialized by applying a baseline SMB (either a time series or a climatology).  ISMIP6 will provide SMB climatologies for the reference period (January 1995 to December 2014) from the same models computing the anomalies.  ISMs can use these these climatologies for spin-up, if desired, but are free to use their own preferred SMB forcing.   
+
Before applying SMB anomalies, ISMs need to be initialized by applying a baseline SMB (either a time series or a climatology).  ISMIP6 provides SMB climatologies for the reference period (January 1960 to December 1989) from the same models computing the anomalies.  This reference period for Greenland SMB was chosen because the ice sheet is assumed to be in steady state with the surrounding climate over that period. ISMs can use these climatologies for spin-up, if desired, but are free to use their own preferred SMB forcing.   
  
 
Elevation feedbacks have been shown to be important in century-scale simulations.  For example, Le clec’h et al. (2019) considered differences between no coupling (SMB independent of z), one-way coupling (i.e., correcting SMB outputs for ISM topography changes), and two-way coupling (allowing ice-sheet topographic changes to feed back on the RCM simulation.)  Their results suggest that one-way coupling is sufficient to represent elevation feedbacks until the end of the 21st century.  For large topographic changes on longer timescales, it might be necessary to incorporate two-way feedbacks, which are beyond the scope of the standalone ISM effort (but will be considered in the ISMIP6 coupled GCM-ISM experiments).
 
Elevation feedbacks have been shown to be important in century-scale simulations.  For example, Le clec’h et al. (2019) considered differences between no coupling (SMB independent of z), one-way coupling (i.e., correcting SMB outputs for ISM topography changes), and two-way coupling (allowing ice-sheet topographic changes to feed back on the RCM simulation.)  Their results suggest that one-way coupling is sufficient to represent elevation feedbacks until the end of the 21st century.  For large topographic changes on longer timescales, it might be necessary to incorporate two-way feedbacks, which are beyond the scope of the standalone ISM effort (but will be considered in the ISMIP6 coupled GCM-ISM experiments).
  
 
=== Provided forcing data sets ===
 
=== Provided forcing data sets ===
ISMIP6 will provide surface forcing datasets for the Greenland ice sheet (GrIS) based on CMIP global climate model (GCM) simulations. The GCM output is  re-interpreted through higher-resolution regional climate models (RCMs). The later allows to capture narrow regions at the periphery of the Greenland ice sheet with large surface mass balance (SMB) gradients, which are not captured by CMIP5 GCMs.  
+
ISMIP6 provides surface forcing datasets for the Greenland ice sheet (GrIS) based on CMIP AOGCM simulations. The AOGCM output is  re-interpreted through higher-resolution regional climate models (RCMs). The later allows to capture narrow regions at the periphery of the Greenland ice sheet with large surface mass balance (SMB) gradients, which are not captured by CMIP5 AOGCMs. For CMIP6, many of the AOGCMs that have indicated participation in ISMIP6 now use multiple elevation classes to downscale SMB to finer grid resolution. Once these models have completed the CMIP6 projections, our goal is to include additional ISMIP6 projections using SMB downscaled via elevation class.   
 
+
For CMIP6, many of the GCMs that have indicated participation in ISMIP6 now use multiple elevation classes to downscale SMB to finer grid resolution. Once these models have completed the CMIP6 projections, our goal is to include additional ISMIP6 projections using SMB downscaled via elevation class.  '''HG:Move to somewhere else?'''
+
 
   
 
   
For the ISMIP6 projections based on CMIP5 GCMs, the  surface forcing datasets were prepared by Xavier Fettweis, using the MAR regional climate model. RCM downscaling can take account of modest future topography/extent changes, and thus cope with the fact that individual ISM runs may not use exactly the same geometry. All RCM runs use a fixed topography, but vertical SMB gradients for each grid cell are derived using the method described by Franco et al. (2012).  At each location the vertical gradient in SMB is found by summing and averaging pairwise differences in the nearest neighbor cells.  This vertical gradient is used to downscale the SMB (15 km) to a finer grid (1 km), allowing resolution of steep topography that is not represented accurately on the coarse grid. The same information will be used to parameterise the SMB-height feedback in the projections. In addition, MAR calculates a potential SMB term for areas that are outside of the observed ice sheet extent, allowing application to ISMs with ice lying outside the MAR ice-sheet mask. However, experience with initMIP has shown that large variations in ice-sheet extent can lead to considerable bias in the projections. We therefore propose a method (see below) for models with large difference from the observed ice sheet extent, to remap the SMB anomaly to the individual modelled ice sheet geometry.     
+
For the ISMIP6 projections based on CMIP5 AOGCMs, the  surface forcing datasets were prepared by Xavier Fettweis, using the MAR regional climate model. RCM downscaling can take account of modest future topography/extent changes, and thus cope with the fact that individual ISM runs may not use exactly the same geometry. All RCM runs use a fixed topography, but vertical SMB gradients for each grid cell are derived using the method described by Franco et al. (2012).  At each location the vertical gradient in SMB is found by summing and averaging pairwise differences in the nearest neighbor cells.  This vertical gradient is used to downscale the SMB (15 km) to a finer grid (1 km), allowing resolution of steep topography that is not represented accurately on the coarse grid. The same information is used to parameterise the SMB-height feedback in the projections. In addition, MAR calculates a potential SMB term for areas that are outside of the observed ice sheet extent, allowing application to ISMs with ice lying outside the MAR ice-sheet mask. However, experience with initMIP has shown that large variations in ice-sheet extent can lead to considerable bias in the projections. We therefore propose a method (see below) for models with large difference from the observed ice sheet extent, to remap the SMB anomaly to the individual modelled ice sheet geometry.     
  
The atmospheric forcings  consist of annual anomalies of SMB and surface temperature, along with a fields to represent the dependence of SMB and surface temperature on elevation (dSMBdz, dTdz).  SMB is needed by ISMs to compute mass changes at the surface, and surface temperature (i.e., the ice temperature at the base of the snow, as distinct from the 2-m air temperature or skin temperature) is used by many ISMs as an upper boundary condition in the ice temperature calculations.  
+
The atmospheric forcings  consist of annual anomalies of SMB (aSMB) and surface temperature, along with a fields to represent the dependence of SMB and surface temperature on elevation (dSMBdz, dTdz).  SMB is needed by ISMs to compute mass changes at the surface, and surface temperature (i.e., the ice temperature at the base of the snow, as distinct from the 2-m air temperature or skin temperature) is used by many ISMs as an upper boundary condition in the ice temperature calculations. ISMIP6 also provides the climatologies that were used to calculate the anomalies: aSMB_AOGCM(x,y,t) = SMB_AOGCM(x,y,t) - SMB_CLIM_AOGCM(x,y), where SMB_AOGCM is the downscaled SMB for a given AOGCM (using MAR RCM) and SMB_CLIM_AOGCM is the corresponding climatology. The SMB_CLIM_AOGCM were computed by taking the mean value of all SMB_AOGCM over the reference period (from January 1960 to December 1989). ISMs can use these climatologies for spin-up, if desired, but are free to use their own preferred SMB forcing for spin-up and historical run. 
  
 
The SMB anomaly aSMB is given in units [kg m-2 s-1] in yearly values, one year per file. It should be applied constant over a full year and step change at the beginning of a new year. To convert to units [m yr-1] typically used in an ice sheet model, multiply the netcdf variable by 31556926 s/yr, 1/1000 m3/kg and by the density ratio rhow/rhoi:
 
The SMB anomaly aSMB is given in units [kg m-2 s-1] in yearly values, one year per file. It should be applied constant over a full year and step change at the beginning of a new year. To convert to units [m yr-1] typically used in an ice sheet model, multiply the netcdf variable by 31556926 s/yr, 1/1000 m3/kg and by the density ratio rhow/rhoi:
Line 159: Line 202:
 
aSMB [m yr-1] = aSMB [kg m-2 s-1] * 31556926 / 1000 * (1000/rhoi), where rhoi is your specific ice density (typically 917.0 or similar).
 
aSMB [m yr-1] = aSMB [kg m-2 s-1] * 31556926 / 1000 * (1000/rhoi), where rhoi is your specific ice density (typically 917.0 or similar).
  
The SMB and its anomalies will be provided on the ISMIP6 1 km standard ice sheet grid for Greenland. ISMs will then horizontally interpolate the anomaly forcing conservatively from the standard grid to their native grids.
+
The SMB climatology and its anomalies are provided on the ISMIP6 1 km standard ice sheet grid for Greenland. ISMs then horizontally interpolate the anomaly forcing conservatively from the standard grid to their native grids.
  
The SMB change with surface elevation dSMBdz is given in units [kg m-2 s-1 m-1] in yearly values, one year per file. To parameterise the SMB-height feedback, the SMB has to be corrected by dSMBdz * h-h_ref, updated every full year, where h and h_ref are the modelled surface elevation the initial modelled ice sheet surface elevation, respectively. 
+
The SMB change with surface elevation dSMBdz is given in units [kg m-2 s-1 m-1] in yearly values, one year per file. To parameterise the SMB-height feedback, the SMB has to be corrected by dSMBdz * h-h_ref, updated every full year, where h is the time evolving modelled ice sheet surface elevation and h_ref is the modelled surface elevation in December 2014.
 
+
The ISMIP6 standard grid has a horizontal resolution of 1 km and uses the projection EPSG:3413. Note that this grid is different than that used in initMIP-Greenland. We have changed to the new grid  because it's projection is used by many observational datasets that are key input for ice sheet models.
+
  
 
=== Implementation ===
 
=== Implementation ===
  
Let SMB_ref(x,y) denote the SMB used to initialize the ISM, and let h_ref(x,y) denote the ice sheet surface elevation at the end of the initialization.  If a time-dependent SMB is used for spin-up, then SMB_ref(x,y) is the average over the reference period.
+
Let SMB_ref(x,y) denote the SMB used to initialize the ISM, and let h_ref(x,y) denote the ice sheet surface elevation at the end of the initialization.  If a time-dependent SMB is used for spin-up, then SMB_ref(x,y) is the average over the SMB reference period (from January 1960 to December 1989 for Greenland).  If an SMB climatology is used in the assimilation or spin-up, then SMB_ref(x,y) is simply the climatology. ISMIP6 recognizes that the use of existing climatologies (or dataset of SMB averaged over many year) used in the initialization may not align with the time period for the SMB reference period. However, it is assumed that the differences between climatologies will be less than the inter annual variability from the SMB resulting from the AOGCMs, and thus changes in aSMB. What is important is that the SMB climatology (or SMB_ref) is computed over many years. This assumption also allows for ISM to work with their favorite SMB.
  
 
Here we propose two different methods for implementation of atmospheric forcing, depending on how close the ice sheet mask at the end of the initialization is to the observed ice sheet mask, the one assumed by the RCM.  
 
Here we propose two different methods for implementation of atmospheric forcing, depending on how close the ice sheet mask at the end of the initialization is to the observed ice sheet mask, the one assumed by the RCM.  
Line 173: Line 214:
 
'''Method 1: when the ice sheet mask is similar to the observed'''  
 
'''Method 1: when the ice sheet mask is similar to the observed'''  
  
This will typically be the case for ice sheet models that use data assimilation in their initialization, but could be the case for other modelling approaches. For those models, ISMIP6 provides aSMB(x,y,t) at h_rcm(x,y), along with dSMBdz(x,y,t).  Here, aSMB is the time-dependent SMB anomaly in a changing climate, computed in an RCM with fixed surface topography h_rcm, and dSMBdz(x,y,t) is the time-dependent vertical gradient of SMB. aSMB and dSMBdz will be provided on an annual basis and should be updated every full year.
+
This is typically the case for ice sheet models that use data assimilation in their initialization, but could be the case for other modelling approaches. For those models, ISMIP6 provides aSMB(x,y,t) at h_rcm(x,y), along with dSMBdz(x,y,t).  Here, aSMB is the time-dependent SMB anomaly in a changing climate, computed in an RCM with fixed surface topography h_rcm, and dSMBdz(x,y,t) is the time-dependent vertical gradient of SMB. aSMB and dSMBdz are provided on an annual basis and should be updated every full year.
  
Given aSMB(x,y,t) and dSMBdz(x,y,t) on the standard grid, the ISM will horizontally interpolate these fields to its local grid.  Then during runtime, the SMB at a given time and location is computed as
+
Given aSMB(x,y,t) and dSMBdz(x,y,t) on the standard grid, the ISM horizontally interpolate these fields to its local grid.  Then during runtime, the SMB at a given time and location is computed as
  
 
SMB(x,y,t) = SMB_ref(x,y) + aSMB(x,y,t) + dSMBdz(x,y) * [h(x,y,t) - h_ref(x,y)],
 
SMB(x,y,t) = SMB_ref(x,y) + aSMB(x,y,t) + dSMBdz(x,y) * [h(x,y,t) - h_ref(x,y)],
  
where h(x,y,t) is the time-dependent surface elevation.  ISMs will likely need to implement code changes to handle the lapse-rate correction. The models will not need h_rcm(x,y) to compute SMB, but it is provided for reference.  Since dSMBdz is computed at h_rcm and is not given as a function of z, this approach may be inaccurate if h_ref is significantly different from h_rcm.
+
where h(x,y,t) is the time-dependent surface elevation.  ISMs will likely need to implement code changes to handle the lapse-rate correction. The models do not need h_rcm(x,y) to compute SMB, but it is provided for reference.  Since dSMBdz is computed at h_rcm and is not given as a function of z, this approach may be inaccurate if h_ref is significantly different from h_rcm.
  
The datasets can be obtained via the ISMIP6 ftp server (email ismip6@gmail.com to obtain the login information).
+
The datasets of SMB (aSMB, dSMBdz), and surface temperature (aST, dSTdz) from 1950 to 2100, along with the climatology (1960-1989) for MIROC5-rcp8.5, MIROC5-rcp2.6, and NorESM1-rcp8.5 can be obtained via the ISMIP6 ftp server (email ismip6@gmail.com to obtain the login information) at
  
A Test data set is available for MAR3.9-MIROC5-RCP85 at ftp://cryoftp1.gsfc.nasa.gov/ISMIP6/Atmosphere_Forcing/Greenland/
+
sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6/Projections/GrIS/Atmosphere_Forcing/aSMB_observed/
  
 +
Unless advised otherwise, always use the latest version of the forcing files (v1, on Feb 25, 2019)
  
 
'''Method 2: when the ice sheet mask is very different from the observed'''  
 
'''Method 2: when the ice sheet mask is very different from the observed'''  
  
This will typically be the case for ice sheet models that use a glacial-interglacial spinup in their initialization, but could also be the case e.g. for other models that fully relax to a suboptimal SMB. For those models, ISMIP6 will generate a time-dependent SMB anomaly, aSMB(x,y,t) and dSMBdz(x,y,t) that are applied as described above. However, the main difference is that the forcing files are specific for the geometry of your modelled initial state.  
+
This is typically the case for ice sheet models that use a glacial-interglacial spinup in their initialization, but could also be the case e.g. for other models that fully relax to a suboptimal SMB. For those models, ISMIP6 generates a time-dependent SMB anomaly, aSMB(x,y,t) and dSMBdz(x,y,t) that are applied as described above. However, the main difference is that the forcing files are specific for the geometry of your modelled initial state.  
  
 
In order to make the forcing applicable for different ice sheet geometries, we first translate a given SMB anomaly field as a function of absolute location, to a function of surface elevation for 25 regional drainage basins. This step exploits the strong elevation dependence of aSMB. We can then remap aSMB (and dSMBdz) to different modelled geometries. This preserves the overall aSMB patterns and reduces unphysical biases. The procedure to generate the remapped aSMB and dSMBdz is described in Goelzer et al. (2019, in prep).  
 
In order to make the forcing applicable for different ice sheet geometries, we first translate a given SMB anomaly field as a function of absolute location, to a function of surface elevation for 25 regional drainage basins. This step exploits the strong elevation dependence of aSMB. We can then remap aSMB (and dSMBdz) to different modelled geometries. This preserves the overall aSMB patterns and reduces unphysical biases. The procedure to generate the remapped aSMB and dSMBdz is described in Goelzer et al. (2019, in prep).  
  
 
To obtain the datasets, modelers simply need to provide ISMIP6 with their modelled initial surface elevation h_ref(x,y) and ice mask sftgif(x,y). ISMIP6 will then compute aSMB for you. Please send an email to ismip6-at-gmail.com when you are ready to upload your initial state (h_ref, sftgif).  
 
To obtain the datasets, modelers simply need to provide ISMIP6 with their modelled initial surface elevation h_ref(x,y) and ice mask sftgif(x,y). ISMIP6 will then compute aSMB for you. Please send an email to ismip6-at-gmail.com when you are ready to upload your initial state (h_ref, sftgif).  
 +
 +
An example data set produced for the observed geometry as h_ref(x,y) is available for MAR3.9-MIROC5-RCP85 at
 +
 +
sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6/Projections/GrIS/Atmosphere_Forcing/aSMB_remapped/
  
 
The decision on which forcing method to use depends on the expected biases inherent to both approaches. Given the initial state (h_ref and sftgif) we can estimate the biases with a simple integration of the SMB anomaly for a static case (i.e. no ice dynamics). If you are not sure which forcing strategy is best suited for you model, please contact us by email to ismip6-at-gmail.com.
 
The decision on which forcing method to use depends on the expected biases inherent to both approaches. Given the initial state (h_ref and sftgif) we can estimate the biases with a simple integration of the SMB anomaly for a static case (i.e. no ice dynamics). If you are not sure which forcing strategy is best suited for you model, please contact us by email to ismip6-at-gmail.com.
  
 
== Oceanic forcing: Calving and frontal melt ==
 
== Oceanic forcing: Calving and frontal melt ==
ISMIP6 provides dataset of runoff and ocean thermal forcing for models that have their own methods for implementing oceanic induced retreat. In addition, modeling groups are encouraged to participate with the ''ISMIP6 Standard approach'' described below. The later is a simple retreat intended to be easily implemented by the majority of ISM taking part in ISMIP6. In addition, for models that wish to implement a more complex oceanic forcing, the ISMIP6 Greenland ocean focus group has developed a second methodology, described in ''ISMIP6 high resolution ocean melt rate approach''.   
+
ISMIP6 provides dataset of runoff and ocean thermal forcing for models that have their own methods for implementing oceanic induced retreat. In addition, modeling groups are expected to participate with one of the two ISMIP6 approaches described below. The  ''ISMIP6 Standard approach'' is a simple retreat intended to be easily implemented by the majority of ISM taking part in ISMIP6. Alternatively, for models that wish to implement a more complex oceanic forcing, the ISMIP6 Greenland ocean focus group has developed a second methodology, described in ''ISMIP6 high resolution ocean melt rate approach''.   
  
'''ISMIP6 Standard approach''' imposes an empirically-derived, sector averaged retreat (Fig 2) as a function of climate forcing. This method was developed for ISMIP6 as a result of the ocean forcing focus group and is described in greater details in Slater et al. (in prep), and in the webinar:
+
'''ISMIP6 Standard approach''' imposes an empirically-derived, sector averaged retreat (Fig 2) as a function of climate forcing. This method was developed for ISMIP6 as a result of the ocean forcing focus group and is described in greater details in Slater et al. (2019a; 2019b), and in the webinar:
ftp://cryoftp1.gsfc.nasa.gov/ISMIP6/Ocean_Forcing/Greenland/Webinar_2018_11_08
+
sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6/Projections/GrIS/Ocean_Forcing/Webinar_2018_11_08/
  
 +
[[File:Oceans_overview_greenland.png|720x480px]]
 
[[File:retreat_rate.png|720x480px]]
 
[[File:retreat_rate.png|720x480px]]
  
Fig 2: Example of the empirically-derived retreat scenarios for the 7 sectors of the Greenland ice sheet, obtained with MIROC5, RCP8.5. '''TO DO: MODIFY FIG so that it shows the 7 ISMIP6 Basins
+
Fig 2: Example of the empirically-derived retreat scenarios for the 7 sectors of the Greenland ice sheet, obtained with MIROC5, RCP8.5.  
'''
+
  
As described in the webinar, retreat will be imposed when the ice sheet geometry and ice front retreat scenario indicates that the land_ice_area_fraction_retreat mask is ice free for a given year. (Note the name was chosen so that it is closely related to the standard name land_ice_area_fraction -corresponding to variable name sftgif- in the ISMIP6 data request. For ease of communication, we use the longer standard name). Implementation for a specific ISM (also illustrated in Fig 3) requires the following steps:
+
As described in the webinar, retreat is imposed when the ice sheet geometry and ice front retreat scenario indicates that the land_ice_area_fraction_retreat mask is ice free for a given year. (Note the name was chosen so that it is closely related to the standard name land_ice_area_fraction -corresponding to variable name sftgif- in the ISMIP6 data request. For ease of communication, we use the longer standard name). Implementation for a specific ISM (also illustrated in Fig 3) requires the following steps:
  
 
1. Identify ice prone to outlet glacier retreat, by interpolating initial ice mask conservatively to 1 km ISMIP6 diagnostic grid
 to obtain a mask for ice fraction:
 
1. Identify ice prone to outlet glacier retreat, by interpolating initial ice mask conservatively to 1 km ISMIP6 diagnostic grid
 to obtain a mask for ice fraction:
Line 233: Line 279:
 
ISMIP6 will generate the land_ice_area_fraction_retreat masks (steps 2-3) for each model.  As retreat is provided as a series of ice fraction masks, ice sheet models with coarse resolution should use a sub-grid approach. A suggestion is to apply an land_ice_area_fraction_retreat that is relative to the reference thickness. Models may have a different strategy for this sub-grid implementation.  
 
ISMIP6 will generate the land_ice_area_fraction_retreat masks (steps 2-3) for each model.  As retreat is provided as a series of ice fraction masks, ice sheet models with coarse resolution should use a sub-grid approach. A suggestion is to apply an land_ice_area_fraction_retreat that is relative to the reference thickness. Models may have a different strategy for this sub-grid implementation.  
  
''question: do we want to give an eq for sub-grid implementation?"
+
''Note on the retreat dataset:'' the retreat rate dataset was calibrated using grounding line position of glaciers that do not have ice shelves. Although the retreat dataset is therefore not optimum for glaciers that have a floating tongue, it is suggested that retreat is imposed at the ice front. Modeling groups that have the capability of computing ice tongue basal melt may use the provided dataset of thermal forcing per basin.
  
To obtain the land_ice_area_fraction_retreat masks appropriate for your model, please submit your land_ice_area_fraction mask (step 1) to ISMIP6 on the 1km ISMIP6 grid (EPSG:3413) as a netcdf file called sftgif.nc. This file should be generated by conservative regridding. If your native grid is regular and on EPSG:3413, please provide your original modelled ice mask with x,y information instead. Upload your file to the directory ISMIP6/Ocean_Forcing/Greenland/Retreat_Implementation/MODELFILES/MODELNAME, where MODELNAME is the name of your model, and let us know via email that your file is uploaded.   
+
To obtain the land_ice_area_fraction_retreat masks appropriate for your model, please submit your land_ice_area_fraction mask (step 1) to ISMIP6. This should be your ice mask from end 2014. We distinguish different cases to reduce interpolation artefacts. 1) If your native grid is regular and on EPSG:3413, please provide your original modelled ice mask with x,y information. 2) If you have a high resolution irregular grid or a grid on a different projection than EPSG:3413, please get in touch so we can find the best solution for you. When ready, upload your file to the directory /ISMIP6/Projections/GrIS/Ocean_Forcing/Retreat_Implementation/MODELFILES/MODELNAME, where MODELNAME is the name of your model, and let us know via email that your file is uploaded.   
  
'''ISMIP6 open approach''' is used sample a larger variety of oceanic forcing parameterizations, as it remains an active field of research. Models are free to continue applying the ocean forcing parameterization they used during the model initialization or their preferred method, but should still rely on the ocean forcing datasets provided by ISMIP6 to simulate future ocean conditions.
+
'''ISMIP6 high resolution ocean melt rate''' allows ice sheet models to specify terminus retreat for each individual marine-terminating outlet glacier, thus allowing glaciers to retreat at different rates. This is accomplished by specifying a glacier terminus melt rate, calculated as a function of subglacial discharge (approximated as surface runoff from each glacier catchment) and ocean thermal forcing (Fig. 4). In addition to the melt rate, a calving rate must be specified to obtain the total frontal ablation rate at each terminus. Additional details can be found in the webinar and the updated presentation slides: sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6/Projections/GrIS/Ocean_Forcing/GrIS_Melt_Parameterization_email_2019Jan.pptx
  
'''ISMIP6 high resolution ocean melt rate''' allows ice sheet models to specify terminus retreat for each individual marine-terminating outlet glacier, thus allowing glaciers to retreat at different rates. This is accomplished by specifying a glacier terminus melt rate, calculated as a function of subglacial discharge (approximated as surface runoff from each glacier catchment) and ocean thermal forcing (Fig. 4). In addition to the melt rate, a calving rate must be specified to obtain the total frontal ablation rate at each terminus. The high resolution ocean melt rate is considered an "open approach". Additional details can be found in the webinar and the updated presentation slides: ftp://cryoftp1.gsfc.nasa.gov/ISMIP6/Ocean_Forcing/Greenland/Webinar_2018_11_08
+
Dataset for thermal forcing and runoff from 1950 to 2100 are available for MIROC5-rcp2.6, MIROC5-rcp8.5, and NorESM1-rcp8.5 at:  
 +
sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6/Projections/GrIS/Ocean_Forcing/Melt_Implementation/
 +
As the dataset were created for an ice sheet geometry corresponding to present day observations from GIMP and BedMachine3, ISM that have differ substantially from this initial geometry should not use this dataset. (ie: if you are using Method 2 for the atmospheric forcing/SMB remapping, then the dataset is not applicable for your model, as the observed drainage basins used to generate the runoff field may not work with your model. You could create a dataset that is appropriate for your model from the remapped runoff and a water routing consistent with your geometry).
  
 
[[File:GrIS_Melt_Parameterization_flowchart.png|650x312px]]
 
[[File:GrIS_Melt_Parameterization_flowchart.png|650x312px]]
Line 245: Line 293:
 
Fig 4: High-resolution melt-rate approach flow chart. Ocean thermal forcing and runoff fields are provided to the ice sheet modelers by the ISMIP6 ocean forcing working group. All parameters in the recommended melt-rate parameterization (Xu et al., 2013; Rignot et al., 2016) are provided, as well. Each ice sheet model must implement the approach and specify their own calving parameterization to calculate a frontal ablation rate for each outlet glacier.
 
Fig 4: High-resolution melt-rate approach flow chart. Ocean thermal forcing and runoff fields are provided to the ice sheet modelers by the ISMIP6 ocean forcing working group. All parameters in the recommended melt-rate parameterization (Xu et al., 2013; Rignot et al., 2016) are provided, as well. Each ice sheet model must implement the approach and specify their own calving parameterization to calculate a frontal ablation rate for each outlet glacier.
  
== Requirements for the standard experiments ==  
+
'''ISMIP6 open approach''' is used sample a larger variety of oceanic forcing parameterizations, as it remains an active field of research. Models are free to continue applying the ocean forcing parameterization they used during the model initialization or their preferred method, but should still rely on the ocean forcing datasets provided by ISMIP6 to simulate future ocean conditions.
 +
 
 +
== Requirements for the projections ==  
  
 
• Participants can and are encouraged to contribute with different models and/or initialisation methods
 
• Participants can and are encouraged to contribute with different models and/or initialisation methods
Line 251: Line 301:
 
• Models have to be able to prescribe a given SMB anomaly
 
• Models have to be able to prescribe a given SMB anomaly
  
• Models have to be able to prescribe a given ice front retreat
+
• Models have to be able to prescribe a given ice front retreat or the high resolution approach for the standard experiments. For the open experiments, models can choose the ocean parameterization of their choice but should use the ocean forcing provided.
  
• Adjustment of SMB due to geometric changes in forward experiments is encouraged.  
+
• Adjustment of SMB due to geometric changes in forward experiments is encouraged using the provided dSMB/dz.  
  
 
• Bedrock adjustment in forward experiment is allowed.
 
• Bedrock adjustment in forward experiment is allowed.
Line 259: Line 309:
 
• The choice of model input data is unconstrained to allow participants the use of their preferred model setup without modification. Modelers without preferred data set choice can have a look at the ISMIP6 [[Datasets]] page for possible options.
 
• The choice of model input data is unconstrained to allow participants the use of their preferred model setup without modification. Modelers without preferred data set choice can have a look at the ISMIP6 [[Datasets]] page for possible options.
  
• To allow for analysis, any modeling choice needs to be well documented. A README file needs to be submitted along the outputs as an integral part of the contribution to the ISMIP6. It may be obtained here ('''need to update the readme file''') or requested by email to ismip6-at-gmail.com.
+
• To allow for analysis, any modeling choice needs to be well documented. A README file needs to be submitted along the outputs as an integral part of the contribution to the ISMIP6. It may be obtained [https://www.dropbox.com/s/7ejbqlhkwild4i4/README_GrIS_Projections.docx?dl=0 hereor requested by email to ismip6-at-gmail.com.
  
== Requirements for the open experiments ==  
+
==References ==
• Participants can and are encouraged to contribute with different models and/or initialisation methods
+
  
• Models have to be able to prescribe a given SMB anomaly
+
Alice Barthel, Cecile Agosta, Christopher M. Little, Tore Hatterman, Nicolas C. Jourdain, Heiko Goelzer, Sophie Nowicki, Helene Seroussi, Fiammetta Straneo, and Thomas J. Bracegirdle: CMIP5 model selection for ISMIP6 ice sheet model forcing: Greenland and Antarctica, The Cryosphere Discuss., https://doi.org/10.5194/tc-2019-191, 2019.
  
• Models can choose the ocean parameterization of their choice but this parameterization should use the ocean forcing provided
+
Franco, B., Fettweis, X., Lang, C., and Erpicum, M.: Impact of spatial resolution on the modelling of the Greenland ice sheet surface mass balance between 1990–2010, using the regional climate model MAR, The Cryosphere, 6, 695-711, https://doi.org/10.5194/tc-6-695-2012, 2012.
  
• Adjustment of SMB due to geometric changes in forward experiments is encouraged.  
+
Goelzer, H., Nowicki, S., Edwards, T., Beckley, M., Abe-Ouchi, A., Aschwanden, A., Calov, R., Gagliardini, O., Gillet-Chaulet, F., Golledge, N. R., Gregory, J., Greve, R., Humbert, A., Huybrechts, P., Kennedy, J. H., Larour, E., Lipscomb, W. H., Le clec'h, S., Lee, V., Morlighem, M., Pattyn, F., Payne, A. J., Rodehacke, C., Rückamp, M., Saito, F., Schlegel, N., Seroussi, H., Shepherd, A., Sun, S., van de Wal, R., and Ziemen, F. A.: Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison, The Cryosphere, 12, 1433-1460, https://doi.org/10.5194/tc-12-1433-2018, 2018.
  
• Bedrock adjustment in forward experiment is allowed.
+
Goelzer, H., Noel, B. P. Y., Edwards, T. L., Fettweis, X., Gregory, J. M., Lipscomb, W. H., van de Wal, R. S. W., and van den Broeke, M. R.: Remapping of Greenland ice sheet surface mass balance anomalies for large ensemble sea-level change projections, The Cryosphere Discuss., 2019, 1-20, https://doi.org/10.5194/tc-2019-188, 2019.
  
• The choice of model input data is unconstrained to allow participants the use of their preferred model setup without modification. Modelers without preferred data set choice can have a look at the ISMIP6 [[Datasets]] page for possible options.
+
Le clec'h, S., Charbit, S., Quiquet, A., Fettweis, X., Dumas, C., Kageyama, M., Wyard, C., and Ritz, C.: Assessment of the Greenland ice sheet–atmosphere feedbacks for the next century with a regional atmospheric model coupled to an ice sheet model, The Cryosphere, 13, 373-395, https://doi.org/10.5194/tc-13-373-2019, 2019.
  
• To allow for analysis, any modeling choice needs to be well documented. A README file needs to be submitted along the outputs as an integral part of the contribution to the ISMIP6. It may be obtained here ('''need to update the readme file''') or requested by email to ismip6-at-gmail.com.
+
Rignot, E., Xu, Y., Menemenlis, D., Mouginot, J., Scheuchl, B., Li, X., et al.: Modeling of ocean-induced ice melt rates of five west Greenland glaciers over the past two decades. Geophysical Research Letters, 43(12). http://doi.org/10.1002/2016GL068784, 2016.
  
==References ==
+
Shannon, S.R., Payne A.J., Bartholomew I.D., Van Den Broeke M.R., Edwards T.L., Fettweis X., Gagliardini O., Gillet-Chaulet F., Goelzer H., Hoffman M.J., Huybrechts P.: Enhanced basal lubrication and the contribution of the Greenland ice sheet to future sea-level rise, Proceedings of the National Academy of Sciences, 110(35):14156-61. https://doi.org/10.1073/pnas.1212647110, 2013.
Franco, B., Fettweis, X., Lang, C., and Erpicum, M.: Impact of spatial resolution on the modelling of the Greenland ice sheet surface mass balance between 1990–2010, using the regional climate model MAR, The Cryosphere, 6, 695-711, https://doi.org/10.5194/tc-6-695-2012, 2012.
+
  
Goelzer, H., Nowicki, S., Edwards, T., Beckley, M., Abe-Ouchi, A., Aschwanden, A., Calov, R., Gagliardini, O., Gillet-Chaulet, F., Golledge, N. R., Gregory, J., Greve, R., Humbert, A., Huybrechts, P., Kennedy, J. H., Larour, E., Lipscomb, W. H., Le clec'h, S., Lee, V., Morlighem, M., Pattyn, F., Payne, A. J., Rodehacke, C., Rückamp, M., Saito, F., Schlegel, N., Seroussi, H., Shepherd, A., Sun, S., van de Wal, R., and Ziemen, F. A. (2018).
+
Slater, D. A., Straneo, F., Felikson, D., Little, C. M., Goelzer, H., Fettweis, X., and Holte, J.: Estimating Greenland tidewater glacier retreat driven by submarine melting, The Cryosphere, 13, 2489-2509, https://doi.org/10.5194/tc-13-2489-2019, 2019a.
[https://www.the-cryosphere.net/12/1433/2018/tc-12-1433-2018.html Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison], The Cryosphere, 12, 1433-1460, doi:10.5194/tc-12-1433-2018.
+
  
Goelzer, H. et al. (2019) Remapping of Greenland ice sheet surface mass balance anomalies for large ensemble sea-level change projections, in prep.
+
Slater, D. A., Felikson, D., Straneo, F., Goelzer, H., Little, C. M., Morlighem, M., Fettweis, X., and Nowicki, S.: 21st century ocean forcing of the Greenland Ice Sheet for modeling of sea level contribution, The Cryosphere Discuss., 2019, 1-34, https://doi.org/10.5194/tc-2019-222, 2019b.  
  
Le clec'h, S., Charbit, S., Quiquet, A., Fettweis, X., Dumas, C., Kageyama, M., Wyard, C., and Ritz, C.: Assessment of the Greenland ice sheet–atmosphere feedbacks for the next century with a regional atmospheric model coupled to an ice sheet model, The Cryosphere, 13, 373-395, https://doi.org/10.5194/tc-13-373-2019, 2019.
+
Xu, Y., Rignot, E., Fenty, I., Menemenlis, D., & Flexas, M. M.: Subaqueous melting of Store Glacier, west Greenland from three-dimensional, high-resolution numerical modeling and ocean observations. Geophysical Research Letters, 40(17). http://doi.org/10.1002/grl.50825, 2013.
  
Shannon, S.R., Payne A.J., Bartholomew I.D., Van Den Broeke M.R., Edwards T.L., Fettweis X., Gagliardini O., Gillet-Chaulet F., Goelzer H., Hoffman M.J., Huybrechts P. (2013) Enhanced basal lubrication and the contribution of the Greenland ice sheet to future sea-level rise, Proceedings of the National Academy of Sciences, 110(35):14156-61.
+
==Acknowledgements==
 +
The experimental protocol and datasets for the ISMIP6-Projections-Greenland standalone ice sheet simulations would not have been possible without the effort of many scientists that have given their time and expertise, and have run models to convert the CMIP5 models output into datasets that standalone ice sheet models can use. ISMIP6 would like to thank the ocean focus group under the leadership of Fiamma Straneo, the atmospheric focus group under the leadership of Bill Lipscomb and Robin Smith, and the CMIP5 model evaluation focus group under the leadership of Alice Barthel. Donald Slater, Denis Felixson, Mathieu Morlinghem and Heiko Goelzer have been instrumental in the development of the ice front retreat and melt parameterization and associated dataset. Xavier Fettweis, Patrick Alexander and Heiko Goelzer prepared the atmospheric dataset. Alice Barthel, Chris Little, Cecile Agosta, and Jamie Holte provided a rigorous analysis of the CMIP5 models against historical data, which allowed the CMIP5 model evaluation group and the ISMIP6 steering committee to select the CMIP5 models used in this effort. Finally, we thank the ISMIP6 ice sheet modelers for their feedback on the design of the protocol and their willingness to participate in ISMIP6.
  
 
==Appendix 1 – Output grid definition and interpolation==
 
==Appendix 1 – Output grid definition and interpolation==
Line 303: Line 351:
 
*  An overview of the regridding process can be found on the [[Regridding]] page.
 
*  An overview of the regridding process can be found on the [[Regridding]] page.
 
*  [[Regridding_with_CDO]] contains tools and tips that have been used by ISMIP6 members
 
*  [[Regridding_with_CDO]] contains tools and tips that have been used by ISMIP6 members
 +
*  ISMIP6 is designing tools to help with the regridding.
 +
*  If you need help with conservative interpolation, please email ismip6-at-gmail.com.
  
 
==Appendix 2 – Naming conventions, upload and model output data.==
 
==Appendix 2 – Naming conventions, upload and model output data.==
Please provide:
+
'''COMING SOON: GUIDANCE FOR CMIP6 MODEL REGISTRATION, ADDITIONAL GUIDANCE FOR FILE PREPARATIONS and REVISED FILENAME CONVENTION. HOWEVER, THE CURRENT FILE STRUCTURE AND GUIDANCE BELOW IS DESIGNED TO FACILITATE THE FINAL FILE PREPARATION (NEEDED BEFORE UPLOADING TO THE CMIP6 ARCHIVE). THE EXTRA INFORMATION/FILE RENAMING WILL BE IMPLEMENTED BY SCRIPTS CURRENTLY WRITTEN BY ISMIP6, WHICH WILL ALSO CHECK FOR CF COMPLIANCE ETC. WE WILL PROVIDE HELP TO MODELERS FOR THE FINAL FILE FORMATING PREPARATION, AND MODELERS SHOULD PROCEED WITH SAVING THEIR RESULTS USING THE INFORMATION BELOW '''
  
• one variable per file for all 2D fields
+
Please provide:
  
all variables in one file for the scalar variables
+
• one variable per file for all 2D fields and scalar variables.
  
 
• a completed readme file
 
• a completed readme file
  
 
===A2.1 File name convention===
 
===A2.1 File name convention===
'''TO DO: REVISE EXP NAME (LOOK AT CMIP6 GUIDANCE, INCLUDING HOW TO REGISTER MODEL) AND CHANGE INFO'''
+
 
File name convention for 2D fields:
+
File name convention for 2D fields and scalar variables:
  
 
<variable>_<IS>_<GROUP>_<MODEL>_<EXP>.nc
 
<variable>_<IS>_<GROUP>_<MODEL>_<EXP>.nc
 
File name convention for scalar variables:
 
 
scalar_<IS>_<GROUP>_<MODEL>_<EXP>.nc
 
  
 
File name convention for readme file:
 
File name convention for readme file:
Line 337: Line 383:
 
<MODEL> = model acronym (all upper case or numbers, no special characters)
 
<MODEL> = model acronym (all upper case or numbers, no special characters)
  
<EXP> = experiment name (init, ctrl or asmb)
+
<EXP> = experiment name  
  
For example, a file containing the scalar variables for the Greenland ice sheet, submitted by group “JPL” with model “ISSM” for experiment “ctrl” would be called:
+
For example, a file containing the variable "orog" for the Greenland ice sheet, submitted by group “JPL” with model “ISSM” for experiment “ctrl” would be called:
scalar_GIS_JPL_ISSM_ctrl.nc
+
orog_GIS_JPL_ISSM_ctrl.nc
  
 
If JPL repeats the experiments with a different version of the model (for example, by changing the sliding law), it could be named ISSM2, and so forth.
 
If JPL repeats the experiments with a different version of the model (for example, by changing the sliding law), it could be named ISSM2, and so forth.
  
===A2.2 Uploading your model output===
+
When uploading your files to the server, please follow the standard directory structure.  
Please upload your model output on the FTP server cryoftp1.gsfc.nasa.gov, and email ismip6@gmail.com for the user name and latest password. Note sftp does not work!
+
  
'''TO DO: CREATE DIRECTORY FOR UPLOAD AND CHANGE INFO'''
+
Example:
  
After log in, go to the ISMIP6/initMIP/output directory via:
+
  group1
 +
      model1
 +
        exp05_05
 +
        expC01_05
  
ftp> cd /ISMIP6/initMIP/output
+
  group2
 +
      model2
 +
        exp05_01
 +
        expC01_01
 +
 
 +
Group, model and experiment names have to be identical in the directory names and file names, except for the resolution suffix (rr=[01, 05, ..] the resolution the submitted data in km).
 +
 
 +
Example for results at 1 km resolution:
 +
  AWI/ISSM1/ctrl_01/acabf_GIS_AWI_ISSM1_ctrl.nc
 +
  xxx yyyyy zzzz rr          xxx yyyyy zzzz
 +
 
 +
The experiment names (exp_id) can be found at https://docs.google.com/spreadsheets/d/1YEiPV3Uc0K8EqD57mz-ZqBIF4XnVb5UbyRDRmox0KeE/edit#gid=1531275537
 +
 
 +
A few examples:
 +
 
 +
{| class="wikitable"
 +
 
 +
| align="center" style="background:#f0f0f0;"|'''exp_id'''
 +
| align="center" style="background:#f0f0f0;"|'''RCP'''
 +
| align="center" style="background:#f0f0f0;"|'''GCM'''
 +
| align="center" style="background:#f0f0f0;"|'''Ocean'''
 +
|-
 +
|| exp05
 +
|| 8.5
 +
|| MIROC5
 +
|| Medium
 +
|-
 +
|| exp06
 +
|| 8.5
 +
|| NorESM
 +
|| Medium
 +
|-
 +
|| expa01
 +
|| 8.5
 +
|| IPSL-CM5-MR
 +
|| Medium
 +
|-
 +
|| expb01
 +
|| ssp585
 +
|| CNRM-CM6
 +
|| Medium
 +
|-
 +
|| expc01
 +
|| 8.5
 +
|| MIROC5 AO
 +
|| Medium
 +
|-
 +
|| expc01
 +
|| 8.5
 +
|| NorESM
 +
|| High
 +
|}
 +
 
 +
 
 +
===A2.2 Retrieving dataset and Uploading your model output===
 +
====A2.2.1 How to gain access to server====
 +
All ISMIP6 datasets are stored at the University of Buffalo’s CCR transfer server (UB CCR).
 +
 
 +
In order to gain access to the ISMIP6 datasets and upload or download files please send a request to ismip6 [at] gmail.com address with your public ssh key, name, and affiliation for account setup.
 +
 
 +
To generate your public ssh key use the command below:
 +
ssh-keygen -t rsa
 +
 
 +
Once your account is set up you will receive a confirmation email with your username.
 +
 
 +
All the current and previous ISMIP6 related information and datasets are located in the
 +
sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6 folder.
 +
 
 +
You can access the server and upload or download datasets on the UB CCR server the following ways below.
 +
 
 +
'''Command line access with sftp:'''
 +
sftp -o "IdentityFile=sshkey" yourusername@transfer.ccr.buffalo.edu
 +
 
 +
Upload files
 +
put /path/to/your/file /path/to/folder/on/server
 +
 +
Upload folder
 +
put -r /path/to/your/folder /path/to/folder/on/server
 +
 +
Download file
 +
get /path/to/file/on/server /path/to/your/folder/
 +
 +
Download whole directory
 +
get -r /path/to/folder/on/server /path/to/your/folder/
 +
 
 +
'''Commands for rsync ssh access:'''
 +
 +
Copy a file from local folder to remote folder:
 +
rsync -avz -e "ssh -i sshkey" /home/localuser/testfile.txt remoteuser@transfer.ccr.buffalo.edu:/home/remoteuser/transfer/
 +
 
 +
Copy a file from remote server into a local folder:
 +
rsync -avz -e "ssh -i sshkey" remoteuser@transfer.ccr.buffalo.edu:/home/remoteuser/transfer/testfile.txt  /home/localuser/
 +
 
 +
'''FileZilla/Cyberduck Setup'''
 +
 +
Please see link below for easy FileZilla setup
 +
https://www.a2hosting.com/kb/getting-started-guide/accessing-your-account/setting-up-sftp-access-for-multiple-users#Configure-an-FTP-client-to-use-the-private-key
 +
 +
Please set the PORT to 22.
 +
 
 +
====A2.2.2 Where to upload your results====
 +
 
 +
After log in, go to the ISMIP6/Projections/GrIS/output directory via:
 +
sftp> cd /ISMIP6/Projections/GrIS/output
  
 
and create a directory named <GROUP> with the following sub-directory structure:
 
and create a directory named <GROUP> with the following sub-directory structure:
 +
sftp> mkdir /Projections/GrIS/output/<GROUP>/<MODEL>/<EXP_RES>
  
initMIP
+
<EXP_RES> should include both the experiment name and the output grid used to simplify the processing (e.g., asmb_08). Only the directory name should include this resolution, unlike the output files.
output/
+
<GROUP>/
+
<MODEL>/
+
init/
+
ctrl/
+
asmb/
+
  
 
Create additional <MODEL> directories when participating with more than one model or model version.
 
Create additional <MODEL> directories when participating with more than one model or model version.
  
An example of model output files can be found in /ISMIP6/initMIP/output/ISMIP6/REF.
+
====A2.2.3 Reducing the size of files====
 +
The size of the model files on higher resolution grid can be largely reduced by file compression which will save space on the storage server. An example command is given below and the results before and after. In the examples that follow we can get a factor of 10 compression and for the masks even more given that contiguous masks are highly compressible because they are repeated data. NetCDF files have been designed with compression in mind. A NetCDF file can be compressed and nothing has to be changed in the way that it is read into Matlab or Python (or any other language that uses standard NetCDF read/write libraries).
 +
 +
The nccopy command copies an input netCDF file to an output netCDF file after compressing the file significantly. The ‘-d’ option stands for the deflation level, from 1 (faster but lower compression) to 9 (slower but more compression) and the ‘-s’ option is the shuffling option to improve compression even more. We recommend using ‘d1’ option since this option seems to accomplish the desired compression.
  
===A2.3 Model output variables and README file===
+
Example of netcdf compression command:
 +
nccopy -d1 -s sftgif_GIS_JPL_ISSMPALEO_historical.nc sftgif_GIS_JPL_ISSMPALEO_historical_c.nc
  
The README file is an important contribution to the ISMIP6 submission. It may be obtained [https://www.dropbox.com/s/tkot25wq5n6ig1b/README_%3CIS%3E_%3CGROUP%3E_%3CMODEL%3E.doc?dl=0| here] or requested by email to ismip6-at-gmail.com
+
Example of compression variant, seems to work better for masks:
 +
nccopy -d1 sftgif_GIS_JPL_ISSMPALEO_historical.nc sftgif_GIS_JPL_ISSMPALEO_historical_c.nc
  
 +
===A2.3 Model output variables and README file===
  
 +
The README file is an important contribution to the ISMIP6 submission. It may be obtained [https://www.dropbox.com/s/7ejbqlhkwild4i4/README_GrIS_Projections.docx?dl=0| here] or requested by email to ismip6-at-gmail.com
 +
 +
====A2.3.1 General guidelines====
 
The variables requested in the table below serve to evaluate and compare the different models and initialization techniques. Some of the variables may not be applicable for your model, in which case they are to be omitted (with explanation in the README file). Also, specify missing values in your netcdf file where needed, and fields should be undefined outside of the ice mask.  
 
The variables requested in the table below serve to evaluate and compare the different models and initialization techniques. Some of the variables may not be applicable for your model, in which case they are to be omitted (with explanation in the README file). Also, specify missing values in your netcdf file where needed, and fields should be undefined outside of the ice mask.  
  
We distinguish between state variables (e.g. ice thickness, temperatures and velocities) and flux variables (e.g. SMB). Flux variables are defined as positive when the process adds mass to the ice sheet and negative otherwise. Note the different treatment for state variables (snapshots) and fluxes (time average). The standard should be averaging over all native time steps for yearly scalar output and for 5 year periods for 2D fields. Please specify how your reported flux data has been averaged over time in the README file.
+
All "missing data" must be assigned the single precision floating point value of 1.e20
  
Example model output files can be found in /ISMIP6/initMIP/output/ISMIP6/REF on the ftp server.
+
We distinguish between state variables "ST" (e.g. ice thickness, temperatures and velocities) and flux variables "FL" (e.g. SMB). State variables should be given as snapshot information at the end of one year for both scalars and 2D variables (for initMIP, 2D variables were only requested over five year periods), while flux variables are to be averaged over the respective periods. Please specify in your README file how your reported flux data has been averaged over time. Ideally, the standard would be go average over all native time steps.
 +
 
 +
Flux variables are defined positive when the process adds mass to the ice sheet and negative otherwise.
 +
 
 +
If you redo the initMIP experiments (because you have a new initial state for the projections), please save the files at a yearly interval instead of the 5 years interval requested as part of the original initMIP. Also upload your iniMIP results in the projections directory.
 +
 
 +
====A2.3.2 How to record time in historical and projection files====
 +
 
 +
In compliance with CMIP6, time should be defined in "days since <basetime>", where <basetime> must be specified by the user, typically in the form year-month-day (e.g., "days since 1800-1-1"). For simulations meant to represent a particular historical period, set the ‘base time’ to the time at the beginning of the simulation.  A historical run initialized with forcing for year 2007 would, for example, have units of “days since 2007-1-1”.  For the future scenario runs, retain the same <basetime> as used in the historical run from which it was initiated. Note the CF definition for years ([http://cfconventions.org/Data/cf-conventions/cf-conventions-1.7/cf-conventions.html#paleoclimate-time-axis-ex. section 4.4]): a common_year is 365 days, a leap_year is 366 days, a Julian_year is 365.25 days,  a Gregorian_year is 365.2425 days, a 360_day has all years with 360 days divided into 30 day months (please see the CF link above for other examples on calendar setting in section 4.4).
 +
 
 +
To illustrate a time recording for the historical file and projections for a typical state variable (ST, eg thickness) and flux variable (FL, eg SMB), we assume that our <basetime> is January 1st 2013, and that we use a calendar = 360_days. '''Other calendars can be used, but you need to indicate the calendar used in the netcdf, and of course if you use a different calendar, the time entries will be different.''' What needs to be recorded is shown in green in the Table below. For state variables, <time> is the day corresponding to the entry that you are saving since the your <basetime>. For flux variables, since these are averaged over a year, <time> is the day since <basetime> corresponding to the middle of the year, while <time_bnds> records the day since <basetime> at the start and end of a year. Note that in CMIP a full year is typically from first of January to the first of January of the following year. We also provide below example of what the  netcdf would look like for our example..
 +
 
 +
[[File:Figure_A1.png|951x471px]]
 +
 
 +
 
 +
For state variables, like thickness for the historical:
 +
dimensions:
 +
    time = UNLIMITED ; // (3 currently)
 +
variables:
 +
    double time(time) ;
 +
        time:units = "days since 1-1-2013" ; // This date correspond to the example basetime
 +
        time:calendar = "360_day" ; // Other calendars can be used... change here to relevant calendar
 +
        time:axis = "T" ;
 +
        time:long_name = "time" ;
 +
        time:standard_name = "time" ;
 +
data:
 +
  time = 0, 360, 720; // If you use a different calendar these values will change
 +
 
 +
and thickness for the projection (note that the full time entries are not shown, only beginning and end) would be:
 +
dimensions:
 +
    time = UNLIMITED ;
 +
variables:
 +
    double time(time) ;
 +
        time:units = "days since 1-1-2013" ; // This date correspond to the example basetime
 +
        time:calendar = "360_day" ; // Other calendars can be used... change here to relevant calendar
 +
        time:axis = "T" ;
 +
        time:long_name = "time" ;
 +
        time:standard_name = "time" ;
 +
data:
 +
  time = 720, 1080, 1440, …, 31320, 31680; // If you use a different calendar these values will change
 +
 
 +
The flux variable, like SMB, would be recorded as the average over a full year, so for the historical:
 +
dimensions:
 +
    time = UNLIMITED ; // (2 currently)
 +
    bnds = 2 ;
 +
variables:
 +
    double time(time) ;
 +
        time:bounds = "time_bnds" ;
 +
        time:units = "days since 1-1-2013 " ; // This date correspond to the example basetime
 +
        time:calendar = "360_day" ; // Other calendars can be used... change here to relevant calendar
 +
        time:axis = "T" ;
 +
        time:long_name = "time" ;
 +
        time:standard_name = "time" ;
 +
    double time_bnds(time, bnds) ;
 +
data:
 +
  time = 180, 540 ; // If you use a different calendar these values will change.
 +
                              //This is the middle of the time_bnds
 +
  time_bnds =
 +
  0, 360, //If you use a different calendar these values will change.
 +
              //These are the day since basetime at the beginning and end of the year
 +
  360, 720 ;
 +
 
 +
and the projection (note that the full time entries are not shown, only beginning and end):
 +
dimensions:
 +
    time = UNLIMITED ; //
 +
    bnds = 2 ;
 +
variables:
 +
    double time(time) ;
 +
        time:bounds = "time_bnds" ;
 +
        time:units = "days since 1-1-2013 " ; // This date correspond to the example basetime
 +
        time:calendar = "360_day" ;  // Other calendars can be used... change here to relevant calendar
 +
        time:axis = "T" ;
 +
        time:long_name = "time" ;
 +
        time:standard_name = "time" ;
 +
    double time_bnds(time, bnds) ;
 +
variables:
 +
    double time(time) ;
 +
        time:bounds = "time_bnds" ;
 +
data:
 +
  time = 900, 1260, 1620, ..., 31140, 31500; // If you use a different calendar these values will change.
 +
                                                      //This is the middle of the time_bnds
 +
  time_bnds =
 +
720, 1080, //If you use a different calendar these values will change.
 +
                  //These are the day since basetime at the beginning and end of the year
 +
1080, 1440,
 +
1440, 1800,
 +
....
 +
30960, 31320,
 +
31320, 31680;
 +
 
 +
====A2.3.3 Table A1: Variable request for ISMIP6====
 +
 
 +
If your quantity does not change with time, then simply save one time entry. An example is geothermal heat flux, which varies in some models but not others.
  
 
{| class="wikitable"
 
{| class="wikitable"
 +
! colspan="7"| Table A1: Variable request for ISMIP6 projections.
 +
 +
Bold names or "alias" indicate a change compared to initMIP, to align the request with the CMIP6 official MIPtable "IyrGre" or names in the CF convention. If possible please use the new names, and if not, the name change will occur when your files are checked for CMIP compliance.
 +
 +
The first entry should be that from which the simulation starts. Fields such as surface mass balance flux should be what was applied as boundary conditions.
 +
 +
|-
 
| align="center" style="background:#f0f0f0;"|'''Variable'''
 
| align="center" style="background:#f0f0f0;"|'''Variable'''
 
| align="center" style="background:#f0f0f0;"|'''Dim'''
 
| align="center" style="background:#f0f0f0;"|'''Dim'''
Line 389: Line 643:
 
|- align="center"
 
|- align="center"
  
|colspan="7"| 2D variables requested every five years, starting at t=0, snapshots for type ST and as five year average for type FL
+
|colspan="7"|''' 2D variables requested yearly as snapshots (end of the year) for type ST and as yearly average for type FL.'''
  
 
|- align="center"
 
|- align="center"
Line 419: Line 673:
 
|| bedrock_altitude
 
|| bedrock_altitude
 
|| m
 
|| m
|| The bedrock topography (unchanged in forward exps.)
+
|| The bedrock topography (may change during the projections)
  
 
|- align="center"
 
|- align="center"
  
 
|| Geothermal heat flux
 
|| Geothermal heat flux
|| x,y
+
|| x,y,t
||C
+
||FL
 
|| hfgeoubed
 
|| hfgeoubed
|| upward_geothermal_heat_flux_at_ground_level  
+
|| '''upward_geothermal_heat_flux_in_land_ice''' alias "upward_geothermal_heat_flux_at_ground_level"
 
|| W m-2
 
|| W m-2
|| Geothermal Heat flux (unchanged in forward exps.)
+
|| Geothermal Heat flux at the ice interface (only needed beneath the grounded ice). If this quantity does not change with time simply enter one timestep
  
 
|- align="center"
 
|- align="center"
Line 443: Line 697:
 
|- align="center"
 
|- align="center"
  
|| Basal mass balance flux
+
|| Basal mass balance flux beneath grounded ice
 
|| x,y,t
 
|| x,y,t
 
||FL
 
||FL
|| libmassbf
+
|| '''libmassbfgr''' alias "libmassbf"
 
|| land_ice_basal_specific_mass_balance_flux
 
|| land_ice_basal_specific_mass_balance_flux
 
|| kg m-2 s-1
 
|| kg m-2 s-1
|| Basal mass balance flux (for areas covered by ice only)
+
|| Basal mass balance flux (only beneath grounded ice)
 +
 
 +
|- align="center"
 +
 
 +
|| Basal mass balance flux beneath floating ice
 +
|| x,y,t
 +
||FL
 +
|| '''libmassbffl''' alias "libmassbf"
 +
|| land_ice_basal_specific_mass_balance_flux
 +
|| kg m-2 s-1
 +
|| Basal mass balance flux (only beneath floating ice)
 +
 
  
 
|- align="center"
 
|- align="center"
Line 466: Line 731:
 
|| x,y,t
 
|| x,y,t
 
||ST
 
||ST
|| uvelsurf
+
|| '''xvelsurf''' alias "uvelsurf"
 
|| land_ice_surface_x_velocity
 
|| land_ice_surface_x_velocity
 
|| m s-1
 
|| m s-1
Line 476: Line 741:
 
|| x,y,t
 
|| x,y,t
 
||ST
 
||ST
|| vvelsurf
+
|| '''yvelsurf''' alias "vvelsurf"
 
|| land_ice_surface_y_velocity
 
|| land_ice_surface_y_velocity
 
|| m s-1
 
|| m s-1
Line 486: Line 751:
 
|| x,y,t
 
|| x,y,t
 
||ST
 
||ST
|| wvelsurf
+
|| '''zvelsurf''' alias "wvelsurf"
 
|| land_ice_surface_upward_velocity  
 
|| land_ice_surface_upward_velocity  
 
|| m s-1
 
|| m s-1
Line 496: Line 761:
 
|| x,y,t
 
|| x,y,t
 
||ST
 
||ST
|| uvelbase
+
|| '''xvelbase''' alias "uvelbase"
 
|| land_ice_basal_x_velocity
 
|| land_ice_basal_x_velocity
 
|| m s-1
 
|| m s-1
Line 506: Line 771:
 
|| x,y,t
 
|| x,y,t
 
||ST
 
||ST
|| vvelbase
+
|| '''yvelbase''' alias "vvelbase"
 
|| land_ice_basal_y_velocity
 
|| land_ice_basal_y_velocity
 
|| m s-1
 
|| m s-1
Line 516: Line 781:
 
|| x,y,t
 
|| x,y,t
 
||ST
 
||ST
|| wvelbase
+
|| '''zvelbase''' alias "wvelbase"
 
|| land_ice_basal_upward_velocity  
 
|| land_ice_basal_upward_velocity  
 
|| m s-1
 
|| m s-1
Line 526: Line 791:
 
|| x,y,t
 
|| x,y,t
 
|| ST
 
|| ST
|| uvelmean
+
|| '''xvelmean''' alias "uvelmean"
 
|| land_ice_vertical_mean_x_velocity
 
|| land_ice_vertical_mean_x_velocity
 
|| m s-1
 
|| m s-1
Line 536: Line 801:
 
|| x,y,t
 
|| x,y,t
 
|| ST
 
|| ST
|| vvelmean
+
|| '''yvelmean''' alias "vvelmean"
 
|| land_ice_vertical_mean_y_velocity
 
|| land_ice_vertical_mean_y_velocity
 
|| m s-1
 
|| m s-1
Line 546: Line 811:
 
|| x,y,t
 
|| x,y,t
 
||ST
 
||ST
|| litempsnic
+
|| '''litemptop''' alias "litempsnic"
|| temperature_at_ground_level_in_snow_or_firn  
+
|| '''temperature_at_top_of_ice_sheet_model''' alias "temperature_at_ground_level_in_snow_or_firn"
 
|| K
 
|| K
 
|| Ice temperature at surface
 
|| Ice temperature at surface
Line 553: Line 818:
 
|- align="center"
 
|- align="center"
  
|| Basal temperature
+
|| Basal temperature beneath grounded ice sheet
 
|| x,y,t
 
|| x,y,t
 
||ST
 
||ST
|| litempbot
+
|| '''litempbotgr''' alias "litempbot"
|| land_ice_basal_temperature
+
|| '''temperature_at_base_of_ice_sheet_model''' alias "land_ice_basal_temperature"
 
|| K
 
|| K
|| Ice temperature at base
+
|| Ice temperature at base of grounded ice sheet
 +
 
 +
|- align="center"
 +
 
 +
|| Basal temperature beneath floating ice shelf
 +
|| x,y,t
 +
||ST
 +
|| '''litempbotfl''' alias "litempbot"
 +
|| '''temperature_at_base_of_ice_sheet_model''' alias "land_ice_basal_temperature"
 +
|| K
 +
|| Ice temperature at base of floating ice shelf
  
 
|- align="center"
 
|- align="center"
Line 567: Line 842:
 
||ST
 
||ST
 
|| strbasemag
 
|| strbasemag
|| magnitude_of_land_ice_basal_drag
+
|| '''land_ice_basal_drag''' alias "magnitude_of_land_ice_basal_drag"
 
|| Pa
 
|| Pa
|| Magnitude of basal drag
+
|| Basal drag
  
 
|- align="center"
 
|- align="center"
Line 581: Line 856:
 
|| Loss of ice mass resulting from iceberg calving. Only for grid cells in contact with ocean
 
|| Loss of ice mass resulting from iceberg calving. Only for grid cells in contact with ocean
  
 +
|- align="center"
 +
 +
|| '''Ice front calving and melt flux'''
 +
|| x,y,t
 +
||FL
 +
|| '''lifmassbf'''
 +
|| '''land_ice_specific_mass_flux_due_to_calving_and_ice_front_melting'''
 +
|| kg m-2 s-1
 +
|| Loss of ice mass resulting from calving and ice front melting. Only for grid cells in contact with ocean
 
|- align="center"
 
|- align="center"
  
Line 607: Line 891:
 
||ST
 
||ST
 
|| sftflf
 
|| sftflf
|| floating_ice_sheet_area_fraction
+
|| '''floating_ice_shelf_area_fraction''' alias "floating_ice_sheet_area_fraction"
 
|| 1
 
|| 1
 
|| Fraction of grid cell covered by ice sheet flowing over seawater
 
|| Fraction of grid cell covered by ice sheet flowing over seawater
Line 613: Line 897:
 
|- align="center"
 
|- align="center"
  
| colspan="7"| Scalar outputs requested every full year, as snapshots for type ST as 1 year averages for type FL. The t=0 value should contain the data of the initialization.
+
| colspan="7"| '''Scalar outputs requested every full year: snapshots for type ST and 1 year averages for type FL. '''
  
 
|- align="center"
 
|- align="center"
Line 640: Line 924:
 
|| t
 
|| t
 
||ST
 
||ST
|| iareag
+
|| '''iareagr''' alias "iareag"
|| grounded_land_ice_area
+
|| grounded_ice_sheet_area
 
|| m^2
 
|| m^2
 
|| spatial integration
 
|| spatial integration
Line 650: Line 934:
 
|| t
 
|| t
 
||ST
 
||ST
|| iareaf
+
|| '''iareafl''' alias "iareaf"
 
|| floating_ice_shelf_area  
 
|| floating_ice_shelf_area  
 
|| m^2
 
|| m^2
Line 674: Line 958:
 
|| kg s-1
 
|| kg s-1
 
|| spatial integration
 
|| spatial integration
 +
 +
|- align="center"
 +
 +
|| '''Total BMB flux beneath floating ice'''
 +
|| t
 +
||FL
 +
|| '''tendlibmassbffl'''
 +
|| tendency_of_land_ice_mass_due_to_basal_mass_balance
 +
|| kg s-1
 +
|| '''spatial integration (computed beneath floating ice only)'''
  
 
|- align="center"
 
|- align="center"
Line 682: Line 976:
 
|| tendlicalvf
 
|| tendlicalvf
 
|| tendency_of_land_ice_mass_due_to_calving
 
|| tendency_of_land_ice_mass_due_to_calving
 +
|| kg s-1
 +
|| spatial integration
 +
 +
|- align="center"
 +
 +
|| '''Total calving and ice front melting flux'''
 +
|| t
 +
||FL
 +
|| '''tendlifmassbf'''
 +
|| '''tendency_of_land_ice_mass_due_to_calving_and_ice_front_melting'''
 
|| kg s-1
 
|| kg s-1
 
|| spatial integration
 
|| spatial integration
Line 696: Line 1,000:
 
| align="center" style="background:#f0f0f0;"|'''Group'''
 
| align="center" style="background:#f0f0f0;"|'''Group'''
  
|}
+
|- align="center"
 +
||[https://www.awi.de/ueber-uns/organisation/mitarbeiter/martin-rueckamp.html Martin Rückamp], [https://www.awi.de/ueber-uns/organisation/mitarbeiter/angelika-humbert.html Angelika Humbert] 
 +
||ISSM
 +
||AWI
 +
||Alfred Wegener Institute for Polar and Marine Research, DE /University of Bremen, DE
  
===Model Characteristics===
+
|- align="center"
 +
||[http://www.bristol.ac.uk/geography/people/victoria-lee/index.html Victoria Lee], [http://www.bristol.ac.uk/geography/people/tony-j-payne/index.html Tony Payne], [https://www.swansea.ac.uk/staff/science/geography/cornfordsl/ Stephen Cornford], [https://crd.lbl.gov/departments/computational-science/anag/about/staff-and-postdocs/martin/ Daniel Martin]
 +
||BISICLES
 +
||BGC
 +
||University of Bristol, Bristol, UK, Department of Geography, Swansea University, UK, Computational Research Division, Lawrence Berkeley National Laboratory, California, USA
  
{| class="wikitable" style="margin: 1em auto 1em auto;"
+
|- align="center"
| align="center" style="background:#f0f0f0;"|'''Model ID'''
+
||[https://science.gsfc.nasa.gov/sed/bio/isabel.j.nias Isabel Nias], [https://science.gsfc.nasa.gov/sed/bio/sophie.nowicki Sophie Nowicki], [https://science.gsfc.nasa.gov/sed/bio/denis.felikson Denis Felikson]
| align="center" style="background:#f0f0f0;"|'''Numerics'''
+
||ISSM
| align="center" style="background:#f0f0f0;"|'''Ice Flow'''
+
||GSFC
| align="center" style="background:#f0f0f0;"|'''Initialization'''
+
||NASA Goddard Space Flight Center, Greenbelt, USA
| align="center" style="background:#f0f0f0;"|'''Initial Year'''
+
| align="center" style="background:#f0f0f0;"|'''Initial SMB'''
+
| align="center" style="background:#f0f0f0;"|'''Velocity'''
+
| align="center" style="background:#f0f0f0;"|'''Bed'''
+
| align="center" style="background:#f0f0f0;"|'''Surface'''
+
| align="center" style="background:#f0f0f0;"|'''GHF'''
+
| align="center" style="background:#f0f0f0;"|'''Res min'''
+
| align="center" style="background:#f0f0f0;"|'''Res max'''
+
  
 
|- align="center"
 
|- align="center"
||
+
||[https://urldefense.proofpoint.com/v2/url?u=http-3A__wwwice.lowtem.hokudai.ac.jp_-7Egreve_&d=DwIDaQ&c=ApwzowJNAKKw3xye91w7BE1XMRKi2LN9kiMk5Csz9Zk&r=DYrfyWvYVW12-MpHGIDM1nFcJ6mmhPanJO8ZdvyeEtk&m=CLB5CLoK8qf5ACikPd_QyvBhIDiveMr4ixXFt4-me1g&s=fXk7_nxexbdnG71yPry2rXEGh8F1DJF-6xAUkD4pntI&e= Ralf Greve], [https://urldefense.proofpoint.com/v2/url?u=https-3A__www.pik-2Dpotsdam.de_members_calov_&d=DwIDaQ&c=ApwzowJNAKKw3xye91w7BE1XMRKi2LN9kiMk5Csz9Zk&r=DYrfyWvYVW12-MpHGIDM1nFcJ6mmhPanJO8ZdvyeEtk&m=CLB5CLoK8qf5ACikPd_QyvBhIDiveMr4ixXFt4-me1g&s=H_M8eXkHJmWhn73Dxg28K2ZF0-QbOjhBpZJxGDgHzoQ&e= Reinhard Calov], [https://urldefense.proofpoint.com/v2/url?u=https-3A__twitter.com_YouStormorg&d=DwICaQ&c=ApwzowJNAKKw3xye91w7BE1XMRKi2LN9kiMk5Csz9Zk&r=DYrfyWvYVW12-MpHGIDM1nFcJ6mmhPanJO8ZdvyeEtk&m=LXh_UBJnfYZGMro-mksV9qqP8lpABelH3XiPJF37n4o&s=yIqXja01f7Rywj4jzx_MZ-cnxlWsFL1zOErYbSsiIFk&e=    Chris Chambers]
||
+
||SICOPOLIS
||
+
||ILTS_PIK
||
+
||Institute of Low Temperature Science, Hokkaido University, Sapporo, JP, Potsdam Institute for Climate Impact Research, Potsdam, DE
||
+
 
||
+
|- align="center"
||
+
||[https://www.uu.nl/staff/HGoelzer/0 Heiko Goelzer], [https://www.uu.nl/staff/RSWvandeWal Roderik van de Wal], [https://www.staff.science.uu.nl/~broek112/ Michiel van den Broeke]
||
+
||IMAUICE
||
+
||IMAU
||
+
||Utrecht University, Institute for Marine and Atmospheric Research (IMAU), Utrecht, NL
||
+
 
||
+
|- align="center"
 +
||[https://science.jpl.nasa.gov/people/HSeroussi/ Helene Seroussi], [https://science.jpl.nasa.gov/people/NSchlegel/ Nicole Schlegel]
 +
||ISSM
 +
||JPL
 +
||NASA Jet Propulsion Laboratory, Pasadena, USA
 +
 
 +
|- align="center"
 +
||[https://science.jpl.nasa.gov/people/Cuzzone/ Joshua K. Cuzzone], [https://science.jpl.nasa.gov/people/NSchlegel/ Nicole Schlegel]
 +
||ISSMPALEO
 +
||JPL
 +
||NASA Jet Propulsion Laboratory, Pasadena, USA
 +
 
 +
|- align="center"
 +
||[https://sites.google.com/site/aurelienquiquet/ Aurélien Quiquet], [https://www.lsce.ipsl.fr/Phocea/Membres/Annuaire/index.php?uid=dumas Christophe Dumas]
 +
||GRISLI
 +
||LSCE
 +
||LSCE/IPSL, Laboratoire des Sciences du Climat et de l’Environnement, CEA-CNRS-UVSQ, Gif-sur-Yvette, FR
 +
 
 +
|- align="center"
 +
||[https://www.physics.mun.ca/~lev/ Lev Tarasov]
 +
||GSM
 +
||MUN
 +
||Dept of Physics and Physical Oceanography, Memorial University of Newfoundland, USA
 +
 
 +
|- align="center"
 +
||[https://staff.ucar.edu/users/lipscomb William Lipscomb], [https://staff.ucar.edu/users/gunterl Gunter Leguy]
 +
||CISM
 +
||NCAR
 +
||National Center for Atmospheric Research, Boulder, CO, USA
 +
 
 +
|- align="center"
 +
||[http://glaciers.gi.alaska.edu/people/aschwanden Andy Aschwanden]
 +
||PISM
 +
||UAF
 +
||Geophysical Institute, University of Alaska Fairbanks, USA
 +
 
 +
|- align="center"
 +
||[https://www.ess.uci.edu/node/11148 Youngmin Choi], [https://science.jpl.nasa.gov/people/HSeroussi/ Helene Seroussi], [https://www.ess.uci.edu/people/mmorligh Mathieu Morlighem]
 +
||ISSM
 +
||UCIJPL
 +
||NASA Jet Propulsion Laboratory, Pasadena, USA / University of California Irvine, Irvine, USA
 +
 
 +
|- align="center"
 +
||[http://we.vub.ac.be/en/sébastien-le-clech Sébastien Le Clec’h] and [https://we.vub.ac.be/en/philippe-huybrechts Philippe Huybrechts]
 +
||GISM
 +
||VUB
 +
||Vrije Universiteit Brussel, Brussels, BE
 +
 
 +
|- align="center"
 +
||[https://www.victoria.ac.nz/antarctic/about/staff/nicholas-golledge Nick Golledge], [https://www.victoria.ac.nz/antarctic/about/grad-students/dan-lowry Dan Lowry]
 +
||PISM
 +
||VUW
 +
||GNS Science, Lower Hutt, NZ / Antarctic Research Centre, Victoria University of Wellington, NZ
  
 
|}
 
|}
 +
 +
===Model Characteristics===
 +
 +
The Model Characteristics table can be found [https://docs.google.com/spreadsheets/d/1YEiPV3Uc0K8EqD57mz-ZqBIF4XnVb5UbyRDRmox0KeE/edit#gid=700377079 here].

Latest revision as of 10:37, 4 February 2020

Overview

NOTE: THIS PAGE IS UNDER CONSTRUCTION. WHEN IT IS FINISHED, THIS STATEMENT WILL BE REMOVED. LAST UPDATE: 19 September 2019

This page describes the experimental protocol for the ISMIP6 projections that target the upcoming IPCC AR6 assessment. Due to the delay in CMIP6 climate simulations, the initial set of ISMIP6 simulations are based on CMIP5 projections. As CMIP6 model output become available, ISMIP6 will include simulations based on these models.

The experimental framework was revised in September 2018 during an ISMIP6 workshop held in Sassenheim (NL). The protocol, summarized in Fig 1, allows for:

• Sampling CMIP scenarios: main focus is on the high emission RCP8.5, but ice sheet evolution in response to low emission RCP2.6 is also investigated.

• Sampling CMIP models: 6 AOGCMs have been selected from the CMIP5 model ensemble. The AOGCMs were identified based on the following steps: 1) present plausible climates near Greenland (evaluated by model biases over the historical period), 2) have the data at the temporal resolution needed for RCM downscaling, 3) sample a diversity of forcing (evaluated by differences in projections and code similarities) and 4) allow only models with RCP8.5 and ideally with RCP2.6. The sampling methodology for the CMIP5 models is described in Barthel et al. (under review). As CMIP6 models became available, ISMIP6 prepared dataset for CNRM-CM6 (ssp585 and ssp126), CESM2 (ssp585, in prep) and UKEM1-CM6 (ssp585). Unlike for the rigorous analysis for the CMIP5 models, the CMIP6 models were selected because of their availability.

• Sampling ice sheet model uncertainty: "standard" and "open" experiments. The "standard" experiments are based on parameterizations developed by the ocean focus group, while "open" experiments utilize the parameterizations already in use by respective ice sheet models. The open experiments are an important contribution as the use of different parameterizations reflect the uncertainty in our current understanding of ice-ocean interactions.

• Sampling ocean forcing uncertainty: the standard experiments include "high", "mid" and "low" parameters.

• Experiment ranking: This experimental framework results in a series of projections (divided into core and targeted experiments), a historical run and a control run. Not every ice sheet model will be able to carry out the full set of experiments, but they are strongly encouraged to participate in the full suite of core experiments (Table 1, which uses both the "open" and "standard" experiments/parameterizations). Given that the standard experiments requires new implementations, it is OK to only participate with only the open experiments. Similarly, it is OK to only participate with the "standard" experiments. Groups are encouraged to work through the lists presented below (see Table 1 for core experiments, and targeted experiments will follow soon), starting from the top, and complete as many experiments as possible. This approach of combining core and targeted experiments is based on Shannon et al. (2013): it ensures that all groups do a subset of identical experiments, while it also allows faster models to explore the targeted experiment space more fully.

Greenland exp design.png

Figure 1: Overview of the Greenland experimental framework

List of Projections

With the help of the atmosphere and ocean focus groups, a number of CMIP5 AOGCMs have been selected for ISMIP6 standalone ice sheet model projections. Table 1 lists the core experiments which are the minimum contribution expected from ISMIP6 models, in addition to the initialization experiments listed in Table 2. All groups are encouraged to contribute to the "standard" experiments. Groups that have their own methods for implementing ocean and atmosphere forcing, are encouraged to do the suite with "open" experiments (1-4), but these are not compulsory. Models that perform the "open" experiments can use the parameterization of their choice to simulate atmospheric and oceanic forcings, but these parameterizations must use the given CMIP5 AOGCM outputs.

Modeling groups that can run many simulations will be encouraged to further explore the ice sheet response using targeted experiments (See Table). These include three additional CMIP5 AOGCMs under RCP8.5, and experiments that explore the ocean forcing uncertainty. Depending on the results of experiments 3 and 7, which consider RCP2.6, additional AOGCMs may be suggested with RCP2.6 for models that are able to do many simulations, but these would be a lower priority than the completing the set of experiments with the 6 AOGCMs for the RCP8.5 scenario. As CMIP6 AOGCMs are becoming available, we are preparing these datasets. The spreadsheet will be updated as new dataset become available. Because there is value in both completing the 6 CMIP5 AOGCMs (to sample the uncertainty in CMIP) and simulations with CMIP6 models, we encourage groups do to as many experiments as possible.

Note: As of August 6, all the datasets needed for the experiments in Table 1, as well as the CMIP5 and CMIP6 AOGCMs listed in the targeted experiments (See Table) are ready.


Table 1: Core Experiments based on MIROC5, NorESM1-M and HadGEM2-ES
Exp RCP AOGCM Std/open Ocean Forcing Unc Note
1 8.5 MIROC5 Open Medium Expected largest response to SMB, median ocean warming
2 8.5 NorESM Open Medium Low atmosphere change, low ocean warming
3 2.6 MIROC5 Open Medium Expected largest response to SMB, median ocean warming
4 8.5 HadGEM2-ES Open Medium Expected median response to SMB, median ocean warming
5 8.5 MIROC5 Standard Medium Expected largest response to SMB, median ocean warming
6 8.5 NorESM Standard Medium Low atmosphere changes, low ocean warming
7 2.6 MIROC5 Standard Medium Expected largest response to SMB, median ocean warming
8 8.5 HadGEM2-ES Standard Medium Expected median response to SMB, median ocean warming
9 8.5 MIROC5 Standard High Ocean Forcing Uncertainty
10 8.5 MIROC5 Standard Low Ocean Forcing Uncertainty

Initial state, control run, historical run and projections set up

The core and targeted experiments all start on January 2015 and end in December 2100. The start date follows the CMIP6 protocol for projections, while the end date is constrained by the availability of forcing. In many cases, modelers will need to run a short historical run to bring their models from the “initialization date” to the “projection start date” of January 2015.

The “initialization date” (or initial state) is left to the modeling groups discretion and can be any time prior to January 2015. The “initialization date” corresponds to the date assigned to the initialization procedure. Groups can reuse their initMIP initialization configuration or generate a new initial state. In the later case, it is important to redo the initMIP schematic experiments ('asmb', Goelzer et al. 2018), as it will help understanding how a novel initial state contributes to the uncertainty in ice sheet evolution.

Note that ISMIP6 has a new standard grid (horizontal resolution of 1 km on projection EPSG:3413), which is different than that used in initMIP-Greenland. The revised datasets for the initMIP 'asmb' forcing for the new grid can be obtained from: sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6/Projections/GrIS/initMIP/dSMB_epsg3413 We have changed to the new grid because its projection is standardised and used by many observational datasets that are key input for ice sheet models.


A control run (‘ctrl’) is also needed to evaluate model drift. As for initMIP, the control run is obtained by running the model forward, keeping the surface mass balance and ocean forcing used in the initialization technique unchanged. The control run starts from the initial state (typically before 2015) and should last a minimum of 100 years, the same duration as the schematic initMIP experiments. The control run should also be sufficiently long to reach 2100. (See examples in Table 2).

Note that in the event that initMIP schematic experiments and control run are redone as part of the projection setup, then consistency with the projections protocol is more important than consistency with the original initMIP setup. For example, in initMIP bedrock was not allowed to evolve. However, an ISM planning to run the projection with evolving bedrock AND planning to redo the initial state, would also rerun initMIP (ctrl, asmb) with bedrock change. Similarly, the original initMIP requested 2D output every 5 years, whereas the projections protocol request for yearly values. Therefore if rerunning initMIP, the schematic experiments should be saved yearly.

A single "historical run" is required for each ice sheet model, from which all the projections will branch off. The historical run starts at the initMIP state and ends in December 2014. Groups are free to choose how to run the “historical run”, using a reanalysis, a historical run from an RCM, a historical run from an AOGCM or combination of multiple datasets. However, we reject multiple historical runs for each individual AOGCM, because this would very much complicate the forcing strategy and interpretation. As a consequence, in case AOGCM data is used, please decide for one. ISMIP6 provides a climatology for the SMB and surface temperature for each of the AOGCMs used to generate the projection dataset, as well as anomalies. Groups that would like to use an ISMIP6 Greenland dataset for initialization and historical run are recommended to use MIROC5. Note that the SMB and temperature climatology corresponds to January 1960-December 1989 because the ice sheet is assumed in steady state with the climate during this period. To give a concrete example of a historical run that seems to work well: The ice sheet is initialised to a steady state in 1960. RCM anomalies (here RACMO) are applied from 1960 - 2014 on top of the SMB assumed during initialisation (here RACMO average 1960-1990). Other SMB products forced with reanalysis data (e.g. MAR) are expected to work equally well.

The "projection control" (ctrl_proj) is an unforced simulation that starts at the same ice sheet state as the projections. It will run until 2100, and is implemented with zero anomalies. It is meant to capture how much drift arises from the historical run.

In most cases the same historical run (and potentially spin-up) can be used for both the standard and open experiments. However, if the implementation of the open and standard experiments requires changes to the ISM that are substantially different (in terms of physics in the ISM), then modelers are allowed to carry out an historical run for the open experiments and an historical run for the standard experiments.

Table 2: Initialization experiments and examples of different initialization start date
Experiment Note start 1 (duration) start 2 (duration) start 3 (duration)
ctrl Unforced control run, needed for model drift evaluation 2015 (100 years) 2005 (100 years) 1980 (120 years)
asmb* initMIP prescribed surface mass balance anomaly 2015 (100 years) 2005 (100 years) 1980 (100 years)
historical needed to bring model from initial state to projection start date N/A (0 years) 2005 (10 years) 1980 (35 years)
ctrl_proj Unforced control run, starting from January 2015, needed for model drift evaluation following historical 2015 (85 years) 2015 (85 years) 2015 (85 years)
*only needed if initial state is different from initMIP

Atmospheric forcing: SMB and temperature anomalies

Introduction

ISMIP6 provides anomalies of SMB and temperature, along with associated climatology, as it allows for an experimental framework that is applicable to a diverse set of ice sheet models.

Before applying SMB anomalies, ISMs need to be initialized by applying a baseline SMB (either a time series or a climatology). ISMIP6 provides SMB climatologies for the reference period (January 1960 to December 1989) from the same models computing the anomalies. This reference period for Greenland SMB was chosen because the ice sheet is assumed to be in steady state with the surrounding climate over that period. ISMs can use these climatologies for spin-up, if desired, but are free to use their own preferred SMB forcing.

Elevation feedbacks have been shown to be important in century-scale simulations. For example, Le clec’h et al. (2019) considered differences between no coupling (SMB independent of z), one-way coupling (i.e., correcting SMB outputs for ISM topography changes), and two-way coupling (allowing ice-sheet topographic changes to feed back on the RCM simulation.) Their results suggest that one-way coupling is sufficient to represent elevation feedbacks until the end of the 21st century. For large topographic changes on longer timescales, it might be necessary to incorporate two-way feedbacks, which are beyond the scope of the standalone ISM effort (but will be considered in the ISMIP6 coupled GCM-ISM experiments).

Provided forcing data sets

ISMIP6 provides surface forcing datasets for the Greenland ice sheet (GrIS) based on CMIP AOGCM simulations. The AOGCM output is re-interpreted through higher-resolution regional climate models (RCMs). The later allows to capture narrow regions at the periphery of the Greenland ice sheet with large surface mass balance (SMB) gradients, which are not captured by CMIP5 AOGCMs. For CMIP6, many of the AOGCMs that have indicated participation in ISMIP6 now use multiple elevation classes to downscale SMB to finer grid resolution. Once these models have completed the CMIP6 projections, our goal is to include additional ISMIP6 projections using SMB downscaled via elevation class.

For the ISMIP6 projections based on CMIP5 AOGCMs, the surface forcing datasets were prepared by Xavier Fettweis, using the MAR regional climate model. RCM downscaling can take account of modest future topography/extent changes, and thus cope with the fact that individual ISM runs may not use exactly the same geometry. All RCM runs use a fixed topography, but vertical SMB gradients for each grid cell are derived using the method described by Franco et al. (2012). At each location the vertical gradient in SMB is found by summing and averaging pairwise differences in the nearest neighbor cells. This vertical gradient is used to downscale the SMB (15 km) to a finer grid (1 km), allowing resolution of steep topography that is not represented accurately on the coarse grid. The same information is used to parameterise the SMB-height feedback in the projections. In addition, MAR calculates a potential SMB term for areas that are outside of the observed ice sheet extent, allowing application to ISMs with ice lying outside the MAR ice-sheet mask. However, experience with initMIP has shown that large variations in ice-sheet extent can lead to considerable bias in the projections. We therefore propose a method (see below) for models with large difference from the observed ice sheet extent, to remap the SMB anomaly to the individual modelled ice sheet geometry.

The atmospheric forcings consist of annual anomalies of SMB (aSMB) and surface temperature, along with a fields to represent the dependence of SMB and surface temperature on elevation (dSMBdz, dTdz). SMB is needed by ISMs to compute mass changes at the surface, and surface temperature (i.e., the ice temperature at the base of the snow, as distinct from the 2-m air temperature or skin temperature) is used by many ISMs as an upper boundary condition in the ice temperature calculations. ISMIP6 also provides the climatologies that were used to calculate the anomalies: aSMB_AOGCM(x,y,t) = SMB_AOGCM(x,y,t) - SMB_CLIM_AOGCM(x,y), where SMB_AOGCM is the downscaled SMB for a given AOGCM (using MAR RCM) and SMB_CLIM_AOGCM is the corresponding climatology. The SMB_CLIM_AOGCM were computed by taking the mean value of all SMB_AOGCM over the reference period (from January 1960 to December 1989). ISMs can use these climatologies for spin-up, if desired, but are free to use their own preferred SMB forcing for spin-up and historical run.

The SMB anomaly aSMB is given in units [kg m-2 s-1] in yearly values, one year per file. It should be applied constant over a full year and step change at the beginning of a new year. To convert to units [m yr-1] typically used in an ice sheet model, multiply the netcdf variable by 31556926 s/yr, 1/1000 m3/kg and by the density ratio rhow/rhoi:

aSMB [m yr-1] = aSMB [kg m-2 s-1] * 31556926 / 1000 * (1000/rhoi), where rhoi is your specific ice density (typically 917.0 or similar).

The SMB climatology and its anomalies are provided on the ISMIP6 1 km standard ice sheet grid for Greenland. ISMs then horizontally interpolate the anomaly forcing conservatively from the standard grid to their native grids.

The SMB change with surface elevation dSMBdz is given in units [kg m-2 s-1 m-1] in yearly values, one year per file. To parameterise the SMB-height feedback, the SMB has to be corrected by dSMBdz * h-h_ref, updated every full year, where h is the time evolving modelled ice sheet surface elevation and h_ref is the modelled surface elevation in December 2014.

Implementation

Let SMB_ref(x,y) denote the SMB used to initialize the ISM, and let h_ref(x,y) denote the ice sheet surface elevation at the end of the initialization. If a time-dependent SMB is used for spin-up, then SMB_ref(x,y) is the average over the SMB reference period (from January 1960 to December 1989 for Greenland). If an SMB climatology is used in the assimilation or spin-up, then SMB_ref(x,y) is simply the climatology. ISMIP6 recognizes that the use of existing climatologies (or dataset of SMB averaged over many year) used in the initialization may not align with the time period for the SMB reference period. However, it is assumed that the differences between climatologies will be less than the inter annual variability from the SMB resulting from the AOGCMs, and thus changes in aSMB. What is important is that the SMB climatology (or SMB_ref) is computed over many years. This assumption also allows for ISM to work with their favorite SMB.

Here we propose two different methods for implementation of atmospheric forcing, depending on how close the ice sheet mask at the end of the initialization is to the observed ice sheet mask, the one assumed by the RCM.

Method 1: when the ice sheet mask is similar to the observed

This is typically the case for ice sheet models that use data assimilation in their initialization, but could be the case for other modelling approaches. For those models, ISMIP6 provides aSMB(x,y,t) at h_rcm(x,y), along with dSMBdz(x,y,t). Here, aSMB is the time-dependent SMB anomaly in a changing climate, computed in an RCM with fixed surface topography h_rcm, and dSMBdz(x,y,t) is the time-dependent vertical gradient of SMB. aSMB and dSMBdz are provided on an annual basis and should be updated every full year.

Given aSMB(x,y,t) and dSMBdz(x,y,t) on the standard grid, the ISM horizontally interpolate these fields to its local grid. Then during runtime, the SMB at a given time and location is computed as

SMB(x,y,t) = SMB_ref(x,y) + aSMB(x,y,t) + dSMBdz(x,y) * [h(x,y,t) - h_ref(x,y)],

where h(x,y,t) is the time-dependent surface elevation. ISMs will likely need to implement code changes to handle the lapse-rate correction. The models do not need h_rcm(x,y) to compute SMB, but it is provided for reference. Since dSMBdz is computed at h_rcm and is not given as a function of z, this approach may be inaccurate if h_ref is significantly different from h_rcm.

The datasets of SMB (aSMB, dSMBdz), and surface temperature (aST, dSTdz) from 1950 to 2100, along with the climatology (1960-1989) for MIROC5-rcp8.5, MIROC5-rcp2.6, and NorESM1-rcp8.5 can be obtained via the ISMIP6 ftp server (email ismip6@gmail.com to obtain the login information) at

sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6/Projections/GrIS/Atmosphere_Forcing/aSMB_observed/

Unless advised otherwise, always use the latest version of the forcing files (v1, on Feb 25, 2019)

Method 2: when the ice sheet mask is very different from the observed

This is typically the case for ice sheet models that use a glacial-interglacial spinup in their initialization, but could also be the case e.g. for other models that fully relax to a suboptimal SMB. For those models, ISMIP6 generates a time-dependent SMB anomaly, aSMB(x,y,t) and dSMBdz(x,y,t) that are applied as described above. However, the main difference is that the forcing files are specific for the geometry of your modelled initial state.

In order to make the forcing applicable for different ice sheet geometries, we first translate a given SMB anomaly field as a function of absolute location, to a function of surface elevation for 25 regional drainage basins. This step exploits the strong elevation dependence of aSMB. We can then remap aSMB (and dSMBdz) to different modelled geometries. This preserves the overall aSMB patterns and reduces unphysical biases. The procedure to generate the remapped aSMB and dSMBdz is described in Goelzer et al. (2019, in prep).

To obtain the datasets, modelers simply need to provide ISMIP6 with their modelled initial surface elevation h_ref(x,y) and ice mask sftgif(x,y). ISMIP6 will then compute aSMB for you. Please send an email to ismip6-at-gmail.com when you are ready to upload your initial state (h_ref, sftgif).

An example data set produced for the observed geometry as h_ref(x,y) is available for MAR3.9-MIROC5-RCP85 at

sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6/Projections/GrIS/Atmosphere_Forcing/aSMB_remapped/

The decision on which forcing method to use depends on the expected biases inherent to both approaches. Given the initial state (h_ref and sftgif) we can estimate the biases with a simple integration of the SMB anomaly for a static case (i.e. no ice dynamics). If you are not sure which forcing strategy is best suited for you model, please contact us by email to ismip6-at-gmail.com.

Oceanic forcing: Calving and frontal melt

ISMIP6 provides dataset of runoff and ocean thermal forcing for models that have their own methods for implementing oceanic induced retreat. In addition, modeling groups are expected to participate with one of the two ISMIP6 approaches described below. The ISMIP6 Standard approach is a simple retreat intended to be easily implemented by the majority of ISM taking part in ISMIP6. Alternatively, for models that wish to implement a more complex oceanic forcing, the ISMIP6 Greenland ocean focus group has developed a second methodology, described in ISMIP6 high resolution ocean melt rate approach.

ISMIP6 Standard approach imposes an empirically-derived, sector averaged retreat (Fig 2) as a function of climate forcing. This method was developed for ISMIP6 as a result of the ocean forcing focus group and is described in greater details in Slater et al. (2019a; 2019b), and in the webinar: sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6/Projections/GrIS/Ocean_Forcing/Webinar_2018_11_08/

Oceans overview greenland.png Retreat rate.png

Fig 2: Example of the empirically-derived retreat scenarios for the 7 sectors of the Greenland ice sheet, obtained with MIROC5, RCP8.5.

As described in the webinar, retreat is imposed when the ice sheet geometry and ice front retreat scenario indicates that the land_ice_area_fraction_retreat mask is ice free for a given year. (Note the name was chosen so that it is closely related to the standard name land_ice_area_fraction -corresponding to variable name sftgif- in the ISMIP6 data request. For ease of communication, we use the longer standard name). Implementation for a specific ISM (also illustrated in Fig 3) requires the following steps:

1. Identify ice prone to outlet glacier retreat, by interpolating initial ice mask conservatively to 1 km ISMIP6 diagnostic grid
 to obtain a mask for ice fraction:

land_ice_area_fraction(x,y) = [0.0, …, 1.0] with 0.0 for ice free and 1.0 for ice covered

2. Calculate the distance to nearest ocean grid cell, or "distance_map", based on model ice fraction and observed mask of ice with potential ocean contact -> distance_map(x,y)

3. Calculate land ice fraction retreat based on distance map and retreat scenario: 


land_ice_area_fraction_retreat(x,y,t) = [0.0, …, 1.0] with 0.0 for ice free and 1.0 for ice covered

4. Interpolate land_ice_area_fraction_retreat mask from diagnostic grid to model grid (conservatively)

5. Apply retreat in forward experiments (sub-grid implementation may be required for models with coarse resolutions to allow for partial retreat):

if land_ice_area_fraction_retreat(x,y,t) = 0.0, apply full retreat

if 0.0 < land_ice_area_fraction_retreat(x,y,t) < 1.0 apply partial retreat

Greenland implementation.png

Fig 3: Illustration of the steps required for the implementation of the oceanic retreat.

ISMIP6 will generate the land_ice_area_fraction_retreat masks (steps 2-3) for each model. As retreat is provided as a series of ice fraction masks, ice sheet models with coarse resolution should use a sub-grid approach. A suggestion is to apply an land_ice_area_fraction_retreat that is relative to the reference thickness. Models may have a different strategy for this sub-grid implementation.

Note on the retreat dataset: the retreat rate dataset was calibrated using grounding line position of glaciers that do not have ice shelves. Although the retreat dataset is therefore not optimum for glaciers that have a floating tongue, it is suggested that retreat is imposed at the ice front. Modeling groups that have the capability of computing ice tongue basal melt may use the provided dataset of thermal forcing per basin.

To obtain the land_ice_area_fraction_retreat masks appropriate for your model, please submit your land_ice_area_fraction mask (step 1) to ISMIP6. This should be your ice mask from end 2014. We distinguish different cases to reduce interpolation artefacts. 1) If your native grid is regular and on EPSG:3413, please provide your original modelled ice mask with x,y information. 2) If you have a high resolution irregular grid or a grid on a different projection than EPSG:3413, please get in touch so we can find the best solution for you. When ready, upload your file to the directory /ISMIP6/Projections/GrIS/Ocean_Forcing/Retreat_Implementation/MODELFILES/MODELNAME, where MODELNAME is the name of your model, and let us know via email that your file is uploaded.

ISMIP6 high resolution ocean melt rate allows ice sheet models to specify terminus retreat for each individual marine-terminating outlet glacier, thus allowing glaciers to retreat at different rates. This is accomplished by specifying a glacier terminus melt rate, calculated as a function of subglacial discharge (approximated as surface runoff from each glacier catchment) and ocean thermal forcing (Fig. 4). In addition to the melt rate, a calving rate must be specified to obtain the total frontal ablation rate at each terminus. Additional details can be found in the webinar and the updated presentation slides: sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6/Projections/GrIS/Ocean_Forcing/GrIS_Melt_Parameterization_email_2019Jan.pptx

Dataset for thermal forcing and runoff from 1950 to 2100 are available for MIROC5-rcp2.6, MIROC5-rcp8.5, and NorESM1-rcp8.5 at: sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6/Projections/GrIS/Ocean_Forcing/Melt_Implementation/ As the dataset were created for an ice sheet geometry corresponding to present day observations from GIMP and BedMachine3, ISM that have differ substantially from this initial geometry should not use this dataset. (ie: if you are using Method 2 for the atmospheric forcing/SMB remapping, then the dataset is not applicable for your model, as the observed drainage basins used to generate the runoff field may not work with your model. You could create a dataset that is appropriate for your model from the remapped runoff and a water routing consistent with your geometry).

GrIS Melt Parameterization flowchart.png

Fig 4: High-resolution melt-rate approach flow chart. Ocean thermal forcing and runoff fields are provided to the ice sheet modelers by the ISMIP6 ocean forcing working group. All parameters in the recommended melt-rate parameterization (Xu et al., 2013; Rignot et al., 2016) are provided, as well. Each ice sheet model must implement the approach and specify their own calving parameterization to calculate a frontal ablation rate for each outlet glacier.

ISMIP6 open approach is used sample a larger variety of oceanic forcing parameterizations, as it remains an active field of research. Models are free to continue applying the ocean forcing parameterization they used during the model initialization or their preferred method, but should still rely on the ocean forcing datasets provided by ISMIP6 to simulate future ocean conditions.

Requirements for the projections

• Participants can and are encouraged to contribute with different models and/or initialisation methods

• Models have to be able to prescribe a given SMB anomaly

• Models have to be able to prescribe a given ice front retreat or the high resolution approach for the standard experiments. For the open experiments, models can choose the ocean parameterization of their choice but should use the ocean forcing provided.

• Adjustment of SMB due to geometric changes in forward experiments is encouraged using the provided dSMB/dz.

• Bedrock adjustment in forward experiment is allowed.

• The choice of model input data is unconstrained to allow participants the use of their preferred model setup without modification. Modelers without preferred data set choice can have a look at the ISMIP6 Datasets page for possible options.

• To allow for analysis, any modeling choice needs to be well documented. A README file needs to be submitted along the outputs as an integral part of the contribution to the ISMIP6. It may be obtained here or requested by email to ismip6-at-gmail.com.

References

Alice Barthel, Cecile Agosta, Christopher M. Little, Tore Hatterman, Nicolas C. Jourdain, Heiko Goelzer, Sophie Nowicki, Helene Seroussi, Fiammetta Straneo, and Thomas J. Bracegirdle: CMIP5 model selection for ISMIP6 ice sheet model forcing: Greenland and Antarctica, The Cryosphere Discuss., https://doi.org/10.5194/tc-2019-191, 2019.

Franco, B., Fettweis, X., Lang, C., and Erpicum, M.: Impact of spatial resolution on the modelling of the Greenland ice sheet surface mass balance between 1990–2010, using the regional climate model MAR, The Cryosphere, 6, 695-711, https://doi.org/10.5194/tc-6-695-2012, 2012.

Goelzer, H., Nowicki, S., Edwards, T., Beckley, M., Abe-Ouchi, A., Aschwanden, A., Calov, R., Gagliardini, O., Gillet-Chaulet, F., Golledge, N. R., Gregory, J., Greve, R., Humbert, A., Huybrechts, P., Kennedy, J. H., Larour, E., Lipscomb, W. H., Le clec'h, S., Lee, V., Morlighem, M., Pattyn, F., Payne, A. J., Rodehacke, C., Rückamp, M., Saito, F., Schlegel, N., Seroussi, H., Shepherd, A., Sun, S., van de Wal, R., and Ziemen, F. A.: Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison, The Cryosphere, 12, 1433-1460, https://doi.org/10.5194/tc-12-1433-2018, 2018.

Goelzer, H., Noel, B. P. Y., Edwards, T. L., Fettweis, X., Gregory, J. M., Lipscomb, W. H., van de Wal, R. S. W., and van den Broeke, M. R.: Remapping of Greenland ice sheet surface mass balance anomalies for large ensemble sea-level change projections, The Cryosphere Discuss., 2019, 1-20, https://doi.org/10.5194/tc-2019-188, 2019.

Le clec'h, S., Charbit, S., Quiquet, A., Fettweis, X., Dumas, C., Kageyama, M., Wyard, C., and Ritz, C.: Assessment of the Greenland ice sheet–atmosphere feedbacks for the next century with a regional atmospheric model coupled to an ice sheet model, The Cryosphere, 13, 373-395, https://doi.org/10.5194/tc-13-373-2019, 2019.

Rignot, E., Xu, Y., Menemenlis, D., Mouginot, J., Scheuchl, B., Li, X., et al.: Modeling of ocean-induced ice melt rates of five west Greenland glaciers over the past two decades. Geophysical Research Letters, 43(12). http://doi.org/10.1002/2016GL068784, 2016.

Shannon, S.R., Payne A.J., Bartholomew I.D., Van Den Broeke M.R., Edwards T.L., Fettweis X., Gagliardini O., Gillet-Chaulet F., Goelzer H., Hoffman M.J., Huybrechts P.: Enhanced basal lubrication and the contribution of the Greenland ice sheet to future sea-level rise, Proceedings of the National Academy of Sciences, 110(35):14156-61. https://doi.org/10.1073/pnas.1212647110, 2013.

Slater, D. A., Straneo, F., Felikson, D., Little, C. M., Goelzer, H., Fettweis, X., and Holte, J.: Estimating Greenland tidewater glacier retreat driven by submarine melting, The Cryosphere, 13, 2489-2509, https://doi.org/10.5194/tc-13-2489-2019, 2019a.

Slater, D. A., Felikson, D., Straneo, F., Goelzer, H., Little, C. M., Morlighem, M., Fettweis, X., and Nowicki, S.: 21st century ocean forcing of the Greenland Ice Sheet for modeling of sea level contribution, The Cryosphere Discuss., 2019, 1-34, https://doi.org/10.5194/tc-2019-222, 2019b.

Xu, Y., Rignot, E., Fenty, I., Menemenlis, D., & Flexas, M. M.: Subaqueous melting of Store Glacier, west Greenland from three-dimensional, high-resolution numerical modeling and ocean observations. Geophysical Research Letters, 40(17). http://doi.org/10.1002/grl.50825, 2013.

Acknowledgements

The experimental protocol and datasets for the ISMIP6-Projections-Greenland standalone ice sheet simulations would not have been possible without the effort of many scientists that have given their time and expertise, and have run models to convert the CMIP5 models output into datasets that standalone ice sheet models can use. ISMIP6 would like to thank the ocean focus group under the leadership of Fiamma Straneo, the atmospheric focus group under the leadership of Bill Lipscomb and Robin Smith, and the CMIP5 model evaluation focus group under the leadership of Alice Barthel. Donald Slater, Denis Felixson, Mathieu Morlinghem and Heiko Goelzer have been instrumental in the development of the ice front retreat and melt parameterization and associated dataset. Xavier Fettweis, Patrick Alexander and Heiko Goelzer prepared the atmospheric dataset. Alice Barthel, Chris Little, Cecile Agosta, and Jamie Holte provided a rigorous analysis of the CMIP5 models against historical data, which allowed the CMIP5 model evaluation group and the ISMIP6 steering committee to select the CMIP5 models used in this effort. Finally, we thank the ISMIP6 ice sheet modelers for their feedback on the design of the protocol and their willingness to participate in ISMIP6.

Appendix 1 – Output grid definition and interpolation

All 2D data is requested on a regular grid with the following description. Polar stereo-graphic projection with standard parallel at 70° N and a central meridian of 45° W (315° E) on datum WGS84 (EPSG3413 projection). The lower left cell center is at (-720000m,-3450000m) with nx=1681 and ny=2881 cells in x and y-direction at full km positions (xmin = -720 km, xmax = +960 km, ymin = -3450 km, ymax = -570 km). The output should be submitted on a resolution adapted to the resolution of the model and can be 20 km, 10 km, 5 km, 2 km or 1 km. The data will be conservatively interpolated to 1 km resolution for archiving and 5 km resolution for diagnostic processing by ISMIP6.

If interpolation is required in order to transform the SMB forcing to your native grid, and transform your model variables to the ISMIP6 output grid (20 km, 10 km, 5 km, 2 km, 1 km), it is required that conservative interpolation is used. The motivation for using a common method for all models is to minimize model to model differences due to the choice of interpolation method.

Note: The previously requested regular grid was in polar stereo-graphic projection with standard parallel at 71° N and a central meridian of 39° W (321° E) on datum WGS84. The lower left corner is at (-800000 m, -3400000 m) and the upper right at (700000 m, -600000 m). This is the same grid (Bamber et al., 2001) used to provide the SMB anomaly forcing previously. This grid was changed to the EPSG3413 projection described above.

A1.1 Regridding Tools and Tips

  • An overview of the regridding process can be found on the Regridding page.
  • Regridding_with_CDO contains tools and tips that have been used by ISMIP6 members
  • ISMIP6 is designing tools to help with the regridding.
  • If you need help with conservative interpolation, please email ismip6-at-gmail.com.

Appendix 2 – Naming conventions, upload and model output data.

COMING SOON: GUIDANCE FOR CMIP6 MODEL REGISTRATION, ADDITIONAL GUIDANCE FOR FILE PREPARATIONS and REVISED FILENAME CONVENTION. HOWEVER, THE CURRENT FILE STRUCTURE AND GUIDANCE BELOW IS DESIGNED TO FACILITATE THE FINAL FILE PREPARATION (NEEDED BEFORE UPLOADING TO THE CMIP6 ARCHIVE). THE EXTRA INFORMATION/FILE RENAMING WILL BE IMPLEMENTED BY SCRIPTS CURRENTLY WRITTEN BY ISMIP6, WHICH WILL ALSO CHECK FOR CF COMPLIANCE ETC. WE WILL PROVIDE HELP TO MODELERS FOR THE FINAL FILE FORMATING PREPARATION, AND MODELERS SHOULD PROCEED WITH SAVING THEIR RESULTS USING THE INFORMATION BELOW

Please provide:

• one variable per file for all 2D fields and scalar variables.

• a completed readme file

A2.1 File name convention

File name convention for 2D fields and scalar variables:

<variable>_<IS>_<GROUP>_<MODEL>_<EXP>.nc

File name convention for readme file:

README_<IS>_<GROUP>_<MODEL>.doc

where

<variable> = netcdf variable name (e.g. lithk)

<IS> = ice sheet (AIS or GIS)

<GROUP> = group acronym (all upper case or numbers, no special characters)

<MODEL> = model acronym (all upper case or numbers, no special characters)

<EXP> = experiment name

For example, a file containing the variable "orog" for the Greenland ice sheet, submitted by group “JPL” with model “ISSM” for experiment “ctrl” would be called: orog_GIS_JPL_ISSM_ctrl.nc

If JPL repeats the experiments with a different version of the model (for example, by changing the sliding law), it could be named ISSM2, and so forth.

When uploading your files to the server, please follow the standard directory structure.

Example:

  group1
     model1
       exp05_05
       expC01_05
  group2
     model2
       exp05_01
       expC01_01

Group, model and experiment names have to be identical in the directory names and file names, except for the resolution suffix (rr=[01, 05, ..] the resolution the submitted data in km).

Example for results at 1 km resolution:

  AWI/ISSM1/ctrl_01/acabf_GIS_AWI_ISSM1_ctrl.nc
  xxx yyyyy zzzz rr           xxx yyyyy zzzz

The experiment names (exp_id) can be found at https://docs.google.com/spreadsheets/d/1YEiPV3Uc0K8EqD57mz-ZqBIF4XnVb5UbyRDRmox0KeE/edit#gid=1531275537

A few examples:

exp_id RCP GCM Ocean
exp05 8.5 MIROC5 Medium
exp06 8.5 NorESM Medium
expa01 8.5 IPSL-CM5-MR Medium
expb01 ssp585 CNRM-CM6 Medium
expc01 8.5 MIROC5 AO Medium
expc01 8.5 NorESM High


A2.2 Retrieving dataset and Uploading your model output

A2.2.1 How to gain access to server

All ISMIP6 datasets are stored at the University of Buffalo’s CCR transfer server (UB CCR).

In order to gain access to the ISMIP6 datasets and upload or download files please send a request to ismip6 [at] gmail.com address with your public ssh key, name, and affiliation for account setup.

To generate your public ssh key use the command below:

ssh-keygen -t rsa

Once your account is set up you will receive a confirmation email with your username.

All the current and previous ISMIP6 related information and datasets are located in the sftp://transfer.ccr.buffalo.edu/projects/grid/ghub/ISMIP6 folder.

You can access the server and upload or download datasets on the UB CCR server the following ways below.

Command line access with sftp:

sftp -o "IdentityFile=sshkey" yourusername@transfer.ccr.buffalo.edu
 

Upload files

put /path/to/your/file /path/to/folder/on/server

Upload folder

put -r /path/to/your/folder /path/to/folder/on/server

Download file

get /path/to/file/on/server /path/to/your/folder/

Download whole directory

get -r /path/to/folder/on/server /path/to/your/folder/

Commands for rsync ssh access:

Copy a file from local folder to remote folder:

rsync -avz -e "ssh -i sshkey" /home/localuser/testfile.txt remoteuser@transfer.ccr.buffalo.edu:/home/remoteuser/transfer/

Copy a file from remote server into a local folder:

rsync -avz -e "ssh -i sshkey" remoteuser@transfer.ccr.buffalo.edu:/home/remoteuser/transfer/testfile.txt  /home/localuser/

FileZilla/Cyberduck Setup

Please see link below for easy FileZilla setup https://www.a2hosting.com/kb/getting-started-guide/accessing-your-account/setting-up-sftp-access-for-multiple-users#Configure-an-FTP-client-to-use-the-private-key

Please set the PORT to 22.

A2.2.2 Where to upload your results

After log in, go to the ISMIP6/Projections/GrIS/output directory via:

sftp> cd /ISMIP6/Projections/GrIS/output

and create a directory named <GROUP> with the following sub-directory structure:

sftp> mkdir /Projections/GrIS/output/<GROUP>/<MODEL>/<EXP_RES>

<EXP_RES> should include both the experiment name and the output grid used to simplify the processing (e.g., asmb_08). Only the directory name should include this resolution, unlike the output files.

Create additional <MODEL> directories when participating with more than one model or model version.

A2.2.3 Reducing the size of files

The size of the model files on higher resolution grid can be largely reduced by file compression which will save space on the storage server. An example command is given below and the results before and after. In the examples that follow we can get a factor of 10 compression and for the masks even more given that contiguous masks are highly compressible because they are repeated data. NetCDF files have been designed with compression in mind. A NetCDF file can be compressed and nothing has to be changed in the way that it is read into Matlab or Python (or any other language that uses standard NetCDF read/write libraries).

The nccopy command copies an input netCDF file to an output netCDF file after compressing the file significantly. The ‘-d’ option stands for the deflation level, from 1 (faster but lower compression) to 9 (slower but more compression) and the ‘-s’ option is the shuffling option to improve compression even more. We recommend using ‘d1’ option since this option seems to accomplish the desired compression.

Example of netcdf compression command:

nccopy -d1 -s sftgif_GIS_JPL_ISSMPALEO_historical.nc sftgif_GIS_JPL_ISSMPALEO_historical_c.nc

Example of compression variant, seems to work better for masks:

nccopy -d1 sftgif_GIS_JPL_ISSMPALEO_historical.nc sftgif_GIS_JPL_ISSMPALEO_historical_c.nc

A2.3 Model output variables and README file

The README file is an important contribution to the ISMIP6 submission. It may be obtained here or requested by email to ismip6-at-gmail.com

A2.3.1 General guidelines

The variables requested in the table below serve to evaluate and compare the different models and initialization techniques. Some of the variables may not be applicable for your model, in which case they are to be omitted (with explanation in the README file). Also, specify missing values in your netcdf file where needed, and fields should be undefined outside of the ice mask.

All "missing data" must be assigned the single precision floating point value of 1.e20

We distinguish between state variables "ST" (e.g. ice thickness, temperatures and velocities) and flux variables "FL" (e.g. SMB). State variables should be given as snapshot information at the end of one year for both scalars and 2D variables (for initMIP, 2D variables were only requested over five year periods), while flux variables are to be averaged over the respective periods. Please specify in your README file how your reported flux data has been averaged over time. Ideally, the standard would be go average over all native time steps.

Flux variables are defined positive when the process adds mass to the ice sheet and negative otherwise.

If you redo the initMIP experiments (because you have a new initial state for the projections), please save the files at a yearly interval instead of the 5 years interval requested as part of the original initMIP. Also upload your iniMIP results in the projections directory.

A2.3.2 How to record time in historical and projection files

In compliance with CMIP6, time should be defined in "days since <basetime>", where <basetime> must be specified by the user, typically in the form year-month-day (e.g., "days since 1800-1-1"). For simulations meant to represent a particular historical period, set the ‘base time’ to the time at the beginning of the simulation. A historical run initialized with forcing for year 2007 would, for example, have units of “days since 2007-1-1”. For the future scenario runs, retain the same <basetime> as used in the historical run from which it was initiated. Note the CF definition for years (section 4.4): a common_year is 365 days, a leap_year is 366 days, a Julian_year is 365.25 days, a Gregorian_year is 365.2425 days, a 360_day has all years with 360 days divided into 30 day months (please see the CF link above for other examples on calendar setting in section 4.4).

To illustrate a time recording for the historical file and projections for a typical state variable (ST, eg thickness) and flux variable (FL, eg SMB), we assume that our <basetime> is January 1st 2013, and that we use a calendar = 360_days. Other calendars can be used, but you need to indicate the calendar used in the netcdf, and of course if you use a different calendar, the time entries will be different. What needs to be recorded is shown in green in the Table below. For state variables,

Figure A1.png


For state variables, like thickness for the historical:

dimensions:
   time = UNLIMITED ; // (3 currently)
variables:
   double time(time) ;
       time:units = "days since 1-1-2013" ; // This date correspond to the example basetime 
       time:calendar = "360_day" ; // Other calendars can be used... change here to relevant calendar
       time:axis = "T" ;
       time:long_name = "time" ;
       time:standard_name = "time" ;
data:
 time = 0, 360, 720; // If you use a different calendar these values will change

and thickness for the projection (note that the full time entries are not shown, only beginning and end) would be:

dimensions:
   time = UNLIMITED ; 
variables:
   double time(time) ;
       time:units = "days since 1-1-2013" ; // This date correspond to the example basetime
       time:calendar = "360_day" ; // Other calendars can be used... change here to relevant calendar
       time:axis = "T" ;
       time:long_name = "time" ;
       time:standard_name = "time" ;
data:
 time = 720, 1080, 1440, …, 31320, 31680; // If you use a different calendar these values will change

The flux variable, like SMB, would be recorded as the average over a full year, so for the historical:

dimensions:
   time = UNLIMITED ; // (2 currently)
   bnds = 2 ;
variables:
   double time(time) ;
       time:bounds = "time_bnds" ;
       time:units = "days since 1-1-2013 " ; // This date correspond to the example basetime
       time:calendar = "360_day" ; // Other calendars can be used... change here to relevant calendar
       time:axis = "T" ;
       time:long_name = "time" ;
       time:standard_name = "time" ;
   double time_bnds(time, bnds) ;
data:
 time = 180, 540 ; // If you use a different calendar these values will change. 
                             //This is the middle of the time_bnds
 time_bnds =
  0, 360, //If you use a different calendar these values will change. 
              //These are the day since basetime at the beginning and end of the year
  360, 720 ;

and the projection (note that the full time entries are not shown, only beginning and end):

dimensions:
   time = UNLIMITED ; // 
   bnds = 2 ;
variables:
   double time(time) ;
       time:bounds = "time_bnds" ;
       time:units = "days since 1-1-2013 " ; // This date correspond to the example basetime
       time:calendar = "360_day" ;  // Other calendars can be used... change here to relevant calendar
       time:axis = "T" ;
       time:long_name = "time" ;
       time:standard_name = "time" ;
   double time_bnds(time, bnds) ;
variables:
   double time(time) ;
       time:bounds = "time_bnds" ;
data:
 time = 900, 1260, 1620, ..., 31140, 31500; // If you use a different calendar these values will change. 
                                                     //This is the middle of the time_bnds
 time_bnds =
720, 1080, //If you use a different calendar these values will change. 
                 //These are the day since basetime at the beginning and end of the year
1080, 1440,
1440, 1800,
....
30960, 31320,
31320, 31680;

A2.3.3 Table A1: Variable request for ISMIP6

If your quantity does not change with time, then simply save one time entry. An example is geothermal heat flux, which varies in some models but not others.

Table A1: Variable request for ISMIP6 projections.

Bold names or "alias" indicate a change compared to initMIP, to align the request with the CMIP6 official MIPtable "IyrGre" or names in the CF convention. If possible please use the new names, and if not, the name change will occur when your files are checked for CMIP compliance.

The first entry should be that from which the simulation starts. Fields such as surface mass balance flux should be what was applied as boundary conditions.

Variable Dim Type Variable Name Standard Name Units Comment
2D variables requested yearly as snapshots (end of the year) for type ST and as yearly average for type FL.
Ice thickness x,y,t ST lithk land_ice_thickness m The thickness of the ice sheet
Surface elevation x,y,t ST orog surface_altitude m The altitude or surface elevation of the ice sheet
Bedrock elevation x,y,t ST topg bedrock_altitude m The bedrock topography (may change during the projections)
Geothermal heat flux x,y,t FL hfgeoubed upward_geothermal_heat_flux_in_land_ice alias "upward_geothermal_heat_flux_at_ground_level" W m-2 Geothermal Heat flux at the ice interface (only needed beneath the grounded ice). If this quantity does not change with time simply enter one timestep
Surface mass balance flux x,y,t FL acabf land_ice_surface_specific_mass_balance_flux kg m-2 s-1 Surface Mass Balance flux (for areas covered by ice only)
Basal mass balance flux beneath grounded ice x,y,t FL libmassbfgr alias "libmassbf" land_ice_basal_specific_mass_balance_flux kg m-2 s-1 Basal mass balance flux (only beneath grounded ice)
Basal mass balance flux beneath floating ice x,y,t FL libmassbffl alias "libmassbf" land_ice_basal_specific_mass_balance_flux kg m-2 s-1 Basal mass balance flux (only beneath floating ice)


Ice thickness imbalance x,y,t FL dlithkdt tendency_of_land_ice_thickness m s-1 dHdt
Surface velocity in x x,y,t ST xvelsurf alias "uvelsurf" land_ice_surface_x_velocity m s-1 u-velocity at land ice surface
Surface velocity in y x,y,t ST yvelsurf alias "vvelsurf" land_ice_surface_y_velocity m s-1 v-velocity at land ice surface
Surface velocity in z x,y,t ST zvelsurf alias "wvelsurf" land_ice_surface_upward_velocity m s-1 w-velocity at land ice surface
Basal velocity in x x,y,t ST xvelbase alias "uvelbase" land_ice_basal_x_velocity m s-1 u-velocity at land ice base
Basal velocity in y x,y,t ST yvelbase alias "vvelbase" land_ice_basal_y_velocity m s-1 v-velocity at land ice base
Basal velocity in z x,y,t ST zvelbase alias "wvelbase" land_ice_basal_upward_velocity m s-1 w-velocity at land ice base
Mean velocity in x x,y,t ST xvelmean alias "uvelmean" land_ice_vertical_mean_x_velocity m s-1 The vertical mean land ice velocity is the average from the bedrock to the surface of the ice
Mean velocity in y x,y,t ST yvelmean alias "vvelmean" land_ice_vertical_mean_y_velocity m s-1 The vertical mean land ice velocity is the average from the bedrock to the surface of the ice
Surface temperature x,y,t ST litemptop alias "litempsnic" temperature_at_top_of_ice_sheet_model alias "temperature_at_ground_level_in_snow_or_firn" K Ice temperature at surface
Basal temperature beneath grounded ice sheet x,y,t ST litempbotgr alias "litempbot" temperature_at_base_of_ice_sheet_model alias "land_ice_basal_temperature" K Ice temperature at base of grounded ice sheet
Basal temperature beneath floating ice shelf x,y,t ST litempbotfl alias "litempbot" temperature_at_base_of_ice_sheet_model alias "land_ice_basal_temperature" K Ice temperature at base of floating ice shelf
Basal drag x,y,t ST strbasemag land_ice_basal_drag alias "magnitude_of_land_ice_basal_drag" Pa Basal drag
Calving flux x,y,t FL licalvf land_ice_specific_mass_flux_due_to_calving kg m-2 s-1 Loss of ice mass resulting from iceberg calving. Only for grid cells in contact with ocean
Ice front calving and melt flux x,y,t FL lifmassbf land_ice_specific_mass_flux_due_to_calving_and_ice_front_melting kg m-2 s-1 Loss of ice mass resulting from calving and ice front melting. Only for grid cells in contact with ocean
Land ice area fraction x,y,t ST sftgif land_ice_area_fraction 1 Fraction of grid cell covered by land ice (ice sheet, ice shelf, ice cap, glacier)
Grounded ice sheet area fraction x,y,t ST sftgrf grounded_ice_sheet_area_fraction 1 Fraction of grid cell covered by grounded ice sheet, where grounded indicates that the quantity correspond to the ice sheet that flows over bedrock
Floating ice sheet area fraction x,y,t ST sftflf floating_ice_shelf_area_fraction alias "floating_ice_sheet_area_fraction" 1 Fraction of grid cell covered by ice sheet flowing over seawater
Scalar outputs requested every full year: snapshots for type ST and 1 year averages for type FL.
Total ice mass t ST lim land_ice_mass kg spatial integration, volume times density
Mass above floatation t ST limnsw land_ice_mass_not_displacing_sea_water kg spatial integration, volume times density
Grounded ice area t ST iareagr alias "iareag" grounded_ice_sheet_area m^2 spatial integration
Floating ice area t ST iareafl alias "iareaf" floating_ice_shelf_area m^2 spatial integration
Total SMB flux t FL tendacabf tendency_of_land_ice_mass_due_to_surface_mass_balance kg s-1 spatial integration
Total BMB flux t FL tendlibmassbf tendency_of_land_ice_mass_due_to_basal_mass_balance kg s-1 spatial integration
Total BMB flux beneath floating ice t FL tendlibmassbffl tendency_of_land_ice_mass_due_to_basal_mass_balance kg s-1 spatial integration (computed beneath floating ice only)
Total calving flux t FL tendlicalvf tendency_of_land_ice_mass_due_to_calving kg s-1 spatial integration
Total calving and ice front melting flux t FL tendlifmassbf tendency_of_land_ice_mass_due_to_calving_and_ice_front_melting kg s-1 spatial integration

Appendix 3 – Participating Models and Characteristics

Greenland Standalone Ice Sheet Modeling

Contributors Model Group ID Group
Martin Rückamp, Angelika Humbert ISSM AWI Alfred Wegener Institute for Polar and Marine Research, DE /University of Bremen, DE
Victoria Lee, Tony Payne, Stephen Cornford, Daniel Martin BISICLES BGC University of Bristol, Bristol, UK, Department of Geography, Swansea University, UK, Computational Research Division, Lawrence Berkeley National Laboratory, California, USA
Isabel Nias, Sophie Nowicki, Denis Felikson ISSM GSFC NASA Goddard Space Flight Center, Greenbelt, USA
Ralf Greve, Reinhard Calov, Chris Chambers SICOPOLIS ILTS_PIK Institute of Low Temperature Science, Hokkaido University, Sapporo, JP, Potsdam Institute for Climate Impact Research, Potsdam, DE
Heiko Goelzer, Roderik van de Wal, Michiel van den Broeke IMAUICE IMAU Utrecht University, Institute for Marine and Atmospheric Research (IMAU), Utrecht, NL
Helene Seroussi, Nicole Schlegel ISSM JPL NASA Jet Propulsion Laboratory, Pasadena, USA
Joshua K. Cuzzone, Nicole Schlegel ISSMPALEO JPL NASA Jet Propulsion Laboratory, Pasadena, USA
Aurélien Quiquet, Christophe Dumas GRISLI LSCE LSCE/IPSL, Laboratoire des Sciences du Climat et de l’Environnement, CEA-CNRS-UVSQ, Gif-sur-Yvette, FR
Lev Tarasov GSM MUN Dept of Physics and Physical Oceanography, Memorial University of Newfoundland, USA
William Lipscomb, Gunter Leguy CISM NCAR National Center for Atmospheric Research, Boulder, CO, USA
Andy Aschwanden PISM UAF Geophysical Institute, University of Alaska Fairbanks, USA
Youngmin Choi, Helene Seroussi, Mathieu Morlighem ISSM UCIJPL NASA Jet Propulsion Laboratory, Pasadena, USA / University of California Irvine, Irvine, USA
Sébastien Le Clec’h and Philippe Huybrechts GISM VUB Vrije Universiteit Brussel, Brussels, BE
Nick Golledge, Dan Lowry PISM VUW GNS Science, Lower Hutt, NZ / Antarctic Research Centre, Victoria University of Wellington, NZ

Model Characteristics

The Model Characteristics table can be found here.