Description: Spanning the northern sector of High Asia, the Altai region contains a rich landform record of glaciation. We report the extent, chronologies, and dynamics of two paleoglaciers on opposite flanks of the Ikh Turgen mountains (In Russian: Chikhacheva Range), straddling the border between Russia and Mongolia, using a combination of remote sensing-based glacial geomorphological mapping, 10Be surface exposure dating, and geomorphometric analysis. On the eastern side (Mongolia), the Turgen-Asgat paleoglacier, with its potential for developing a large accumulation area (∼257 km2), expanded 40 km down valley, and mean ages from a latero-frontal moraine indicate deglaciation during marine oxygen isotope stage (MIS) 3 (45.1 ± 1.8 ka, n = 4) and MIS 2 (22.8 ± 3.3 ka, n = 5). These minimum age constraints are consistent with other 10Be glacial chronologies and paleoclimate records from the region, which indicates glacier culmination during cold and wet conditions coinciding with MIS 3 (piedmont-style glaciation; inferred for a few sites across the region) and glacier culmination during cold and dry conditions coinciding with MIS 2 (mainly valley-style glaciation; inferred from several sites across the region). On the western side (Russia), the Boguty paleoglacier had a smaller accumulation area (∼222 km2), and advanced 30 km down valley across a low gradient forefield. Surface exposure ages from two moraine complexes on this side of the mountains exhibit wide scatter (∼14–53 ka, n = 8), making paleoclimate inferences and comparison to other proxies difficult. Ice surface profile reconstructions imply that the two paleoglaciers likely shared an ice divide.

Description: The shallow ice approximation (SIA) is commonly used in ice-sheet models to simplify the force balance equations within the ice. However, the SIA cannot adequately reproduce the dynamics of the fast flowing ice streams usually found at the margins of ice sheets. To overcome this limitation, recent studies have introduced heuristic hybrid combinations of the SIA and the shelfy stream approximation. Here, we implement four different hybrid schemes into a model of the Antarctic Ice Sheet in order to compare their performance under present-day conditions. For each scheme, the model is calibrated using an iterative technique to infer the spatial variability in basal sliding parameters. Model results are validated against topographic and velocity data. Our analysis shows that the iterative technique compensates for the differences between the schemes, producing similar ice-sheet configurations through quantitatively different results of the sliding coefficient calibration. Despite this we observe a robust agreement in the reconstructed patterns of basal sliding parameters. We exchange the calibrated sliding parameter distributions between the schemes to demonstrate that the results of the model calibration cannot be straightforwardly transferred to models based on different approximations of ice dynamics. However, easily adaptable calibration techniques for the potential distribution of basal sliding coefficients can be implemented into ice models to overcome such incompatibility, as shown in this study

Description: The orientation and tectonic regime of the observed crustal/lithospheric stress field contribute to our knowledge of different deformation processes occurring within the Earth's crust and lithosphere. In this study, we analyze the influence of the thermal and density structure of the upper mantle on the lithospheric stress field and topography. We use a 3-D lithosphere–asthenosphere numerical model with power-law rheology, coupled to a spectral mantle flow code at 300 km depth. Our results are validated against the World Stress Map 2016 (WSM2016) and the observation-based residual topography. We derive the upper mantle thermal structure from either a heat flow model combined with a seafloor age model (TM1) or a global S-wave velocity model (TM2). We show that lateral density heterogeneities in the upper 300 km have a limited influence on the modeled horizontal stress field as opposed to the resulting dynamic topography that appears more sensitive to such heterogeneities. The modeled stress field directions, using only the mantle heterogeneities below 300 km, are not perturbed much when the effects of lithosphere and crust above 300 km are added. In contrast, modeled stress magnitudes and dynamic topography are to a greater extent controlled by the upper mantle density structure. After correction for the chemical depletion of continents, the TM2 model leads to a much better fit with the observed residual topography giving a good correlation of 0.51 in continents, but this correction leads to no significant improvement of the fit between the WSM2016 and the resulting lithosphere stresses. In continental regions with abundant heat flow data, TM1 results in relatively small angular misfits. For example, in western Europe the misfit between the modeled and observation-based stress is 18.3°. Our findings emphasize that the relative contributions coming from shallow and deep mantle dynamic forces are quite different for the lithospheric stress field and dynamic topography.

Description: Ice-shelf basal melting is the largest contributor to the negative mass balance of the Antarctic ice sheet. However, current implementations of ice/ocean interactions in ice-sheet models disagree with the distribution of sub-shelf melt and freezing rates revealed by recent observational studies. Here we present a novel combination of a continental-scale ice flow model and a calibration technique to derive the spatial distribution of basal melting and freezing rates for the whole Antarctic ice-shelf system. The modelled ice-sheet equilibrium state is evaluated against topographic and velocity observations. Our high-resolution (10-km spacing) simulation predicts an equilibrium ice-shelf basal mass balance of −1648.7 Gt a−1 that increases to −1917.0 Gt a−1 when the observed ice-shelf thinning rates are taken into account. Our estimates reproduce the complexity of the basal mass balance of Antarctic ice shelves, providing a reference for parameterisations of sub-shelf ocean/ice interactions in continental ice-sheet models. We perform a sensitivity analysis to assess the effects of variations in the model set-up, showing that the retrieved estimates of basal melting and freezing rates are largely insensitive to changes in the internal model parameters, but respond strongly to a reduction of model resolution and the uncertainty in the input datasets.

