scholarly journals Local stresses, dyke arrest and surface deformation in volcanic edificesand rift zones

2009 ◽  
Vol 47 (4) ◽  
Author(s):  
A. Gudmundsson ◽  
L. S. Brenner
Keyword(s):  
Surface Deformation ◽  
Numerical Models ◽  
Local Stress ◽  
Volcanic Eruptions ◽  
Field Studies ◽  
Normal Faults ◽  
Small Surface ◽  
Tensile Stresses ◽  
Surface Stresses ◽  
Rift Zones

Field studies indicate that nearly all eruptions in volcanic edifices and rift zones are supplied with magma through fractures (dykes) that are opened by magmatic overpressure. While (inferred) dyke injections are frequent during unrest periods, volcanic eruptions are, in comparison, infrequent, suggesting that most dykes become arrested at certain depths in the crust, in agreement with field studies. The frequency of dyke arrest can be partly explained by the numerical models presented here which indicate that volcanic edifices and rift zones consisting of rocks of contrasting mechanical properties, such as soft pyroclastic layers and stiff lava flows, commonly develop local stress fields that encourage dyke arrest. During unrest, surface deformation studies are routinely used to infer the geometries of arrested dykes, and some models (using homogeneous, isotropic half-spaces) infer large grabens to be induced by such dykes. Our results, however, show that the dyke-tip tensile stresses are normally much greater than the induced surface stresses, making it difficult to explain how a dyke can induce surface stresses in excess of the tensile (or shear) strength while the same strength is not exceeded at the (arrested) dyke tip. Also, arrested dyke tips in eroded or active rift zones are normally not associated with dyke-induced grabens or normal faults, and some dykes arrested within a few metres of the surface do not generate faults or grabens. The numerical models show that abrupt changes in Young's moduli(stiffnesses), layers with relatively high dyke-normal compressive stresses (stress barriers), and weak horizontal contacts may make the dyke-induced surface tensile stresses too small for significant fault or graben formation to occur in rift zones or volcanic edifices. Also, these small surface stresses may have no simple relation to the dyke geometry or the depth to its tip. Thus, for a layered crust with weak contacts, straightforward inversion of surface geodetic data may lead to unreliable geometries of arrested dykes in active rift zones and volcanic edifices.

Forecasting the propagation paths of fluid-driven fractures, particularly dikes and inclined sheets

2020 ◽  
Author(s):  
Agust Gudmundsson ◽  
Kyriaki Drymoni ◽  
Mohsen Bazargan ◽  
Kayode Adeoye-Akinde
Keyword(s):  
Numerical Models ◽  
Volcanic Eruptions ◽  
Pyroclastic Flows ◽  
Propagation Path ◽  
Normal Faults ◽  
Hydraulic Fractures ◽  
Pure Mode ◽  
Field Observations ◽  
Least Action ◽  
Propagation Paths

<p>It is of great importance in many fields to be able to forecast the likely propagation paths of fluid-driven factures. These include mineral veins, human-made hydraulic fractures, and dikes/inclined sheets. The physical principles that control the propagation of all fluid-driven fractures are the same. Here the focus is on dikes and inclined sheets where the selected path determines whether, where, and when a particular dike/sheet reaches the surface to erupt. Here we provide analytical and numerical models on dike/sheet paths in crustal segments (including volcanoes) that include layers of various types (lava flows, pyroclastic flows, tuff layers, soil layers, etc) as well as mechanically weak contacts and faults. The modelling results are then compared with, and tested on, actual data of two types. (a) Seismic data on the paths of dikes/sheets as well as human-made hydraulic fractures, and (b) field data on the actual propagation paths of dikes/sheets in layered and faulted rocks</p><p>The numerical results show that, particularly in stratovolcanoes, the paths are likely to be complex with common deflections along layer contacts, in agreement with field observations.  Also, some dikes/sheets may use existing faults as parts of their paths, primarily steeply dipping and recently active normal faults. The propagation path is thus not entirely in pure mode I but rather partly in a mixed mode. The energy required to propagate the dike/sheet is mainly the surface energy needed to rupture the rock, to form two new surfaces and move them apart as the fracture propagates. The energy available to drive the fracture is the stored elastic energy in the hosting crustal segment.</p><p>From its point of initiation in the magma-chamber roof, a dike/sheet can, theoretically, select any one of an infinite number of paths to follow to its point of arrest or eruption. It is shown that the eventual path selected is the one of least action, that is, the path along which the time integral of the difference between the kinetic and potential energies is an extremum (normally a minimum) relative to all other possible paths with the same endpoints. If the kinetic energy is omitted, and there are no constraints, then least action becomes the minimum potential energy, which was postulated as a basis for understanding dike propagation by Gudmundsson (1986). Here it is shown how this theoretical framework can help us make reliable forecasts of dike/sheet paths and associated volcanic eruptions.</p><p>Gudmundsson, A., 1986. Formation of dykes, feeder-dykes, and the intrusion of dykes from magma chambers. Bulletin of Volcanology, 47, 537-550.</p><p>Gudmundsson, A., 2020. Volcanotectonics: Understanding the Structure, Deformation, and Dynamics of Volcanoes. Cambridge University Press, Cambridge.</p><p>Drymoni, K., Browning, J. Gudmundsson, A., 2020. Dyke-arrest scenarios in extensional regimes: insights from field observations and numerical models, Santorini, Greece. Journal of Volcanology and Geothermal Research (in press).</p><p> </p>

