The Mesozoic and Cenozoic motion of Adria (central Mediterranean): a review of constraints and limitations

This paper presents kinematic analysis on the motion of Adria, which is the continental mass that bridges Africa and Europe in the central Mediterranean. Palaeomagnetic data show a general coherence between the motion of Adria and Africa since the Late Paleozoic. This mutual motion, for the period from 120 Ma and the present, is verified by comparing inferred palaeolatitudes from relatively stable parts of Adria (Apulia, Gargano, Istria, and the Southern Alps) and the Hyblean Plateau, with latitudinal changes that are calculated from the motion of Africa with respect to hotspots. Additional constraints on the motion of Adria are provided from the Late Paleozoic-Early Mesozoic passive margin of Adria in the Ionian Sea. The seismic structure of the floor of the Ionian Sea resembles the structure of the oceanic crust in marginal back-arc basins, suggesting that it formed as a small ocean basin. Furthermore, the Ionian lithosphere in the Calabrian arc has been subjected to rapid rollback, which commonly occurs only when the subducting slab is made of oceanic lithosphere. This oceanic domain marks the Pennian-Triassic to Jurassic plate boundary between Adria and Africa, suggesting that a small amount of independent motion between Adria and Africa took place at that time. Since the Jurassic, Adria and Africa have shared a relatively coherent motion path. (C) 2004 Lavoisier SAS. All rights reserved.

Numerical modelling of chemical effects of magma solidification problems in porous rocks

The solidification of intruded magma in porous rocks can result in the following two consequences: (1) the heat release due to the solidification of the interface between the rock and intruded magma and (2) the mass release of the volatile fluids in the region where the intruded magma is solidified into the rock. Traditionally, the intruded magma solidification problem is treated as a moving interface (i.e. the solidification interface between the rock and intruded magma) problem to consider these consequences in conventional numerical methods. This paper presents an alternative new approach to simulate thermal and chemical consequences/effects of magma intrusion in geological systems, which are composed of porous rocks. In the proposed new approach and algorithm, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with the proposed mass source and physically equivalent heat source. The major advantage in using the proposed equivalent algorithm is that a fixed mesh of finite elements with a variable integration time-step can be employed to simulate the consequences and effects of the intruded magma solidification using the conventional finite element method. The correctness and usefulness of the proposed equivalent algorithm have been demonstrated by a benchmark magma solidification problem. Copyright (c) 2005 John Wiley & Sons, Ltd.

Fluid-mobile trace element constraints on the role of slab melting and implications for Archaean crustal growth models

We use published and new trace element data to identify element ratios which discriminate between arc magmas from the supra-subduction zone mantle wedge and those formed by direct melting of subducted crust (i.e. adakites). The clearest distinction is obtained with those element ratios which are strongly fractionated during refertilisation of the depleted mantle wedge, ultimately reflecting slab dehydration. Hence, adakites have significantly lower Pb/Nd and B/Be but higher Nb/Ta than typical arc magmas and continental crust as a whole. Although Li and Be are also overenriched in continental crust, behaviour of Li/Yb and Be/Nd is more complex and these ratios do not provide unique signatures of slab melting. Archaean tonalite-trondhjemite-granodiorites (TTGs) strongly resemble ordinary mantle wedge-derived arc magmas in terms of fluid-mobile trace element content, implying that they-did not form by slab melting but that they originated from mantle which was hydrated and enriched in elements lost from slabs during prograde dehydration. We suggest that Archaean TTGs formed by extensive fractional crystallisation from a mafic precursor. It is widely claimed that the time between the creation and subduction of oceanic lithosphere was significantly shorter in the Archaean (i.e. 20 Ma) than it is today. This difference was seen as an attractive explanation for the presumed preponderance of adakitic magmas during the first half of Earth's history. However, when we consider the effects of a higher potential mantle temperature on the thickness of oceanic crust, it follows that the mean age of oceanic lithosphere has remained virtually constant. Formation of adakites has therefore always depended on local plate geometry and not on potential mantle temperature.

Stratification of the earth beneath the Azores from P and S receiver functions

