Geochemistry of the gabbroic lithosphere in IODP Holes 304-U1309D and 305-U1309D

IODP Hole U1309D (Atlantis Massif, Mid-Atlantic Ridge 30°N) is the second deepest hole drilled into slow spread gabbroic lithosphere. It comprises 5.4% of olivine-rich troctolites (~ > 70% olivine), possibly the most primitive gabbroic rocks ever drilled at mid-ocean ridges. We present the result of an in situ trace element study carried out on a series of olivine-rich troctolites, and neighbouring troctolites and gabbros, from olivine-rich intervals in Hole U1309D. Olivine-rich troctolites display poikilitic textures; coarse-grained subhedral to medium-grained rounded olivine crystals are included into large undeformed clinopyroxene and plagioclase poikiloblasts. In contrast, gabbros and troctolites have irregularly seriate textures, with highly variable grain sizes, and locally poikilitic clinopyroxene oikocrysts in troctolites. Clinopyroxene is high Mg# augite (Mg# 87 in olivine-rich troctolites to 82 in gabbros), and plagioclase has anorthite contents ranging from 77 in olivine-rich troctolites to 68 in gabbros. Olivine has high forsterite contents (82-88 in olivine-rich troctolites, to 78-83 in gabbros) and is in Mg-Fe equilibrium with clinopyroxene. Clinopyroxene cores and plagioclase are depleted in trace elements (e.g., Ybcpx ~ 5-11 * Chondrite), they are in equilibrium with the same MORB-type melt in all studied rock-types. These compositions are not consistent with the progressively more trace element enriched (evolved) compositions expected from olivine rich primitive products to gabbros in a MORB cumulate sequence. They indicate that clinopyroxene and plagioclase crystallized concurrently, after melts having the same trace element composition, consistent with crystallization in an open system with a buffered magma composition. The slight trace element enrichments and lower Cr contents observed in clinopyroxene rims and interstitial grains results from crystallization of late-stage differentiated melts, probably indicating the closure of the magmatic system. In contrast to clinopyroxene and plagioclase, olivine is not in equilibrium with MORB, but with a highly fractionated depleted melt, similar to that in equilibrium with refractory oceanic peridotites, thus possibly indicating a mantle origin. In addition, textural relationships suggest that olivine was in part assimilated by the basaltic melts after which clinopyroxene and plagioclase crystallized (impregnation). These observations suggest a complex crystallization history in an open system involving impregnation by MORB-type melt(s) of an olivine-rich rock or mush. The documented magmatic processes suggest that olivine-rich troctolites were formed in a zone with large magmatic transfer and accumulation, similar to the mantle-crust transition zone documented in ophiolites and at fast spreading ridges.

Strength of the Iberian intraplate lithosphere: Cenozoic deformations and seismicity

Previous studies about the strength of the lithosphere in the center of Iberia fail to resolve the depth of earthquakes because of the rheological uncertainties. Therefore, new contributions are considered (the crustal structure from a density model) and several parameters (tectonic regime, mantle rheology, strain rate) are checked in this paper to properly examine the role of lithospheric strength in the intraplate seismicity and the Cenozoic evolution. The strength distribution with depth, the integrated strength, the effective elastic thickness and the seismogenic thickness have been calculated by a finite element modelling of the lithosphere across the Central System mountain range and the bordering Duero and Madrid sedimentary basins. Only a dry mantle under strike-slip/extension and a strain rate of 10-15 s-1, or under extension and 10-16 s-1, causes a strong lithosphere. The integrated strength and the elastic thickness are lower in the mountain chain than in the basins. This heterogeneity has been maintained since the Cenozoic and determine the mountain uplift and the biharmonic folding of the Iberian lithosphere during the Alpine deformations. The seismogenic thickness bounds the seismic activity in the upper–middle crust, and the decreasing crustal strength from the Duero Basin towards the Madrid Basin is related to a parallel increase in Plio–Quaternary deformations and seismicity. However, elasto–plastic modelling shows that current African–Eurasian convergence is resolved elastically or ductilely, which accounts for the low seismicity recorded in this region.

Strength of the Iberian intraplate lithosphere: Cenozoic deformations and seismicity

