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.     

Finite-frequency traveltime tomography of San Francisco Bay region crustal velocity structure

Seismic velocity structure of the San Francisco Bay region crust is derived using measurements of finite-frequency traveltimes. A total of 57 801 relative traveltimes are measured by cross-correlation over the frequency range 0.5–1.5 Hz. From these are derived 4862 'summary' traveltimes, which are used to derive 3-D P -wave velocity structure over a 341 × 140 km² area from the surface to 25 km depth. The seismic tomography is based on sensitivity kernels calculated on a spherically symmetric reference model. Robust elements of the derived P -wave velocity structure are: a pronounced velocity contrast across the San Andreas fault in the south Bay region (west side faster); a moderate velocity contrast across the Hayward fault (west side faster); moderately low velocity crust around the Quien Sabe volcanic field and the Sacramento River delta; very low velocity crust around Lake Berryessa. These features are generally explicable with surface rock types being extrapolated to depth ∼10 km in the upper crust. Generally high mid-lower crust velocity and high inferred Poisson's ratio suggest a mafic lower crust.     

Determination of the three-dimensional seismic structure of the crust and upper mantle in the Central Midlands of England

Summary. Teleseismic P -wave residuals relative to CWF, a permanent shortperiod seismic station on Charnwood Forest in the Central Midlands of England, have been determined for two small aperture arrays deployed over the Precambrian block of Charnwood and its surrounding Phanerozoic sediments. The data have been inverted to produce a block model of the P -wave velocity variations in the crust and upper mantle beneath the study region. The results are consistent with significant variations penetrating to a depth of at least 50 km. Low velocities are associated with two upper crustal intrusive bodies, the Caledonian Mountsorrel granodiorite and the South Leicestershire diorites. A longer-wavelength variation at lower crustal/upper mantle depths could arise from the Moho dipping to the south-west beneath the study region, and whose strike sub-parallels the dominant Charnian trend of the major basement structures in this part of Central England.     

Crustal and uppermost mantle structure in southern Africa revealed from ambient noise and teleseismic tomography

Rayleigh wave phase velocity maps in southern Africa are obtained at periods from 6 to 40 s using seismic ambient noise tomography applied to data from the Southern Africa Seismic Experiment (SASE) deployed between 1997 and 1999. These phase velocity maps are combined with those from 45 to 143 s period which were determined previously using a two-plane-wave method by Li & Burke. In the period range of overlap (25–40 s), the ambient noise and two-plane-wave methods yield similar phase velocity maps. Dispersion curves from 6 to 143 s period were used to estimate the 3-D shear wave structure of the crust and uppermost mantle on an 1°× 1° grid beneath southern Africa to a depth of about 100 km. Average shear wave velocity in the crust is found to vary from 3.6 km s^–1 at 0–10 km depths to 3.86 km s^–1 from 20 to 40 km, and velocity anomalies in these layers correlate with known tectonic features. Shear wave velocity in the lower crust is on average low in the Kaapvaal and Zimbabwe cratons and higher in the surrounding Proterozoic terranes, such as the Limpopo and the Namaqua-Natal belts, which suggests that the lower crust underlying the Archean cratons is probably less mafic than beneath the Proterozoic terranes. Crustal thickness estimates agree well with a previous receiver function study of Nair et al. . Archean crust is relatively thin and light and underlain by a fast uppermost mantle, whereas the Proterozoic crust is thick and dense with a slower underlying mantle. These observations are consistent with the southern African Archean cratons having been formed by the accretion of island arcs with the convective removal of the dense lower crust, if the foundering process became less vigorous in arc environments during the Proterozoic.     

A crustal seismic velocity model for the UK, Ireland and surrounding seas

A regional model of the 3-D variation in seismic P -wave velocity structure in the crust of NW Europe has been compiled from wide-angle reflection/refraction profiles. Along each 2-D profile a velocity–depth function has been digitised at 5 km intervals. These 1-D velocity functions were mapped into three dimensions using ordinary kriging with weights determined to minimise the difference between digitised and interpolated values. An analysis of variograms of the digitised data suggested a radial isotropic weighting scheme was most appropriate. Horizontal dimensions of the model cells are optimised at 40 × 40 km and the vertical dimension at 1 km. The resulting model provides a higher resolution image of the 3-D variation in seismic velocity structure of the UK, Ireland and surrounding areas than existing models. The construction of the model through kriging allows the uncertainty in the velocity structure to be assessed. This uncertainty indicates the high density of data required to confidently interpolate the crustal velocity structure, and shows that for this region the velocity is poorly constrained for large areas away from the input data.     

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.     

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.     

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.     

Calibration of the Yellowknife Seismic Array with First Zone Explosions

