The Fennoscandian uplift and glacial isostasy

Various hypotheses have been put forth in relation to the land uplift of Fennoscandia, which is well documented both by geological and geodetic observations. Most modern authors attribute the present uplift to an isostatic rebound of the earth after the last deglaciation. Recent information on the gravity field, both from the satellite data and land survey measurements is examined to ascertain whether the Fennoscandian uplift is associated with a gravity minimum and a mass deficit. Free air anomalies correlate well with the central area of uplift and predict a remaining uplift of about 100 m. Results of secular gravity measurements are inconclusive. Seismicity of Fennoscandia does not show a close association with the area of maximum uplift. Different rheological models proposed for the mantle below the Fennoscandian shield are reviewed and it is shown that the available data on the rates of uplift for the last 9000 years are more compatible with a low-viscosity (10²⁰ P) asthenosphere of 100–200 km thickness.     

The present rate of uplift of Fennoscandia implies a low-viscosity asthenosphere

The rate of displacement in Fennoscandia has been intensively discussed for many years. It is now widely accepted to be an isostatic response of the glacial history of the area. The Earth's present response to deglaciation in Fennoscandia is simulated using a three-dimensional (3D), viscoelastic model in which the asthenosphere and mantle viscosity are allowed to vary so that the maximum rate of present uplift matches its observed value. The deglaciation history is considered to be known, and the C¹⁴-datings are converted to sidereal years. The pattern of the present uplift gives a firm match with the observed data when a low-viscosity asthenosphere is introduced. Assuming a 15,000 years load cycle, i.e. the glacier was applied to the surface for 15,000 years before the melting started, the best fitting earth viscosity model is a 10²⁴ Nm lithosphere overlying a 75 km-thick 2.0 × 10¹⁹ Pas asthenosphere and a 1.2 × 10²¹ Pas mantle. The simulations suggest a remaining maximum uplift of 40 m.     

Late Quaternary glacio- and hydro-isostasy, on a layered earth

Isostatic response of the Earth to changes in Quaternary Times of ice and water loads is partly elastic, and partly involves viscous mantle flow. The relaxation spectrum of the Earth, critical for estimation of the mantle flow component, is estimated from published determinations of Fennoscandian and Laurentide rebound, and of the nontidal acceleration of the Earth's rotation. The spectrum is consistent with an asthenosphere viscosity around 10²¹P, and a viscosity around 10²³P below 400 km depth. Calculation of relaxation effects is done by convoluting the load history with the response function in spherical harmonics for global effects, and in rectangular or cylindrical transforms for smaller regional effects. Broad-scale deformation of the globe, resulting from the last deglaciation and sea level rise, is calculated to have involved an average depression of ocean basins of about 8 m, and mean upward movement of continents of about 16 m, relative to the center of the Earth, in the last 7000 yr. Deflection in the ocean margin “hinge zone” varies with continental shelf geometry and rigidity of the underlying lithosphere: predictions are made for different model cases. The computational methods is checked by predicting Fennoscandian and Laurentide postglacial warping, from published estimates of icecap histories, with good results. The depth variations of shorelines formed around 17,000 BP (e.g., North America, 90–130 m; Australia, 130–170 m), are largely explainable in terms of combined elastic and relaxation isostasy. Differences between Holocene eustatic records from oceanic islands (Micronesia, Bermuda), and continental coasts (eastern North America, Australia), are largely but not entirely explained in the same terms.     

Uplift and active tectonics of southern Albania inferred from incision of alluvial terraces

