Na depletion in modern adakites via melt / rock reaction within the sub-arc mantle

Interaction between slab-derived melt and mantle peridotite and the role of slab melt as a metasomatizing agent in the sub-arc mantle is being increasingly recognized. Adakite, the slab melt erupted on the surface, usually exhibits anomalously high MgO, CaO, Cr and Ni contents that indicate interaction with mantle peridotitite. Here we note that Cenozoic adakites have Na₂O contents below 5.8 wt.% with ∼95% samples lower than 5.0 wt.%, and are generally depleted in this component relative to experimental basalt partial melts (mostly beyond 5.0 wt.% and up to 9.0 wt.% Na₂O) produced under 1.5-3.0 GPa conditions that are most relevant to adakite production. We interpret the adakite Na depletion to be also a consequence of the melt / rock reaction that takes place within the hot mantle wedge. During ascent and reaction with mantle peridotite, primary adakite melts gain mantle components MgO, CaO, Cr and Ni but lose Na₂O, SiO₂ and perhaps K₂O to the mantle, leading to Na-rich mantle metasomatism. Selective assimilation of predominately mantle clinopyroxene, some spinel and minor olivine at high T/P has been considered to be an important process in producing high-Mg adakites from primary low-Mg slab melts [Killian, R., Stern, C. R., 2002. Constraints on the interaction between slab melts and the mantle wedge from adakitic glass in peridotite xenoliths. Eur. J. Mineral. 14, 25-36]. In such a process, Na depletion in the assimilated melt is the result of dilution due to the increase in melt mass. Phase relationships in the reaction system siliceous melt + peridotite and quantitative calculation suggest that assimilation of mantle clinopyroxene, olivine and spinel and fractional crystallization of sodic amphibole and orthopyroxene, under conditions of moderate T/P and increasing melt mass, is also an important process that modifies the composition of adakites and causes the Na depletion.     

Mantle plumes from top to bottom

Hotspots include midplate features like Hawaii and on-axis features like Iceland. Mantle plumes are a well-posed hypothesis for their formation. Starting plume heads provide an explanation of brief episodes of flood basalts, mafic intrusions, and radial dike swarms. Yet the essence of the hypothesis hides deep in the mantle. Tests independent of surface geology and geochemistry to date have been at best tantalizing. It is productive to bare the current ignorance, rather than to dump the plume hypothesis. One finds potentially fruitful lines of inquiry using simple dynamics and observations. Ancient lithospheric xenoliths may reveal heating by plumes and subsequent thermal equilibration in the past. The effect at the base of the chemical layer is modest 50-100 K for transient heating by plume heads. Thinning of nonbuoyant platform lithosphere is readily observed but not directly attributable to plumes. The plume history in Antarctica is ill constrained because of poor geological exposure. This locality provides a worst case on what is known about surface evidence of hotspots. Direct detection of plume tail conduits in the mid-mantle is now at the edge of seismic resolution. Seismology does not provide adequate resolution of the deep mantle. We do not know the extent of a chemically dense dregs layer or whether superplume regions are cooler or hotter than an adiabat in equilibrium with the asthenosphere. Overall, mid-mantle seismology is most likely to give definitive results as plume conduits are the guts of the dynamic hypothesis. Finding them would bring unresolved deep and shallow processes into place.     

Fracture analysis in the south-western Corinth rift (Greece) and implications on fault hydraulic behavior

This paper reviews the data concerning the fracture network and the hydraulic characteristics of faults in an active zone of the Gulf of Corinth. Pressure gap measured through fault planes shows that in this area the active normal faults (Aigion, Helike) act, at least temporarily and locally, as transversal seal. The analysis of the carbonate cements in the fractures on both the hangingwall and the footwall of the faults also suggests that they have acted as local seals during the whole fault zone evolution. However, the pressure and the characteristics of the water samples measured in the wells indicate that meteoric water circulates from the highest part of the relief to the coast, which means it goes through the fault zones. Field quantitative analysis and core studies from the AIG-10 well have been performed to define both regional and fault-related fracture networks. Then laboratory thin section observations have been done to recognize the different fault rocks characterizing the fault zone components. These two kinds of approach give information on the permeability characteristics of the fault zone. To synthesize the data, a schematic conceptual 3D fluid flow modeling has been performed taking into account fault zone permeability architecture, sedimentation, fluid flow, fault vertical offset and meteoric water influx, as well as compaction water flow. This modeling allows us to fit all the data with a model where the fault segments act as a seal whereas the relays between these segments allow for the regional flow from the Peloponnese topographic highs to the coast.     

