Magnetotelluric studies in the Market Weighton area of eastern England

Summary Magnetotelluric measurements at periods from 30 to 1000s were made at eight locations in the Market Weighton (MW) area, along an east—west profile across gravity and magnetic anomalies. Dimensional parameters were developed for assessing the structural dimensionality of the electrical conductivity of the Earth from the data. One-dimensional inversion modelling techniques were employed to interpret the data at each site, and four-layer models were obtained to explain the main structure of the crust in the area studied. If it is assumed that all strata are unmagnetized then the results show that there is a highly resistive layer in the crust, the thickness of the highly resistive layer ranges from 12 km in the east to 44 km in the west with a large change in the middle near the MW site. A structural boundary lying north—south near MW was also indicated by the principal directions of rotated apparent resistivities and transfer functions. Both electrical conductivity and magnetic permeability contrast in the ground were considered in an attempt to interpret the observed variations in apparent resistivity at different periods.     

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.     

西南极菲尔德斯半岛长城站地区深部电性结构

本文分析了1992年12月长城站附近三个测点的大地电磁测深资料(1s到4096s的低频信号数据),得出菲尔德斯半岛风谷断裂电性主轴为北东110°,半岛地壳厚度为22.3km,壳内含四个主要电性层,厚度分别为:1.3km,6.7km,1.2km和13.3km,其中上地壳为9.2km,下地壳为13.1km。     

Electrical conductivity anomalies and their relationship with the tectonics of South Australia

Summary. Geomagnetic variation studies have been conducted in the Gawler Craton and Adelaide Geosyncline of South Australia. The magnetometer stations extend from the coast up to the southern edge of the 1970 array of Gough, McElhinny & Lilley. The coast effect is the dominant feature of the data but use is made of the hypothetical event technique to identify two zones of telluric current concentration. Both of them appear to be associated with linear zones of enhanced electrical conductivity within the crust. The Southern Eyre Peninsula anomaly lies within the Gawler Craton and may identify a major fracture or shear in the upper crust. The conductivity anomaly within the Adelaide Geosyncline appears to be the continuation of the Flinders anomaly discovered by the 1970 array study. It correlates well with the arcuate fold pattern of the Southern Flinders Zone of the Geosyncline and with the local pattern of seismicity. In both anomalies the enhanced conductivity is probably caused by saline waters within fractured crustal rocks.     

南极区和北极区高空电流体系的比较研究

本文用ＩＧＹ／ＩＧＣ期间全球地磁台网的资料计算出地磁太阳日变化（Ｓ）和太阴日变化（Ｌ）的电流体系，对比分析了南极区与北极区电流体系的特点。分析表明：（１）两极区的外源电流体系存在明显差别，这反映了产生该电流系的发电机过程（对Ｓ和Ｌ）和场向电流（对Ｓ）的不同。两极区磁场结构的特征可能是导致这一差异的根本原因。（２）两极区内源电流存在明显差异，这一方面归因于外源施感场（电流）的差异，另一方面也反映了两极区地下电导率的不同。分析表明，从总体来看，南极区地下电导率高于北极区     

Resistivity studies over the Flinders conductivity anomaly, South Australia

Summary. Seven Schlumberger resistivity soundings with maximum current electrode spacings of 20 km have been conducted south of Lake Frome in South Australia. These experiments were done partly to test new electrical sounding equipment and partly to investigate a large conductivity anomaly previously delineated by other workers using magnetometer array and MT methods (the 'Flinders'anomaly). These previous studies left some doubt as to the depth to the conductive region responsible for the anomaly.
The electrical soundings did not detect a buried conductive zone, which constrains it to lie deeper than 5–7 km. However, the study did show the surface sediments of the region to be very conductive; resistivities of 2–9 μm were measured over thicknesses of 50–400 m, with sediment thickness inferred to be up to 2 km to the north of the studied area. This raises the question of whether current channelling in the surface sediments could have been responsible for the earlier results. Simple modelling and application of the criteria given by Jones suggest this may be so.
The equipment used for this study is a low power (200 W), computer controlled system which employs synchronous stacking and other signal processing to achieve signal to noise improvement ratios of up to 1000.     

East African earthquakes below 20 km depth and their implications for crustal structure