Seismic recordings of IRIS/IDA/GSN station CMLA and of several temporary stations in the Azores archipelago are processed with P and S receiver function (PRF and SRF) techniques. Contrary to regional seismic tomography these methods provide estimates of the absolute velocities and of the Vp/Vs ratio up to a depth of similar to 300 km. Joint inversion of PRFs and SRFs for a few data sets consistently reveals a division of the subsurface medium into four zones with a distinctly different Vp/Vs ratio: the crust similar to 20 km thick with a ratio of similar to 1.9 in the lower crust, the high-Vs mantle lid with a strongly reduced VpNs velocity ratio relative to the standard 1.8, the low-velocity zone (LVZ) with a velocity ratio of similar to 2.0, and the underlying upper-mantle layer with a standard velocity ratio. Our estimates of crustal thickness greatly exceed previous estimates (similar to 10 km). The base of the high-Vs lid (the Gutenberg discontinuity) is at a depth of-SO km. The LVZ with a reduction of S velocity of similar to 15% relative to the standard (IASP91) model is terminated at a depth of similar to 200 km. The average thickness of the mantle transition zone (TZ) is evaluated from the time difference between the S410p and SKS660p, seismic phases that are robustly detected in the S and SKS receiver functions. This thickness is practically similar to the standard IASP91 value of 250 km. and is characteristic of a large region of the North Atlantic outside the Azores plateau. Our data are indicative of a reduction of the S-wave velocity of several percent relative to the standard velocity in a depth interval from 460 to 500 km. This reduction is found in the nearest vicinities of the Azores, in the region sampled by the PRFs, but, as evidenced by SRFs, it is missing at a distance of a few hundred kilometers from the islands. We speculate that this anomaly may correspond to the source of a plume which generated the Azores hotspot. Previously, a low S velocity in this depth range was found with SRF techniques beneath a few other hotspots.

Mantle source heterogeneity, magma generation and magmatic evolution at Terceira Island (Azoresarchipelago): Constraints from elemental and isotopic (Sr, Nd, Hf, and Pb) data

This work addresses the present-day (<100 ka) mantle heterogeneity in the Azores region through the study of two active volcanic systems from Terceira Island. Our study shows that mantle heterogeneities are detectable even when "coeval" volcanic systems (Santa Barbara and Fissural) erupted less than 10 km away. These volcanic systems, respectively, reflect the influence of the Terceira and D. Joao de Castro Bank end-members defined by Beier et at (2008) for the Terceira Rift Santa Barbara magmas are interpreted to be the result of mixing between a HIMU-type component, carried to the upper mantle by the Azores plume, and the regional depleted MORB magmas/source. Fissural lavas are characterized by higher Ba/Nb and Nb/U ratios and less radiogenic Pb-206/Pb-204, Nd-143/Nd-144 and Hf-176/Hf-177, requiring the small contribution of delaminated sub-continental lithospheric mantle residing in the upper mantle. Published noble gas data on lavas from both volcanic systems also indicate the presence of a relatively undegassed component, which is interpreted as inherited from a lower mantle reservoir sampled by the ascending Azores plume. As inferred from trace and major elements, melting began in the garnet stability field, while magma extraction occurred within the spinel zone. The intra-volcanic system's chemical heterogeneity is mainly explained by variable proportions of the above-mentioned local end-members and by crystal fractionation processes. (C) 2011 Elsevier By. All rights reserved.

Deep structure of crust and mantle beneath Iberia and Western Mediterranean from P and S Receiver functions and SKS Waveforms

We have conducted a P and S receiver functions [PRFs and SRFs] analysis for 19 seismic stations on the Iberia and western Mediterranean. In the transition zone [TZ] the PRFs analysis reveals a band [from Gibraltar to Balearic] increased by 10-20 km relative to the standard 250 km. The TZ thickness variations are strongly correlated with the P660s times in PRFs. We interpret the variable depth of the 660-km discontinuity as an effect of subduction. Over the anomalous TZ we found a reduced velocity zone in the upper mantle. Joint inversion of PRFs and SRFs reveals a subcrustal high S velocity lid and an underlying LVZ. A reduction of the S velocity in the LVZ is less than 10%. The Gutenberg discontinuity is located at 65±5 km, but in several models sampling the Mediterranean, the lid is missing or its thickness is reduced to ~30 km. In the Gibraltar and North Africa this boundary is located at ~100 km. The lid Vp/Vs beneath Betics is reduced relative to the standard 1.8. Another evidence of the Vp/Vs anomaly is provided by S410p phase late arrivals in the SRFs. The azimuthal anisotropy analysis with a new technology was conducted at 5 stations and at 2 groups of stations. The fast direction in the uppermost mantle layer is ~90º in Iberian Massif. In Balearic is in the azimuth of ~120º. At a depth of ~60 km the direction becomes 90º. Anisotropy in the upper layer can be interpreted as frozen, whereas anisotropy in the lower layer is active, corresponding to the present-day or recent flow. The effect of the asthenosphere in the SKS splitting is much larger than the effect of the subcrustal lithosphere.

Age constraints on the late cretaceous alkaline magmatism on the West Iberian margin

Cretaceous Research 30 (2009) 575–586

