Syn-eruptive breakdown of pyrrhotite: a record of magma fragmentation,air entrainment,and oxidation

Air entrainment in fragmented magmas controls the dynamics of volcanic eruptions. Pyroclast oxidation kinetics may be applied to quantify the degree of magma–air interaction. Pyrrhotite (Po) in volcanic rocks is often oxidized to form magnetite (Mt) and hematite (Hm), and its reaction mechanisms are well constrained. To test utilizing Po oxidation as a marker for magma–air interactions, we compared the occurrence of Po oxidation products from three different eruption styles during the Sakurajima 1914–1915 eruption. Pumices from the Plinian eruption include columnar-type Fe oxides (Mt with subordinate width of Hm) often accompanied by relict Po. This columnar type is also found in clastogenic lava, where it is almost completely oxidized to Hm. The effusive lava contains framboidal aggregates of subhedral to anhedral Mt crystals without Hm. The formation mechanisms of columnar and framboidal Fe oxides were estimated. The columnar type Fe oxides were formed syn-eruptively through gaseous reactions, as opposed to the melt in a magma chamber, as demonstrated by the Ti-free nature of the columnar Mt and its synchronous oxidation to Hm. By contrast, the framboidal type was formed in a melt with decreasing fS₂. The calculation of Hm growth in a conductively cooling pumice clast constrains the surface temperature of pumice in the eruption column. The paragenesis and oxidation degree of Po and Fe oxides are consistent with the eruption processes in terms of magma fragmentation, air entrainment, and welding, and can, therefore, be a responsive marker for the magma–air interaction.     

The ground layer of Ata pyroclastic flow deposit,southwestern Japan — evidence for the capture of lithic fragments

Lithic fragments in the ground layer of the Ata pyroclastic flow deposit, southwestern Japan, were supplied from two different sources. One is the eruptive vent and the other is the basement rock exposed underneath the path of flow. Lithic fragments captured at the eruptive vent gradually decrease in size with distance from the source. Local increases of ML or Md are proportional to increased amounts of captured lithic fragments. The pyroclastic flow eroded basement formations on slopes dipping away from the source, and deposited the lithics within the ground layer on slopes dipping towards the source. The ground layer was found only in the western half of the Ata pyroclastic flow deposit. The absence of the ground layer in the eastern half of the pyroclastic flow deposit is interpreted to result from a selective loss of lithics when the flow traversed a bay or a lake located just east from the vent.     

Mixing efficiency in decaying stably stratified turbulence

Mixing efficiency in stratified flows is a measure of the proportion of turbulent kinetic energy that goes into increasing the potential energy of the fluid by irreversible mixing. In this research direct numerical simulations (DNS) and rapid distortion theory (RDT) calculations of transient turbulent mixing events are carried out in order to study this aspect of mixing. In particular, DNS and RDT of decaying, homogeneous, stably-stratified turbulence are used to determine the mixing efficiency as a function of the initial turbulence Richardson number R_it0=(N_L0/_u0)²$R{i}_{t0}={(N{L}_0/u_0)}^2$, where N   is the buoyancy frequency and _L0$L_0$ and _u0$u_0$ are initial length and velocity scales of the turbulence. The results show that the mixing efficiency increases with increasing R_it0$R{i}_{t0}$ for small R_it0$R{i}_{t0}$, but for larger R_it0$R{i}_{t0}$ the mixing efficiency becomes approximately constant. These results are compared with data from towed grid experiments. There is qualitative agreement between the DNS results and the available experimental data, but significant quantitative discrepancies. The grid turbulence experiments suggest a maximum mixing efficiency (at large R_it0$R{i}_{t0}$) of about 6%, while the DNS and RDT results give about 30%. We consider two possible reasons for this discrepancy: Prandtl number effects and non-matching initial conditions. We conclude that the main source of the disagreement probably is due to inaccuracy in determining the initial turbulence energy input in the case of the grid turbulence experiments.