scholarly journals Syn-Eruptive Conditions of the AD 1530 Sub-Plinian Eruption of La Soufrière of Guadeloupe (Lesser Antilles)

2021 ◽  
Vol 9 ◽  
Author(s):  
Caroline Martel ◽  
Michel Pichavant ◽  
Hélène Balcone-Boissard ◽  
Georges Boudon

A prerequisite in refining volcanic hazard at explosive volcanoes is a better quantification of the timescales of the syn-eruptive processes, such as magma degassing and crystallization prior to eruption. To this aim, new data on the matrices (microlites, residual glasses, and bubbles) of pumice, scoria, and dense clasts erupted during the AD 1530 andesitic eruption of La Soufrière of Guadeloupe are combined with published data from phase-equilibrium and kinetic experiments, in order to estimate pressures of microlite crystallization and magma ascent rates. From the timescale data, we infer that the AD 1530 eruption started with phreatomagmatic explosions tapping magmas that ascended during about 1 month (decompression rate of ∼50 Pa/s) from the coldest parts of the reservoir (∼825°C and a 74 wt% SiO2 melt). These magmas continuously crystallized microlites (∼25 vol% plagioclase, orthopyroxene, clinopyroxene, magnetite, quartz, and apatite), but did not outgas likely due to sealed conduit margins. The flank collapse (preexisting “cold” edifice) that followed the phreatomagmatic phase triggered a sub-Plinian eruption that progressively tapped the hotter main reservoir (∼875°C and 71 wt% SiO2 interstitial melt), emitting banded and homogeneous pumice. The banded pumice did not significantly outgas and mostly lack decompression-induced microlites, suggesting short ascent durations of the order of 0.5–1 day (decompression rates of 1,400–4,000 Pa/s). The following Strombolian phase emitted dark scoria that did not significantly outgas and only crystallized rare microlites, suggesting ascent duration of the order of 2 days (decompression rates of ∼550 Pa/s). The terminal lava dome growth involved fully outgassed magmas ascended during more than 1 month, giving time for microlite crystallization (∼40 vol% plagioclase, orthopyroxene, clinopyroxene, magnetite, and cristobalite). The detection of any shallow new magmatic intrusion is crucial, since it can trigger a sequence of conduit processes leading to an eruption marked by a succession of different and disastrous eruptive styles, following a scenario similar to the AD 1530 eruption. Overall, we provide a combined approach of petrological, geochemical, and experimental data that may be used to infer ascent conditions and rates at other volcanic systems.

Solid Earth ◽  
2013 ◽  
Vol 4 (2) ◽  
pp. 201-213 ◽  
Author(s):  
S. B. Mueller ◽  
N. R. Varley ◽  
U. Kueppers ◽  
P. Lesage ◽  
G. Á. Reyes Davila ◽  
...  

Abstract. The most recent eruptive phase of Volcán de Colima, Mexico, started in 1998 and was characterized by dome growth with a variable effusion rate, interrupted intermittently by explosive eruptions. Between November 2009 and June 2011, activity at the dome was mostly limited to a lobe on the western side where it had previously started overflowing the crater rim, leading to the generation of rockfall events. As a consequence of this, no significant increase in dome volume was perceivable and the rate of magma ascent, a crucial parameter for volcano monitoring and hazard assessment could no longer be quantified via measurements of the dome's dimensions. Here, we present alternative approaches to quantify the magma ascent rate. We estimate the volume of individual rockfalls through the detailed analysis of sets of photographs (before and after individual rockfall events). The relationship between volume and infrared images of the freshly exposed dome surface and the seismic signals related to the rockfall events were then investigated. Larger rockfall events exhibited a correlation between its previously estimated volume and the surface temperature of the freshly exposed dome surface, as well as the mean temperature of rockfall mass distributed over the slope. We showed that for larger events, the volume of the rockfall correlates with the maximum temperature of the newly exposed lava dome as well as a proxy for seismic energy. It was therefore possible to calibrate the seismic signals using the volumes estimated from photographs and the count of rockfalls over a certain period was used to estimate the magma extrusion flux for the period investigated. Over the course of the measurement period, significant changes were observed in number of rockfalls, rockfall volume and hence averaged extrusion rate. The extrusion rate was not constant: it increased from 0.008 ± 0.003 to 0.02 ± 0.007 m3 s−1 during 2010 and dropped down to 0.008 ± 0.003 m3 s−1 again in March 2011. In June 2011, magma extrusion had come to a halt. The methodology presented represents a reliable tool to constrain the growth rate of domes that are repeatedly affected by partial collapses. There is a good correlation between thermal and seismic energies and rockfall volume. Thus it is possible to calibrate the seismic records associated with the rockfalls (a continuous monitoring tool) to improve volcano monitoring at volcanoes with active dome growth.


