Combining geomorphological and geodetic data to determine postglacial tilting in Manitoba

Estimates of postglacial rebound in central North America from Laurentide deglaciation to the present time are uncertain as a result of lack of data from the continental interior. A more precise knowledge of postglacial tilt history will assist studies of the evolution of the major lakes in Manitoba and will facilitate the engineering and environmental management of the present-day hydrological system. This paper explores the benefits of combining geomorphological data with high-precision, real-time geodetic data (GPS positioning and absolute gravity) and lake-gauge tilt data now being collected for postglacial rebound studies in Manitoba and adjacent regions in the US. Presently-available data sets representing these data types are (1) tilting of the 9.5 kyr B.P. Campbell strand line south and west of Lake Winnipeg, (2) the rate of decrease in absolute gravity values measured from 1987 to 1995 at Churchill, Manitoba, and (3) the present-day regional tilt rate derived from water-level gauges in southern Manitoba lakes. These data are compared to theoretical predictions based on the published ICE-3G loading history and on a model of Earth rheology characterized by a 1066B density and elastic structure, an upper-mantle viscosity of 10 ²¹Pa s, a lower-mantle viscosity of 2 × 10 ²¹Pa s, and a lithosphere thickness of 120 km (Tushingham & Peltier, 1991). All three data types show disagreement in Manitoba with ICE-3G and the standard Earth model. ICE-4G does better but could not be investigated in any detail. The constraints on model parameters provided by the different data types were investigated by varying, one at a time, three key parameters, (1) the thickness of the lithosphere in excess of 120 km, (2) the lower mantle viscosity, and (3) the thickness of Laurentide ice over the Prairies, to obtain better fits to the data. The present data do not appear to constrain lithosphere thickness in excess of 120 km very well. While both the Campbell strand line data and the Churchill absolute gravity data are consistent with an increase in lower-mantle viscosity, the present-day, lake-gauge data are not. All three data types are consistent with a thinning of the Laurentide ice-sheet over the Prairies relative to the ICE-3G model. Simultaneous adjustment of model parameters with the advantage of anticipated new data in Manitoba and adjacent regions in the US will lead to better understanding of the trade-offs between Earth rheology and ice sheet history and hence to an improved Laurentide postglacial rebound model.     

Spectral-finite-element approach to two-dimensional electromagnetic induction in a spherical earth

We present a spectral-finite-element approach to the 2-D forward problem for electromagnetic induction in a spherical earth. It represents an alternative to a variety of numerical methods for 2-D global electromagnetic modelling introduced recently (e.g. the perturbation expansion approach, the finite difference scheme). It may be used to estimate the effect of a possible axisymmetric structure of electrical conductivity of the mantle on surface observations, or it may serve as a tool for testing methods and codes for 3-D global electromagnetic modelling. The ultimate goal of these electromagnetic studies is to learn about the Earth's 3-D electrical structure.
Since the spectral-finite-element approach comes from the variational formulation, we formulate the 2-D electromagnetic induction problem in a variational sense. The boundary data used in this formulation consist of the horizontal components of the total magnetic intensity measured on the Earth's surface. In this the variational approach differs from other methods, which usually use spherical harmonic coefficients of external magnetic sources as input data. We verify the assumptions of the Lax-Milgram theorem and show that the variational solution exists and is unique. The spectral-finite-element approach then means that the problem is parametrized by spherical harmonics in the angular direction, whereas finite elements span the radial direction. The solution is searched for by the Galerkin method, which leads to the solving of a system of linear algebraic equations. The method and code have been tested for Everett & Schultz's (1995) model of two eccentrically nested spheres, and good agreement has been obtained.     

Paleoshoreline evidence for postglacial tilting in Southern Manitoba

