Rheology of the Earth and a thermoconvective mechanism for sedimentary basin formation

A power-law non-Newtonian fluid is usually assumed to model slow flows in the mantle and, in particular, convective flows. However, the power-law fluid has no memory, in contrast to a real material. A new non-linear integral (having a memory) model is proposed to describe the rheology of rocks. The model is consistent with the theory of simple fluids with fading memory and with laboratory studies of rock creep. The proposed model reduces to the power-law fluid model for stationary flows and to the Andrade model for flows associated with small strains. Stationary convection beneath continents has been studied by Fleitout & Yuen (1984 ), who used the power-law fluid model and obtained the cold immobile boundary layer (continental lithosphere). In a stability analysis of this layer, the Andrade model must be used. The analysis shows that the lithosphere is overstable (the period of oscillation is about 200  Ma). In the present study, it is suggested that these thermoconvective oscillations of the lithosphere are a mechanism for sedimentary basin formation. The vertical crustal movement in sedimentary basins can be considered as a slow subsidence on which small-amplitude oscillations are superimposed. The longest period of oscillatory crustal movement is of the same order of magnitude as the period of convective oscillation of the lithosphere found in the stability analysis. Taking into account the difference between depositional and erosional transport rates we can explain the permanent subsidence as well as the oscillations.     

Three-dimensional gravity inversion for Moho depth at rifted continental margins incorporating a lithosphere thermal gravity anomaly correction

This paper describes a method for determining Moho depth, lithosphere thinning factor (γ= 1 − 1/β) and the location of the ocean–continent transition at rifted continental margins using 3-D gravity inversion which includes a correction for the large negative lithosphere thermal gravity anomaly within continental margin lithosphere. The lateral density changes caused by the elevated geotherm in thinned continental margin and adjacent ocean basin lithosphere produce a significant lithosphere thermal gravity anomaly which may be in excess of −100 mGal, and for which a correction must be made in order to determine Moho depth accurately from gravity inversion. We describe a method of iteratively calculating the lithosphere thermal gravity anomaly using a lithosphere thermal model to give the present-day temperature field from which we calculate the lithosphere thermal density and gravity anomalies. For continental margin lithosphere, the lithosphere thermal perturbation is calculated from the lithosphere thinning factor (γ= 1 − 1/β) obtained from crustal thinning determined by gravity inversion and breakup age for thermal re-equilibration time. For oceanic lithosphere, the lithosphere thermal model used to predict the lithosphere thermal gravity anomaly may be conditioned using ocean isochrons from plate reconstruction models to provide the age and location of oceanic lithosphere. A correction is made for crustal melt addition due to decompression melting during continental breakup and seafloor spreading. We investigate the sensitivity of the lithosphere thermal gravity anomaly and the predicted Moho depth from gravity inversion at continental rifted margins to the methods used to calculate and condition the lithosphere thermal model using both synthetic models and examples from the North Atlantic.     

Lithospheric thickness beneath the Dabie Shan, central eastern China from S receiver functions

P and S receiver functions obtained from a portable array of 34 broad-band stations in east central China provide a detailed image of the crust–mantle and lithosphere–asthenosphere boundaries (LAB) in the Dabie Shan and its adjacent areas. Clear S -to- P converted waves produced at the LAB show a thin lithosphere beneath the whole study area. Based on our results, the thickest lithosphere of 72 km is observed beneath the southern part of the area within the Yangtze craton, whereas beneath the North-China platform, the lithosphere is only 60 km thick. S receiver functions also reveal, in good agreement with P receiver functions, a maximum depth of the Moho beneath the Dabie Shan orogen at approximately 40 km. Furthermore, we interpret the structural difference at 32° latitude as the probable location of the mantle suture formed between the Yangtze and the Sino-Korean cratons.     

Crustal thickness estimation beneath the southern central Andes at 30°S and 36°S from S wave receiver function analysis

