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.

Segmental dynamics in entangled linear polymer melts

We present molecular dynamics (MD) and slip-springs model simulations of the chain segmental dynamics in entangled linear polymer melts. The time-dependent behavior of the segmental orientation autocorrelation functions and mean-square segmental displacements are analyzed for both flexible and semiflexible chains, with particular attention paid to the scaling relations among these dynamic quantities. Effective combination of the two simulation methods at different coarse-graining levels allows us to explore the chain dynamics for chain lengths ranging from Z ≈ 2 to 90 entanglements. For a given chain length of Z ≈ 15, the time scales accessed span for more than 10 decades, covering all of the interesting relaxation regimes. The obtained time dependence of the monomer mean square displacements, g1(t), is in good agreement with the tube theory predictions. Results on the first- and second-order segmental orientation autocorrelation functions, C1(t) and C2(t), demonstrate a clear power law relationship of C2(t) C1(t)m with m = 3, 2, and 1 in the initial, free Rouse, and entangled (constrained Rouse) regimes, respectively. The return-to-origin hypothesis, which leads to inverse proportionality between the segmental orientation autocorrelation functions and g1(t) in the entangled regime, is convincingly verified by the simulation result of C1(t) g1(t)−1 t–1/4 in the constrained Rouse regime, where for well-entangled chains both C1(t) and g1(t) are rather insensitive to the constraint release effects. However, the second-order correlation function, C2(t), shows much stronger sensitivity to the constraint release effects and experiences a protracted crossover from the free Rouse to entangled regime. This crossover region extends for at least one decade in time longer than that of C1(t). The predicted time scaling behavior of C2(t) t–1/4 is observed in slip-springs simulations only at chain length of 90 entanglements, whereas shorter chains show higher scaling exponents. The reported simulation work can be applied to understand the observations of the NMR experiments.

Dissipation of shear-free turbulence near boundaries

The rapid-distortion model of Hunt & Graham (1978) for the initial distortion of turbulence by a flat boundary is extended to account fully for viscous processes. Two types of boundary are considered: a solid wall and a free surface. The model is shown to be formally valid provided two conditions are satisfied. The first condition is that time is short compared with the decorrelation time of the energy-containing eddies, so that nonlinear processes can be neglected. The second condition is that the viscous layer near the boundary, where tangential motions adjust to the boundary condition, is thin compared with the scales of the smallest eddies. The viscous layer can then be treated using thin-boundary-layer methods. Given these conditions, the distorted turbulence near the boundary is related to the undistorted turbulence, and thence profiles of turbulence dissipation rate near the two types of boundary are calculated and shown to agree extremely well with profiles obtained by Perot & Moin (1993) by direct numerical simulation. The dissipation rates are higher near a solid wall than in the bulk of the flow because the no-slip boundary condition leads to large velocity gradients across the viscous layer. In contrast, the weaker constraint of no stress at a free surface leads to the dissipation rate close to a free surface actually being smaller than in the bulk of the flow. This explains why tangential velocity fluctuations parallel to a free surface are so large. In addition we show that it is the adjustment of the large energy-containing eddies across the viscous layer that controls the dissipation rate, which explains why rapid-distortion theory can give quantitatively accurate values for the dissipation rate. We also find that the dissipation rate obtained from the model evaluated at the time when the model is expected to fail actually yields useful estimates of the dissipation obtained from the direct numerical simulation at times when the nonlinear processes are significant. We conclude that the main role of nonlinear processes is to arrest growth by linear processes of the viscous layer after about one large-eddy turnover time.

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.

Rheology of discrete failure regimes of anisotropic sea ice

