Anatolian surface wave evaluated at GEOFON Station ISP Isparta, Turkey

We have studied seismic surface waves of 255 shallow regional earthquakes recently recorded at GEOFON station ISP (Isparta, Turkey) and have selected these 52 recordings with high signal-to-noise ratio for further analysis. An attempt was made by the simultaneous use of the Rayleigh and Love surface wave data to interpret the planar crust and uppermost mantle velocity structure beneath the Anatolian plate using a differential least-square inversion technique. The shear-wave velocities near the surface show a gradational change from approximately 2.2 to 3.6 km s^− 1 in the depth range 0–10 km. The mid-crustal depth range indicating a weakly developed low velocity zone has shear-wave velocities around 3.55 km s^− 1. The Moho discontinuity characterizing the crust–mantle velocity transition appears somewhat gradual between the depth range  25–45 km. The surface waves approaching from the northern Anatolia are estimated to travel a crustal thickness of  33 km whilst those from the southwestern Anatolia and part of east Mediterranean Sea indicate a thicker crust at  37 km. The eastern Anatolia events traveled even thicker crust at  41 km. A low sub-Moho velocity is estimated at  4.27 km s^− 1, although consistent with other similar studies in the region. The current velocities are considerably slower than indicated by the Preliminary Reference Earth Model (PREM) in almost all depth ranges.     

Formation, subduction, and exhumation of Penninic oceanic crust in the Eastern Alps: time constraints from Ar/Ar geochronology

New ⁴⁰Ar/³⁹Ar geochronology places time constraints on several stages of the evolution of the Penninic realm in the Eastern Alps. A 186±2 Ma age for seafloor hydrothermal metamorphic biotite from the Reckner Ophiolite Complex of the Pennine–Austroalpine transition suggests that Penninic ocean spreading occurred in the Eastern Alps as early as the Toarcian (late Early Jurassic). A 57±3 Ma amphibole from the Penninic subduction–accretion Rechnitz Complex dates high-pressure metamorphism and records a snapshot in the evolution of the Penninic accretionary wedge. High-pressure amphibole, phengite, and phengite+paragonite mixtures from the Penninic Eclogite Zone of the Tauern Window document exhumation through ≤15 kbar and >500 °C at 42 Ma to 10 kbar and 400 °C at 39 Ma. The Tauern Eclogite Zone pressure–temperature path shows isothermal decompression at mantle depths and rapid cooling in the crust, suggesting rapid exhumation. Assuming exhumation rates slower or equal to high-pressure–ultrahigh-pressure terrains in the Western Alps, Tauern Eclogite Zone peak pressures were reached not long before our high-pressure amphibole age, probably at ≤45 Ma, in accordance with dates from the Western Alps. A late-stage thermal overprint, common to the entire Penninic thrust system, occurred within the Tauern Eclogite Zone rocks at 35 Ma. The high-pressure peak and switch from burial to exhumation of the Tauern Eclogite Zone is likely to date slab breakoff in the Alpine orogen. This is in contrast to the long-lasting and foreland-propagating Franciscan-style subduction–accretion processes that are recorded in the Rechnitz Complex.     

TRANSALP—deep crustal Vibroseis and explosive seismic profiling in the Eastern Alps

The TRANSALP consortium, comprising institutions from Italy, Austria and Germany, carried out deep seismic reflection measurements in the Eastern Alps between Munich and Venice in 1998, 1999 and 2001. In order to complement each other in resolution and depth range, the Vibroseis technique was combined with simultaneous explosive source measurements. Additionally, passive cross-line recording provided three-dimensional control and alternative north–south sections. Profits were obtained by the combination of the three methods in sectors or depths where one method alone was less successful.The TRANSALP sections clearly image a thin-skinned wedge of tectonic nappes at the northern Alpine front zone, unexpected graben or half-graben structures within the European basement, and, thick-skinned back-thrusting in the southern frontal zone beneath the Dolomite Mountains. A bi-vergent structure at crustal scale is directed from the Alpine axis to the external parts. The Tauern Window obviously forms the hanging wall ramp anticline above a southward dipping, deep reaching reflection pattern interpreted as a tectonic ramp along which the Penninic units of the Tauern Window have been up-thrusted.The upper crystalline crust appears generally transparent. The lower crust in the European domain is characterized by a 6–7 km thick laminated structure. On the Adriatic side the lower crust displays a much thicker or twofold reflective pattern. The crustal root at about 55 km depth is shifted around 50 km to the south with respect to the main Alpine crest.     