In Albania, the Osum and Vjoje rivers cross the active graben system and the active frontal thrust system of the Albanides. The effects of climatic and geodynamic forcing on the development of these two rivers were investigated by the means of field mapping, topographic surveying and absolute exposure-age dating. We established the chronology of terraces abandonment from the compilation of new dating (¹⁴C and in situ produced ¹⁰Be) and previously published data. We identified nine fluvial terraces units developed since Marine Isotope Stage 6 up to historic times. From this reconstituted history, we quantified the vertical uplift on a time scale shorter than the glacial climatic cycle. Regional bulging produces an overall increase of the incision rate from the west to the east that reaches a maximum value of 2.8 m/ka in the hinterland. Local pulses of incision are generated by activation of normal faults. The most active faults have a SW–NE trend and a vertical slip rate ranging from 1.8 to 2.2 m/ka. This study outlines the geodynamic control in the development of rivers flowing through the Albanides on the scale of 10³–10⁵ka.     

Structure of the thermal lithosphere and asthenosphere beneath oceans and continents

The lithosphere and asthenosphere make up a common geodynamic system but are characterized by different physical parameters. The former has a temperature of 1200–1300°C, a density of 3.3 g/cm³, and a viscosity of 10²² poise, while the latter has a density of 3.23 g/cm³, a viscosity in the range 10²¹-10^18–19 poise, and a temperature from 1200–1300°C to 1600–1700°C. The asthenosphere is distinguished by a great variability of its physical state in the lateral and vertical directions. This circumstance necessitates the recognition of the different types of the asthenosphere: seismic (LVZ zone), electrical, thermal, and seismological. The structure and the physical state of the thermal asthenosphere is considered in this paper on the basis of P-T parameters. Its state normally fits viscous Newtonian liquid beneath the continents and provides partial (5–20%) melting in spreading zones and along continental margins. No partial melting is detected beneath the main portion of the continents. The interaction between the asthenosphere and lithosphere is characterized by spatiotemporal migration of partial melting zones and asthenosphere upwelling, and such interaction determines the entire range of geodynamic processes from spreading and rifting to collision and vertical motions of different senses.     

Stress and viscosity in the asthenosphere

Stresses and effective viscosities in the asthenosphere to a depth of 400 km are calculated on the basis of Weertmans “temperature method” i.e., on relating viscosity to the ratio of the temperature to the melting point (=homologous temperature). Some oceanic and continental geotherms and two melting point—depth curves, the dry pyrolite solidus and the forsterite₉₀ melting curve are used for the conversion of the homologous temperature to the effective viscosity. Two creep laws are considered, the linear, grain-size-dependent Nabarro—Herring (NH) creep law, and a power creep law, in which the creep rate is proportional to the third power of the stress. A plate tectonic model yields creep rates of 2 · 10⁻¹⁴ s⁻¹ for the oceanic and 3 · 10⁻¹⁵ s⁻¹ for the continental asthenosphere. These values are held constant for the calculations and may be valid for regions inside plates.The dry pyrolite mantle model results in high homologous temperatures in the asthenosphere below oceans (0.9), very low stresses (a few bars and lower) and shows a low viscosity “layer” of about 200-km thickness. Below continental shields the homologous temperature has a maximum value of 0.73, stresses are around 5–20 bar and the low-viscosity region is thicker and less pronounced than in the oceanic case. The Fo₉₀ mantle model generally gives lower homologous temperatures (maximum value below oceans beside active ridges 0.75). The stresses in the asthenosphere beneath oceans vary from a few bars to about 50 bar and below continents to about 100 bar. The low-viscosity region seems to reach great depths without forming a “channel”. The Figs. 1 and 2 show the approximate viscosity—depth distribution for the two mantle models under study.Assuming a completely dry mantle and a mean grain size of 5 mm, power law creep will be the dominating creep process in the asthenosphere. However, grains may grow in a high-temperature—low-stress regime (i.e., below younger oceans), an effect which will further diminish the influence of NH creep. In the upper 100–150 km of the earth some fluid phases may affect considerably creep processes.     

On the P-wave velocity and plate-tectonics implications for the Tyrrhenian deep-earthquake zone