We have used the S wave receiver function (SRF) technique to investigate the crustal thickness beneath two seismic profiles from the CHARGE project in the southern central Andes. A previous study employing the P wave receiver function method has observed the Moho interface beneath much of the profiles. They found, however, that the amplitude of the P to S conversion was diminished in the western part of the profiles and have attributed it to a reduction of the impedance contrast at the Moho due to lower crustal ecologitization. With SRF, we have successfully detected S to P converted waves from the Moho as well as possible conversions from other lithospheric boundaries. The continental South American crust reaches its maximum thickness of ∼70 km (along 30°S between 70°W and 68.5°W) beneath the Principal Cordillera and the Famatina system and becomes thinner towards the Sierras Pampeanas with a thickness of ∼40 km. Negative phases, possibly related to the base of the continental and oceanic lithosphere, can be recognized in the summation traces at different depths. By comparing our results with data obtained from previous investigations, we are able to further constrain the thickness of the crust and lithosphere beneath the central Andes.     

A three-dimensional radially anisotropic model of shear velocity in the whole mantle

We present a 3-D radially anisotropic S velocity model of the whole mantle (SAW642AN), obtained using a large three component surface and body waveform data set and an iterative inversion for structure and source parameters based on Non-linear Asymptotic Coupling Theory (NACT). The model is parametrized in level 4 spherical splines, which have a spacing of ∼ 8°. The model shows a link between mantle flow and anisotropy in a variety of depth ranges. In the uppermost mantle, we confirm observations of regions with   V_SH > V_SV   starting at ∼80 km under oceanic regions and ∼200 km under stable continental lithosphere, suggesting horizontal flow beneath the lithosphere. We also observe a   V_SV > V_SH   signature at ∼150–300 km depth beneath major ridge systems with amplitude correlated with spreading rate for fast-spreading segments. In the transition zone (400–700 km depth), regions of subducted slab material are associated with   V_SV > V_SH   , while the ridge signal decreases. While the mid-mantle has lower amplitude anisotropy (<1 per cent), we also confirm the observation of radially symmetric   V_SH > V_SV   in the lowermost 300 km, which appears to be a robust conclusion, despite an error in our previous paper which has been corrected here. The 3-D deviations from this signature are associated with the large-scale low-velocity superplumes under the central Pacific and Africa, suggesting that   V_SH > V_SV   is generated in the predominant horizontal flow of a mechanical boundary layer, with a change in signature related to transition to upwelling at the superplumes.     

Effect of 3-D viscoelastic structure on post-seismic relaxation from the 2004 M= 9.2 Sumatra earthquake

The 2004 M = 9.2 Sumatra–Andaman earthquake profoundly altered the state of stress in a large volume surrounding the ∼1400 km long rupture. Induced mantle flow fields and coupled surface deformation are sensitive to the 3-D rheology structure. To predict the post-seismic motions from this earthquake, relaxation of a 3-D spherical viscoelastic earth model is simulated using the theory of coupled normal modes. The quasi-static deformation basis set and solution on the 3-D model is constructed using: a spherically stratified viscoelastic earth model with a linear stress–strain relation; an aspherical perturbation in viscoelastic structure; a 'static' mode basis set consisting of Earth's spheroidal and toroidal free oscillations; a "viscoelastic" mode basis set; and interaction kernels that describe the coupling among viscoelastic and static modes. Application to the 2004 Sumatra–Andaman earthquake illustrates the profound modification of the post-seismic flow field at depth by a slab structure and similarly large effects on the near-field post-seismic deformation field at Earth's surface. Comparison with post-seismic GPS observations illustrates the extent to which viscoelastic relaxation contributes to the regional post-seismic deformation.     

Small-scale instabilities below the cooling oceanic lithosphere

Interpretation of satellite altimetry data as well as ship bathymetry data revealed strongly elongated anomalies roughly perpendicular to the mid-ocean ridges in the Indian and east Pacific oceans. A spectral analysis of gravity altimetry data along profiles parallel to the East Pacific Rise indicated wavelengths of about 150–180  km close to the ridge and about 250  km further away. A simple model of Rayleigh–Taylor instabilities developing at the base of the cooling lithosphere is discussed and applied to the data. By considering thermal diffusion and comparing Rayleigh–Taylor growth rates to the velocity of the thermal front in the cooling lithosphere, we are able to explain the observed anomalies by instabilities developing below the lithosphere in a layer with a viscosity of about 10¹⁹  Pa  s above an asthenospheric layer with a viscosity reduction of 2–3 orders of magnitude.     