Surface kinematics in the Alpine–Carpathian–Dinaric and Balkan region inferred from a new multi-network GPS combination solution

The understanding of the intraplate tectonics of Central Europe requires a detailed picture of how stress is transferred from the interaction of the Eurasian, Nubian and Anatolian plates to the Alpine, Carpathian, Pannonian and Dinaric regions. Recent strain distribution is controlled by the Adria horizontal push, by the Vrancea vertical slab pull and associated horizontal displacements, and by the Aegean/Anatolia extension and slab-roll back. We present a horizontal velocity field for the Alpine-Carpathian-Pannonic-Dinaric and Balkan regions resulting from a new combination of seven different GPS networks formed from permanent and campaign stations. Dedicated velocity profiles in two specific regions are studied in detail. One is the Alpine Pannonian region, with a detailed picture of the NS indentation of the Adria microplate into the Southern Alps, in NE Italy, the deformation in the Tauern Window and the eastwards kinematics of a Pannonian plate fragment. The second study region includes Transylvania, the Southern Carpathians up to the Aegean sea and Albania, where a major right lateral shear deformation exists as a consequence of the NE convergence of the Apulia platform towards the Dinarids, and the SSW motion of Macedonia, Western Bulgaria and Rumania, related to the Hellenic arc dynamics in the Eastern Mediterranean. The profiles in the Alpine–Pannonian area indicate that a velocity drop of 2.5 +/− 0.4 mm/yr associated with the Adria indentation concentrates on a segment of some 50 km south of the Periadriatic fault. The deformation becomes extensional by a similar amount just north of the Periadriatic fault, in the Tauern Window, where the updoming of the Tauern Window implies vertical motion which could well be associated with surface extension. In the EW profile, we observe a sudden velocity change of 1.5 +/− 0.2 mm/yr in 20 km, corresponding to the right lateral Lavant fault, which seems to mark the border between dominant indentation kinematics to the West and dominant extrusion kinematics to the East.Three profiles are considered in Southern and Eastern Europe: one across the lower Adriatic sea from Apulia in Italy to the southern Dinarides, which enables it to constrain the velocity drop associated with the subduction of the Adria microplate into the Dinarides to 3.2 +/− 0.5 mm/yr in 140 km. The second profile is longitudinal and constrains the velocity inversion of 7.4 +/− 1.0 mm/yr in 350 km associated with right lateral shear faults in Albania. The third profile crosses the Transylvania region with a shortening of 2.3 +/− 1.0 mm/yr in 220 km, and the Wallachian–Moesian region up to the Chalcidic peninsula in N Greece. This lower part of the profile implies an extensional stretch of the upper crust of 3.2 +/− 0.9 mm/yr in 440 km, culminating in the Hellenic arc. Strain rate maps are presented in this regional scale, showing the excellent agreement between fault plane solutions of crustal earthquakes and the eigenvectors of the GPS derived two dimensional strain rate tensor.Three profiles are considered in the Balkan and SE Carpathians: one across the lower Adriatic sea from Apulia in Italy to the southern Dinarides, which enables to constrain the velocity drop associated to the subduction of the Adria microplate into the Dinarides to 3.2 +/− 0.5 mm/yr in 140 km. The second profile is longitudinal and constrains the velocity inversion of 7.4 +/− 1.0 mm/yr in 350 km associated to right lateral shear faults in Macedonia, a highly seismic region. The third profile crosses the Transylvania with a shortening2.3 +/− 1.0 mm/yr in 220 km, and the Wallachian–Moesian region up to the Chalcidic peninsula in N Greece. This lower part of the profile implies an extensional stretch of the upper crust of 3.2 +/− 0.9 mm/yr in 440 km, culminating in the Hellenic arc.     