2010 ◽  
Vol 10 (6) ◽  
pp. 1281-1292 ◽  
Author(s):  
B. Poisson ◽  
R. Pedreros

Abstract. Two historical landslide-induced tsunamis that reached the coasts of the French Lesser Antilles are studied. First, the Martinique coast was hit by a tsunami down the western flank of Montagne Pelée at the beginning of the big eruption of May 1902. More recently, the northeastern coast of Guadeloupe was affected by a tsunami that had been generated around Montserrat by pyroclastic flows entering the sea, during the July 2003 eruption of the Soufrière Hills volcano. We use a modified version of the GEOWAVE model to compute numerical simulations of both events. Two source hypotheses are considered for each tsunami. The comparison of the simulation results with reported tsunami height data helps to discriminate between the tested source decriptions. In the Martinique case, we obtain a better fit to data when considering three successive lahars entering the sea, as a simplified single source leads to an overstimation of the tsunami wave heights at the coast. In the Montserrat case, the best model uses a unique source which volume corresponds to published data concerning the peak volume flow. These findings emphasize the importance of an accurate description of the relevant volume as well as the timing sequence of the source event in landslide-generated tsunami modelling. They also show that considering far-field effects in addition to near-field effects may significantly improve tsunami modelling.


CFD letters ◽  
2021 ◽  
Vol 13 (9) ◽  
pp. 13-27
Author(s):  
Mohamad Lutfi Samsudin ◽  
Hasril Hasini

Meshing of domain in CFD is an important step to ensure accuracy of the solution. In the past, hexahedral or tetrahedral mesh systems were commonly used, and both have their merits and demerits. For large and complex geometry, polyhedral is another option but its accuracy is claimed to be lacking. In this paper, the use of polyhedral mesh system by past researchers are reviewed. Evaluation on the application of polyhedral mesh system for the study of the vortex formation with a simple single pump sump model is made. Validation was made through the comparison of the results from hexahedral, tetrahedral and polyhedral mesh sizes and the experimental data from published data. The polyhedral mesh system was found to perform satisfactorily and was able to match the results from the hexahedral mesh system as well as the experimental data.


2008 ◽  
Vol 74 (3) ◽  
pp. 391-429 ◽  
Author(s):  
G. HERDRICH ◽  
D. PETKOW

AbstractThe development of the inductively driven plasma wind tunnel PWK3, which enables the electrodeless generation of high-enthalpy plasmas for the development of heat shield materials required for space vehicles performing entry manoeuvres in the atmospheres of Venus, Earth and Mars, is described. The facility with its modular inductive plasma generators allows operation with gases such as carbon dioxide, air, oxygen and nitrogen and was qualified for thermal plasma powers up to 60 kW. Previously developed models for determining plasma properties and plasma source related characteristics enable a maximum plasma power in combination with long operational periods using different operational gases and gas mixtures. This is achieved by an optimization using the optimum operational frequency, a minimization of field losses using very thin plasma tube wall thicknesses and the successful application of MHD effects. Based on the solved cylinder problem for ICPs, a one-dimensional model for radial Lorentz forces and magnetic pressure has been developed. Here, a synthesis of previously published data and works is made where the new algebraic model for the calculation of Lorentz forces and magnetic pressures in an ICP was used and applied to experimental data. In addition, results from the model using the experimental data are shown to be consistent and, in addition, a comparison with a simpler model based on the well-known exponential approach for ICPs showed that the simpler model is covered without fail by the new model. The new model also states that there is a maximum of the Lorentz forces over the damping parameter d/δ (plasma diameter divided by skin depth) which almost corresponds with the position of the maximum plasma power of the cylindric model for ICPs. For the magnetic pressure the position of the maximum pressure is identical to the value for d/δ for the maximum plasma power.


2021 ◽  
Author(s):  
Conner Sarich ◽  
Adam Hope ◽  
Jim Rule

