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.     

Conductivity anomalies in the Baltic Shield in Finland

In 1981–84 six arrays of 30 or 31 magnetometers were operated on the Baitic Shield in central and southern Finland to deduce information about the electrical conductivity within the Earth's crust. The magnetometer sites cover the area with 20–60 km spacing allowing the identification of large-scale electrical structures within the crust. Each array recorded the three orthogonal magnetic field components on analogue film for about two months simultaneously at all the stations. Altogether 17 magnetic variation events of 2–6hr length have been digitized and analysed.
This paper includes results of the last array in south-western Finland, a comparison of two methods of induction vector estimation and combination of results from all the arrays. Four zones of anomalously high electrical conductivity have been observed. They are all interpreted as lying at depths greater than 5 km. The 'Oulu Anomaly' is wider than the others and rather short, about 100 km. The 'Southern Finland Anomaly' is probably a 500 km-long narrow zone. This and the 'Ostrobothnian Anomaly' surround three sides of the large Central Finland Granite area. The 'Outokumpu Anomaly' is probably deeper lying than the others. Blocks separated by these anomalies are presented. The blocks have different average resistivities.     

Magnetometer array studies and deep Schlumberger soundings in the Damara orogenic belt, South West Africa

Summary. A zone of concentrated induced electric currents crossing parts of Zimbabwe, Botswana and South West Africa was discovered during a magnetovariational study conducted in 1972. In 1977, a second study was made with 27 recording magnetometers distributed across the width of South West Africa between latitudes 19 and 22°S. Several geomagnetic disturbances were recorded with high recording efficiencies. Three of these time sequences were digitized for analysis. Magnetograms and Fourier transform amplitude and phase maps in the period range 22–128min were used to delineate the westward continuation of the conductive structure revealed by the earlier investigation. The conductive zone runs approximately east-west from the Botswana border (21°E) to 17°E longitude. From here to the Atlantic coast it trends in a NE—SW direction. Anomalous fields, normalized to the horizontal field at a station recording the normal field, were used to obtain maximum depth estimates of around 45 km for the induced currents. Several deep Schlumberger soundings were done over the anomalous zone and the results showed that the conductive structure is, in places, only 3 km from the surface and that it has a resistivity of less than 20 Ωm. The resistivity of the upper crust outside the structure ranges from 5000 to more than 20000 Ωm. Some 14 post-Karoo alkaline igneous complexes occur along the course of the resistivity anomaly. These intrusive complexes represent the youngest igneous activity in the Damara Orogenic Belt and were most probably emplaced along a line of weakness in the lithosphere. The resistivity anomaly would seem to delineate this line of weakness.     

Crustal structure of Iraq from receiver functions and surface wave dispersion: implications for understanding the deformation history of the Arabian–Eurasian collision

We report the crustal structure for two locations in Iraq estimated by joint inversion of P -wave receiver functions (RFs) and surface (Rayleigh) wave group velocity dispersion. RFs were computed from teleseismic recordings at two temporary broad-band seismic stations located in Mosul (MSL) in the Zagros Fold Belt and Baghdad (BHD) in the Mesopotamian Foredeep. Group velocity dispersion curves at the sites were derived from continental-scale tomography. The inversion results show that the crustal thicknesses are 39 km at MSL and 43 km at BHD. We observe a strong Ps _Moho at BHD consistent with a sharp Moho discontinuity. However, at MSL we observe a weak Ps _Moho suggesting a transitional Moho where crustal thickening is likely to be occurring in the deep crust. Both sites reveal low velocity surface layers consistent with sedimentary thickness of about 3 km at station MSL and 7 km at BHD and agreeing well with the previous reports. Ignoring the sediments, the crystalline crustal velocities and thicknesses are remarkably similar at both stations. The similarity of crustal structure suggests that the crust of the northeastern proto-Arabian Platform was uniform before subsidence and deposition of the sediments in the Cenozoic. If crystalline crustal structure is uniform across the northern Arabian Platform then crustal thickness variations in the Zagros Fold Belt and Thrust Zone should reveal the history of deformation and crustal shortening in the Arabian–Eurasian collision zone and not reflect pre-existing crustal thickness variations in the Arabian Plate.     

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.     

