Skip to main content
zenodoopen

Model data for probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under high ice-shelf basal melt conditions

<p>Model data for probabilistic projections of the Amery Ice Shelf catchment,&nbsp;<br>Antarctica, under high ice-shelf basal melt conditions<br>======================================================================</p> <p>This archive contains ice-sheet model output and statistical models&nbsp;<br>used in the manuscript:<br>"Probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under<br>high ice-shelf basal melt conditions" by:<br>Sanket Jantre, Matthew Hoffman, Nathan M. Urban, Trevor Hillebrand, Mauro<br>Perego, Stephen Price, John D. Jakeman</p> <p>OVERVIEW<br>--------</p> <p>Statistical models and corresponding datasets are archived in the file<br>Amery_UQ_Study_Statistical_Models.zip. &nbsp;That file contains a README that<br>describes its contents.<br>Contact for statistical models:<br>Sanket Jantre, Brookhaven National Laboratory, sjantre@bnl.gov</p> <p>The remaining files in this archive contain output from ice-sheet model&nbsp;<br>simulations using the MPAS-Albany Land Ice (MALI) model&nbsp;<br>(Hoffman et al. 2018, https://github.com/MALI-Dev/E3SM)<br>applied to a regional domain of the Amery Ice Shelf catchment<br>of Antarctica with mesh resolution varying from 4 to 20 km.<br>The contents of these files are described below.<br>Contact for MALI simulations:<br>Matthew Hoffman, Los Alamos National Laboratory, mhoffman@lanl.gov</p> <p>ENSEMBLES<br>---------</p> <p>The ensembles consist of all or a subset of 200 model runs with values for 6<br>parameters sampled from a Sobol' sequence over a specified range.</p> <p>Parameters and Sampled Ranges:<br>1. Ice stiffness scaling factor, C&phi;: (0.8, 1.2)<br>2. Basal friction scaling factor, C&mu;: (0.8, 1.2)<br>3. Basal slip exponent, q: (0.1, 0.333)<br>4. Calving yield stress, &sigma;max: (80, 180) kPa<br>5. Ice-shelf melt coefficient, &gamma;0: (9620, 471000) m yr&minus;1<br>6. Ice-shelf basal melt rate, m: (12, 58) Gt yr&minus;1</p> <p>SCENARIOS<br>---------</p> <p>The archive contains output from 4 scenarios described in the manuscript:</p> <p>Historical relaxation (RELX): For each ensemble member, we conducted a 50 year<br>relaxation from the initial condition using historical climate forcing to<br>integrate out fast transient behavior. For surface mass balance, we applied a<br>1995&ndash;2017 climatological average from RACMO2.3p1 (Van Wessem et al., 2014; van<br>den Broeke, 2019). The ocean thermal forcing was the observation-based<br>climatology compiled for ISMIP6-Antarctica, which uses data from 1995&ndash;2018<br>(Jourdain et al., 2020; Nowicki et al., 2020). The 50 year relaxation duration<br>was chosen as the most rapid adjustments occur in the first few decades of<br>integration, while the long term adjustment to a fully steady state takes<br>thousands of years. Relaxation to full steady state would require<br>substantially more computing resources than our entire set of ensembles and<br>also leads to the complication of different runs having potentially very<br>different initial states. Future improvements to model initialization that<br>account for surface elevation change (Perego et al., 2014) may reduce model<br>drift and adjust this requirement. The final model state in each run at the<br>end of RELX was given the nominal date of January 1, 2015, and all three<br>projection ensembles were branched from these states.</p> <p>Control projection (CTRL): The CTRL projection ensemble were an extension of<br>the RELX configurations, continuing the same surface mass balance and ocean<br>thermal forcing from January 1, 2015, to January 1, 2300. This ensemble was<br>used to assess model drift relative to the forced response of the climate<br>scenarios.&nbsp;</p> <p>SSP1-2.6 projection (SSP1): Our SSP1 projection used annual surface mass<br>balance and ocean thermal forcing derived from a UKESM SSP1-2.6 climate<br>scenario (expAE10 from Seroussi et al., 2024). Surface mass balance and ocean<br>thermal forcing were applied as anomalies relative to the climatological mean<br>forcings in RELX/CTRL to avoid issues related to climate model bias and abrupt<br>changes in forcing. This ensemble was also run from January 1, 2015, to<br>January 1, 2300.</p> <p>SSP5-8.5 projection (SSP5) Our SSP5 projection used the UKESM SSP5-8.5<br>projection forcings (expAE05 from Seroussi et al., 2024), again applied as<br>anomalies and from 2015 to 2300.</p> <p><br>OUTPUT<br>------</p> <p>Each ensemble directory contains at the base level:<br>* branch_ensemble.cfg - configuration file for the ensemble_generator test<br>&nbsp; case (https://mpas-dev.github.io/compass/latest/users_guide/landice/test_groups/ensemble_generator.html)&nbsp;<br>&nbsp; of COMPASS (Configuration Of Model for Prediction Across Scales Setups, https://github.com/MPAS-Dev/compass)<br>&nbsp; used to set up the ensemble.<br>* mesh_vars.nc - netCDF file containing the mesh variables. &nbsp;This file is the<br>&nbsp; same for all ensembles. &nbsp;These fields need to be appended to an output file to<br>&nbsp; visualize in, e.g., Paraview. &nbsp;This can be done with the command:<br>&nbsp; ncks -A mesh_vars.nc output.nc</p> <p>Each ensemble contains subdirectories for each run numbered 000-199. &nbsp;Each run<br>directory contains the files:<br>* globalStats.nc - scalar metrics at every time step<br>* output.nc - spatial fields at 10 year intervals<br>* run_info.cfg - summary of parameter values</p> <p>Notes:<br>* The SSP1 and SSP5 ensembles contain a subset of the full 200 runs<br>&nbsp; because some runs were filtered out as being unnecessary during model<br>&nbsp; calibration (see manuscript).<br>* The RELX ensemble is run for 50 years with an arbitrary start year of 2000.<br>&nbsp; The CTRL, SSP1, and SSP5 ensembles are started in 2015 from the final state<br>&nbsp; of the corresponding RELX run (nominally 2050). &nbsp;Because spatial data has<br>&nbsp; been saved in this archive at 10 year intervals, the first output for the<br>&nbsp; three projection ensembles is 2020. &nbsp;To get the initial condition for the<br>&nbsp; projection ensembles, use the 2050 state from RELX for the corresponding<br>&nbsp; model run.<br>* Spatial field output at 1 year intervals, as well as run input, log, and<br>&nbsp; restart files are available from the corresponding contact. &nbsp;The complete<br>&nbsp; run data is about 2 TB.</p> <p><br>CITATION<br>--------<br>If you find these data useful, please cite the DOI for this&nbsp;<br>archive (10.5281/zenodo.11166628), as well as our manuscript<br>submitted to The Cryosphere.</p> <p>This archive is approved by Los Alamos National Laboratory<br>for public release under LA-UR-24-24969; distribution is unlimited.</p> <p>&nbsp;</p> <p>REFERENCES<br>----------</p> <p>Hoffman, M. J., Perego, M., Price, S. F., Lipscomb, W. H., Zhang, T.,<br>Jacobsen, D., Tezaur, I., Salinger, A. G., Tuminaro, R., and Bertagna, L.:<br>MPAS-Albany Land Ice (MALI): a variable-resolution ice sheet model for Earth<br>system modeling using Voronoi grids, Geoscientific Model Development, 11,<br>3747&ndash;3780, 2018.</p> <p>Jourdain, N. C., Asay-Davis, X., Hattermann, T., Straneo, F., Seroussi, H.,<br>Little, C. M., and Nowicki, S.: A protocol for calculating basal melt rates in<br>the ISMIP6 Antarctic ice sheet projections, The Cryosphere, 14, 3111&ndash;3134,<br>https://doi.org/10.5194/tc-14-3111-2020, 2020.</p> <p>Nowicki, S., Goelzer, H., Seroussi, H., Payne, A. J., Lipscomb, W. H.,<br>Abe-Ouchi, A., Agosta, C., Alexander, P., Asay-Davis, X. S., Barthel, A.,<br>Bracegirdle, T. J., Cullather, R., Felikson, D., Fettweis, X., Gregory, J. M.,<br>Hattermann, T., Jourdain, N. C., Kuipers Munneke, P., Larour, E., Little, C.<br>M., Morlighem, M., Nias, I., Shepherd, A., Simon, E., Slater, D., Smith, R.<br>S., Straneo, F., Trusel, L. D., van den Broeke, M. R., and van de Wal, R.:<br>Experimental protocol for sea level projections from ISMIP6 stand-alone ice<br>sheet models, The Cryosphere, 14, 2331&ndash;2368,<br>https://doi.org/10.5194/tc-14-2331-2020, 2020.</p> <p>Perego, M., Price, S., and Stadler, G.: Optimal Initial Conditions for<br>Coupling Ice Sheet Models to Earth System Models, Journal of Geophysical<br>Research: Earth Surface, 119, 1894&ndash;1917, https://doi.org/10.1002/2014JF003181,<br>2014.</p> <p>Seroussi, H., et al. 2024. ISMIP6 Projections 2300 Antarctica Protocol.<br>https://theghub.org/groups/ismip6/wiki/ISMIP6-Projections2300-Antarctica.&nbsp;</p> <p>van den Broeke, M.: RACMO2.3p1 annual surface mass balance Antarctica<br>(1979-2014), PANGAEA - Data Publisher for Earth &amp; Environmental Science,<br>https://doi.org/10.1594/PANGAEA.896940, 2019.</p> <p>van Wessem, J., Reijmer, C., Morlighem, M., Mouginot, J., Rignot, E., Medley,<br>B., Joughin, I., Wouters, B., Depoorter, M., Bamber, J., Lenaerts, J., Van De<br>Berg, W., Van Den Broeke, M., and Van Meijgaard, E.: Improved Representation<br>of East Antarctic Surface Mass Balance in a Regional Atmospheric Climate<br>Model, Journal of Glaciology, 60, 761&ndash;770,<br>https://doi.org/10.3189/2014JoG14J051, 2014.</p>

ShareScore

36/100

Overall dataset sharing score

Score breakdown

These five areas show where the dataset supports — or may limit — practical reuse.

Stewardship
8
Harmonization
4
Access
16
Reuse readiness
8
Engagement
0

Topics