Description: Lithospheric plates move over the low‐viscosity asthenosphere balancing several forces, which generate plate motions. We use a global 3‐D lithosphere‐asthenosphere model (SLIM3D) with visco‐elasto‐plastic rheology coupled to a spectral model of mantle flow at 300 km depth to quantify the influence of intra‐plate friction and asthenospheric viscosity on plate velocities. We account for the brittle‐ductile deformation at plate boundaries (yield stress) using a plate boundary friction coefficient to predict the present‐day plate motion and net rotation of the lithospheric plates. Previous modeling studies have suggested that small friction coefficients ( urn:x-wiley:15252027:media:ggge21498:ggge21498-math-0001, yield stress urn:x-wiley:15252027:media:ggge21498:ggge21498-math-0002 MPa) can lead to plate tectonics in models of mantle convection. Here we show that in order to match the observed present‐day plate motion and net rotation, the frictional parameter must be less than 0.05. We obtain a good fit with the magnitude and orientation of the observed plate velocities (NUVEL‐1A) in a no‐net‐rotation (NNR) reference frame with urn:x-wiley:15252027:media:ggge21498:ggge21498-math-0003 and a minimum asthenosphere viscosity of urn:x-wiley:15252027:media:ggge21498:ggge21498-math-0004 Pas to 1020 Pas. Our estimates of net rotation (NR) of the lithosphere suggest that amplitudes urn:x-wiley:15252027:media:ggge21498:ggge21498-math-0005 ( urn:x-wiley:15252027:media:ggge21498:ggge21498-math-0006/Ma), similar to most observation‐based estimates, can be obtained with asthenosphere viscosity cutoff values of urn:x-wiley:15252027:media:ggge21498:ggge21498-math-0007 Pas to urn:x-wiley:15252027:media:ggge21498:ggge21498-math-0008 Pas and friction coefficients urn:x-wiley:15252027:media:ggge21498:ggge21498-math-0009.

Description: The timing of past glaciation across the Tian Shan provides a proxy for past climate change in this critical area. Correlating glacial stages across the region is difficult but cosmogenic exposure ages have considerable potential. A drawback is the large observed scatter in 10Be surface exposure data. To quantify the robustness of the dating, we compile, recalculate, and perform statistical analyses on sets of 10Be surface exposure ages from 25 moraines, consisting of 114 new and previously published ages. We assess boulder age scatter by dividing boulder groups into quality classes and rejecting boulder groups of poor quality. This allows us to distinguish and correlate robustly dated glacier limits, resulting in a more conservative chronology than advanced in previous publications. Our analysis shows that only one regional glacial stage can be reliably correlated across the Tian Shan, with glacier expansions occurring between 15 and 28 ka during marine oxygen isotope stage (MIS) 2. However, there are examples of older more extensive indicators of glacial stages between MIS 3 and MIS 6. Paleoglacier extent during MIS 2 was mainly restricted to valley glaciation. Local deviations occur: in the central Kyrgyz Tian Shan paleoglaciers were more extensive and we propose that the topographic context explains this pattern. Correlation between glacial stages prior to late MIS 2 is less reliable, because of the low number of samples and/or the poor resolution of the dating. With the current resolution and spatial coverage of robustly-dated glacier limits we advise that paleoclimatic implications for the Tian Shan glacial chronology beyond MIS 2 are speculative and that continued work toward robust glacial chronologies is needed to resolve questions regarding drivers of past glaciation in the Tian Shan and Central Asia.