3-D magnetotelluric inversion and model validation with gravity data for the investigation of flood basalts and associated volcanic rifted margins

20 magnetotelluric (MT) soundings were collected on the Isle of Skye, Scotland to provide a high-resolution three-dimensional (3-D) electrical resistivity model of a volcanic province within the framework of a project jointly interpreting gravity, seismic, geological and MT data. The full 3-D inversion of the MT data jointly interpreted with gravity data reveals upper crustal structure. The main features of the model are interpreted in conjunction with previous geological mapping and borehole data. Our model extends to 13 km depth, several kilometres below the top of the Lewisian basement. The top of the Lewisian basement is at approximately 7–8 km depth and the topography of its surface was controlled by Precambrian rifting, during which a 4.5 km thick sequence of Torridonian sediments was deposited. The Mesozoic sediments above, which can reach up to 2.2 km thick, have small-scale depocentres and are covered by up to 600 m of Tertiary lava flows. The interpretation of the resistivity model shows that 3-D MT inversion is an appropriate tool to image sedimentary structures beneath extrusive basalt units, where conventional seismic reflection methods may fail.     

Three-dimensional image of the Moho undulations beneath the Gulf of Bothnia using wide-angle seismic data

The BABEL marine seismic experiment has been carried out to investigate the lithospheric structure and antecedent tectonic signatures of the Baltic Shield, including the Archaean-Proterozoic collisional structure in the northern part of the Gulf of Bothnia.
Lithospheric seismic-reflection streamer data and simultaneously recorded wide-angle reflection and refraction data collected in the Gulf of Bothnia as part of the BABEL project have been used for 3-D modelling. The distribution of land stations around the Gulf provides a good 3-D ray coverage of the PMP reflection data recorded at the eight stations in the area and allows an estimation of strikes and dips of the Moho boundary in the area. The traveltimes of reflected phases are calculated using a method that utilizes the finite-difference solution of the eikonal equation. The Moho wide-angle-reflection (PMP) traveltimes are modelled using an inversion method. A 2-D model from the Gulf of Bothnia extended into the third dimension is used as an initial model. During the inversion the velocity is kept constant and only the Moho boundary is allowed to vary. To estimate the strike of the Moho boundary and the stability of the inversion, two initial models with different strikes are examined.
The results indicate that the Moho depth in the Gulf of Bothnia undulates and has a maximum depth of 55 km in the south, rising to 42 km in the north. The Moho depth variations seem to be step-like. This change in the Moho depth coincides with the location of the presumed fossil subduction zone in the area. The crustal-thickness variations seem to be well approximated by a nearly 2-D structure striking parallel to a postulated subduction zone immediately to the south of the Skellefte area. The presence of the step at the crust/mantle boundary can be interpreted as a result of a plate-collision event at about 2 Ga.     

Deep Schlumberger sounding and the crustal resistivity structure of central Australia

Summary. Three 200 km Schlumberger resistivity soundings have been conducted over the central Australian shield, using telephone lines to obtain the large electrode spacings. These represent the first crustal scale controlled source electrical study to be carried out in this continent. A computer controlled data acquisition system was used which allowed precise measurements to be made with only modest emission currents (0.1–0.5 A).
The three soundings, centred on the towns of Renner Springs, Wauchope and Aileron, showed the southern part of the study area (the Arunta Block) to be an order of magnitude more resistive than the more northerly section (the Tennant Creek Block). This difference correlates with the higher heat flow of the Tennant Creek Block. A lowering of apparent resistivity at large electrode spacings for one sounding (Wauchope) is taken to indicate the presence of a low resistivity layer in the middle crust, at a depth less than 20 km. However, the effect of the highly conductive overburden characteristic of inland Australia, combined with the large transverse resistance of the crust, prevented the other two soundings from detecting such a layer. Without support from these two soundings, it is impossible to be sure that the lowered resistivity at Wauchope is not caused merely by lateral variations in near-surface resistivity.
The data also show that crustal resistivities are much lower than the expected values for dry rock, whether or not a low resistivity layer is included in the model. This implies a widespread occurrence of free water in the crust, with greater amounts occurring at depth if the low resistivity zone exists.     

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

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

