Geophysical imaging of buried volcanic structures within a continental back-arc basin: The Central Volcanic Region, North Island, New Zealand

Hidden beneath the ~ 2 km thick low-velocity volcaniclastics on the western margin of the Central Volcanic Region, North Island, New Zealand, are two structures that represent the early history of volcanic activity in a continental back-arc. These ~ 20 × 20 km structures, at Tokoroa and Mangakino, form an adjacent gravity high and low, respectively. Interpretations from seismic refraction arrivals and gravity modelling indicate the − 65 mgal Mangakino residual gravity anomaly can be modelled, in part, by two low-density bodies that reach depths of ~ 6.5 km, whereas the Tokoroa gravity anomaly is due to a higher density rock coming, at most, to within ~ 650 m of the surface. The Mangakino anomaly is interpreted to be due to the remnants of magma chambers that fed large ignimbrite eruptions from about 1.2 Ma. An andesite volcano or complex volcanic structure is the preferred interpretation for the Tokoroa gravity high. The size of the putative volcanic structure is comparable to the presently active Tongariro Volcanic Complex in the centre of North Island.     

Progressive cooling of the hyaloclastite ridge at Gjálp, Iceland, 1996–2005

In the subglacial eruption at Gjálp in October 1996 a 6 km long and 500 m high subglacial hyaloclastite ridge was formed while large volumes of ice were melted by extremely fast heat transfer from magma to ice. Repeated surveying of ice surface geometry, measurement of inflow of ice, and a full Stokes 2-D ice flow model have been combined to estimate the heat output from Gjálp for the period 1996–2005. The very high heat output of order 10⁶ MW during the eruption was followed by rapid decline, dropping to  2500 MW by mid 1997. It remained similar until mid 1999 but declined to 700 MW in 1999–2001. Since 2001 heat output has been insignificant, probably of order 10 MW. The total heat carried with the 1.2 × 10¹² kg of basaltic andesite erupted (0.45 km³ DRE) is estimated to have been 1.5 × 10¹⁸ J. About two thirds of the thermal energy released from the 0.7 km³ edifice in Gjálp occurred during the 13-day long eruption, 20% was released from end of eruption until mid 1997, a further 10% in 1997–2001, and from mid 2001 to present, only a small fraction remained. The post-eruption heat output history can be reconciled with the gradual release of 5 × 10¹⁷ J thermal energy remaining in the Gjálp ridge after the eruption, assuming single phase liquid convection in the cooling edifice. The average temperature of the edifice is found to have been approximately 240 °C at the end of the eruption, dropping to  110 °C after 9 months and reaching  40 °C in 2001. Although an initial period of several months of very high permeability is possible, the most probable value of the permeability from 1997 onwards is of order 10^− 12 m². This is consistent with consolidated/palagonitized hyaloclastite but incompatible with unconsolidated tephra. This may indicate that palagonitization had advanced sufficiently in the first 1–2 years to form a consolidated hyaloclastite ridge, resistant to erosion. No ice flow traversing the Gjálp ridge has been observed, suggesting that it has effectively been shielded from glacial erosion in its first 10 years of existence.     

Eocene melting of subducting continental crust and early uplifting of central Tibet: Evidence from central-western Qiangtang high-K calc-alkaline andesites, dacites and rhyolites

