Crustal model of the Bransfield Rift, West Antarctica, from detailed OBS refraction experiments

The first detailed deep seismic refraction study in the Bransfield Strait, West Antarctica, using sensitive OBSs (ocean bottom seismographs) was carried out successfully during the Antarctic summer of 1990/1991. The experiment focused on the deep crustal structure beneath the axis of the Bransfield Rift. Seismic profile DSS-20 was located exactly in the Bransfield Trough, which is suspected to be a young rift system. Along the profile, five OBSs were deployed at spacings of 50-70 km. 51 shots were fired along the 310 km profile. This paper gives the first presentation of the results. A detailed model of the crustal structure was obtained by modelling the observed traveltimes and amplitudes using a 2-D ray-tracing technique. The uppermost (sedimentary?) cover, with velocities of 2.0-5.5 km s⁻¹, reaches a depth of up to 8 km. Below this, a complex with velocities of 6.4-6.8 km s⁻¹ is observed. The presence of a high-velocity body, with V _p= 7.3-7.7 km s⁻¹, was detected in the 14-32 km depth range in the central part of the profile. These inhomogeneities can be interpreted as a stage of back-arc spreading and stretching of the continental crust, coinciding with the Deception-Bridgeman volcanic line. Velocities of 8.1 km s⁻¹, characteristic of the Moho, are observed along the profile at a depth of 30-32 km.     

Seismic anisotropy within the uppermost mantle of southern Germany

This paper presents an updated interpretation of seismic anisotropy within the uppermost mantle of southern Germany. The dense network of reversed and crossing refraction profiles in this area made it possible to observe almost 900 traveltimes of the P_n phase that could be effectively used in a time-term analysis to determine horizontal velocity distribution immediately below the Moho. For 12 crossing profiles, amplitude ratios of the P_n phase compared to the dominant crustal phase were utilized to resolve azimuthally dependent velocity gradients with depth. A P -wave anisotropy of 3–4 per cent in a horizontal plane immediately below the Moho at a depth of 30 km, increasing to 11 per cent at a depth of 40 km, was determined. For the axis of the highest velocity of about 8.03 km s⁻¹ at a depth of 30 km a direction of N31°F was obtained. The azimuthal dependence of the observed P_n amplitude is explained by an azimuth-dependent sub-Moho velocity gradient decreasing from 0.06 s⁻¹ in the fast direction to 0 s⁻¹ in the slow direction of horizontal P -wave velocity. From the seismic results in this study a petrological model suggesting a change of modal composition and percentage of oriented olivine with depth was derived.     

Crustal and upper mantle structure of the northwestern North Island, New Zealand, from seismic refraction data

The crustal and upper mantle structure of the northwestern North Island of New Zealand is derived from the results of a seismic refraction experiment; shots were fired at the ends and middle of a 575 km-long line extending from Lake Taupo to Cape Reinga. The principal finding from the experiment is that the crust is 25 ± 2 km thick, and is underlain by what is interpreted to be an upper mantle of seismic velocity 7.6 ± 0.1 km s⁻¹, that increases to 7.9 km s⁻¹ at a depth of about 45 km. Crustal seismic velocities vary between 5.3 and 6.36 km s⁻¹ with an average value of 6.04 km s⁻¹. There are close geophysical and geological similarities between the north-western North Island of New Zealand and the Basin and Range province of the western United States. In particular, the conditions of low upper-mantle seismic velocities, thin crust with respect to surface elevation, and high heat-flow (70–100 mW m⁻²) observed in these two areas can be ascribed to their respective positions behind an active convergent margin for about the past 20 Myr.     

Upper crustal shear velocity models from higher mode Rayleigh wave dispersion in Scotland