P-wave velocities in the Tyrrhenian mantle have been determined for the 230–480 km depth range. Analysis of P-wave travel times for a set of Tyrrhenian deep earthquakes gives a velocity-distribution law which shows different behaviours in the 230–300 km and 300–480 km depth intervals. For the first interval the velocity gradient is 0.64 · 10⁻² sec⁻¹ and for the second one it is 0.59 · 10⁻² sec⁻¹. At a depth of 300 km the velocity decreases rapidly from 8.75 to 8.43 km/sec.The results have been analyzed in the framework of a Tyrrhenian structural model characterized by a lithospheric slab dipping 55–60° in the WNW direction.It is also pointed out that the analysis of some geodynamic features of the slabs of Pacific island arcs carried out by Oliver et al. (1973) and Sleep (1973) can be applied to the Tyrrhenian mantle geodynamic features.     

Isostasie flexurale et érosion différentielle : modélisation numérique appliquée au seuil du PoitouFlexural isostasy and differential erosion: numerical modelling applied at the ‘Seuil du Poitou’

One uses the principle of flexural isostasy to represent the uplift of the ‘Seuil du Poitou’ in response to erosion. The model indicates that a different uplift of on the last five millions years resulted from an erosion paradoxically stronger in the crystalline basement than in the limestones bedrock. To cite this article: J.-C. Maurin, K. Renaud, C. R. Geoscience 334 (2002) 1149–1155.     

Holocene relative sea-level change, isostatic subsidence and the radial viscosity structure of the mantle of northwest Europe (Belgium, the Netherlands, Germany, southern North Sea)

A comprehensive observational database of Holocene relative sea-level (RSL) index points from northwest Europe (Belgium, the Netherlands, northwest Germany, southern North Sea) has been compiled in order to compare and reassess the data collected from the different countries/regions and by different workers on a common time–depth scale. RSL rise varies in magnitude and form between these regions, revealing a complex pattern of differential crustal movement which cannot be solely attributed to tectonic activity. It clearly contains a non-linear, glacio- and/or hydro-isostatic subsidence component, which is only small on the Belgian coastal plain but increases significantly to a value of ca 7.5 m relative to Belgium since 8 cal. ka BP along the northwest German coast. The subsidence is at least in part related to the Post-Glacial collapse of the so-called peripheral forebulge which developed around the Fennoscandian centre of ice loading during the Last Glacial Maximum. The RSL data have been compared to geodynamic Earth models in order to infer the radial viscosity structure of the Earth's mantle underneath NW Europe (lithosphere thickness, upper- and lower-mantle viscosity), and conversely to predict RSL in regions where we have only few observational data (e.g. in the southern North Sea). A very broad range of Earth parameters fit the Belgian RSL data, suggesting that glacial isostatic adjustment (GIA) only had a minor effect on Belgian crustal dynamics during and after the Last Ice Age. In contrast, a narrow range of Earth parameters define the southern North Sea region, reflecting the greater influence of GIA on these deeper/older samples. Modelled RSL data suggest that the zone of maximum forebulge subsidence runs in a relatively narrow, WNW–ESE trending band connecting the German federal state of Lower Saxony with the Dogger Bank area in the southern North Sea. Identification of the effects of local-scale factors such as past changes in tidal range or tectonic activity on the spatial and temporal variations of sea-level index points based on model-data comparisons is possible but is still complicated by the relatively large range of Earth model parameters fitting each RSL curve, emphasising the need for more high-quality observational data.     

Tectonics and global change - inferences from Late Cenozoic subsidence and uplift patterns in the Atlantic/ Mediterranean region

