Stress field control of eruption dynamics at a rift volcano: Nyamuragira, D.R.Congo

Using the record of 30 flank eruptions over the last 110 years at Nyamuragira, we have tested the relationship between the eruption dynamics and the local stress field. There are two groups of eruptions based on their duration (< 80days >) that are also clustered in space and time. We find that the eruptions fed by dykes parallel to the East African Rift Valley have longer durations (and larger volumes) than those eruptions fed by dykes with other orientations. This is compatible with a model for compressible magma transported through an elastic-walled dyke in a differential stress field from an over-pressured reservoir (Woods et al., 2006). The observed pattern of eruptive fissures is consistent with a local stress field modified by a northwest-trending, right lateral slip fault that is part of the northern transfer zone of the Kivu Basin rift segment. We have also re-tested with new data the stochastic eruption models for Nyamuragira of Burt et al. (1994). The time-predictable, pressure-threshold model remains the best fit and is consistent with the typically observed declining rate of sulphur dioxide emission during the first few days of eruption with lava emission from a depressurising, closed, crustal reservoir. The 2.4-fold increase in long-term eruption rate that occurred after 1977 is confirmed in the new analysis. Since that change, the record has been dominated by short-duration eruptions fed by dykes perpendicular to the Rift. We suggest that the intrusion of a major dyke during the 1977 volcano-tectonic event at neighbouring Nyiragongo volcano inhibited subsequent dyke formation on the southern flanks of Nyamuragira and this may also have resulted in more dykes reaching the surface elsewhere. Thus that sudden change in output was a result of a changed stress field that forced more of the deep magma supply to the surface. Another volcano-tectonic event in 2002 may also have changed the magma output rate at Nyamuragira.

Cyclic extrusion of a lava dome based on a stick-slip mechanism

Lava dome eruptions are sometimes characterised by large periodic fluctuations in extrusion rate over periods of hours that may be accompanied by Vulcanian explosions and pyroclastic flows. We consider a simple system of nonlinear equations describing a 1D flow of lava extrusion through a deep elastic dyke feeding a shallower cylindrical conduit in order to simulate this short-period cyclicity. Stick-slip conditions depending on a critical shear stress are assumed at the wall boundary of the cylindrical conduit. By analogy with the behaviour of industrial polymers in a plastic extruder, the elastic dyke acts like a barrel and the shallower cylindrical portion of the conduit as a die for the flow of magma acting as a polymer. When we applied the model to the Soufrière Hills Volcano, Montserrat, for which the key parameters have been evaluated from previous studies, cyclic extrusions with periods from 3 to 30 h were readily simulated, matching observations. The model also reproduces the reduced period of cycles observed when a major unloading event occurs due to lava dome collapse.

Pulsatory andesite lava flow at Bagana Volcano

Using a time series of TerraSAR-X spaceborne radar images we have measured the pulsatory motion of an andesite lava flow over a 14-month period at Bagana volcano, Papua New Guinea. Between October 2010 and December 2011, lava flowed continuously down the western flank of the volcano forming a 3 km-long blocky lava flow with a channel, levees, overflows and branches. We captured four successive pulses of lava advancing down the channel system, the first such behaviour of an andesite flow to be recorded using radar. Each pulse had a volume of the order of 107 m3 emplaced over many weeks. The average extrusion rate estimated from the radar data was 0.92 ± 0.35 m3 s-1 , and varied between 0.3 and 1.8 m3 s-1, with higher rates occurring earlier in each pulse. This, together with observations of sulphur dioxide emissions, explosions and incandescence suggest a variable supply rate of magma through Bagana’s conduit as the most likely source of the pulsatory behaviour.

Co-eruptive and inter-eruptive surface deformation measured by satellite radar interferometry at Nyamuragira volcano, D.R. Congo, 1996 to 2010.