Summary. Group velocities for first and second higher mode Rayleigh waves, in the frequency range 0.8–4.8 Hz, generated from a local earthquake of magnitude 3.7 M _L in western Scotland, are measured at stations along the 1974 LISPB line. These provide detailed information about the crustal structure west of the line. The data divide the region into seven apparently homogeneous provinces. Averaged higher mode velocity dispersion curves for each province are analysed simultaneously using a linearized inversion technique, yielding regionalized shear velocity profiles down to a depth of 17 km into the upper crust. Shear wave velocity is between 3.0 and 3.4 km s⁻¹ in the upper 2 km, with a slow increase to around 3.8 km s⁻¹. P -wave models computed using these results agree with profiles from the LISPB and LUST refraction experiments.     

Crust and upper mantle structure of the central Iberian Meseta (Spain)

Summary. Quarry blasts recorded along three lines on the central Iberian Meseta are used in an attempt to interpret the crustal structure. The results of the interpretation of the data, together with published surface wave and earthquake data, suggest a layered structure of the crust having the following features: the basement, in some areas covered by up to 4 km of sediments, has a P -velocity of 6.1 km s⁻¹; a low-velocity layer, between 7 and 11 km depth, seems to exist on the basis of both P and S interpretation of seismic data; a thick middle crust of 12 km has a P -velocity of 6.4 km s⁻¹ and overlies a lower crust with a mean P -velocity of 6.9 km s⁻¹ and a possible slight negative gradient; the mean v _p/ v _s ratio for the crust is about 1.75; the Moho is reached at about 31 km depth and consists of a transition zone at least 1.5 km thick. The P -velocity of the upper mantle is close to 8.1 km s⁻¹ and the S -velocity about 4.5 km s⁻¹, which gives a v _p /v _s ratio of 1.8 for the uppermost mantle. A tentative petrological interpretation of the velocities and composition of the layers is given.     

Velocity structure of upper-mantle transition zones beneath central eurasia from seismic inversion using genetic algorithms

We present velocity constraints for the upper-mantle transition zones beneath Central Siberia based on observations of the 1982 RIFT Deep Seismic Sounding (DSS) profile. The data consist of seismic recordings of a nuclear explosion in north-western Siberia along a 2600 km long seismic profile extending from the Yamal Peninsula to Lake Baikal. We invert seismic data from the mantle transition zones using a non-linear inversion scheme using a genetic algorithm for optimization and the WKBJ method to compute the synthetic seismograms. A statistical error analysis using a graph-binning technique was performed to provide uncertainty values in the velocity models.
Our best model for the upper-mantle velocity discontinuity near 410 km depth has a two-stage velocity-gradient structure, with velocities increasing from 8.70–9.25 km s⁻¹ over a depth range of 400–415 km, a gradient of 0.0433 s⁻¹, and from 9.25–9.60 km s⁻¹ over a depth range of 415–435 km, a gradient of 0.0175 s⁻¹. This derived model is consistent with other seismological observations and mineral-physics models. The model for the velocity discontinuity near 660 km depth is simple, sharp and includes velocities increasing from 10.15 km s⁻¹ at 655 km depth to 10.70 km s⁻¹ at 660 km depth, a gradient of 0.055 s⁻¹.     

Axial Intrusion Zone beneath the Median Valley of the Mid-Atlantic Ridge at 37°N Detected by Explosion Seismology

Summary. Four seismic refraction lines, three of which had shots every 250 m, were shot across, along and parallel to the median valley of the Mid-Atlantic Ridge at 37° N. A method has been developed for calculating the effect on the travel times of the rough sea-floor relief beneath the profiles and has been used to correct all the travel times for this effect. Most arrivals were from a main refractor of apparent velocity 5·4 to 6·3 km s⁻¹; only beyond 35 km were faster arrivals observed from an 8·09 ± 36 km s⁻¹ refractor. The main refractor corresponds in depth, at least approximately, to the top of Layer 3 of the ocean basins but its velocity is significantly less than normal for Layer 3, perhaps due to dip. A study of time residuals along two profiles across the median valley indicates the presence of a 2 to 3 km wide low velocity zone (about 3·2 km s⁻¹) beneath the median valley floor. This zone extends over the upper 2·5 km of the crust and is believed to represent a zone of intrusion through which magma passes on its way to the sea floor.     