Crustal folds alter local stress fields as demonstrated by magma sheet – fold interactions in the Central Andes

2021 ◽  
Author(s):  
Matías Clunes ◽  
John Browning ◽  
José Cembrano ◽  
Carlos Marquardt ◽  
Agust Gudmundsson
Keyword(s):  
Surface Deformation ◽  
Local Stress ◽  
Volcanic Eruptions ◽  
Central Andes ◽  
Stress Fields ◽  
Field Observations ◽  
Magma Chambers ◽  
The Core ◽  
Principal Compressive Stress ◽  
The Andes

<p>For magma chambers to form or volcanic eruptions to occur magma must propagate through the crust as dikes, inclined sheets and sills. The vast majority of models that investigate magma paths assume the crust to be either homogeneous or horizontally layered, often composed of rocks of contrasting mechanical properties. In subduction regions that have experienced orogenesis, like the Andes, the crust has been deformed over several million years, resulting in rock layers that are commonly folded and steeply dipping. The assumption of homogeneous properties or horizontal layering then does not capture all of the potential magma path crustal interactions. Here we tackle this problem by determining the effect of a crust made of steeply inclined layers in which sills and inclined sheets are emplaced. We combine field observations from a sill emplaced in the core of an anticlinal fold at El Juncal in the Chilean Central Andes, such as lithologies, sill and fold limbs attitude, sill length and layers and sill thickness, with a suite of finite element method models to explore the mechanical interactions between inclined layers and magma paths. Our results demonstrate that the properties of the host rock layers as well as the contacts between the layers and the crustal geometry all play an important role on magma propagation and emplacement at shallow levels. Sill propagation and emplacement through heterogeneous and anisotropic crustal segments changes the crustal stress field promoting sill arrest, deflection or propagation. Specifically, sills are more likely to be deflected when encountering shallow dipping layers rather than steeply dipping layers of a fold. Mechanically weak contacts encourage sill deflection due to the related rotation of the maximum principal compressive stress and this effect is attenuated when the fold layers are more steeply dipping. This processes may change the amount and style of surface deformation recorded, with significant implications for monitoring of active volcanoes.</p>

Optimizing Ensemble-Based Inversions for Non-unique Volcanic Systems

2020 ◽  
Author(s):  
John Albright ◽  
Patricia Gregg
Keyword(s):  
Data Assimilation ◽  
Surface Deformation ◽  
Mechanical Stability ◽  
Numerical Models ◽  
Random Noise ◽  
Source Parameters ◽  
Volcanic Eruptions ◽  
Data Availability ◽  
Analytical Models ◽  
Experimental Conditions

