Subduction beneath western South America: evidence from converted phases

Summary. There are two conflicting models for the dip of the subduction zone beneath central Peru and beneath central Chile. One model, based primarily on the distribution of hypocentres thought to be most reliably located, postulates a shallow (∼ 10°) dip for the downgoing plate. The second model, re-examined in this paper and based chiefly on an analysis of ScSp converted phases, postulates a normal (∼ 30°) dip for the subduction.
A detailed examination of ScSp data for central Peru and central Chile shows that ScSp arrivals can be identified on all seismograms on which the predicted signal/noise ratio is greater than 1, and the measured amplitude ratios of ScSp to ScS imply a normal dip for the conversion interfaces. The characteristics of the ScSp arrivals in these regions are virtually identical to those for ScSp arrivals in southern Peru, where a well-defined Benioff zone and the calculated ScSp -conversion region both imply a normal dip for the subduction. Hence ScSp observations in South America support the model of an approximately constant dip of subduction from central Peru to central Chile to a depth of at least 120 km, and impose a constraint that must be satisfied by any satisfactory tectonic model for these regions.     

New insight into the crust and upper mantle structure under Alaska

To better understand the seismic structure of the subducting Pacific plate under Alaska, we determined the three-dimensional P-wave velocity structure to a depth of approximately 200 km beneath Alaska using 438,146 P-wave arrival times from 10,900 earthquakes. In this study an irregular grid parameterization was adopted to express the velocity structure under Alaska. The number of grid nodes increases from north to south in the study area so that the spacing between grid nodes is approximately the same in the longitude direction. Our results suggest that the subducting Pacific slab under Alaska can be divided into three different parts based on its geometry and velocity structure. The western part has features similar to those in other subduction zones. In the central part a thick low-velocity zone is imaged at the top of the subducting Pacific slab beneath north of the Kenai Peninsula, which is believed to be most likely the oceanic crust plus an overlying serpentinized zone and the coupled Yakutat terrane subducted with the Pacific slab. In the eastern part, significant high-velocity anomalies are visible to 60–90 km depth, suggesting that the Pacific slab has only subducted down to that depth.     

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.     

Geodetically determined strain across the southern end of the Tonga-Kermadec-Hikurangi subduction zone

Summary. A method of simultaneous reduction is presented for determining strain rates from multiple triangulation surveys where common triangulation stations have been used, but the angles of the old survey have not necessarily been reobserved. This method is applied to triangulation in the northern South Island, at the southern end of the Tonga-Kermadec-Hikurangi subduction zone. From a profile of shear strain across the Indian-Pacific plate boundary, the displacement of the Indian plate relative to the Pacific is calculated to be 54 ± 9 mm yr⁻¹ at an azimuth of 84°± 10°, in remarkable agreement with the motion predicted by global plate tectonic models. Most of this motion occurs within a 150 km wide zone bounded on the east by the Hikurangi Trough. Within this zone the motion is partitioned: near the Hikurangi Trough no slip is occurring at the upper surface of the subducting Pacific plate (the subduction thrust) and motion is predominantly thrusting normal to the trough axis: to the west is a region of predominantly dextral strike slip faulting. This pattern is consistent with Fitch's model of oblique subduction. To the south of the profile, a change is observed in the azimuth of the faulting along a line which marks the southern extent of the subduction slab, indicating the end of the partitioned motion.     

Strain accumulation along an oblique plate boundary: the Reykjanes Peninsula, southwest Iceland

We use annual GPS observations on the Reykjanes Peninsula (RP) from 2000 to 2006 to generate maps of surface velocities and strain rates across the active plate boundary. We find that the surface deformation on the RP is consistent with oblique plate boundary motion on a regional scale, although considerable temporal and spatial strain rate variations are observed within the plate boundary zone. A small, but consistent increase in eastward velocity is observed at several stations on the southern part of the peninsula, compared to the 1993–1998 time period. The 2000–2006 velocities can be modelled by approximating the plate boundary as a series of vertical dislocations with left-lateral motion and opening. For the RP plate boundary we estimate left-lateral motion  18⁺⁴₋₃ mm yr⁻¹  and opening of  7⁺³₋₂ mm yr⁻¹  below a locking depth of  7⁺¹₋₂ km  . The resulting deep motion of  20⁺⁴₋₃ mm yr⁻¹  in the direction of  N(100⁺⁸₋₆)°E  agrees well with the predicted relative North America–Eurasia rate. We calculate the areal and shear strain rates using velocities from two periods: 1993–1998 and 2000–2006. The deep motion along the plate boundary results in left-lateral shear strain rates, which are perturbed by shallow deformation due to the 1994–1998 inflation and elevated seismicity in the Hengill–Hrómundartindur volcanic system, geothermal fluid extraction at the Svartsengi power plant, and possibly earthquake activity on the central part of the peninsula.     

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