Detailed air photo interpretation and four seasons of field mapping and surveying in southern Manitoba have revealed that the once-level paleoshorelines of Lake Winnipegosis and Dauphin Lake and the Burnside shoreline of former Lake Agassiz have been tilted up to the northeast by postglacial differential rebound. Our investigation has also revealed that Lake Winnipegosis has the best preserved paleoshoreline record of any of the large lakes in southern Manitoba, including lakes Winnipeg and Manitoba. This is because northeasterly uptilting shifts the region's lakes to the southwest. Lakes with southern outlets, like Lake Winnipegosis, undergo general regression as the outlet is lowered relative to the rest of the basin. Lakes with northern outlets, like lakes Winnipeg and Manitoba, undergo general transgression as northeasterly uptilting raises the outlet relative to the rest of the basin. Along the northeastern shore of Lake Winnipegosis a staircase of at least 32 abandoned Winnipegosis shorelines exists that is consistent with northeasterly tilting. The Dawson level represents the major mid-Holocene highstand on Lake Winnipegosis. It persisted for about 500 years, peaking at 5290 ¹⁴C yr B.P. (early Dawson) and then falling about 3 m by 4740 ¹⁴C yr B.P. (late Dawson). The early Dawson shoreline is tilted at 13.5 cm km^-1 in a direction N24.3°E. Three other shorelines informally named shoreline 4, shoreline 3, and shoreline 2 are also tilted up to the northeast. Their radiocarbon ages (and slopes in cm km^-1) are 3330 yr B.P. (2.2), 1510 yr B.P. (1.3), and 1080 yr B.P. (0.7), respectively. On Dauphin Lake shoreline IV is the oldest level mapped for this study. It has a ¹⁴C age of 7910 yr B.P. and is tilted at 21.7 cm km^-1 in a direction N44.4°E. The Id shoreline marks the major mid-Holocene highstand for Dauphin Lake. It peaked at 4640 ¹⁴C yr B.P. followed by a rapid decline of about 1 m to the Ib shoreline, which is dated at 4320 ¹⁴C yr B.P. Id is tilted up at 8.8 cm km^-1 in a direction N53.4°E. The next major shoreline is Ia3 which has a ¹⁴C age of 3020 yr B.P. and is tilted up at 5.3 cm km^-1 in a direction N62.3°E. Tilt directions are significantly more easterly for the Dauphin Lake shorelines than those from Lake Winnipegosis or any of the much older Lake Agassiz shorelines. Taken together, the Winnipegosis and Dauphin isobases indicate that the direction of tilt in southern Manitoba is more complex than a simple uni-directional pattern. The observed pattern of tilting for paleoshorelines in southern Manitoba agrees better with predictions derived from the recently revised loading history model ICE-4G than with those from its predecessor ICE-3G. In general, the calculated tilt based on the ICE-3G load tends to exceed the observed tilt, while ICE-4G tends to underestimate it. Both ice load models appear to disagree most with our observed tilts in this region during the interval before about 9000 cal yr B.P., when deglaciation was proceeding rapidly and the large water load associated with Lake Agassiz covered the region. Because both of these ice load models have been estimated mainly from a global data set of relative sea level curves from marine coast sites, it is not unexpected that model tilts derived from them do not agree well with observations in the North American continental interior. The pattern of postglacial crustal deformation for southern Manitoba described in this paper could be used to further refine ice load models for the North American continental interior.     

Glacial isostatic adjustment in Fennoscandia for a laterally heterogeneous earth

