Crustal structure of the Whipple Mountains, southeastern California: a refraction study across a region of large continental extension

Summary. Closely spaced refraction profiling across the Whipple Mountains metamorphic core complex in southeastern California yields a complex picture of crustal structure in this region of large continental extension. A NE-directed profile, parallel to the extension direction, reveals a high-velocity mid-crustal layer (6.6–6.8 km s⁻¹) at 16-18 km depth, bounded above and below by laterally discontinuous low-velocity zones (<6.0 km s⁻¹). In marked contrast, a NW-directed profile shows a more uniform 6.0 km s⁻¹ crust down to the crust-mantle boundary. The apparent contrast between these two perpendicular profiles may be related not only to a more complex geologic structure in the NW-SE direction, but also to velocity anisotropy associated with mid-crustal mylonites. Despite the differences between the two refraction profiles, both define a flat Moho at 26-27 km depth with an associated upper mantle-velocity of 7.8 km s⁻¹. This observation is significant as it suggests that, although the amount of extension has been highly variable regionally, the crust is no thinner beneath the Whipple Mountains (where extension has been extreme) than the surrounding mountain ranges. Such an observation requires either that the crust was considerably thicker prior to extension, or that lateral flow in the lower crust and/or inflation of the crust via magmatism occurred contemporaneous with extension.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

The crustal structure of the eastern Marmara region, Turkey by teleseismic receiver functions

We study the crustal structure of eastern Marmara region by applying the receiver function method to the data obtained from the 11 broad-band stations that have been in operation since the 1999 İzmit earthquake. The stacked single-event receiver functions were modelled by an inversion algorithm based on a five-layered crustal velocity model to reveal the first-order shear-velocity discontinuities with a minimum degree of trade-off. We observe crustal thickening from west (29–32 km) to east (34–35 km) along the North Anatolian Fault Zone (NAFZ), but we observe no obvious crustal thickness variation from north to south while crossing the NAFZ. The crust is thinnest beneath station TER (29 km), located near the Black Sea coast in the west and thickest beneath station TAR (35 km), located inland in the southeast. The average crustal thickness and S -wave velocity for the whole regions are  31 ± 2  km and  3.64 ± 0.15 km s⁻¹  , respectively. The eastern Marmara region with its average crustal thickness, high heat flow value (101 ± 11 mW m⁻²) and with its remarkable extensional features seems to have a Basin and Range type characteristics, but the higher average shear velocities (∼3.64 km s⁻¹) and crustal thickening from 29 to 35 km towards the easternmost stations indicate that the crustal structure shows a transitional tectonic regime. Therefore, we conclude that the eastern Marmara region seems to be a transition zone between the Marmara Sea extensional domain and the continental Anatolian inland region.     

The Hatton Bank continental margin—I. Shallow structure from two-ship expanding spread seismic profiles

Summary. The Hatton Bank passive continental margin exhibits thick seaward dipping reflector sequences which consist of basalts extruded during rifting between Greenland and Rockall Plateau. Multichannel seismic reflection profiling across the margin reveals three reflector wedges with a maximum thickness near 7 km, extending from beneath the upper continental slope to the deep ocean basin. We present results of the velocity structure within the dipping reflector sequences at eight locations across the margin, interpreted by synthetic seismogram modelling a set of multichannel expanding spread profiles parallel to the margin. At the top of some reflector sequences, we observe a series of 100 m thick high- and low-velocity zones, which are interpreted as basalt flows alternating with sediments or weathered and rubble layers. At the profile locations, the base of the dipping reflectors correlates with P -wave velocities near 6.5 km s⁻¹. However, elsewhere the reflectors appear to extend significantly deeper than the inferred 6.5 km s⁻¹ velocity contour, indicating that the velocity structure may not be controlled solely by lithological boundaries but also by metamorphic effects. Shear-waves were observed on two lines, permitting the calculation of Poisson's ratio. The decrease in Poisson's ratio from 0.28 to near 0.25 in the upper 5 km of crust may also indicate the effect of metamorphism on seismic properties, or alternatively may be explained by crack closure under load.     

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.     

Crustal structure in the Western Somali Basin