<p>In recent years, the advent of ensemble-based methods in volcanology has greatly facilitated the use of numerical models within data assimilation frameworks that had previously been limited, either computationally or mathematically, to simpler analytical models. Because numerical models can simulate stress conditions throughout the model space, recent inversions based on assimilated volcanic deformation data are able to track not only the basic parameters of a magma reservoir, but also how those parameters affect the overall mechanical stability of the system. Although this approach has produced successful forecasts and hind-casts of volcanic eruptions, much work remains to be done in assessing its full capabilities and limitations. In particular, non-uniqueness in how source parameters are reflected in surface deformation can significantly impair the inversion’s ability to resolve the magma system’s true state and, by extension, the likelihood of eruption. While this problem is nearly intractable for deep reservoirs, for which changes in pressure and size are indistinguishable from deformation alone, preliminary synthetic tests at shallower systems have demonstrated a limited ability to resolve the main inflation mechanism. In this study, we investigate how the performance of an Ensemble Kalman Filter (EnKF) data assimilation framework varies under a wider range of experimental conditions than used in these initial investigations. In particular, we test how different mathematical implementations of the filter and how different levels of data availability affect the EnKF’s ability to distinguish inflation drivers and to accurately resolve reservoir parameters. To implement this experiment, two time series of synthetic GPS and InSAR data are generated, one in which deformation is driven by excess pressure and another in which it is driven by lateral expansion of the reservoir. For each filter implementation these datasets are down-sampled and given random noise prior to inversion, and after assimilation the resulting model is compared to the original synthetic conditions. We find that newer deterministic formulations of the EnKF are more accurate and consistent than the original stochastic implementation, although the improvement is relatively small. Moreover, some amount of parameter inflation is required to avoid model collapse, but more sophisticated adaptive inflation schemes do not produce better results than more basic formulations. Finally, we show that while increased data sampling does improve performance, this effect is subject to diminishing returns. In particular, data resolution near the center of inflation is more important than overall range of coverage. As new inversion techniques are developed or adapted from other fields, rigorous testing as demonstrated here will be a key step in being able to interpret future results and develop new forecasting frameworks for volcanic eruptions.</p>

Influence of viscous layers on the growth of normal faults: insights from experimental and numerical models

2003 ◽  
Vol 25 (9) ◽  
pp. 1471-1485 ◽  
Author(s):  
Nicolas Bellahsen ◽  
Jean-Marc Daniel ◽  
Laurent Bollinger ◽  
Evgenii Burov
Keyword(s):  
Numerical Models ◽  
Normal Faults

Ensemble-based data assimilation of volcanic aerosols using FALL3D+PDAF

2021 ◽  
Author(s):  
Leonardo Mingari ◽  
Andrew Prata ◽  
Federica Pardini
Keyword(s):  
Data Assimilation ◽  
Information Transfer ◽  
High Performance ◽  
Volcanic Ash ◽  
Numerical Models ◽  
Volcanic Eruptions ◽  
Column Height ◽  
Eruption Column ◽  
Operational Environment ◽  
Assimilation System

<p>Modelling atmospheric dispersion and deposition of volcanic ash is becoming increasingly valuable for understanding the potential impacts of explosive volcanic eruptions on infrastructures, air quality and aviation. The generation of high-resolution forecasts depends on the accuracy and reliability of the input data for models. Uncertainties in key parameters such as eruption column height injection, physical properties of particles or meteorological fields, represent a major source of error in forecasting airborne volcanic ash. The availability of nearly real time geostationary satellite observations with high spatial and temporal resolutions provides the opportunity to improve forecasts in an operational context. Data assimilation (DA) is one of the most effective ways to reduce the error associated with the forecasts through the incorporation of available observations into numerical models. Here we present a new implementation of an ensemble-based data assimilation system based on the coupling between the FALL3D dispersal model and the Parallel Data Assimilation Framework (PDAF). The implementation is based on the last version release of FALL3D (versions 8.x) tailored to the extreme-scale computing requirements, which has been redesigned and rewritten from scratch in the framework of the EU Center of Excellence for Exascale in Solid Earth (ChEESE). The proposed methodology can be efficiently implemented in an operational environment by exploiting high-performance computing (HPC) resources. The FALL3D+PDAF system can be run in parallel and supports online-coupled DA, which allows an efficient information transfer through parallel communication. Satellite-retrieved data from recent volcanic eruptions were considered as input observations for the assimilation system.</p>

Modelling of interactions between dykes, inclined sheets and faults at Santorini volcano

2021 ◽  
Author(s):  
Kyriaki Drymoni ◽  
John Browning ◽  
Agust Gudmundsson
Keyword(s):  
Tensile Strength ◽  
Fault Zone ◽  
Numerical Models ◽  
Damage Zone ◽  
Sheet Thickness ◽  
Analytical Models ◽  
Energy Conditions ◽  
Dyke Swarm ◽  
Normal Faults ◽  
Santorini Volcano