Glaciation and deglaciation in Fennoscandia during the last glacial cycles has significantly perturbed the Earth's equilibrium figure. Changes in the Earth's solid and geoidal surfaces due to external and internal mass redistributions are recorded in sequences of ancient coastlines, now either submerged or uplifted, and are still visible in observations of present‐day motions of the surface and glacially induced anomalies in the Earth's gravitational field. These observations become increasingly sophisticated with the availability of GPS measurements and new satellite gravity missions.
Observational evidence of the mass changes is widely used to constrain the radial viscosity structure of the Earth's mantle. However, lateral changes in earth model properties are usually not taken into account, as most global models of glacial isostatic adjustment assume radial symmetry for the earth model. This simplifying assumption contrasts with seismological evidence of significant lateral variations in the Earth's crust and upper mantle throughout the Fennoscandian region.
We compare predictions of glacial isostatic adjustment based on an ice model over the Fennoscandian region for the last glacial cycle for both radially symmetric and fully 3‐D earth models. Our results clearly reveal the importance of lateral variations in lithospheric thickness and asthenospheric viscosity for glacially induced model predictions. Relative sea‐level predictions can differ by up to 10–20 m, uplift rate predictions by 1–3 mm yr⁻¹ and free‐air gravity anomaly predictions by 2–4 mGal when a realistic 3‐D earth structure as proposed by seismic modelling is taken into account.     

Iceocean mass balance during the Late Pleistocene glacial cycles in view of CHAMP and GRACE satellite missions

During the last glacial cycles, global sea level dropped several times by about 120 m and large ice sheets covered North America, northern Europe and Antarctica during the glacial stages. The changes in the iceocean mass balance have displaced mantle material mainly via viscous flow, and the perturbation of the equilibrium figure of the Earth by glacial isostatic adjustment is still observable today in timedependent changes of gravitational and rotational observations. Contemporary iceocean mass balance from volume changes of polar ice caps also contributes to secular variations of the Earth's gravitational field.
In the near future, several satellite gravity missions will significantly improve the accuracy of the observed timedependent gravitational field. In view of the expected improvements in the observations, we predict glacially induced perturbations of the gravitational field, induced by Late Pleistocene and contemporary ice volume changes, for a variety of radial mantle viscosity profiles. We assess the degree of uncertainty for the glacially induced contributions to gravitational and rotational parameters, both in the spectral and the spatial domain.
Predictions of power spectra for the glacially induced freeair gravity and geoid anomalies are about one order of magnitude lower than the observed values, and uncertainties arising from different plausible viscosity profiles are around 0.150.4 mGal and 0.21.5 m, respectively. Uncertainties from different ice models are of secondary importance for the predicted power spectra. Predicted secular changes in geoid anomalies in formerly glaciated areas are mainly controlled by the viscosity profile and contemporary ice volume changes. We also show that the simple threelayer viscosity profiles currently employed for the majority of postglacial rebound studies represent a limited subset for model predictions of the timedependent gravitational field.     

Spectral analysis using orthonormal functions with a case study on the sea surface topography

b
Spherical harmonics are orthonormalized using the Gram-Schmidt process in a function space. The problem of linear dependence of spherical harmonics over the oceans is studied using the Gram matrices and consequently three sets of orthonormal (ON) functions have been constructed. For the process an efficient formula for computing inner products of spherical harmonics has been developed. Important spectral properties of the ON functions are addressed. The ON functions may be used for representing the sea surface topography (SST) in the analysis of satellite altimeter data. The geoid error can be transformed to a representation by the ON functions and hence the comparison of powers of the geoid error and the SST signal only over the oceans is possible, leading to a better way of determining the cut-off frequency of the SST in the simultaneous solution using satellite altimeter data. As a case study, the modified Levitus SST is expanded into the ON functions. The results show that 99.90 per cent of that signal's energy is contained within degree 24 of the orthonormal functions. Such expansions also render better spectral behaviour of oceanic signals as compared to that from spherical harmonic expansions. The study shows that these generalized Fourier functions are suitable for spectral analyses of oceanic signals and they can be applied to future altimetric mission where the geoid and the SST are to be recovered.     

Two-dimensional modelling of subduction zone anisotropy with application to southwestern Japan