Effect of the viscoelastic lithosphere on polar wander speed caused by the Late Pleistocene glacial cycles

Previous studies of the wander of the rotation pole associated with the Late Pleistocene glacial cycles indicate that the predicted polar wander speed is sensitive to the density jump at the 670 km discontinuity, the thickness of the elastic lithosphere, and the lower mantle viscosity. In particular, the M1 mode related to the density jump at 670 km depth has been shown to contribute a dominant portion of predicted polar wander speed for sufficiently small lower mantle viscosities. In this study, we examine the sensitivity of polar wander to variations in the viscosity of the viscoelastic lithosphere using simplified compressible Maxwell viscoelastic earth models. Model calculations for earth models with a viscoelastic lithosphere of finite viscosity indicate that the contribution of the M1 mode is similar to those associated with the density discontinuity at the core–mantle boundary (C0 mode) and the lithosphere (L0 mode). We speculate that this is due to the interaction between the M1 mode and the transient mode associated with the viscoelastic lithosphere, which reduces the magnitude of polar wander rates. Therefore, the M1 mode does not contribute a dominant portion of the predicted polar wander speed for earth models with a viscoelastic lithosphere of finite viscosity. In this case, predictions of polar wander speed as a function of lower mantle viscosity exhibit the qualitative form of an 'inverted parabola', as predicted for the J ˙₂ curve. We caution, however, that these results are obtained for simplified earth models, and the results for seismological earth models such as PREM may be complicated by the interaction between the M1 mode and the large set of transient modes.     

Polarization anomaly of Love waves caused by lateral heterogeneity

We calculate surface waves propagating in a laterally heterogeneous structure beneath the Kuril trench, where significant Love-wave polarization anomalies, called quasi-Love waves, are generated. Since 3-D wave propagation in the two-dimensionally heterogeneous structure can be assumed, we apply the 2.5-D finite difference method to the surface-wave calculations. The calculations show that a velocity contrast of 7 per cent at depths of less than 210 km beneath the Kuril trench cannot generate quasi-Love waves, and that an unlikely contrast of 20 per cent is required to generate clear quasi-Love waves. The possible cause of the quasi-Love waves inferred from previous studies on coupled free oscillations is a lateral variation in azimuthal anisotropy. The lateral variation in azimuthal anisotropy beneath the Kuril trench suggests a change in the mantle flow induced by the subducting slab.     

Post-seismic reloading and temporal clustering on a single fault

Geological studies show evidence for temporal clustering of large earthquakes on individual fault systems. Since post-seismic deformation due to the inelastic rheology of the lithosphere may result in a variable loading rate on a fault throughout the interseismic period, it is reasonable to expect that the rheology of the non-seismogenic lower crust and mantle lithosphere may play a role in controlling earthquake recurrence times. We study this phenomenon using a 2-D, finite element method continuum model of the lithosphere containing a single strike-slip fault. This model builds on a previous study using a 1-D spring-dashpot-slider analogue of a single fault system to study the role of Maxwell viscoelastic relaxation in producing non-periodic earthquakes. In our 2-D model, the seismogenic portion of the fault slips when a predetermined yield stress is exceeded; stress accumulated on the seismogenic fault is shed to the viscoelastic layers below and recycled back to the seismogenic fault through viscoelastic relaxation. We find that random variation of the fault yield stress from one earthquake to the next can cause the earthquake sequence to be clustered; the amount of clustering depends on a non-dimensional number, W , called the Wallace number defined as the standard deviation of the randomly varied fault yield stress divided by the effective viscosity of the system times the tectonic loading rate. A new clustering metric based on the bimodal distribution of interseismic intervals allows us to investigate clustering behaviour of systems over a wide range of model parameters and those with multiple viscoelastic layers. For models with   W ≥ 1  clustering increases with increasing W , while those with   W ≤ 1  are unclustered, or quasi-periodic.     

