The surface geometry of the Last Glacial Maximum ice sheet in the Andøya-Skånland region, northern Norway, constrained by surface exposure dating and clay mineralogy

Blockfields, weathering boundaries and marginal moraines have been mapped along a longitudinal transect from northern Andøya to Skånland in northern Norway. The degree of rock-surface weathering above and below glacial trimlines, clay-mineral assemblages and surface exposure dating based on in situ cosmogenic ¹⁰Be have been used to reconstruct the vertical dimensions and timing of the Last Glacial Maximum (LGM) of the Scandinavian Ice Sheet in this region. The cosmogenic exposure dates suggest that the lower blockfield boundary/trimline along the Andøya-Skånland transect represents the upper limit of the Late Weichselian ice sheet, with an average surface gradient of c . 9.5 m/km. The surface exposure dates from Andøya pre-date the LGM, suggesting that the LGM ice sheet did not reach mountain plateaux at northwest Andøya. The results thus support evidence from lake sediment records that the northern tip of Andøya was not covered by the Scandinavian Ice Sheet during the LGM.     

A preliminary history of Holocene colluvial (debris-flow) activity,Leirdalen, Jotunheimen,Norway

A stratigraphic succession of alternating peat and minerogenic sediments at the foot of a steep mountain slope provides the basis for the reconstruction of a preliminary colluvial history from the alpine zone of Jotunheimen, southern Norway. Layers of silty sand and sandy silt, typically 5–10 cm thick and interpreted as distal debris-flow facies, are separated by layers of peat that have been radiocarbon dated. Deposition from at least 7500 to about 3800 ¹⁴C yr BP of predominantly minerogenic material suggests relatively infrequent but large-magnitude debris-flow events in an environment warmer and/or drier than today. Particularly low colluvial activity between about 6500 and 3900 ¹⁴C yr BP was terminated by a succession of major debris-flow events between about 3800 and 3400 ¹⁴C yr BP. Unhumified peats indicative of higher water tables, separate six debris-flows that occurred between about 3300 and 2300 ¹⁴C yr BP and signify a continuing high frequency of colluvial activity. Uninterrupted peat accumulation between about 2400 and 1600 ¹⁴C yr BP indicates reduced debris-flow activity; subsequent renewed activity appears to have culuminated in the ‘Little Ice Age’ after about 600 ¹⁴C yr BP. This pattern of colluvial deposition demonstrates a long history of natural Holocene low-alpine landscape instability, suggests an increase first in the magnitude and then in the frequency of debris-flow activity coincident with late Holocene climatic deterioration, and points to the potential of debris-flow records as a unique source of palaeoclimatic information related to extreme rainfall events. © 1997 John Wiley & Sons, Ltd.     

Nunataks of the last ice sheet in northwest Scotland

High-level weathering limits separating ice-scoured topography from an upper zone of frost-weathered detritus were identified on 17 mountains in NW Scotland at altitudes of <600 m to< 900 m, and a further 6 peaks were found to support evidence of ice scouring to summit level. Weathering limits are most clearly defined on Torridon Sandstone, which is resistant to frost shattering, but can also be mapped on Cambrian Quartzite, Lewisian Gneiss and Moine Schist. Contrasts in degree of rock surface weathering above and below the weathering limits were evaluated using measurements of joint depth and rock surface hardness, and through X-ray diffraction analyses of clay mineral assemblages. The results indicate significantly more advanced rock and soil weathering above the weathering limits. Widespread persistence of gibbsite above the weathering limits suggests that they represent the upper limit of Late Devensian glacial erosion, and the regularity of the decline in weathering limit altitude along former flowlines eliminates the possibility that it represents a former thermal boundary between protective cold-based and erosive warm-based ice. The weathering limits are therefore interpreted as periglacial trimlines defining the maximum surface altitude of the last ice sheet around former nunataks. Calculated basal shear stresses of 50–78 kPa are consistent with this interpretation. The altitude of the trimlines implies that the former ice shed lay at 900–930 m in the Fannich Mountains and descended gently northwards, and that the ice surface descended NW from the ice shed to >500 m over the extreme NW tip of Scotland and to 700–730 m at the head of Little Loch Broom.     

