Here we show a brief selection of animations, plots, and illustrations directly related to our research topics.

  • This animation is part of an article published in Nature Climate Change (Robinson et al. 2012). The results of the simulations carried out in this work have been reflected in a video that has obtained the first prize of Scientific Dissemination in the category of videos of the Campus of International Excellence CEI Moncloa:
Robinson, Alexander; Calov, Reinhard; Ganopolski, Andrey. 2012. Multistability and critical thresholds of the Greenland ice sheet Nature Climate Change. Nature Publishing Group. 2-4, pp.429-432.
  • Animation of the Laurentide and Greenland ice sheets. Hudson Bay/Strait ice stream velocities are illustrated in colors. The time series at bottom show the evolution of the GRIP d18O core, the subsurface oceanic index, and the simulated calving of icebergs into the Labrador Sea.
Alvarez-Solas, J., Robinson, A., Montoya, M. and Ritz, C., 2013. Iceberg discharges of the last glacial period driven by oceanic circulation changes. Proceedings of the National Academy of Sciences, 110(41), pp.16350-16354.
  • Difference in Greenland Ice Sheet response to millennial-scale and orbital-only oceanic variability. Basal melt at the marine margins (upper side figure; millennial-scale signal in blue, orbital-only signal in dashed black); millennial minus orbital-only ice thickness (lower left side figure; green for mass gained, brown for mass loss); millennial minus orbital-only ice speed (lower right side figure; red for a speed up for the run forced with millennial-scale variability, blue for a speed up for the orbital-only run). Grounding line positions are also shown in Greenland maps (millennial in solid black, orbital-only in dashed black).
Tabone et al., 2019: Impact of millennial-scale oceanic variability on the Greenland ice-sheet evolution throughout the last glacial period, CP, 15, 593–609, 2019,
  • Simulated response of the Antarctic Ice Sheet to glacial millennial scale variability for a high oceanic sensitivity (kappa = 15 m/yr K). The left panel represents the surface elevation in m; the right panel shows the surface velocities in m/yr. The middle panel shows the applied oceanic millennial forcing (extracted from Dome C) in K; the bottom panel the total grounded ice volume in 10^6 km^3. No refreezing is allowed.
Blasco, J., Tabone, I., Alvarez-Solas, J., Robinson, A., and Montoya, M.: The Antarctic Ice Sheet response to glacial millennial-scale variability, Clim. Past, 15, 121–133,, 2019.
  • The NEGIS past response to submarine melt for the last glacial period. NEGIS grounding line distance from its present position (upper figure) simulated by GRISLI-UCM (dashed red for the run with no submarine melt applied at the marine margins, solid blue for the run forced by submarine melt with peaks of 3 m/yr at around 45 kyr ago) compared to observations (dashed black curve). Ice velocities at the Northeast Greenland (lower figure) are shown respectively in the left hand side and right hand side figures. Solid black lines on maps stand for the NEGIS simulated grounding line, while dashed black lines represent its margin position reconstructed for the last glacial maximum.
Tabone et al., 2019: Submarine melt as a potential trigger of the North East Greenland Ice Stream margin retreat during Marine Isotope Stage 3, TC, 13, 1911–1923, 2019,