[1] A combined geophysical-petrological methodology to study the thermal, compositional, density, and seismological structure of lithospheric/sublithospheric domains is presented. A new finite-element code (LitMod) is used to produce 2-D forward models from the surface to the 410-km discontinuity. The code combines data from petrology, mineral physics, and geophysical observables within a self-consistent framework. The final result is a lithospheric/sublithospheric model that simultaneously fits all geophysical observables and consequently reduces the uncertainties associated with the modeling of these observables alone or in pairs, as is commonly done. The method is illustrated by applying it to both oceanic and continental domains. We show that anelastic attenuation and uncertainties in seismic data make it unfeasible to identify compositional variations in the lithospheric mantle from seismic studies only. In the case of oceanic lithosphere, plates with thermal thicknesses of 105 ± 5 km satisfy geophysical and petrological constraints. We find that Vp are more sensitive to phase transitions than Vs, particularly in the case of the spinel-garnet transition. A low-velocity zone with absolute velocities and gradients comparable to those observed below ocean basins is an invariable output of our oceanic models, even when no melt effects are included. In the case of the Archean subcontinental lithospheric mantle, we show that ''typical'' depleted compositions (and their spatial distribution) previously thought to be representative of these mantle sections are compatible neither with geophysical nor with petrological data. A cratonic keel model consisting of (1) strongly depleted material (i.e., dunitic/harzburgitic) in the first 100-160 km depth and (2) less depleted (approximately isopycnic) lower section extending down to 220-300 km depth is necessary to satisfy elevation, geoid, SHF, seismic velocities, and petrological constraints. This highly depleted (viscous) upper layer, and its chemical isolation, may play a key role in the longevity and stability of cratons.

No abstract

We discuss the implications of a lithospheric model of the Moroccan Atlas Mountains based on topography, heat flow, gravity and geoid anomalies, taking into account the regional geology. The NW African cratonic lithosphere, some 160-180 km thick, thins to c. 80 km beneath the Atlas fold-thrust belts, in contrast with the shortening regime prevailing there since the early Cenozoic. This fact explains several geological and geophysical features as high topography with modest tectonic shortening, the occurrence of alkaline magmatism contemporaneous to compression, the absence of large crustal roots to support elevation, the scarce development of foreland basins, and a marked geoid high. The modelled lithosphere thinning is related to a thermal upwelling constrained between the Iberia-Africa convergent plate boundary and the Saharan craton.

We present a two‐dimensional lithospheric thermal and density model along a transect running from the southwestern Iberian Peninsula to the northwestern Sahara. The main goal is to investigate the lithosphere structure underneath the Gulf of Cadiz and the Atlas Mountains. The model is based on the assumption of topography in local isostatic equilibrium and is constrained by surface heat flow, gravity anomalies, geoid, and topography data. The crustal structure has been constrained by seismic and geological data where available. Mantle density is supposed to vary linearly with temperature, providing the link between thermal and density‐related data. The lithospheric thickness varies strongly along the profile, going from near 100 km under the Iberian Peninsula to at least 160–190 km under the Gulf of Cadiz and the Gharb foreland basin in Morocco and to 70 km underneath the Atlas Mountains, coinciding with a region of Neogene volcanism. The thickening of the lithosphere is interpreted as a SW trending lithospheric slab extending from the western Betics to the Gulf of Cadiz and the Gharb Basin, whereas the thin lithosphere underneath the Atlas may be interpreted as plume‐like asthenospheric upwelling similar to those observed in the west European Alpine foreland or as a side effect of a slab penetrating the less viscous asthenosphere.

Abstract. A three-dimensional gravity modeling combined with integrated heat flow and elevation modeling is conducted to map out the crustal and lithospheric mantle thickness in the Alboran Basin, in the westernmost Mediterranean. A "sediment"-corrected Bouguer anomaly has been derived using a depth-to-the-basement map and densities determined from well logs and seismic data. The gravity effect of the base of the lithosphere has been removed from the sediment-corrected Bouguer anomaly to obtain a "crustal" Bouguer anomaly, which has been inverted for crustal thickness. The resulting lithospheric structure is further constrained by elevation data under the assumption of local isostasy. The low residual elevation anomalies obtained (_+100 m in average) suggest that the area is in local isostasy, particularly the medium-and long-wavelength topography features.

No abstract

This paper presents a new southern North Atlantic plate model from Late Cretaceous to present, with the aim of constraining the kinematics of the Iberian plate during the last 83.5 Myr. This model is presented along with a detailed isochron map generated through the analysis of 3 aeromagnetic tracks and ~400 ship tracks from the National Centers for Environmental Information database. We present a new technique to obtain well‐constrained estimates of the Iberia‐North America plate motions from magnetic anomalies, overcoming the scarcity of large‐offset fracture zones and transform faults. We build an integrated kinematic model for NW Africa, Morocco, Iberia, Europe, and North America, which shows that the deformation is partitioned between Pyrenees and Betic‐Rif orogenic domain during the Late Cretaceous‐Oligocene time interval. In the Eastern Betics domain, the calculated amount of NW Africa‐Iberia convergence is ~80 km between 83.5 and 34 Ma, followed by ~150 km since the Oligocene. The motion of Iberia relative to Europe in the Central Pyrenees is characterized by overall NE directed transpressional motion during the Campanian and the Paleocene, followed by NW directed transpressional movement until the Lutetian and overall NNW directed convergence from Bartonian to Chattian. This motion occurs along the axis of the Bay of Biscay from the Santonian–Campanian boundary to the middle Priabonian, subsequently jumping to King's Trough at Anomaly 17 (36.62 Ma).

A two‐dimensional algorithm to determine the steady state thermal structure of the lithosphere that integrates thermal, gravity, and local isostasy analyses is presented. Gravity analyses together with seismic data are used to constrain spatial variations in density and crustal structure, while absolute elevation is used to determine the lithospheric mantle thickness. The calculation is performed using a finite element technique that links the different physical equations. The program optionally calculates the temperature at any material boundary and, with given rheological parameters, the strength distribution and the total lithospheric strength in selected columns. We apply the algorithm to the Northeastern Spanish Geotransect which extends from the Pyrenees to the Balearic Promontory and along which a strong variation in crustal and lithospheric thickness is evident. We assess the use of two different inferred density models for the lithospheric mantle: The first assumes a linear decrease in density with increasing temperature using the asthenospheric density as a reference; the second model assumes a constant density for the whole lithospheric mantle. Although conceptually the two hypotheses differ substantially, the results obtained do not show significant differences. Lithospheric thicknesses of 120–130 km below the Pyrenees, 60–65 km in the Valencia Trough, and 65–75 km below the Balearic Promontory are deduced. In all cases the mean lithospheric mantle density has to be 40–60 kg m−3 higher than the asthenospheric density. The algorithm is shown to be a powerful tool in lithospheric thermal modeling especially in areas where surface heat flow is poorly constrained because of the temperature‐density‐elevation relationship.