A rheological model of sea ice is presented that incorporates the orientational distribution of ice thickness in leads embedded in isotropic floe ice. Sea ice internal stress is determined by coulombic, ridging and tensile failure at orientations where corresponding failure criteria are satisfied at minimum stresses. Because sea ice traction increases in thinner leads and cohesion is finite, such failure line angles are determined by the orientational distribution of sea ice thickness relative to the imposed stresses. In contrast to the isotropic case, sea ice thickness anisotropy results in these failure lines becoming dependent on the stress magnitude. Although generally a given failure criteria type can be satisfied at many directions, only two at most are considered. The strain rate is determined by shearing along slip lines accompanied by dilatancy and closing or opening across orientations affected by ridging or tensile failure. The rheology is illustrated by a yield curve determined by combining coulombic and ridging failure for the case of two pairs of isotropically formed leads of different thicknesses rotated with regard to each other, which models two events of coulombic failure followed by dilatancy and refreezing. The yield curve consists of linear segments describing coulombic and ridging yield as failure switches from one lead to another as the stress grows. Because sliding along slip lines is accompanied by dilatancy, at typical Arctic sea ice deformation rates a one-day-long deformation event produces enough open water that these freshly formed slip lines are preferential places of ridging failure.

Modeling Coulombic failure of sea ice with leads

[1] Sea ice failure under low-confinement compression is modeled with a linear Coulombic criterion that can describe either fractural failure or frictional granular yield along slip lines. To study the effect of anisotropy we consider a simplified anisotropic sea ice model where the sea ice thickness depends on orientation. Accommodation of arbitrary deformation requires failure along at least two intersecting slip lines, which are determined by finding two maxima of the yield criterion. Due to the anisotropy these slip lines generally differ from the standard, Coulombic slip lines that are symmetrically positioned around the compression direction, and therefore different tractions along these slip lines give rise to a nonsymmetric stress tensor. We assume that the skewsymmetric part of this tensor is counterbalanced by an additional elastic stress in the sea ice field that suppresses floe spin. We consider the case of two leads initially formed in an isotropic ice cover under compression, and address the question of whether these leads will remain active or new slip lines will form under a rotation of the principal compression direction. Decoupled and coupled models of leads are considered and it is shown that for this particular case they both predict lead reactivation in almost the same way. The coupled model must, however, be used in determining the stress as the decoupled model does not resolve the stress asymmetry properly when failure occurs in one lead and at a new slip line.

Modelling the rheology of sea ice as a collection of diamond-shaped floes

In polar oceans, seawater freezes to form a layer of sea ice of several metres thickness that can cover up to 8% of the Earth’s surface. The modelled sea ice cover state is described by thickness and orientational distribution of interlocking, anisotropic diamond-shaped ice floes delineated by slip lines, as supported by observation. The purpose of this study is to develop a set of equations describing the mean-field sea ice stresses that result from interactions between the ice floes and the evolution of the ice floe orientation, which are simple enough to be incorporated into a climate model. The sea ice stress caused by a deformation of the ice cover is determined by employing an existing kinematic model of ice floe motion, which enables us to calculate the forces acting on the ice floes due to crushing into and sliding past each other, and then by averaging over all possible floe orientations. We describe the orientational floe distribution with a structure tensor and propose an evolution equation for this tensor that accounts for rigid body rotation of the floes, their apparent re-orientation due to new slip line formation, and change of shape of the floes due to freezing and melting. The form of the evolution equation proposed is motivated by laboratory observations of sea ice failure under controlled conditions. Finally, we present simulations of the evolution of sea ice stress and floe orientation for several imposed flow types. Although evidence to test the simulations against is lacking, the simulations seem physically reasonable.

Studying the contact point and interface moving in a sinusoidal tube with lattice Boltzmann method

The multicomponent nonideal gas lattice Boltzmann model by Shan and Chen (S-C) is used to study the immiscible displacement in a sinusoidal tube. The movement of interface and the contact point (contact line in three-dimension) is studied. Due to the roughness of the boundary, the contact point shows "stick-slip" mechanics. The "stick-slip" effect decreases as the speed of the interface increases. For fluids that are nonwetting, the interface is almost perpendicular to the boundaries at most time, although its shapes at different position of the tube are rather different. When the tube becomes narrow, the interface turns a complex curves rather than remains simple menisci. The velocity is found to vary considerably between the neighbor nodes close to the contact point, consistent with the experimental observation that the velocity is multi-values on the contact line. Finally, the effect of three boundary conditions is discussed. The average speed is found different for different boundary conditions. The simple bounce-back rule makes the contact point move fastest. Both the simple bounce-back and the no-slip bounce-back rules are more sensitive to the roughness of the boundary in comparison with the half-way bounce-back rule. The simulation results suggest that the S-C model may be a promising tool in simulating the displacement behaviour of two immiscible fluids in complex geometry.