We present a series of 2-D numerical models of viscous flow in the mantle wedge induced by a subducting lithospheric plate. We use a kinematically defined slab geometry approximating the subduction of the Philippine Sea plate beneath Eurasia. Through finite element modelling we explore the effects of different rheological and thermal constraints (e.g. a low-viscosity region in the wedge corner, power law versus Newtonian rheology, the inclusion of thermal buoyancy forces and a temperature-dependent viscosity law) on the velocity and finite strain field in the mantle wedge. From the numerical flow models we construct models of anisotropy in the wedge by calculating the evolution of the finite strain ellipse and combining its geometry with appropriate elastic constants for effective transversely isotropic mantle material. We then predict shear wave splitting for stations located above the model domain using expressions derived from anisotropic perturbation theory, and compare the predictions to ∼500 previously published shear wave splitting measurements from seventeen stations of the broad-band F-net array located in southwestern Japan. Although the use of different model parameters can have a substantial effect on the character of the finite strain field, the effect on the average predicted splitting parameters is small. However, the variations with backazimuth and ray parameter of individual splitting intensity measurements at a given station for different models are often different, and rigorous analysis of details in the splitting patterns allows us to discriminate among different rheological models for flow in the mantle wedge. The splitting observed in southwestern Japan agrees well with the predictions of trench-perpendicular flow in the mantle wedge along with B-type olivine fabric dominating in a region from the wedge corner to about 125 km from the trench.     

Time dependent piezomagnetic fields in viscoelastic medium

We investigated time dependent piezomagnetic fields due to volcanic sources embedded in a viscoelastic, homogeneous half-space. Especially in volcanic areas, the presence of inhomogeneous materials and high temperatures produce a lower effective viscosity of the Earth's crust that calls for considering anelastic properties of the medium. Piezomagnetic properties are carried by grains of titano-magnetite, which occupy only a small fraction of ordinary rock volume and are supposed to be elastic, while the non-magnetic surrounding matrix is assumed to be viscoelastic. From all the possible rheological models, we investigated two cases in which the bulk modulus is purely elastic and the shear modulus relaxes as: (i) a Maxwell solid and (ii) a standard linear solid (SLS). We applied the Correspondence Principle to the analytical elastic solutions for pressurized spherical sources and dislocation sources in order to determine the time dependent piezomagnetic fields in a viscoelastic medium. The piezomagnetic field completely vanishes after the relaxation process for a Maxwell rheology, whereas it is found to decrease over time and reach some finite offset value for a SLS rheology. These different behaviours provide helpful hints in understanding the temporal evolution of piezomagnetic anomalies in volcanic regions.     

Dynamics of the last glaciation in eastern Svalbard as inferred from glacier-movement indicators

Glacial striae and other ice movement indicators such as roche moutonées, glacial erratics, till fabric and glaciotectonic deformation have been used to reconstruct the Late Weichselian ice movements in the region of eastern Svalbard and the northern Barents Sea. The ice movement pattern may be divided into three main phases: (1) a maximum phase when ice flowed out of a centre east or southeast of Kong Karls Land. At this time the southern part of Spitsbergen was overrun by glacial ice from the Barents Sea; (2) the phase of deglaciation of the Barents Sea Ice Sheet, when an ice cap was centred between Kong Karls Land and Nordaustlandet. At the same time ice flowed southwards along Storfjorden; and (3) the last phase of the Late Weichselian glaciation in eastern Svalbard is represented by local ice caps on Spitsbergen, Nordaustlandet, Barentsoya and Edgeøya.
The reconstructed ice flow pattern during maximum glaciation is compatible with a centre of uplift in the northern Barents Sea as shown by isobase reconstructions and suggested by isostatic modelling.     

The onset of Pleistocene glaciation in the Barents Sea: implications for glacial isostatic adjustment