Changes in oceanic O–Sr isotopic compositions and global cooling beginning in the Eocene are considered to have been caused by the uplift of the Tibetan Plateau. The specific timing and uplift mechanism, however, have long been subjects of debate. We investigated the Duogecuoren lavas of the central-western Qiangtang Block, which form the largest outcrops among Cenozoic lavas in northern-central Tibet and have widely been considered as shoshonitic. Our study demonstrates, however, that most of these lavas are high-K calc-alkaline andesites, dacites and rhyolites. Moreover, they are characterized by high Sr (367–2472 ppm) and Al₂O₃ (14.55–16.86 wt.%) and low Y (3.05–16.9 ppm) and Yb (0.31–1.48 ppm) contents and high La/Yb (27–100) and Sr/Y (48–240) ratios, similar to adakitic rocks derived by partial melting of an eclogitic source. They can be further classified as either peraluminous and metaluminous subtypes. The peraluminous rocks have relatively high SiO₂ (> 66 wt.%) contents, and low MgO (< 1.0 wt.%), Cr (4.94–23.3 ppm) and Ni (2.33–17.0 ppm) contents and Mg^# (20–50) values, while the metaluminous rocks exhibit relatively low SiO₂ (55–69 wt.%) contents, and high MgO (1.41–6.34), Cr (25.7–383 ppm), Ni (14.13–183 ppm) and Mg^# (46–69) values, similar to magnesian andesites. ⁴⁰Ar/³⁹Ar and SHRIMP zircon U–Pb dating reveal that both peraluminous and metaluminous adakitic rocks erupted in the Eocene (46–38 Ma). Paleocene–Early Miocene thrust faults and associated syn-contractional basin deposits in the Qiangtang Block suggest that this region was undergoing crustal shortening within a continent during the Eocene. The low ε_Nd (− 2.81 to − 6.91) and high ⁸⁷Sr/⁸⁶Sr (0.7057–0.7097), Th (11.2–32.3 ppm) and Th/La (0.23–0.88) values in the Duogecuoren adakitic rocks further indicate that they were not derived by partial melting of subducted oceanic crust. Taking into account tectonic and geophysical data and the compositions of xenoliths in Cenozoic lava in northern-central Tibet, we suggest that the peraluminous adakitic rocks were most probably derived by partial melting of subducted sediment-dominated continent of the Songpan-Ganzi Block along the Jinsha suture to the north at a relatively shallow position (the hornblende + garnet stability field), but the metaluminous adakitic rocks likely originated from the interaction between peraluminous adakitic melts generated at greater depths (the garnet + rutile stability field) and mantle. Because the Duogecuoren adakitic rocks must have originated from a garnet-bearing (namely, eclogite facies) source, Eocene continental subduction along the Jinsha suture caused the thickening of the Qiangtang crust. Given that crustal thickening generally equates with elevation, the uplift of the Central Tibetan Plateau probably began as early as 45–38 Ma, which provides important evidence for tectonically driven models of oceanic O–Sr isotope evolution during global cooling and Asian continental aridification beginning in the Eocene.     

Anomalous subsidence on the rifted volcanic margin of Pakistan: No influence from Deccan plume

The role of hotter than ambient plume mantle in the formation of a rifted volcanic margin in the northern Arabian Sea is investigated using subsidence analysis of a drill site located on the seismically defined Somnath volcanic ridge. The ridge has experienced > 4 km of subsidence since 65 Ma and lies within oceanic lithosphere. We estimate crustal thickness to be 9.5–11.5 km. Curiously < 400 m of the thermal subsidence occurred prior to 37 Ma, when subsidence rates would normally be at a maximum. We reject the hypothesis that this was caused by increasing plume dynamic support after continental break-up because the size of the thermal anomalies required are unrealistic (> 600 °C), especially considering the rapid northward drift of India relative to the Deccan-Réunion hotspot. We suggest that this reflects very slow lithospheric growth, possibly caused by vigorous asthenospheric convection lasting > 28 m.y., and induced by the steep continent–ocean boundary. Post-rift slow subsidence is also recognized on volcanic margins in the NE Atlantic and SE Newfoundland and cannot be used as a unique indicator of plume mantle involvement in continental break-up.     

Structure of Réunion Island (Indian Ocean) inferred from the interpretation of gravity anomalies