Statistical evaluation of methods for the calculation of static formation temperatures in geothermal and oil wells using an extension of the error propagation theory

The feasibility of using the more sophisticated weighted least-squares (WLS) model, as opposed to the traditional ordinary least-squares (OLS), in linear regressions of BHT data to estimate the static formation temperatures (SFT) was investigated. The most commonly used analytical methods (line-source; spherical and radial heat flow; and cylindrical heat source) were applied. Error propagation equations were derived to calculate errors in the time function of each method. These errors were combined with the BHT measurement errors to compute weighting factors for applying the WLS. Intercept uncertainties were estimated for all regressions using sets of synthetic and actual borehole logs taken from geothermal and oil applications. SFT computed with the spherical and radial heat flow method were generally greater than those from the other two methods.     

Temperature–composition–time (T–X–t) data from authigenic K-feldspar: An integrated methodology for dating fluid flow events

Integration of fluid inclusion analysis with high spatial resolution Ar–Ar dating of K-feldspar cements has been used to resolve and reconstruct palaeo-fluid flow. Fluid inclusion analysis allows discrimination of distinct cement phases, thereby identifying discrete episodes of fluid flow. Ar–Ar dating of the same cements via high spatial resolution laserprobe establishes absolute age constraints on the framework previously constructed. Integration of these two datasets yields temperature–composition–time data.     

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.     

Upper mantle structure of the Northern Eurasia from peaceful nuclear explosion data

Several long-range seismic profiles were carried out in Russia with Peaceful Nuclear Explosions (PNE). The data from 25 PNEs recorded along these profiles were used to compile a 3-D upper mantle velocity model for the central part of the Northern Eurasia. 2-D crust and upper mantle models were also constructed for all profiles using a common methodology for wavefield interpretation. Five basic boundaries were traced over the study area: N1 boundary (velocity level, V = 8.35 km/s; depth interval, D = 60–130 km), N2 (V = 8.4 km/s; D = 100–140 km), L (V = 8.5 km/s; D = 180–240 km) and H (V = 8.6 km/s; D = 300–330 km) and structural maps were compiled for each boundary. Together these boundaries describe a 3-D upper mantle model for northern Eurasia. A map characterised the velocity distribution in the uppermost mantle down to a depth of 60 km is also presented. Mostly horizontal inhomogeneity is observed in the uppermost mantle, and the velocities range from the average 8.0–8.1 km/s to 8.3–8.4 km/s in some blocks of the Siberian Craton. At a depth of 100–200 km, the local high velocity blocks disappear and only three large anomalies are observed: lower velocities in West Siberia and higher velocities in the East-European platform and in the central part of the Siberian Craton. In contrast, the depths to the H boundary are greater beneath the craton and lower beneath in the West Siberian Platform. A correlation between tectonics, geophysical fields and crustal structure is observed. In general, the old and cold cratons have higher velocities in the mantle than the young platforms with higher heat flows.Structural peculiarities of the upper mantle are difficult to describe in form of classical lithosphere–asthenosphere system. The asthenosphere cannot be traced from the seismic data; in contrary the lithosphere is suggested to be rheologically stratified. All the lithospheric boundaries are not simple discontinuities, they are heterogeneous (thin layering) zones which generate multiphase reflections. Many of them may be a result of fluids concentrated at some critical P–T conditions which produce rheologically weak zones. The most visible rheological variations are observed at depths of around 100 and 250 km.     

Seismic lamination and anisotropy of the Lower Continental Crust

