Volume and mass changes of the Greenland ice sheet inferred from ICESat and GRACE

This study examines the recent evolution of the Greenland ice sheet and its six major drainage basins. Based on laser altimetry data acquired by the Ice, Cloud and Land Elevation Satellite (ICESat), covering the period September–November 2003 to February–March 2008, ice surface height changes and their temporal variations were inferred. Our refined repeat track analysis is solely based on ICESat data and is independent of external elevation models, since it accounts for both ice height changes and the local topography. From the high resolution ice height change pattern we infer an overall mean surface height trend of −0.12 ± 0.006 m yr⁻¹. Furthermore, the largest changes could be identified at coastal margins of the ice sheet, exhibiting rates of more than −2 m yr⁻¹. The total ice volume change of the entire ice sheet amounts to −205.4 ± 10.6 km³ yr⁻¹. In addition, we assessed mass changes from 78 monthly Gravity Recovery and Climate Experiment (GRACE) solutions. The Release-04 gravity field solutions of GeoForschungsZentrum Potsdam cover the period between August 2002 and June 2009. We applied an adjusted regional integration approach in order to minimize the leakage effects. Attention was paid to an optimized filtering which reduces error effects from different sources. The overall error assessment accounts for GRACE errors as well as for errors due to imperfect model reductions. In particular, errors caused by uncertainties in the glacial isostatic adjustment models could be identified as the largest source of errors. Finally, we determined both seasonal and long-term mass change rates. The latter amounts to an overall ice mass change of −191.2 ± 20.9 Gt yr⁻¹ corresponding to 0.53 ± 0.06 mm yr⁻¹ equivalent eustatic sea level rise. From the combination of the volume and mass change estimates we determined a mean density of the lost mass to be 930 ± 11 kg m⁻³. This value supports our applied density assumption 900 ± 30 kg m⁻³ which was used to perform the volume–mass-conversion of our ICESat results. Hence, mass change estimates from two independent observation techniques were inferred and are generally in good agreement.     

The effects of climatic variability on estimates of recharge from temperature profiles

Using heat as a tracer allows for estimation of ground water recharge rates based on subsurface temperature measurements. While possible in theory, it may be difficult in practice to discriminate the effects of climate from the effects of ground water advection. This study uses synthetic simulations to determine the influence of variability of ground surface temperature (GST) on the ability to estimate vertical specific discharge from temperature profiles. Results suggest that in cases where temperature measurements are sufficiently deep and specific discharge is sufficiently high, estimates of specific discharges will be reasonably accurate. Increasing the number of times temperatures are measured, or producing models that incorporate variations in GST, will increase the reliability of any studies using temperatures to estimate specific discharge. Furthermore, inversions of temperature measurements should be combined with other methods of estimating recharge rates to improve the reliability of recharge estimates.     

Regional geoid modeling in the area of subglacial Lake Vostok,Antarctica

We present a geoid model for the area of Lake Vostok, Antarctica, from a combination of local airborne gravity, ice-surface and ice-thickness data and a lake bathymetry model. The topography data are used for residual terrain modeling (RTM) in a remove–restore approach together with GOCE satellite data. The quasigeoid is predicted by least-squares collocation (LSC) and subsequently converted to geoid heights. Special aspects of that method in presence of an ice sheet are discussed.It is well known that a body freely floating in water is in a state of hydrostatic equilibrium (HE). This usually applies, e.g., to ice shelves or sea ice. However, it has been shown that this is valid also for the ice sheet covering the subglacial Lake Vostok. Thus, we demonstrate the use of such a refined regional geoid model for glaciological and geophysical applications by means of the HE surface of that lake. The mean quadratic residual geoid signal (0.56 m) w.r.t. the GOCE background model exceeds the residual variations of the estimated apparent lake level (ALL) (0.26 m) within the central part of the lake. An approach considering the actual geopotential at the ALL has been derived and subsequently applied. In this context, downward continuation of the potential field within the ice sheet as well as the latitudinal tilt of off-geoid equipotential surfaces are discussed. In view of the accuracy of the ice-thickness measurements that dominate the total error budget of the estimated ALL these effects are negligible. Thus, the HE surface of subglacial lakes may safely be described by a constant height bias in small-scale regional applications. However, field continuation is significant with respect to the formal uncertainty of the quasigeoid, which is at the level of 5 cm given that accurate airborne gravity data (±2 mGal) are available.     