Réunion is a volcanic edifice whose origin is related to a hot spot in the Indian Ocean. Only 3% of its volume is emergent. Many geological and geophysical studies were carried out on Réunion Island during the 1980's but few of them allow study of the internal structure of the edifice. Several gravity surveys have been carried out on the island since 1976 and we have compiled the available data set. The lack of data on the western side of the island led us to conduct a regional survey in 1993 to obtain a more homogeneous distribution of the stations. Computation of Bouguer anomalies for different correction densities accounts for the variable density of the rocks constituting the edifice and provides a distribution of gravity anomalies interpreted as dense bodies of intrusive rocks inside the edifice. Two very large intrusive complexes can be unambiguously recognised: one beneath Piton des Neiges and one beneath the Grand Brûlé area. Both have been penetrated by geothermal exploration drill holes and the first is also known from outcrop observations. 2.5D simple models were constructed to reveal the geometry and extent of the buried intrusives. They are deeply rooted, extending several kilometres below sea level, and extensive (20–25 km long and 10–13 km wide for the Piton des Neiges complex, 12–15 km long and some kilometres wide for the Grand Brûlé complex). The development of such complexes implies that the activity of the two volcanic centres was long lasting and remained stable while the volcanoes were growing. The Grand Brûlé complex has been interpreted as relics of an old volcano named Alizés Volcano. The interpretation of the gravity maps suggests the presence of a ridge of dense rocks to the North of the axis joining the centres of Piton des Neiges and Piton de la Fournaise volcanoes. By analogy with the other structures, 2.5D models show that this structure would culminate between 0 and 1 km below sea level and be 15 km wide. This complex induces a maximum anomaly in Takamaka Valley and we thus propose to name it Takamaka Volcano. No geological evidence of the nature of these dense rocks is available but the ridge coincides with structures revealed by magnetic and seismic data. Interpretation of the Bouguer anomaly maps suggests that the inner gravity structure of Piton de la Fournaise is not characterised by the presence of a voluminous dense body but probably by more restricted concentrations of dense rocks. Some structures can be recognised: along the present NE and SE rift zones and in the previous central part of Piton de la Fournaise to the West of the present summit. The recent eastward migration of the centre of activity of Piton de la Fournaise accounts for the lack of a large positive anomaly beneath the active craters.     

Seismic evidence for a metastable olivine wedge in the subducting Pacific slab under Japan Sea

We apply a forward-modeling approach to high-quality arrival time data from 23 deep earthquakes greater than 400 km depth to investigate the detailed structure of the subducting Pacific slab beneath the Japan Sea. Our results show that a finger-like anomaly exists within the subducting Pacific slab below 400 km depth, which has a P-wave velocity 5% lower than the surrounding slab velocity (or 3% lower than that of the normal mantle), suggesting the existence of a metastable olivine wedge (MOW) in the slab. The MOW top and bottom depths are 400 and 560 km, respectively. The MOW is estimated to be about 50 km wide at 400 km depth and close to the slab upper boundary. At 560 km depth the MOW is located at about 25 km below the slab upper boundary. Most of the deep earthquakes are located in the MOW. Our results favor transformational faulting as the mechanism for deep earthquakes.     

Crustal structure and origin of the Cape Verde Rise

The Cape Verde Islands are located on a mid-plate topographic swell and are thought to have formed above a deep mantle plume. Wide-angle seismic data have been used to determine the crustal and uppermost mantle structure along a ~ 440 km long transect of the archipelago. Modelling shows that ‘normal’ oceanic crust, ~ 7 km in thickness, exists between the islands and is gently flexed due to volcano loading. There is no direct evidence for high density bodies in the lower crust or for an anomalously low density upper mantle. The observed flexure and free-air gravity anomaly can be explained by volcano loading of a plate with an effective elastic thickness of 30 km and a load and infill density of 2600 kg m^− 3. The origin of the Cape Verde swell is poorly understood. An elastic thickness of 30 km is expected for the ~ 125 Ma old oceanic lithosphere beneath the islands, suggesting that the observed height of the swell and the elevated heat flow cannot be attributed to thermal reheating of the lithosphere. The lack of evidence for high densities and velocities in the lower crust and low densities and velocities in the upper mantle, suggests that neither a crustal underplate or a depleted swell root are the cause of the shallower than expected bathymetry and that, instead, the swell is supported by dynamic uplift associated with the underlying plume.     

Effusive history of the Grande Découverte Volcanic Complex, southern Basse-Terre (Guadeloupe, French West Indies) from new K–Ar Cassignol–Gillot ages