In this paper the effect of a delayed onset of glaciation in the Barents Sea on glacial isostatic adjustment is investigated. The model calculations solve the sea-level equation governing the total mass redistributions associated with the last glaciation cycle on a spherically symmetric, linear, Maxwell viscoelastic earth for two different scenarios for the growth phase of the Barents Sea ice sheet. In the first ice model a linear growing history is used for the Barents Sea ice sheet, which closely relates its development to the build-up of other major Late Pleistocene ice sheets. In the second ice model the accumulation of the Barents Sea ice sheet is restricted to the last 6 ka prior to the last glacial maximum.
The calculations predict relative sea levels, present-day radial velocities, and gravity anomalies for the area formerly covered by the Weichselian ice sheet. The results show that observed relative sea levels in the Barents Sea are appropriate for distinguishing between the different glaciation histories. In particular, present-day observables such as the free-air gravity anomaly over the Barents Sea, and the present-day radial velocities are sensitive to changes in the glaciation history on this scale.
A palaeobathymetry derived from relative sea-level predictions before the last glacial maximum based on the second ice model essentially agrees with a palaeobathymetry derived by Lambeck (1995). The additional emerged areas provide centres for the build-up of an ice sheet and thus support the theory of Hald, Danielsen & Lorentzen (1990) and Mangerud et al. (1992) that the Barents Sea was an essentially marine environment shortly before the last glacial maximum.     

RESPONSE OF A TEMPERATE ALPINE VALLEY GLACIER TO CLIMATE CHANGE AT THE DECADAL SCALE

Glacier advance and recession are considered key indicators of climate change. Understanding the relationship between climatic variations and glacial responses is crucial. Here, we apply archival digital photogrammetry to reconstruct the decadal scale glacial history of an unmonitored Alpine valley glacier, the Haut Glacier d'Arolla, Switzerland, and we use the data generated to explore the linkages between glacier recession and climate forcing. High precision digital elevation models were derived. They show continual recession of the glacier since 1967, associated with long‐term climatic amelioration but only a weak reaction to shorter‐term climatic deterioration. Glacier surface velocity estimates obtained using surface particle tracking showed that, unlike for most Swiss glaciers during the late 1970s and early 1980s, ice mass flux from the accumulation zone was too low to compensate for the effects of glacier thinning and subsequent snout recession, especially during the rapid warming that occurred through the 1980s. The results emphasise the dangers of inferring glacier response to climate forcing from measurements of the terminus position only and the importance of using remote sensing methods as an alternative, especially where historical imagery is available.     

Parametrizing surface wave tomographic models with harmonic spherical splines

We present a mathematical framework and a new methodology for the parametrization of surface wave phase-speed models, based on traveltime data. Our method is neither purely local, like block-based approaches, nor is it purely global, like those based on spherical harmonic basis functions. Rather, it combines the well-known theory and practical utility of the spherical harmonics with the spatial localization properties of spline basis functions. We derive the theoretical foundations for the application of harmonic spherical splines to surface wave tomography and summarize the results of numerous numerical tests illustrating the performance of a practical inversion scheme based upon them. Our presentation is based on the notion of reproducing-kernel Hilbert spaces, which lends itself to the parametrization of fully 3-D tomographic earth models that include body waves as well.     

Semi-analytical solution for viscous Stokes flow in two eccentrically nested spheres

Applying the infinite Prandtl number approximation, a semi-analytical solution for computing 2-D axisymmetric viscous Stokes flow in a model consisting of two eccentrically nested spheres of different viscosities is derived. Since numerical codes based on spectral or finite techniques for modelling mantle flow in a spherical geometry in the presence of lateral viscosity variation are becoming more and more popular, reliable examples for testing and validating such codes are extremely useful. The eccentrically nested sphere solution was used to test a numerical algorithm based on a mixed spherical-harmonic finite-element formulation of the Stokes problem, and good agreement was obtained.     

Quantifying the erosional impact of the Fennoscandian ice sheet in the Torneträsk–Narvik corridor, northern Sweden, based on cosmogenic radionuclide data

