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.     

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.     

Source mechanism of the deep Colombian earthquake of 1970 July 31 from the free oscillation data

Summary. The method proposed by Mendiguren to determine the source parameters from free oscillation data is applied to the 1970 July 31 deep Colombian earthquake. The results indicate a source propagating horizontally for about 150 km along the lithosphere and cutting across its width. The slab behaves as a guide for source propagation. The horizontal propagation velocity is determined as 3.8 km/s. The intensity of the source grew proportionately to the second power of the propagation distance. This rate of source intensity growth may be interpreted either by a fan-shaped fault model or by a cone-shaped volume source. The average slip and stress drop are estimated as 360 cm and 300 bar for the fault model. For the volume source model the transformational shear strain and stress are estimated as 11 × 10⁻⁵ and 160 bar. There is no evidence of a double couple radiation preceding the P origin time. It is shown that the isotropic and deviatoric components of the moment tensor cannot be uniquely resolved when only observations of a single mode are available. It is observed that, from a statistical basis, the available ₀ S_n data for Colombian shock can be equally well explained by a pure deviatoric source model or by a source model including an isotropic component. Numerical experiments indicate that the inclusion of higher mode data does not change this situation. But, on the other hand, numerical experiments show that the available data and the scheme used for the inversion would not result in a solution including an artificial implosive component if the actual source were pure deviatoric. If the departure from a pure deviatoric source is produced by noise, it has to be non-random, as it could be produced by lateral heterogeneities not included in the inversion scheme.     

A three-dimensional radially anisotropic model of shear velocity in the whole mantle

We present a 3-D radially anisotropic S velocity model of the whole mantle (SAW642AN), obtained using a large three component surface and body waveform data set and an iterative inversion for structure and source parameters based on Non-linear Asymptotic Coupling Theory (NACT). The model is parametrized in level 4 spherical splines, which have a spacing of ∼ 8°. The model shows a link between mantle flow and anisotropy in a variety of depth ranges. In the uppermost mantle, we confirm observations of regions with   V_SH > V_SV   starting at ∼80 km under oceanic regions and ∼200 km under stable continental lithosphere, suggesting horizontal flow beneath the lithosphere. We also observe a   V_SV > V_SH   signature at ∼150–300 km depth beneath major ridge systems with amplitude correlated with spreading rate for fast-spreading segments. In the transition zone (400–700 km depth), regions of subducted slab material are associated with   V_SV > V_SH   , while the ridge signal decreases. While the mid-mantle has lower amplitude anisotropy (<1 per cent), we also confirm the observation of radially symmetric   V_SH > V_SV   in the lowermost 300 km, which appears to be a robust conclusion, despite an error in our previous paper which has been corrected here. The 3-D deviations from this signature are associated with the large-scale low-velocity superplumes under the central Pacific and Africa, suggesting that   V_SH > V_SV   is generated in the predominant horizontal flow of a mechanical boundary layer, with a change in signature related to transition to upwelling at the superplumes.     

Mantle transition zone thickness beneath Afar: implications for the origin of the Afar hotspot

The velocity spectrum stacking method is applied to receiver functions from stations ATD and AAE to image P -to- S converted phases originating at the 410 and 660 km discontinuities beneath Afar. A transition zone thickness of 244 ± 19 km is obtained, similar to the global average transition zone thickness. This result suggests that any broad thermal anomaly beneath Afar probably does not extend as far down as the transition zone. However, because of the 19 km uncertainty in the thickness estimate, a small thermal anomaly of ~100–150 K at mantle transition zone depths cannot be ruled out.     

2001 August earthquake swarm at Shadwan Island, Gulf of Suez, Egypt