Results from regional vibroseis profiling: Appalachian ultra-deep core hole site study

Summary. In 1985, 180 km of regional vibroseis profiles were acquired in the Carolinas and Georgia, southeastern United States, as part of the Appalachian Ultra-Deep Core Hole (ADCOH) Site Study. The data quality is excellent, with large-amplitude reflections from faults and crystalline rocks, lower Palaeozoic shelf strata and from within autochthonous Grenville basement. The profiles image the subsurface more clearly than other available data and allow the possibility of alternative interpretations of important elements of the tectonic framework of the southern Appalachians.
The major points in the interpretation are: 1) The Blue Ridge master decollement is at a depth of 2-3 km beneath the Blue Ridge. This thrust increases in dip just NW of the Brevard fault zone. 2) The Brevard fault zone appears to splay from the master decollement at 6 km (2.2 s) near Westminster, S.C., and defines the base of the crystalline Inner Piedmont allochthon. 3) Below the Blue Ridge thrust sheet are images of duplex and imbricate structures ("duplex tuning wedges") connected by other thrust faults that duplicate shelf strata to a thickness of 4–5 km. 4) Subhorizontal reflections from depths of 6 to 9 km may be from relatively undisturbed lower Palaeozoic strata as suggested by others. 5) Eocambrian-Cambrian(?) rift basins in the Grenville basement are also imaged.
The ADCOH data were originally recorded with 14–56 Hz bandwidth and 8 s length, but an extended Vibroseis correlation was used to produce 17 s data length revealing reflections from within the upper crust. Below 8 s, reflections from within the Grenville basement become weak, but are observable as late as 13 s; however, these Moho (?) reflections are generally short segments.     

Results of magnetotelluric and gravimetric measurements in western Nicaragua, Central America

Magnetotelluric and gravity data have been collected within a ca. 170 km long traverse running from the Pacific coast of Nicaragua in the west to the Nicaraguan Highland in the east. This part of Nicaragua is characterized by sedimentary rocks of the Pacific Coastal Plain, separated from the Tertiary volcanic rocks of the Highland by the NW-SE-trending Nicaraguan Depression. 2-D interpretation of the magnetotelluric (MT) data, collected at 13 stations, indicates four regions of high electrical conductivity in addition to the conductive coastal region. Two of these are associated with conducting sediments and pyroclastics in the upper part of the crust. Two other conductive structures have been defined at depth around 20 km and the one best defined is located below the depression. From the distribution of seismic events, volcanic activity in the depression and the similarity in geophysical characteristics with areas such as the Rio Grande Rift, this conductor is interpreted as a melt layer or a complex of magma chambers. Models of the upper lithosphere, constrained by the MT model, vertical electrical sounding (VES) data, seismic data and densities, have been tested using gravity data. A model that passes this test shows a gradual thickening of the crust eastwards from the Pacific coast. An anomaly centred over the depression is interpreted to have its origin in a thinning of the crust. In this model the melt layer is situated on top of the bulge of the lower lithosphere. A change in the composition of the crust, from the Pacific Coastal Plain to the Highland, is indicated from the change in character of the MT response and from the density distribution in the gravity model. This may support the hypothesis that the Pacific region is an accreted terrane. MT and gravity data indicate a depth to a resistive and high-density basement in the depression of ca. 2 km. On the basis of this, the vertical setting in the depression is estimated to be of the order of 2.5 km.     

Application of Magnetotelluric (MT) Resistivity to Imaging of Regional Three-Dimensional Geologic Structures and Groundwater Systems