Kinematics of the Inntal shear zone–sub-Tauern ramp fault system and the interpretation of the TRANSALP seismic section, Eastern Alps, Austria

A new interpretation of the Inntal–Tauern sector of the TRANSALP seismic section is presented. One of the most prominent contrasts in reflectivity in the TRANSALP seismic section is the contact between the Bajuvaric unit in the footwall and the overlying Tirolic unit and its basement across a moderately south-dipping interface. We trace this contact from the surface at the southern margin of the Inn valley to a depth of 5 km. There, the contact is deformed or cut by the Tauern Window northern margin. We define the contact between Bajuvaric and Tirolic units as Brixlegg thrust, which is older than Miocene Tauern window exhumation and has a Paleogene age. The sub-Tauern ramp connects with the Inntal fault system at the surface and roots below the Tauern window. Oblique thrust movements across this fault system in the Miocene caused exhumation of the hanging wall, where the fault has a ramp geometry, which is in the area of the TRANSALP cross section and west of it. East of the TRANSALP cross section, the fault system merges with Alpine basal thrust, which is a flat. No Miocene exhumation occurred above the flat.     

Crustal velocity inhomogeneities along the Hirapur–Mandla profile, central India and its tectonic implications

Wide-angle seismic and gravity data across the Narmada-Son lineament (NSL) in central India are analyzed to determine crustal structure, velocity inhomogeneities and hence constrain the tectonics of the lineament. We present the 2-D crustal velocity structure from deep wide-angle reflection data by using a ray-trace inverse approach. The main result of the study is the delineation of fault-bounded horst raised to a subsurface depth (1.5 km) and the Moho upwarp beneath the NSL. The crust below the basement consists of three layers with velocities of 6.45–6.7, 6.2–6.5 and 6.7–6.95 km/s and interface depths of about 5.5–8.7, 14–17 and 18–23 km along the profile. The low-velocity (6.2–6.5 km/s) layer goes up to a depth of 5 km and becomes the thickest part (13 km), while the overlying high-velocity (6.45–6.7 km/s) layer becomes the thinnest (3 km) and upper boundary lies at a depth of 1.5 km beneath the NSL. The overall uncertainties of various velocity and boundary nodes are of the order of ±0.12 km/s and ±1.40 km, respectively. The up-lifted crustal block and the up-warping Moho beneath the NSL indicate that the north and south faults bounding the NSL are deeply penetrated through which mafic materials from upper mantle have been intruded into the upper crust. Gravity modeling was also undertaken to assess the seismically derived crustal features and to fill the seismic data gap. The lateral and vertical heterogeneous nature of the structure and velocity inhomogeneities in the crust cause instability to the crustal blocks and played an important role in reactivation of the Narmada south fault during the 1997 Jabalpur earthquake.     

Influence of the precollisional stage on subduction dynamics and the buried crust thermal state: Insights from numerical simulations