Absolute S-velocity estimation from receiver functions

We present a novel method to recover absolute S velocities from receiver functions.
For a homogeneous half-space the S velocity can be calculated from the horizontal slowness and the angle of surface particle motion for an incident P wave. Generally, the calculated S velocity is an apparent half-space value which depends on model inhomogeneity and P -waveform. We therefore, suggest to calculate such apparent half-space S velocities from low-pass filtered (smoothed) receiver functions using a suite of filter-parameters, T . The use of receiver functions neutralize the influence of the P -waveform, and the successive low-pass filterings emphasize the variation of S velocity with depth.
We apply this   V _{S ,app.}( T )  technique to teleseismic data from three stations: FUR, BFO and SUM, situated on thick sediments, bedrock and the Greenland ice cap, respectively. The observed   V _{S ,app.}( T )  curves indicate the absolute S velocities from the near surface to the uppermost mantle beneath each station, clearly revealing the different geological environments. Application of linearized, iterative inversion quantify these observations into   V _S ( z )  models, practically independent of the S -velocity starting model. The obtained models show high consistency with independent geoscientific results. These cases provide also a general validation of the   V _{S ,app.}( T )  method.
We propose the computation of   V _{S ,app.}( T )  curves for individual three-component broad-band stations, both for direct indication of the S velocities and for inverse modelling.     

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 reflection measurements behind the Hikurangi convergent margin, southern North Island, New Zealand.

Summary. The active Australian-Pacific plate boundary passes through New Zealand. In the north, the Pacific plate subducts beneath the Australian plate with an accretionary wedge forming the eastern continental (Hikurangi) margin of the North Island. The structure of the region behind the Hikurangi margin changes from the extensional back-arc basin under central North Island to a postulated crustal downwarp under the southern North Island. A 100 km long multichannel seismic reflection profile was recorded across the region of crustal downwarp. The data show discontinuous coherent reflectors dipping westwards at the east end of the profile, and east dipping reflectors at the west end, from depths of 9 to 15 s two way time. Simple hand migration of these events indicate that the east dipping reflectors, interpreted as the base of the Australian plate crust, abut against the west dipping reflectors which are interpreted as marking the top of the subducted Pacific plate. Detailed earthquake hypocentre locations in the area show a dipping zone of high seismicity, the top of which coincides closely with the west dipping events, thus supporting this interpretation.     

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.     

Travel times of Pn -waves in the Aegean and surrounding area

Summary. Based on accurately located 23 very shallow earthquakes ( h = 1–14 km) in northern and central Greece by portable networks of seismic stations and by the joint epicentre method, the travel times of the P_n -waves from the foci of these earthquakes to the sites of 54 permanent stations in the Balkan region have been determined. The travel times of P_n -waves in the central and eastern part of the area (eastern Greece, south-eastern Yugoslavia, the Aegean Sea, Bulgaria, southern Romania, western Turkey) fit a straight line very well with the P_n velocity equal to 7.9 ± 0.1 km s^-1. On the contrary, the travel times of P_n -waves to stations in the western part of the area (Albania, western Greece) do not fit this curve because the P_n -waves travelling to these stations are delayed by more than 1 s due to the thicker crust under the Dinarides–Hellenides mountain range. Time delays for P_n -waves have been calculated for each permanent station in the Balkan area with respect to the mean travel-time curve of these waves in the central and eastern part of the area. Corrections of the travel times for these delays contribute very much to the improvement of the accuracy in the location of the shallow earthquakes in the Aegean and surrounding area.     

