Timing of deformation in the Norseman-Wiluna Belt, Yilgarn Craton, Western Australia

Establishing relative and absolute time frameworks for the sedimentary, magmatic, tectonic and gold mineralisation events in the Norseman-Wiluna Belt of the Archean Yilgarn Craton of Western Australia, has long been the main aim of research efforts. Recently published constraints on the timing of sedimentation and absolute granite ages have emphasized the shortcomings of the established rationale used for interpreting the timing of deformation events. In this paper the assumptions underlying this rationale are scrutinized, and it is shown that they are the source of significant misinterpretations. A revised time chart for the deformation events of the belt is established. The first shortening phase to affect the belt, D₁, was preceded by an extensional event D_1e and accompanied by a change from volcanic-dominated to plutonic-dominated magmatism at approximately 2685–2675 Ma. Later extension (D_2e) controlled deposition of the ca 2655 Ma Kurrawang Sequence and was followed by D₂, a major shortening event, which folded this sequence. D₂ must therefore have started after 2655 Ma—at least 20 Ma later than previously thought and after the voluminous 2670–2655 Ma high-Ca granite intrusion. Younger transcurrent deformation, D₃–D₄, waned at around 2630 Ma, suggesting that the crustal shortening deformation cycle D₂–D₄ lasted approximately 20–30 Ma, contemporaneous with low-volume 2650–2630 Ma low-Ca granites and alkaline intrusions. Time constraints on gold deposits suggest a late mineralisation event between 2640–2630 Ma. Thus, D₂–D₄ deformation cycle and late felsic magmatism define a 20–30 Ma long tectonothermal event, which culminated with gold mineralisation. The finding that D₂ folding took place after voluminous high-Ca granite intrusion led to research into the role of competent bodies during folding by means of numerical models. Results suggest that buoyancy-driven doming of pre-tectonic competent bodies trigger growth of antiforms, whereas non-buoyant, competent granite bodies trigger growth of synforms. The conspicuous presence of pre-folding granites in the cores of anticlines may be a result from active buoyancy doming during folding.     

Magnetic anisotropy produced by magma flow: theoretical model and experimental data from Ferrar dolerite sills (Antarctica)

Volcanic rocks forming sills, dykes or lava flows may display a magnetic anisotropy derived from the viscous flow during their emplacement. We model a sill as a steady-state flow of a Bingham fluid, driven by a pressure gradient in a horizontal conduit. The magma velocity as a function of depth is calculated from the motion and constitutive equations. Vorticity and strain rate are determined for a reference system moving with the fluid. The angular velocity and the orientation of an ellipsoidal magnetic grain immersed in the fluid are calculated as functions of time or strain. Magnetic susceptibility is then calculated for a large number of grains with a uniform distribution of initial orientations. It is shown that the magnetic lineation oscillates in the vertical plane through the magma flow direction, and that the magnetic foliation plane changes periodically from horizontal to vertical. The results are compared with the magnetic fabric of Ferrar dolerite sills (Victoria Land, East Antarctica) derived from low-field susceptibility measurements.     

Quantitative imaging of complex structures from dense wide-aperture seismic data by multiscale traveltime and waveform inversions: a case study

An integrated multiscale seismic imaging flow is applied to dense onshore wide‐aperture seismic data recorded in a complex geological setting (thrust belt). An initial P‐wave velocity macromodel is first developed by first‐arrival traveltime tomography. This model is used as an initial guess for subsequent full‐waveform tomography, which leads to greatly improved spatial resolution of the P‐wave velocity model. However, the application of full‐waveform tomography to the high‐frequency part of the source bandwidth is difficult, due to the non‐linearity of this kind of method. Moreover, it is computationally expensive at high frequencies since a finite‐difference method is used to model the wave propagation. Hence, full‐waveform tomography was complemented by asymptotic prestack depth migration to process the full‐source bandwidth and develop a sharp image of the short wavelengths. The final traveltime tomography model and two smoothed versions of the final full‐waveform tomography model were used as a macromodel for the prestack depth migration. In this study, wide‐aperture multifold seismic data are used. After specific preprocessing of the data, 16 frequency components ranging from 5.4 Hz to 20 Hz were inverted in cascade by the full‐waveform tomography algorithm. The full‐waveform tomography successfully imaged SW‐dipping structures previously identified as high‐resistivity bodies. The relevance of the full‐waveform tomography models is demonstrated locally by comparison with a coincident vertical seismic profiling (VSP) log available on the profile. The prestack depth‐migrated images, inferred from the traveltime, and the smoothed full‐waveform tomography macromodels are shown to be, on the whole, consistent with the final full‐waveform tomography model. A more detailed analysis, based on common‐image gather computations, and local comparison with the VSP log revealed that the most accurate migrated sections are those obtained from the full‐waveform tomography macromodels. A resolution analysis suggests that the asymptotic prestack depth migration successfully migrated the wide‐aperture components of the data, allowing medium wavelengths in addition to the short wavelengths of the structure to be imaged. The processing flow that we applied to dense wide‐aperture seismic data is shown to provide a promising approach, complementary to more classical seismic reflection data processing, to quantitative imaging of complex geological structures.     

