Insights into the lithospheric structure and tectonic setting of the Barents Sea region from isostatic considerations

We study the tectonic setting and lithospheric structure of the greater Barents Sea region by investigating its isostatic state and its gravity field. 3-D forward density modelling utilizing available information from seismic data and boreholes shows an apparent shift between the level of observed and modelled gravity anomalies. This difference cannot be solely explained by changes in crustal density. Furthermore, isostatic calculations show that the present crustal thickness of 35–37 km in the Eastern Barents Sea is greater than required to isostatically balance the deep basins of the area (>19 km). To isostatically compensate the missing masses from the thick crust and deep basins and to adequately explain the gravity field, high-density material (3300–3350 kg m⁻³) in the lithospheric mantle below the Eastern Barents Sea is needed. The distribution of mantle densities shows a regional division between the Western and Eastern Barents and Kara Seas. In addition, a band of high-densities is observed in the lower crust along the transition zone from the Eastern to Western Barents Sea. The distribution of high-density material in the crust and mantle suggests a connection to the Neoproterozoic Timanide orogen and argues against the presence of a Caledonian suture in the Eastern Barents Sea. Furthermore, the results indicate that the basins of the Western Barents Sea are mainly affected by rifting, while the Eastern Barents Sea basins are located on a stable continental platform.     

Crustal and upper-mantle structure in the Eastern Mediterranean from the analysis of surface wave dispersion curves

The dispersive properties of surface waves are used to infer earth structure in the Eastern Mediterranean region. Using group velocity maps for Rayleigh and Love waves from 7 to 100 s, we invert for the best 1-D crust and upper-mantle structure at a regular series of points. Assembling the results produces a 3-D lithospheric model, along with corresponding maps of sediment and crustal thickness. A comparison of our results to other studies finds the uncertainties of the Moho estimates to be about 5 km. We find thick sediments beneath most of the Eastern Mediterranean basin, in the Hellenic subduction zone and the Cyprus arc. The Ionian Sea is more characteristic of oceanic crust than the rest of the Eastern Mediterranean region as demonstrated, in particular, by the crustal thickness. We also find significant crustal thinning in the Aegean Sea portion of the backarc, particularly towards the south. Notably slower S -wave velocities are found in the upper mantle, especially in the northern Red Sea and Dead Sea Rift, central Turkey, and along the subduction zone. The low velocities in the upper mantle that span from North Africa to Crete, in the Libyan Sea, might be an indication of serpentinized mantle from the subducting African lithosphere. We also find evidence of a strong reverse correlation between sediment and crustal thickness which, while previously demonstrated for extensional regions, also seems applicable for this convergence zone.     

Small-scale mantle flows and induced lithospheric stress near island arcs

Summary. This paper explores the middle ground between complex thermally-coupled viscous flow models and simple corner flow models of island arc environments. The calculation retains the density-driven nature of convection and relaxes the geometrical constraints of corner flow, yet still provides semianalytical solutions for velocity and stress. A novel aspect of the procedure is its allowance for a coupled elastic lithosphere on top of a Newtonian viscous mantle. Initially, simple box-like density drivers illustrate how vertical and horizontal forces are transmitted through the mantle and how the lithosphere responds by trench formation. The flexural strength of the lithosphere spatially broadens the surface topography and gravity anomalies relative to the functional form of the vertical flow stresses applied to the plate base. I find that drivers in the form of inclined subducting slabs cannot induce self-driven parallel flow; however, the necessary flow can be provided by supplying a basal drag of 1–5 MPa to the mantle from the oceanic lithosphere. These basal drag forces create regional lithospheric stress and they should be quantifiable through seismic observations of the neutral surface. The existence of a shallow elevated phase transition is suggested in two slab models of 300 km length where a maximum excess density of 0.2 g cm⁻³ was needed to generate an acceptable mantle flow. A North New Hebrides subduction model which satisfies flow requirements and reproduces general features of topography and gravity contains a high shear stress zone (75 MPa) around the upper slab surface to a depth of 150 km and a deviatoric tensional stress in the back arc to a depth of 70 km. The lithospheric stress state of this model suggests that slab detachment is possible through whole plate fracture.     

