ISMIP6 Logo
  • Overview
  • Members
  • Wiki14
  • Resources
  • Forum
  • Blog
  • Wish List
  • Usage
  • Calendar
  • Announcements
  • Projects
  • Files
  • Activity
  • Discoverability Visible
  • Join Policy Restricted
  • Created 01 Jul 2022

ISMIP6

►

Wiki

ISMIP6 Projections 2300 Antarctica

Version 34
by (unknown)
Version 35
by (unknown)

Deletions or items before changed

Additions or items after changed

1 [[Image(ice_sheet_image_thin.png)]]
2
3 == Overview ==
4 This page describes the experimental protocol for the ISMIP6 2300 projections that focus on simulations of the Antarctic Ice Sheet (AIS) extended to year 2300.
5
6 These simulations are based on CMIP5 and CMIP6 climate model outputs, and are a follow-on to the simulations to 2100 described in [ISMIP6-Projections-Antarctica] and the papers by [https://tc.copernicus.org/articles/14/2331/2020/ Nowicki et al. (2020)] and [https://tc.copernicus.org/articles/14/3033/2020/ Seroussi et al. (2020)].
7
8 Some experiments use climate forcing from coupled global models that were run until 2300 under CMIP forcing scenarios, while other experiments use repeated forcing from the 2080-2100 period, sampled randomly between 2100 and 2300.
9
10 -
This dataset is only available for ISMIP6 ice sheet modeling groups and can be obtained via a request to ismip6 - at- gmail.com. Please see also the bottom of the page for guidance on usage and how to refer this dataset.
+
This experiment is now closed and the dataset is available via Ghub/Globus datasets. Please see also the bottom of the page for guidance on usage and how to refer this dataset.
11
12 ----
13 == List of Projections ==
14 Two CMIP5 models (CCSM4 and !HadGEM2) and two CMIP6 models (CESM2 and UKESM) were run with extended high CO,,2,, forcing to 2300 (RCP8.5 and ssp5-85, respectively, for CMIP5 and CMIP6) and were selected for long-term projections.
15 No !NorESM1-M extension to 2300 is available, so the extended experiments with !NorESM forcing use repeat forcing only.
16
17 '''There are 14 experiments in all'''. Each experiment ID begins with the prefix ''''AE'''' to signify ''''Antarctic extension''''. '''Table 1''' lists six '''''Tier-1''''' experiments that we ask each group to run if possible. This selection includes one run for each of the four climate models with extended forcing to 2300; a low-forcing run for comparison to the high-forcing runs; and one run with repeat forcing from the late 21^st^ century forcing for comparison to runs with extended forcing. If a group is unable to run all six '''''Tier-1''''' experiments, we ask that they choose a subset.
18
19 '''Table 2''' lists an additional eight '''''Tier-2''''' experiments that we encourage each group to run if resources allow. These include three more experiments with repeat late 21^st^ century forcing, for comparison to the runs with extended forcing. There is one experiment with low emissions (ssp1-26) verextended to 2300, to complement the RCP2.6 experiment in '''''Tier-1'''''. Finally, there are four experiments with prescribed ice-shelf collapse driven by hydrofracture, to compare to the runs without shelf collapse. We excluded shelf-collapse runs from '''''Tier-1''''' because not all ice sheet models might have this capability.
20
21 Modeling groups that can run many simulations are encouraged to further explore the ice sheet response using ''targeted experiments''.
22 For these experiments, groups should repeat the runs in '''Tables 1''' and '''2''' with ''either'' atmospheric or ocean forcing enabled, but not both. This will help us analyze the relative contributions of atmospheric and ocean forcing, especially for the runs with more extreme warming. The experiment numbers are as in '''Tables 1''' and '''2''', but with a lower-case '''a''' or '''o''' appended to show which forcing is applied.
23
24 For the original ISMIP6 projections, each group was asked to run with a ''standard'' sub-shelf melting parameterization based on the method described by [https://tc.copernicus.org/articles/14/3111/2020/ Jourdain et al. (2020)],
25 and optionally an ''open'' parameterization chosen by the group. For the extended Antarctic projections, we no longer prescribe a standard method. Instead, sub-shelf melting schemes are left to the discretion of each group. This change reflects our desire to fully sample the methods in use. As in the original projections, we ask each group to use the thermal forcing data provided by ISMIP6, so that different ice sheet model responses can be attributed to differences in the models rather than the forcing.
26
27 '''''Note:''' All datasets needed for the '''''Tier-2''''', '''''Tier-2''''', and targeted experiments are available on Ghub's [https://theghub.org/dataset-listing Browse Datasets]. The same datasets are used for the additional targeted experiments. (See [https://docs.google.com/spreadsheets/d/1J1AjPrnOLqXu-n8BJeYjZqdzXKv6g3wGRWmZt0tGuNc/edit#gid=0 Table])''
28
29
30 ||||||||||||'''Table 1: Tier-1 Experiments'''||
31 ||'''Exp''' ||'''Model'''||'''Scenario'''||'''Forcing'''||'''Collapse'''||'''Notes'''||
32 ||'''AE01'''||!NorESM1-M||RCP2.6||Repeat||No||Low warming scenario||
33 ||'''AE02'''||CCSM4||RCP8.5||To 2300||No||Extended high-emissions CMIP5 scenario||
34 ||'''AE03'''||!HadGEM2||RCP8.5||To 2300||No||Extended high-emissions CMIP5 scenario||
35 ||'''AE04'''||CESM2||ssp5-85||To 2300||No||Extended high-emissions CMIP6 scenario||
36 ||'''AE05'''||UKESM||ssp5-85||To 2300||No||Extended high-emissions CMIP6 scenario||
37 ||'''AE06'''||UKESM||ssp5-85||Repeat||No||Repeat forcing for comparison to AE05||
38 ||||||||||||'''Table 2: Tier-2 Experiments'''||
39 ||'''Exp''' ||'''Model'''||'''Scenario'''||'''Forcing'''||'''Collapse'''||'''Notes'''||
40 ||'''AE07'''||!NorESM1-M||RCP8.5||Repeat||No||Extension (with repeat forcing) of ISMIP6-Antarctica Exp 1||
41 ||'''AE08'''||!HadGEM2||RCP8.5||Repeat||No||Repeat forcing for comparison to AE03||
42 ||'''AE09'''||CESM2||ssp5-85||Repeat||No||Repeat forcing for comparison to AE04||
43 ||'''AE10'''||UKESM||ssp1-26||To 2300||No||Extended low warming scenario||
44 ||'''AE11'''||CCSM4||RCP8.5||To 2300||Yes||Collapse experiment for comparison to AE02||
45 ||'''AE12'''||!HadGEM||RCP8.5||To 2300||Yes||Collapse experiment for comparison to AE03||
46 ||'''AE13'''||CESM2||ssp5-85||To 2300||Yes||Collapse experiment for comparison to AE04||
47 ||'''AE14'''||UKESM||ssp5-85||To 2300||Yes||Collapse experiment for comparison to AE05 ||
48
49 ----
50 [[Image(ice_sheet_image_thin_2.png)]]
51
52 == Initialization, historical run, control run, and projection runs ==
53
54 All projection experiments start on 1 January 2015 and end on 31 December 2300. The start date follows the CMIP6 protocol for projections, while the end date is constrained by the availability of forcing.
55
56 The '''''initialization date''''' (or initial state) is left to each group's discretion and can be any time before January 2015. The initialization date corresponds to the date assigned to the initialization procedure.
57
58 In many cases, a short '''''historical run''''' will be needed to bring the models from the initialization date (say, 1990) to the '''''projection start date''''' of January 2015.
59 Each model configuration should have a '''single''' historical run, from which all the projections will branch. Groups are free to choose the forcing for the historical run – for example, using a reanalysis, historical forcing from an RCM or AOGCM, or a combination of multiple datasets. Groups should not carry out a separate historical run for each AOGCM experiment, because this would
60 complicate the forcing strategy and interpretation. Models without a historical run, e.g. that start directly in January 2015, should report their initial conditions as the historical run.
61
62 In addition to the projection runs, each model configuration should have a ''projection control run'' (`ctrl_proj`). This is an unforced simulation that starts in January 2015 in the same ice sheet state as the projections, and continues through 2300. It is meant to capture how much drift arises from the historical run. The projection control run is implemented with ''zero anomalies'' relative to the 20-year AR6 reference period of January 1995 through December 2014. Below, we offer guidance on choosing SMB and ocean forcing climatologies for the reference period.
63
64 Some experiments (#1, 2, 3, 5, 6, 7, 8, 11, 12 and 14 in '''Tables 1''' and '''2''') use the same forcing for years 2015–2100 as in the original Antarctic projections.
65
66 However, the CESM2 ssp5-85 forcing (#4, 9 and 13) has changed, coming from a run with a different atmosphere component (WACCM instead of CAM). The low-emission UKESM forcing (#10) was not part of the original projections. For the experiments with identical 21^st^ century forcing, groups that already ran their models through 2100 can simply continue from 1 January 2101, ''provided their ice sheet models have not changed''. If the models have changed, we ask that groups repeat their initialization and historical run, and then start the projections from 2015.
67
68 For groups choosing AOGCM forcing for the initialization or historical run, ISMIP6 provides an SMB and surface temperature climatology, along with anomalies, for each AOGCM used to generate a projection dataset. For Antarctica, the SMB and temperature climatology corresponds to 1995–2014, to align with the AR6 reference period. Antarctic SMB and temperature anomalies are available from 1950. For the Southern Ocean, the datasets start from 1850, and the climate model climatology corresponds to 1995–2014. Groups using an Antarctica dataset provided by ISMIP6 are recommended to use the !NorESM1-M climatology and anomalies for SMB and surface temperature (in the directory Atmosphere_Forcing/noresm1-m_rcp8.5). See '''Appendix A2.2''' for accessing the forcing data and directories. For the ocean, modelers can use observational climatology (ISMIP6 2100 dataset forcing directory via /AIS/Ocean_Forcing/climatology_from_obs_1995-2017 directory) and/or anomalies (Ocean_Forcing/noresm1-m_rcp8.5/1850-1994 directory). The climate model ocean climatologies (Ocean_Forcing/noresm1-m_rcp8.5/climatology_1995-2014 directory) are not intended for use by modelers, but are provided so that users can see what was subtracted during the dataset preparation, as the ocean forcing data (Ocean_Forcing/noresm1-m_rcp8.5/1850-1994 directory) is the sum of the observational climatology and model anomalies.
69
70 To better sample uncertainties, we encourage groups to submit results using more than one model configuration – for example, from a model run at two or more grid resolutions, or with substantially different physics options. Each configuration would be associated with a separate suite of up to 14 extension experiments. In this case, it is appropriate to do an independent '''initialization''', '''historical run''', and '''projection control run''' for each configuration.
71
72 ----
73 == Atmospheric forcing: SMB and temperature anomalies ==
74 ISMIP6 provides surface forcing datasets for the AIS based on CMIP AOGCM simulations.
75
76 '''Two approaches''' are possible:
77 # using AOGCM output directly, or
78 # re-interpreting the GCM climates through higher-resolution regional climate models (RCMs).
79
80 The latter approach, which better captures large surface mass balance (SMB) gradients regions near the periphery of the ice sheet, has been used for ISMIP6 '''Greenland experiments''' ('''Goelzer et al.''', 2020).
81
82 For the '''Antarctic experiments''', RCMs are not used, so SMB anomalies based on AOGCM output are applied directly.
83 Several AOGCMs now use multiple elevation classes to downscale SMB to finer grid resolution over ice sheets. In the future, we may add ISMIP6 experiments using SMB downscaled via elevation class.
84
85 For the ISMIP6 projections based on CMIP5 and CMIP6 AOGCMs, the surface forcing consists of anomalies in SMB and surface temperature (illustrated in '''Fig. 2'''). 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. The following remarks refer mostly to SMB, but the same comments would generally apply to surface temperature as well.
86
87 ISMIP6 provides yearly averaged surface mass balance anomalies, `aSMB(x,y,t)`, along with its components (precipitation, evaporation and runoff), as well as the SMB climatologies used to compute the anomalies:
88
89 `aSMB_AOGCM(x,y,t) = SMB_AOGCM(x,y,t) - SMB_CLIM_AOGCM(x,y)`
90
91
92 '''where''' `SMB_AOGCM` is the SMB for a given AOGCM and `SMB_CLIM_AOGCM` is the climatology for that AOGCM. The `SMB_CLIM_AOGCM` were computed by taking the temporal average of all `SMB_AOGCM` over the reference period (from January 1995 to December 2014). ISMs can use these climatologies for initialization runs, if desired, but are free to use their preferred SMB forcing for these runs.
93
94 During the projection run, modelers need to reintroduce the climatology that best fits their simulations. SMB is computed as:
95
96 `SMB(x,y,t) = SMB_ref(x,y) + aSMB(x,y,t)`
97
98
99 '''where''' `SMB_ref` is the SMB that the ice sheet model would have used over the reference period (from January 1995 to December 2014) and should be the same for all experiments. If a time-dependent SMB is used, then `SMB_ref(x,y)` is the average over the reference period. If an SMB climatology is used, then `SMB_ref(x,y)` is simply the climatology. ISMIP6 accepts that existing climatologies (or datasets of SMB averaged over many years) may not align exactly with the AR6 reference period.
100
101 However, we assume that the differences between climatologies will be less than the inter-annual variability of the SMB derived from the AOGCMs, and thus changes in `aSMB`. What is important is that `SMB_ref` is computed over many years. This assumption allows each ISM to use their preferred SMB.
102
103 `aSMB(x,y,t)` is constant over the entire year and changes stepwise at the beginning of the following year. SMB climatologies and anomalies are given in units of kg m^-2^ s^-1^ (water equivalent), and surface temperature in units of deg K. To convert `aSMB` to units [!m yr^-1^ ice] typically used in an ice sheet model, multiply the `netcdf` variable by 31,556,926 s/yr, 1/1000 m^3^/kg and by the density ratio ρ,,w,,/ρ,,i,,:
104
105 `aSMB [m yr^-1^] = aSMB [kg m^-2^ s^-1^] * 31,556,926 / 1000 * (1000/ρ^i^)`
106
107
108 '''where''' ρ,,w,, = 1000 kg m^3^ is the density of water, and ρi is your specific ice density (typically 917.0 kg m^3^ or similar).
109
110 The datasets can be obtained from the ISMIP6 2300 Forcing Globus endpoint AIS/Atmosphere_Forcing/. See 2300 Projections Antarctic files at the end of this wiki to set up a Globus account. Files are provided for several resolutions (1 km, 2 km, 4 km, 8 km, 16 km, and 32 km). Modeling groups should use the resolution closest to their native grid to conservatively interpolate the data to the model grid (see '''Appendix 1''', below).
111
112
113 [[Image(SmbAnomaly.png)]]
114
115 '''Figure 2a''': SMB anomaly (Gt/yr) timeseries, 1950-2300, for CCSM4, CESM2, !HadGEM2, UKESM, and !NorESM1. All forcings are from high emission scenarios except !NorESM1 RCP2.6. !NorESM1 anomalies are repeat forcings (see overview).
116
117 ----
118 [[Image(Smb2300All.png, px=750)]]
119
120 '''Figure 2b''': Change in SMB between the projection start and end date (2300 minus 2015) in units of meters ice equivalent for the AOGCMs shown in '''Fig. 2a'''. '''Note''' the uneven colorbar.
121
122 ----
123 == Oceanic forcing: temperature, salinity, thermal forcing and melt rate parameterization==
124 ISMIP6 provides datasets of extrapolated ocean "ambient" temperature ('''T'''), salinity ('''S''') and thermal forcing ('''TF''') from 1850—2300 that are appropriate for present and future ice-shelf cavities. These datasets originate from CMIP models and have been extrapolated under ice shelves by Xylar Asay-Davis, using rules that account for sills and troughs ('''Fig. 3'''). The datasets are on the ISMIP6 8-km Antarctic grid. For more information on how the datasets were produced, please see '''Nowicki et al.''' (2020), '''Jourdain et al.''' (2020), and the [https://github.com/xylar/ismip6-ocean-forcing ISMIP6 Antarctic Ocean Forcing GitHub].
125
126
127 [[Image(Ocean_overview_antarctica_basin.png, 850px)]]
128 ----
129
130 '''Figure 3''': Bathymetry and IMBIE2 basins (left) used in the sub-ice shelf extrapolation of ocean temperature (right).
131
132 '''Modeling groups are free to use their sub-ice-shelf melt rate parameterization of choice''', provided that the parameterization uses the ocean forcing datasets (T, S, and TF) provided by ISMIP6. The temperature, salinity and thermal forcing data provided for CMIP5 and CMIP6 models are the anomalies of each model with respect to its January 1995—December 2014 average, added to an observational climatology (based on WOA, EN4 and MEOP datasets). Thus, the datasets can be used directly by models without computing anomalies or selecting reference observations of your own. However, groups might need to compute anomalies. In such cases, anomalies should be computed with respect to the January 1995—December 2014 average, as this period was used to anomalize the CMIP5 model input (and is slightly different from the time period, 1995—2017, spanned by the observations). '''Groups with multiple submissions using different melt rate parameterizations''' should carry out all the '''''Tier-1''''' experiments ('''Table 1'''), and optionally the '''''Tier-2''''' experiments ('''Table 2'''), for each submission. Each submission will have an independent initialization and historical run.
133
134 Groups may use the '''quadratic-dependence melt rate parameterization''' that was developed by the Antarctic ocean focus group and was the standard melt parameterization in ISMIP6 Antarctica ('''Seroussi et al.''', 2020). This parameterization is described in detail in '''Nowicki et al.''' (2020), '''Jourdain et al.''' (2020), [ISMIP6-Projections-Antarctica], and '''Favier et al.''' (2019). Modeling groups are free to use either median, 5th percentile, or 95th percentile values of the γ,,,0,,, and ΔT parameters (provided in /Ocean_Forcing/parameterizations directory), so long as the parameter choice remains consistent. We encourage groups to explore the sensitivity of the melt parameterizations to the values of γ,,,0,,, and ΔT. In this case, groups should submit independent initialization and historical runs for each configuration.
135
136 The datasets can be obtained from the ISMIP6 2300 Forcing Globus endpoint AIS/Ocean_Forcing/ .
137
138 == Antarctic ice shelf collapse ==
139 Surface melting can trigger ice shelf collapse (for example, the Larsen B ice shelf in the Antarctic Peninsula). This mechanism is distinct from, but a precursor to, cliff collapse. Although the mechanisms for Larsen B-style ice shelf collapse are not well understood, ISMIP6 provides datasets for ice shelf collapse in the form of a time-dependent mask ('''Fig. 4'''). These datasets were derived from AOGCM near-surface air temperature (tas) using the method described in '''Trusel et al.''' (2015) to compute annual surface melt.
140
141 For ISMIP6, Luke Trussel prepared the bias-corrected annual surface melt, which was used to generate the masks. Ice shelves are assumed to collapse following a 10-year period with average surface melt above 725 mm/yr. Some experiments (see '''Table 2''') require modeling ice shelf collapse using the ISMIP6 masks. Ice shelf collapse is included only in '''''Tier-2''''' experiments and is not a requirement for ISMIP6-2300 participation.
142
143 Models are free to decide on the appropriate method to simulate tributary glaciers' behavior following ice shelf collapse. Since the masks were derived from observations, the observed ice shelf may not always corresponds to an ice shelf in the ISM. In the event that the ice shelf collapse mask corresponds to a region which the ISM considers to be part of the grounded ice sheet, the collapse should not be imposed. Similarly, in the event that applying the mask results in "icebergs" (i.e., regions of floating ice that are now detached from the main ice shelf), these floating regions should be removed.
144
145 The datasets can be obtained from the ISMIP6 2300 Forcing Globus endpoint AIS/Ice_Shelf_Fracture/ directory.
146
147 [[Image(CollapseMaskCcsm.png)]]
148
149 '''Figure 4''': Ice shelf collapse mask for CCSM4 under RCP8.5
150
151 ----
152 == Requirements and options for the projections ==
153 We encourage participants to contribute results using different models, grid resolutions, physics options, and/or initialization methods.
154
155 * Models must be able to prescribe a given SMB anomaly.
156 * Models must be able to use an ocean melt parameterization based on ocean thermal forcing evolving over time, such as the one used for the '''''Tier-1''''' experiments of the ISMIP6-Antarctica Projections.
157 * The adjustment of the SMB due to geometric ice-sheet changes in forward experiments is encouraged.
158 * Bedrock adjustment in forward experiments is allowed.
159 * The choice of model input data is unconstrained, to allow participants to use their preferred model setup without modification. Modelers without a preferred dataset can look at the ISMIP6 [https://theghub.org/dataset-listing Browse Datasets] page for possible options.
160 * Participants must submit a '''README''' file along with the model outputs as an integral part of the contribution to ISMIP6. The README template may be obtained [https://www.dropbox.com/s/6snnco54zo76s02/README_AIS_Projections.docx?dl=0|| here] or requested by email to ismip6@gmail.com. To allow for analysis, models must be well documented.
161
162 [[Image(ice_sheet_image_thin_2.png)]]
163 ----
164 ==Appendix 1 – Output grid definition and interpolation==
165 All '''2D data''' is requested on a regular grid with the following description: polar stereographic projection with standard parallel at 71° S and a central meridian of 0° W on datum WGS84. The lower left corner is at (-3,040,000 m, -3,040,000 m) and the upper right at (3,040,000 m, 3,040,000 m). This is the same grid used to provide the 8-km SMB and oceanic forcings. The output should be submitted at a resolution similar to the native model resolution and can be 32 km, 16 km, 8 km, 4 km, 2 km, or 1 km. The data will be stored at this resolution for archiving and conservatively interpolated to the 8-km grid for diagnostic processing by ISMIP6. Output should be provided with single precision.
166
167 If interpolation is required to transform the SMB forcing to your native grid, and to transform your model variables to a standard ISMIP6 output grid (32 km, 16 km, 8 km, 4 km, 2 km, or 1 km), conservative interpolation must be used. This requirement minimizes model-to-model differences due to the choice of interpolation method.
168
169 ===A1.1 Regridding Tools and Tips ===
170 An overview of the regridding process can be found on the following pages:
171 * [https://theghub.org/groups/ismip6/wiki/RegriddingwithCDO Regridding_with_CDO] contains tools and tips that have been used by ISMIP6 members.
172 * [https://theghub.org/groups/ismip6/wiki/RegriddingBISICLESoutputwithESMFandNCO Regridding with BISICLES with ESMF and NCO ]
173
174 ISMIP6 is designing tools to help with the regridding. If you need help with conservative interpolation, please email ismip6@gmail.com.
175
176 ----
177 ==Appendix 2 – Naming conventions, accessing, uploading and modeling output data==
178
179 '''Please provide:'''
180 * one variable per file for all '''2D fields''' and '''scalar variables'''
181 * a completed readme file
182 * single precision should be used for all output
183
184 ----
185 ===A2.1 File name convention===
186
187 File name convention for '''2D fields''' and '''scalar variables''':
188
189 `____.nc`
190
191 File name convention for '''readme file''':
192
193 `README___.doc`
194
195 '''where'''
196
197 `` = variable name (e.g. lithk)
198
199 `` = ice sheet (AIS or GIS)
200
201 `` = group acronym (all upper case or numbers, no special characters)
202
203 `` = model acronym (all upper case or numbers, no special characters)
204
205 `` = experiment name (from the experiment list, e.g. AE01)
206
207 For example, a file containing the variable "orog" for the Antarctic ice sheet, submitted by group “JPL” with model “ISSM” for experiment 'AE05' would be called:
208
209 `orog_AIS_JPL_ISSM_AE05.nc`
210
211 If JPL repeats the experiments with a different version of the model (for example, by changing the sliding law), the model could be named ISSM2, and so forth.
212
213 ----
214 ===A2.2 Accessing ISMIP6 datasets and accessing the forcing data ===
215
216 '''ISMIP6 datasets and directories''' are distributed via the Ghub Globus web application. Public datasets can be found in Ghub's [https://theghub.org/dataset-listing Browse Data] page. ISMIP6-specific initMIP Antarctic (and initMIP Greenland and projection data) can be accessed through the '''Ghub endpoints via Globus UI'''. To access and download data, one must create a Ghub account and register with Globus. Instructions to create accounts can be referenced in the '''2300 AIS ISMIP6 Globus Instructions (v. Feb 2023)''' instruction document at the end of this wiki.
217
218 This documentation provides instructions on how to use Globus to download Ghub data in general, including the ISMIP6 datasets distributed via Ghub. These datasets are from earlier ISMIP6 activities, such as the initMIP, ABUMIP or projections to 2100. ISMIP6 and GHub is partnered with UB CCR to provide access to large datasets. These datasets are described in detail on our [https://theghub.org/dataset-listing Browse Data] page.
219
220 If you have any questions or issues, please contact us by email at ismip6@gmail.com. See more details on Ghub's [https://theghub.org/accessing-data-with-globus Accessing Data] wiki. Please also check the suggested text to acknowledge the many scientists and organizations that made the ISMIP6 data possible.
221
222 ===A2.2.1 Where to upload your results===
223 Model results should also be uploaded on GHub. Once ready to upload your results, you should send an email to ismip6@gmail.com and ask that a new directory be created for your model results. You will then be able to upload your results in this directory using Globus.
224
225 ===A2.2.2 Reducing the size of files===
226 The size of the model files on higher-resolution grids can be reduced by file compression, which will save space on the storage server. Example commands are given below. Using these commands, we can get 10x compression, and for the masks even more given that contiguous masks contain repeated data. !NetCDF files have been designed with compression in mind. A !NetCDF file can be compressed without changing the way it is read into Matlab or Python (or any other language that uses standard !NetCDF read/write libraries).
227
228 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). The `‘-s’` option is the shuffling option to improve compression even more. We recommend using the `‘d1’` option, which seems to accomplish the desired compression.
229
230 Example of '''netcdf compression command''':
231
232 {{{
233 nccopy -d1 -s sftgif_GIS_JPL_ISSMPALEO_historical.nc sftgif_GIS_JPL_ISSMPALEO_historical_c.nc
234 }}}
235
236 Example of '''compression variant''', seems to work better for masks:
237
238 {{{
239 nccopy -d1 sftgif_GIS_JPL_ISSMPALEO_historical.nc sftgif_GIS_JPL_ISSMPALEO_historical_c.nc
240 }}}
241
242 ----
243 ===A2.3 Model output variables and README file===
244
245 The '''README''' file is an important contribution to the ISMIP6 submission. A template may be obtained [https://drive.google.com/open?id=1HOtbZhVUttXwsYPdJp0dsUqVZDcO6yN4&authuser=helene.l.seroussi%40dartmouth.edu&usp=drive_fs|| here] or requested by email to ismip6@gmail.com.
246
247
248 ===A2.3.1 General guidelines===
249
250 The variables requested in '''Table A1''' serve to evaluate and compare the different models and initialization techniques. Some variables may not apply to your model, in which case they can be omitted (with an explanation in the README file).
251
252 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 each year for both scalars and 2D variables (for initMIP, 2D variables were requested only over five-year periods), while flux variables should be averaged over the respective period. Please specify in your '''README''' file how your reported flux data has been averaged over time. Ideally, fluxes are averaged over all native model time steps.
253
254 Flux variables are defined as positive when the process adds mass to the ice sheet and negative otherwise.
255
256 All "missing data" must be assigned the single-precision floating point value of 1e^20^. Fields should be undefined outside the ice mask.
257
258
259 ===A2.3.2 How to record time in historical and projection files===
260
261 In compliance with CMIP6, time should be defined in "days since ", where `` 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 `` as used in the historical run from which it was initiated.
262
263 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]):
264 * `common_year` is 365 days,
265 * `leap_year` is 366 days
266 * `Julian_year` is 365.25 days,
267 * `Gregorian_year` is 365.2425 days,
268 * `360_day` year has 360 days divided into twelve 30-day months
269 * Please see the CF link above for other examples on calendar setting in section 4.4.
270
271 To illustrate a time recording for the historical file and projections for a typical state variable ('''ST''', e.g. thickness) and flux variable ('''FL''', e.g. SMB), we assume that our `` is January 1st 2013, and that we use a `calendar = 360_days`.
272
273 Other calendars can be used, but you need to indicate the calendar used in the netcdf file, 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,
274
275 '''''Note:'''''
276 * In CMIP a full year is typically from first of January to the first of January of the following year. We also provide below an example of what the `netcdf` would look like for our example.
277 * The end of the historical run should not be included at the beginning of the projections to avoid repetitions of the same values. This allows to merge together the different periods without having repeated entries. Models that do not have a historical period and start directly in January 2015 should still submit a historical run with one time step containing the model's initial conditions.
278
279
280 [[Image(ISMIP6_time2300_ISM.png)]]
281
282 ----
283 =For '''state variables''', like thickness for the historical:=
284
285 {{{
286 dimensions:
287 time = UNLIMITED ; // (3 currently)
288 variables:
289 double time(time) ;
290 time:units = "days since 1-1-2013" ; // This date correspond to the example basetime
291 time:calendar = "360_day" ; // Other calendars can be used... change here to relevant calendar
292 time:axis = "T" ;
293 time:long_name = "time" ;
294 time:standard_name = "time" ;
295 data:
296 time = 0, 360, 720; // If you use a different calendar these values will change
297 }}}
298
299 =For '''thickness for the projection'''=
300 '''''note:''''' that the full time entries are not shown, only beginning and end) would be:
301
302 {{{
303 dimensions:
304 time = UNLIMITED ;
305 variables:
306 double time(time) ;
307 time:units = "days since 1-1-2013" ; // This date correspond to the example basetime
308 time:calendar = "360_day" ; // Other calendars can be used... change here to relevant calendar
309 time:axis = "T" ;
310 time:long_name = "time" ;
311 time:standard_name = "time" ;
312 data:
313 time = 720, 1080, 1440, …, 103320, 103680; // If you use a different calendar these values will change
314 }}}
315
316 =The '''flux variable''', like SMB, would be recorded as the average over a full year, so for the historical: =
317
318 {{{
319 dimensions:
320 time = UNLIMITED ; // (2 currently)
321 bnds = 2 ;
322 variables:
323 double time(time) ;
324 time:bounds = "time_bnds" ;
325 time:units = "days since 1-1-2013 " ; // This date correspond to the example basetime
326 time:calendar = "360_day" ; // Other calendars can be used... change here to relevant calendar
327 time:axis = "T" ;
328 time:long_name = "time" ;
329 time:standard_name = "time" ;
330 double time_bnds(time, bnds) ;
331 data:
332 time = 180, 540 ; // If you use a different calendar these values will change.
333 //This is the middle of the time_bnds
334 time_bnds =
335 0, 360, //If you use a different calendar these values will change.
336 //These are the day since basetime at the beginning and end of the year
337 360, 720 ;
338 }}}
339
340
341 =For the '''projection=
342 '''''note:''''' that the full time entries are not shown, only beginning and end):
343
344 {{{
345 dimensions:
346 time = UNLIMITED ; //
347 bnds = 2 ;
348 variables:
349 double time(time) ;
350 time:bounds = "time_bnds" ;
351 time:units = "days since 1-1-2013 " ; // This date correspond to the example basetime
352 time:calendar = "360_day" ; // Other calendars can be used... change here to relevant calendar
353 time:axis = "T" ;
354 time:long_name = "time" ;
355 time:standard_name = "time" ;
356 double time_bnds(time, bnds) ;
357 variables:
358 double time(time) ;
359 time:bounds = "time_bnds" ;
360 data:
361 time = 900, 1260, 1620, ..., 103140, 103500; // If you use a different calendar these values will change.
362 //This is the middle of the time_bnds
363 time_bnds =
364 720, 1080, //If you use a different calendar these values will change.
365 //These are the day since basetime at the beginning and end of the year
366 1080, 1440,
367 1440, 1800,
368 ....
369 102960, 103320,
370 103320, 103680;
371 }}}
372
373
374 [[Image(ice_sheet_image_thin.png)]]
375
376 ===A2.3.3 Table A1: Variable request for ISMIP6===
377
378 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.
379
380 '''Model Characteristics'''
381 The Model Characteristics table can be found [https://docs.google.com/spreadsheets/d/1YEiPV3Uc0K8EqD57mz-ZqBIF4XnVb5UbyRDRmox0KeE/edit#gid=1534572380 here].
382
383
384 ||||||||||||||'''Table A1: Variable request for ISMIP6 projections.'''
385
386 '''Bold names''' or "alias" indicate a change compared to initMIP, to align the request with the CMIP6 official MIPtable "!IyrAnt" 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.
387
388 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. ||
389 ||'''Variable'''||'''Dim'''||'''Type'''||'''Variable Name '''||'''Standard Name '''||'''Units ''' ||'''Comment '''||
390 ||||||||||||||'''''2D variables requested yearly as snapshots (end of the year) for type ST and as yearly average for type FL.'''''||
391 ||Ice thickness||x,y,t||ST||lithk||land_ice_thickness||m||The thickness of the ice sheet||
392 ||Surface elevation||x,y,t||ST||orog||surface_altitude||m||The altitude or surface elevation of the ice sheet||
393 ||Base elevation||x,y,t||ST||base||base_altitude||m||The altitude of the lower ice surface elevation of the ice sheet||
394 ||Bedrock elevation||x,y,t||ST||topg||bedrock_altitude||m||The bedrock topography (may change during the projections)||
395 ||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 land ice interface (only needed beneath the grounded ice). If this quantity does not change with time, then a single entry is sufficient ||
396 ||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||
397 ||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)||
398 ||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)||
399 ||Ice thickness imbalance|| x,y,t||FL||dlithkdt||tendency_of_land_ice_thickness||m s^-1^||dHdt||
400 ||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||
401 ||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||
402 ||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||
403 ||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||
404 ||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||
405 ||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||
406 ||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||
407 ||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||
408 ||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||
409 ||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||
410 ||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||
411 ||Basal drag||x,y,t||ST||strbasemag||'''land_ice_basal_drag''' alias "magnitude_of_land_ice_basal_drag"||Pa||Basal drag||
412 ||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||
413 ||'''Ice front melt and calving 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 ice front melting and calving. Only for grid cells in contact with ocean||
414 ||Grounding line flux||x,y,t||FL||ligroundf|| land_ice_specific_mass_flux_at_grounding_line ||kg m^-2^ s^-1^||Loss of grounded ice mass resulting at grounding line. Only for grid cells in contact with grounding line||
415 ||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)||
416 ||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||
417 ||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||
418 ||||||||||||||'''''Scalar outputs requested every full year: snapshots for type ST and 1 year averages for type FL.'''''||
419 ||Total ice mass||t||ST||lim||land_ice_mass||kg||spatial integration, volume times density||
420 ||Mass above floatation||t||ST||limnsw||land_ice_mass_not_displacing_sea_water||kg||spatial integration, volume times density||
421 ||Grounded ice area||t||ST||'''iareagr''' alias "iareag"||'''grounded_ice_sheet_area''' alias "grounded_land_ice_area"||m^2^||spatial integration||
422 ||Floating ice area||t||ST||'''iareafl''' alias "iareaf"||floating_ice_shelf_area ||m^2^||spatial integration||
423 ||Total SMB flux||t||FL||tendacabf||tendency_of_land_ice_mass_due_to_surface_mass_balance||kg s^-1^||spatial integration||
424 ||Total BMB flux||t||FL||tendlibmassbf||tendency_of_land_ice_mass_due_to_basal_mass_balance||kg s-1||spatial integration||
425 ||'''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)'''||
426 ||Total calving flux||t||FL||tendlicalvf||tendency_of_land_ice_mass_due_to_calving||kg s^-1^||spatial integration||
427 ||'''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||
428 ||Total grounding line flux||t||FL||tendligroundf||tendency_of_grounded_ice_mass||kg s^-1^||spatial integration||
429
430 ----
431 ==Appendix 3 – Participating Models and Characteristics==
432
433 ===Antarctica Standalone Ice Sheet Modeling for 2300 Projections===
434
435 The list of participating models will be added as simulations are submitted.
436
437
438 ----
439 [[Image(ice_sheet_image_thin_2.png)]]
440
441 ==References ==
442
443 * '''Barthel, A.''', Agosta, C., Little, C.M., Hattermann, T., Jourdain, N.C., Goelzer, H., Nowicki, S., Seroussi, H., Straneo, F. and Bracegirdle, T.J. (2020): [https://www.the-cryosphere.net/14/855/2020/ CMIP5 model selection for ISMIP6 ice sheet model forcing: Greenland and Antarctica], The Cryosphere, 14(3), 855–879, https://doi.org/10.5194/tc-14-855-2020.
444 * '''Favier, L.''', Jourdain, N. C., Jenkins, A., Merino, N., Durand, G., Gagliardini, O., Gillet-Chaulet, F., and Mathiot, P. (2019): [https://www.geosci-model-dev-discuss.net/gmd-2019-26/ Assessment of Sub-Shelf Melting Parameterisations Using the Ocean-Ice Sheet Coupled Model NEMO(v3.6)-Elmer/Ice(v8.3)], Geosci. Model Dev., https://doi.org/10.5194/gmd-12-2255-2019.
445 * '''Jourdain''', N.C., Asay-Davis, X., Hattermann, T., Straneo, F., Seroussi, H., Little, C.M. and Nowicki, S. (2020): [https://tc.copernicus.org/articles/14/3111/2020/tc-14-3111-2020.pdf A protocol for calculating basal melt rates in the ISMIP6 Antarctic ice sheet projections]. The Cryosphere, 14(9), 3111-3134. https://doi.org/10.5194/tc-14-3111-2020.
446 * '''Nowicki, S.''', Goelzer, H., Seroussi, H., Payne, A. J., Lipscomb, W. H., Abe-Ouchi, A., Agosta, C., Alexander, P., Asay-Davis, X. S., Barthel, A., Bracegirdle, T. J., Cullather, R., Felikson, D., Fettweis, X., Gregory, J. M., Hattermann, T., Jourdain, N. C., Kuipers Munneke, P., Larour, E., Little, C. M., Morlighem, M., Nias, I., Shepherd, A., Simon, E., Slater, D., Smith, R. S., Straneo, F., Trusel, L. D., van den Broeke, M. R., and van de Wal, R. (2020): [https://tc.copernicus.org/articles/14/2331/2020/ Experimental protocol for sea level projections from ISMIP6 stand-alone ice sheet models], The Cryosphere, 14, 2331–2368, https://doi.org/10.5194/tc-14-2331-2020.
447 * '''Seroussi, H.''', Nowicki, S., Simon, E., Abe Ouchi, A., Albrecht, T., Brondex, J., Cornford, S., Dumas, C., Gillet-Chaulet, F., Goelzer, H., Golledge, N. R., Gregory, J. M., Greve, R., Hoffman, M. J., Humbert, A., Huybrechts, P., Kleiner, T., Larour, E., Leguy, G., Lipscomb, W. H., Lowry, D., Mengel, M., Morlighem, M., Pattyn, F., Payne, A. J., Pollard, D., Price, S., Quiquet, A., Reerink, T., Reese, R., Rodehacke, C. B., Schlegel, N.-J., Shepherd, A., Sun, S., Sutter, J., Van Breedam, J., van de Wal, R. S. W., Winkelmann, R., and Zhang, T. (2019): [https://www.the-cryosphere.net/13/1441/2019/ initMIP-Antarctica: An ice sheet model initialization experiment of ISMIP6], The Cryosphere., 13, 1441-1471, https://doi.org/10.5194/tc-13-1441-2019.
448 * '''Seroussi, H.''', Nowicki, S., Payne, A. J., Goelzer, H., Lipscomb, W. H., Abe-Ouchi, A., Agosta, C., Albrecht, T., Asay-Davis, X., Barthel, A., Calov, R., Cullather, R., Dumas, C., Galton-Fenzi, B. K., Gladstone, R., Golledge, N. R., Gregory, J. M., Greve, R., Hattermann, T., Hoffman, M. J., Humbert, A., Huybrechts, P., Jourdain, N. C., Kleiner, T., Larour, E., Leguy, G. R., Lowry, D. P., Little, C. M., Morlighem, M., Pattyn, F., Pelle, T., Price, S. F., Quiquet, A., Reese, R., Schlegel, N.-J., Shepherd, A., Simon, E., Smith, R. S., Straneo, F., Sun, S., Trusel, L. D., Van Breedam, J., van de Wal, R. S. W., Winkelmann, R., Zhao, C., Zhang, T., and Zwinger, T. (2020): [https://tc.copernicus.org/articles/14/3033/2020/ ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21^st^ century], The Cryosphere, 14, 3033–3070, https://doi.org/10.5194/tc-14-3033-2020.
449
450 ----
451 ==Acknowledgements==
452 The experimental protocol and datasets for the ISMIP6-Projections2300-Antarctica standalone ice sheet simulations would not have been possible without the effort of many scientists who have given their time and expertise, and have run models to convert the CMIP5 and CMIP6 model outputs into datasets that standalone ice sheet models can use.
453
454 ISMIP6 would like to thank the ocean focus group under the leadership of Fiamma Straneo, the atmospheric focus group under the leadership of William Lipscomb and Robin Smith, and the CMIP5 model evaluation focus group under the leadership of Alice Barthel. Xylar Asay-Davis, Nicolas Jourdain, Tore Hattermann, Chris Little, and Helene Seroussi were instrumental in the development of the ice shelf basal melt rate parameterization and associated datasets. Erika Simon, Richard Cullather and Sophie Nowicki prepared the atmospheric dataset. Luke Trusel and Helene Seroussi prepared the ice shelf fracture dataset. Alice Barthel, Chris Little, Cecile Agosta, Nicolas Jourdain, and Tore Hattermann 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.
455
456 Finally, we thank the ISMIP6 ice sheet modelers for their feedback on the design of the protocol and their willingness to participate in ISMIP6.
457
458 ----
459 =How to obtain the 2300 AIS ISMIP6 dataset=
460 -
The dataset for this suite of experiments is currently distributed upon request via Globus. Please see the word document below. The data will be made public and distributed via Ghub once the community paper has been published.
+
The dataset for this suite of experiments is currently distributed via Ghub/Globus in a manner similar to the other ISMIP6 datasets.
461 -
+
If you use the dataset, please reference the ISMIP6 protocol paper (Nowicki et al. 2020), the model selection (Barthel et al., 2020) and the Antarctic basal melt (Jourdain et al. 2020) and the community publication (Seroussi et al., 2024). The work from CMIP6 should also be recognized and we suggest text to include in any acknowledgement section in our publication page. Please also consider asking for an ISMIP6 contribution number, so that we can include your work on our publication page.
462 -
If you use the dataset, please reference the ISMIP6 protocol paper (Nowicki et al. 2020), the model selection (Barthel et al., 2020) and the Antarctic basal melt (Jourdain et al. 2020). As soon as the community publication is available, this should be cited too. The work from CMIP6 should also be recognized and we suggest text to include in any acknowledgement section in our publication page. Please also consider asking for an ISMIP6 contribution number, so that we can include your work on our publication page.
+
463 -
+
464 -
[[File(Globus_Instructions_2300AIS_Feb2023.docx)]]
+
465
466 [[Image(ice_sheet_image_thin.png)]]