Interferometric Synthetic Aperture Radar (InSAR) measurements of surface deformation at Nyamuragira Volcano between 1996 and 2010 reveal a variety of co-eruptive and inter-eruptive signals. During 7 of the 8 eruptions in this period deformation was measured that is consistent with the emplacement of shallow near-vertical dykes feeding the eruptive fissures and associated with a NNW-trending fissure zone that traverses the summit caldera. Between eruptions the caldera and the summit part of this fissure zone subsided gradually (b3–5 cm/year). We also find evidence of post-eruption subsidence around the sites of the main vents of some flank eruptions (2002, 2004, 2006, and 2010). In the 6 months prior to the 2010 eruption a10-km wide zone centred on the caldera inflated by 1–2 cm. The low magnitude of this signal suggests that the presumed magma reservoir at 3–8 km depth contains highly compressible magma with little stored elastic strain energy. To the north of the caldera the fissure zone splits into WNW and NE branches around a zone that has a distinct InSAR signal. We interpret this zone to represent an elevated, 'stable' block of basement rocks buried by lavas within the Rift Zone.

Coupled subdaily and multiweek cycles during the lava dome erupt of Soufriere Hills Volcano, Montserrat

Observations of volcanoes extruding andesitic lava to produce lava domes often reveal cyclic behaviour. At Soufriere Hills Volcano, Montserrat, cycles with sub-daily and multi-week periods have been recognised on many occasions. These two types of cycle have been modelled separately as stick-slip magma flow at the junction between a dyke and an overlying cylindrical conduit (Costa et al. 2012), and as the filling and discharge of magma through the elastic-walled dyke (Costa et al., 2007a) respectively. Here, we couple these two models to simulate the behaviour over a period of well-observed multi-week cycles, with accompanying sub-daily cycles, from 13 May to 21 September 1997. The coupled model captures well the asymmetrical first-order behaviour: the first 40% of the multi-week cycle consists of high rates of lava extrusion during short period/high amplitude sub-daily cycles as the dyke reservoir discharges itself. The remainder of the cycle involves increasing pressurization as more magma is stored, and extrusion rate falls, followed by a gradual increase in the period of the sub-daily cycles.

Emplacement mechanism of the main granite pluton of the Lavras do Sul intrusive complex, South Brazil, determined by magnetic anisotropies

Magnetic fabric and rock magnetism studies were performed on apparently isotropic granite facies from the main intrusion of the Lavras do Sul Intrusive Complex pluton (LSIC, Rio Grande do Sul, South Brazil). This intrusion is roughly circular (similar to 12 x 13.5 km), composed of alkali-calcic and alkaline granitoids, with the latter occupying the margin of the pluton. Magnetic fabrics were determined by applying both anisotropy of low-field magnetic susceptibility (AMS) and anisotropy of anhysteretic remanent magnetization (AARM). The two fabrics are coaxial. The parallelism between AMS and AARM tensors excludes the presence of a single domain (SD) effect on the AMS fabric of the granites. Several rock-magnetism experiments performed in one specimen from each sampled site show that for all sites the magnetic susceptibility is dominantly carried by ferromagnetic minerals, while mainly magnetite carries the magnetic fabrics. Lineations and foliations in the granite facies were successful determined by applying magnetic methods. Magnetic lineations are gently plunging and roughly parallel to the boundaries of the pluton facies, except at the few sites in the central facies which have a radial orientation pattern. In contrast, the magnetic foliations tend to follow the contacts between the different granite facies. They are gently outerward-dipping inside the pluton, and become either steeply southwesterly dipping or vertical towards its margin. The lack of solid-state and subsolidus deformations at outcrop scale and in thin sections precludes deformation after full crystallization of the pluton. This evidence allows us to interpret the observed magnetic fabrics as primary in origin (magmatic) acquired when the rocks were solidified as a result of processes reflecting magma flow. The foliation pattern displays a dome-shaped form for the main LSIC-pluton. However, the alkaline granites which outcrop in the southern part of the studied area have an inward-dipping foliation, and the steeply plunging magnetic lineation suggests that this area could be part of a feeder zone. The magma ascent probably occurred due to ring-diking. (C) 2008 Elsevier B.V. All rights reserved.