At continental subduction initiation, the continental crust buoyancy may induce, first, a convergence slowdown, and second, a compressive stress increase that could lead to the forearc lithosphere rupture. Both processes could influence the slab surface P–T conditions, favoring on one side crust partial melting or on the opposite the formation of ultra-high pressure/low temperature (UHP-LT) mineral. We quantify these two effects by performing numerical simulations of subduction. Water transfers are computed as a function of slab dehydration/overlying mantle hydration reactions, and a strength decrease is imposed for hydrated mantle rocks. The model starts with an old oceanic plate ( 100 Ma) subducting for 145.5 Myr with a 5 cm/yr convergence rate. The arc lithosphere is thermally thinned between 100 km and 310 km away from the trench, due to small-scale convection occuring in the water-saturated mantle wedge. We test the influence of convergence slowdown by carrying on subduction with a decreased convergence rate (≤ 2 cm/yr). Surprisingly, the subduction slowdown yields not only a strong slab warming at great depth (> 80 km), but also a significant cooling of the forearc lithosphere at shallower depth. The convergence slowdown increases the subducted crust temperature at 90 km depth to 705 ± 62 °C, depending on the convergence rate reduction, and might thus favor the oceanic crust partial melting in presence of water. For subduction velocities ≤ 1 cm/yr, slab breakoff is triggered 20–32 Myr after slowdown onset, due to a drastic slab thermal weakening in the vicinity of the interplate plane base. At last, the rupture of the weakened forearc is simulated by imposing in the thinnest part of the overlying lithosphere a dipping weakness plane. For convergence with rates ≥ 1 cm/yr, the thinned forearc first shortens, then starts subducting along the slab surface. The forearc lithosphere subduction stops the slab surface warming by hot asthenosphere corner flow, and decreases in a first stage the slab surface temperature to 630 ± 20 °C at 80 km depth, in agreement with P–T range inferred from natural records of UHP-LT metamorphism. The subducted crust temperature is further reduced to 405 ± 10 °C for the crust directly buried below the subducting forearc. Such a cold thermal state at great depth has never been sampled in collision zones, suggesting that forearc subduction might not be always required to explain UHP-LT metamorphsim.     

Crustal structure at the easternmost termination of the Variscan belt based on CELEBRATION 2000 and ALP 2002 data

The eastern margin of the Variscan belt in Europe comprises plate boundaries between continental blocks and terranes formed during different tectonic events. The crustal structure of that complicated area was studied using the data of the international refraction experiments CELEBRATION 2000 and ALP 2002. The seismic data were acquired along SW–NE oriented refraction and wide-angle reflection profiles CEL10 and ALP04 starting in the Eastern Alps, passing through the Moravo-Silesian zone of the Bohemian Massif and the Fore-Sudetic Monocline, and terminating in the TESZ in Poland. The data were interpreted by seismic tomographic inversion and by 2-D trial-and-error forward modelling of the P waves. Velocity models determine different types of the crust–mantle transition, reflecting variable crustal thickness and delimiting contacts of tectonic units in depth. In the Alpine area, few km thick LVZ with the Vp of 5.1 km s^− 1 dipping to the SW and outcropping at the surface represents the Molasse and Helvetic Flysch sediments overthrust by the Northern Calcareous Alps with higher velocities. In the Bohemian Massif, lower velocities in the range of 5.0–5.6 km s^− 1 down to a depth of 5 km might represent the SE termination of the Elbe Fault Zone. The Fore-Sudetic Monocline and the TESZ are covered by sediments with the velocities in the range of 3.6–5.5 km s^− 1 to the maximum depth of 15 km beneath the Mid-Polish Trough. The Moho in the Eastern Alps is dipping to the SW reaching the depth of 43–45 km. The lower crust at the eastern margin of the Bohemian Massif is characterized by elevated velocities and high Vp gradient, which seems to be a characteristic feature of the Moravo-Silesian. Slightly different properties in the Moravian and Silesian units might be attributed to varying distances of the profile from the Moldanubian Thrust front as well as a different type of contact of the Brunia with the Moldanubian and its northern root sector. The Moho beneath the Fore-Sudetic Monocline is the most pronounced and is interpreted as the first-order discontinuity at a depth of 30 km.     

Late Miocene increasing exhumation rates in the eastern part of the Alps – implications from low temperature thermochronology

A new set of apatite fission‐track and apatite (U–Th)/He data reveals a hitherto undated late Miocene exhumation pulse in the eastern part of the Eastern Alps. While distinct parts of the study area, including the Seckauer Tauern, have been at near surface conditions (<100 °C) since the Eocene, the neighbouring Niedere Tauern experienced enhanced cooling and exhumation in the middle Miocene and again at the late Miocene/Pliocene boundary. Middle Miocene exhumation is interpreted as a result of tectonic escape and convergence that operated simultaneously during lateral extrusion of the Eastern Alps. As the higher late Miocene/Pliocene exhumation rates are restricted to a single tectonic block, namely the Niedere Tauern, we infer a tectonic trigger that is probably related to a change in the external stress field that affected the Alps during this time.     