The Grande Découverte Volcanic Complex (GDVC), active since at least 0.2 Ma, is the most recent volcanic complex of the Basse-Terre Island (Guadeloupe, Lesser Antilles Arc). A detailed geochronological study using the K–Ar Cassignol–Gillot technique has been undertaken in order to reconstruct the history of effusive activity of this long-lived volcanic system. Twenty new ages permit to suggest that the GDVC experienced at least six main effusive stages, from 200 ka to present time. To the north of the GDVC, the GDS (Grande Découverte–Soufrière volcano) has been active since at least 200 ka, and to the south, the TRMF (Trois-Rivières–Madeleine Field), started to be emplaced 100 ka. Morphological investigations suggest that the whole TRMF volcanism was emitted from vents distinct from the GDS, most probably a large E–W fissure network linked to the Marie-Galante rift. The mean age of 62 ± 5 ka, obtained for the E–W Madeleine–Le Palmiste alignment suggests that a fissure-opening event occurred at that time. However, whole-rock major and trace element signatures are similar for both systems, suggesting that a common complex magma-plumbing system has fed the overall GDVC. We report very young ages for lava flows from the TRMF, which implies that < 10 ka volcanic activity is now identified for both massifs. Although hazards associated with such effusive volcanism are much lower than those associated with potential flank-collapse of the Soufrière lava dome or a magmatic dome eruption with explosive phases within the GDS, the emplacement of relatively large Holocene age lava flows (3–1 × 10⁸ m³) suggests that a revised integrated volcanic hazard assessment for Southern Basse-Terre should now consider the potential for renewed future activity from two Holocene volcanic centers including the TRMF.     

Geoelectric structure of the Ein Gedi sinkhole occurrence site at the Dead Sea shore in Israel

Hazardous sinkholes started to appear in alluvial fans and unconsolidated sediments along the western Dead Sea coast in 1990. Since then hundreds of sinkholes have appeared from north to south along the shoreline. The Electrical Resistivity Tomography (ERT) method was used to achieve a better understanding of the subsurface geoelectric structure at the sinkhole development sites, taking into account that electric parameters (such as resistivity or conductivity) are very sensitive to formation properties and their variations in time. Fifteen image lines were surveyed at the Ein Gedi area during a period of active sinkhole development (in 2001–2002) over an area of 300 × 550 m². Resistivity cross-sections and maps were constructed from 2-D linear surveys. The process of sinkhole formation in the surveyed area is located in a strip 50–70 m wide and 300–500 m long, extending approximately in a north–south direction. The sinkholes are arranged along a tortuous line within this strip. On resistivity maps and sections this U-shaped zone appears as an alternation of high resistivity anomalies of 350–1000 Ωm (at sinkhole group locations) with narrow background resistivity zones of 50–100 Ωm. The large size of resistivity anomalies (250 × 300 m²), which are considerably greater than those of the sinkholes, form one of the features of the sinkhole sites in the Ein Gedi area. The anomalies continue down to the water table or even deeper (maximum of 25–35 m depth). A low resistivity layer of 1–8 Ωm underlies them. The combined analysis of the image results and other geophysical data shows that high resistivity anomalies are associated with the decompaction of the soil mass at the sinkhole development sites and surrounding areas. Recent studies have shown that sinkholes in the Ein Gedi area are developing along the salt western edge located at a depth of 50 m. The subsurface high resistivity anomaly conforms to the sinkhole line (and salt boundary). They are presumably located above the great dissolution caverns at the salt edge. The heterogeneity of the resistivity structure within the high resistivity anomaly (seen in both lateral and vertical planes) confirms that a disintegration of internal formation structure takes place. Away from the sinkhole sites the subsurface resistivity distribution is homogeneous.     

Phase relations of iron and iron–nickel alloys up to 300 GPa: Implications for composition and structure of the Earth's inner core