Tectonic fabric revealed by AARM of the proterozoic mafic dike swarm in the Salvador city (Bahia state): Sao Francisco Craton, NE Brazil

Magnetic fabric and rock magnetism studies were performed on 25 unmetamorphosed mafic dikes of the Meso-Late Proterozoic (similar to 1.02 Ga) dike swarm from Salvador (Bahia State, NE Brazil). This area lies in the north-eastern part of the Sao Francisco Craton, which was dominantly formed/reworked during the Transamazonian orogeny (2.14-1.94 Ga). The dikes crop out along the beaches and in quarries around Salvador city, and cut across both amphibolite dikes and granulites. Their widths range from a few centimeters up to 30 m with an average of similar to 4 m, and show two main trends N 140-190 and N 100-120 with vertical dips. Magnetic fabrics were determined using both anisotropy of low-field magnetic susceptibility (AMS) and anisotropy of anhysteretic remanent magnetization (AARM). The magnetic mineralogy was investigated by many experiments including remanent magnetization measurements at variable low temperatures (10-300 K), Mossbauer spectroscopy, high temperature magnetization curves (25-700 degrees C) and scanning electron microscopy (SEM). The rock magnetism study suggests pseudo-single-domain magnetite grains carrying the bulk magnetic susceptibility and AARM fabrics. The magnetite grains found in these dikes are large and we discard the presence of single-domain grains. Its composition is close to stoichiometric with low Ti substitution, and its Verwey transition occurs around 120 K. The main AMS fabric recognized in the swarm is so-called normal, in which the K(max)-K(int) plane is parallel to the dike plane and the magnetic foliation pole K(min)) is perpendicular to it. This fabric is interpreted as due to magma flow, and analysis of the K m inclination permitted to infer that approximately 80% of the dikes were fed by horizontal or sub-horizontal flows (K(max) < 30 degrees). This interpretation is supported by structural field evidence found in five dikes. In addition, based on the plunge of K(max), two mantle sources could be inferred; one of them which fed about 80% of the swarm would be located in the southern part of the region, and the other underlied the Valeria quarry. However, for all dikes the AARM tensors are not coaxial with AMS fabrics and show a magnetic lineation (AARM(max)) oriented to N30-60E, suggesting that magnetite grains were rotated clockwise from dike plane. The orientation of AARM lineation is similar to the orientation of a system of faults in which the Salvador normal fault is the most important. These faults were formed during Cretaceous rifting in the Reconcavo-Tucano-jatoba assemblage that corresponds to an aborted intra-continental rift formed during the opening of the South Atlantic. Therefore, the AARM fabric found for the Salvador dikes is probably tectonic in origin and suggests that the dike swarm was affected by the important tectonic event responsible for the break-up of the Gondwanaland. (C) 2008 Elsevier B.V. All rights reserved.

Precambrian history of the Zona Transversal of the Borborema Province, NE Brazil: Insights from Sm-Nd and U-Pb geochronology