Quantitative subsidence analysis for a number of rifted basins in the northern Atlantic/Mediterranean region provides evidence for rapid phases of PlioQuaternary subsidence. The observed acceleration in tectonic subsidence occurs after a phase of general quiescence in subsidence in these basins and deviates from predictions of stretching models. The latter indicate a decay of subsidence with time after Mesozoic-Tertiary basin formatioh and a slow tectonic subsidence in Plio-Quaternary times. A possible explanation for the observed patterns of anomalous subsidence could be an increase in the level of intraplate compression in the northern Atlantic region. Intraplate stress changes in the Plio-Quaternary are related to the dynamics of African/ Eurasian collision processes and a reorganization of spreading directions in the AtlAntic, possibly reflecting a plate reorganization of global nature.
It seems that the Plio-Quaternary record reflects a period of increased levels of neotectonic activity, interplaying with periods of (de)glaciation. Stress-induced topography in the onshore parts of continental margins, coupled with the stress-induced subsidence in the offshore deeper parts of the basins, could have contributed to recent phases of uplift in Fennoscandia, augmenting the uplift induced by glacial unloading. Estimates of ice thicknesses are directly inferred from the observed uplift ignoring other driving mechanisms whereas topography plays a crucial role in the dynamics of glaciation. It is, therefore, important to quantdy the interplay of rapid tectonic uplift and subsidence phases with climatic effects during the Plio-Quaternary.     

A possible mechanism of inversion of the vertical lithosphere movements in the back-arc basins of the West Pacific region

We present results of the study of a three-layer tectonosphere model of the West Pacific Transition Zone based on modeling of a piecewise inhomogeneous medium caused by local density reduction of the asthenosphere, whose viscosity decreases due to fluid accumulation. We used the viscous liquid motion equation in the Stokes approximation. It was shown that the anomalous asthenosphere in the back-arc basins can move as a convective cell with an uprising flow in the zone of maximum density reduction and extension of the lithosphere above it. At the initial stages, this process causes formation of the central crustal uplift, which is transformed into a system of depressions as the asthenosphere viscosity decreases to values of the order of 4.0 × 10¹⁹ Pa s and lower. The modeling results satisfactorily agree with the reconstructions of the Cretaceous Cenozoic lithotectonic evolution of the Okhotsk Sea region.     

The past valley glacier network in the Himalayas and the Tibetan ice sheet during the last glacial period and its glacial-isostatic, eustatic and climatic consequences

Since 1973 new data were obtained on the maximum extent of glaciation in High Asia. Evidence for an ice sheet covering Tibet during the Last Glacial Period means a radical rethinking about glaciation in the Northern Hemisphere. The ice sheet's subtropical latitude, vast size (2.4 million km²) and high elevation (6000 m asl) are supposed to have resulted in a substantial, albedo-induced cooling of the Earth's atmosphere and the disruption of summer monsoon circulation. Moraines were found to reach down to 460 m asl on the southern flank of the Himalayas and to 2300 m asl on the northern slope of the Tibetan Plateau, in the Qilian Shan region. On the northern slopes of the Karakoram, Aghil and Kuen-Lun mountains, moraines occur as far down as 1900 m asl. In southern Tibet radiographic analyses of erratics suggest a former ice thickness of at least 1200 m. Glacial polish and roches moutonnées in the Himalayas and Karakoram suggest former glaciers as thick as 1200–2700 m. On the basis of this evidence, a 1100–1600 m lower equilibrium line (ELA) has been reconstructed, resulting in an ice sheet of 2.4 million km², covering almost all of Tibet. Radiometric ages, obtained by different methods, classify this glaciation as isotope stage 3–2 in age (Würmian = last glacial period). With the help of 13 climate measuring stations, radiation- and radiation balance measurements have been carried out between 3800 and 6650 m asl in Tibet. They indicate that the subtropical global radiation reaches its highest energies on the High Plateau, thus making Tibet today's most important heating surface of the atmosphere. At glacial times 70% of those energies were reflected into space by the snow and firn of the 2.4 million km² extended glacier area covering the upland. As a result, 32% of the entire global cooling during the ice ages, determined by the albedo, were brought about by this area — now the most significant cooling surface. The uplift of Tibet to a high altitude about 2.75 Ma ago, coincides with the commencement of the Quaternary Ice Ages. When the Plateau was lifted above the snowline (= ELA) and glaciated, this cooling effect gave rise to the global depression of the snowline and to the first Ice Age. The interglacial periods are explained by the glacial-isostatic lowering of Tibet by 650 m, having the effect that the initial Tibet ice – which had evoked the build-up of the much more extended lowland ices – could completely melt away in a period of positive radiation anomalies. The next ice age begins, when – because of the glacial-isostatic reverse uplift – the surface of the Plateau has again reached the snowline. This explains, why the orbital variations (Milankovic-theory) could only have a modifying effect on the Quaternary climate dynamic, but were not primarily time-giving: as long as Tibet does not glaciate automatically by rising above the snowline, the depression in temperature is not sufficient for initiating a worldwide ice age; if Tibet is glaciated, but not yet lowered isostatically, a warming-up by 4 °C might be able to cause an important loss in surface but no deglaciation, so that its cooling effect remains in a maximum intensity. Only a glaciation of the Plateau lowered by isostasy, can be removed through a sufficiently strong warming phase, so that interglacial climate conditions are prevailing until a renewed uplift of Tibet sets in up to the altitude of glaciation.An average ice thickness for all of Tibet of approximately 1000 m would imply that 2.2 million km³ of water were stored in the Tibetan ice sheet. This would correspond to a lowering in sea level of about 5.4 m.     