A seismic study of deep geological structure in the Bristol Channel area, SW Britain

Summary. Results from eight seismic refraction lines, 35–90 km long, in the Bristol Channel area are presented. The data, mostly land recordings of marine shots, have been interpreted mainly by ray-tracing and time-term modelling. Upper layer velocities through Palaeozoic rocks usually fall within the range 4.8–5.2 km s⁻¹. Below the Carboniferous Limestone with a normal velocity of 5.1–5.2 kms⁻¹, the Old Red Sandstone with a velocity of 4.7–4.8 kms⁻¹ acts as a low velocity layer, as do parts of the underlying Lower Palaeozoic succession. In the central South Wales/Bristol Channel area and the Mendips, a 5.4–5.5 km s⁻¹ refractor is correlated with a horizon at or near the top of the Lower Palaeozoic succession. Under the whole area, except for north Devon, a 6.0–6.2 km s⁻¹ basal refractor has been located and is correlated with Precambrian crystalline basement rocks. In general, this refractor deepens southwards from a series of basement highs, which existed before the major movements of the Variscan orogeny in South Wales, resulting in a southerly thickening of the pre Upper Carboniferous supra-basement sequence. In north Devon, a 6.2 km s⁻¹ refractor at shallow depth, interpreted as a horizon in the Devonian or Lower Palaeozoic succession, overlies a deep reflector that may represent the Precambrian crystalline basement.     

Detailed S-wave structure in the Dublin Basin and its northern margin

Summary. The seismic structure has been measured to a depth of about 3 km along a 30 km seismic profile in east central Ireland. This profile is unusual in that it is the S -wave velocity—depth structure that has been measured to a degree of precision more normally associated with P -wave results. One reason for this is that the sources used were quarry blasts which generated strong S -waves and short-period surface waves but rather weak P -waves.
The results show a layer of Carboniferous limestone with shear velocity 2.65 km⁻¹ s overlying a layer with a velocity of 3.06 km s⁻¹. This second layer was interpreted as Lower Palaeozoic strata (Silurian/Ordovician) since this velocity was evident in an inlier seen at the surface at the northern end of the line. A third refraction horizon, shear velocity 3.45 km s⁻¹ and displaying a basinal structure, was also recognized. This may be Cambrian or Precambrian basement.     

Crustal structure of the Mid-Atlantic ridge crest at 37° N

Summary. A structural model of the Mid-Atlantic Ridge at 37° N is proposed on the basis of travel-time data and synthetic seismograms. At the ridge axis the crust is only 3 km thick and overlies material with an anomalously low'upper mantle'velocity of 7.2 km s⁻¹. Crustal thickening and the formation of layer 3 and a layer with velocity 7.2–7.3 km s⁻¹ takes place within a few kilometres of the axis, producing a 6–7 km thick crust by less than 10 km from the axis. A normal upper mantle velocity of 8.1 km s⁻¹ exists within 10 km of the axis. Shear waves propagate across the axis, thus precluding the existence of any sizeable magma chamber at shallow depth.     

Ivrea seismic array: a study of continental crust and upper mantle

A seismic-array study of the continental crust and upper mantle in the Ivrea-Yerbano and Strona-Ceneri zones (northwestern Italy) is presented. A short-period network is used to define crustal P - and S -wave velocity models from earthquakes. The analysis of the seismic-refraction profile LOND of the CROP-ECORS project provided independent information and control on the array-data interpretation.
Apparent-velocity measurements from both local and regional earthquakes, and time-term analysis are used to estimate the velocity in the lower crust and in the upper mantle. The geometry of the upper-lower crust and Moho boundaries is determined from the station delay times.
We have obtained a three-layer crustal seismic model. The P -wave velocity in the upper crust, lower crust and upper mantle is 6.1±0.2 km s⁻¹, 6.5±0.3 km s⁻¹ and 7.8±0.3 km s⁻¹ respectively. Pronounced low-velocity zones in the upper and lower crust are not observed. A clear change in the velocity structure between the upper and lower crust is documented, constraining the petrological interpretation of the Ivrea-type reflective lower continental crust derived from small-scale petrophysical data. Moreover, we found a V _P/ V _S ratio of 1.69±0.04 for the upper crust and 1.82±0.08 for the lower crust and upper mantle. This is consistent with the structural and petrophysical differences between a compositionally uniform and seismically transparent upper crust and a layered and reflective lower crust. The thickness of the lower crust ranges from about 8 km in front of the Ivrea body (ARVO, Arvonio station) in the northern part of the array to a maximum of about 15 km in the southern part of the array. The lower crust reaches a minimum depth of 5 km below the PROV (Provola) station.     