The gravity field of a dipping plate in Greece

Summary. The surface gravity anomalies for a three-dimensional density model of a dipping lithospheric plate under the Aegean island arc and the Aegean Sea in Greece in the eastern Mediterranean Sea have been calculated. Such a dipping plate has been suggested in geophysical investigations. The calculated gravity maximum over the dipping plate is located in the area of the largest observed Bouguer anomalies inside the island arc. It is suggested that this should be taken into account in studies of crustal structure and gravity in the area.     

Travel-time anomalies in the mantle under the North Atlantic

Summary. Lateral heterogeneity exists in the Earth's mantle, and may result in seismic velocity anomalies up to several per cent. If convection cells and plumes extend down to the core, then these features may be associated with local inhomogeneities observed in the lower mantle.
Published data for direct and core-reflected P -wave residuals are used to delineate velocity anomalies in the middle—lower mantle under the North Atlantic. Differential ( PcP — P ) residuals indicate travel-time anomalies near the core—mantle transition, and may be due to core topography or lateral variations in velocity. It is assumed that the anomalies occur near the midpoints of the ray paths. The main source of error in the data set may arise from phases which have been identified incorrectly. Hence trend-surfaces are fitted to the residual data to show only the large-scale trends in anomaly values, with wavelengths of the order of 1000 km.
The Azores and Colorado hot spots occur in a region covered by the data. A possible interpretation of the trend maps is that an anomalous zone extends from a relatively fast region at the core boundary at 35° N, 50° W up to these hot spots, at about 30 degrees from the vertical. This may agree with the suggestion of Anderson that plumes are chemical rather than thermal in origin. If inclⁱned plumes do exist, the deviation from the ideal vertical plume or convection cell boundary may imply that lateral shear or other distortion effects exist in the mantle.     

Tectonic evolution of the North Sea basin: crustal stretching and subsidence

Summary. The lithospheric stretching model for the formation of sedimentary basins was tested in the central North Sea by a combined study of crustal thinning and basement subsidence patterns. A profile of crustal structure was obtained by shooting a long-range seismic experiment across the Central Graben, the main axis of subsidence. A seabed array of 12 seismometers in the graben was used to record shots fired in a line 530 km long across the basin. The data collected during the experiment were interpreted by modelling synthetic seismograms from a laterally varying structure, and the final model showed substantial crustal thinning beneath the graben. Subsidence data from 19 exploration wells were analysed to obtain subsidence patterns in the central North Sea since Jurassic times. Changes in water depth were quantified using foraminiferal assemblages where possible, and observed basement subsidence paths were corrected for sediment loading, compaction and changes in water depth through time. The seismic model is shown to be compatible with the observed gravity field, and the small size of observed gravity anomalies is used to argue that the basin is in local isostatic equilibrium. Both crustal thinning and basement subsidence studies indicate about 70 km of stretching across the Central Graben during the mid-Jurassic to early Cretaceous extensional event. This extension appears to have occurred over crust already slightly thinned beneath the graben, and the seismic data suggest that total extension since the early Permian may have been more than 100km. The data presented here may all be explained using a simple model of uniform extension of the lithosphere.     

Gravity studies of the Rockall and Exmouth Plateaux using SEASAT altimetry

Abstract SEASAT altimetric measurements are used to determine the gravity anomalies across two passive continental margins: the western margin of the Rockall Plateau, UK, and the Exmouth Plateau off north-west Australia. The small gravity anomalies observed over the starved western margin of the Rockall Plateau require the existence of a major density contrast within the crust, as well as the Moho, and show that the elastic thickness is less than 5 km at the time of rifting. The gravity anomaly over the Exmouth Plateau is compared with the gravity anomaly calculated from the sediment loading of a thin elastic plate, taking account of the variation in crustal thickness. This comparison shows that the Exmouth Plateau also has a small effective elastic thickness of 5 km, even for loads emplaced between 60 and 120 Myr after rifting. Elastic thicknesses of about 5 km have also been reported for other sedimentary basins, and are to be expected if the rheological properties of the crust and mantle depend on the ratio of the present temperature to the melting temperature. Flexural effects are therefore likely to be of minor importance in sedimentary basins.     