Constraints from Moho geometry and crustal thickness on the geodynamic origin of the Vrancea Seismogenic Zone (Romania)

Reprocessing of industry deep seismic reflection data (Ramnicu Sarat and Braila profiles) from the SE Carpathian foreland of Romania provides important new constraints on geodynamic models for the origin of the intermediate depth Vrancea Seismogenic Zone (VSZ). Mantle (70–200 km) earthquakes of the VSZ are characterized by high magnitudes (greater than 6.5), frequent occurrence rates (approximately 25 years), and confinement in a very narrow (30 × 70 × 200 km³) near vertical zone atypical for a Wadati–Benioff plane, located in front of the orogen. These two deep (20 s) seismic reflection profiles (70 km length across the foreland) reveal (1) a high-amplitude, gently east-dipping reflection across most of the section from what we interpret to be the Moho at  15 s (40–42 km) on the Ramnicu Sarat line to  16 s (47–48 km) on the Braila line, (2) a thick sedimentary cover increasing in thickness from east (1 s;  800 m) to west (7.5 s; 14 km), (3) an eastward increase in crustal thickness from 38 km (near VSZ) to  45 km, (4) seismic and topographic evidence for a newly imaged, possibly seismically active basement fault with a surface offset of 30 m observed on the Ramnicu Sarat line, (5) a lack of notable west-dipping structures in the crust and across the Moho, and (6) variable displacements on Peceneaga–Camena Fault of  5 km at Moho and  200 m at the basement–sedimentary cover contact.These observations appear to argue against recent models for west-dipping subduction of oceanic lithosphere at or in the vicinity of the Vrancea Seismogenic Zone given the lack of west-dipping fabrics in the lower crust and across the crust–mantle boundary. Consequently, one possible explanation for the geodynamic origin of VSZ could be partial delamination of the continental lithosphere in an intra-plate setting along a sub-horizontal lithospheric interface in the Carpathian hinterland that likely involves remnant lithospheric coupling between the crust and uppermost mantle in the foreland.     

Shear wave velocity and attenuation structure for the shallow crust of the southern Korean peninsula from short period Rayleigh waves

We analyzed the short period Rayleigh waves from the first crustal-scale seismic refraction experiment in the Korean peninsula, KCRUST2002, to determine the shear wave velocity and attenuation structure of the uppermost 1 km of the crust in different tectonic zones of the Korean peninsula and to examine if this can be related to the surface geology of the study area. The experiment was conducted with two large explosive sources along a 300-km long profile in 2002. The seismic traces, recorded on 170 vertical-component, 2-Hz portable seismometers, show distinct Rayleigh waves in the period range between 0.2 s and 1.2 s, which are easily recognizable up to 30–60 km from the sources. The seismic profiles, which traverse three tectonic regions (Gyeonggi massif, Okcheon fold belt and Yeongnam massif), were divided into five subsections based on tectonic boundaries as well as lithology. Group and phase velocities for the five subsections obtained by a continuous wavelet transform method and a slant stack method, respectively, were inverted for the shear wave models. We obtained shear wave velocity models up to a depth of 1.0 km. Overall, the shear wave velocity of the Okcheon fold belt is lower than that of the Gyeonggi and Yeongnam massifs by  0.4 km/s in the shallowmost 0.2 km and by 0.2 km/s at depths below 0.2 km. Attenuation coefficients, determined from the decay of the fundamental mode Rayleigh waves, were used to obtain the shear wave attenuation structures for three subsections (one for each of the three different tectonic regions). We obtained an average value of Q_β^− 1 in the upper 0.5 km for each subsection. Q_β^− 1 for the Okcheon fold belt ( 0.026) is approximately three times larger than Q_β^− 1 for the massif areas ( 0.008). The low shear wave velocity in the Okcheon fold belt is consistent with the high attenuation in this region.     

Geothermobarometry and fluid inclusions of dioritic rocks in Bangladesh: Implications for emplacement depth and exhumation rate