Despite spectacular landform evidence of a dominant role for glacial action in shaping landscapes under former northern hemisphere ice sheets, there is little quantitative evidence for rates and patterns of erosion associated with specific glaciations. Here we use cosmogenic nuclide data to assess rates of subglacial erosion underneath the Fennoscandian ice sheet. By testing whether there are remnant nuclide concentrations in samples taken from sites that include both relict areas and features and landscapes typically associated with vigorous glacial erosion, we can constrain the level and pattern of modification that resulted from the last glaciation. Cosmogenic ¹⁰Be and ³⁶Cl data from the Torneträsk region confirm the temporal and spatial variability of glacial erosion suggested by geomorphological mapping. At some sites, glacial erosion estimates in what appear to be heavily scoured areas indicate erosion of only c. 2 ± 0.4 m of bedrock, based on cosmogenic nuclide inheritance. This implies that the generation of severely scoured terrain in this study area required multiple glaciations. The overall modification produced by ice sheets along glacial corridors may be more restricted than previously thought, or may have occurred preferentially during earlier Quaternary glacial periods.     

Quantifying the erosional impact of the Fennoscandian ice sheet in the Torneträsk–Narvik corridor, northern Sweden, based on cosmogenic radionuclide data

Despite spectacular landform evidence of a dominant role for glacial action in shaping landscapes under former northern hemisphere ice sheets, there is little quantitative evidence for rates and patterns of erosion associated with specific glaciations. Here we use cosmogenic nuclide data to assess rates of subglacial erosion underneath the Fennoscandian ice sheet. By testing whether there are remnant nuclide concentrations in samples taken from sites that include both relict areas and features and landscapes typically associated with vigorous glacial erosion, we can constrain the level and pattern of modification that resulted from the last glaciation. Cosmogenic ¹⁰Be and ³⁶Cl data from the Torneträsk region confirm the temporal and spatial variability of glacial erosion suggested by geomorphological mapping. At some sites, glacial erosion estimates in what appear to be heavily scoured areas indicate erosion of only c. 2 ± 0.4 m of bedrock, based on cosmogenic nuclide inheritance. This implies that the generation of severely scoured terrain in this study area required multiple glaciations. The overall modification produced by ice sheets along glacial corridors may be more restricted than previously thought, or may have occurred preferentially during earlier Quaternary glacial periods.     

THE LAST GLACIATION OF SHETLAND, NORTH ATLANTIC

Evidence relating to the extent, dynamics, and relative chronology of the last glaciation of the Shetland Islands, North Atlantic, is presented here, in an attempt to better illuminate some of the controversies that still surround the glacial history of the archipelago. We appraise previous interpretations and compare these earlier results with new evidence gleaned from the interpretation of a high resolution digital terrain model and from field reconnaissance. By employing a landsystems approach, we identify and describe three quite different assemblages of landscape features across the main islands of Mainland, Yell and Unst. Using the spatial interrelationship of these landsystems, an assessment of their constituent elements, and comparisons with similar features in other glaciated environments, we propose a simple model for the last glaciation of Shetland. During an early glacial phase, a coalescent British and Scandinavian ice sheet flowed approximately east to west across Shetland. The terrestrial land‐forms created by this ice sheet in the north of Shetland suggest that it had corridors of relatively fast‐flowing ice that were partially directed by bed topography, and that subsequent deglaciation was interrupted by at least one major stillstand. Evidence in the south of Shetland indicates the growth of a local ice cap of restricted extent that fed numerous radial outlet glaciers during, or after, ice‐sheet deglaciation. Whilst the absolute age of these three landsystems remains uncertain, these new geo‐morphological and palaeoglaciological insights reconcile many of the ideas of earlier workers, and allow wider speculation regarding the dynamics of the former British ice sheet.     

Rapid and rhythmic ice sheet fluctuations in western Scandinavia 15-40 Kya–a review