Description: Numerical ice sheet models constrained by theory and refined by comparisons with observational data are a central component of work to address the interactions between the cryosphere and changing climate, at a wide range of scales. Although there continue to be significant advances in modelling, major challenges still exist, in particular in terms of downscaling global climate model output to estimate regional and local climate patterns that are critical controls for the dynamics of glaciers and ice sheets. Ice sheet models are tested and refined by comparing model predictions of past ice geometries with field-based reconstructions from geological, geomorphological, and ice core data. However, on the East Antarctic Ice sheet, there is a critical gap in the empirical data required to reconstruct changes in ice sheet geometry in the Dronning Maud Land (DML) region. In addition, there is poor control on the regional climate history of the ice sheet margin, because ice core locations, where detailed reconstructions of climate history exist, are located on high inland domes. This leaves numerical models of regional glaciation history in this near-coastal area largely unconstrained. MAGIC-DML is an ongoing Swedish-US-Norwegian-German-UK collaboration with a focus on improving ice sheet models by combining advances in modeling with filling critical data gaps that exist in our knowledge of the timing and pattern of ice surface changes on the western Dronning Maud Land margin. A combination of geomorphological mapping using remote sensing data, field investigations, cosmogenic nuclide surface exposure dating, and numerical ice-sheet modelling are being used in an iterative manner to produce a comprehensive reconstruction of the glacial history of western Dronning Maud Land. We present an overview of the project, as well as results of the initial mapping and modelling that has been used to identify high potential sites for field sampling in 2016/17 and 2017/18.

Description: Given current concern about the stability of ice sheets, and potential sea level rise, it is imperative that we are able to reconstruct and predict the response of ice sheets to climate change. The Intergovernmental Panel on Climate Change (IPCC), amongst others, have highlighted that our current ability to do so is limited. Numerical ice sheet models are a central component of the work to address this challenge. An unresolved key issue in this work concerns the volume and rate of ice mass loss needed to explain the large difference between late glacial and interglacial global sea levels. Some 20% of observed sea level rise since the Last Glacial Maximum (LGM) cannot be attributed to any known former ice mass, indicating that this inconsistency arises from the deficiencies in modelled reconstructions of ice sheet volumes and postglacial rebound. Ice sheet models are tested and refined by comparing model predictions of past ice geometries with field-based reconstructions from geological, geomorphological and ice core data. However, on the East Antarctic Ice sheet, Dronning Maud Land (DML) presents a critical gap in the empirical data required to reconstruct changes in ice sheet geometry. In addition, there is poor control on regional climate histories of ice sheet margins, because ice core locations, where detailed reconstructions of climate history exist, are located on high inland domes. This leaves numerical models of regional glaciation history largely unconstrained. MAGIC-DML is a Swedish-US-Norwegian-German-UK collaboration with a focus on filling the critical data gaps that exist in our knowledge of the timing and pattern of ice surface changes on the western Dronning Maud Land margin. Here we describe a series of high-resolution modelling experiments to help identify those areas across western Dronning Maud Land that are the most sensitive to uncertainties in the regional climate history and the choice of model parameters. For this we employ a wide range of climate and ocean histories combining published outputs of 18 general circulation models for the LGM and mid-Holocene with ice core records. The modelling results together with remote sensing mapping of glacial landforms is informing and guiding cosmogenic nuclide sampling campaigns in western Dronning Maud Land starting 2016/17. Successful integration of numerical modelling and field investigations in an iterative manner is key to achieving the anticipated outcome of the MAGIC-DML project, a reconstruction of the long-term pattern and timing of vertical changes in ice surface elevation since the mid-Pliocene warm period, which will provide the missing empirical data required to constrain numerical ice sheet models.

Description: In the context of future climate change, understanding the nature and behaviour of ice sheets during warm intervals in Earth history is of fundamental importance. The late Pliocene warm period (also known as the PRISM interval: 3.264 to 3.025 million years before present) can serve as a potential analogue for projected future climates. Although Pliocene ice locations and extents are still poorly constrained, a significant contribution to sea-level rise should be expected from both the Greenland ice sheet and the West and East Antarctic ice sheets based on palaeo sea-level reconstructions. Here, we present results from simulations of the Antarctic ice sheet by means of an international Pliocene Ice Sheet Modeling Intercomparison Project (PLISMIP-ANT). For the experiments, ice-sheet models including the shallow ice and shelf approximations have been used to simulate the complete Antarctic domain (including grounded and floating ice). We compare the performance of six existing numerical ice-sheet models in simulating modern control and Pliocene ice sheets by a suite of five sensitivity experiments. We include an overview of the different ice-sheet models used and how specific model configurations influence the resulting Pliocene Antarctic ice sheet. The six ice-sheet models simulate a comparable present-day ice sheet, considering the models are set up with their own parameter settings. For the Pliocene, the results demonstrate the difficulty of all six models used here to simulate a significant retreat or re-advance of the East Antarctic ice grounding line, which is thought to have happened during the Pliocene for the Wilkes and Aurora basins. The specific sea-level contribution of the Antarctic ice sheet at this point cannot be conclusively determined, whereas improved grounding line physics could be essential for a correct representation of the migration of the grounding-line of the Antarctic ice sheet during the Pliocene.

