The Meso-Cenozoic thermo-tectonic evolution of the Eastern Pyrenees: an <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar fission track and (U–Th)/He thermochronological study of the Canigou and Mont-Louis massifs

Vertical displacements on the SW–NE Têt fault (Eastern Pyrenees Axial Zone, France), which separates the Variscan Canigou-Carança and Mont-Louis massifs, were constrained using a thermochronologic multi-method approach. ⁴⁰Ar/³⁹Ar data from the granitic Mont-Louis massif record its Variscan cooling history and reveal no ages younger than Early Cretaceous, while the Canigou-Carança gneiss massif records systematically younger ⁴⁰Ar/³⁹Ar ages. These younger ⁴⁰Ar/³⁹Ar ages in the Canigou-Carança gneiss massif are the result of partial to total rejuvenation of argon isotopic systems related to a thermal flow coeval with the Cretaceous HT-BP metamorphism in the North Pyrenean Zone. Only the deepest rocks from the Canigou-Carança suffered this extensive Mid-Cretaceous thermal overprint probably due to differential burial around 4 km at that time. The post Mid-Cretaceous vertical displacements along the Têt fault are recorded by “low” temperature thermochronology using K-feldspar ⁴⁰Ar/³⁹Ar, zircon and apatite fission track and (U–Th)/He datings. The Mont-Louis granite samples experienced a long period of protracted cooling reflecting a lack of thermo-tectonic activity in this area from Late Palaeozoic to Early Cenozoic, followed by cooling from 55–60 Ma to Late Eocene at a mean rate of 15–20°C/Ma in the final stage. This cooling stage corresponds to Têt fault reactivation with a reversed component, promoting exhumation of the Mont-Louis roof zone contemporaneously with the south-vergent Pyrenean thrusting. In the Canigou-Carança massif, the main cooling event occurred from 32 to 18 Ma at a maximum rate of 30°C/Ma during Early Oligocene followed by a more moderate rate of 3°C/Ma from Late Oligocene to Early Burdigalian, coeval with the normal reactivation of the Têt fault in brittle conditions that accommodated the final exhumation of the massif during the opening of the Gulf of Lion.     

Provenance of late Palaeozoic metasediments of the Patagonian proto-Pacific margin (southernmost Chile and Argentina)

In this provenance study of late Palaeozoic metasediments of the Eastern Andean Metamorphic Complex (EAMC) along the south Patagonian proto-Pacific margin of Gondwana, the palaeogeological setting of the continental margin in Devonian–Carboniferous and Permian times is reconstructed. The study is based on detrital heavy mineral contents, chemical compositions of tourmaline grains, and whole rock element and Nd-Sr isotopic compositions. Element and isotopic compositions reveal that Devonian–Carboniferous metaturbidites deposited before the development of a Late Carboniferous–Permian magmatic arc along the margin were mainly fed from felsic, recycled, old continental rocks. The last recycling phase involved erosion of metasediments that were exposed in Patagonia. Feeder systems to the basin cut either through epidote-rich or garnet-rich metasediments. In Permian time, EAMC metaturbidites were deposited next to the evolving magmatic arc and were derived from felsic, crustal rocks. Two provenance domains are recognised. The metasediments of the northern one are chemically similar to those of the Devonian–Carboniferous metasediments. This domain was fed from the metasedimentary host rocks of the magmatic arc. The southern domain probably was fed from the arc proper, as indicated mainly by the dominance of metaplutonic lithic fragments, abundant detrital biotite, and the major element composition of the metasediments.     

High heat flow anomalies on an old oceanic plate observed seaward of the Japan Trench

Thirty-three new measurements on the seaward slope and outer rise of the Japan Trench along a parallel of 38°45′N revealed the existence of high heat flow anomalies on the subducting Pacific plate, where the seafloor age is about 135 m.y.. The most prominent anomaly with the highest value of 114 mW/m² is associated with a small mound on the outer rise, which was reported to be a kind of mud volcano. On the seaward slope of the trench, heat flow is variable: high (70–90 mW/m²) at some locations and normal for the seafloor age (about 50 mW/m²) at others. The spatial variation of heat flow may be related to development of normal faults and horst/graben structures due to bending of the Pacific plate before subduction, with fluid flow along the fault zones enhancing the vertical heat transfer. Possible heat sources of the high heat flow anomalies are intra-plate volcanism in the last several million years like that discovered recently on the Pacific plate east of the Japan Trench.     

Estimates of heat flow and heat production and a thermal model of the São Francisco craton

An updated analysis of geothermal data from the highland area of eastern Brazil has been carried out and the characteristics of regional variations in geothermal gradients and heat flow examined. The database employed includes results of geothermal measurements at 45 localities. The results indicate that the Salvador craton and the adjacent metamorphic fold belts northeastern parts of the study area are characterized by geothermal gradients in the range of 6–17°C/km. The estimated heat flow values fall in the range of 28–53 mW/m², with low values in the cratonic area relative to the fold belts. On the other hand, the São Francisco craton and the intracratonic São Francisco sedimentary basin in the southwestern parts are characterized by relatively higher gradient values, in the range of 14–42°C/km, with the corresponding heat flow values falling in the range of 36–89 mW/m². Maps of regional variations indicate that high heat flow anomaly in the São Francisco craton is limited to areas of sedimentary cover, to the west of the Espinhaço mountain belt. Crustal thermal models have been developed to examine the implications of the observed intracratonic variations in heat flow. The thermal models take into consideration variation of thermal conductivity with temperature as well as change of radiogenic heat generation with depth. Vertical distributions of seismic velocities were used in obtaining estimates of radiogenic heat production in crustal layers. Crustal temperatures are calculated based on a procedure that makes simultaneous use of the Kirchoff and Generalized Integral Transforms, providing thereby analytical solutions in 2D and 3D geometry. The results point to temperature variations of up to 300°C at the Moho depth, between the northern Salvador and southern São Francisco cratons. There are indications that differences in rheological properties, related to thermal field, are responsible for the contrasting styles of deformation patterns in the adjacent metamorphic fold belts.     

