Cordilleran ice sheet improved bedrock simulations continuous variables

These data contain a subset of time-dependent glacier model output variables. The ghf70 data files are an update on the reference below, fixing significant problems affecting the computation of the bedrock deformation in response to ice load (PISM Github issues #370 and #377) and the computation of ice temperature (PISM Github issue #371). The other data files additionally include spatially-variable geothermal heat flux (dav13, gou11comb, gou11simi, sha04), different lithospheric rigidity (eet30km) or mantle viscosity (num1e21), and higher horizontal resolution (3km).

Reference:

• Seguinot, J., Rogozhina, I., Stroeven, A. P.,  Margold, M. and Kleman, J.: Numerical simulations of the Cordilleran ice sheet through the last glacial cycle, The Cryosphere, 10(2), 639–664, doi:10.5194/tc-10-639-2016, 2016.

File names:

cisbed.{res}.{forcing}.{ex.100a|ts.10a}.{ghf}.{props}.nc

• Horizontal resolution:
  • 10km: 10 km horizontal resolution
  • 5km: 5 km horizontal resolution
  • 3km: 3 km horizontal resolution
• Temperature forcing:
  • epica: EPICA ice core temperature forcing
  • grip: GRIP ice core temperature forcing
• Variable types:
  • ex.100a: spatial diagnostics every hundred years
  • ts.10a: scalar time-series every ten years
• Geothermal heat flow:
  • ghf70: constant 70 mW m-2 heat flow
  • dav13: Davies (2013) geothermal heat flow map
  • gou11comb: Goutorbe et al. (2011) best combination method
  • gou11simi: Goutorbe et al. (2011) similarity method
  • sha04: Shapiro and Ritzwoller (2004) heat flow map
• Bedrock properties
  • eet30km: lithosphere elastic thickness of 30 km
  • num1e21: astenosphere viscosity of 1e21 Pa s

Data format:

The data use compressed netCDF format. For quick inspection I recommend ncview. Spatial diagnostics (*.ex.100a.nc) can be converted to GeoTIFF (and other GIS formats) e.g. using GDAL:

gdal_translate NETCDF:filename.nc:variable -b band filename.variable.band.tif

The list of variables (subdatasets) can be obtained from ncdump or gdalinfo. Band information can be displayed with:

gdalinfo NETCDF:filename.nc:variable

Variable long names, units, PISM configuration parametres and additional information are contained within the netCDF metadata.

Funding:

Swiss National Supercomputing Centre (CSCS) grants s573 and sm13 to J. Seguinot, Swiss National Science Foundation (SNSF) grants no.~200020-169558 and 200021-153179/1 to M. Funk, and Research Foundation – Flanders (FWO) Odysseus Type II project G0DCA23N 'GlaciersMD' to H. Zekollari.

Changelog:

• Version 1:
  • Initial version.