Previous studies about the strength of the lithosphere in the Iberia centre fail to resolve the depth of earthquakes because of the rheological uncertainties. Therefore, new contributions are considered (the crustal structure from a density model) and several parameters (tectonic regime, mantle rheology, strain rate) are checked in this paper to properly examine the role of lithospheric strength in the intraplate seismicity and the Cenozoic evolution. The strength distribution with depth, the integrated strength, the effective elastic thickness and the seismogenic thickness have been calculated by a finite element modelling of the lithosphere across the Central System mountain range and the bordering Duero and Madrid sedimentary basins. Only a dry mantle under strike-slip/extension and a strain rate of 10-15 s-1, or under extension and 10-16 s-1, causes a strong lithosphere. The integrated strength and the elastic thickness are lower in the mountain chain than in the basins. These anisotropies have been maintained since the Cenozoic and determine the mountain uplift and the biharmonic folding of the Iberian lithosphere during the Alpine deformations. The seismogenic thickness bounds the seismic activity in the upper–middle crust, and the decreasing crustal strength from the Duero Basin towards the Madrid Basin is related to a parallel increase in Plio–Quaternary deformations and seismicity. However, elasto–plastic modelling shows that current African–Eurasian convergence is resolved elastically or ductilely, which accounts for the low seismicity recorded in this region.

Thermal models, lithosphere thickness and heat flow in South Portugal. Some comments about the subject

Reinterpretation of old heat flow data or use of new data and new techniques of detection of the temperature under the surface have conducted to new heat flow density values in some regions of the globe. The problem of ice melting in Greenland and Antarctica caught the public's attention to the importance of knowledge on heat flow values and thermal structure of the globe. In the last years, several models were presented trying to obtain lithosphere and Moho thickness of the Iberia Peninsula. The work we intend to present is related with the SW part of the Iberia Peninsula ( south of the Ossa Morena zone, South Portuguese Zone and Algarve). The results obtained show a decrease in the thickness of the crust and the lithosphere in this region. Density anomalies in the crust are also referred. I intend to make the connection between the results of these models and the heat flow thermal conductivity, heat production and geological data available for the region, trying to explain the results of heat flow density data obtained.

Lithosphere Thickness, Heat Flow and Moho Depth in the South of Portugal

Abstract In the last years, several models have been presented trying to obtain lithosphere and Moho thickness in the Iberian Peninsula, using data related to geoid elevation and topography, gravity, seismicity and thermal analysis. The results obtained show a decrease in the thickness of the crust and the lithosphere in the SW part of the Iberian Peninsula. Density anomalies in the crust are also referred. Data obtained in the region was collected and deviations from average values used were detected. In this work, models were made taking into account the specific characteristics of the region. Heat flow, thermal conductivity, heat production, topography, gravity, seismic and geological data available for the region, were used to adjust the model. The results show that this region is different from other parts of the Iberian Peninsula and a special attention must be given to it. This work shows the importance of trying to know and understand the thermal structure of the region.

Lithosphere Thickness, Heat Flow and Moho Depth in the South of Portugal

ABSTRACT In the last years, several models were presented trying to obtain lithosphere and Moho thickness in the Iberian Peninsula, using data related to geoid elevation and topography, gravity, seismicity and thermal analysis. The results obtained show a decrease in the thickness of the crust and the lithosphere in the SW part of the Iberian Peninsula. Density anomalies in the crust are also referred. The work I intend to present is related with the south of the Ossa Morena Zone, the South Portuguese Zone and the Algarve, in the south of Portugal. Data obtained in the region was collected and deviations from average values used were detected. Models were made taking into account the specific characteristics of the region. Heat flow, thermal conductivity, heat production, topography, gravity, seismic and geological data available for the region, are used to adapt the models. A special attention will be given to the spatial variation of heat flow values and to Moho depth in the region. The results show that this region is different from other parts of the Iberian Peninsula and a special attention must be given to it. The different values obtained using seismic, gravity, and geoid height data, and the results obtained with models using thermal data shows the importance of trying to know and understand the thermal structure of the regions. Problems related with the use of average values will be focused.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, 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 a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.