The earthquake swarm that struck Shadwan Island at the entrance of the Gulf of Suez in 2001 August included 408 events. Almost all of these events (94 per cent) were microearthquakes and only 6 per cent had small measurable magnitudes  (5.0 > M _L≥ 3.0)  . Most of the earthquakes were weak and followed each other so closely in time that they could not be identified at more distant stations. The fault plane solutions of the strongest events of the swarm show almost identical focal mechanisms, predominately normal faulting with a significant sinistral strike-slip component for nodal planes trending NW–SE. A comparison with the mechanisms of the 1969 and 1972 events which took place 20 km north of the swarm epicentral region shows similarities in faulting type and orientation of nodal planes. The azimuths of T -axes determined from focal mechanisms in this study are oriented in the NNE–SSW direction. This direction is consistent with the present-day stress field derived from borehole breakouts in the southern Gulf of Suez and the last phase of stress field changes in the Late Pleistocene, as well as with recent GPS results.
The source parameters of the largest  ( M _L≥ 3.0)  events of the 2001 August Shadwan swarm have been estimated from the P -wave spectra of the Egyptian National Seismograph Network (ENSN). Averaging of the values obtained at different stations shows relatively similar source parameters, including a fault length of  0.65 ≤ L ≤ 2 km  , a seismic moment of  7.1 × 10¹²≤ Mo ≤ 3.0 × 10¹⁴ N m  and a stress drop of  0.4 ≤Δσ≤ 10  bar.     

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).     

A wide-angle seismic traverse through the Variscan of southwest Ireland

A wide-angle seismic profile across the western peninsulas of SW Ireland was performed. This region corresponds to the northernmost Variscan thrust and fold deformation. The dense set of 13 shots and 109 stations along the 120  km long profile provides a detailed velocity model of the crust.
  The seismic velocity model, obtained by forward and inverse modelling, defines a five-layer crust. A sedimentary layer, 5–8  km thick, is underlain by an upper-crustal layer of variable thickness, with a base generally at a depth of 10–12  km. Two mid-crustal layers are defined, and a lower-crustal layer below 22  km. The Moho lies at a depth of 30–32  km. A low-velocity zone, which coincides with a well-defined gravity low, is observed in the central part of the region and is modelled as a Caledonian granite which intruded upper-crustal basement. The granite may have acted as a buffer to northward-directed Variscan thrusting. The Dingle–Dungarvan Line (DDL) marks a major change in sedimentary and crustal velocity and structure. It lies immediately to the north of the velocity and gravity low, and shows thickness and velocity differences in many of the underlying crustal layers and even in the Moho. This suggests a deep, pre-Variscan control of the structural development of this area. The model is compatible with thin-skinned tectonics, which terminated at the DDL and which incorporated thrusts involving the sedimentary and upper-crustal layers.     

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.     

State of stress in the Southern Tyrrhenian subduction zone from fault-plane solutions

In this paper we present revised locations and original focal mechanisms computed for intermediate and deep earthquakes that occurred within the Southern Tyrrhenian subduction zone between 1988 and 1994, in order to improve our knowledge of the state of stress for this compressional margin. In particular, we define the stress distribution within a large portion of the descending slab, between 40 and about 450 km depth. The seismicity distribution reveals a continuous 40–50 km thick slab that abruptly increases its dip from subhorizontal in the Ionian Sea to a constant 70° dip in the Tyrrhenian. We computed focal mechanisms for events with magnitudes ranging from 2.7 and 5.7, obtaining the distribution of P - and T -axes for many events for which centroid moment tensor (CMT) solutions are not available, thus enabling the sampling of a larger depth range compared to previous studies. We define three portions of the slab characterized by different distributions of P - and T -axes. A general down-dip compression is found between 165 and 370 km depth, whereas in the upper part of the slab (40–165 km depth) the fault-plane solutions are strongly heterogeneous. Below 370 km the P -axes of the few deep events located further to the north have a shallower dip and are not aligned with the 70° dipping slab, possibly suggesting that they belong to a separated piece of subducted lithosphere. There is a good correspondence between the depth range in which the P -axes plunge closer to the slab dip (∼ 70°) and the interval characterized by the highest seismic energy release (190–370 km).     

Magnetotelluric image of the crust and upper mantle in the backarc of the northwestern Argentinean Andes