Application of hornblende thermobarometry and fluid inclusion studies to the Palaeoproterozoic (1.7 Ga) basement rocks from Maddhapara, NW Bangladesh, provide information on the pressure and temperature (P–T) conditions of crystallization, the emplacement depth and composition of magmatic fluid. The basement rocks are predominantly diorite or quartz diorite with a mineral assemblage of plagioclase, hornblende, biotite, quartz, K-feldspar, titanite, and secondary epidote and chlorite. The calculated P–T conditions of the dioritic rocks are 680–725 °C and 4.9–6.4 kbar, which probably correspond to crystallization conditions. Fluid inclusion studies suggest that low- to medium-salinity (0–6.4 wt.% NaCl_eq) H₂O-rich fluids are trapped during the crystallization of quartz and plagioclase. The isochore range calculated for primary aqueous inclusions is consistent with the P–T condition obtained by geothermobarometry. The basement rocks likely crystallized at a depth of 17–22 km, with a minimum average exhumation rate of 12–15 m/Ma during Palaeoproterozoic to Lopingian time. Such slow exhumation indicates low relief continental shield surface during this period.     

Seismic image of the deep crust at the eastern margin of the Alps (Austria): indications for crustal extension in a convergent orogen

A 39-km-long deep seismic reflection profile recorded during two field campaigns in 1996 and 2002 provides a first detailed image of the deep crust at the eastern margin of the Eastern Alps (Austria). The ESE–WNW-trending, low-fold seismic line crosses Austroalpine basement units and extends approximately from 20 km west of the Penninic window group of Rechnitz to 60 km SSE of the Alpine thrust front.The explosive-source seismic data reveals a transparent shallow crust down to 5 km depth, a complexly reflective upper crust and a highly reflective lowermost crust. The upper crust is dominated by three prominent west-dipping packages of high-amplitude subparallel reflections. The upper two of these prominent packages commence at the eastern end of the profile at about 5 and 10 km depth and are interpreted as low-angle normal shear zones related to the Miocene exhumation of the Rechnitz metamorphic core complex. In the western portion of the upper crust, east-dipping and less significant reflections prevail. The lowermost package of these reflections is suggested to represent the overall top of the European crystalline basement.Along the western portion of the line, the lower crust is characterised by a 6–8-km-thick band of high-amplitude reflection lamellae, typically observed in extensional provinces. The Moho can be clearly defined at the base of this band, at approximately 32.5 km depth. Due to insufficient signal penetration, outstanding reflections are missing in the central and eastern portion of the lower crust. We speculate that the result of accompanying gravity measurements and lower crustal sporadic reflections can be interpreted as an indication for a shallower Moho in the east, preferable at about 30.5 km depth.The high reflectivity of the lowermost part of the lower crust and prominent reflection packages in the upper crust, the latter interpreted to represent broad extensional mylonite zones, emphasises the latest extensional processes in accordance with eastward extrusion.     

Seismic velocity structure of a large mafic intrusion in the crust of central Denmark from project ESTRID

The origin of regional sedimentary basins is being investigated by the ESTRID project (Explosion Seismic Transects around a Rift In Denmark). This project investigates the mechanisms of the formation of wide, regional basins and their interrelation to previous rifting processes in the Danish–Norwegian Basin in the North Sea region. In May 2004 a 143 km long refraction seismic profile was acquired along the strike direction of a suspected major mafic intrusion in the crust in central Denmark. The data confirms the presence of a body with high seismic velocity (> 6.5 km/s) extending from a depth of  10–12 km depth into the lower crust. There is a remarkable Moho relief between 27 and 34 km depth along this new along-strike profile as based on ray-tracing modelling of PmP reflections. The lack of PmP reflections at a zone of very high velocity in the lowest crust (7.3–7.5 km/s) suggests a possible location of a feeder channel to the batholith. The presence of volcanic rocks of Carboniferous–Permian age above the intrusion (mafic batholith) suggests a similar age of the intrusion. An older obliquely crossing profile and two new fan profiles deployed perpendicular to the main ESTRID profile, show that the batholith is about 30–40 km wide. The existence of this large mafic batholith supports the hypothesis that the origin of the Danish–Norwegian Basin is related to cooling and contraction after intrusion of large amounts of mafic melts into the crust during the late Carboniferous and early Permian. The data and interpretations from project ESTRID will form the basis for subsidence modelling. Tentatively, we interpret the formation of the Danish–Norwegian Basin as a thermal subsidence basin, which developed after widespread rifting of the region.     