We report source parameters for eight earthquakes in East Africa obtained using a number of techniques, including (1) inversion of long-period P and SH waves for moment tensors and source-time functions, (2) forward modelling of first-motion polarities and P and pP amplitudes on short-period seismograms, and (3) determination of pP-P and sP-P differential traveltimes from short-period records. The foci of these earthquakes lie between depths of 24 and 34 km in Archean and Proterozoic lithosphere, and all but one fault-plane solution indicates normal faulting (primarily E-W extension), consistent with the regional stress regime in East Africa. Because many of these earthquakes occurred in areas where the crust may have been thinned by rifting, it is difficult to ascertain whether or not their foci lie within the lower crust or upper mantle. Some of them, however, occurred away from rift structures in Proterozoic crust that is possibly 35–40 km thick or thicker, and thus they probably nucleated within the lower crust. Strength profile calculations suggest that in order to account for seismogenic (i.e. brittle) behaviour at sufficient depths to explain lower crustal earthquakes in East Africa, the lower crust must not only be composed of mafic lithologies, as suggested by previous investigators, but also that significantly more heat (∼100 per cent) must come from the upper crust than predicted by the crustal heat source distribution obtained from a 1-D interpretation of the linear relationship between heat flow and heat production observed in Proterozoic terrains within eastern and southern Africa. Precambrian mafic dike swarms throughout East Africa provide evidence for magmatic events which could have delivered large amounts of mafic material to the lower crust over a very broad area, thus explaining why the lower crust in East Africa might be mafic away from the volcanogenic rift valleys.     

Separation of local and regional information in distorted GDS response functions by hypothetical event analysis

Electromagnetic investigations are usually intended to examine regional structures where induction takes place at a given period range. However, the regional information is often distorted by galvanic effects at local conductivity boundaries. Bahr (1985) and Groom & Bailey (1989) developed a physical distortion model for decomposing the MT impedance tensor, based upon local galvanic distortion of a regional 2-D electromagnetic field. We have extended their method to predict the magnetic variation fields created at an array of sites. The magnetic response functions at periods around 1000 s may be distorted by large-scale inhomogeneities in the upper or middle crust. In this period range, the data measured by a magnetometer array contain common information that can be extracted if the data set is treated as a unit, for example by using hypothetical event analysis. With this technique it is always possible to recover the regional strike direction from distorted data, even if a strong, spatially varying regional vertical field component is present in the data set. The determination of the regional impedance phases, on the other hand, is far more sensitive to deviations from the physical distortion model.
The approach has been used to investigate the Iapetus data set. For the array, which covers an area of 200  km × 300  km in northern England/southern Scotland, the technique revealed a common regional strike azimuth of ca . N125° E in the period range 500–2000  s. This direction differs from the strike indicated by the induction arrows, which seem influenced mainly by local current concentrations along the east–west-striking Northumberland Trough and a NE–SW-striking mid-crustal conductor. Both impedance phases are positive and differ by ca . 10°, which supports the assumptions of distortion fields in the data set and that the regional structure is 2-D.     

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.     

Electromagnetic induction studies in the Eyre Peninsula, South Australia

Magnetic field fluctuations have been recorded by an array of portable three-component magnetometers at 60 sites across the Eyre Peninsula in South Australia between December 1993 and March 1995. An additional 54 magnetometer data records, collected prior to 1989 and described by Milligan (1989) and Milligan, White & Chamalaun (1989), were included in the analysis. A major conductive feature in the crust, first noted by White & Milligan (1984) as the Eyre Peninsula Anomaly (EPA), is re-examined to assess its continuity to the north of the original arrays and to investigate its relationship with major tectonic features.
Magnetic-field time-series were converted to induction arrows in the frequency domain. These induction arrows were initially inverted using the minimum-structure 2-D Occam approach to estimate the electrical conductance of the crust. Following this, thin-sheet forward modelling was used to examine the relationship between the conductance and the dominant tectonic features. The principal results of the modelling are that a narrow conductive feature extends inland from the coast about 160 km before terminating, and the conductance is in the range 3000 to 10 000 S, which decreases inland.
A strong correlation exists between the electrical conductance of the Eyre Peninsula and Bouguer gravity anomalies, and in particular the EPA is coincident with a significant Bouguer gravity gradient. There is also good agreement between the locations of the foci of earthquakes of magnitude greater than 4.0 and the EPA. We believe that the anomaly is associated with a geological fracture in the Precambrian upper crust as a result of crustal extension prior to the rifting of Australia from Antarctica in the Jurassic (160 Ma).     

Magnetometer array study in a tectonically active region of Quebec, Canada

Summary. An array of 26 magnetometers deployed in a tectonically active area of Quebec has mapped a boundary in terrestrial electrical conductivity for 200 km along strike. The contrast in conductivity across the boundary, from previous magnetotelluric soundings, is about one order of magnitude. Anomalous variation fields associated with electric currents flowing along the boundary are readily detected at pulsation periods only when the horizontal field is polarized transverse to the structure (the E -polarization case). The anomaly is hardly visible in transfer functions from substorms, for a number of reasons: a predominant H -polarization orientation of the substorm fields, the single order-of-magnitude contrast in conductivity, and the probably small depth extent of the structure. Attempts were made to estimate the response of a one-dimensional earth via the inductive scale length with gradients evaluated from polynomial surfaces fitted to the smoothly varying substorm fields. The results were inconsistent, owing to vertical fields with strong external components and to horizontal fields with scale lengths too small relative to their penetration distances.     