Summary Recordings from a crustal seismic experiment, which was conducted in the Yellowknife area in 1966, were used for calibration of the Yellow-knife seismic array. In the immediate vicinity of the array the crust is found to be very uniform. A superficial layer with an intercept time of 0–172 ± 0–012s and unknown velocity is underlain by a crust with a P wave velocity of 6.04 ± 0–01 km s^-1 near the top: assuming this velocity constant throughout the second layer, the total thickness of the crust is about 34 ± 2 km. The Mohorovicic discontinuity is horizontal under the array within the resolution of this experiment and the apparent P_n velocity is 8.15 km s^-1. At a distance of a few tens of kilometres the crustal uniformity breaks down. The distances are such that, for most teleseismic signals, the effect of these in homogeneities should be negligible.     

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.     

A long-range seismic profile in Southeastern Australia

Summary. Nine portable seismic stations deployed across the Western Plains of New South Wales recorded signals in the distance range 250–1000 km from large timed explosions at both ends of the line. A velocity—depth model derived from the travel-time data has the following features: a two-layer crust with a thickness of 35 km; a sub-Moho velocity of 7.98 km/s; an abrupt increase to 8.36 km/s at 100 km depth; a further step to 8.72 km/s at 190 km depth, with a low-velocity channel immediately above the discontinuity. The model has several features in common with others derived from long-range profiles in Australia and elsewhere. The data, however, provide the first suggestion of a low P -velocity channel in Eastern Australia.     

Pre-drilling reflection survey of the Black Forest, SW Germany

Summary. The unified seismic exploration program, consisting of 345 km of deep reflection profiling, a 200 km refraction profile, an expanding spread profile and near-surface high resolution reflection meaasurements, revealed a strongly differentiated crust beneath the Black Forest. The highly reflective lower crust contains numerous horizontal and dipping reflectors at depths of 13-14 km down to the crust-mantle boundary (Moho). The Moho appears as a flat horizontal first order discontinuity at a relatively shallow level of 25–27 km above a transparent upper mantle. From modelling of synthetic near-vertical and wide-angle seismograms using the reflectivity method the lower crust is supposed to be composed of laminae with an average thickness of about 100 m and velocity differences of greater than 10% increasing from top to bottom. The upper crust is characterised by mostly dipping reflectors, associated with bivergent underthrusting and accretion tectonics of Variscan age and with extensional faults of Mesozoic age. A bright spot at 9.5 km depth is characterised by low velocity material suggesting a fluid trap. It appears on all of the three profiles in the centre of the intersection region. The upper crust seems to be decoupled from the lowest crust by a relatively transparent zone which is' also identified as a low-velocity zone. This low velocity channel is situated directly above the laminated lower crust. The laminae in the Rhinegraben area are displaced vertically to greater depths indicating an origin before Tertiary rift formation and a subsidence of the whole graben wedge.     

Estimation of seasonal changes in the flow of Shirase Glacier using JERS-1/SAR image correlation

This paper presents estimates of detailed seasonal variations in ice-flow velocity for Shirase Glacier calculated using data obtained by Japanese Earth Resources Satellite-1 (JERS-1) synthetic aperture radar (SAR). We used 12 pairs of images (44-day repeat cycle) over the interval from 30 April 1996 to 1 July 1998 to estimate ice-flow fields using an image correlation method. Geometric registration was performed with reference to the RADARSAT Antarctic Mapping Project (RAMP) image dataset. Error analysis based on feature mismatch indicated an absolute error of ±0.30 km/a and relative error of ±0.04 km/a in the estimated flow velocity. The obtained ice-flow velocity increases rapidly from the upstream region (1.18 km/a) to the grounding line, where it becomes stagnant (2.32 km/a), before accelerating gradually to 2.62–2.82 km/a in the downstream region and then increasing to 3.05–3.50 km/a at the terminus of the floating ice tongue. The ice-flow velocities in the downstream region are highly variable, depending on both the distance from the grounding line and the observed epoch (season). Most of the obtained seasonal variations in ice-flow velocity at the floating ice tongue are within the range of the associated error estimate, but the annual difference between 1997 (3.11 km/a) and 1998 (3.50 km/a) is significant, reflecting a possible acceleration in the ice-flow velocity in association with the disappearance of the floating ice tongue between April and May of 1998. In terms of the summer–winter difference in averaged air temperature, the large difference recorded in 1997 (17.0 °C) relative to 1996 (13.9 °C) corresponds to a reduced ice-flow velocity in 1997 (approximately 0.20 km/a) relative to that in 1996 (approximately 0.30 km/a), indicating interactions between air, sea ice, and glacier flow in Lützow-Holm Bay.     

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.     