Description: Climate change in Siberia is currently receiving a lot of attention because large permafrost-covered areas could provide a strong positive feedback to global warming through the release of carbon that has been sequestered there on glacial–interglacial timescales. Geological evidence and climate model experiments show that the Siberian region also played an exceptional role during glacial periods. The region that is currently known for its harsh cold climate did not experience major glaciations during the last ice age, including its severest stages around the Last Glacial Maximum (LGM). On the contrary, it is thought that glacial summer temperatures were comparable to the present day. However, evidence of glaciation has been found for several older glacial periods. We combine LGM experiments from the second and third phases of the Paleoclimate Modelling Intercomparison Project (PMIP2 and PMIP3) with sensitivity experiments using the Community Earth System Model (CESM). Together, these climate model experiments reveal that the intermodel spread in LGM summer temperatures in Siberia is much larger than in any other region of the globe and suggest that temperatures in Siberia are highly susceptible to changes in the imposed glacial boundary conditions, the included feedbacks and processes, and to the model physics of the different components of the climate model. We find that changes in the circumpolar atmospheric stationary wave pattern and associated northward heat transport drive strong local snow and vegetation feedbacks and that this combination explains the susceptibility of LGM summer temperatures in Siberia. This suggests that a small difference between two glacial periods in terms of climate, ice buildup or their respective evolution towards maximum glacial conditions can lead to strongly divergent summer temperatures in Siberia, allowing for the buildup of an ice sheet during some glacial periods, while during others, above-freezing summer temperatures preclude a multi-year snowpack from forming.

Description: Ice-penetrating radar1,2,3 and ice core drilling4 have shown that large parts of the north-central Greenland ice sheet are melting from below. It has been argued that basal ice melt is due to the anomalously high geothermal flux1,4 that has also influenced the development of the longest ice stream in Greenland1. Here we estimate the geothermal flux beneath the Greenland ice sheet and identify a 1,200-km-long and 400-km-wide geothermal anomaly beneath the thick ice cover. We suggest that this anomaly explains the observed melting of the ice sheet’s base, which drives the vigorous subglacial hydrology3 and controls the position of the head of the enigmatic 750-km-long northeastern Greenland ice stream5. Our combined analysis of independent seismic, gravity and tectonic data6,7,8,9 implies that the geothermal anomaly, which crosses Greenland from west to east, was formed by Greenland’s passage over the Iceland mantle plume between roughly 80 and 35 million years ago. We conclude that the complexity of the present-day subglacial hydrology and dynamic features of the north-central Greenland ice sheet originated in tectonic events that pre-date the onset of glaciation in Greenland by many tens of millions of years.

Description: After more than a century of geological research, the Cordilleran ice sheet of North America remains among the least understood in terms of its former extent, volume, and dynamics. Because of the mountainous topography on which the ice sheet formed, geological studies have often had only local or regional relevance and shown such a complexity that ice-sheet-wide spatial reconstructions of advance and retreat patterns are lacking. Here we use a numerical ice sheet model calibrated against field-based evidence to attempt a quantitative reconstruction of the Cordilleran ice sheet history through the last glacial cycle. A series of simulations is driven by time-dependent temperature offsets from six proxy records located around the globe. Although this approach reveals large variations in model response to evolving climate forcing, all simulations produce two major glaciations during marine oxygen isotope stages 4 (62.2–56.9 ka) and 2 (23.2–16.9 ka). The timing of glaciation is better reproduced using temperature reconstructions from Greenland and Antarctic ice cores than from regional oceanic sediment cores. During most of the last glacial cycle, the modelled ice cover is discontinuous and restricted to high mountain areas. However, widespread precipitation over the Skeena Mountains favours the persistence of a central ice dome throughout the glacial cycle. It acts as a nucleation centre before the Last Glacial Maximum and hosts the last remains of Cordilleran ice until the middle Holocene (6.7 ka).

Description: Paleoglaciological reconstructions based on glacial geological and geomorphological traces are the basis for testing and constraining numerical models of ice sheet extent and dynamics. In MAGIC-DML (“Mapping, Measuring and Modelling Antarctic Geomorphology and Ice Change, in Dronning Maud Land”) we aim to reconstruct the timing and pattern of ice surface elevation changes since the mid-Pliocene across western Dronning Maud Land, East Antarctica. The choice of study area is motivated by an absence of empirical field data and considerable ice sheet model uncertainties. We present a remote sensing-based mapping of glacial geomorphology and structures on the ice sheet surface for a coastal-inland transect including Ahlmannryggen-Borgmassivet-Kirwanveggen using high-resolution WorldView imagery. The primary aim of the study is to map traces of a thicker ice sheet on nunatak slopes that were formerly partly or entirely covered during surface highstands. Panchromatic and multispectral images were analysed in a multi-step procedure using ArcGIS, including image processing and mosaicking, visual feature recognition, and mapping. The identification of crucial landforms (such as till veneers and erratic boulders) required the adoption of some assumptions to differentiate, for example, till from regolith. Where patterned ground was mapped, we infer a presence of till rather than regolith since subglacial erosion is more likely to produce fine material than subaerial weathering. Huge boulders on plateau surfaces are mapped as erratics because they could not have been delivered by slope processes to local highpoints. Sediment veneers with ridges were mapped as till because the ridges are inferred to originate from reworking by ice. However, because of the necessity to invoke assumptions and the absence of crosscutting relationships between landforms, we are unsure of the reliability of derived paleo-ice sheet reconstructions. At face value, the data indicates thicker ice in the past from the altitude above present ice surface of till cover and erratics on some nunataks. Mapping and landform interpretations will be verified during the MAGIC-DML 2017-18 field season. The mapping of structures on the ice sheet surface is used to (i) infer ice flow characteristics; This was possible by assessing the distribution of primarily blue ice areas, crevasse fields, and supraglacial moraines, and by analysing their connection to wind directions and ablation rates; and (ii) yield target field routes to potential cosmogenic nuclide (CN) sampling locations along vertical transects on nunataks. The timing derived from CN dating will permit the delineation of ice sheet surface elevations as targets for ice sheet modeling.