Magnetotelluric data from the backarc of the Central Andes in NW Argentinawere re-examined by employing impedance tensor decomposition and 2-D inversion and modelling techniques. The data in the period range of 50–15 000 s were collected on a profile of 220 km length reaching from the Eastern Cordillera across the Santa Barbara System to the Andean foreland of the Argentinean Chaco.
After a dimensionality analysis, data from most sites were treated as regional 2-D. The exception was the eastern section of the profile, where the magnetotelluric transfer functions for periods ≤ 1000 s reflect a 3-D earth. Application of two tensor decomposition schemes yielded a regional strike direction of N–S, which is the azimuth of the Central Andean mountain chains. Several 2-D models were obtained by pseudo- and full 2-D Occam inversion schemes. Special emphasis was placed on the inversion of phase data to reduce the influence of static shifts in the apparent resistivity data. The smooth inversion models all show a good conductor at depth. A final model was then calculated using a finite element forward algorithm.
The most prominent feature of the resulting model is a conductor which rises from depths of 180 km below the Chaco region to 80 km beneath the Santa Barbara System and the Eastern Cordillera. Its interpretation as a rise of the electrical asthenosphere is supported by seismic attenuation studies. Magnetotelluric results, surface heat-flow distribution in the area, and the electrical properties of crustal and mantle rocks suggest that the upper mantle is predominantly ductile beneath the Eastern Cordillera and the western Santa Barbara System. This generally agrees with anelastic seismic attenuation models of the area and is useful in discriminating between models of Q quality factor distribution.     

Shear-wave velocity structure for a spherically averaged earth

Summary. The shear-wave velocity distribution in a spherically averaged Earth is estimated statistically from previously published short-period S travel-time measurements (Uhrhammer). An algorithm is defined for integral inversion techniques which allows estimation of the variance of the velocity distribution from the uncertainties in the S slowness model. Comparisons are made between the resulting S -velocity solution and other solutions in common use. There are significant differences (at the 95 per cent confidence level) between the 5-velocity model determined here and the Jeffreys-Bullen model over the depth ranges of 150–550 km and 2100–2350 km. The 95 per cent confidence level in the present velocity distribution ranges from ± 0.025 km/s at 625 km to ±0.32km/s at 2766 km and averages about ±0.063 or ±1 percent.
Correlations between azimuthally dependent source and station adjustments (which were previously determined (Uhrhammer)) indicate widespread lateral inhomogeneities (up to 3.4 per cent) to depths of approximately 700 km. Up to three-quarters of the source adjustments are due to lateral velocity variations in the source regions. Station adjustments for differential 5 minus P times are significantly correlated with elevation and crustal age, but not with station instrumental magnification.     

Seismotectonics and analysis of earthquakes from the Hindukush region

Summary. Over 80 earthquakes, exclusively from the Hindukush focal region, which were recorded at the Gauribidanur seismic array (GBA) have been used in this study. These events have similar epicentral distances and a narrow azimuthal range from GBA but varying focal depths from 10 to 240 km. A fault plane dipping steeply (75°) in the north-west direction and striking N 66° E has been investigated on the basis of the spatial distribution of earthquakes in two vertical planes through 68° E and 32° N. Short period P -wave recordings up to 30 s were processed using the adaptive cross-correlation filtering technique. Slowness and azimuthal anomalies were obtained for first arrivals. These anomalies show positive as well as negative bias and are attributed to a steep velocity gradient in the upper mantle between the 400–700 km depth range where the seismic rays have their maximum penetration. Relative time residuals between the stations of GBA owe their origin very near to the surface beneath the array. A search of the signals across the array revealed that most of the events occurring at shallower depths had complex signatures as compared to the deeper events. The structure near the source region, complicated source functions and the scattering confined to the crust—upper mantle near source are mainly responsible for the complexity of the Hindukush earthquakes as the transmission zone of the ray tubes from turning point to the recording station is practically the same.     

Source history of the 1905 great Mongolian earthquakes (Tsetserleg, Bolnay)