From the Bay of Biscay to the High Atlas: completing the anisotropic characterization of the upper mantle beneath the westernmost Mediterranean region

The knowledge of the anisotropic properties beneath the Iberian Peninsula and Northern Morocco has been dramatically improved since late 2007 with the analysis of the data provided by the dense TopoIberia broadband seismic network, the increasing number of permanent stations operating in Morocco, Portugal and Spain, and the contribution of smaller scale/higher resolution experiments. Results from the two first TopoIberia deployments have evidenced a spectacular rotation of the fast polarization direction (FPD) along the Gibraltar Arc, interpreted as an evidence of mantle flow deflected around the high velocity slab beneath the Alboran Sea, and a rather uniform N100 degrees E FPD beneath the central Iberian Variscan Massif, consistent with global mantle flow models taking into account contributions of surface plate motion, density variations and net lithosphere rotation. The results from the last Iberarray deployment presented here, covering the northern part of the Iberian Peninsula, also show a rather uniform FPD orientation close to N100 degrees E, thus confirming the previous interpretation globally relating the anisotropic parameters to the LPO of mantle minerals generated by mantle flow at asthenospheric depths. However, the degree of anisotropy varies significantly, from delay time values of around 0.5 s beneath NW Iberia to values reaching 2.0 sin its NE comer. The anisotropic parameters retrieved from single events providing high quality data also show significant differences for stations located in the Variscan units of NW Iberia, suggesting that the region includes multiple anisotropic layers or complex anisotropy systems. These results allow to complete the map of the anisotropic properties of the westernmost Mediterranean region, which can now be considered as one of best constrained regions worldwide, with more than 300 sites investigated over an area extending from the Bay of Biscay to the Sahara platform. (C) 2015 Elsevier B.V. All rights reserved.

Mantle beneath the Gibraltar Arc from receiver functions

P and S receiver functions (PRF and SRF) from 19 seismograph stations in the Gibraltar Arc and the Iberian Massif reveal new details of the regional deep structure. Within the high-velocity mantle body below southern Spain the 660-km discontinuity is depressed by at least 20 km. The Ps phase from the 410-km discontinuity is missing at most stations in the Gibraltar Arc. A thin (similar to 50 km) low-S-velocity layer atop the 410-km discontinuity is found under the Atlantic margin. At most stations the S410p phase in the SRFs arrives 1.0-2.5 s earlier than predicted by IASP91 model, but, for the propagation paths through the upper mantle below southern Spain, the arrivals of S410p are delayed by up to +1.5 s. The early arrivals can be explained by elevated Vp/Vs ratio in the upper mantle or by a depressed 410-km discontinuity. The positive residuals are indicative of a low (similar to 1.7 versus similar to 1.8 in IASP91) Vp/Vs ratio. Previously, the low ratio was found in depleted lithosphere of Precambrian cratons. From simultaneous inversion of the PRFs and SRFs we recognize two types of the mantle: 'continental' and 'oceanic'. In the 'continental' upper mantle the S-wave velocity in the high-velocity lid is 4.4-4.5 km s(-1), the S-velocity contrast between the lid and the underlying mantle is often near the limit of resolution (0.1 km s(-1)), and the bottom of the lid is at a depth reaching 90 100 km. In the 'oceanic' domain, the S-wave velocities in the lid and the underlying mantle are typically 4.2-4.3 and similar to 4.0 km s(-1), respectively. The bottom of the lid is at a shallow depth (around 50 km), and at some locations the lid is replaced by a low S-wave velocity layer. The narrow S-N-oriented band of earthquakes at depths from 70 to 120 km in the Alboran Sea is in the 'continental' domain, near the boundary between the 'continental' and 'oceanic' domains, and the intermediate seismicity may be an effect of ongoing destruction of the continental lithosphere.

Magma flow pattern in dykes of the Azores revealed by ansiotropy of magnetic susceptibility

The localization of magma melting areas at the lithosphere bottom in extensional volcanic domains is poorly understood. Large polygenetic volcanoes of long duration and their associated magma chambers suggest that melting at depth may be focused at specific points within the mantle. To validate the hypothesis that the magma feeding a mafic crust, comes from permanent localized crustal reservoirs, it is necessary to map the fossilized magma flow within the crustal planar intrusions. Using the AMS, we obtain magmatic flow vectors from 34 alkaline basaltic dykes from São Jorge, São Miguel and Santa Maria islands in the Azores Archipelago, a hot-spot related triple junction. The dykes contain titanomagnetite showing a wide spectrum of solid solution ranging from Ti-rich to Ti-poor compositions with vestiges of maghemitization. Most of the dykes exhibit a normal magnetic fabric. The orientation of the magnetic lineation k1 axis is more variable than that of the k3 axis, which is generally well grouped. The dykes of São Jorge and São Miguel show a predominance of subhorizontal magmatic flows. In Santa Maria the deduced flow pattern is less systematic changing from subhorizontal in the southern part of the island to oblique in north. These results suggest that the ascent of magma beneath the islands of Azores is predominantly over localized melting sources and then collected within shallow magma chambers. According to this concept, dykes in the upper levels of the crust propagate laterally away from these magma chambers thus feeding the lava flows observed at the surface.