A tentative 2D thermal model of central India across the Narmada-Son Lineament (NSL)

This work deals with 2D thermal modeling in order to delineate the crustal thermal structure of central India along two Deep Seismic Sounding (DSS) profiles, namely Khajuriakalan–Pulgaon and Ujjan–Mahan, traversing the Narmada-Son-Lineament (NSL) in an almost north–south direction. Knowledge of the crustal structure and P-wave velocity distribution up to the Moho, obtained from DSS studies, has been used for the development of the thermal model. Numerical results reveal that the Moho temperature in this region of central India varies between 500 and 580 °C. The estimated heat flow density value is found to vary between 46 and 49 mW/m². The Curie depth varies between 40 and 42 km and is in close agreement with the Curie depth (40±4 km) estimated from the analysis of MAGSAT data. Based on the present work and previous work, it is suggested that the major part of peninsular India consisting of the Wardha–Pranhita Godavari graben/basin, Bastar craton and the adjoining region of the Narmada Son Lineament between profiles I and III towards the north and northwest of the Bastar craton are characterized with a similar mantle heat flow density value equal to 23 mW/m². Variation in surface heat flow density values in these regions are caused by variation in the radioactive heat production and fluid circulation in the upper crustal layer.     

Folding of Europa's icy lithosphere: an analysis of viscous-plastic buckling and subsequent topographic relaxation

Potential contractional folds on Jupiter's icy moon Europa have been identified. The best example is at the extensional band Astypalaea Linea, where a series of subtle topographic undulations, 25 km in wavelength, possess parasitic tectonic features that support a folding origin. A scenario has been qualitatively proposed, whereby folds form via unstable contraction of the icy lithosphere, compensate for extension elsewhere on Europa, and then subsequently relax. Here, we quantitatively address this scenario, applying a model for viscous-plastic buckling of planetary lithospheres and finite element simulations of topographic relaxation. Our results suggest that the lithosphere of Europa could indeed be unstable, but the low required surface temperatures limit fold formation to higher latitudes, and the high required driving stresses (9–10 MPa) are difficult to achieve on the satellite. The depth to the brittle–ductile transition is well constrained, and high thermal gradients are indicated, implying heat flows near 100 mW m⁻². In addition, topographic relaxation progresses so slowly even at these heat flows that it is not a viable mechanism to eliminate such features over the age of Europa's surface. Given the paucity of identified folds, we conclude that the necessary conditions for their formation are rare and that lithospheric folding is a minor mechanism for compensating the large amounts of extension seen elsewhere on Europa.     

Heat flow, heat production and fission track data from the Hercynian basement around the Provençal Basin (Western Mediterranean)

In the complex structural framework of the Western Mediterranean. Hercynian areas are expected to be thermally preserved from the recent tectonic evolution. The thermal regime of these areas is studied using heat flow, heat production and fission track data. The surface heat flow is significantly higher in Corsica (76 ± 10 mW m⁻²) than in the Maures and Estérel (58 ± 2 mW m⁻²). Neither heat production nor erosion subsequent to the Alpine orogeny in Corsica can explain such a difference. It is suggested that a deep thermal source related to the asymmetric evolution of the Provençal basin could explain the higher heat flow in Corsica. A model of thermal structure based on the present day thermal regime of the Maures and Estérei is proposed for the stable Hercynian crust in this area. The mantle heat flow is 20–25 mW m⁻² and the temperature at Moho level is 375–500°C, depending on the thermal parameter distribution with depth.     

Thermal evolution of the Indian shield and its subjacent mantle

The estimate of the paleothermal state (at −3 Ga) of the Indian shield is deduced using pressure and temperature estimates from the mineral assemblages combined with models of the heat transport mechanism and depth distribution of radiogenic heat in the crust. We find that at 3 Ga ago, the mean thermal gradient within the crust, reduced heat flow and mean mantle temperature were 28°C/km 1.08 h.f.u. and 2350°C, respectively.     