We have investigated the phase relations of iron and iron–nickel alloys with 18 to 50 wt.% Ni up to over 300 GPa using a laser-heated diamond-anvil cell. The synchrotron X-ray diffraction measurements show the wide stability of hcp-iron up to 301 GPa and 2000 K and 319 GPa and 300 K without phase transition to dhcp, orthorhombic, or bcc phases. On the other hand, the incorporation of nickel has a remarkable effect on expanding the stability field of fcc phase. The geometry of the temperature–composition phase diagram of iron–nickel alloys suggests that the hcp–fcc–liquid triple point is located at 10 to 20 wt.% Ni at the pressure of the inner core boundary. The fcc phase could crystallize depending on the nickel and silicon contents in the Earth's core, both of which are fcc stabilizer.     

Sedimentological constraints on hydrometeor-enhanced particle deposition: 1992 Eruptions of Crater Peak, Alaska

Water is a dominant component of volcanic clouds and has fundamental control on very fine particle deposition. Particle size characteristics of distal tephra-fall (100s km from source volcano) have a higher proportion of very fine particles compared to predictions based on single particle settling rates. In this study, sedimentological analyses of fallout from for the 18 August and 16–17 September 1992 eruptions of Crater Peak, Alaska, are combined with satellite observations, and cloud trajectory and microphysics modeling to investigate meteorological influences on particle sedimentation. Total grain size distributions of tephra fallout were reconstructed for both Crater Peak eruptions and indicate a predominance of fine particles < 125 μm. Polymodal analysis of the deposits has identified a particle subpopulation with mode ~ 15–18 μm involved in particle aggregation. Accounting for the magmatic water source only, calculated ice water content of the 3.7 hour old September 1992 Spurr cloud was ~ 4.5 × 10^− 2 g m^− 3 (based on an estimated cloud thickness of ~ 1000 m from trajectory modeling). Hydrometeor formation on particles in the volcanic cloud and subsequent sublimation may induce a cloud base instability that leads to rapid bulk (en masse) sedimentation of very fine particles through a mammatus-like mechanism.     

Magnetic survey of topsoils in Windsor–Essex County, Canada

Windsor–Essex County is a major cross-border truck and transportation route, with significant localized industrialization as well as rural and farming areas. Magnetic property measurements (in-field and laboratory susceptibility, frequency-dependent susceptibility, hysteresis properties, thermomagnetic and thermosusceptibility curves, anhysteretic and isothermal magnetizations) were made in order to determine the potential for using such variables to distinguish between natural and anthropogenic pollutants. In-field magnetic susceptibility measured on 324 soil sampling sites on a 0.5–2 km grid spacing through Windsor–Essex County ranged from 3.7 × 10^− 6 to 305.2 × 10^− 6 SI (average 36.2 ± 35.8 × 10^− 6 SI), and showed that high magnetic susceptibility values were obtained on soil sampling sites in and around the cities/towns of Windsor, Harrow, Olinda and Oakland and near the beaches of Point Pelee National Park (PPNP) and Deerbrook, whereas lower susceptibility values were observed in near the towns of Lakeshore and Essex. On this grid spacing, Highway 401 (the major truck route) did not show anomalous susceptibility values; however, closer (1–3 m) sampling on other roads did show anomalously high values, suggesting that the coarser grid spacing may have missed anomalies. Laboratory measurements indicated that the dominant magnetic mineral in the Windsor–Essex County soils is magnetite; however, the grain size is variable. Pseudo-single domain (PSD)–multidomain (MD) magnetite is generally found on beaches and in PPNP, whereas single domain (SD)–PSD magnetite has been found near the City of Windsor and other towns. While certain correlations exist between some anthropogenic activities and the measured magnetic susceptibility and magnetic property values, no overall correlation can be made. A variety of geologic and anthropogenic factors must be considered when interpreting the origin of the magnetic signal in a particular area.     

Geochemistry and geochronology of Late Triassic volcanic rocks in the Chiang Khong region, northern Thailand