Climate change inferred from analysis of borehole temperatures: First results from Alberta Basin,Canada

High quality temperature measurements have been made to depths of 30 to 220 m at 42 sites in 62 observational hydrogeological wells in Alberta. The temperature profiles commonly show near-surface inversions with a minimum temperature at depths of 30 to 50 m. Thermal modelling suggests a surface temperature history with warming reaching 2°C over the past 30 to 60 years. Recent climate warming evident from the analysis of the air temperature data in the region seems to provide at least a partial explanation of the increased ground temperatures. A sudden increase of the surface ground temperature caused by land clearing may be the other explanation, although modelling of such a sudden increase can only explain the observed temperature-depth data if the onset of such warming is 20–30 years old, which is in disagreement with the history of land development in the studied area. The effect of near-surface inversions of the temperature profiles also has been observed in the forested areas. The above support the climate based effect. The superposition of the climatic effect and man-made activity effect upon the ground warming is a very complicated process calling for considerably more research.     

Decoupled warming and monsoon precipitation in East Asia over the last deglaciation

Our understanding of the continental climate development in East Asia is mainly based on loess–paleosol sequences and summer monsoon precipitation reconstructions based on oxygen isotopes (δ¹⁸O) of stalagmites from several Chinese caves. Based on these records, it is thought that East Asian Summer Monsoon (EASM) precipitation generally follows Northern Hemisphere (NH) summer insolation. However, not much is known about the magnitude and timing of deglacial warming on the East Asian continent. In this study we reconstruct continental air temperatures for central China covering the last 34,000 yr, based on the distribution of fossil branched tetraether membrane lipids of soil bacteria in a loess–paleosol sequence from the Mangshan loess plateau. The results indicate that air temperature varied in phase with NH summer insolation, and that the onset of deglacial warming at ~ 19 kyr BP is parallel in timing with other continental records from e.g. Antarctica, southern Africa and South-America. The air temperature increased from ~ 15 °C at the onset of the warming to a maximum of ~ 27 °C in the early Holocene (~ 12 kyr BP), in agreement with the temperature increase inferred from e.g. pollen and phytolith data, and permafrost limits in central China.Comparison of the tetraether membrane lipid-derived temperature record with loess–paleosol proxy records and stalagmite δ¹⁸O records shows that the strengthening of EASM precipitation lagged that of deglacial warming by ca. 3 kyr. Moreover, intense soil formation in the loess deposits, caused by substantial increases in summer monsoon precipitation, only started around 12 kyr BP (ca. 7 kyr lag). Our results thus show that the intensification of EASM precipitation unambiguously lagged deglacial warming and NH summer insolation, and may contribute to a better understanding of the mechanisms controlling ice age terminations.     

Atmospheric temperature profiling by joint Raman,Rayleigh and Fe Boltzmann lidar measurements

Simultaneous and complete temperature profiles from near ground to about 100 km are essential for studying the dynamical coupling between different atmospheric layers. They are acquired by combining three different lidar techniques at Wuhan, China (30.5°N, 114.4°E). The atmospheric temperatures from about 3 to 25 km are calculated from the nitrogen molecule density profiles obtained from the N₂ vibrational Raman backscatter, while the atmospheric temperatures between 30 and ∼75 km are calculated by the standard Rayleigh scattering method. The temperatures in the 80–100 km altitude region are derived from the Fe Boltzmann technique. The temperature profiles measured by our lidar systems exhibit good agreement when compared with the radiosonde and satellite data, as well as the model. A Lomb–Scargle spectral analysis of the normalized temperature perturbations in the altitude range from 4 to 60 km shows that the spectral slopes of the vertical wave number spectra tended to −3 for large vertical wave numbers. This is consistent with the model predictions of saturated gravity wave spectra.     