Abstract Precipitation kinetics were investigated in select Fe, Ni, and Al alloys using a CALPHAD based precipitation model based on Langer-Schwartz theory. Thermodynamic and kinetic data are taken from commercially available CALPHAD software, but reliable interfacial energy data for precipitates needed for the calculations is often lacking. While models exist to approximate these interfacial energies, this study has focused on deriving more reliable estimates by comparison with experimental data. By performing simulations with thermal histories, nucleation sites, and precipitate morphologies that closely replicate experimental data found in literature, the interfacial energies were optimized until volume fraction and mean radius values closely matched the published data. Using this technique, interfacial energy values have been determined for carbides in Grade 22 low alloy steels, delta phase in Ni 625 and 718, SPhase in Al 2024, and Q’ and β’’ in Al 6111, and can be used for future predictive precipitation simulations.


2002 ◽  
Vol 76 (4) ◽  
pp. 349-353 ◽  
Author(s):  
K.N. Mouritsen

AbstractThe release ofMaritrema subdolumcercariae (Digenea: Microphallidae) from the marine mud snailHydrobia ulvaeis significantly affected by temperature, salinity, light and exudates from the second intermediate amphipod host. Based on (i) previously published data on temperature–salinity dependentH. ulvaeactivity, (ii) new experimental data onH. ulvaeactivity in light and darkness as well as in the presence and absence of host exudates, and (iii) the cercarial emergence rate from free moving snails and snails prevented from crawling, the present analysis indicates that emergence ofM. subdolumlarvae is regulated mainly by host activity as the ultimate factor for release. The adaptive significance of such an emergence strategy is emphasized.


Materials ◽  
2019 ◽  
Vol 12 (17) ◽  
pp. 2728 ◽  
Author(s):  
Sergey V. Ushakov ◽  
Alexandra Navrotsky ◽  
Qi-Jun Hong ◽  
Axel van de Walle

Among transition metal carbides and nitrides, zirconium, and hafnium compounds are the most stable and have the highest melting temperatures. Here we review published data on phases and phase equilibria in Hf-Zr-C-N-O system, from experiment and ab initio computations with focus on rocksalt Zr and Hf carbides and nitrides, their solid solutions and oxygen solubility limits. The systematic experimental studies on phase equilibria and thermodynamics were performed mainly 40–60 years ago, mostly for binary systems of Zr and Hf with C and N. Since then, synthesis of several oxynitrides was reported in the fluorite-derivative type of structures, of orthorhombic and cubic higher nitrides Zr3N4 and Hf3N4. An ever-increasing stream of data is provided by ab initio computations, and one of the testable predictions is that the rocksalt HfC0.75N0.22 phase would have the highest known melting temperature. Experimental data on melting temperatures of hafnium carbonitrides are absent, but minimum in heat capacity and maximum in hardness were reported for Hf(C,N) solid solutions. New methods, such as electrical pulse heating and laser melting, can fill the gaps in experimental data and validate ab initio predictions.


Author(s):  
Robert A. Basterfield ◽  
Chris J. Lawrence ◽  
Michael J. Adams

Pastes occur as intermediates or final product forms in many industrially important manufacturing sectors. The use of computer simulation techniques, such as the finite element method, is becoming more common in the design of paste processing operations. A major problem in the application of this approach is the development of sufficiently representative materials models. It has been established that pastes may be described as elasto-viscoplastic materials with the plastic flow being governed by the Herschel-Bulkley relationship. This paper describes the development of analytical and numerical models that can be used as a basis for deriving the material parameters from experimental data obtained using extrusion, compression and bending procedures. Measurements have also been carried out on a model paste and the derived material parameters are compared with published data for the same paste. The merits of the three experimental methods are compared on this basis.


2010 ◽  
Vol 133 (3) ◽  
Author(s):  
J. Michael Owen

Ingress of hot gas through the rim seals of gas turbines can be modeled theoretically using the so-called orifice equations. In Part I of this two-part paper, the orifice equations were derived for compressible and incompressible swirling flows, and the incompressible equations were solved for axisymmetric rotationally induced (RI) ingress. In Part II, the incompressible equations are solved for nonaxisymmetric externally induced (EI) ingress and for combined EI and RI ingress. The solutions show how the nondimensional ingress and egress flow rates vary with Θ0, the ratio of the flow rate of sealing air to the flow rate necessary to prevent ingress. For EI ingress, a “saw-tooth model” is used for the circumferential variation of pressure in the external annulus, and it is shown that ε, the sealing effectiveness, depends principally on Θ0; the theoretical variation of ε with Θ0 is similar to that found in Part I for RI ingress. For combined ingress, the solution of the orifice equations shows the transition from RI to EI ingress as the amplitude of the circumferential variation of pressure increases. The predicted values of ε for EI ingress are in good agreement with the available experimental data, but there are insufficient published data to validate the theory for combined ingress.


Sign in / Sign up

Export Citation Format

Share Document