The dammed Hikurangi Trough: a channel-fed trench blocked by subducting seamounts and their wake avalanches (New Zealand–France GeodyNZ Project)

The Hikurangi Trough, off eastern New Zealand, is at the southern end of the Tonga–Kermadec–Hikurangi subduction system, which merges into a zone of intracontinental transform. The trough is mainly a turbidite-filled structural trench but includes an oblique-collision, foredeep basin. Its northern end has a sharp boundary with the deep, sediment-starved, Kermadec Trench. Swath-mapping, sampling and seismic surveys show modern sediment input is mainly via Kaikoura Canyon, which intercepts littoral drift at the southern, intracontinental apex of the trough, with minor input from seep gullies. Glacial age input was via many canyons and about an order of magnitude greater. Beyond a narrow, gravelly, intracontinental foredeep, the southern trench-basin is characterized by a channel meandering around the seaward edge of mainly Plio-Pleistocene, overbank deposits that reach 5 km in thickness. The aggrading channel has sandy turbidites, but low-backscatter, and long-wavelength bedforms indicating thick flows. Levées on both sides are capped by tangentially aligned mudwaves on the outsides of bends, indicating centrifugal overflow from heads of dense, fast-moving, autosuspension flows. The higher, left-bank levée also has levée-parallel mudwaves, indicating Coriolis and/or boundary currents effects on dilute flows or tail plumes. In the northern trough, basin-fill is generally less than 2 km thick and includes widespread overbank turbidites, a massive, blocky, avalanche deposit and an extensive, buried, debris flow deposit. A line of low seamounts on the subducting plate acts as a dam preventing modern turbidity currents from reaching the Kermadec Trench. Major margin collapse probably occurred in the wake of a large subducting seamount; this seamount and its wake debris flow probably dammed the trench from 2 Ma to 0.5 Ma. Before this, similar dams may have re-routed turbidity currents across the plateau.     

Reflection mapping across the convergent margin of western Canada

Summary. Five marine multichannel seismic reflection profiles totalling 520 km were recorded across the western Canada convergent margin where the Juan de Fuca plate is subducting beneath North America. The data extend the results of LITHOPROBE on Vancouver Island. The primary objectives are definition of the offshore accretionary structures and clarification of the convergent interaction between the two plates. The main features of this preliminary interpretation are: (1) the subduction deformation front is complex with evidence of sediments being accreted and subducted; (2) the top of the oceanic crust and the Mono are imaged below the deep water sedimentary basin; (3) the top of the subducting plate is clearly imaged below the shelf; (4) beneath the inner shelf, one band of high reflectivity underlain by a zone of lesser reflectivity lies above the plate; (5) alternative interpretations place the present zone of decoupling at the base of the reflective band or the top of the plate; (6) the San Juan and Leech River faults that bound small accreted terranes are imaged as thrusts that merge at depth.     

Palaeomagnetism of the continental sector of the Cameroon Volcanic Line, West Africa

Measurement of samples from 154 sites in the continental sector of the Cameroon Volcanic Line yielded six palaeomagnetic poles, at 243.6°E, 84.6°N, α ₉₅ = 6.8°; 224.3°E, 81.2°N, α ₉₅ = 8.4°; 176.1°E, 82.0°N, α ₉₅ = 8.5°; 164.3°E, 86.4°N, α ₉₅ = 3.4°; 169.4°E, 82.6°N, α ₉₅ = 4.6° and 174.7°E, 72.8°N, α ₉₅ = 9.5°, belonging to rocks which have been dated by the K–Ar method at 0.4–0.9  Ma, 2.6  Ma, 6.5–11  Ma, 12–17  Ma, 20–24  Ma and 28–31  Ma, respectively. The results are in general agreement with other palaeomagnetic poles from Oligocene to Recent formations in Africa.
  The first three poles for rocks formed between 0.4 and 11  Ma are not significantly different from the present geographical pole. Together with other African poles for the same period, this suggests that the African continent has moved very little relative to the pole since 11  Ma. The other three poles for rocks dated between 12 and 31  Ma are significantly different from the present geographical pole, showing a 5° polar deviation from the present pole in the Miocene and 13° in the Middle Oligocene.     