Description: 3000 m of ice sheet thickness has ensured that central Greenland has kept it geothermal heat flow (GHF) distribution enigmatic. Some few direct ice temperature measurements from deep ice cores reveal a GHF of 50 to 60 mW/m2 in the Summit region and this is noticeably above what would be expected for the underlying Early Proterozoic lithosphere. In addition, indirect estimates from zones of rapid basal melting suggest extreme anomalies 15 to 30 times continental background. Subglacial topography indicates caldera like topographic features in the zones hinting at possible volcanic activity in the past [1], and all of these observations combined hint at an anomalous lithospheric structure. Further supporting this comes from new high-resolution P-wave tomography, which shows a strong thermal anomaly in the lithosphere crossing Greenland from east to west [2]. Rock outcrops at the eastern and western end of this zone indicate significant former magmatic activity, older in the east and younger in the west. Additionally, plate modelling studies suggest that the Greenland plate passed over the mantle plume that is currently under Iceland from late Cretaceous to Neogene times, consistent with the evidence from age of magmatism. Evidence of rapid basal melt revealed by ice penetrating radar along the hypocentre of the putative plume track indicates that it continues to affect the Greenland continental geotherm today. We analyse plume-induced thermal disturbance of the present-day lithosphere and their effects on the central Greenland ice sheet by using a novel evolutionary model of the climate-ice-lithosphere-upper mantle system. Our results indicate that mantle plume-induced erosion of the lithosphere has occurred, explaining caldera-type volcanic structures, the GHF anomaly, and requiring dyke intrusion into the crust during the early Cenozoic. The residual thermo-mechanical effect of the mantle plume has raised deep-sourced heat flow by over 25 mW/m2 since 60 Ma and explains the high basal melting rates of the Greenland ice sheet observed in the study area. [1] Fahnestock, M., Abdalati, W., Joughin, I., Brozena, J., Gogineni, P., 2001. High geothermal heat flow, Basal melt, and the origin of rapid ice flow in central Greenland. Science (New York, N.Y.). 294, 2338–2342. [2] Jakovlev, A.V., Bushenkova, N.A., Koulakov, I.Y., Dobretsov, N.L., 2012. Structure of the upper mantle in the Circum-Arctic region from regional seismic tomography. Russian Geology and Geophysics. 53, 963–971.