A comparison of the isostatic response of bathymetric features in the north Pacific Ocean and Philippine Sea

Summary. This paper concerns the calculation and analysis of admittance functions from large and uniform data sets of gravity and topography in four regions of the northern and western Pacific Ocean. The purpose is to separate and describe possible differences in isostatic compensation between several 'type' regions of oceanic crust: a mid-ocean ridge (Juan de Fuca), a mid-plate seamount chain (Hawaiian Ridge), fracture zone topography on old crust (north of Hawaii) and a marginal basin (Philippine Sea). Results suggest that there are significant differences in the degree to which long wavelength topography has been compensated which can be distinguished between regions. These differences are set in the perspective of three simple compensation mechanisms. Two of these consider local Airy models in which raised topography is compensated at depth either by crustal roots or low density mantle. A third considers the effects of an elastic plate of variable thickness supporting crustal variations. Conclusions are that: (a) a thick plate possibly in excess of 30 km supports the Hawaiian Ridge; (b) a much thinner plate of 5 to 15 km existed when the fracture zone topography was formed; (c) the Juan de Fuca Ridge is compensated either regionally by a plate 5 to 10 km thick or locally by sub-crustal low densities at depths of 15 to 20 km; and (d) the Philippine Sea shows no evidence for regional support: ridges are compensated locally by differences in crustal thickness whereas the basins are underlain by density variations at depths comparable to those of the much younger Juan de Fuca Ridge. The major difference between admittance functions for the Philippine Sea and comparably aged regions of the north Pacific Ocean adds further new evidence of possible evolutionary differences between it and normal ocean basins.     

Gravity constraints on the dynamics of the crust–mantle system during Calabrian subduction

The thermomechanic evolution of the lithosphere–upper mantle system during Calabrian subduction is analysed using a 2-D finite element approach, in which the lithosphere is compositionally stratified into crust and mantle. Gravity and topography predictions are cross-checked with observed gravity and topography patterns of the Calabrian region. Modelling results indicate that the gravity pattern in the arc-trench region is shaped by the sinking of light material, belonging to both the overriding and subduction plates. The sinking of light crustal material, up to depths of the order of 100–150 km is the ultimate responsible for the peculiar gravity signature of subduction, characterized by a minimum of gravity anomaly located at the trench, bounded by two highs located on the overriding and subducting plates, with a variation in magnitude of the order of 200 mGal along a wavelength of 200 km, in agreement with the isostatically compensated component of gravity anomaly observed along a transect crossing the Calabrian Arc, from the Tyrrhenian to the Ionian Seas. The striking agreement between the geodetic retrieved profiles and the modelled ones in the trench region confirms the crucial role of compositional stratification of the lithosphere in the subduction process and the correctness of the kinematic hypotheses considered in our modelling, that the present-day configuration of crust–mantle system below the Calabrian arc results from trench's retreat at a rate of about 3 cm yr⁻¹, followed by gravitational sinking of the subducted slab in the last 5 Myr.     

Crustal structure and isostatic compensation near the Kane fracture zone from topography and gravity measurements—I. Spectral analysis approach