Intraplate earthquake swarm in Belo Jardim, NE Brazil: reactivation of a major Neoproterozoic shear zone (Pernambuco Lineament)

Intraplate earthquakes in stable continental areas have been explained basically by reactivation of pre-existing zones of weakness, stress concentration, or both. Zones of weakness are usually identified as sites of the last major orogeny, provinces of recent alkaline intrusions, or stretched crust in ancient rifts. However, it is difficult to identify specific zones of weakness and intraplate fault zones are not always easily correlated with known geological features. Although Northeastern Brazil is one of the most seismically active areas in the country (magnitudes 5 roughly every 5 yr), with hypocentral depths shallower than similar to 10 km and seismic zones as long as 30-40 km, no clear relationship with the known surface geology can be usually established with confidence, and a clear identification of zones of weakness has not yet been possible. Here we present the first clear case of seismic activity occurring as reactivation of an old structure in Brazil: the Pernambuco Lineament, a major Neoproterozoic shear zone. The 2004 earthquake swarm of Belo Jardim (magnitudes up to 3.1) and the recurrent activities in the nearby towns of Sao Caetano and Caruaru (magnitudes up to 4.0 and 3.8), show that the Pernambuco Lineament is a weak zone. A local seismic network showed that the Belo Jardim swarm of 2004 November occurred by normal faulting on a North dipping, E-W oriented fault plane in close agreement with the E-W trending structures within the Pernambuco Lineament. The Belo Jardim activity was concentrated in a 1.5 km (E-W) by 2 km (downdip) fault area, and average depth of 4.5 km. The nearby Caruaru activity occurs as both strike-slip and normal faulting, also consistent with local structures of the Pernambuco Lineament. The focal mechanisms of Belo Jardim, Caruaru and S. Caetano, indicate E-W compressional and N-S extensional principal stresses. The NS extension of this stress field is larger than that predicted by numerical models such as those of Coblentz & Richardson and we propose that additional factors such as flexural stresses from the nearby Sergipe-Alagoas marginal basin could also affect the current stress field in the Pernambuco Lineament.

The intraplate Porto dos Gauchos seismic zone in the Amazon craton - Brazil

The largest earthquake observed in the stable continental interior of the South American plate occurred in Serra do Tombador, Mato Grosso state - Brazil, on January 31,1955 with a magnitude of 6.2 m(b). Since then no other earthquake has been located near the 1955 epicentre. However, in Porto dos Gauchos, 100 km northeast of Serra do Tombador, a recurrent seismicity has been observed since 1959. Both Serra do Tombador and Porto dos Gauchos are located in the Phanerozoic Parecis basin. Two magnitude 5 earthquakes occurred in Porto dos Gauchos, in 1998 and 2005, with intensities up to VI and V, respectively. These two main shocks were followed by aftershock sequences lasting more than three years each. Local seismic stations have been deployed by the Seismological Observatory of the University of Brasilia since 1998 to study the ""Porto dos Gauchos"" seismic zone (PGSZ). A local seismic refraction survey was carried out with two explosions to help define the seismic velocity model. Both the 1998 and 2005 earthquake sequences occurred in the same WSW-ENE oriented fault zone with right-lateral strike-slip mechanisms. The epicentral zone is in the Parecis basin, near its northern border where there are buried grabens, generally trending WNW-ESE, such as the deep Mesoproterozoic Caiabis graben which lies partly beneath the Parecis basin. However, the epicentral distribution indicates that the 1998 and 2005 sequences are related to a N60 degrees E fault which probably crosses the entire Caiabis graben. The 1955 earthquake, despite the uncertainty in its epicentre, does not seem to be directly related to any buried graben either. The seismicity in the Porto dos Gauchos seismic zone, therefore, is not directly related to rifted crust. The probable direction of the maximum horizontal stress near Porto dos Gauchos is roughly E-W, consistent with other focal mechanisms further south in the Pantanal basin and Paraguay. but seems to be different from the NW-SE direction observed further north in the Amazon basin. The recurrent seismicity observed in Porto dos Gauchos, and the large 1955 earthquake nearby, make this area of the Parecis basin one of the most important seismic zones of Brazil. (C) 2009 Elsevier B.V. All rights reserved.