The Central Scandinavian Dolerite Group—Protracted hotspot activity or back-arc magmatism?: Constraints from U–Pb baddeleyite geochronology and Hf isotopic data

The Central Scandinavian Dolerite Group (CSDG) occurs in five separate complexes in central Sweden and SW Finland. U–Pb baddeleyite ages of dolerite dikes and sills fall into three age intervals: 1264–1271 (the Dalarna complex), 1256–1259 (the Västerbotten-Ulvö-Satakunta complexes) and 1247 Ma (the Jämtland complex). Timing and spatial distribution of CSDG are unlike expressions of the voluminous and short-lived magmatism which characterises plume-associated large igneous provinces (LIPs). Protracted mafic magmatism in association with mantle plume tail (hotspot) activity beneath the Fennoscandian lithosphere or discrete events of extension behind an active margin (subduction) are considered more plausible tectonic settings. Both settings are consistent with timing, relative magma volumes between complexes and vertical ascent of individual magma pulses through the crust, as inferred from seismic sections [Korja, A., Heikkinen, P., Aaro, S., 2001. Crustal structure of the northern Baltic Sea palaeorift. Teconophysics 331, 341–358]. In the hotspot model, the lack of a linear track of intrusions can be explained by an almost stationary position of Fennoscandia relative to the hotspot, in agreement with palaeomagnetic data [Elming, S.-Å., Mattsson, H., 2001. Post Jotnian basic intrusion in the Fennoscandian Shield, and the break up of Baltica from Laurentia: a palaeomagnetic and AMS study. Precambrian Res. 108, 215–236]). Together with geological evidence, dolerite sill complexes and dike swarms in Labrador (Canada), S Greenland and central Scandinavia in the range 1234–1284 Ma are best explained by long-lived subduction along a continuous Laurentia-Baltica margin preceding Rodinia formation. There is no support for the hypothesis that CSDG was fed by magma derived from a distal mantle plume located between Baltica and Greenland and, hence, for rifting between the cratons at 1.26 Ga.The epsilon-Hf in various members of the CSDG varies between 4.7 and 10.3, which are overall higher than both older and younger Mesoproterozoic mafic intrusions in central Fennoscandia. Magma generated from a hotspot mantle source that was mixed to highly variable degrees with an enriched subcontinental lithospheric mantle could account for the wide range in Hf isotope composition. In the course of Hf isotope development work during this project we have analysed four fragments of the Geostandard 91500 reference zircon and after evaluating the existing ICPMS and TIMS data we calculate a mean ¹⁷⁶Hf/¹⁷⁷Hf value of 0.282303 ± 0.000003 (2σ).     