Summary. Six gravity and bathymetry profiles perpendicular to the Kane fracture zone, each more than 300 km long, were gathered to study the variation in crustal structure in the vicinity of a major fracture zone and the gravitational edge effect at the contact between lithosphere of two different ages. A spectral analysis of the gravity and bathymetric series as a function of wavelength shows that the gravitational edge effect is only significant at the longest wavelengths. For remaining wavelengths the admittance, the ratio of the amplitude of the gravity anomaly to the amplitude of the bathymetry, is best explained by a model of isostasy in which topographic loads are partially supported by the flexural rigidity of an elastic plate, about 6 km in thickness. After subtracting the gravitational attraction of the bathymetry and its compensation, substantial isostatic anomalies remain. We interpret these anomalies as being caused by variations in crustal thickness which have little correlation with surface topography, except at very long wavelengths. The apparent crustal thickness varies by as much as a factor of 2, but there is no evidence indicating systematic thinning of the crust beneath the fracture zone. Our data do suggest that such density variations within the plate are also compensated by the isostatic response of an elastic plate but with very different effect from those at the surface. This indicates that there are two different modes of crustal formation with different gravity and topographic signatures: effusive volcanism which loads the surface of the elastic plate producing both topographic relief and coherent gravity anomalies, and intrusive volcanism or underplating producing gravity anomalies but little topographic relief.     

Lithospheric structure of the Arabian and Eurasian collision zone in eastern Turkey from S-wave receiver functions

Crustal and upper-mantle seismic discontinuities beneath eastern Turkey are imaged using teleseismic S -to- P converted phases. Three crustal phases are observed: the Moho with depth ranging between 30 and 55 km, indicating variable tectonic regimes within this continental collision zone; an upper-crustal discontinuity at approximately 10 km depth; and various crustal low-velocity zones, possibly associated with recent Quaternary volcanism. Imaging of the upper mantle is complicated by the 3-D geometry of the region, in particular due to the Bitlis–Zagros suture zone. However, several upper-mantle S -to- P converted phase are identified as being the signature of the lithosphere–asthenosphere boundary (LAB). The inferred LAB for the Eastern Anatolian Accretionary Complex indicates that eastern Turkey has an anomalously thin (between ∼60 and 80 km) lithosphere which is consistent with an oceanic slab detachment model. The observed LAB phases for the Arabian shield and Iranian plateau indicate that lithospheric thickness for these stable regions is on the order of 100 to 125 km thick, which is typical of continental margins.     

A Teleseismic Delay Time Study Across the Central African Shear Zone In the Adamawa Region of Cameroon, West Africa

Summary. P -wave relative teleseismic residuals were measured for a network of seismological stations along a 300 km profile across the Adamawa Plateau and the Central African Shear Zone of central Cameroon, to determine the variation in crust and upper mantle velocity associated with these structures. A plot of the mean relative residuals for the stations shows a long wavelength (> 300 km) variation of amplitude 0.45 s. the slowest arrivals are located over and just to the north, of the faulted northern margin of the Adamawa Plateau. the residuals do not correlate with topography, surface geology or the previously determined crustal structure, in any simple way.
The Aki inversion technique has been used to invert the relative residuals into a 3-D model of velocity perturbations from a mean earth model. the results show the region is divided roughly into three blocks by two subvertical boundaries, striking ENE and traversing both the crust and upper mantle down to depths greater than 190km. the central block, which is 2 per cent slower than the adjacent blocks, roughly corresponds to the Central African Shear Zone. the Adamawa Plateau, as an individual uplifted area, is explained by the interaction of a regional anomalous upper mantle associated with the West African Rift System, and the Central African Shear Zone, which provided a conduit for heat flow to the surface.     

Seismic structure of the upper mantle beneath the United States by three-dimensional inversion of body wave arrival times

Summary. Teleseismic P and S arrival times to North American stations are obtained from the ISC bulletins for the 10-yr period 1964–73, and relative travel-time delays are calculated with respect to standard tables. Station anomalies as well as variations of the delays with azimuth and epicentral distance from station are analysed, and the location of the velocity anomalies responsible for them is discussed. Inversion of the P delays to infer upper mantle velocity structure down to a depth of 700 km is obtained using three-dimensional blocks, as proposed by Aki, Christofferson & Husebye. Three layers can be resolved in this depth range. It is found that the heterogeneities responsible for the travel-time delays are primarily located in the first 250 km of the upper mantle, and that they correlate with surface features. Significant heterogeneities subsist to depths of at least 700 km and their broad scale pattern also correlates with the surface features: in the third layer (500 to 700 km depth) there is an increase of velocity from the West to the East of the United States, while the second layer (250 to 450 km depth) exhibits a reversed pattern. A tentative interpretation of these deeper anomalies is made, as being due mainly to topography of the major upper mantle discontinuities, near 400 and 650 km depth.     

