Showing papers by "British Geological Survey published in 2003"

Groundwater: the processes and global significance of aquifer degradation

Abstract: The exploitation of groundwater resources for human use dates from the earliest civilizations, but massive resource development has been largely restricted to the past 50 years. Although global in scope, the emphasis of this paper is on groundwater-based economies in a developing nation context, where accelerated resource development has brought major social and economic benefits over the past 20 years. This results from groundwater's significant role in urban water supply and in rural livelihoods, including irrigated agriculture. However, little of the economic benefit of resource development has been reinvested in groundwater management, and concerns about aquifer degradation and resource sustainability began to arise. A general review, for a broad-based audience, is given of the mechanisms and significance of three semi-independent facets of aquifer degradation. These are (i) depletion of aquifer storage and its effects on groundwater availability, terrestrial and aquatic ecosystems; (ii) groundwater salinization arising from various different processes of induced hydraulic disturbance and soil fractionation; and (iii) vulnerability of aquifers to pollution from land-use and effluent discharge practices related to both urban development and agricultural intensification. Globally, data with which to assess the status of aquifer degradation are of questionable reliability, inadequate coverage and poor compilation. Recourse has to be made to 'type examples' and assumptions about the extension of similar hydrogeological settings likely to be experiencing similar conditions of groundwater demand and subsurface contaminant load. It is concluded that (i) aquifer degradation is much more than a localized problem because the sustainability of the resource base for much of the rapid socio-economic development of the second half of the twentieth century is threatened on quite a widespread geographical basis; and (ii) major (and long overdue) investments in groundwater resource and quality protection are urgently needed. These investments include appropriate institutional provisions, demand-side management, supply-side enhancement and pollution control.

Common-Pb corrected in situ U–Pb accessory mineral geochronology by LA-MC-ICP-MS

Abstract: LA-MC-ICP-MS is shown to be a rapid, precise and accurate method for determination of U–Pb ages of accessory minerals. For the protocol described, total analysis time is <3 min with a main acquisition sequence of only 30 s. Using a raster ablation protocol, within-run inter-element fractionation can be effectively eliminated and an external ablation standard used to quantify an overall error for the analysis. Reproducibilities of 206Pb/238U = ca. 3% and 207Pb/206Pb = <1% (2 σ) are achieved, with the resulting age accurate to within 1% as determined using in-house samples previously characterised by TIMS. A key control on the Pb/Pb reproducibility is shown to be the size of the 207Pb peak and an error propagation curve is determined for the accurate representation of this data. Propagation of these errors allows each individual sample analysis to be considered a stand-alone result, removing the need for statistical averaging of multiple data points. Simultaneous collection of flat-topped peaks enables precise measurement and correction of isobaric interference from 204Hg and a procedure for the consistent correction of common-Pb using 204Pb is described. Determination and correction of the common-Pb component is shown to be critical to the reliable interpretation of the data for certain minerals including those phases where a correction is often deemed unnecessary. Combined with time-resolved analysis of the data, this allows the Pb-loss history and nature of discordance within individual crystal domains to be ascertained. Successful analyses of zircons using a non-matrix matched (monazite) standard are also demonstrated suggesting that particle size distribution, ionisation efficiency and plasma loading, are more important issues in controlling inter-element fractionation in the plasma than absolute matrix matching.

Landslide tsunami case studies using a Boussinesq model and a fully nonlinear tsunami generation model

Abstract: . Case studies of landslide tsunamis require integration of marine geology data and interpretations into numerical simulations of tsunami attack. Many landslide tsunami generation and propagation models have been proposed in recent time, further motivated by the 1998 Papua New Guinea event. However, few of these models have proven capable of integrating the best available marine geology data and interpretations into successful case studies that reproduce all available tsunami observations and records. We show that nonlinear and dispersive tsunami propagation models may be necessary for many landslide tsunami case studies. GEOWAVE is a comprehensive tsunami simulation model formed in part by combining the Tsunami Open and Progressive Initial Conditions System (TOPICS) with the fully non-linear Boussinesq water wave model FUNWAVE. TOPICS uses curve fits of numerical results from a fully nonlinear potential flow model to provide approximate landslide tsunami sources for tsunami propagation models, based on marine geology data and interpretations. In this work, we validate GEOWAVE with successful case studies of the 1946 Unimak, Alaska, the 1994 Skagway, Alaska, and the 1998 Papua New Guinea events. GEOWAVE simulates accurate runup and inundation at the same time, with no additional user interference or effort, using a slot technique. Wave breaking, if it occurs during shoaling or runup, is also accounted for with a dissipative breaking model acting on the wave front. The success of our case studies depends on the combination of accurate tsunami sources and an advanced tsunami propagation and inundation model.