Teleseismic delay-time tomography of the upper mantle beneath southeastern India: imprint of Indo–Antarctica rifting

A 3-D P -velocity map of the crust and upper mantle beneath the southeastern part of India has been reconstructed through the inversion of teleseismic traveltimes. Salient geological features in the study region include the Archean Dharwar Craton and Eastern Ghat metamorphic belt (EGMB), and the Proterozoic Cuddapah and Godavari basins. The Krishna–Godavari basin, on the eastern coastal margin, evolved in response to the Indo–Antarctica breakup. A 24-station temporary network provided 1161 traveltimes, which were used to model 3-D P -velocity variation. The velocity model accounts of 80 per cent of the observed data variance. The velocity picture to a depth of 120 km shows two patterns: a high velocity beneath the interior domain (Dharwar craton and Cuddapah basin), and a lower velocity beneath the eastern margin region (EGMB and coastal basin). Across the array velocity variations of 7–10 per cent in the crust (0–40 km) and 3–5 per cent in the uppermost mantle (40–120 km) are observed. At deeper levels (120–210 km) the upper-mantle velocity differences are insignificant among different geological units. The presence of such a low velocity along the eastern margin suggests significantly thin lithosphere (<100 km) beneath it compared to a thick lithosphere (>200 km) beneath the eastern Dharwar craton. Such lithospheric thinning could be a consequence of Indo–Antarctica break-up.     

Impact of regional mantle flow on subducting plate geometry and interplate stress: insights from physical modelling

Physical models of subduction investigate the impact of regional mantle flow on the structure of the subducted slab and deformation of the downgoing and overriding plates. The initial mantle flow direction beneath the overriding plate can be horizontal or vertical, depending on its location with respect to the asthenospheric flow field. Imposed mantle flow produces either over or underpressure on the lower surface of the slab depending on the initial mantle flow pattern (horizontal or vertical, respectively). Overpressure promotes shallow dip subduction while underpressure tends to steepen the slab. Horizontal mantle flow with rates of 1–10 cm yr⁻¹ provides sufficient overpressure on a dense subducting lithosphere to obtain a subduction angle of  ∼60°  , while the same lithospheric slab sinks vertically when no flow is imposed. Vertical drag force (due to downward mantle flow) exerted on a slab can result in steep subduction if the slab is neutrally buoyant but fails to produce steep subduction of buoyant oceanic lithosphere. The strain regime in the overriding plate due to the asthenospheric drag force depends largely on slab geometry. When the slab dip is steeper than the interplate zone, the drag force produces negative additional normal stress on the interplate zone and tensile horizontal stress in the overriding plate. When the slab dip is shallower than the interplate zone, an additional positive normal stress is produced on the interplate zone and the overriding plate experiences additional horizontal compressive stress. However, the impact of the mantle drag force on interplate pressure is small compared to the influence of the slab pull force since these stress variations can only be observed when the slab is dense and interplate pressure is low.     

The Canary Islands swell: a coherence analysis of bathymetry and gravity

The Canary Archipelago is an intraplate volcanic chain, located near the West African continental margin, emplaced on old oceanic lithosphere of Jurassic age, with an extended volcanic activity since Middle Miocene. The adjacent seafloor does not show the broad oceanic swell usually observed in hotspot-generated oceanic islands. However, the observation of a noticeable depth anomaly in the basement west of the Canaries might indicate that the swell is masked by a thick sedimentary cover and the influence of the Canarian volcanism. We use a spectral approach, based on coherence analysis, to determine the swell and its compensation mechanism. The coherence between gravity and topography indicates that the swell is caused by a subsurface load correlated with the surface volcanic load. The residual gravity/geoid anomaly indicates that the subsurface load extends 600 km SSW and 800 km N and NNE of the islands. We used computed depth anomalies from available deep seismic profiles to constrain the extent and amplitude of the basement uplift caused by a relatively low-density anomaly within the lithospheric mantle, and coherence analysis to constrain the elastic thickness of the lithosphere ( T_e ) and the compensation depth of the swell. Depth anomalies and coherence are well simulated with T_e =28–36 km, compensation depth of 40–65 km, and a negative density contrast within the lithosphere of ∼33 kg m⁻³. The density contrast corresponds to a temperature increment of ∼325°C, which we interpret to be partially maintained by a low-viscosity convective layer in the lowermost lithosphere, and which probably involves the shallower parts of the asthenosphere. This interpretation does not require a significant rejuvenation of the mechanical properties of the lithosphere.     