The seismic structure of thinned continental crust in the northern Bay of Biscay

Summary. The stretching and thinning of the continental crust, which occurs during the formation of passive continental margins, may cause important changes in the velocity structure of such crust. Further, crust attenuated to a few kilometres' thickness, can be found underlying 'oceanic' water depths. This paper poses the question of whether thinned continental crust can be distinguished seismically from normal oceanic crust of about the same thickness. A single seismic refraction line shot over thinned continental crust as part of the North Biscay margin transect in 1979 was studied in detail. Tau— p inversion suggested that there are differences between oceanic and continental crust in the lower crustal structure. This was confirmed when synthetic seismograms were calculated. The thinned continental crust (β± 7.0) exhibits a two-gradient structure in the non-sedimentary crust with velocities between 5.9 and 7.4 km s⁻¹; an upper 0.8 s⁻¹ layer overlies a 0.4 s⁻¹ layer. No layer comparable to oceanic layer 3 was detected. The uppermost mantle also contains a low-velocity zone.     

Deep seismic soundings along Hirapur-Mandla profile, central India

Summary. The crustal depth section along Hirapur-Mandla profile has been computed in two steps from Deep Seismic Sounding (DSS) data. The shallow section up to the crystalline basement is derived by inverting first arrival refraction travel times. The upper Vindhyan sediments (velocity 4.5 km s⁻¹) have a maximum thickness of about 1.5 km at Bakshaho. The lower Vindhyan sediments (velocity 5.4 km s⁻¹) were deposited north of Narmada-Son lineament between Katangi and Narsinghgarh in a graben developed in crystalline basement. The thickness of the lower Vindhyans increases from north to south towards Katangi and the depth to the basement reaches 5.5 km near Jabera. The depth to the Moho boundary varies from 39.5 km near Tikaria to 45 km at Narsinghgarh. The narrow block between Katangi and Jabalpur forms a horst feature which represents the Narmada-Son lineament forming the southern boundary of the Vindhyan basin. Two-dimensional ray tracing was performed generating travel time curves from various shot points which were matched with observed travel time data.     

An anomalous high-velocity layer at shallow crustal depths in the Narmada zone, India

The Narmada zone in central India is a zone of weakness that separates the region of Vindhyan (Meso-Neoproterozoic) deposition to the north from Gondwana (Permo-Carboniferous–lower Cretaceous) deposits to the south. The reinterpretation of analogue seismic refraction data, acquired during the early 1980s, using 2-D ray-tracing techniques reveals a basement (velocity 5.8–6.0 km s⁻¹ ) topography suggesting that the Narmada zone, bounded by the Narmada North and Narmada South faults is a region of basement uplift. A layer of anomalously high velocity (6.5–6.7 km s⁻¹ ) at depths between 1.5 and 9.0 km appears to be present in the entire region. Within the Narmada zone this layer occurs at shallower depths than outside the Narmada zone. At two places within the Narmada zone this layer is at a depth of about 1.5 km. This layer cannot be considered as the top of the lower crust because in this case it should have produced large positive gravity anomalies at the shallowest parts. Instead, these parts correspond to Bouguer gravity lows. Furthermore, lower crust at such shallow depths has not been reported from any other part of the Indian shield. Therefore, this layer is likely to represent the top of a high-velocity mafic body that has different thicknesses in different places.     