The O and H stable isotope composition of freshwaters in the British Isles. 2. Surface waters and groundwater

Abstract: . The utility of stable isotopes as tracers of the water molecule has a long pedigree. The study reported here is part of an attempt to establish a comprehensive isotopic "baseline" for the British Isles as background data for a range of applications. Part 1 of this study (Darling and Talbot, 2003) considered the isotopic composition of rainfall in Britain and Ireland. The present paper is concerned with the composition of surface waters and groundwater. In isotopic terms, surface waters (other than some upland streams) are poorly characterised in the British Isles; their potential variability has yet to be widely used as an aid in hydrological research. In what may be the first study of a major British river, a monthly isotopic record of the upper River Thames during 1998 was obtained. This shows high damping of the isotopic variation compared to that in rainfall over most of the year, though significant fluctuations were seen for the autumn months. Smaller rivers such as the Stour and Darent show a more subdued response to the balance between runoff and baseflow. The relationship between the isotopic composition of rainfall and groundwater is also considered. From a limited database, it appears that whereas Chalk groundwater is a representative mixture of weighted average annual rainfall, for Triassic sandstone groundwater there is a seasonal selection of rainfall biased towards isotopically-depleted winter recharge. This may be primarily the result of physical differences between the infiltration characteristics of rock types, though other factors (vegetation, glacial history) could be involved. In the main, however, groundwaters appear to be representative of bulk rainfall within an error band of 0.5‰ δ18O. Contour maps of the δ18O and δ2H content of recent groundwaters in the British Isles show a fundamental SW-NE depletion effect modified by topography. The range of measured values, while much smaller than those for rainfall, still covers some ‰ for δ18O and 30‰ for δ2H. Over lowland areas the "altitude effect" is of little significance, but in upland areas is consistent with a range of –0.2 to –0.3‰ per 100 m increase in altitude. Groundwaters dating from the late Pleistocene are usually modified in δ18O and δ2H owing to the effects of climate change on the isotopic composition of rainfall and thus of recharge. Contour maps of isotopic variability prior to 10 ka BP, based on the relatively limited information available from the British Isles, allow a first comparison between groundwaters now and at the end of the last Ice Age. The position of the British Isles in the context of the stable isotope systematics of NW Europe is reviewed briefly. Keywords: Stable isotopes, surfacewaters, groundwater, British Isles

Exhumation of the Main Central Thrust from Lower Crustal Depths, Eastern Bhutan Himalaya

Abstract: Geothermometry and mineral assemblages show an increase of temperature structurally upwards across the Main Central Thrust (MCT); however, peak metamorphic pressures are similar across the boundary, and correspond to depths of 35-45 km. Garnet-bearing samples from the uppermost Lesser Himalayan sequence (LHS) yield metamorphic conditions of 650-675 ! C and 9-13 kbar. Staurolite-kyanite schists, about 30 m above the MCT, yield P-T conditions near 650 ! C, 8-10 kbar. Kyanite-bearing migmatites from the Greater Himalayan sequence (GHS) yield pressures of 10-14 kbar at 750-800 ! C. Top-to-the- south shearing is synchronous with, and postdates peak metamorphic mineral growth. Metamorphic monazite from a deformed and metamorphosed Proterozoic gneiss within the upper LHS yield U⁄Pb ages of 20-18 Ma. Staurolite-kyanite schists within the GHS, a few metres above the MCT, yield monazite ages of c. 22 ± 1 Ma. We interpret these ages to reflect that prograde meta- morphism and deformation within the Main Central Thrust Zone (MCTZ) was underway by c. 23 Ma. U⁄Pb crystallization ages of monazite and xenotime in a deformed kyanite-bearing leucogranite and kyanite-garnet migmatites about 2 km above the MCT suggest crystallization of partial melts at 18- 16 Ma. Higher in the hanging wall, south-verging shear bands filled with leucogranite and pegmatite yield U⁄Pb crystallization ages for monazite and xenotime of 14-15 Ma, and a 1-2 km thick leuco- granite sill is 13.4 ± 0.2 Ma. Thus, metamorphism, plutonism and deformation within the GHS con- tinued until at least 13 Ma. P-T conditions at this time are estimated to be 500-600 ! C and near 5 kbar. From these data we infer that the exhumation of the MCT zone from 35 to 45 km to around 18 km, occurred from 18 to 16 to c. 13 Ma, yielding an average exhumation rate of 3-9 mm year )1 . This process of exhumation may reflect the ductile extrusion (by channel flow) of the MCTZ from between the overlying Tibetan Plateau and the underthrusting Indian plate, coupled with rapid erosion.