Constraining models of the tectonic setting of the giant Olympic Dam iron oxide–copper–gold deposit, South Australia, using deep seismic reflection data

In the Gawler Craton, the completeness of cover concealing the crystalline basement in the region of the giant Olympic Dam Cu–Au deposit has impeded any sufficient understanding of the crustal architecture and tectonic setting of its IOCG mineral-system. To circumvent this problem, deep seismic reflection data were recently acquired from  250 line-km of two intersecting traverses, centered on the Olympic Dam deposit. The data were recorded to 18 s TWT ( 55 km). The crust consists of Neoproterozoic cover, in places more than 5 km thick, over crystalline basement with the Moho at depths of 13–14 s TWT ( 40–42 km). The Olympic Dam deposit lies on the boundary between two distinct pieces of crust, one interpreted as the Archean–Paleoproterozoic core to the craton, the other as a Meso–Neoproterozoic mobile belt. The host to the deposit, a member of the  1590 Ma Hiltaba Suite of granites, is situated above a zone of reduced impedance contrast in the lower crust, which we interpret to be source-region for its  1000 °C magma. The crystalline basement is dominated by thrusts. This contrasts with widely held models for the tectonic setting of Olympic Dam, which predict extension associated with heat from the mantle producing the high temperatures required to generate the Hiltaba Suite granites implicated in mineralization. We use the seismic data to test four hypotheses for this heat-source: mantle underplating, a mantle-plume, lithospheric extension, and radioactive heating in the lower crust. We reject the first three hypotheses. The data cannot be used to reject or confirm the fourth hypothesis.     

Paleomagnetic constraints on the paleogeography and oroclinal bending of the Devonian volcanic arc in Kazakhstan

Numerical modelling, incorporating coupling between surface processes and induced flow in the lower continental crust, is used to address the Quaternary evolution of the Gulf of Corinth region in central Greece. The post-Early Pleistocene marine depocentre beneath this Gulf overlies the northern margin of an older (Early Pleistocene and earlier) lacustrine basin, the Proto Gulf of Corinth Basin or PGCB. In the late Early Pleistocene, relief in this region was minimal but, subsequently, dramatic relief has developed, involving the creation of  900 m of bathymetry within the Gulf and the uplift by many hundreds of metres of the part of the PGCB, south of the modern Gulf, which forms the Gulf's main sediment supply. It is assumed that, as a result of climate change around 0.9 Ma, erosion of this sediment source region and re-deposition of this material within the Gulf began, both processes occurring at spatial average rates of  0.2 mm a^− 1. Modelling of the resulting isostatic response indicates that the local effective viscosity of the lower crust is  4 × 10¹⁹ Pa s, indicating a Moho temperature of  560 °C. It predicts that the  10 mm a^− 1 of extension across this  70 km wide model region, at an extensional strain rate of  0.15 Ma^− 1, is partitioned with  3 mm a^− 1 across the sediment source,  2 mm a^− 1 across the depocentre, and  5 mm a^− 1 across the ‘hinge zone’ in between, the latter value being an estimate of the extension rate on normal faults forming the major topographic escarpment at the southern margin of the Gulf. This modelling confirms the view, suggested previously, that coupling between this depocentre and sediment source by lower-crustal flow can explain the dramatic development in local relief since the late Early Pleistocene. The effective viscosity of the lower crust in this region is not particularly low; the strong coupling interpreted between the sediment source and depocentre results instead from their close proximity. In detail, the effective viscosity of the lower crust is expected to decrease northward across this model region, due to the northward increase in exposure of the base of the continental lithosphere to the asthenosphere; in the south the two are separated by the subducting Hellenic slab. The isostatic consequences of such a lateral variation in viscosity provide a natural explanation for why, since  0.9 Ma, the modern Gulf has developed asymmetrically over the northern part of the PGCB, leaving the rest of the PGCB to act as its sediment source.