Stability of massive ground ice bodies in University Valley,McMurdo Dry Valleys of Antarctica: Using stable O–H isotope as tracers of sublimation in hyper-arid regions

To date, studies of the stability of subsurface ice in the McMurdo Dry Valleys of Antarctica have been mainly based on climate-based vapor diffusion models. In University Valley (1800 m), a small glacier is found at the base of the head of the valley, and adjacent to the glacier, a buried body of massive ice was uncovered beneath 20–40 cm of loose cryotic sediments and sandstone boulders. This study assesses the origin and stability of the buried body of massive ice by measuring the geochemistry and stable O–H isotope composition of the ice and applies a sublimation and molecular diffusion model that accounts for the observed trends. The results indicate that the buried massive ice body represents an extension of the adjacent glacier that was buried by a rock avalanche during a cold climate period. The contrasting δ¹⁸O profiles and regression slope values between the uppermost 6 cm of the buried massive ice (upward convex δ¹⁸O profile and S_D-18O = 5.1) and that below it (progressive increase in δ¹⁸O and S_D-18O = 6.4) suggest independent post-depositional processes affected the isotope composition of the ice. The upward convex δ¹⁸O profile in the uppermost 6 cm is consistent with the ice undergoing sublimation. Using a sublimation and molecular diffusion model, and assuming that diffusion occurred through solid ice, the sublimation rate needed to fit the measured δ¹⁸O profile is 0.2 ? 10^? 3 mm yr^? 1, a value that is more similar to net ice removal rates derived from ³He data from cobbles in Beacon Valley till (7.0 ? 10^? 3 mm yr^? 1) than sublimation rates computed based on current climate (0.1–0.2 mm yr^?1). We suggest that the climate-based sublimation rates are offset due to potential ice recharge mechanisms or to missing parameters, particularly the nature and thermo-physical properties of the overlying sediments (i.e., temperature, humidity, pore structure and ice content, grain size).     

Ground response analysis of Kanpur soil along Indo-Gangetic Plains

The present work deals with 1D and 2D ground response analysis and liquefaction analysis of alluvial soil deposits from Kanpur region along Indo-Gangetic plains. Standard penetration tests and seismic down hole tests have been conducted at four locations namely IITK, Nankari village, Mandhana and Bithoor at 1.5 m interval up to a depth of 30 m below the ground surface to find the variation of penetration blows and the shear wave velocity along the depth. From the selected sites undisturbed as well as representative soil samples have been collected for detailed soil classification. The soil profiles from four sites have been considered for 1D and 2D ground response analysis by applying the free field motions of three Himalayan earthquakes namely Chamba earthquake (M_w—5.1), Chamoli earthquake (M_w—6.4) and Uttarkashi earthquake (M_w—6.5). An average value of Peak Ground Acceleration (PGA) obtained from 1D and 2D analysis is considered for liquefaction analysis and post-liquefaction settlement. The excess pore water pressure ratio is greater than 0.8 at a depth of 24 m from ground surface for IITK, Nankari village, Bithoor sites. More than 50% of post liquefaction settlement is contributed by layers from 21–30 m for all sites. In general, the soil deposits in Kanpur region have silty sand and sand deposits and are prone to liquefaction hazards due to drastic decrease of cyclic resistance ratio (CRR) at four chosen sites in Kanpur.     

Quantifying and specifying the solar influence on terrestrial surface temperature

This investigation is a follow-up of a paper in which we showed that both major magnetic components of the solar dynamo, viz. the toroidal and the poloidal ones, are correlated with average terrestrial surface temperatures. Here, we quantify, improve and specify that result and search for their causes.We studied seven recent temperature files. They were smoothed in order to eliminate the Schwabe-type (11 years) variations. While the total temperature gradient over the period of investigation (1610–1970) is 0.087 °C/century; a gradient of 0.077 °C/century is correlated with the equatorial (toroidal) magnetic field component. Half of it is explained by the increase of the Total Solar Irradiance over the period of investigation, while the other half is due to feedback by evaporated water vapour. A yet unexplained gradient of ?0.040 °C/century is correlated with the polar (poloidal) magnetic field. The residual temperature increase over that period, not correlated with solar variability, is 0.051 °C/century. It is ascribed to climatologic forcings and internal modes of variation.We used these results to study present terrestrial surface warming. By subtracting the above-mentioned components from the observed temperatures we found a residual excess of 0.31° in 1999, this being the triangularly weighted residual over the period 1990–2008.We show that solar forcing of the ground temperature associated with significant feedback is a regularly occurring feature, by describing some well observed events during the Holocene.     