The Agulhas Plateau: structure and evolution of a Large Igneous Province

Large Igneous Provinces (LIP) are of great interest due to their role in crustal generation, magmatic processes and environmental impact. The Agulhas Plateau in the southwest Indian Ocean off South Africa has played a controversial role in this discussion due to unclear evidence for its continental or oceanic crustal affinity. With new geophysical data from seismic refraction and reflection profiling, we are able to present improved evidence for its crustal structure and composition. The velocity–depth model reveals a mean crustal thickness of 20 km with a maximum of 24 km, where three major units can be identified in the crust. In our seismic reflection records, evidence for volcanic flows on the Agulhas Plateau can be observed. The middle crust is thickened by magmatic intrusions. The up to 10 km thick lower crustal body is characterized by high seismic velocities of 7.0–7.6 km s⁻¹. The velocity–depth distribution suggests that the plateau consists of overthickened oceanic crust similar to other oceanic LIPs such as the Ontong-Java Plateau or the northern Kerguelen Plateau. The total volume of the Agulhas Plateau was estimated to be 4 × 10⁶ km³ of which about 10 per cent consists of extruded igneous material. We use this information to obtain a first estimate on carbon dioxide and sulphur dioxide emission caused by degassing from this material. The Agulhas Plateau was formed as part of a larger LIP consisting of the Agulhas Plateau itself, Northeast Georgia Rise and Maud Rise. The formation time of this LIP can be estimated between 100 and 94 (± 5) Ma.     

Shallow velocity structure along the Hirapur–Mandla profile using traveltime inversion of wide-angle seismic data, and its tectonic implications

In order to investigate the velocity structure, and hence shed light on the related tectonics, across the Narmada–Son lineament, traveltimes of wide-angle seismic data along the 240 km long Hirapur–Mandla profile in central India have been inverted. A blocky, laterally heterogeneous, three-layer velocity model down to a depth of 10 km has been derived. The first layer shows a maximum thickness of the upper Vindhyans (4.5 km s⁻¹ ) of about 1.35 km and rests on top of normal crystalline basement, represented by the 5.9 km s⁻¹ velocity layer. The anomalous feature of the study is the absence of normal granitic basement in the great Vindhyan Graben, where lower Vindhyan sediments (5.3 km s⁻¹ ) were deposited during the Precambrian on high-velocity (6.3 km s⁻¹ ) metamorphic rock. The block beneath the Narmada–Son lineament represents a horst feature in which high-velocity (6.5 km s⁻¹ ) lower crustal material has risen to a depth of less than 2 km. South of the lineament, the Deccan Traps were deposited on normal basement during the upper Cretaceous period and attained a maximum thickness of about 800 m.     

Crust and upper mantle shear velocities from controlled sources

b
A two ship refraction profile was undertaken on the Australian continental shelf during the Banda Sea geophysical program, carried out by the Woods Hole Oceanographic Institution, the Scripps Institution of Oceanography and the Geological Survey of Indonesia. S waves originating close to the sea bottom were observed to distances of up to 1150 km at an array of stations in northern Australia.
These observations are interpreted as implying S mantle velocities of 4.60 km s^-1 from a depth of 45 km to a depth of 76 km and 4.72 km s^-1 below a depth of 76 km.
Ratios of the P and S travel times (V_p/V_s) have been determined to be 1.74 in the crust rising to a value of greater than 1.79 below a velocity discontinuity at a depth of 200 km. It is inferred that this high value arises because the effect of temperature is greater for S than for P .
Using the data from this and other studies in the shield region of Northern Australia it has been found that the _S travel times are significantly less than predicted by the Jeffreys—Bullen tables.     

Seismic moment release during slab rupture beneath the Banda Sea