Early Triassic palaeomagnetic results from the Ryeongnam Block, Korean Peninsula: the eastern extension of the North China Block

Greenish sandstones in the Early Triassic Nogam Formation of the Ryeongnam Block, Korean Peninsula were collected at 23 sites for palaeomagnetic study. A high-temperature magnetization component with unblocking temperatures of 670–690 °C was isolated from seven sites and yielded a positive fold test at the 95 per cent confidence level. The high-temperature component is interpreted to be of primary origin because the folding age is Middle Triassic. The Early Triassic palaeomagnetic direction for the Ryeongnam Block after tilt correction is D =347.1°, I =23.8° ( α ₉₅=5.5°). The palaeomagnetic pole (62.5°N, 336.8°E, A ₉₅ = 4.7°) shows good agreement with the coeval pole for the North China Block, suggesting that the Ryeongnam Block has been part of the North China Block at least since Early Triassic times. A tectonic history of the Korean Peninsula includes obduction of the eastern part of the South China Block onto the central part of the Korean Peninsula in the Permian, with the Ryeongnam Block geographically isolated from the main part of the North China Block. Collision of the North and South China blocks commenced initially at the Korean Peninsula, and suturing of the two blocks progressed westwards.     

Double-planed deep seismic zone and upper-mantle structure in the Northeastern Japan Arc

Summary. The ScSp wave converted from the ScS wave at the boundary between the descending lithospheric slab and the mantle above it was clearly observed from a nearby deep earthquake with magnitude 7.7 at some stations of the seismic network of Tohoku University which covers the Tohoku District, the northeastern part of Honshu, Japan. By applying the three-dimensional seismic-ray tracing method, the location of this boundary was determined from the difference in arrival time between the ScS and ScSp waves. The result shows that the upper boundary of the descending slab lies exactly on the upper plane of the double-planed deep seismic zone found in the Northeastern Japan Arc.
There is an additional evidence that the boundary is located on the upper plane of the double-planed deep seismic zone. The hypocentre distribution of intermediate-depth earthquakes located by the small-scale seismic-array observation is extremely different from that obtained by the relatively large-scale seismic network. The discrepancy in the distribution of hypocentres of the same earthquake independently located is well explained by the inclined lithospheric slab model derived from the difference in arrival time between the ScS and ScSp waves.
The earthquakes with reverse faulting or with down-dip compressional stresses occur at the upper boundary of the descending slab. Within the descending slab, the earthquakes with down-dip extensional stresses also occur in a very narrow zone from 30 to 40 km below the dipping boundary in the depth range from 50 to about 200 km, and these shocks form the lower plane of the double-planed deep seismic zone.     

Run-up from impact tsunami

We report on calculations of the on-shore run-up of waves that might be generated by the impact of subkilometre asteroids into the deep ocean. The calculations were done with the COULWAVE code, which models the propagation and shore-interaction of non-linear moderate- to long-wavelength waves  ( kh < π)  using the extended Boussinesq approximation. We carried out run-up calculations for several different situations: (1) laboratory-scale monochromatic wave trains onto simple slopes; (2) 10–100 m monochromatic wave trains onto simple slopes; (3) 10–100 m monochromatic wave trains onto a compound slope representing a typical bathymetric profile of the Pacific coast of North America; (4) time-variable scaled trains generated by the collapse of an impact cavity in deep water onto simple slopes and (5) full-amplitude trains onto the Pacific coast profile. For the last case, we also investigated the effects of bottom friction on the run-up. For all cases, we compare our results with the so-called 'Irribaren scaling': The relative run-up   R / H ₀=ξ= s ( H ₀/ L ₀)^−1/2  , where the run-up is   R , H ₀  is the deep-water waveheight, L ₀ is the deep-water wavelength, s is the slope and ξ is a dimensionless quantity known as the Irribaren number. Our results suggest that Irribaren scaling breaks down for shallow slopes   s ≤ 0.01  when  ξ < 0.1 − 0.2  , below which   R / H ₀  is approximately constant. This regime corresponds to steep waves and very shallow slopes, which are the most relevant for impact tsunami, but also the most difficult to access experimentally.     