The Borborema Province has three major subprovinces. The northern subprovince lies north of the Patos shear zone and is comprised of Paleoproterozoic cratonic basement with Archean nuclei, plus overlying Neoproterozoic supracrustal rocks and Brasiliano plutonic rocks. The central subprovince occurs between the Patos and Pernambuco shear zones and is mainly comprised of the Zona Transversal. The southern subprovince occurs between the Pernamabuco shear zone and the Sao Francisco craton and is comprised of a tectonic collage of various blocks, terranes, or domains ranging in age from Archean to Neoproterozoic. This report focuses on the Zona Transversal, especially on Brasiliano rocks for which we have the most new information. Paleoproterozoic gneisses with ages of 2.0-2.2 Ga occur discontinuously throughout the Zona Transversal. The Cariris Velhos suite consists of metavolcanic, metasedimentary, and metaplutonic rocks yielding U-Pb zircon ages of 995-960 Ma. This suite is mainly confined to a 100 km wide belt that extends for more than 700 km within the Alto Pajeu terrane. Sm-Nd model ages in metaigneous rocks cluster about 1.3-1.6 Ga, indicating that older crust was involved in genesis of their magmas. Brasiliano supracrustal rocks dominate the Pianco-Alto Brigida terrane, and they probably also constitute significant parts of the Alto Pajeu and Rio Capibaribe terranes. They are only slightly older than early stages of Brasiliano plutonism, with detrital zircon ages at least as young as 620 Ma; most T(DM) ages range from 1.2 to 1.6 Ga. Brasiliano plutons range from ca. 640 to 540 Ma, and their T(DM) ages range from 1.2 to 2.5 Ga. Previous workers have shown significant correlations among U-Pb ages, Sm-Nd model ages, petrology, and geochemistry, and we are able to reinforce and extend these correlations. Stage I plutons formed 640 -610 Ma and have T(DM) ages less than 1.5 Ga. Stage 11 (610-590 Ma) contains few plutons, but coincides with the peak of compressional deformation, metamorphism, and formation of migmatites. Stage III plutons (590 to ca. 575 Ma) have older T(DM) ages (ca. 1.8-2.0 Ga), as do Stage IV plutons (575 to ca. 550 Ma; T(DM) from 1.9 to 2.4 Ga). Stage III plutons formed during the transition from compressional to transcurrent deformation, while Stage IV plutons are mainly post-tectonic. Stage V plutons (550-530 Ma) are commonly undeformed (except along younger shear zones) and have A-type geochemistry. The five stages have distinct geochemical properties, which suggest that the tectonic settings evolved from early, arc-related magma-genesis (Stage I) to within-plate magma-genesis (Stage V), with perhaps some intermediate phases of extensional environments. (C) 2011 Elsevier Ltd. All rights reserved.

Evaluation of petrogenetic models for intermediate and silicic plutonic rocks from the Sierra de Valle Fertil-La Huerta, Argentina: Petrologic constraints on the origin of igneous rocks in the Ordovician Famatinian-Puna paleoarc

The whole Valle Fertil-La Huerta section appears as a calc-alkaline plutonic suite typical of a destructive plate margin. New Sr and Nd isotopic whole-rock data and published whole-rock geochemistry suggest that the less-evolved intermediate (dioritic) rocks can be derived by magmatic differentiation, mainly by hornblende + plagioclase +/- Fe-Ti oxide fractional crystallization, from mafic (gabbroic) igneous precursors. Closed-system differentiation, however, cannot produce the typical intermediate (tonalitic) and silicic (granodioritic) plutonic rocks, which requires a preponderant contribution of crustal components. Intermediate and silicic plutonic rocks from Valle Fertil-La Huerta section have formed in a plate subduction setting where the thermal and material input of mantle-derived magmas promoted fusion of fertile metasedimentary rocks and favored mixing of gabbroic or dioritic magmas with crustal granitic melts. Magma mixing is observable in the field and evident in variations of chemical elemental parameters and isotopic ratios, revealing that hybridization coupled with fractionation of magmas took place in the crust. Consideration of the whole-rock geochemical and isotopic data in the context of the Famatinian-Puna magmatic belt as a whole demonstrates that the petrologic model postulated for the Sierra Valle Fertil-La Huerta section has the potential to explain the generation of plutonic and volcanic rocks across the Early Ordovician paleoarc from central and northwestern Argentina. As the petrologic model does not require the intervention of old Precambrian continental crust, the nature of the basement on which thick accretionary turbiditic sequences were deposited remains a puzzling aspect. Discussion in this paper provides insights into the nature of magmatic source rocks and mechanisms of magma generation in Cordilleran-type volcano-plutonic arcs of destructive plate margins. (C) 2010 Elsevier Ltd. All rights reserved.

Two- and three-dimensional shape fabric analysis by the intercept method in grey levels