P- and S-velocity images of the lithosphere–asthenosphere system in the Central Andes from local-source tomographic inversion

About 50 000 P and S arrival times and 25 000 values of t * recorded at seismic arrays operated in the Central Andes between 20°S and 25°S in the time period from 1994 to 1997 have been used for locating more than 1500 deep and crustal earthquakes and creating 3-D P , S velocity and Qp models. The study volume in the reference model is subdivided into three domains: slab, continental crust and mantle wedge. A starting velocity distribution in each domain is set from a priori information: in the crust it is based on the controlled sources seismic studies; in slab and mantle wedge it is defined using relations between P and S velocities, temperature and composition given by mineral physics. Each iteration of tomographic inversion consists of the following steps: (1) absolute location of sources in 3-D velocity model using P and S arrival times; (2) double-difference relocation of the sources and (3) simultaneous determination of P and S velocity anomalies, P and S station corrections and source parameters by inverting one matrix. Velocity parameters are computed in a mesh with the density of nodes proportional to the ray density with double-sided nodes at the domain boundaries. The next iteration is repeated with the updated velocity model and source parameters obtained at the previous step. Different tests aimed at checking the reliability of the obtained velocity models are presented. In addition, we present the results of inversion for Vp and Vp/Vs parameters, which appear to be practically equivalent to Vp and Vs inversion. A separate inversion for Qp has been performed using the ray paths and source locations in the final velocity model. The resulting Vp , Vs and Qp distributions show complicated, essentially 3-D structure in the lithosphere and asthenosphere. P and S velocities appear to be well correlated, suggesting the important role of variations of composition, temperature, water content and degree of partial melting.     

The shear-wave velocity structure in the upper mantle beneath Eurasia

We develop an approach that allows us to invert for the mantle velocity structure within a finely parametrized region as a perturbation with respect to a low-resolution, global tomographic model. We implement this technique to investigate the upper-mantle structure beneath Eurasia and present a new model of shear wave velocity, parametrized laterally using spherical splines with ∼2.9° spacing in Eurasia and ∼11.5° spacing elsewhere. The model is obtained from a combined data set of surface wave phase velocities, long-period waveforms and body-wave traveltimes. We identify many features as narrow as few hundred kilometres in diameter, such as subducting slabs in eastern Eurasia and slow-velocity anomalies beneath tectonically active regions. In contrast to regional studies in which these features have been identified, our model encompasses the structure of the entire Eurasian continent. Furthermore, including mantle- and body-wave waveforms helped us constrain structures at depths larger than 250 km, which are poorly resolved in earlier models. We find that up to +9 per cent faster-than-average anomalies within the uppermost ∼200 km of the mantle beneath cratons and some orogenic regions are separated by a sharp gradient zone from deeper, +1 to +2 per cent anomalies. We speculate that this gradient zone may represent a boundary separating the lithosphere from the continental root, which might be compositionally distinct from the overlying lithosphere and remain stable either due to its compositional buoyancy or due to higher viscosity compared with the suboceanic mantle. Our regional model of anisotropy is not significantly different from the global one.     

Deglacial land emergence and lateral upper-mantle heterogeneity in the Svalbard Archipelago—II. Extended results for high-resolution load models