Widespread definition of a groundwater system in three dimensions is necessary for the management and maintenance of groundwater resources. A magnetotelluric (MT) survey can be an effective geophysical prospecting method for imaging regional geological structures by measuring both shallow and deep resistivity. To demonstrate the capability of an MT survey to characterize a groundwater system, the Kumamoto area of central Kyushu in southwestern Japan was selected as a case study site because of its rich groundwater resources. Three-dimensional (3D) MT resistivity structure to a depth of 5?km was modeled by 1D inversion analysis of raw MT data and 3D interpolation of the resultant resistivity column data by the optimization principle method. Consequently, both deep and shallow aquifers were detected. A high-resistivity zone appears at depths between 500 and 2,000?m between the Futagawa?CHinagu faults and the Usuki?CYatsushiro tectonic line, which supports the existence of an aquiclude under the aquifer. The most important characteristic inferred from the 3D resistivity model is that the deep groundwater system below a depth of 1,000?m has two main flow paths. One path is likely to be through porous rocks because the low resistivity zone is regarded as tuff with sand and gravel, and the other flow path is interpreted to be through fractured zones along the Hinagu faults. Hence, the path and direction of the groundwater flows are probably controlled by geologic structures and the configuration of the active faults. These findings support the effectiveness of the MT method for investigating groundwater systems.     

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.     

Geophysical Models for the Cu-Dominated VHMS Mineralization in Katta District, Western Ethiopia

Volcanic hosted massive sulfides (VHMS) have been recognized for many years within the Arabian–Nubian Shield of western Ethiopia. The Shield is a collision belt formed when East and West Gondwana collided during the Neoproterozoic time. It covers NE Africa and the Middle East and is known to have hosted several economic mineral deposits. Situated within the Shield, the Katta VHMS mineralization is primarily made up of chalcopyrite, sphalerite, pyrite, pyrrhotite, and magnetite. These minerals are known to have contrasting values in resistivity, chargeability, and magnetic susceptibility with the enclosing host rocks. These physical attributes make resistivity, IP, and magnetic methods ideal to study and detect the mineralization and possibly delineate the ore body. Dipole–dipole DC resistivity and induced polarization data collected along parallel profiles are inverted to generate both 2D and 3D models of resistivity and chargeability. The chargeability model successfully delineates the zone of sulfide mineralization and the resistivity model gives information on the stratigraphic horizons and the structures. Magnetic data collected along more closely spaced profiles were also inverted to recover the 3D magnetic susceptibility model. The recovered model shows a linear magnetic anomaly that could be associated with the sulfides. This anomaly shows good spatial continuity and is in good agreement with the result obtained in the chargeability model. The geophysical inversion results indicate a NNW oriented mineralized zone of considerable strike length and extending to a depth of about 380 m. This study suggests the presence of exploitable VHMS presence in Katta district and that the inversion methodology is an important interpretation tool in guiding the exploration and exploitation of sulfide mineralization.     

The determination of resistivity and porosity of the sediment and fractured basalt layers near the Juan de Fuca Ridge

Summary. In the summer of 1984 an electrical survey using magnetometric off-shore electrical sounding (MOSES) was conducted at two sites in Middle Valley, part of the northern Juan de Fuca Ridge complex. MOSES has been designed to minimize the difficulties inherent in electrical surveys of the crust below the electrically conductive sea layer. Site 1, at 48°32N, 128°42W, is in the central part of the turbidite-filled basin. Using a two-layer model of conductive sediments overlying a fractured basalt basement, the sediment resistivity and thickness were found to be 0.82 ± 0.06 Ωm and 1800 ± 300 m, respectively. The basement resistivity, although not well constrained by the data is consistent with the results obtained at site 2.
Site 2, located at 48°10N, 128°50W, has a thinner sediment layer, which appears to vary with position. The sediment conductivity—thickness product is the parameter determined by the data. If the sediment resistivity were the same as at site 1, the sediment thickness would be 140 ± 30 m to the SE of site 2, and 240 ± 55 m to the NW. The fractured basalt basement has a resistivity of 8.5 ± 3.4 Ωm and is at least 1000 m thick.
Using temperature-corrected pore fluid resistivity, the calculated porosity is found to vary from 62 per cent at the top to 21 per cent at the base of the sediments and is 8 per cent in the basement. These values are in good agreement with estimates from seismic velocities for a thick turbidite sequence in a nearby sediment-filled basin and determined for layer 2A/B basalts in DSDP hole 504B, respectively.     