The onshore record of Middle to Late Weichselian sediments and glacial history in Norway indicates a succession of four major ice advances alternating with rapid, considerable ice recessions and interstadial conditions. During all the glacial advances the ice sheet expanded from onshore/inland positions to the shelf areas. The basis for visualizing these variations in glaciation curves constructed along nine transects from inland to shelf, and for interpretation of the palaeoclimatic history, is the regional Quaternary stratigraphy, more than 300 datings, fossil content and some palaeomagnetic data. The methods applied in recent years for AMS radiocarbon dating of glacial sediments with low organic carbon content have given promising results with respect to accuracy and precision, and the results of such datings were an important tool for our reconstructions and for timing of the ice oscillations. The rapid and rhythmic ice fluctuations, as reconstructed in our new model, have been fairly synchronous in most parts of Norway. Ice advances commenced and culminated at 40, 30-28, 24-21 and 18-15 (¹⁴C) Kya. We describe three intervening interstadials from inland sites: Hattfjelldal I, Hattfjelldal II and Trofors. Our stratigraphical record also includes many indications of high, pre-Holocene, relative sea levels, suggesting a considerable glacioisostatic depression of western Scandinavia during the interstadials. In our glaciation model we suggest that, in addition to precipitation, the mountainous fjord and valley topography, glacial isostasy and relative sea level changes were probably more important for the size of the glacial fluctuations than were air temperature changes.     

Material versus isobaric internal boundaries in the Earth and their influence on postglacial rebound

Most previous earth models used to calculate viscoelastic relaxation after the removal of the Late Pleistocene ice loads implicitly assume that there is no exchange of mass across the mantle density discontinuities on periods of tens of thousands of years (the material boundary formulation). In the present study, simple incompressible models are used to determine the Earth's behaviour in the case where the density discontinuity remains at a constant pressure rather than deforming with the material (the isobaric boundary formulation). The calculation of the movement of the boundary is more rigorous than in earlier studies and uses the local incremental pressure calculated at the depth of the boundary and allows for the vertical deformation caused by the change in volume as material changes phase. It is shown that the buoyancy modes associated with the density discontinuities decrease in strength and increase in relaxation time analogous to what results when the density contrast is reduced. Also, two viscoelastic modes arise from an isobaric boundary, which is also predicted when there is a contrast in rigidity or viscosity across a material boundary. The difference in predicted radial deformation between the isobaric boundary model and the material boundary model is largest for long-wavelength loads for which the material incremental pressure at depth is largest. If the isobaric boundary model is appropriate for the treatment of the mineral phase changes in the mantle on glacial rebound timescales, then previous inferences of the deep-mantle to shallow-mantle viscosity ratio based on large-scale deformation (spherical harmonic degree < 10) of the Earth and including data from the early part of the glacio-isostatic uplift are too small.     

Drumlin Formation Time: Evidence from Northern and Central Sweden

Large‐scale drumlins occur abundantly throughout central and northern Sweden. Whereas many drumlins in the north are an integral part of a relict glacial landscape >100,000 years old, those to the south are generally interpreted as of last deglaciation age. Typically, the latter ones have not been overprinted by younger glacial landforms. Despite this apparent difference in formation history, drumlins in both regions have similar directional and morphological characteristics. A systematic analysis of >3000 drumlins in (i) areas within relict landscapes, (ii) areas with an ambiguous deglaciation age assignment, and (iii) areas within deglacial landscapes, indicates that these latter deglaciation drumlins differ clearly in both shape and size from drumlins in the other two types of landscapes. In addition, numerical modelling indicates that basal melting conditions, a prerequisite for drumlin formation, prevailed only for a very limited time over much of northern Sweden during the last deglaciation, but lasted for longer periods of time during earlier stages of the Weichselian. A reconnaissance radionuclide bedrock exposure date from the crag of a large drumlin in the relict landscape indicates that glacial erosion, and presumably drumlin formation, at this location predated Marine Isotope Stage 7. We conclude, therefore, that the large‐scale drumlins of central and northern Sweden did not form during the last deglaciation, or during any other specific ice flow event. Instead, we suggest that they were formed by successive phases of erosion and deposition by ice sheets of similar magnitude and configuration.