Storage of fossil fuel-derived carbon dioxide beneath the surface of the earth

Abstract: ▪ Abstract Underground storage in porous and permeable reservoir rocks is a technically feasible way to dispose of industrial quantities of carbon dioxide such as are produced by a fossil fuel–fired power plant. All the necessary steps are commercially proven and in use today. Extensive, naturally occurring CO2 accumulations indicate that under favorable conditions CO2 can be retained in underground reservoirs for millions of years. Large-scale commercial underground CO2 sequestration has begun at the Sleipner West gas field in the North Sea. Some of the major issues to be addressed if this technology is to make an impact on CO2 emissions to the atmosphere are cost of CO2 capture, safety and security of storage, and public acceptability.

Modelling frequency-dependent seismic anisotropy in fluid-saturated rock with aligned fractures: implication of fracture size estimation from anisotropic measurements

Abstract: Measurements of seismic anisotropy in fractured rock are used at present to deduce information about the fracture orientation and the spatial distribution of fracture intensity. Analysis of the data is based upon equivalent-medium theories that describe the elastic response of a rock containing cracks or fractures in the long-wavelength limit. Conventional models assume frequency independence and cannot distinguish between microcracks and macrofractures. The latter, however, control the fluid flow in many subsurface reservoirs. Therefore, the fracture size is essential information for reservoir engineers. In this study we apply a new equivalent-medium theory that models frequency-dependent anisotropy and is sensitive to the length scale of fractures. The model considers velocity dispersion and attenuation due to a squirt-flow mechanism at two different scales: the grain scale (microcracks and equant matrix porosity) and formation-scale fractures. The theory is first tested and calibrated against published laboratory data. Then we present the analysis and modelling of frequency-dependent shear-wave splitting in multicomponent VSP data from a tight gas reservoir. We invert for fracture density and fracture size from the frequency dependence of the time delay between split shear waves. The derived fracture length matches independent observations from borehole data.

Petrogenesis of Cenozoic Basalts from Mongolia: Evidence for the Role of Asthenospheric versus Metasomatized Lithospheric Mantle Sources