<p>Dykes and inclined sheets are known occasionally to exploit faults as parts of their paths, but the conditions that allow this to happen are still not fully understood. Here we report field observations from a well-exposed dyke swarm of the Santorini volcano, Greece, that show dykes and inclined sheets deflected into faults and the results of analytical and numerical models to explain the conditions for deflection. The deflected dykes and sheets belong to a local swarm of 91 dyke/sheet segments that was emplaced in a highly heterogeneous and anisotropic host rock and partially cut by some regional faults and a series of historic caldera collapses, the caldera walls providing, excellent exposures of the structures. The numerical models focus on a normal-fault dipping 65° with a damage zone composed of parallel layers or zones of progressively more compliant rocks with increasing distance from the fault rupture plane. We model sheet-intrusions dipping from 0˚ to 90˚ and with overpressures of alternatively 1 MPa and 5 MPa, approaching the fault. We further tested the effects of changing (1) the sheet thickness, (2) the fault-zone thickness, (3) the fault-zone dip-dimension (height), and (4) the loading by, alternatively, regional extension and compression. We find that the stiffness of the fault core, where a compliant core characterises recently active fault zones, has pronounced effects on the orientation and magnitudes of the local stresses and, thereby, on the likelihood of dyke/sheet deflection into the fault zone. Similarly, the analytical models, focusing on the fault-zone tensile strength and energy conditions for dyke/sheet deflection, indicate that dykes/sheets are most likely to be deflected into and use steeply dipping recently active (zero tensile-strength) normal faults as parts of their paths.</p>

A model for Archean tectonism. Part 2. Numerical models of vertical tectonism in greenstone belts

1980 ◽  
Vol 17 (1) ◽  
pp. 60-71 ◽  
Author(s):  
Jean-Claude Mareschal ◽  
Gordon F. West
Keyword(s):  
Numerical Models ◽  
Volcanic Eruptions ◽  
Crustal Rocks ◽  
Tectonic Model ◽  
Mechanical Evolution ◽  
Accumulated Strain ◽  
Lava Sequence ◽  
Order Of Magnitude ◽  
Greenstone Belts ◽  
Horizontal Temperature Gradients

A tectonic model that attempts to explain common features of Archean geology is investigated. The model supposes the accumulation, by volcanic eruptions, of a thick basaltic pile on a granitoid crust. The thermal blanketing effect of this lava raises the temperature of the granitic crust and eventually softens it enough that gravitational slumping and downfolding of the lava follows.Numerical models of the thermal and mechanical evolution of a granitoid crust covered with a thick lava sequence indicate that such an evolution is possible when reasonable assumptions are made about the temperature dependence of the viscosity in crustal rocks. These models show the lava sinking in relatively narrow regions while wider granite diapirs appear in between. The convection produces strong horizontal temperature gradients that may cause lateral changes in metamoprhic facies. A one order of magnitude drop in accumulated strain occurs between the granite–basalt interface and the center of the granite diaper at a depth of 10–15 km.

scholarly journals Looking beyond kinematics: 3D thermo-mechanical modelling reveals the dynamic of transform margins

2021 ◽  
Author(s):  
Anthony Jourdon ◽  
Charlie Kergaravat ◽  
Guillaume Duclaux ◽  
Caroline Huguen
Keyword(s):  
Boundary Conditions ◽  
Strain Localization ◽  
Numerical Models ◽  
Local Stress ◽  
Shear Zones ◽  
Strike Slip ◽  
Convergent Margins ◽  
Mechanical Modelling ◽  
Kinematic Models ◽  
Extension Direction