The Chiang Khong segment of the Chiang Khong–Lampang–Tak Volcanic Belt is composed of three broadly meridional sub‐belts of mafic to felsic volcanic, volcaniclastic, and associated intrusive rocks. Associated sedimentary rocks are largely non‐marine red beds and conglomerates. Three representative Chiang Khong lavas have Late Triassic (223–220 Ma) laser ablation inductively coupled mass‐spectroscopy U–Pb zircon ages. Felsic‐dominated sequences in the Chiang Khong Western and Central Sub‐belts are high‐K calc–alkaline rocks that range from basaltic to dominant felsic lavas with rare mafic dykes. The Western Sub‐belt lavas have slightly lower high field strength element contents at all fractionation levels than equivalent rocks from the Central Sub‐belt. In contrast, the Eastern Sub‐belt is dominated by mafic lavas and dykes with compositions transitional between E‐mid‐oceanic ridge basalt and back‐arc basin basalts. The Eastern Sub‐belt rocks have higher FeO* and TiO₂ and less light rare earth element enrichment than basalts in the high‐K sequences. Basaltic and doleritic dykes in the Western and Central sub‐belts match the composition of the Eastern Sub‐belt lavas and dykes. A recent geochemical study of the Chiang Khong rocks concluded that they were erupted in a continental margin volcanic arc setting. However, based on the dominance of felsic lavas and the mainly non‐marine associated sediments, we propose an alternative origin, in a post‐collisional extensional setting. A major late Middle to early Late Triassic collisional orogenic event is well documented in northern Thailand and Yunnan. We believe that the paucity of radiometric dates for arc‐like lavas in the Chiang Khong–Lampang–Tak Volcanic Belt that precede this orogenic event, coupled with the geochemistry of the Chiang Khong rocks, and strong compositional analogies with other post‐collisional magmatic suites, are features that are more typical of volcanic belts formed in a rapidly evolving post‐collisional, basin‐and range‐type extensional setting.     

Pre-eruptive and syn-eruptive conditions in the Black Butte, California dacite: Insight into crystallization kinetics in a silicic magma system

A series of experiments and petrographic analyses have been run to determine the pre-eruption phase equilibria and ascent dynamics of dacitic lavas composing Black Butte, a dome complex on the flank of Mount Shasta, California. Major and trace element analyses indicate that the Black Butte magma shared a common parent with contemporaneously erupted magmas at the Shasta summit. The Black Butte lava phenocryst phase assemblage (20 v.%) consists of amphibole, plagioclase (core An_77.5), and Fe–Ti oxides in a fine-grained (< 0.5 mm) groundmass of plagioclase, pyroxene, Fe–Ti oxides, amphibole, and cristobalite. The phenocryst assemblage and crystal compositions are reproduced experimentally between 890 °C and 910 °C, ≥ 300 MPa, X_H2O = 1, and oxygen fugacity = NNO + 1. This study has quantified the extent of three crystallization processes occurring in the Black Butte dacite that can be used to discern ascent processes. Magma ascent rate was quantified using the widths of amphibole breakdown rims in natural samples, using an experimental calibration of rim development in a similar magma at relevant conditions. The majority of rims are 34 ± 10 μm thick, suggesting a time-integrated magma ascent rate of 0.004–0.006 m/s among all four dome lobes. This is comparable to values for effusive samples from the 1980 Mount St. Helens eruption and slightly faster than those estimated at Montserrat. A gap between the compositions of plagioclase phenocryst cores and microlites suggests that while phenocryst growth was continuous throughout ascent, microlite formation did not occur until significantly into ascent. The duration of crystallization is estimated using the magma reservoir depth and ascent rate, as determined from phase equilibria and amphibole rim widths, respectively. Textural analysis of the natural plagioclase crystals yields maximum growth rates of plagioclase phenocryst rims and groundmass microlites of 8.7 × 10^− 8 and 2.5 × 10^− 8 mm/s, respectively. These rates are comparable to values determined from time-sequenced samples of dacite erupted effusively from Mount St. Helens during 1980 and indicate that syneruptive crystallization processes were important during the Black Butte eruptive cycle.     

Seismic crustal study of the Oslo Rift