Seismic studies of the Brasília fold belt at the western border of the São Francisco Craton, Central Brazil, using receiver function, surface-wave dispersion and teleseismic tomography

The Tocantins Province in Central Brazil is composed of a series of SSW–NNE trending terranes of mainly Proterozoic ages, which stabilized in the Neoproterozoic in the final collision between the Amazon and São Francisco cratons. No previous information on crustal seismic properties was available for this region. Several broadband stations were used to study the regional patterns of crustal and upper mantle structure, extending the results of a recent E–W seismic refraction profile. Receiver functions and surface wave dispersion showed a thin crust (33–37 km) in the Neoproterozoic Magmatic Arc terrane. High average crustal Vp/Vs ratios (1.74–1.76) were consistently observed in this unit. The foreland domain of the Brasília foldbelt, on the other hand, is characterized by thicker crust (42–43 km). Low Vp/Vs ratios (1.70–1.72) were observed in the low-grade foreland fold and thrust zone of the Brasília belt adjacent to the São Francisco craton. Teleseismic P-wave tomography shows that the lithospheric upper mantle has lower velocities beneath the Magmatic Arc and Goiás Massif compared with the foreland zone of the belt and São Francisco craton. The variations in crustal thickness and upper mantle velocities observed with the broadband stations correlate well with the measurements along the seismic refraction profile. The integration of all seismic observations and gravity data indicates a strong lithospheric contrast between the Goiás Massif and the foreland domain of the Brasília belt, whereas little variation was found across the foldbelt/craton surface boundary. These results support the hypothesis that the Brasília foreland domain and the São Francisco craton were part of a larger São Francisco-Congo continental plate in the final collision with the Amazon plate.     

Historical measurements of the Earth's magnetic field compared with remanence directions from lava flows in Italy over the last four centuries

Direct measurements of the Earth's magnetic field in Italy since 1640 a.d. have been used to check the remanence directions derived from historically dated volcanic rocks of Etna and Vesuvius. Direct measurements consist of the records of L’Aquila and Pola geomagnetic observatories, the repeat stations of the Italian Magnetic Network and the data base of the Historical Italian Geomagnetic Data Catalogue. All have been relocated to the same reference site (Viterbo — lat. 42.45°N, long. 12.03°E) in order to draw a reference secular variation (SV) curve. The direction of the Earth's field at Viterbo has also been calculated from the historical records (2000-1600) of ref. [Jackson, A., Jonkers, A.R.T., Walker, M.R., 2000. Four centuries of geomagnetic secular variation from historical records. Phil. Trans. R. Soc. London, Ser. A 358, 957-990] database. The remanence directions from Etna show a general agreement with the trend of the SV curve, although their inclination is usually lower than that from the direct measurement. The directions from Vesuvius are more scattered. Large discrepancies occur at both volcanoes and in some cases have been ascribed in the literature to poor geographic information, making it difficult to identify the flows actually emplaced during the eruptions reported in the chronicles. Closer examination shows that the great majority of the best-defined remanence directions (semi-angle of confidence α₉₅ < 2.5°) deviate significantly from the geomagnetic direction measured at the time of the emplacement, the angle between the two directions being larger than the α₉₅ value. The value of 2.5-3.0° can thus be regarded as a conservative evaluation of the error when dealing with dating Etna and Vesuvius lava flows older than 17th century, even when the accuracy attained in remanence measurements is higher. In default of a SV curve for Italy derived from archaeological artefacts, a further error in dating is introduced when reference is made to SV curves of other countries, even if well-established, as these are from regions too far from Italy (>600 km) to confidently relocate magnetic directions.     

A complex, Subplinian-type eruption from low-viscosity, phonolitic to tephri-phonolitic magma: the AD 472 (Pollena) eruption of Somma-Vesuvius, Italy

The combined use of field investigation and laboratory analyses allowed the detailed stratigraphic reconstruction of the Pollena eruption (472 AD) of Somma-Vesuvius. Three main eruptive phases were recognized, related either to changes in the eruptive processes and/or to relative changes of melt composition. The eruption shows a pulsating behavior with deposition of pyroclastic fall beds and generation of dilute and dense pyroclastic density currents (PDC). The eruptive mechanisms and transportation dynamics were reconstructed for the whole eruption. Column heights were between 12 and 20 km, corresponding to mass discharge rates (MDR) of 7×10⁶ kg/s and 3.4×10⁷ kg/s. Eruptive dynamics were driven by magmatic fragmentation of a phono-tephritic to tephri-phonolitic magma during Phases I and II, whereas phreatomagmatic fragmentation dominated Phase III. Magma composition varies between phonolitic and tephritic-phonolitic, with melt viscosity likely not in excess of 10³ Pa s. The volume of the pyroclastic fall deposits, calculated by using of proximal isopachs, is 0.44 km³. This increases to 1.38 km³ if ash volumes are extrapolated on a log thickness vs. square root area diagram using one distal isopach and column height.Editorial responsibility: R Cioni