Arctic rockwall retreat rates estimated using laboratory‐calibrated ERT measurements of talus cones in Longyeardalen,Svalbard

Holocene rockwall retreat rates quantify integral values of rock slope erosion and talus cone evolution. Here we investigate Holocene rockwall retreat of exposed arctic sandstone cliffs in Longyeardalen, central Svalbard and apply laboratory‐calibrated electrical resistivity tomography (ERT) to determine talus sediment thickness. Temperature–resistivity functions of two sandstone samples are measured in the laboratory and compared with borehole temperatures from the talus slope. The resistivity of the higher and lower‐porosity sandstone at relevant borehole permafrost temperatures defines a threshold range that accounts for the lithological variability of the dominant bedrock and debris material. This helps to estimate the depth of the transition from higher resistivities of ice‐rich debris to lower resistivities of frozen bedrock in the six ERT transects. The depth of the debris–bedrock transition in ERT profiles is confirmed by a pronounced apparent resistivity gradient in the raw data plotted versus depth of investigation. High‐resolution LiDAR‐scanning and ERT subsurface information were collated in a GIS to interpolate the bedrock surface and to calculate the sediment volume of the talus cones. The resulting volumes were referenced to source areas to calculate rockwall retreat rates. The rock mass strength was estimated for the source areas. The integral rockwall retreat rates range from 0.33 to 1.96 mm yr^–1, and are among the highest rockwall retreat rates measured in arctic environments, presumably modulated by harsh environmental forcing on a porous sandstone rock cliff with a comparatively low rock mass strength. Here, we show the potential of laboratory‐calibrated ERT to provide accurate estimates of rockwall retreat rates even in ice‐rich permafrost talus slopes. Copyright © 2012 John Wiley & Sons, Ltd.     

Cosmogenic 36Cl glacial chronologies of the Late Quaternary glaciers on Mount Geyikdağ in the Eastern Mediterranean

We report the timing of glaciations during the Late Quaternary in the central Taurus Mountains of Turkey in the Eastern Mediterranean. Forty moraine samples from three glacial valleys on Mount Geyikdağ (36.53°N, 32.10^°E, 2877 m), near the Eastern Mediterranean coast of Turkey, were dated with in-situ cosmogenic ³⁶Cl. These glacial valleys are located on the southern flank of the mountain and were filled with few km long glaciers that terminated at elevations of about 1750 m above sea level. Three glacial retreats/advances were determined in this study. During the Last Glacial Maximum (LGM), glaciers reached their maximum positions at 20.6 ± 0.6 ka ago (±1σ). This date is in accordance with the timing of local maximum ice extent, represented by piedmont glaciers in the northern side of the mountain. Glaciers started to retreat after the LGM and shortly stabilized or re-advanced two times before they completely vanished out. The first stage ended before 13.7 ± 0.8 ka ago during the Late-glacial. The last glaciation occurred during the Holocene and ended between 9.6 ± 1.4 ka and 5.9 ± 0.5 ka ago. Later, glaciers mostly vanished from the study area, but a few rock glaciers developed during the Late Holocene. Glacial chronologies of Mount Geyikdağ are mostly comparable with the globally observed advances elsewhere.     

The penultimate deglaciation: High-resolution paleoceanographic evidence from a north–south transect along the eastern Nordic Seas