Abstract. Transform margins represent ~30 % of the non-convergent margins worldwide. Their formation and evolution have long been addressed through kinematic models that do not account for the mechanical behaviour of the lithosphere. In this study, we use high resolution 3D numerical thermo-mechanical modelling to simulate and investigate the evolution of the intra-continental strain localization under oblique extension. The obliquity is set through velocity boundary conditions that range from 15° (high obliquity) to 75° (low obliquity) every 15° for strong and weak lower continental crust rheologies. Numerical models show that the formation of localized strike-slip shear zones leading to transform continental margins always follows a thinning phase during which the lithosphere is thermally and mechanically weakened. For low (75°) to intermediate (45°) obliquity cases, the strike-slip faults are not parallel to the extension direction but form an angle of 20° to 40° with the plates' motion while for higher obliquities (30° to 15°) the strike-slip faults develop parallel to the extension direction. Numerical models also show that during the thinning of the lithosphere, the stress and strain re-orient while boundary conditions are kept constant. This evolution, due to the weakening of the lithosphere, leads to a strain localization process in three major phases: (1) strain initiates in a rigid plate where structures are sub-perpendicular to the extension direction; (2) distributed deformation with local stress field variations and formation of transtensional and strike-slip structures; (3) formation of highly localized plates boundaries stopping the intra-continental deformation. Our results call for a thorough re-evaluation of the kinematic approach to studying transform margins.

A study of microseismicity in Northern Baja California, Mexico

1976 ◽  
Vol 66 (6) ◽  
pp. 1921-1929 ◽  
Author(s):  
Tracy L. Johnson ◽  
Juan Madrid ◽  
Theodore Koczynski
Keyword(s):  
Fault Zone ◽  
Focal Mechanism ◽  
Baja California ◽  
Surface Deformation ◽  
Normal Faults ◽  
East Side ◽  
Composite Focal Mechanism ◽  
Agua Blanca Fault ◽  
Eastern Edge ◽  
Low Activity

abstract Five microearthquake instruments were operated for 2 months in 1974 in a small mobile array deployed at various sites near the Agua Blanca and San Miguel faults. An 80-km-long dection of the San Miguel fault zone is presently active seismically, producing the vast majority of recorded earthquakes. Very low activity was recorded on the Agua Blanca fault. Events were also located near normal faults forming the eastern edge of the Sierra Juarez suggesting that these faults are active. Hypocenters on the San Miguel fault range in depth from 0 to 20 km although two-thirds are in the upper 10 km. A composite focal mechanism showing a mixture of right-lateral and dip slip, east side up, is similar to a solution obtained for the 1956 San Miguel earthquake which proved consistent with observed surface deformation.

scholarly journals Aeolian Remobilisation of Volcanic Ash: Outcomes of a Workshop in the Argentinian Patagonia

2020 ◽  
Vol 8 ◽  
Author(s):  
Paul A. Jarvis ◽  
Costanza Bonadonna ◽  
Lucia Dominguez ◽  
Pablo Forte ◽  
Corine Frischknecht ◽  
...  
Keyword(s):  
Volcanic Ash ◽  
Friction Velocity ◽  
Critical Infrastructure ◽  
Dust Emission ◽  
Source Area ◽  
Field Measurements ◽  
Volcanic Eruptions ◽  
Field Studies ◽  
Sediment Traps ◽  
Economic Sectors

During explosive volcanic eruptions, large quantities of tephra can be dispersed and deposited over wide areas. Following deposition, subsequent aeolian remobilisation of ash can potentially exacerbate primary impacts on timescales of months to millennia. Recent ash remobilisation events (e.g., following eruptions of Cordón Caulle 2011; Chile, and Eyjafjallajökull 2010, Iceland) have highlighted this to be a recurring phenomenon with consequences for human health, economic sectors, and critical infrastructure. Consequently, scientists from observatories and Volcanic Ash Advisory Centers (VAACs), as well as researchers from fields including volcanology, aeolian processes and soil sciences, convened at the San Carlos de Bariloche headquarters of the Argentinian National Institute of Agricultural Technology to discuss the “state of the art” for field studies of remobilised deposits as well as monitoring, modeling and understanding ash remobilisation. In this article, we identify practices for field characterisation of deposits and active processes, including mapping, particle characterisation and sediment traps. Furthermore, since forecast models currently rely on poorly-constrained dust emission schemes, we call for laboratory and field measurements to better parameterise the flux of volcanic ash as a function of friction velocity. While source area location and extent are currently the primary inputs for dispersion models, once emission schemes become more sophisticated and better constrained, other parameters will also become important (e.g., source material volume and properties, effective precipitation, type and distribution of vegetation cover, friction velocity). Thus, aeolian ash remobilisation hazard and associated impact assessment require systematic monitoring, including the development of a regularly-updated spatial database of resuspension source areas.

Sign in / Sign up

Export Citation Format

Share Document