In Paper I (Breuer & Wolf 1995), a preliminary interpretation of the postglacial land emergence observed at a restricted set of six locations in the Svalbard Archipelago was given. The study was based on a simple model of the Barents Sea ice sheet and suggested increases in lithosphere thickness and asthenosphere viscosity with increasing distance from the continental margin.
In the present paper, the newly developed high-resolution load model. BARENTS-2, and land-uplift observations from an extended set of 25 locations are used to study further the possibility of resolving lateral heterogeneity in the upper mantle below the northern Barents Sea. A comparison of the calculated and observed uplift values shows that the lithosphere thickness is not well resolved by the observations, although values above 110 km are most common for this parameter. In contrast to this, there are indications of a lateral variation of asthenosphere viscosity. Whereas values in the range 10¹⁸-10²⁰Pas are inferred for locations close to the continental margin, 10²⁰-10²¹ Pa s are suggested further away from the margin.
A study of the sensitivity of the values found for lithosphere thickness and asthenosphere viscosity to modifications of load model BARENTS-2 shows that such modifications can be largely accommodated by appropriate changes in lithosphere thickness, whereas the suggested lateral variation of asthenosphere viscosity is essentially unaffected. An estimate of the influence of the Fennoscandian. ice sheet leads to the conclusion that its neglect results in an underestimation of the thickness of the Barents Sea ice sheet by about 10 per cent.     

A perturbation theory for Love waves in anelastic media

Summary. We give a systematic formulation and a rigorous justification of a perturbation technique for the computation of the eigenvalues and eigenfunctions of Love waves (and toroidal oscillations by an appropriate change for variables) in an anelastic medium with a constitutive law modelling geophysical media of current interest such as the Kelvin—Voigt Solid, the Maxwell Solid, the Standard Linear Solid, and the Standard Linear Solid with a continuous spectrum of relaxation times. We develop expressions relating the eigenvalues of eigenfunctions for Love waves in a continuously varying vertically stratified anelastic half-space to the corresponding elastic eigenvalues and eigenfunctions. Analytically, our correspondence principle has the form of a regular perturbation expansion in terms of a parameter for both the eigenvalues and eigenfunctions. The identification of ɛ is motivated by the dissipativity principle of viscoelasticity theory. Moreover, we show that our correspondence principle applies respectively only in the high and low frequency range for the Maxwell and Kelvin—Voigt Solids. Outside of the applicable range of frequencies, our correspondence principle yields no useful information. For the family of Standard Linear Solids it is uniformly applicable for all non-zero frequencies.
We also derive an explicit formula to estimate the radius of convergence of our perturbation expansions. This estimate of the radius of convergence for each eigenvalue and eigenfunction is functionally defined by the constitutive model for the anelastic medium. The estimate is frequency dependent and depends on the separation distance between the eigenvalue and the remainder of the spectrum of the corresponding elastic problem.     

Late Jurassic subsidence and passive margin evolution in the Vulcan Sub-basin, north-west Australia: constraints from basin modelling

The Vulcan Sub-basin, located in the Timor Sea, north-west Australia, developed during the Late Jurassic extension which ultimately led to Gondwanan plate breakup and the development of the present-day passive continental margin. This paper describes the evolution of upper crustal extension and the development of Late Jurassic depocentres in this subbasin, via the use of forward modelling techniques. The results suggest that a lateral variation in structural style exists. The south of the basin is characterized by relatively large, discrete normal faults which have generated deep sub-basins, whereas more distributed, small-scale faulting further north reflects a collapse of the early basin margin, with the development of a broader, 'sagged' basin geometry. By combining forward and reverse modelling techniques, the degree of associated lithosphere stretching can be quantified. Upper crustal faulting, which represents up to 10% extension, is not balanced by extension in the deeper, ductile lithosphere; the magnitude of this deeper extension is evidenced by the amount of post-Valanginian thermal subsidence. Reverse modelling shows that the lithosphere stretching
factor has a magnitude of up to β=1.55 in the southern Vulcan Sub-basin, decreasing to β=1.2 in the northern Vulcan Sub-basin. It is proposed that during plate breakup, deformation in the Vulcan Sub-basin consisted of depth-dependent lithosphere extension. This additional component of lower crustal and lithosphere stretching is considered to reflect long-wavelength partitioning of strain associated with continental breakup, which may have extended 300–500 km landward of the continent–ocean boundary.     