NON-LINEAR AVO INVERSION FOR A STACK OF ANELASTIC LAYERS1

Parameters in a stack of homogeneous anelastic layers are estimated from seismic data, using the amplitude versus offset (AVO) variations and the travel-times. The unknown parameters in each layer are the layer thickness, the P-wave velocity, the S-wave velocity, the density and the quality factor. Dynamic ray tracing is used to solve the forward problem. Multiple reflections are included, but wave-mode conversions are not considered. The S-wave velocities are estimated from the PP reflection and transmission coefficients. The inverse problem is solved using a stabilized least-squares procedure. The Gauss-Newton approximation to the Hessian matrix is used, and the derivatives of the dynamic ray-tracing equation are calculated analytically for each iteration. A conventional velocity analysis, the common mid-point (CMP) stack and a set of CMP gathers are used to identify the number of layers and to establish initial estimates for the P-wave velocities and the layer thicknesses. The inversion is carried out globally for all parameters simultaneously or by a stepwise approach where a smaller number of parameters is considered in each step. We discuss several practical problems related to inversion of real data. The performance of the algorithm is tested on one synthetic and two real data sets. For the real data inversion, we explained up to 90% of the energy in the data. However, the reliability of the parameter estimates must at this stage be considered as uncertain.     

SEISMIC REFLECTION AMPLITUDES1

Seismic amplitude variations with offset contain information about the elastic parameters. Prestack amplitude analysis seeks to extract this information by using the variations of the reflection coefficients as functions of angle of incidence. Normally, an approximate formula is used for the reflection coefficients, and variations with offset of the geometrical spreading and the anelastic attenuation are often ignored. Using angle of incidence as the dependent variable is also computationally inefficient since the data are recorded as a function of offset. Improved approximations have been derived for the elastic reflection and transmission coefficients, the geometrical spreading and the complex travel-time (including anelastic attenuation). For a 1 D medium, these approximations are combined to produce seismic reflection amplitudes (P-wave, S-wave or converted wave) as a Taylor series in the offset coordinate. The coefficients of the Taylor series are computed directly from the parameters of the medium, without using the ray parameter. For primary reflected P-waves, dynamic ray tracing has been used to compute the offset variations of the transmission coefficients, the reflection coefficient, the geometrical spreading and the anelastic attenuation. The offset variation of the transmission factor is small, while the variations in the geometrical spreading, absorption and reflection coefficient are all significant. The new approximations have been used for seismic modelling without ray tracing. The amplitude was approximated by a fourth-order polynomial in offset, the traveltime by the normal square-root approximation and the absorption factor by a similar expression. This approximate modelling was compared to dynamic ray tracing, and the results are the same for zero offset and very close for offsets less than the reflector depth.     

Cirque glaciers as morphological evidence for a thin Younger Dryas ice sheet in east-central southern Norway

Three localities with marginal moraines deposited by former cirque glaciers are investigated in east-central southern Norway. The wet-based (erosive) cirque glaciers with aspects towards S-SW and N-NE are mapped at altitudes above 1100 m, and have a mean equilibrium-line altitude of 1275 m. With a suggested mean annual winter precipitation close to the average for the modern accumulation season (1 October-30 April) when the cirque glaciers existed, the mean air-temperature depression during the ablation season (1 May-30 September) is calculated to be 6–7°C lower than at present. The high-altitude cirques of central Rondane were still covered by ice when the low-altitude cirque glaciers developed in distal position for this massif in eastern Rondane and on isolated mountains. Hence, the cirque glaciers are suggested to have existed during the deglaciation after the Late Weichselian maximum, and most likely during the Younger Dryas (11000–10000 BP). The cirque glaciers indicate a downwasting ice-sheet surface well below an altitude of 1100 m prior to the Younger Dryas, and this supports a limited (small) vertical extent for the Late Weichselian ice sheet in this region. With the contemporaneous level for instantaneous glacierization (glaciation threshold) just below the highest elevated peaks in east-central southern Norway, this fits with the idea of a continuous downwasting of the Late Weichselian ice sheet since the 'first' nunataks appeared. The occurrence of the cirque glaciers indicates a multidomed Scandinavian ice-sheet geometry during the Late Weichselian.     