Induced seismicity in the Castarthao reservoir, NE Brazil - Preliminary results

The Castanhao reservoir was built in the state of Ceara, a dry region in Northeastern Brazil, to regulate the flow of the Jaguaribe River, for irrigation, and for power generation. It is an earth-filled dam, 60 m high, with a water capacity of 4.5 x 10(9) m(3). The seismicity in the area has been monitored since 1998, with a few interruptions, using one analog or one digital station and, during a few periods, a three-station network. The first earthquakes likely to be induced events were detected in 2003, when the water level was about 20 in high. In early 2004 a very heavy rainfall season quickly filled the reservoir. Shortly after, an increase in the seismic activity occurred and many micro-earthquakes were recorded. We suggest that this activity resulted from an increase in pore pressure due to undrained response. Therefore, we may classify this cluster of microearthquakes as ""initial seismicity."" We deployed a network with four analog stations in the area, following this activity, to determine the epicentral zone. At least three epicentral areas under the reservoir were detected. The spatio-temporal analysis of the available data revealed that the seismicity occurs in clusters and that these were activated at different periods. We identified four sets of faults (N-S-, E-W-, NW-SE-, and NE-SW-oriented), some of which moved in shallow crustal levels and as recently as the Quaternary (1.8 Ma). Under the present-day stress regime, the last two sets moved as strike-slip structures. We suggest a possible correlation between dormant faults and the observed induced seismicity. (c) 2008 Elsevier B.V. All rights reserved.

Late Quaternary sedimentation in the Paraiba Basin, Northeastern Brazil: Landform, sea level and tectonics in Eastern South America passive margin

Late Quaternary deposits in the northeastern Brazil have been scarcely investigated, despite their relevance to the discussion of the post-rift evolution of the South American passive margin within the context of landform, sea level and tectonic deformation. Sedimentological, stratigraphic and morphological characterization of these deposits, referred as Post-Barreiras Sediments, led to their distinction from underlying Early/Middle Miocene strata. Based on optically stimulated luminescence (OSL) dating, two sedimentary units (PB1 and PB2) were recognized and related to the time intervals between 74.8 +/- 9.3 and 30.8 +/- 6.9 ka, and 8.8 +/- 0.9 and 1.8 +/- 0.2 ka, respectively. Unit PB1 consists of indurated sandstones and breccias either with massive bedding or complex types of soft sediment deformation structures generated by contemporaneous seismic activity. Unit PB2 is composed of massive sands or sands related to structures developed by dissipation of dunes. The present work, focusing on the Post-Barreiras Sediments, discusses landform, sea level and tectonics of the eastern South American passive margin during the latest Quaternary. Non-deposition and sub-aerial exposure related to the Tortonian worldwide low sea level combined with tectonic quiescence followed the Miocene transgression. Tectonic deformation in the latest Pleistocene created space to accommodate unit PB1 in downthrown faulted blocks and, perhaps, also synclines produced by strike-slip deformation. Although deposition of this unit was simultaneous with the progressive fall in sea level that followed the Last Interglacial Maximum, punctuated rises combined with land subsidence led to marine deposition close to the modern coastline. Renewed subsidence in the Holocene gave rise to accommodation of the Post-Barreiras Sediments. Most of unit PB2 was deposited during the Holocene Transgression, but it is not composed of marine sediments, which suggests either an insignificant rise in relative sea level or aeolian reworking of thin transgressive sands. The data presented here lead to a review of the evolution of the South American passive margin based on assumptions of uniform sedimentation and undeformed planation surfaces over a wide coastal area of the northeastern Brazil. (C) 2011 Elsevier B.V. All rights reserved.

Quartz recrystallization regimes, c-axis texture transitions and fluid inclusion reequilibration in a prograde greenschist to amphibolite facies mylonite zone (Ribeira Shear Zone, SE Brazil)