Depleted arc volcanism in the Alboran Sea and shoshonitic volcanism in Morocco: geochemical and isotopic constraints on Neogene tectonic processes

Samples of volcanic rocks from Alboran Island, the Alboran Sea floor and from the Gourougou volcanic centre in northern Morocco have been analyzed for major and trace elements and Sr-Nd isotopes to test current theories on the tectonic geodynamic evolution of the Alboran Sea. The Alboran Island samples are low-K tholeiitic basaltic andesites whose depleted contents of HFS elements (similar to0.5xN-MORB), especially Nb (similar to0.2xN-MORB), show marked geochemical parallels with volcanics from immature intra-oceanic arcs and back-arc basins. Several of the submarine samples have similar compositions, one showing low-Ca boninite affinity. Nd-143/Nd-144 ratios fall in the same range as many island-arc and back-arc basin samples, whereas Sr-87/Sr-86 ratios (on leached samples) are somewhat more radiogenic. Our data point to active subduction taking place beneath the Alboran region in Miocene times, and imply the presence of an associated back-arc spreading centre. Our sea floor suite includes a few more evolved dacite and rhyolite samples with (Sr-87/Sr-86)(0) up to 0.717 that probably represent varying degrees of crustal melting. The shoshonite and high-K basaltic andesite lavas from Gourougou have comparable normalized incompatible-element enrichment diagrams and Ce/Y ratios to shoshonitic volcanics from oceanic island arcs, though they have less pronounced Nb deficits. They are much less LIL- and LREE-enriched than continental arc analogues and post-collisional shoshonites from Tibet. The magmas probably originated by melting in subcontinental lithospheric mantle that had experienced negligible subduction input. Sr-Nd isotope compositions point to significant crustal contamination which appears to account for the small Nb anomalies. The unmistakable supra-subduction zone (SSZ) signature shown by our Alboran basalts and basaltic andesite samples refutes geodynamic models that attribute all Neogene volcanism in the Alboran domain to decompression melting of upwelling asthenosphere arising from convective thinning of over-thickened lithosphere. Our data support recent models in which subsidence is caused by westward rollback of an eastward-dipping subduction zone beneath the westemmost Mediterranean. Moreover, severance of the lithosphere at the edges of the rolling-back slab provides opportunities for locally melting lithospheric mantle, providing a possible explanation for the shoshonitic volcanism seen in northern Morocco and more sporadically in SE Spain. (C) 2004 Elsevier B.V. All rights reserved.

Characteristic reflection patterns in the southeast Canadian cordillera, northern Appalachians and Swiss Alps

There are some striking similarities and some differences between the seismic reflection sections recorded across the fold and thrust belts of the southeast Canadian Cordillera, Quebec-Maine Appalachians and Swiss Alps. In the fold and thrust belts of all three mountain ranges, seismic reflection surveys have yielded high-quality images of. (1) nappes (thin thrust sheets) stacked on top of ancient continental margins; (2) ramp anticlines in the hanging walls of faults that have ramp-flat or listric geometries; (3) back thrusts and back folds that developed during the terminal phases of orogeny; and (4) tectonic wedges and regional decollements. A principal result of the Cordilleran and Appalachian deep crustal studies has been the recognition of master decollements along which continental margin strata have been transported long distances, whereas a principal result of the Swiss Alpine deep crustal program has been the identification of the Adriatic indenter, a crustal-scale wedge that caused delamination of the European lithosphere. Significant crustal roots are observed beneath the fold and thrust belts of the Alps, southeast Canadian Cordillera and parts of the southern Appalachians, but such structures beneath the northern Appalachians have probably been removed by post-orogenic collapse and/or crustal attenuation associated with the Mesozoic opening of the Atlantic Ocean.

Predicting community and ecosystem outcomes of mycorrhizal responses to global change.