Two great Mongolian earthquakes, Tsetserleg and Bolnay, occurred on 1905 July 9 and 23. We determined the source history of these events using body waveform inversion. The Tsetserleg rupture (azimuth N60°) correspond to a N60° oriented branch of the long EW oriented Bolnay fault.
Historical seismograms recorded by Wiechert instruments are digitized and corrected for the geometrical deformation due to the recording system. We use predictive filters to recover the signals lost at the minute marks.
The total rupture length for the Tsetserleg earthquake may reach up to 190 km, in order to explain the width of the recorded body waves. This implies adding 60 km to the previously mapped fault. The rupture propagation is mainly eastward. It starts at the southwest of the central subsegment, showing a left lateral strike-slip with a reverse component. The total duration of the modelled source function is 65 s. The seismic moment deduced from the inversion is 10²¹ N m, giving a magnitude   M _w = 8  .
The nucleation of the Bolnay earthquake was at the intersection between the main fault (375 km left lateral strike-slip) and the Teregtiin fault (N160°, 80 km long right lateral strike-slip with a vertical component near the main fault). The rupture was bilateral along the main fault: 100 km to the west and 275 km to east. It also propagated 80 km to the southeast along the Teregtiin fault. The source duration was 115 s. The moment magnitude Mw varies between 8.3 and 8.5.
The nucleation and rupture depths remain uncertain. We tested three cases: (1) nucleation and rupture depth limited to the seismogenic zone; (2) nucleation in the seismogenic zone and rupture propagation going to the base of the crust and (3) nucleation within the crust–upper mantle interface and rupture propagation within the upper mantle.     

The crustal structure beneath the northwest fjords, Iceland, from receiver functions and surface waves

Five broad-band seismic stations were operated in the northwest fjords area of Iceland from 1996 to 1998 as part of the Iceland Hotspot project. The structures of the upper 35  km or so beneath these stations were determined by the modelling and joint inversion of receiver functions and regional surface wave phase velocities. More than 40 teleseismic events and a few regional events containing high-quality surface wave trains were used. Although the middle period passband of the seismograms is corrupted by oceanic microseismic noise, which hinders the interpretation of structural details, the inversions reveal the overall features. Many profiles obtained exhibit large velocity gradients in the upper 5  km or so, smaller zero gradients below this, and, at ~23  km depth, a zone 2–4  km thick with higher velocity gradients. The two shallower intervals are fairly consistent with the 'upper' and 'lower' crust, defined by Flovenz (1980 ). The deep zone of enhanced velocity gradient seems to correspond to the sharp reflector first reported by Bjarnason et al . (1993 ) and identified by them as the 'Moho'. However, this type of structure is not ubiquitous beneath the northwest fjords area. The distinctiveness of the three intervals is variable, and in some cases a structure with velocity gradient increasing smoothly with depth is observed. We term these two end-members structures of the first and second types respectively. Structures of the second type correlate with older areas. Substantial variation in fundamental structure is to be expected in Iceland because of the great geological heterogeneity there.     

Eigendecomposition of the two-point correlation tensor for optimal characterization of mantle convection

The two-point correlation tensor provides complete information on mantle convection accurate up to second-order statistics. Unfortunately, the two-point spatial correlation tensor is in general a data-intensive quantity. In the case of mantle convection, a simplified representation of the two-point spatial correlation tensor can be obtained by using spherical symmetry. The two-point correlation can be expressed in terms of a planar correlation tensor, which reduces the correlation's dependence to only three independent variables: the radial locations of the two points and their angular separation. The eigendecomposition of the planar correlation tensor provides a rational methodology for further representing the second-order statistics contained within the two-point correlation in a compact manner. As an illustration, results on the planar correlation are presented for the thermal anomaly obtained from the tomographic model of Su, Woodward & Dziewonski (1994 ) and the corresponding velocity field obtained from a simple constant-viscosity convection model Zhang & Christensen 1993 ). The first 10 most energetic eigensolutions of the planar correlation, which constitute an almost three orders of magnitude reduction in the data, capture the two-point correlation to 97 per cent accuracy. Furthermore, the energetic eigenfunctions efficiently characterize the thermal and flow structures of the mantle. The signature of the transition zone is clearly evident in the most energetic temperature eigenfunction, which clearly shows a reversal of thermal fluctuations at a depth of around 830  km. In addition, a local peak in the thermal fluctuations can be observed around a depth of 600  km. In contrast, due to the simplicity of the convection model employed, the velocity eigenfunctions exhibit a simple cellular structure that extends over the entire depth of the mantle and do not exhibit transition-zone signatures.     

A test of the precursory accelerating moment release model on some recent New Zealand earthquakes