Low-frequency shear and structural relaxation in rhyolite melt

The frequency-dependence of the shear viscosity and modulus of rhyolite melt has been determined over a 10–14 log₁₀ Pa s viscosity range and a 0.03–63 rad s^–1 angular frequency range. The frequency-dependent viscosity determined at high frequencies is 5 orders of magnitude lower than the Newtonian viscosity. At lower frequencies, a frequency-independent viscosity identical with the Newtonian viscosity is observed. The measured shear modulus increases from zero to 30.5 ± 2.5 GPa with increasing frequency. The viscoelastic regime consists of a maximum in viscous loss centered on the Max-well relaxation time. The width and height of the loss modulus as a function of frequency is inconsistent with a single relaxation-time. The frequency-dependent shear modulus is best described by a distribution of relaxation-times with a sharp cutoff at times slightly longer than the Maxwell relaxation time, and a long tail at shorter times extending up to 5 orders of magnitude less than the Maxwell relaxation time. This distribution of relaxation-times is in contrast with the single-relaxation-time behavior observed in low viscosity silicate melts.     

Exposure ages from relict lateral moraines overridden by the Fennoscandian ice sheet

Lateral moraines constructed along west to east sloping outlet glaciers from mountain centred, pre-last glacial maximum (LGM) ice fields of limited extent remain largely preserved in the northern Swedish landscape despite overriding by continental ice sheets, most recently during the last glacial. From field evidence, including geomorphological relationships and a detailed weathering profile including a buried soil, we have identified seven such lateral moraines that were overridden by the expansion and growth of the Fennoscandian ice sheet. Cosmogenic ¹⁰Be and ²⁶Al exposure ages of 19 boulders from the crests of these moraines, combined with the field evidence, are correlated to episodes of moraine stabilisation, Pleistocene surface weathering, and glacial overriding. The last deglaciation event dominates the exposure ages, with ¹⁰Be and ²⁶Al data derived from 15 moraine boulders indicating regional deglaciation 9600 ± 200 yr ago. This is the most robust numerical age for the final deglaciation of the Fennoscandian ice sheet. The older apparent exposure ages of the remaining boulders (14,600-26,400 yr) can be explained by cosmogenic nuclide inheritance from previous exposure of the moraine crests during the last glacial cycle. Their potential exposure history, based on local glacial chronologies, indicates that the current moraine morphologies formed at the latest during marine oxygen isotope stage 5. Although numerous deglaciation ages were obtained, this study demonstrates that numerical ages need to be treated with caution and assessed in light of the geomorphological evidence indicating moraines are not necessarily formed by the event that dominates the cosmogenic nuclide data.     

Tectonics and crustal structures related to Bhuj earthquake of January 26, 2001: based on gravity and magnetic surveys constrained from seismic and seismological studies