The penultimate termination has been studied with focus on oceanographic changes in the eastern Nordic Seas and the influence of these changes on the surrounding ice sheets and vice versa. Repeatedly, major changes in the strength of the Atlantic Meridional Overturning Circulation (AMOC) occurred during the studied interval. Times of strong overturning and increased heat transport northwards were of importance in triggering one minor and two major disintegration events. The two major disintegration events were separated by a deglacial pause, characterized by a strong AMOC. The same variability is seen throughout the eastern Nordic Seas, from the Faeroe–Shetland Channel in the south to the Fram Strait in the north. Some of the oceanographic changes occurring during the penultimate termination are comparable with changes seen through the last deglaciation. Reduced winter moisture flux and increased summer melting due to the present insolation forcing further amplified the rate of ice sheet disintegration. Calculated sea-level change through TII shows a mean change of 121 ± 4 m, 41 ± 16 m in the first step and 80 ± 13 m in the last step.     

A 3D resistivity model derived from the transient electromagnetic data observed on the Araba fault,Jordan

72 inloop transient electromagnetic soundings were carried out on two 2 km long profiles perpendicular and two 1 km and two 500 m long profiles parallel to the strike direction of the Araba fault in Jordan which is the southern part of the Dead Sea transform fault indicating the boundary between the African and Arabian continental plates. The distance between the stations was on average 50 m.The late time apparent resistivities derived from the induced voltages show clear differences between the stations located at the eastern and at the western part of the Araba fault. The fault appears as a boundary between the resistive western (ca. 100 Ωm) and the conductive eastern part (ca. 10 Ωm) of the survey area. On profiles parallel to the strike late time apparent resistivities were almost constant as well in the time dependence as in lateral extension at different stations, indicating a 2D resistivity structure of the investigated area.After having been processed, the data were interpreted by conventional 1D Occam and Marquardt inversion. The study using 2D synthetic model data showed, however, that 1D inversions of stations close to the fault resulted in fictitious layers in the subsurface thus producing large interpretation errors. Therefore, the data were interpreted by a 2D forward resistivity modeling which was then extended to a 3D resistivity model. This 3D model explains satisfactorily the time dependences of the observed transients at nearly all stations.     

Application of the Brewster angle to quantify the dielectric properties of ground ice formations

Previous studies of ground ice using moveout type ground-penetrating radar (GPR) surveys indicate that the dielectric permittivity can constrain the type of ground ice present in the subsurface. Due to the high-loss nature of the active layer over permafrost targets, however, the signal strength of GPR signals is often insufficient to resolve the basal boundary required for determining the dielectric permittivity of an underlying unit. We apply a non-conventional antenna orientation and post-processing method to determine the dielectric permittivity of the unit underlying the lowest resolvable boundary. We conduct moveout surveys using a 450 MHz GPR with collinear parallel oriented antennas on two adjacent ground ice formations in the region of Thomas Lee Inlet, Devon Island, Nunavut. We exploit the Brewster angle to calculate the approximate dielectric permittivity of ground ice formations below the active layer. The results agree within 1 dielectric unit with on-ice permittivity measurements made during a complementary study of the site.     

Calculating basal temperatures in ice sheets: an Excel spreadsheet method

The flow of ice sheets and their geomorphological impact is greatly influenced by their basal thermal regime. Calculations of basal temperatures in ice sheets are therefore fundamental in evaluating glacier dynamics and in determining the spatial distribution of zones of erosion and deposition beneath ice masses. Calculations of basal temperatures are not frequently attempted, however, primarily because of the techniques required to solve the heat conduction equation between the ice surface and the base. This paper describes a new Excel spreadsheet method of solving this equation that can readily be applied to both former and contemporary ice sheets. The application of the spreadsheet is illustrated with two examples. The first provides a calculation of basal thermal regime beneath the north eastern part of the Scottish ice sheet during the last glacial maximum; the second shows how basal ice temperatures can be calculated beneath the modern Antarctic ice sheet. Copyright © 2002 John Wiley & Sons, Ltd.     

Resolving sea level contributions by identifying fingerprints in time-variable gravity and altimetry