Evolution of the intracratonic Officer Basin, central Australia: implications from subsidence analysis and gravity modelling

ABSTRACT The intracratonic basins of central Australia are distinguished by their large negative Bouguer gravity anomalies, despite the absence of any significant topography. Over the Neoproterozoic to Palaeozoic Officer Basin, the anomalies attain a peak negative amplitude in excess of 150 mGal, amongst the largest of continental anomalies observed on Earth. Using well data from the Officer and Amadeus basins and a data grid of sedimentary thicknesses from the eastern Officer Basin, we have assessed the evolution of these intracratonic basins. One-dimensional backstripping analysis reveals that Officer and Amadeus basin tectonic subsidence was not entirely synchronous. This implies that the basins evolved as discrete geological features once the Centralian Superbasin was dismembered into its constituent basins. Two- and three-dimensional backstripping and gravity modelling suggest that the eastern Officer Basin evolved from a broad continental sag into a region of intracratonic flexural subsidence from the latest Neoproterozoic, when flexure of the lithosphere deepened the northern basin. The results from gravity modelling improve when the crust is thickened beneath the northern margin of the basin and thinned at the southern margin, as has been suggested by recent deep seismic data. The crustal thickening beneath the basin's northern margin abuts the region of greatest topographic relief and is consistent with the observed structure at the edges of many orogenic belts. If the Officer Basin evolved as a foreland-type basin from the late Proterozoic and has retained those features to the present, then one implication is that in the absence of any significant topography, cratonic lithosphere must be able to support stresses over very long periods of geological time.     

Upper mantle heterogeneities beneath Central Europe

Summary. An inversion of ISC travel-time data from selected earthquakes in the distance range 30°-90° to 53 stations in Central Europe has been used to model velocity down to 600 km depth. The model explains 0.1–0.2s of the residuals, as for other array studies, leaving 0.5 s unexplained as noise. The uppermost 100 km of the mantle and crust contains inhomogeneities that correlate remarkably well with the geology. This may be due to deep-seated thermal anomalies or, in some areas, to delays introduced by passage of the rays through sedimentary cover. The deeper anomalies are smaller and unrelated to those in the lithosphere, which suggests that the asthenosphere is decoupled from the rigid lithosphere. The structure at 600 km depth is again quite inhomogeneous and might be due to undulations of the 650 km discontinuity. The models show some suggestion of a high velocity slab trending from east to west beneath the Alps.     

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.     

Global mapping of upper mantle reflectors from long-period SS precursors

Long-period precursors to SS resulting from underside reflections off upper mantle discontinuities ( SdS where d is the discontinuity depth) can be used to map the global distribution and depth of these reflectors. We analyse 5,884 long-period seismograms from the Global Digital Seismograph Network (1976-1987, shallow sources, transverse component) in order to identify SdS arrivals. Corrections for velocity dispersion, topography and crustal thickness at the SS bounce point, and lateral variation in mantle velocity are critical for obtaining accurate estimates of discontinuity depths. The 410 and 660 km discontinuities are observed at average depths of 413 and 653 km, and exhibit large-scale coherent patterns of topography with depth variations up to 40 km. These patterns are roughly correlated with recent tomographic models, with fast anomalies in the transition zone associated with highs in the 410 km discontinuity and lows in the 660 km discontinuity, a result consistent with laboratory measurements of Clapeyron slopes for the appropriate phase changes. The best resolved feature in these maps is a trough in the 660 km discontinuity in the northwest Pacific, which appears to be associated with the subduction zones in this region. Amplitude variations in SdS arrivals are not correlated with discontinuity depths and probably result from focusing and defocusing effects along the ray paths. The SdS arrivals suggest the presence of regional reflectors in the upper mantle above 400 km. However, only the strongest of these features are above probable noise levels due to sampling inadequacies.     