Gravity and magnetic data of the Kachchh basin and surrounding regions have delineated major E–W and NW–SE oriented lineaments and faults, which are even extending up to plate boundaries in the north Arabian Sea and western boundary of the Indian plate, respectively. The epicentral zone of Bhuj earthquake and its aftershocks is located over the junction of Rann of Kachchh and median uplifts viz. Kachchh mainland and Wagad uplifts, which are separated by thrust faults. Gravity data with constraints from the results of the seismic studies along a profile suggest that the basement is uplifted towards the north along thrust faults dipping 40–60° south. Similarly gravity and magnetic modeling along a profile across Wagad uplift suggest south dipping (50–60°) basement contacts separating rocks of high susceptibility and density towards the north. One of these contacts coincides with the fault plane of the Bhuj earthquake as inferred from seismological studies and its projection on the surface coincides with the E–W oriented north Wagad thrust fault. A circular gravity high in contact with the fault in northern part of the Wagad uplift along with high amplitude magnetic anomaly suggests plug type mafic intrusive in this region. Several such gravity anomalies are observed over the island belt in the Rann of Kachchh indicating their association with mafic intrusions. The contact of these intrusives with the country rock demarcates shallow crustal inhomogeneities, which provides excellent sites for the accumulation of regional stress. A regional gravity anomaly map based on the concept of isostasy presents two centers of gravity lows of −11 to −13 mGal (10⁻⁵ m/s²) representing mass deficiency in the epicentral region. Their best-fit model constrained from the receiver function analysis and seismic refraction studies suggest crustal root of 7–8 km (deep crustal inhomogeneity) under them for a standard density contrast of −400 kg/m³. It is, therefore, suggested that significant amount of stress get concentrated in this region due to (a) buoyant crustal root, (b) regional stress due to plate tectonic forces, and (c) mafic intrusives as stress concentrators and the same might be responsible for the frequent and large magnitude earthquakes in this region including the Bhuj earthquake of January 26, 2001.     

Petrogenesis of the Neogene alkaline volcanics with implications for post-collisional lithospheric thinning of the Eastern Pontides, NE Turkey

Whole-rock geochemistry, Sr–Nd–Pb isotopes and K–Ar data are reported for alkaline samples collected from the Neogene alkaline volcanics (NAVs) in the Eastern Pontides, northeastern Turkey, in order to investigate their source and petrogenesis and geodynamic evaluation of the region. The NAVs were made of three groups that comprise of basanite–tephrite (feldspar-free; Group A), tephrite–tephriphonolite (feldspar and feldspathoid-bearing; Group B) and alkaline basalt–rhyolite (feldspathoid-free; Group C) series. These rocks cover a broad compositional range from silica-undersaturated to silica-oversaturated types, almost all of which are potassic in character. They show enrichment of LREE and LILE and depletion of HFSE, without a Eu anomaly in most of the mafic samples. Textural features and calculated pressures based on the Cpx-barometer in each series indicate that the alkaline magma equilibrated at shallow crustal depths under a pressure of about 3–4.5 kbar and approximating a crystallization depth of 9–14 km. The NAVs are slightly depleted in isotopic composition, with respect to ⁸⁷Sr/⁸⁶Sr (ranging from 0.705018 to 0.705643) and ¹⁴³Nd/¹⁴⁴Nd (ranging from 0.512662 to 0.512714) that indicate young Nd model ages (0.51–059 Ga). This may indicate that the parent melts tapped a homogeneous and young lithospheric mantle source which was metasomatized by subduction-derived sediments during the Late Mesozoic. Pb isotopic compositions (²⁰⁶Pb/²⁰⁴Pb = 18.85–18.95; ²⁰⁷Pb/²⁰⁴Pb = 15.60–15.74; ²⁰⁸Pb/²⁰⁴Pb = 38.82–39.25) may also be consistent with a model for an enriched subcontinental lithospheric mantle source. Lithospheric thinning and resultant upwelling of asthenosphere induced by lithospheric delamination may have favoured partial melting of chemically enriched, young lithospheric mantle beneath the Eastern Pontides. Then, the melt subsequently underwent a fractional crystallization process along with or without minor amounts of crustal assimilation, generating a wide variety of rock types in a post-collision extensional regime in the Eastern Pontides during the Neogene.     

Tectonic significance of the young mineral dates and the rates of cooling and uplift in the Himalaya