A seismic refraction investigation across the southern part of the Oslo Rift has been made, based on quarry blasts at three localities. The study shows a three-layered crust with the followingP-wave velocities: . the upper mantleP-wave celocity, is 8.07 km/s. The velocity-depth relationship for the uppermost crust, obtained by solving the Wiechert-Herglotz integral equation numerically, shows a continuously decreasing velocity gradient in the region of the Oslo Rift which approaches zero at a depth of 9 km, the corresponding increase in theP-wave velocity being from 5.55 km/s to 6.34 km/s. The interface separating the subsurface layer ( =6.60 km/s) from the uppermost layer , interpreted as the Conrad discontinuity, is essentially horizontal in the investigated part of the Oslo Rift at a depth of approximately 15 km. A deep crustal layer with aP-wave velocity of 7.10 km/s appears to be related to the rift, though the top of this layer extends somewhat eastwards beneath the Precambrian rocks from the southern part of the rift at a depth of approximately 20 km. The Moho discontinuity is elevated beneath the Oslo Region compared with the surrounding area. A broad regional gravity high of about 45 mgal is observed along the entire rift zone. It is suggested that this anomaly is caused by the elevation of the sub-Conrad and Moho discontinuities during the rifting processes.     

Subsurface structure across the axis of the Tongariro Volcanic Centre,New Zealand