Upper-mantle structure of the Baltic Shield below the Swedish National Seismological Network (SNSN) resolved by teleseismic tomography

Upper-mantle structure under the Baltic Shield is studied using non-linear high resolution teleseismic P -phase tomography. Observed relative arrival-time residuals from 52 teleseismic earthquakes recorded by the Swedish National Seismological Network (SNSN) are inverted to delineate the structure of the upper mantle. The network consists of 47 (currently working) three-component broad-band stations located in an area about 450 km wide and 1450 km long. In order to reduce complications due to possible significant three-dimensionality of Earth structure, events chosen for this study lay close to in-line with the long-axis of the array  (±30°)  . Results indicate P -wave velocity perturbations of ±3 per cent down to at least 470 km below the network. The size of the array allows inversion for structures even at greater depths, and lateral variations of velocity at depths of up to 680 km appear to be resolved. Below the central part of the array (60°–64° N), where ray coverage is best, the data reveals a large region of relatively low velocity at depths of over about 300 km. At depths less than about 250–300 km, the models include a number of features, including an apparent slab-like structure dipping gently towards the north.     

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.     

Crustal root beneath the highlands of southern Norway resolved by teleseismic receiver functions

Teleseismic data have been collected with temporary seismograph stations on two profiles in southern Norway. Including the permanent arrays NORSAR and Hagfors the profiles are 400 and 500 km long and extend from the Atlantic coast across regions of high topography and the Oslo Rift. A total of 1071 teleseismic waveforms recorded by 24 temporary and 8 permanent stations are analysed. The depth-migrated receiver functions show a well-resolved Moho for both profiles with Moho depths that are generally accurate within ±2 km.
For the northern profile across Jotunheimen we obtain Moho depths between 32 and 43 km (below sea level). On the southern profile across Hardangervidda, the Moho depths range from 29 km at the Atlantic coast to 41 km below the highland plateau. Generally the depth of Moho is close to or above 40 km beneath areas of high mean topography (>1 km), whereas in the Oslo Rift the crust locally thins down to 32 km. At the east end of the profiles we observe a deepening Moho beneath low topography. Beneath the highlands the obtained Moho depths are 4–5 km deeper than previous estimates. Our results are supported by the fact that west of the Oslo Rift a deep Moho correlates very well with low Bouguer gravity which also correlates well with high mean topography.
The presented results reveal a ca . 10–12 km thick Airy-type crustal root beneath the highlands of southern Norway, which leaves little room for additional buoyancy-effects below Moho. These observations do not seem consistent with the mechanisms of substantial buoyancy presently suggested to explain a significant Cenozoic uplift widely believed to be the cause of the high topography in present-day southern Norway.     

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.     

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

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

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.     

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.     

Relocation of M ≥ 2 events of the 1989 Dôbi seismic sequence in Afar: evidence for earthquake migration

The Ethiopian side of central Afar was struck in August 1989 by the largest seismic sequence (three 6.1 ≤ M _s ≤ 6.3 events, 15 with M _s or m _b ≥ 5.0) since that of Serdo in 1969. Using the Djibouti seismological network, we relocated 297 of the events of that sequence. As most of the large events took place outside the network, we assessed the accuracy and stability of earthquake relocations by using three different velocity models and two relocation codes to try to relate individual shocks to distinct faults and surface breaks. A majority of the events apparently occurred underneath the floor of the Dôbi graben, an area about 45  km long and 15  km wide, rupturing boundary and inner floor faults, in agreement with the surface cracks and scarps that we mapped in the area. The relocation shows that the principal events propagated about 50  km northwestwards along the graben in the first 40  hr. A day and a half after the beginning of the sequence, smaller events ( M ≤ 4) started to propagate more than 55  km eastwards, towards Asal Lake. Using two three-component stations installed near the Ethiopian border, we could determine reliable depths for 21 events. The depths are compatible with a seismogenic crust about 14  km thick in the Dôbi and Hanle graben area. Although the Dôbi sequence ruptured about 50  km of the fault array extending from Serdo to Asal, where the regional stress was released by earthquakes in 1969 and 1978, respectively, a seismic gap about 50  km long still subsists along the northern part of the Gaggade region (Der'êla half-graben).     

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.