More than two-thirds of the published K-Ar, Rb-Sr and fission-track mineral dates from the Himalaya lie in the 5–75 m.y. range as a result of metamorphic overprint, uplift and cooling during the Late Cretaceous—Tertiary Himalayan orogeny. In contrast, the few but almost invariably old, Rb-Sr whole-rock ages reveal pre-Tertiary magmaticmetamorphic events.The pattern of distribution of these young dates vis-á-vis geological evidence reveals three phases, of the Himalayan orogeny, viz.: (1) folding and metamorphism (50–75 m.y.); (2) uplift (25–40 m.y.); and (3) major uplift, thrusting, formation of nappe structures, mylonitization and regional retrogression. The maximum concentration of dates in the 10–25 m.y. period marks this paroxysmal phase of the Himalayan orogeny.The Rb-Sr dates of co-existing muscovites and biotites have been used to compute the rates of cooling and uplift. Thus, slow cooling at the rate of about 4°C m.y.⁻¹ from 50 to 25 m.y. and rapid cooling at the rate of 19°-21°C m.y.⁻ from 25 m.y. to present have been inferred. The high rate of cooling over the past 25 m.y. is the result of major uplift at the rate of 0.7–0.8 mm yr⁻¹, which is in conformity with the current rate of uplift obtained from geodetic survey.     

Glacial rebound and relative sea levels in Europe from tide-gauge records

Relative sea levels recorded by tide gauges during the past century in northern Europe are dominated by isostatic readjustment of the land following the latest deglaciation of Scandinavia and Scotland. Maximum relative uplift of the land is centered near the northern Gulf of Bothnia (at a rate of 6–7 mm/yr), with a smaller secondary maximum over Scotland (also at a rate of 6–7 mm/yr). Although there probably is a relaxing peripheral bulge surrounding the regions of maximum uplift, such a former bulge is poorly defined by coastal tide gauges; in the North Sea evidence for sinking of a former peripheral bulge of glacial origin is complicated by post-Carboniferous basin deepening with sediment loading and possible rejuvenation associated with glaciation. Other data (gravity, radiocarbon, geomorphology) support the interpretation that glacial isostasy controls the structure of relative sea-level change. Included in this pattern of relative rise of land is a eustatic signal that biases the estimates of glacial rebound. Such a eustatic signal could not be isolated from the isostatic signal using the present data, but glacial isostasy clearly is a major control for relative sea levels of the region.

Absence of significant higher frequency (2–50 yr) cycles in mean annual sea levels of northern Europe reflects the complex hydrologic/oceanographic forces to which sea levels must respond. Whereas other coastal regions show significant higher frequency peaks in the energy spectra of relative sea levels, the many marginal seas in northern Europe preclude a clear relation between hydrologic/oceanographic forcing and relative sea levels, although this relation must exist on a more local scale.     

Cosmogenic nuclide chronology of pre-last glacial maximum moraines at Lago Buenos Aires, 46°S, Argentina

At Lago Buenos Aires, Argentina, ¹⁰Be, ²⁶Al, and ⁴⁰Ar/³⁹Ar ages range from 190,000 to 109,000 yr for two moraines deposited prior to the last glaciation, 23,000–16,000 yr ago. Two approaches, maximum boulder ages assuming no erosion, and the average age of all boulders and an erosion rate of 1.4 mm/10³ yr, both yield a common estimate age of 150,000–140,000 yr for the two moraines. The erosion rate estimate derives from ¹⁰Be and ²⁶Al concentrations in old erratics, deposited on moraines that are >760,000 yr old on the basis of interbedded ⁴⁰Ar/³⁹Ar dated lavas. The new cosmogenic ages indicate that a major glaciation during marine oxygen isotope stage 6 occurred in the mid-latitude Andes. The next five youngest moraines correspond to stage 2. There is no preserved record of a glacial advance during stage 4. The distribution of dated boulders and their ages suggest that at least one major glaciation occurred between 760,000 and >200,000 yr ago. The mid-latitude Patagonian glacial record, which is well preserved because of low erosion rates, indicates that during the last two glacial cycles major glaciations in the southern Andes have been in phase with growth and decay of Northern Hemisphere ice sheets, especially at the 100,000 yr periodicity. Thus, glacial maxima are global in nature and are ultimately paced by small changes in Northern Hemisphere insolation.