THE LAST ICE SHEET IN NORTH-WEST SCOTLAND: RECONSTRUCTION AND IMPLICATIONS

Recent models of the last Scottish ice sheet suggest that nunataks remained above the ice surface in areas peripheral to the main centres of accumulation. This proposition has been investigated on 140 mountains over an area of 10,000 km² in NW Scotland. Outside the limits of the later Loch Lomond Readvance in this area there is evidence for a single high-level weathering limit that separates glacially eroded terrain from higher areas of in situ frost debris. This limit occurs at altitudes ranging from 425 to 450 m in the Outer Hebrides to >950 m on the mainland, and is best developed on lithologies that resisted breakdown after ice-sheet downwastage. Interpretation of this weathering limit as a periglacial trimline cut by the last ice sheet at its maximum thickness is supported by: (1) joint-depth and Schmidt hammer measurements that indicate significantly more advanced rock breakdown above the weathering limit; (2) a much greater representation of gibbsite (a pre-Late Devensian weathering product) in the clay fraction of soils above the limit; (3) cosmogenic isotope dating of the exposure ages of rock outcrops above and below the limit; (4) the sharpness of the limit at some sites and its regular decline along former ice flowlines; and (5) shear stress calculations based on the inferred altitude and gradient of the former ice surface. Reconstruction of the ice surface based on trimline evidence indicates that the mainland ice shed lay near or slightly east of the present watershed and descended northwards from >900 m to ca. 550 m at the north coast. Independent dispersion centres fed broad ice streams that occupied major troughs. On Skye an ice dome >800 m deflected the northwestwards movement of mainland ice, but the mountains of Rum were over-ridden by mainland ice up to an altitude of ca. 700 m. The Outer Hebrides supported an independent ice cap that was confluent with mainland ice in the Minches. Extrapolation of the trimline evidence indicates that most reconstructions of ice extent are too conservative, and suggests that low-gradient ice streams extended across the Hebridean Shelf offshore. Wider implications of this research are: (1) that blockfields and other periglacial weathering covers are not all of the same age or significance, depending on the resistance of different lithologies to frost weathering; (2) that the contrasting degree of glacial modification in the Western and Eastern Highlands of Scotland may reflect a former cover of predominantly warm-based ice in the former and predominantly cold-based ice in the latter; and (3) that the approach and techniques developed in this study have potential application for constraining ice-sheet models, not only in areas peripheral to the main centres of ice accumulation in Britain and Ireland, but also in other mountain areas where nunataks protruded through warm-based Late Pleistocene ice masses.     

The porosity and permeability of chondritic meteorites and interplanetary dust particles

Abstract— We have investigated the porosity of a large number of chondritic interplanetary dust particles (IDPs) and meteorites by three techniques: standard liquid/gas flow techniques, a new, noninvasive ultrasonic technique, and image processing of backscattered images. The latter technique is obviously best-suited to sub-kilogram sized samples. We have also measured the gas and liquid permeabilities of some chondrites by two techniques: standard liquid/gas flow techniques, and a new, nondestructive pressure release technique. We find that chondritic IDPs have a somewhat bimodal porosity distribution. Peaks are present at 0 and 4% porosity; a tail then extends to 53%. Type 1–3 chondrite matrix porosities range up to 30%, with a peak at 2%. The bulk porosities for type 1–3 chondrites have the same approximate range as exhibited by the matrix, which indicates that other components of the bulk meteorites (including chondrules and aggregates) have the same average porosity as the matrix. These results reveal that the porosities of primitive materials at scales ranging from nanogram to kilogram are similar, which implies that similar accretion dynamics operated through 12 orders of size magnitude. Permeabilities of the investigated chondrites vary by several orders of magnitude, and there appears to be no simple dependence of permeability with degree of aqueous alteration, chondrite type or porosity.