The highest intermediate depth moment release rates in Indonesia occur in the slab beneath the largely submerged segment of the Banda arc in the Banda Sea to the east of Roma, termed the Damar Zone. The most active, western-part of this zone is characterized by downdip extension, with moment release rates (∼10¹⁸ Nm yr^–1 per 50 km strike length) implying the slab is stretching at ∼10⁻¹⁴ s⁻¹ consistent with near complete slab decoupling across the 100–200 km depth range. Differential vertical stretching along the length of the Damar Zone is consistent with a slab rupture front at ∼100–200 km depth beneath Roma propagating eastwards at ∼100 km Myr^–1. Complexities in the slab deformation field are revealed by a narrow zone of anomalous in-plane P -axis trends beneath Damar, where subhorizontal constriction suggests extreme stress concentrations ∼100 km ahead of the slab rupture front. Such stress concentrations may explain the anomalously deep ocean gateways in this region, in which case ongoing slab rupture may have played a key role in modulating the Indonesian throughflow in the Banda Sea over the last few million years.     

Crustal structure of the Faeroe–Shetiand Channel

Summary. The deep structure of the Faeroe–Shetland Channel has been investigated as part of the North Atlantic Seismic Project. Shot lines were fired along and across the axis of the Channel, with recording stations both at sea and on adjacent land areas. At 61°N, 1.7 km of Tertiary sediments overlies a 3.9–4.5 km s^-1 basement interpreted as the top of early Tertiary volcanics. A main 6.0–6.6 km s^-1 crustal refractor interpreted as old oceanic crust occurs at about 9 km depth. The Moho (8.0 ° 0.2 km s^-1) is at about 15–17 km depth. There is evidence that P _n may be anisotropic beneath the Faeroe–Shetland Channel. Arrivals recorded at land stations show characteristics best explained by scattering at an intervening boundary which may be the continent–ocean crustal contact or the edge of the volcanics.
The Moho delay times at the shot points, determined by time-term analysis, show considerable variation along the axis of the Channel. They correlate with the basement topography, and the greatest delays occur over the buried extension of the Faeroe Ridge at about 60° 15'N, where they are nearly 1 s more than the delays at 61°N after correction for the sediments. The large delays are attributed to thickening of the early Tertiary volcanic layer with isostatic downsagging of the underlying crust and uppermost mantle in response to the load, rather than to thickening of the main crustal ayer.
The new evidence is consistent with deeply buried oceanic crust beneath the Faeroe–Shetland Channel, forming a northern extension of Rockall Trough. The seabed morphology has been grossly modified by the thick and laterally variable pile of early Tertiary volcanic rocks which swamped the region, accounting for the anomalous shallow bathymetry, the transverse ridges and the present narrowness of the Channel.     

Forward modelling receiver functions for crustal structure beneath station TBZ (Trabzon, Turkey)

We use teleseismic three-component digital data from the Trabzon, Turkey broadband seismic station TBZ to model the crustal structure by the receiver function method. The station is located at a structural transition from continental northeastern Anatolia to the oceanic Black Sea basin. Rocks in the region are of volcanic origin covered by young sediments. By forward modelling the radial receiver functions, we construct 1-D crustal shear velocity models that include a lower crustal low-velocity zone, indicating a partial melt mechanism which may be the source of surfacing magmatic rocks and regional volcanism. Within the top 5 km, velocities increase sharply from about 1.5 to 3.5 km s⁻¹. Such near-surface low velocities are caused by sedimentation, extending from the Black Sea basin. Velocities at around 20 km depth have mantle-like values (about 4.25 km s⁻¹ ), which easily correlate to magmatic rocks cropping out on the surface. At 25 km depth there is a thin low-velocity layer of about 4.0 km s⁻¹. The average Moho velocity is about 4.6 km s⁻¹, and its depth changes from 32 to 40 km. Arrivals on the tangential components indicate that the Moho discontinuity dips approximately southwards, in agreement with the crustal thickening to the south. We searched for the solution of receiver functions around the regional surface wave group velocity inversion results, which helped alleviate the multiple solution problem frequently encountered in receiver function modelling.
Station TBZ is a recently deployed broadband seismic station, and the aim of this study is to report on the analysis of new receiver function data. The analysis of new data in such a structurally complex region provides constraining starting models for future structural studies in the region.