Geochemistry of amphibolitized eclogites and cross-cutting tonalitic–trondhjemitic dykes in the Metamorphic Kimi Complex in East Rhodope (N.E. Greece): implications for partial melting at the base of a thickened crust

In the ultra-high pressure Metamorphic Kimi Complex widespread tonalitic–trondhjemitic dykes, with an intrusion age ca. 65–63 Ma, cross-cut boudins and layers of amphibolitized eclogites. Geochemical investigation proclaims the tied genetic relationship of the amphibolitized eclogites and the associated tonalitic–trondhjemitic dykes. The major and trace element contents and rare earth element patterns of the amphibolitized eclogites indicate formation of their protoliths by fractional crystallization of tholeiitic magmas in a back-arc environment. The tonalites and trondhjemites are characterized by moderate to high Sr contents (>130 ppm), and low Y (<8.2 ppm) and heavy rare earth element contents (Yb content of 0.19–0.88 ppm). The chemical composition of the tonalitic and trondhjemitic dykes are best explained by partial melting of a tholeiitic source like the amphibolitized eclogites with residual garnet and amphibole, at the base of a thickened crust during Early Tertiary subduction/accretion at the southern margins of the European continent.     

Evaluation of change in subsurface thermal environment due to groundwater flow in the Tokyo Lowland,Japan

Information on the distribution of subsurface temperature and hydraulic heads at 24 observation wells in and around the Tokyo Lowland, the eastern part of the Tokyo Metropolitan area, were examined to make clear the relationship between groundwater and the subsurface thermal environment in the urban area. Minimums in temperature–depth logs due to subsurface temperature increasing at shallow parts were recognized in 21 wells. This fact shows subsurface temperature is affected by ground surface warming in almost all of this area. Deeper than minimums, where the effects of surface warming became relatively small, regional variation is observed as follows: high temperatures are shown in the central part to the southern part, and low temperatures shown in the inland to eastern part. The high temperature area corresponds to an area where the lower boundary of groundwater flow is relatively shallow. This area corresponds also to an area with severe land subsidence resulting from excessive groundwater pumping. It is considered that this high temperature area is formed by the effects of upward groundwater flow affected by hydrogeological conditions and pumping. On the other hand, a comparison between past data (1956–1967) and present data (2001–2003) revealed widespread decreasing temperature in the inland area. This is explained by downward groundwater flow based on an analysis of temperature–depth logs. This fact suggests that subsurface temperature is not only increasing from the effects of surface warming but also decreasing from the effects of groundwater environment change due to pumping.     

Hydrogeologic assessment of escalating groundwater exploitation in the Indus Basin, Pakistan

Groundwater development has contributed significantly to food security and reduction in poverty in Pakistan. Due to rapid population growth there has been a dramatic increase in the intensity of groundwater exploitation leading to declining water tables and deteriorating groundwater quality. In such prevailing conditions, the hydrogeological appraisal of escalating groundwater exploitation has become of paramount importance. Keeping this in view, a surface water–groundwater quantity and quality model was developed to assess future groundwater trends in the Rechna Doab (RD), a sub-catchment of the Indus River Basin. Scenario analysis shows that if dry conditions persist, there will be an overall decline in groundwater levels of around 10 m for the whole of RD during the next 25 years. The lower parts of RD with limited surface water supplies will undergo the highest decline in groundwater levels (10 to 20 m), which will make groundwater pumping very expensive for farmers. There is a high risk of groundwater salinization due to vertical upconing and lateral movement of highly saline groundwater into the fresh shallow aquifers in the upper parts of RD. If groundwater pumping is allowed to increase at the current rate, there will be an overall decline in groundwater salinity for the lower and middle parts of RD because of enhanced river leakage.     

Strain localization due to structural in-homogeneities in the Central European Basin System

The large-scale crustal deformations observed in the Central European Basin System (CEBS) are the result of the interplay between several controlling factors, among which lateral rheological heterogeneities play a key role. We present a finite-element integral thin sheet model of stress and strain distribution within the CEBS. Unlike many previous models, this study is based on thermo-mechanical data to quantify the impact of lateral contrasts on the tectonic deformation. Elasto-plastic material behaviour is used for both the mantle and the crust, and the effects of the sedimentary fill are also investigated. The consistency of model results is ensured through comparisons with observed data. The results resemble the present-day dynamics and kinematics when: (1) a weak granite-like lower crust below the Elbe Fault System is modelled in contrast to a stronger lower crust in the area extending north of the Elbe Line throughout the Baltic region; and (2) a transition domain in the upper mantle is considered between the shallow mantle of the Variscan domain and the deep mantle beneath the East European Craton (EEC), extending from the Elbe Line in the south till the Tornquist Zone. The strain localizations observed along these structural contrasts strongly enhance the dominant role played by large structural domains in stiffening the propagation of tectonic deformation and in controlling the basin formation and the evolution in the CEBS.