We have studied the ability of the GRACE gravimetry mission and Jason-1 altimetry to resolve ice and glacier induced contributions to sea level rise, by means of a fingerprint method. Here, the signals from ice sheet and land glacier changes, steric changes, glacial isostatic adjustment and terrestrial hydrology are assumed to have fixed spatial patterns. In a joint inversion using GRACE and Jason-1 data the unknown temporal components can then be estimated by least-squares. In total, we estimate temporal components for up to ∼ 80 individual patterns. From a propagation of the full error-covariance from GRACE and a diagonal error-covariance from Jason-1 altimetry we find that: (1) GRACE almost entirely explains the mass related parameters in the joint inversion, (2) an inversion using only Jason-1 data has a marginal ability to estimate the mass related parameters, while the steric parameters have much better formal accuracy. In terms of mean sea level rise the steric patterns have a maximum formal accuracy of 0.01 mm for an 11 week running mean. In general, strong negative error correlations (ρ < − 0.9) exists between the high and low elevation parts of the ice sheet drainage basins, when those are estimated independently. The largest formal errors found are in the order of 40 Gton for small high elevation subbasins in the southern Greenland ice sheet, which are difficult to separate. In a simplified joint inversion, merging high and low elevation basins, we have investigated the ability of the GRACE and Jason-1 data to separate the geocenter motion into a present-day contribution and a contribution from glacial isostatic adjustment (GIA). We find that the GIA related signal is larger than the present-day component with a maximum of −0.71 mm/year in the Z direction. Total geocenter motion rates are found to be −0.28, 0.43, −1.08 mm/year for the X, Y and Z components, respectively. The inversion results have been propagated to the Jason-1 along-track measurements. Over the time period considered, we see that a large part of the variability in the Pacific, Atlantic and Indian ocean can be explained by our inversion results. The applied inversion method therefore seems a feasible way to separate steric from mass induced sea level changes. At the same time, the joint inversion would benefit from more advanced parameterizations, which may aid in fitting remaining signal from altimetry.     

Postglacial isostatic adjustment in a self-gravitating spherical earth with power-law rheology

Since microphysics cannot say definitively whether the rheology of the mantle is linear or non-linear, the aim of this paper is to constrain mantle rheology from observations related to the glacial isostatic adjustment (GIA) process—namely relative sea-levels (RSLs), land uplift rate from GPS and gravity-rate-of-change from GRACE. We consider three earth model types that can have power-law rheology (n = 3 or 4) in the upper mantle, the lower mantle or throughout the mantle. For each model type, a range of A parameter in the creep law will be explored and the predicted GIA responses will be compared to the observations to see which value of A has the potential to explain all the data simultaneously. The coupled Laplace finite-element (CLFE) method is used to calculate the response of a 3D spherical self-gravitating viscoelastic Earth to forcing by the ICE-4G ice history model with ocean loads in self-gravitating oceans. Results show that ice thickness in Laurentide needs to increase significantly or delayed by 2 ka, otherwise the predicted uplift rate, gravity rate-of-change and the amplitude of the RSL for sites inside the ice margin of Laurentide are too low to be able to explain the observations. However, the ice thickness elsewhere outside Laurentide needs to be slightly modified in order to explain the global RSL data outside Laurentide. If the ice model is modified in this way, then the results of this paper indicate that models with power-law rheology in the lower mantle (with A ∼ 10⁻³⁵ Pa⁻³ s⁻¹ for n = 3) have the highest potential to simultaneously explain all the observed RSL, uplift rate and gravity rate-of-change data than the other model types.     

A quantitative examination of modern sedimentary lithofacies formation on the glacially influenced Gulf of Alaska continental shelf