Accurate hypocentre determination in the seismogenic zone of the subducting Nazca Plate in northern Chile using a combined on-/offshore network

The coupled plate interface of subduction zones—commonly called the seismogenic zone—has been recognized as the origin of fatal earthquakes. A subset of the after-shock series of the great Antofagasta thrust-type event (1995 July 30; M _w = 8.0) has been used to study the extent of the seismogenic zone in northern Chile. To achieve reliable and precise hypocentre locations we applied the concept of the minimum 1-D model, which incorporates iterative simultaneous inversion of velocity and hypocentre parameters. The minimum 1-D model is complemented by station corrections which are influenced by near-surface velocity heterogeneity and by the individual station elevations. By relocating mine blasts, which were not included in the inversion, we obtain absolute location errors of 1  km in epicentre and 2  km in focal depth. A study of the resolution parameters ALE and DSPR documents the importance of offshore stations on location accuracy for offshore events. Based on precisely determined hypo-centres we calculate a depth of 46  km for the lower limit of the seismogenic zone, which is in good agreement with previous studies for this area. For the upper limit we found a depth of 20  km. Our results of an aseismic zone between the upper limit of the seismogenic zone and the surface correlates with a detachment zone proposed by other studies; the results are also in agreement with thermal studies for the Antofagasta forearc region.     

Receiver functions at the stations of the German Regional Seismic Network (GRSN)

Several years of broad-band teleseismic data from the GRSN stations have been analysed for crustal structure using P -to- S converted waves at the crustal discontinuities. An inversion technique was developed which applies the Thomson-Haskell formalism for plane waves without slowness integration. The main phases observed are Moho conversions, their multiples in the crust, and conversions at the base of the sediments. The crustal thickness derived from these data is in good agreement with results from other studies. For the Gräfenberg stations, we have made a more detailed comparison of our model with a previously published model obtained from refraction seismic experiments. The refraction seismic model contains boundaries with strong velocity contrasts and a significant low-velocity zone, resulting in teleseismic waveforms that are too complicated as compared to the observed simple waveforms. The comparison suggests that a significant low-velocity zone is not required and that internal crustal boundaries are rather smooth.     

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

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

Magnetotelluric Soundings In the Canadian Rocky Mountains

Summary. Magnetotelluric soundings have been made at 25 stations in the Rocky Mountain Trench (RMT) and Main Ranges near 53° N, close to the centre of a major conductivity anomaly which had been mapped in a magnetovariation array study. Most stations covered the frequency range 0.01–500 Hz and three stations 0.0002–500 Hz. the resistivity tensor shows low to moderate anisotropy in the RMT, but is strongly 2-D or 3-D in the Rocky Mountains. Apparent resistivities as a function of frequency are displayed in pseudosections along the Trench and along a transverse profile across the RMT and into the Main Ranges. In preparation for 2-D modelling, 1-D inversions have been used to construct resistivity-depth sections satisfying both magnitudes and phases of the MT responses. These show very low resistivities, in the range 1—10Ωm, in the upper crust under the RMT and even lower values under the Main Ranges. the latter values give strong confirmation of the Northern Rockies conductor reported by Bingham, Cough & Ingham and are in agreement with models of the conductors fitted to long-period magnetovariation fields by Ingham, Gough & Parkinson. the MT results here reported add some essential depth and resistivity information. It is suggested that the conductors beneath the Rocky Mountains Main Ranges and Trench constitute a thickening at the edge of the Canadian Cordilleran Regional (CCR) conductor. Gough has argued that a wide variety of geophysical and geological parameters indicate high temperatures and partial melting in the mantle under the CCR conductor. At the upper crustal depths penetrated in this magnetotelluric study, it is considered more probable that the high conductivity is caused by hot, saline water of mantle origin rather than silicate melt. the CCR in general may have two layers of fluid producing its high conductivity, silicate melt below and saline hot water above.     

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.