The proposal that the moment release rate increases in a systematic way in a large region around a forthcoming large earthquake is tested using three recent, large New Zealand events. The three events, 1993–1995, magnitudes 6.7–7.0, occurred in varied tectonic settings. For all three events, a circular precursory region can be found such that the moment release rate of the included seismicity is modelled significantly better by the proposed accelerating model than by a linear moment release model, although in one case the result is dubious. The 'best' such regions have radii from 122 to 167 km, roughly in accord with previous observations world-wide, but are offset by 50–60 km from the associated main shock epicentre. A grid-search procedure is used to test whether these three earthquakes could have been forecast using the accelerating moment release model. For two of the earthquakes the result is positive in terms of location, but the main shock times are only loosely constrained.     

Multimode surface waveform tomography of the Pacific Ocean: a closer look at the lithospheric cooling signature

We present a regional surface waveform tomography of the Pacific upper mantle, obtained using an automated multimode surface waveform inversion technique on fundamental and higher mode Rayleigh waves, to constrain the   V_SV   structure down to ∼400 km depth. We have improved on previous implementations of this technique by robustly accounting for the effects of uncertainties in earthquake source parameters in the tomographic inversion. We have furthermore improved path coverage in the South Pacific region by including Rayleigh wave observations from the French Polynesian Pacific Lithosphere and Upper Mantle Experiment deployment. This improvement has led to imaging of vertical low-velocity structures associated with hotspots within the South Pacific Super-Swell region. We have produced an age-dependent average cross-section for the Pacific Ocean lithosphere and found that the increase in   V_SV   with age is broadly compatible with a half-space cooling model of oceanic lithosphere formation. We cannot confirm evidence for a Pacific-wide reheating event. Our synthetic tests show that detailed interpretation of average   V_SV   trends across the Pacific Ocean may be misleading unless lateral resolution and amplitude recovery are uniform across the region, a condition that is difficult to achieve in such a large oceanic basin with current seismic stations.     

Earthquake processes of the Himalayan collision zone in eastern Nepal and the southern Tibetan Plateau

Focal mechanisms determined from moment tensor inversion and first motion polarities of the Himalayan Nepal Tibet Seismic Experiment (HIMNT) coupled with previously published solutions show the Himalayan continental collision zone near eastern Nepal is deforming by a variety of styles of deformation. These styles include strike-slip, thrust and normal faulting in the upper and lower crust, but mostly strike-slip faulting near or below the crust–mantle boundary (Moho). One normal faulting earthquake from this experiment accommodates east–west extension beneath the Main Himalayan Thrust of the Lesser Himalaya while three upper crustal normal events on the southern Tibetan Plateau are consistent with east–west extension of the Tibetan crust. Strike-slip earthquakes near the Himalayan Moho at depths >60 km also absorb this continental collision. Shallow plunging P -axes and shallow plunging EW trending T -axes, proxies for the predominant strain orientations, show active shearing at focal depths ∼60–90 km beneath the High Himalaya and southern Tibetan Plateau. Beneath the southern Tibetan Plateau the plunge of the P -axes shift from vertical in the upper crust to mostly horizontal near the crust–mantle boundary, indicating that body forces may play larger role at shallower depths than at deeper depths where plate boundary forces may dominate.     

Model of the geoelectrical structure of the mid- and lower mantle in the Europe–Asia region

The conductivity structure of the Earth's mantle was estimated using the induction method down to 2100  km depth for the Europe–Asia region. For this purpose, the responses obtained at seven geomagnetic observatories (IRT, KIV, MOS, NVS, HLP, WIT and NGK) were analysed, together with reliable published results for 11  yr variations. 1-D spherical modelling has shown that, beneath the mid-mantle conductive layer (600–800  km), the conductivity increases slowly from about 1  S  m⁻¹ at 1000  km depth to 10  S  m⁻¹ at 1900  km, while further down (1900–2100  km) this increase is faster. Published models of the lower mantle conductivity obtained using the secular, 30–60  yr variations were also considered, in order to estimate the conductivity at depths down to the core. The new regional model of the lower mantle conductivity does not contradict most modern geoelectrical sounding results. This model supports the idea that the mantle base, situated below 2100  km depth, has a very high conductivity.