The relationship between structure and volcanism in the Tongariro Volcanic Centre, New Zealand, is largely masked by a mantle of young volcanic deposits. Here we report the results of an integrated geophysical investigation (using gravity, multi-level aeromagnetic and magnetotelluric methods) of subsurface deposits and structures in the upper 1–2 km across the axis of the Tongariro Volcanic Centre. Modelling of these data across the Tama Lakes saddle shows that the outcropping volcanic deposits are up to 800 m thick, underlain by Tertiary sediments (of a few 10's to a few 100 m in thickness) and in turn lying above a basement of probable Mesozoic greywacke. Basement faulting is shown to be concentrated in the centre of the rift, which is 18 km wide at this location, but no vertical offset is resolved at the rift axis. Vertical displacements on basement faults of 250–300 m are modelled giving a minimum total basement subsidence of 650 m. A 5 km-wide, deep low resistivity zone occurs at the axis of the rift which is interpreted as either resulting from extensive fracturing and/or hydrothermal alteration within the basement. Steep-sided volcanic bodies with a high proportion of lavas/dykes coincide with the Waihi fault and the rift axis. Coincidence with the Waihi Fault suggests that this fault system may have provided magma pathways to the surface and a focus for dyke emplacement, which could have contributed to rift extension. The lack of offset at the rift axis may reflect the juvenile nature of faulting at this location, which is consistent with the notion of a migration of faulting towards the centre of the graben, alternatively, rifting may have been entirely accommodated by dyke emplacement.     

The geoid in southeastern Brazil and adjacent regions: new constraints on density distribution and thermal state of the lithosphere

Least-squares collocation technique was used to process regional gravity data of the SE South American lithospheric plate in order to map intermediate (10–2000 km) wavelength geoid anomalies. The area between 35–10° S and 60–25° W includes the Paraná CFB Province, the Southern São Francisco Craton and its marginal fold/thrust belts, the Brazilian continental margin and oceanic basins. The main features in the geoid anomaly map are: (a) Paraná CFB Province is characterized by a 1000 km long and 500 km wide, NE-trending, 9 m-amplitude negative anomaly which correlates with the distribution of sediments and basalts within the Paraná basin. (b) A circular (600–800 km in diameter) positive, 8 m-amplitude geoid anomaly is located in the southern S. Francisco craton and extends into the northeastern border of the Paraná CFB Province. This anomaly partially correlates with Alto Paranaíba Igneous Province (APIP), where alkalic volcanism and tholeiitic dikes of ages younger than 80 Ma are found and where a low-velocity zone in the mantle has been mapped using seismic tomography. This positive geoid anomaly extends towards the continental margin at latitude 21° S and joins a linear sequence of short wavelength positive geoid anomalies associated with Vitoria–Trindade seamounts. (c) A NE-trending, 1000 km long and 800 km wide, 4 m-amplitude, positive geoid anomaly, which is located along the southeastern coast of Brazil, from latitude 24 to 35° S. The northern part of this anomaly correlates with the Ponta Grossa Arch and Florianopolis dyke swarm provinces. The age of this intrusive volcanism is 130–120 Ma. (d) A circular positive anomaly with 9 m of amplitude, located over the Rio Grande and Uruguay shields and offshore Pelotas basin. Few alkaline intrusives with ages between 65 and 80 Ma are found in the region and apatite fission track ages in basement rocks indicates cooling at around 30 Ma. A semi-quantitative analysis of the observed geoid anomalies using isostatic considerations suggests that the mechanism which generated Paraná CFB Province did not change, in a significant manner, the lithospheric thermal structure, since the same geoid pattern observed within this province continues northward over the Neoproterozoic fold/thrust belts systems separating the São Francisco and Amazon cratons. Therefore, this observation favours Anderson’s idea of rapid basaltic outpouring through a pull-apart mechanism along a major suture zone. A thermal component may still be present in the Southern São Francisco Craton and in the Rio Grande Shield and contiguous continental margins, sites of Tertiary thermal and magmatic reactivations.     

An Archean oceanic felsic dyke swarm in a nascent arc: the Hunter Mine Group, Abitibi Greenstone Belt, Canada

The 2730-Ma-old Hunter Mine Group (HMG), a dominantly felsic subaqueous volcanic sequence, was formed during early arc construction in the Abitibi greenstone belt (Quebec, Canada). The western part of the HMG contains a felsic dyke swarm up to 1.5 km wide and traceable up-section for 2.5 km. Five distinct generations were identified: (1) aphanitic to feldspar-phyric dykes; (2) quartz-feldspar-phyric dykes with < 5% quartz phenocrysts; (3) quartz-feldspar-phyric dykes with 10–25% quartz phenocrysts; (4) dacitic feldspar-phyric dykes; and (5) mafic dykes. The felsic dykes collectively constitute more than 90% of the dyke swarm. Geochemically, they resemble modern calc-alkaline dacites and rhyolites. Their mantle-normalized incompatible trace-element patterns display a moderate enrichment of Th and light REE relative to HFSE and heavy REE as well as negative Nb, Ta, Eu and Ti anomalies. Most of the major- and trace-element abundance variations in these rocks can be explained by crystallization of feldspars. Geochemical data including depleted mantle-like _Nd values suggest that an older sialic substrate was not involved in their genesis. We infer that the felsic rocks were generated by melting of mafic oceanic crust. The swarm was emplaced during nascent oceanic island-arc development and was related to rifting of the arc. The conformably overlying MORB-like basalts and basaltic komatiites of the Stoughton-Roquemaure Group used the same conduits and further indicate splitting of the arc. HMG and associated parts of the Abitibi greenstone belts bear a strong resemblance to modern rifted intraoceanic arcs of the western Pacific.     

Caldera-scale inflation of the Lazufre volcanic area, South America: Evidence from InSAR

Collapsed calderas are the structural surface expression of the largest volcanic eruptions on Earth and may reach diameters of tens of kilometres while erupting volumes larger than 1000 km³. Remnants of collapse calderas can be found along the South American volcanic arc and are thought to be inactive. However, this study shows that systems of such dimension may become active in a relatively short period of time without attracting much attention. Using satellite-based InSAR data, a 45 km wide elongated area of ground deformation was observed in the Lazufre volcanic region (Chile), where no deformation was detected 10 years ago. The deformation signal shows an uplift of up to ~ 3 cm yr^− 1 during 2003–2006, affecting an area of about 1100 km², comparable in size to super-volcanoes such as Yellowstone or Long Valley. This deformation signal can be explained by an inflating magma body at about 10 km depth, expanding and propagating laterally at a velocity of up to 4 km per year. Although it is not clear whether this intrusion will lead to an eruption, its dimensions and the rapid deformation rate insinuate that a potentially large volcanic system is forming.