Imaging the Hikurangi subduction zone, New Zealand, using teleseismic receiver functions: crustal fluids above the forearc mantle wedge

We image the Hikurangi subduction zone using receiver functions derived from teleseismic earthquakes. Migrated receiver functions show a northwest dipping low shear wave feature down to 60 km depth, which we associate with the crust of the subducted Pacific Plate. Receiver functions (RF) at several stations also show a pair of negative and positive polarity phases with associated conversion depths of ∼20–26 km, where the subducted Pacific Plate is at a depth of ∼40–50 km beneath the overlying Australian Plate. RF inversion solutions model these phases with a thin low S -wave velocity zone less than 4 km thick, and an S -wave velocity contrast of more than ∼0.5 km s⁻¹ with the overlying crust. We interpret this phase pair as representing fluids near the base of the lower crust of the Australian Plate, directly overlying the forearc mantle wedge.     

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.     

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.     

Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis – II. Crustal and upper-mantle structure

We determine the 3-D shear wave speed variations in the crust and upper mantle in the southeastern borderland of the Tibetan Plateau, SW China, with data from 25 temporary broad-band stations and one permanent station. Interstation Rayleigh wave (phase velocity) dispersion curves were obtained at periods from 10 to 50 s from empirical Green's function (EGF) derived from (ambient noise) interferometry and from 20 to 150 s from traditional two-station (TS) analysis. Here, we use these measurements to construct phase velocity maps (from 10 to 150 s, using the average interstation dispersion from the EGF and TS methods between 20 and 50 s) and estimate from them (with the Neighbourhood Algorithm) the 3-D wave speed variations and their uncertainty. The crust structure, parametrized in three layers, can be well resolved with a horizontal resolution about of 100 km or less. Because of the possible effect of mechanically weak layers on regional deformation, of particular interest is the existence and geometry of low (shear) velocity layers (LVLs). In some regions prominent LVLs occur in the middle crust, in others they may appear in the lower crust. In some cases the lateral transition of shear wave speed coincides with major fault zones. The spatial variation in strength and depth of crustal LVLs suggests that the 3-D geometry of weak layers is complex and that unhindered crustal flow over large regions may not occur. Consideration of such complexity may be the key to a better understanding of relative block motion and patterns of seismicity.     

Crustal structure beneath exposed accreted terranes of Southern Alaska

Summary. The crustal structure beneath the exposed terranes of southern Alaska has been explored using coincident seismic refraction and reflection profiling. A wide-angle reflector at 8–9 km depth, at the base of an inferred low-velocity zone, underlies the Peninsular and Chugach terranes, appears to truncate their boundary, and may represent a horizontal decollement beneath the terranes. The crust beneath the Chugach terrane is characterized by a series of north-dipping paired layers having low and high velocities that may represent subducted slices of oceanic crust and mantle. This layered series may continue northward under the Peninsular terrane. Earthquake locations in the Wrangell Benioff zone indicate that at least the upper two low-high velocity layer pairs are tectonically inactive and that they appear to have been accreted to the base of the continental crust. The refraction data suggest that the Contact fault between two similar terranes, the Chugach and Prince William terranes, is a deeply penetrating feature that separates lower crust (deeper than 10 km) with paired dipping reflectors, from crust without such reflectors.     

Upper-crustal structure of the Benevento area (southern Italy): fault heterogeneities and potential for large earthquakes

The Benevento region is part of the southern Apennines seismogenic belt, which experienced large destructive seismic events both in historical and in recent times. The study area lies at the northern end of the Irpinia fault, which ruptured in 1980 with a M_s = 6.9 normal faulting event, which caused about 3000 casualties. The aims of this paper are to image lateral heterogeneities in the upper crust of the Benevento region, and to try to identify the fault segments that are expected to generate such large earthquakes. This work is motivated by the recognition that lithological heterogeneities along major fault zones, inferred from velocity anomalies, reflect the presence of fault patches that behave differently during large rupture episodes. In this paper, we define the crustal structure of the Benevento region by using the background seismicity recorded during 1991 and 1992 by a local seismic array. These data offer a unique opportunity to investigate the presence of structural discontinuities of a major seismogenic zone before the occurrence of the next large earthquake. The main result that we obtained is the delineation of two NW-trending high-velocity zones (HVZs) in the upper crust beneath the Matese limestone massif. These high velocities are interpreted as high-strength regions that extend for 30-40 km down to at least 12 km depth. The correspondence of these HVZs with the maximum intensity regions of historical earthquakes (1688 AD, 1805 AD) suggests that these anomalies delineate the extent of two fault segments of the southern Apenninic belt capable of generating M = 6.5⁻⁷ earthquakes. The lateral offset observed between the two segments from tomographic results and isoseismal areas is possibly related to transverse right-lateral faults.