Seismic lamination in the lower crust associated with marked anisotropy has been observed at various locations. Three of these locations were investigated by specially designed experiments in the near vertical and in the wide-angle range, that is the Urach and the Black Forrest area, both belonging to the Moldanubian, a collapsed Variscan terrane in southern Germany, and in the Donbas Basin, a rift inside the East European (Ukrainian) craton. In these three cases, a firm relationship between lower crust seismic lamination and anisotropy is found. There are more cases of lower-crustal lamination and anisotropy, e.g. from the Basin and Range province (western US) and from central Tibet, not revealed by seismic wide-angle measurements, but by teleseismic receiver function studies with a P–S conversion at the Moho. Other cases of lamination and anisotropy are from exhumed lower crustal rocks in Calabria (southern Italy), and Val Sesia and Val Strona (Ivrea area, Northern Italy). We demonstrate that rocks in the lower continental crust, apart from differing in composition, differ from the upper mantle both in terms of seismic lamination (observed in the near-vertical range) and in the type of anisotropy. Compared to upper mantle rocks exhibiting mainly orthorhombic symmetry, the symmetry of the rocks constituting the lower crust is either axial or orthorhombic and basically a result of preferred crystallographic orientation of major minerals (biotite, muscovite, hornblende). We argue that the generation of seismic lamination and anisotropy in the lower crust is a consequence of the same tectonic process, that is, ductile deformation in a warm and low-viscosity lower crust. This process takes place preferably in areas of extension. Heterogeneous rock units are formed that are generally felsic in composition, but that contain intercalations of mafic intrusions. The latter have acted as heat sources and provide the necessary seismic impedance contrasts. The observed seismic anisotropy is attributed to lattice preferred orientation (LPO) of major minerals, in particular of mica and hornblende, but also of olivine. A transversely isotropic symmetry system, such as expected for sub-horizontal layering, is found in only half of the field studies. Azimuthal anisotropy is encountered in the rest of the cases. This indicates differences in the horizontal components of tectonic strain, which finally give rise to differences in the evolution of the rock fabric.     

An investigation of upper mantle heterogeneity beneath the Archaean and Proterozoic crust of western Canada from Lithoprobe controlled-source seismic experiments

Observations of upper mantle reflectivity at numerous locations around the world have been linked to the presence of a heterogeneous distribution of rock types within a broad layer of the upper mantle. This phenomenon is observed in wide-angle reflection data from Lithoprobe's Alberta Basement Transect [the SAREX and Deep Probe experiments of 1995] and Trans-Hudson Orogen Transect [the THoRE experiment of 1993]. SAREX and Deep Probe image the Archaean lithosphere of the Hearne and Wyoming Provinces, whereas THoRE images the Archaean and Proterozoic lithosphere of the Trans-Hudson Orogen and neighbouring areas.Finite-difference synthetic seismograms are used to constrain the position and physical properties of the reflective layer. SAREX/Deep Probe modelling uses a 2-D visco-elastic finite-difference routine; THoRE modelling uses a pseudospectral algorithm. In both cases, the upper mantle is parameterized in terms of two media. One medium is the background matrix; the other is statistically distributed within the first as a series of elliptical bodies. Such a scheme is suitable for modelling: (1) variations in lithology (e.g., a peridotite matrix with eclogite lenses) or (2) variations in rheology (e.g., lenses of increased strain within a less strained background).The synthetic seismograms show that the properties of heterogeneities in the upper mantle do not change significantly between the two Lithoprobe transects. Beneath the Trans-Hudson Orogen in Saskatchewan, the layer is best modelled to lie at depths between 80 and 150 km. Based on observations from perpendicular profiles, anisotropy of the heterogeneities is inferred. Beneath the Precambrian domains of Alberta, 400 km to the west, upper mantle heterogeneities are modelled to occur between depths of 90 and 140 km. In both cases the heterogeneous bodies within the model have cross-sectional lengths of tens of kilometers, vertical thicknesses less than 1 km, and velocity contrasts from the background of − 0.3 to − 0.4 km/s. Based on consistency with complementary data and other results, the heterogeneous layer is inferred to be part of the continental lithosphere and may have formed through lateral flow or deformation within the upper mantle.     

Modèles géostatistiques de concentrations ou de débits le long des cours d'eau

Estimating concentrations or flow rates along a stream network requires specific models. Two classes of models, recently proposed in the literature, are generalized, to the intrinsic case in particular. We present a global construction by ‘streams’, i.e. on the whole set of paths between sources and outlet. Combining stationary or intrinsic one-dimensional random functions leads to stationary or intrinsic models on segments, with discontinuities at the forks. A construction from outlet to sources, leads to stationary or intrinsic models on each stream, without any discontinuity at the forks. The linear variogram is found as a particular case. The extension to the linear model of coregionalization is immediate, allowing a multivariate modelling of concentrations. To cite this article: C. de Fouquet, C. Bernard-Michel, C. R. Geoscience 338 (2006).