The count intercept is a robust method for the numerical analysis of fabrics Launeau and Robin (1996). It counts the number of intersections between a set of parallel scan lines and a mineral phase, which must be identified on a digital image. However, the method is only sensitive to boundaries and therefore supposes the user has some knowledge about their significance. The aim of this paper is to show that a proper grey level detection of boundaries along scan lines is sufficient to calculate the two-dimensional anisotropy of grain or crystal distributions without any particular image processing. Populations of grains and crystals usually display elliptical anisotropies in rocks. When confirmed by the intercept analysis, a combination of a minimum of 3 mean length intercept roses, taken on 3 more or less perpendicular sections, allows the calculation of 3-dimensional ellipsoids and the determination of their standard deviation with direction and intensity in 3 dimensions as well. The feasibility of this quick method is attested by numerous examples on theoretical objects deformed by active and passive deformation, on BSE images of synthetic magma flow, on drawing or direct analysis of thin section pictures of sandstones and on digital images of granites directly taken and measured in the field. (C) 2010 Elsevier B.V. All rights reserved.

Deformation and geochronology of syntectonic granitoids emplaced in the Major Gercino Shear Zone, southeastern South America

The Major Gercino Shear Zone is one of the NE-SW lineaments that separate the Neoproterozoic Dom Feliciano Belt, of Brazil and Uruguay, into two different domains: a northwestern supracrustal domain from a southeastern granitoid domain. The shear zone, striking NE, is composed of protomylonites to ultramylonites with mainly dextral kinematic indicators. In Santa Catarina State, southern Brazil, the shear zone is composed of two mylonite belts. The mylonites have mineral orientations produced under greenschist fades conditions at a high strain rate. Strong flattening and coaxial deformation indicate the transpressive character, while the role of pure shear is emphasized by the orientation of the mylonite belts in relation to the inferred stress field component. The quartz microstructures point out that different dynamic recrystallization regimes and crystal plasticity were the dominant mechanisms of deformation during the mylonitization process. Additionally, the fabrics suggest that the glide systems are activated for deformation conditions compatible with the metamorphism in the middle greenschist facies. Elongated granitoid intrusions belonging to two petrographically, geochemically and isotopically distinct rock associations occur between the two mylonite belts. The structures observed in the granites result from a deformation range from magmatic to solid-state conditions points to a continuum of magma straining during and just after its crystallization. Conventional U-Pb analysis of multi-crystal zircon fractions yielded essentially identical ages of 609 +/- 16 Ma and 614 +/- 2 Ma for the two granitic associations, and constrain the transpressive phase of the shear zone. K-Ar ages of biotites between 585 and 560 Ma record the slow cooling and uplift of the intrusions. Some K-Ar ages of micas in regional mylonites are similar, suggesting that thermo-tectonic activity was intense up to this time, probably related to the agglutination of the granite belt to the supracrustal belt NW of the MGSZ. (C) 2009 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

Petrogenesis of Upper Paleozoic post-collisional peraluminous leucogranites, Sierra de Ancasti, northwest Argentina

The Santa Rosa and Sauce Guacho plutons are two post-collisional peraluminous Late Devonian to Early Carboniferous leucogranites that intruded the banded schists of the Ancasti Formation. The leucogranites are composed of microcline phenocrysts along with quartz, plagioclase, muscovite, biotite, ilmenite, tourmaline, apatite, monazite and zircon. Their geochemical composition is consistent with S-type granites and mineralogically they belong to MPG granites (muscovite-peraluminous granites). It is proposed that granite magma generation was related to shear zones that concentrated fluids in the metasedimentary crust during a collision or transcurrent tectonics. U-Pb analyses on monazite gave an age of 369.8 +/- 5.3 Ma, while Sm/Nd isotopic data yield epsilon(Nd(t)) values of -5.3 for Sauce Guacho and -5.7 for Santa Rosa indicating crustal provenance. Nd model ages between 1,544 and 1,571 Ma are within the range of magmatic rocks from the Lower Ordovician Famatinian Arc in the Central Sierras Pampeanas.