Lithospheric structure of the continent–continent collision zone: eastern Turkey

We infer the lithospheric structure in eastern Turkey using teleseismic and regional events recorded by 29 broad-band stations from the Eastern Turkey Seismic Experiment (ETSE). We combine the surface wave group velocities (Rayleigh and Love) with telesesimic receiver functions to jointly invert for the S -wave velocity structure, Moho depth and mantle-lid (lithospheric mantle) thickness. We also estimated the transverse anisotropy due to Love and Rayleigh velocity discrepancies. We found anomalously low shear wave velocities underneath the Anatolian Plateau. Average crustal thickness is 36 km in the Arabian Plate, 44 km in Anatolian Block and 48 km in the Anatolian Plateau. We observe very low shear wave velocities at the crustal portion (30–38 km) of the northeastern part of the Anatolian Plateau. The lithospheric mantle thickness is either not thick enough to resolve it or it is completely removed underneath the Anatolian Plateau. The shear velocities and anisotropy down to 100 km depth suggest that the average lithosphere–asthenosphere boundary in the Arabian Plate is about 90 and 70 km in Anatolian block. Adding the surface waves to the receiver functions is necessary to constrain the trade-off between velocity and the thickness. We find slower velocities than with the receiver function data alone. The study reveals three different lithospheric structures in eastern Turkey: the Anatolian plateau (east of Karliova Triple Junction), the Anatolian block and the northernmost portion of the Arabian plate. The boundary of lithospheric structure differences coincides with the major tectonic boundaries.     

Modelling rates and distribution of subsidence due to dynamic topography over subducting slabs: is it possible to identify dynamic topography from ancient strata?

Dynamic topography formed over subducting oceanic lithosphere has been proposed as a mechanism to explain certain otherwise anomalous long-wavelength patterns of subsidence inferred from ancient strata. Forward modelling of mantle flow in response to a subducting slab predicts amplitudes and distributions of dynamic topography that may occur due to various subducting slab geometries and histories. Plotting calculated dynamic topographies at a point against time produces tectonic subsidence curves. These subsidence curves show features such as evolution from convex to concave shape, amplitudes up to ~2000 m, subsidence rates up to ~220 m Myr⁻¹, and a general decrease in subsidence amplitude away from the subduction zone, over a distance of ~2000 km. On many convergent continental margins, dynamic topography is likely to be superimposed on other subsidence mechanisms. In back-arc basins, subsidence due to dynamic topography should be distinguishable from that due to extensional tectonics based simply on the temporal subsidence evolution expressed in the subsidence curve shapes. In a foreland basin setting, comparing dynamic topography models with forward models of flexural loading suggest the two processes can generate similar temporal subsidence patterns, but that dynamic topography causes subsidence over significantly greater wavelengths. Matches between calculated subsidence due to dynamic topography and backstripped subsidence patterns from Upper Cretaceous strata of the Western Interior Basin, USA, support the hypothesis that a long-wavelength 'background subsidence' was caused by dynamic topography.     

A technique for computing geoid height anomalies over two-dimensional earth models

Summary A technique is presented for calculating geoid height anomalies over two-dimensional models of Earth structure. The method consists of convolving gravity anomalies over the structure with filters which take into account the finite size of the structure in the third dimension and the curvature of the Earth. Similar filters are also developed for a flat earth case. The method is applied to a sea-surface gravity profile crossing the Tonga-Kermadec trench and is found to give good agreement with a Geos-3 radar altimetry profile in the same region. The example demonstrates that introducing arbitrary offsets in computing gravity anomalies can result in spurious long-wavelength effects in the computed geoid. Comparison of the results obtained using flat earth and spherical earth filters suggests that the effects of the curvature of the Earth only become significant for wavelengths in the gravity field greater than about 1000 km.