Description: Heat entering large ice sheets and glaciers from the underlying bedrock controls the thermal regime of the basal ice layers in terms of basal melting, the formation of basal temperate ice layers and its influence on sliding processes. In this study, we estimate the effect of the uncertainties in geothermal heat flux (GHF) on the topography, dynamics and basal thermal conditions of simulations of the present-day Greenland Ice Sheet (GIS), using the dynamic polythermal ice-sheet model SICOPOLIS. Today’s GHF at the base of the GIS has been measured at only the locations of the deepest ice cores. It should therefore be considered as highly uncertain forcing component for the dynamic modeling of the major part of the GIS. Three GHF models are applied as lithosphere boundary conditions driving a series of paleoclimatic simulations. The first is based on a 3-D global seismic model of the crust and upper mantle, constrained in regions where heat-flow measurements are available. The other two models are based on tectonic regionalization and a combination of tectonic regionalization and deep ice-core measurements. Although the tectonic regionalization model was originally taken as a starting point for both the seismic mantle/crust and tectonic/ice-core GHF models, the resulting heat flow distributions across the extent of the GIS differ significantly, not only by the mean values but also by their spatial patterns. As a consequence, the present-day GIS driven by the seismic GHF model is 450 meters thicker in the northern and central parts of the GIS, and 450 meters thinner in the southern part relative to the GIS forced by the tectonic/ice-core GHF model. In addition, the horizontal ice transport simulated using the seismic GHF model is 10-50% higher in the northern and western GIS, and 10-30% lower in the southern and eastern GIS than the transport resulting from the simulation driven by the tectonic/ice-core GHF model. These alterations in the GIS topography and horizontal transport induced by the two latter simulations with respect to the simulation driven by the tectonic regionalization model have the opposite tendencies. Finally, we discuss different scenarios of the evolution of the temperate base in the ice-sheet simulations driven by the seismic-based, tectonic regionalization and tectonic/ice-core GHF models. The modeled present-day temperate ice areas with non-zero basal melting cover 34%, 40.6% and 56.8% of the ice-covered area, respectively, for each of the GHF models examined. The minimum coverage of the temperate ice is reached shortly after the last glacial maximum and amounts to less than 20% of the ice-covered area, as determined using the seismic and tectonic GHF models, and about 40% when employing the tectonic/ice-core model. Although the response of the modeled ice-sheet to changing basal conditions may be influenced by the limitations in the classical polythermal approach used in SICOPOLIS (neglecting water diffusivity), we may conclude that the modeled present-day thermodynamical state of the GIS is largely sensitive to the GHF model employed

Description: Since the International Geophysical Year (1957), a view has prevailed that East Antarctica has a relatively homogeneous lithospheric structure, consisting of a craton-like mosaic of Precambrian terranes, stable since the Pan-African orogeny 500 million years ago (e.g. Ferracioli et al. 2011). Recent recognition of a continental-scale rift system cutting the East Antarctic interior has crystallised an alternative view of much more recent geological activity with important implications. The newly defined East Antarctic Rift System (EARS) (Ferraccioli et al. 2011) appears to extend from at least the South Pole to the continental margin at the Lambert Rift, a distance of 2500 km. This is comparable in scale to the well-studied East African rift system. New analysis of RadarSat data by Golynsky & Golynsky (2009) indicates that further rift zones may form widely distributed extension zones within the continent. A pilot study (Vaughan et al. 2012), using a newly developed gravity inversion technique (Chappell & Kusznir 2008) with existing public domain satellite data, shows distinct crustal thickness provinces with overall high average thickness separated by thinner, possibly rifted, crust. Understanding the nature of crustal thickness in East Antarctica is critical because: 1) this is poorly known along the ocean–continent transition, but is necessary to improve the plate reconstruction fit between Antarctica, Australia and India in Gondwana, which will also better define how and when these continents separated; 2) lateral variation in crustal thickness can be used to test supercontinent reconstructions and assess the effects of crystalline basement architecture and mechanical properties on rifting; 3) rift zone trajectories through East Antarctica will define the geometry of zones of crustal and lithospheric thinning at plate-scale; 4) it is not clear why or when the crust of East Antarctica became so thick and elevated, but knowing this can be used to test models of Cenozoic ice sheet formation and stability. References Chappell, A.R. & Kusznir, N.J. 2008. Three-dimensional gravity inversion for Moho depth at rifted continental margins incorporating a lithosphere thermal gravity anomaly correction. Geophysical Journal International, 174 (1), 1–13. Ferraccioli, F., Finn, C.A., Jordan, T.A., Bell, R.E., Anderson, L.M. & Damaske, D. 2011. East Antarctic rifting triggers uplift of the Gamburtsev Mountains Nature, 479, 388–392. Golynsky, A.V. & Golynsky, D.A. 2009. Rifts in the tectonic structure of East Antarctica (in Russian). Russian Earth Science Research in Antarctica, 2, 132–162. Vaughan, A.P.M., Kusznir, N.J., Ferraccioli, F. & Jordan, T.A.R.M. 2012. Regional heat-[U+FB02]ow prediction for Antarctica using gravity inversion mapping of crustal thickness and lithosphere thinning. Geophysical Research Abstracts, 14, EGU2012–8095.