Thermochronology of central Ribeira Fold Belt, SE Brazil: Petrological and geochronological evidence for long-term high temperature maintenance during Western Gondwana amalgamation

The studied sector of the central Ribeira Fold Belt (SE Brazil) comprises metatexites, diatexites, charnockites and blastomylonites. This study integrates petrological and thermochronological data in order to constrain the thermotectonic and geodynamic evolution of this Neoproterozoic-Ordovician mobile belt during Western Gondwana amalgamation. New data indicate that after an earlier collision stage at similar to 610 Ma (zircon, U-Pb age), peak metamorphism and lower crust partial melting, coeval with the main regional high grade D(1) thrust deformation, occurred at 572-562 Ma (zircon, U-Pb ages). The overall average cooling rate was low (<5 degrees C/Ma) from 750 to 250 degrees C (at similar to 455 Ma; biotite-WR Rb-Sr age), but disparate cooling paths indicate differential uplift between distinct lithotypes: (a) metatexites and blastomylonites show a overall stable 3-5 degrees C/Ma cooling rate; (b) charnockites and associated rocks remained at T>650 degrees C during sub-horizontal D(2) shearing until similar to 510-470 Ma (garnet-WR Sm-Nd ages) (1-2 degrees C/Ma), being then rapidly exhumed/cooled (8-30 degrees C/Ma) during post-orogenic D(3) deformation with late granite emplacement at similar to 490 Ma (zircon, U-Pb age). Cooling rates based on garnet-biotite Fe-Mg diffusion are broadly consistent with the geochronological cooling rates: (a) metatexites were cooled faster at high temperatures (6 degrees C/Ma) and slowly at low temperatures (0.1 degrees C/Ma), decreasing cooling rates with time; (b) charnockites show low cooling rates (2 degrees C/Ma) near metamorphic peak conditions and high cooling rates (120 degrees C/Ma) at lower temperatures, increasing cooling rates during retrogression. The charnockite thermal evolution and the extensive production of granitoid melts in the area imply that high geothermal gradients were sustained fora long period of time (50-90 Ma). This thermal anomaly most likely reflects upwelling of asthenospheric mantle and magma underplating coupled with long-term generation of high HPE (heat producing elements) granitoids. These factors must have sustained elevated crustal geotherms for similar to 100 Ma, promoting widespread charnockite generation at middle to lower crustal levels. (C) 2010 Elsevier B.V. All rights reserved.

Petrogenesis of Early Cretaceous silicic volcanism in SE Uruguay: The role of mantle and crustal sources

Early Cretaceous (similar to 129 Ma) silicic rocks crop out in SE Uruguay between the Laguna Merin and Santa Lucia basins in the Lascano, Sierra Sao Miguel. Salamanca and Minas areas They are mostly rhyolites with minor quartz-trachytes and are nearly contemporaneous with the Parana-Etendeka igneous province and with the first stages of South Atlantic Ocean opening A strong geochemical variability (particularly evident from Rb/Nb, Nb/Y trace element ratios) and a wide range of Sr-Nd isotopic ratios ((143)Nd/(144)Nd((129)) = 0.51178-0.51209, (87)Sr/(86)Sr((129)) = 0.70840-0.72417) characterize these rocks Geochemistry allows to distiniguish two compositional groups, corresponding to the north-eastern (Lascano and Sierra Sao Miguel, emplaced on the Neo-Proterozoic southern sector of the Dom Feliciano mobile belt) and south-eastern localities (Salamanca, Minas, emplace on the much older (Archean) Nico Perez teriane or on the boundary between the Dom Feliciano and Nico Perez termites) These compositional differences between the two groups are explained by variable mantle source and crust contributions. The origin of the silicic magmas is best explained by complex processes involving assimilation and fractional crystallization and mixing of a basaltic magma with upper crustal lithologies, for Lascano and Sierra Sao Miguel rhyolites. In the Salamanea and Minas rocks genesis, a stronger contribution from lower crust is indicated.