Post-seismic relaxation following the great 2004 Sumatra-Andaman earthquake on a compressible self-gravitating Earth

The   M _w γ 9.0  2004 December 26 Sumatra-Andaman and   M _w = 8.7  2005 March 28 Nias earthquakes, which collectively ruptured approximately 1800 km of the Andaman and Sunda subduction zones, are expected to be followed by vigorous viscoelastic relaxation involving both the upper and lower mantle. Because of these large spatial dimensions it is desirable to fully account for gravitational coupling effects in the relaxation process. We present a stable method of computing relaxation of a spherically-stratified, compressible and self-gravitating viscoelastic Earth following an impulsive moment release event. The solution is cast in terms of a spherical harmonic expansion of viscoelastic normal modes. For simple layered viscoelastic models, which include a low-viscosity oceanic asthenosphere, we predict substantial post-seismic effects over a region several 100s of km wide surrounding the eastern Indian Ocean. We compare observed GPS time-series from ten regional sites (mostly in Thailand and Indonesia), beginning in 2004 December, with synthetic time-series that include the coseismic and post-seismic effects of the 2004 December 26 and 2005 March 28 earthquakes. A viscosity structure involving a biviscous (Burgers body) rheology in the asthenosphere explains the pattern and amplitude of post-seismic offsets remarkably well.     

The transition from subduction to continental collision: crustal structure in the North Canterbury region, New Zealand

The North Canterbury region marks the transition from Pacific plate subduction to continental collision in the South Island of New Zealand. Details of the seismicity, structure and tectonics of this region have been revealed by an 11-week microearthquake survey using 24 portable digital seismographs. Arrival time data from a well-recorded subset of microearthquakes have been combined with those from three explosions at the corners of the microearthquake network in a simultaneous inversion for both hypocentres and velocity structure. The velocity structure is consistent with the crust in North Canterbury being an extension of the converging Chatham Rise. The crust is about 27 km thick, and consists of an 11 km thick seismic upper crust and 7 km thick seismic lower crust, with the middle part of the crust being relatively aseismic. Seismic velocities are consistent with the upper and middle crust being composed of greywacke and schist respectively, while several lines of evidence suggest that the lower crust is the lower part of the old oceanic crust on which the overlying rocks were originally deposited.
The distribution of relocated earthquakes deeper than 15 km indicates that the seismic lower crust changes dip markedly near 43S. To the south-west it is subhorizontal, while to the north-east it dips north-west at about 10. Fault-plane solutions for these earthquakes also change near 43S. For events to the south, P -axes trend approximately normal to the plate boundary (reflecting continental collision), while for events to the north, T -axes are aligned down the dip of the subducted plate (reflecting slab pull). While lithospheric subduction is continuous across the transition, it is not clear whether the lower crust near 43S is flexed or torn.     

Thermochronometry and palaeomagnetism of the Archaean Nelshoogte Pluton, South Africa

²⁰⁷Pb/²⁰⁶Pb single-grain zircon, ⁴⁰Ar/³⁹Ar single-grain hornblende and biotite, and ⁴⁰Ar/³⁹Ar bulk-sample muscovite and biotite ages from the Nelshoogte trondhjemite pluton located in eastern Transvaal, South Africa, show that this granitoid had a protracted thermal history spanning 3213±4  Ma to about 3000  Ma. Whole-rock ⁴⁰Ar/³⁹Ar ages from cross-cutting dolerite dykes indicate that these were intruded at about 1900  Ma. There is no evidence of this or other, later events significantly affecting the argon systematics of the minerals from the pluton dated by the ⁴⁰Ar/³⁹Ar method.
  The pluton has a well-defined palaeomagnetic pole which is dated at 3179±18 (2 σ ) Ma by ⁴⁰Ar/³⁹Ar dating of hornblende. This pole (18°N, 310°E, A ₉₅=9°) yields a palaeolatitude of 0°, significantly different from other Archaean poles from the Kaapvaal Craton. The palaeolatitude difference implies that there was significant apparent polar wander during the Archaean. A second, overprinting magnetization seen in the pluton is also seen in the lower-Proterozoic dolerite dykes, and is consistent with other lower-Proterozoic (2150–1950  Ma) poles for southern Africa.