Continental-shelf lithofacies are described from a series of cores collected in the northern Gulf of Alaska, a high-energy paraglacial shelf experiencing rapid rates of sediment accumulation. Short-lived tracers (²³⁴Th and chlorophyll-a) indicate that during the annual peak in fluvial sediment input (summer), biologic sediment mixing coefficients in the surficial seabed are generally lower than other coastal environments (<20 cm² yr⁻¹) and mixing extends downward <10 cm.²¹⁰Pb geochronology indicates that sediment accumulation rates (time scales of 10–100 yr) are 0.1–3 cm yr⁻¹. The measured bioturbation and accumulation rates lead to predictions of moderate to bioturbated lithofacies, as observed. Primary depositional fabric is preferentially preserved where sediment accumulation rates >2 cm yr⁻¹ and non-steady sediment deposition occurs. Depositional fabric is also observed in strata at 50–100 m water depths and is similar in appearance to beds that may form through deposition of wave-induced fluid-mud flows, which have been observed forming on other shelves with moderate to high wave energy. Five general lithofacies can be identified for the study area: inner-shelf sand facies, interbedded sandy mud facies, moderate-to-well-bioturbated mud facies, gravelly mud facies, and Tertiary bedrock facies. The moderate-to-well-bioturbated mud facies is areally dominant, representing over 50% of the shelf area, although roughly equal volumes (∼0.4 km³) of strata with some preservation of primary fabric are annually accumulating. Lithofacies on this paraglacial shelf generally resemble mid- and low-latitude allochthonous shelf strata to a much greater degree than Holocene glacimarine strata formed on shelves dominated by icebergs and floating ice shelves. Paraglacial strata may be differentiated from non-glacial shelf strata by lower organic carbon concentrations, a relatively lower degree of bioturbation, and increased preservation of primary depositional fabric.     

Chronology of Lateglacial ice flow reorganization and deglaciation in the Gotthard Pass area,Central Swiss Alps,based on cosmogenic 10Be and in situ 14C

We reconstruct the timing of ice flow reconfiguration and deglaciation of the Central Alpine Gotthard Pass, Switzerland, using cosmogenic ¹⁰Be and in situ ¹⁴C surface exposure dating. Combined with mapping of glacial erosional markers, exposure ages of bedrock surfaces reveal progressive glacier downwasting from the maximum LGM ice volume and a gradual reorganization of the paleoflow pattern with a southward migration of the ice divide. Exposure ages of ∼16–14 ka (snow corrected) give evidence for continuous early Lateglacial ice cover and indicate that the first deglaciation was contemporaneous with the decay of the large Gschnitz glacier system. In agreement with published ages from other Alpine passes, these data support the concept of large transection glaciers that persisted in the high Alps after the breakdown of the LGM ice masses in the foreland and possibly decayed as late as the onset of the Bølling warming. A younger group of ages around ∼12–13 ka records the timing of deglaciation following local glacier readvance during the Egesen stadial. Glacial erosional features and the distribution of exposure ages consistently imply that Egesen glaciers were of comparatively small volume and were following a topographically controlled paleoflow pattern. Dating of a boulder close to the pass elevation gives a minimum age of 11.1 ± 0.4 ka for final deglaciation by the end of the Younger Dryas. In situ ¹⁴C data are overall in good agreement with the ¹⁰Be ages and confirm continuous exposure throughout the Holocene. However, in situ ¹⁴C demonstrates that partial surface shielding, e.g. by snow, has to be incorporated in the exposure age calculations and the model of deglaciation.     

Modeled steric and mass-driven sea level change caused by Greenland Ice Sheet melting

Meltwater from the Greenland Ice Sheet (GIS) has been a major contributor to sea level change in the recent past. Global and regional sea level variations caused by melting of the GIS are investigated with the finite element sea-ice ocean model (FESOM). We consider changes of local density (steric effects), mass inflow into the ocean, redistribution of mass, and gravitational effects. Five melting scenarios are simulated, where mass losses of 100, 200, 500, and 1000 Gt/yr are converted to a continuous volume flux that is homogeneously distributed along the coast of Greenland south of 75°N. In addition, a scenario of regional melt rates is calculated from daily ice melt characteristics. The global mean sea level modeled with FESOM increases by about 0.3 mm/yr if 100 Gt/yr of ice melts, which includes eustatic and steric sea level change. In the global mean the steric contribution is one order of magnitude smaller than the eustatic contribution. Regionally, especially in the North Atlantic, the steric contribution leads to strong deviations from the global mean sea level change. The modeled pattern mainly reflects the structure of temperature and salinity change in the upper ocean. Additionally, small steric variations occur due to local variability in the heat exchange between the atmosphere and the ocean. The mass loss has also affects on the gravitational attraction by the ice sheet, causing spatially varying sea level change mainly near the GIS, but also at greater distances. This effect is accounted for by using Green's functions.