The crystal-plastic behavior of quartz mylonites from the Ribeira Shear Zone (SE Brazil), a major strike-slip structure that was active during a prograde metamorphic phase related to the Neoproterozoic Brasiliano-Pan African Orogeny, was investigated using a multi-method approach. Geothermobarometry results indicate deformational conditions ranging from similar to 300 to similar to 630 degrees C and 500-700 MPa. A strong correlation between mapped metamorphic zones and a dominance of different dynamic recrystallization mechanisms of quartz occurs within the mylonite zone. Bulging recrystallization (BLG) dominates within the chlorite zone between 300 and 410 degrees C, subgrain rotation recrystallization (SGR) operates within the biotite zone from 410 to 520 degrees C, and grain boundary migration recrystallization (GBM) dominates in the garnet zone above 520 degrees C. The development of quartz c-axis textures is mainly governed by temperature and dynamic recrystallization mechanisms. Textures from BLG zone mylonites are characterized by maxima around Z; SGR zone mylonites display single girdles or asymmetric type I crossed girdles; and GBM zone mylonites comprise maxima around Y and intermediate between X and Z. The scarcity or absence of water-bearing fluid inclusions in quartz mylonites from the SGR and GBM zones, which are dominated by carbonic inclusions, suggests water-deficient conditions, whereas BLG zone mylonites are dominated by water-bearing inclusions. This evidence indicates that water was available in the protoliths but has been eliminated with increasing deformation and deformation temperature. No effect of the water content variation on the quartz microstructural and recrystallized grain size evolution was detected, and little influence on c-axis texture development was observed. Most of the fluid inclusion densities were reequilibrated during the shear zone exhumation history, recording a decompression in the range of 300-500 MPa, while microstructural reequilibration effects related to the prograde metamorphism are largely preserved. Fluid inclusion microstructures and densities from two SGR zone samples preserved evidence for a near isothermal compression within the interior of the Ribeira Shear Zone during the prograde metamorphism. (C) 2009 Elsevier B.V. All rights reserved.

Quaternary tectonics in a passive margin: Marajo Island, northern Brazil

Marajo Island is located in a passive continental margin that evolved from rifting associated with the opening of the Equatorial South Atlantic Ocean in the Late Jurassic/Early Cretaceous period. This study, based on remote sensing integrated with sedimentology, as well as subsurface and seismographic data available from the literature, allows discussion of the significance of tectonics during the Quaternary history of marginal basins. Results show that eastern Marajo Island contains channels with evidence of tectonic control. Mapping of straight channels defined four main groups of lineaments (i.e. NNE-SSW, NE-SW, NW-SE and E-W) that parallel main normal and strike-slip fault zones recorded for the Amazon region. Additionally, sedimentological studies of late Quaternary and Holocene deposits indicate numerous ductile and brittle structures within stratigraphic horizons bounded by undeformed strata, related to seismogenic deformation during or shortly after sediment deposition. This conclusion is consistent with subsurface Bouguer mapping suggestive of eastern Marajo Island being still part of the Marajo graben system, where important fault reactivation is recorded up to the Quaternary. Together with the recognition of several phases of fault reactivation, these data suggest that faults developed in association with rift basins might remain active in passive margins, imposing important control on development of depositional systems. Copyright (C) 2007 John Wiley & Sons, Ltd.

Hybrid photocatalysts for the degradation of trichloroethylene in air

Hybrid photocatalysts based on an adsorbent SiMgOx and a photocatalyst TiO(2) were developed in a plate shape. The ceramic surface was coated with TiO(2) by the slip-casting technique. The effect of the support in the photocatalytic degradation of trichloroethylene (TCE) was analyzed by modifying TiO(2) loading and the layer thickness. Photocatalysts were characterised by N(2) adsorption-desorption, mercury intrusion porosimetry, SEM, UV-vis spectroscopy and XRD. A direct relationship between the TiO(2) content and the photocatalytic activity was observed up to three layers of TiO(2) (0.66 wt.%). Our results indicate that intermediate species generated on the TiO(2) layer can migrate through relatively long distances to react with the OH(-) surface groups of the support. By increasing the TiO(2) loading of the photocatalyst two effects were observed: trichloroethylene conversion is enhanced, while the efficiency of the oxidation process is decreased at expenses of increasing the selectivity to COCl(2) and dichloroacetylchloride (DCAC). The results are discussed in terms of the layer thickness, TiO(2) amount, TCE conversion and CO(2), and COCl(2) selectivity. (C) 2009 Elsevier B.V. All rights reserved.