Abstract: Diffuse Cenozoic volcanism in Mongolia forms part of a widespread tectono-magmatic province that extends from NE China to Lake Baikal, Siberia. Mafic lavas from the Gobi Altai, southern Mongolia (33 Ma) and Hangai, central Mongolia ( 70 km. Three isotopic end-members can explain the heterogeneity: (1) is similar to bulk silicate Earth with 206Pb/204Pb >17·8 and is asthenospheric; (2) is EM1-like, characterized by low 206Pb/204Pb (>17·062), and may represent mobilized ancient lithospheric mantle; (3) also lithospheric, is characterized by low 143Nd/144Nd (>0·512292) and shows similarities to EM2, although decoupling of isotopic systems suggests a complex enrichment process. The timing of lithospheric enrichment is unconstrained, but may be related to Mesozoic magmatic events and/or melts mobilized during the Cenozoic responding to higher than ambient potential temperature mantle. Published geophysical studies suggest anomalous material at the base of the lithospheric mantle; however, there is no evidence to suggest a high heat flux mantle plume. Volcanism is likely to occur where localized extensional conditions are favourable.

Source of the Dalradian Supergroup constrained by U-Pb dating of detrital zircon and implications for the East Laurentian margin

Abstract: Detrital zircons in siliciclastic units of the Dalradian Supergroup yield U–Pb ages ranging from 3.2 to 0.5 Ga. Detrital zircons from the sub-Grampian Group basement and the Grampian Group are predominantly of Palaeoproterozoic and Mesoproterozoic ages with Archaean grains absent or rare. In contrast, the overlying Appin, Argyll and Southern Highland groups contain a significant contribution of Archaean detrital zircon grains, some of which locally preserve evidence for late Palaeoproterozoic overgrowths dated at c . 1.8 Ga. In addition, on concordia plots Archaean grains are slightly discordant with a lower intercept at c . 1.8 Ga suggesting they were affected by a tectonothermal event at this time. Late Palaeoproterozoic and Mesoproterozoic grains also show evidence for overprinting by a tectonothermal event around 1.0–0.9 Ga. These tectonothermal events occurred in the source region before accumulation of the siliciclastic detritus. The analysed samples contained no definitive evidence for having been affected by a late Neoproterozoic (Knoydartian) event or of containing detritus derived from a source showing evidence for this event. The overall age range of detritus, combined with sparse palaeocurrent data, is consistent with derivation from the Laurentian foreland, especially the Labrador–Greenland region. Archaean detritus overlaps with that of Archaean cratons, notably the Superior, whereas Palaeoproterozoic detritus corresponds to the timing of suturing of Archaean cratons by a series of orogenic belts (Ketilidian–Makkovik, New Quebec, Nagssugtoqidian, Torngat belts). Mesoproterozoic detritus is consistent with derivation from the Grenville Orogen. The presence of a series of detrital age peaks in the late Palaeoproterozoic and early Mesoproterozoic (1.8–1.5 Ga), the paucity of mid-Mesoproterozoic detritus (1.4–1.2 Ga), and evidence for a tectonothermal event between 1.0 and 0.9 Ga is typical of the geological history of the Labrador–Greenland region of Laurentia.

Hydrogeochemical and isotopic evolution of water in the Complexe Terminal aquifer in the Algerian Sahara

Abstract: The hydrogeochemical and isotopic evolution of groundwaters in the Mio-Pliocene sands of the Complexe Terminal (CT) aquifer in central Algeria are described. The CT aquifer is located in the large sedimentary basin of the Great Oriental Erg. Down- gradient groundwater evolution is considered along the main representative aquifer cross section (south-north), from the southern recharge area (Tinrhert Plateau and Great Oriental Erg) over about 700 km. Groundwater mineralisation increases along the flow line, from 1.5 to 8gl �1 , primarily as a result of dissolution of evaporite minerals, as shown by Br/Cl and strontium isotope ratios. Trends in both major and trace elements demonstrate a progressive evolution along the flow path. Redox reac- tions are important and the persistence of oxidising conditions favours the increase in some trace elements (e.g. Cr) and also NO3 � , which reaches concentrations of 16.8 mg l �1 NO3-N. The range in 14 C, 0-8.4 pmc in the deeper groundwaters, corresponds with late Pleistocene recharge, although there then follows a hiatus in the data with no results in the range 10-20 pmc, interpreted as a gap in recharge coincident with hyper-arid but cool conditions across the Sahara; groundwater in the range 24.7-38.9 pmc signifies a distinct period of Holocene

Use of low frequencies for sub‐basalt imaging

Abstract: Many prospective passive ocean margins are covered by large areas of basalts. These basalts are often extremely heterogeneous and scatter the seismic energy of the conventional seismic reflection system so that it becomes difficult to obtain information on deeper reflectors. Since high frequencies are scattered more than low frequencies, we argue that the acquisition system for sub-basalt targets should be modified to emphasize the low frequencies, using much larger airguns, and towing the source and receivers at about 20 m depth. In the summer of 2001 we obtained seismic reflection data over basalt in the northeast Atlantic using a system modified to enhance the low-frequency energy. These new data show deep reflections that are not visible on lines shot in the same places with a conventional system.

The O and H stable isotope composition of freshwaters in the British Isles. 1. Rainfall

Abstract: . An understanding of the hydrological cycle in stable isotopic terms requires the characterisation of rainfall. This paper reviews existing and new data for the British Isles. Rainfall at the Wallingford (Oxfordshire) collection station was collected daily from November 1979 to October 1980. Large variations in isotopic content were noted, sometimes from day to day. Winter rainfall was similar to summer in amount, and only slightly depleted isotopically. Amount and temperature correlations with δ18O were generally low, only the autumn and winter temperature relationships being significant. A 20-year monthly dataset from 1982 to 2001 for Wallingford gives the following regression: δ2H = 7.0δ18O + 1.2, a slope somewhat below the world meteoric line but consistent with the those from other long-term stations in NW Europe. The data showed uncorrelated maxima and minima for each year, but rather more consistent amount-weighted averages. Although there is only a small difference in gradient between summer and winter rainfall values, when plotted against the month of the year there are clear changes in the values of both isotopes, and the δ2H-δ18O relationship as demonstrated by the d-excess parameter. The isotope-amount correlation is low but significant, with summer months appearing to be well-correlated when considered in terms of month of the year. On this same seasonal basis temperature has a strong correlation throughout the year, giving a positive δ18O-temperature relationship of 0.25 ‰ per °C change. The Wallingford monthly record is compared with data from Keyworth (Nottinghamshire) and the Valentia station of the GNIP (IAEA-WMO Global Network for Isotopes in Precipitation) in SW Ireland. While not large, differences between the stations are broadly attributable to the balance between maritime and continental influences. Over the period September 1981 to August 1982 the maximum number of monthly collection stations was operating across the British Isles. While a comparison of the sites serves mostly to illustrate the variability of British weather in space and time, there is clear isotopic evidence for the predominance of frontal rainfall in winter and convective rainfall in summer. The effect of altitude on isotopic content was measured within a high-relief stream catchment in Scotland. The best correlations occurred during winter, when an average relationship of approximately –0.30 ‰ δ18O per 100 m increase in altitude was observed. It is well established that rainfall isotopic composition changes in response to alterations in climate. However these changes are difficult to detect isotopically in the short term, even when the changes are indexed, e.g. in the form of the NAO (North Atlantic Oscillation). The brief duration of rainfall isotope records is a further hindrance; for the British Isles proxies such as tree-ring cellulose may have some value in extending the record back. Keywords: stable isotopes, rainfall,British Isles

Geostatistical analysis of arsenic concentration in groundwater in Bangladesh using disjunctive kriging

Abstract: The National Hydrochemical Survey of Bangladesh sampled the water from 3,534 tube wells for arsenic throughout most of Bangladesh. It showed that 27% of the shallow tube wells (less than 150 m deep) and 1% of the deep tube wells (more than 150 m deep) exceeded the Bangladesh standard for arsenic in drinking water (50 µg L−1). Statistical analyses revealed the main characteristics of the arsenic distribution. Concentrations ranged from less than the detection limit (≈0.5 µg L−1), to as much as 1,600 µg L−1, though with a very skewed distribution, and with spatial dependence extending to some 180 km. Disjunctive kriging was used to estimate concentrations of arsenic in the shallow ground water and to map the probability that the national limit for arsenic in drinking water was exceeded for most of the country (the Chittagong Hill Tracts and the southern coastal region were excluded). A clear regional pattern was identified, with large probabilities in the south of the country and small probabilities in much of the north including the Pleistocene Tracts. Using these probabilities, it was estimated that approximately 35 million people are exposed to arsenic concentrations in groundwater exceeding 50 µg L−1 and 57 million people are exposed to concentrations exceeding 10 µg L−1 (the WHO guideline value).

Geochemical distribution of arsenic in waters, sediments and weathered gold mineralized rocks from Iron Quadrangle, Brazil

Abstract: The Iron Quadrangle has been the scenery of the most important gold production in Brazil. It is estimated that during the three centuries of gold mining in the Iron Quadrangle, at least 390,000 t of arsenic was discharged into the drainage system. This study presents geochemical data for the three river basins in the region, with focus on surface water and stream sediment monitoring. Samples of primary and oxidized sulfide ores as well as of tailings and groundwater from the major gold mines were also studied. The highest As concentrations in water and stream sediments occur in the vicinity of mining areas. In surface water, up to 300 g As/l were found whereas the As contents in stream sediments were in the range of 20 to 4,000 mg/kg. The As3+/As5+ concentration ratios obtained for some water samples range from 1.10−1 to 4.10−2. The As mobility associated with ore-deposit weathering could be traced in some closed gold mines by observation of in-situ pyrite and arsenopyrite oxidation, precipitation of scorodite and gippsite, As adsorption onto goethite, and final liberation of As into underground and surface waters. This process is likely to produce large volumes of mine effluents containing total As and trivalent As up to 1,960 and 60 g/l, respectively. River sediments and tailings pile samples were submitted to a leaching procedure showing maximal arsenic release from 1 to 4% of the original total As in the samples. There are potential risks for As hazards in some areas induced by, for instance, the dispersion of old tailings by flooding, occupation of poisoned soils for settlements, and occasional consumption of contaminated surface and groundwater.

Sr isotope evidence for population movement within the Hebridean Norse community of NW Scotland.

Abstract: The excavation at Cnip, Isle of Lewis, Scotland, of the largest, and only known family cemetery from the early Norse period in the Hebrides, provided a unique opportunity to use Sr isotope analysis to examine the origins of people who may have been Norwegian Vikings. Sr isotope analysis permits direct investigation of a person9s place of origin rather than indirectly through acquired cultural and artefactual affiliations. Sr isotope data suggest that the Norse group at Cnip was of mixed origins. The majority were consistent with indigenous origins but two individuals, of middle-age and different sex, were immigrants. They were, however, not from Norway but were raised separately, most probably on Tertiary volcanic rocks (e.g. the Inner Hebrides or NE Ireland) or, for the female, on marine carbonate rocks.

Numerical simulation of wave propagation in media with discrete distributions of fractures: effects of fracture sizes and spatial distributions

Abstract: SUMMARY We model seismic wave propagation in media with discrete distributions of fractures using the pseudospectral method. The implementation of fractures with a vanishing width in the 2-D finite-difference grids is done using an effective medium theory (that is, the Coates and Schoenberg method). Fractures are treated as highly compliant interfaces inside a solid rock mass. For the physical representation of the fractures the concept of linear slip deformation or the displacement discontinuity method is used. According to this model, the effective compliance of a rock mass with one or several fracture sets can be found as the sum of the compliances of the host (background) rock and those of all the fractures. To first order, the background rock and fracture parameters can be related to the effective anisotropic coefficients, which govern the influence of anisotropy on various seismic signatures. We test the validity of the method and examine the accuracy of the synthetic seismograms by a comparison with theoretical ray traveltimes. We present three numerical examples to show the effects of different fracture distributions. The first example shows that different spatial distributions of the same fractures produce different wavefield characteristics. The second example examines the effects of variation of fracture scale length (size) compared with the wavelength. The final example examines the case of fractures with a power-law (fractal) distribution of sizes and shows how that affects the wavefield propagation in fractured rock. We conclude that characterization of fractured rock based on the concept of seismic anisotropy using effective medium theories must be used with caution. Scale length and the spatial distributions of fractures, which are not properly treated in such theories, have a strong influence on the characteristics of wave propagation.

Evidence for Neoproterozoic orogenesis and early high temperature Scandian deformation events in the southern East Greenland Caledonides

Abstract: Integrated field structural studies and SHRIMP U–Pb zircon and monazite dating have been undertaken in Renland, west of Scoresby Sund district in the southern part of the East Greenland Caledonides. Southwest Renland is dominated by metasedimentary rocks correlated with the Krummedal supracrustal succession of East Greenland and which on Renland were intruded by augen granites. Krummedal psammite from Renland yielded a spectrum of Mesoproterozoic to Palaeoproterozoic detrital U–Pb zircon dates, the youngest of which indicate deposition of the psammite occurred c. 1000 Ma ago, thus post-dating Grenvillian continent–continent collision in North American Laurentia. These Krummedal metasediments were deformed into regional nappe-scale folds prior to metamorphism, crustal anatexis and genesis of augen granites; an example of the latter has been dated at 915±18 Ma (U–Pb zircon). This demonstrates early Neoproterozoic high-temperature tectono-metamorphism affecting rocks within the southern East Greenland Caledonides, broadly contemporaneous with similar rocks farther north in East Greenland and with Sveconorwegian events on Baltica. Still in southwestern Renland, a later thermal event led to development of uppermost amphibolite to granulite facies metamorphic assemblages, veins and patches of in situ garnetiferous melt-bearing neosome in both metasediments (432±6 Ma, U–Pb zircon) and in the augen granites, and contemporaneous biotite-bearing granite sheets in top-down-to-the-E extensional shear zones (434±5 Ma, U–Pb zircon). Monazites from southwestern Renland record Caledonian thermal events as late as 410−400 Ma. In contrast, southeastern Renland is dominated by quartzofeldspathic migmatites with a strongly Caledonian signature but enclosing relicts of augen granite and retrogressed granulite facies psammitic and pelitic metasediment. There is also a sequence of Caledonian granitoid intrusions. Two samples from a hypersthene monzonite intrusion yielded U–Pb zircon dates of 424±8 Ma and 424±6 Ma. This pluton shows the marginal effects of the regional migmatization and was intruded early in the sequence of granitoid emplacement. An amphibolite facies migmatite, textural evidence from which suggests that it had never hosted granulite facies assemblages, records zircon growth at 423±6 Ma, and closure of monazite by 402±10 Ma. High grade metamorphism, and the protracted sequence of granitoid emplacement and still younger thermal events which together span the period between 430 and 400 Ma may, in part, reflect complicated lithospheric dynamics associated with subduction outboard of the Laurentian margin. Crustal segments carrying the relict evidence of Neoproterozoic and early Caledonian events must then quickly have been thrust northwestwards in foreland-propagating, northwesterly directed thrusts over Cambro-Ordovician platformal sequences on the Laurentian margin. This records the final closure of Iapetus, encroachment of Baltica and continent–continent collision from late Llandovery times (425–430 Ma).