Summary. As part of integrated marine geophysical studies in the Western Somali Basin, we performed 118 sonobuoy experiments to define better the crustal structure of the margins and basin created by the separation of Madagascar and Africa. After using T ²/ X ², conventional slope-intercept methods, and slant-stacked t-p techniques to analyse the data, we combined our solutions with all previous velocity information for the area. Velocity functions were derived for the sediment coiumn, and we detected a high-velocity (4.58 ± 0.29 km s^–1) sediment layer overlying acoustic basement. We confirmed that the crust is indeed seismically oceanic, and that it may be considered either in terms of a layered model – layers 2B (5.42 ± 0.19 km s^–1), 2C (6.23 ± 0.22 km s^–1), 3 (7.03 ± 0.25 km s^–1), and mantle (7.85 ± 0.32 km s^–1) were identified – or a more complex gradient model in which layer 2 is marked by a steeper velocity gradient than underlying layer 3. Integrated igneous crustal thicknesses (1.62 ± 0.22 s, 5.22 ± 0.64 km) are significantly less than what is considered normal. We present a revised seismic transect across the East African margin, as well as total sediment thickness, depth to basement and crustal thickness maps.     

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.     

Seismic and electromagnetic evidence of dehydration as a free water source in the reactivated crust

According to recent estimates, the continental mid-crust contains 35–40 per cent amphibolites. Heating of the crust by an underlying mantle plume, for example beneath continental rifts, high plateaus, and areas of intraplate volcanic activity, releases water. Dehydration of amphibole-bearing rocks at depths of 20–40  km occurs mainly in the temperature range 650–700 °C, and this releases about 0.4  wt per cent of water.
  Seismic tomography studies of the crust in the Kirgyz Tien Shan Range, where the age of the tectonic activity is less than 30  Ma, revealed a low-velocity zone in the mid-crust. The velocity of P waves was 0.4  km  s⁻¹ lower than in normal crust. MT sounding data in the region show the existence of a low-resistivity layer with an average resistivity of about 25  Ω  m at the depth of the low-velocity layer. The spatial correlation of the observed anomalous layers and calculated effect of fluid phase on seismic and electric parameters of rocks suggests the presence of aqueous fluids released by the heating of the mid-crust.     

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⁻¹.     

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.     

Seismic refraction profiles between Cyprus and Israel and their interpretation

Summary. A long seismic refraction profile was carried out between southern Israel and Cyprus. The seismic energy was generated by 33 sea shots each of 0.8 t explosives and was recorded by land stations in Israel and Cyprus and by ocean bottom seismographs deployed along the profile.
The results showed that the continental crust of southern Israel thins towards the Mediterranean underneath a northward thickening sedimentary cover. Cyprus is underlain by a 35 km thick continental crust thinning south-wards and extending to Mt Eratosthenes. Between Mt Eratosthenes and the Israel continental shelf the crystalline crust is composed of high velocity (6.5 km s^-1)material and is about 8 km thick. It is covered by 12–14 km of sediments and may represent a fossil oceanic crust.     

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.     

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.     

Transitional continental–oceanic structure beneath the Norwegian Sea from inversion of surface wave group velocity data

We have analysed the fundamental mode of Love and Rayleigh waves generated by 12 earthquakes located in the mid-Atlantic ridge and Jan Mayen fracture zone. Using the multiple filter analysis technique, we isolated the Rayleigh and Love wave group velocities for periods between 10 and 50  s. The surface wave propagation paths were divided into five groups, and average group velocities calculated for each group. The average group velocities were inverted and produced shear wave velocity models that correspond to a quasi-continental oceanic structure in the Greenland–Norwegian Sea region. Although resolution is poor at shallow depth, we obtained crustal thickness values of about 18  km in the Norwegian Sea area and 9  km in the region between Svalbard and Iceland. The abnormally thick crust in the Norwegian Sea area is ascribed to magmatic underplating and the thermal blanketing effect of sedimentary layers. Maximum crustal shear velocities vary between 3.5 and 3.9  km  s⁻¹ for most paths. An average lithospheric thickness of 60  km was observed, which is lower than expected for oceanic-type structure of similar age. We also observed low shear wave velocities in the lower crust and upper mantle. We suggest that high heat flow extending to depths of about 30  km beneath the surface can account for the thin lithosphere and observed low velocities. Anisotropy coefficients of 1–5 per cent in the shallow layers and >7 per cent in the upper mantle point to the existence of polarization anisotropy in the region.