Mycorrhizal symbioses link the biosphere with the lithosphere by mediating nutrient cycles and energy flow though terrestrial ecosystems. A more mechanistic understanding of these plant-fungal associations may help ameliorate anthropogenic changes to C and N cycles and biotic communities. We explore three interacting principles: (1) optimal allocation, (2) biotic context and (3) fungal adaptability that may help predict mycorrhizal responses to carbon dioxide enrichment, nitrogen eutrophication, invasive species and land-use changes. Plant-microbial feedbacks and thresholds are discussed in light of these principles with the goal of generating testable hypotheses. Ideas to develop large-scale collaborative research efforts are presented. It is our hope that mycorrhizal symbioses can be effectively integrated into global change models and eventually their ecology will be understood well enough so that they can be managed to help offset some of the detrimental effects of anthropogenic environmental change.

Liquid line of descent of a basanitic liquid at 1.5 GPa: Constraints for metasomatic vein formation

The metasomatism observed in the oceanic and continental lithosphere is generally interpreted to represent a continuous differentiation process forming anhydrous and hydrous veins plus a cryptic enrichment in the surrounding peridotite. In order to constrain the mechanisms of vein formation and potentially clarify the nature and origin of the initial metasomatic agent, we performed a series of high-pressure experiments simulating the liquid line of descent of a basanitic magma differentiating within continental or mature oceanic lithosphere. This series of experiments has been conducted in an end-loaded piston cylinder apparatus starting from an initial hydrous ne-normative basanite at 1.5 GPa and temperature varying between 1,250 and 980°C. Near-pure fractional crystallization process was achieved in a stepwise manner in 30°C temperature steps and starting compositions corresponding to the liquid composition of the previous, higher-temperature glass composition. Liquids evolve progressively from basanite to peralkaline, aluminum-rich compositions without significant SiO2 variation. The resulting cumulates are characterized by an anhydrous clinopyroxene + olivine assemblage at high temperature (1,250-1,160°C), while at lower temperature (1,130-980°C), hydrous cumulates with dominantly amphibole + minor clinopyroxene, spinel, ilmenite, titanomagnetite and apatite (1,130-980°C) are formed. This new data set supports the interpretation that anhydrous and hydrous metasomatic veins could be produced during continuous differentiation processes of primary, hydrous alkaline magmas at high pressure. However, the comparison between the cumulates generated by the fractional crystallization from an initial ne-normative liquid or from hy-normative initial compositions (hawaiite or picrobasalt) indicates that for all hydrous liquids, the different phases formed upon differentiation are mostly similar even though the proportions of hydrous versus anhydrous minerals could vary significantly. This suggests that the formation of amphibole-bearing metasomatic veins observed in the lithospheric mantle could be linked to the differentiation of initial liquids ranging from ne-normative to hy-normative in composition. The present study does not resolve the question whether the metasomatism observed in lithospheric mantle is a precursor or a consequence of alkaline magmatism; however, it confirms that the percolation and differentiation of a liquid produced by a low degree of partial melting of a source similar or slightly more enriched than depleted MORB mantle could generate hydrous metasomatic veins interpreted as a potential source for alkaline magmatism by various authors.

On the origin of the ultradeep East Barents Sea basin

Very large subsidence, with up to 20 km thick sediment layers, is observed in the East Barents Sea basin. Subsidence started in early Paleozoic, accelerated in Permo-Triassic times, finished during the middle Cretaceous, and was followed by moderate uplift in Cenozoic times. The observed gravity signal suggests that the East Barents Sea is at present in isostatic balance and indicates that a mass excess is required in the lithosphere to produce the observed large subsidence. Several origins have been proposed for the mass excess. We use 1-D thermokinematic modeling and 2-D isostatic density models of continental lithosphere to evaluate these competing hypotheses. The crustal density in 2-D thermokinematic models resulting from pressure-, temperature-, and composition-dependent phase change models is computed along transects crossing the East Barents Sea. The results indicate the following. (1) Extension can only explain the observed subsidence provided that a 10 km thick serpentinized mantle lens beneath the basin center is present. We conclude that this is unlikely given that this highly serpentinized layer should be formed below a sedimentary basin with more than 10 km of sediments and crust at least 10 km thick. (2) Phase changes in a compositionally homogeneous crust do not provide enough mass excess to explain the present-day basin geometry. (3) Phase change induced densification of a preexisting lower crustal gabbroic body, interpreted as a mafic magmatic underplate, can explain the basin geometry and observed gravity anomalies. The following model is proposed for the formation of the East Barents Sea basin: (1) Devonian rifting and extension related magmatism resulted in moderate thinning of the crust and a mafic underplate below the central basin area explaining initial late Paleozoic subsidence. (2) East-west shortening during the Permian and Triassic resulted in densification of the previously emplaced mafic underplated body and enhanced subsidence dramatically, explaining the present-day deep basin geometry.