Seismic anisotropy in the core–mantle transition zone

Split S waves observed at Hockley, Texas from events in the Tonga–Fiji region of the southwest Pacific show predominantly vertically polarized shear-wave ( SV  ) energy arriving earlier than horizontally polarized ( SH ) energy for rays propagating horizontally through D" . After corrections are made for the effects of upper-mantle anisotropy beneath Hockley, a time lag of 1.5 to 2.0  s remains for the furthest events (93.9°–100.6° ), while the time lags of the nearer observations (90.5°–92.9° ) nearly disappear. At closer distances, the S waves from these same events do not penetrate as deeply into the lower mantle, and are not split. These observations suggest that a patch of D" beneath the central Pacific is anisotropic, while the mantle immediately above the patch is isotropic. The thickness of the anisotropic zone appears to be of the order of 100–200  km.
  Observations of shear-wave splitting have previously been made for paths that traverse D" under the Caribbean and under Alaska. SH leads SV , the reverse of the Hockley observations, but in these areas the fact that SV  leads SH in the HKT data shown here suggests a different sort of anisotropy under the central Pacific from that under Alaska and the Caribbean. The case of SH travelling faster than SV  is consistent with transverse isotropy with a vertical axis of symmetry (VTI) and does not require variations with azimuth. The case of SV  leading SH is consistent with transverse isotropy with a horizontal axis of symmetry (HTI), an azimuthally anisotropic medium, and with a VTI medium formed by a hexagonal crystal. Given that (Mg,Fe)SiO₃ perovskite appears unlikely to form anisotropic fabrics on a large scale, the presence of anisotropy may point to chemical heterogeneity in the lowermost mantle, possibly due to mantle–core interactions.     

Sediment blanketing and the flexural strength of extended continental lithosphere

The flexural rigidity of the oceanic lithosphere is strongly dependent on its temperature structure at the time of loading. It is commonly assumed that the depth to the 450°C isotherm defines the effective elastic thickness T_e of the lithosphere. However, recent gravity studies across the Baltimore Canyon and Nova Scotian margins suggest that temperature may play a more complicated role in controlling the mechanical strength of extended continental lithosphere. For example, the flexural strength of the Baltimore Canyon margin (with sediment thicknesses of ? 15 km) appears to be controlled by the depth to the 150°C isotherm whereas the strength of the Nova Scotian margin (with sediment thicknesses cf ? 10 km) is controlled by the depth to the 250°C isotherm. The apparent correlation between sediment thickness and controlling isotherm suggests that sediment blanketing may play a role in modifying the flexural strength of extended continental lithosphere. This hypothesis was investigated by simulating the sedimentation history of a margin as a Gaussian function in which sedimentation peak and rate are determined by the mean and standard deviation of the function. The temperature structure of the lithosphere is continually modified as sediments are deposited on, and incorporated into the temperature structure of, the underlying lithosphere. Given a ‘starting’ value of T_e defined by the degree of extension of the lithosphere, the modification of T_e appears to be directly proportional to the sedimentation rate and cumulative sediment thickness, and inversely proportional to the time at which the sedimentation rate is a maximum. The first-order consequence of sediment blanketing is to reduce the cooling rate of the lithosphere relative to cooling in the absence of sediments. At thermal equilibrium, the initial value of T_e is reduced by the cumulative sediment thickness. Local isostatic conditions (i. e. T_e? 0) can only be approached when the sedimentation rate is unrealistically high (> 1000 m/Myr) during the rift or early post-rift phase of basin development. However, while these early loads may be locally compensated, any subsequent loads will be regionally compensated. Thus, it is unlikely that the low present-day flexural strengths interpreted from the Baltimore Canyon and Nova Scotian passive continental margins are a consequence of sediment blanketing.