Deep submarine pyroclastic eruptions: theory and predicted landforms and deposits

Submarine pyroclastic eruptions at depths greater than a few hundred meters are generally considered to be rare or absent because the pressure of the overlying water column is sufficient to suppress juvenile gas exsolution so that magmatic disruption and pyroclastic activity do not occur. Consideration of detailed models of the ascent and eruption of magma in a range of sea floor environments shows, however, that significant pyroclastic activity can occur even at depths in excess of 3000 m. In order to document and illustrate the full range of submarine eruption styles, we model several possible scenarios for the ascent and eruption of magma feeding submarine eruptions: (1) no gas exsolution; (2) gas exsolution but no magma disruption; (3) gas exsolution, magma disruption, and hawaiian-style fountaining; (4) volatile content builds up in the magma reservoir leading to hawaiian eruptions resulting from foam collapse; (5) magma volatile content insufficient to cause fragmentation normally but low rise speed results in strombolian activity; and (6) volatile content builds up in the top of a dike leading to vulcanian eruptions. We also examine the role of bulk-interaction steam explosivity and contact-surface steam explosivity as processes contributing to volcaniclastic formation in these environments. We concur with most earlier workers that for magma compositions typical of spreading centers and their vicinities, the most likely circumstance is the quiet effusion of magma with minor gas exsolution, and the production of somewhat vesicular pillow lavas or sheet flows, depending on effusion rate. The amounts by which magma would overshoot the vent in these types of eruptions would be insufficient to cause any magma disruption. The most likely mechanism of production of pyroclastic deposits in this environment is strombolian activity, due to the localized concentration of volatiles in magma that has a low rise rate; magmatic gas collects by bubble coalescence, and ascends in large isolated bubbles which disrupt the magma surface in the vent, producing localized blocks, bombs, and pyroclastic deposits. Another possible mode of occurrence of pyroclastic deposits results from vulcanian eruptions; these deposits, being characterized by the dominance of angular blocks of country rocks deposited in the vicinity of a crater, should be easily distinguishable from strombolian and hawaiian eruptions. However, we stress that a special case of the hawaiian eruption style is likely to occur in the submarine environment if magmatic gas buildup occurs in a magma reservoir by the upward drift of gas bubbles. In this case, a layer of foam will build up at the top of the reservoir in a sufficient concentration to exceed the volatile content necessary for disruption and hawaiian-style activity; the deposits and landforms are predicted to be somewhat different from those of a typical primary magmatic volatile-induced hawaiian eruption. Specifically, typical pyroclast sizes might be smaller; fountain heights may exceed those expected for the purely magmatic hawaiian case; cooling of descending pyroclasts would be more efficient, leading to different types of proximal deposits; and runout distances for density flows would be greater, potentially leading to submarine pyroclastic deposits surrounding vents out to distances of tens of meters to a kilometer. In addition, flows emerging after the evacuation of the foam layer would tend to be very depleted in volatiles, and thus extremely poor in vesicles relative to typical flows associated with hawaiian-style eruptions in the primary magmatic gas case. We examine several cases of reported submarine volcaniclastic deposits found at depths as great as 3000 m and conclude that submarine hawaiian and strombolian eruptions are much more common than previously suspected at mid-ocean ridges. Furthermore, the latter stages of development of volcanic edifices (seamounts) formed in submarine environments are excellent candidates for a wide range of submarine pyroclastic activity due not just to the effects of decreasing water depth, but also to: (1) the presence of a summit magma reservoir, which favors the buildup of magmatic foams (enhancing hawaiian-style activity) and episodic dike emplacement (which favors strombolian-style eruptions); and (2) the common occurrence of alkalic basalts, the CO₂ contents of which favor submarine explosive eruptions at depths greater than tholeiitic basalts. These models and predictions can be tested with future sampling and analysis programs and we provide a checklist of key observations to help distinguish among the eruption styles.     

Intra-oceanic settings of the western Mexico Late Jurassic-Early Cretaceous arc sequences. Implications for the Pacific-Tethys geodynamic relationships during the Cretaceous

Abstract

The Guerrero suspect terrane composed of Late Jurassic-Early Cretaceous sequences, extends from Baja California up to Acapulco and is considered to be coeval with the Late Mesozoic igneous and sedimentary arc sequences of the Greater Antilles, Venezuela and Western Cordillera of Colombia. New geological, petrological and geochemical data from central and southern Mexico, led us to propose a new model for the building of the Alisitos-Teloloapan arc. This arc, partly built on the Pacific oceanic lithosphere and partly on continental fragments, could be related to the subduction of an oceanic basin - the Arperos basin - under the Paleo-Pacific plate. This subduction was dipping southwest.

At the beginning of the magmatic activity of the oceanic segment of this arc, depleted tholeiitic basalts were emitted in a submarine environnement below the CCD. While subduction was going on, the arc magmas evolved from LREE depleted tholeiites to slightly LREE enriched tholeiites and then, to calc-alkaline basalts and andesites enriched in LREE and HFSE. Concurrently, the arc sedimentary environment changed from deep oceanic to neritic with the deposition of Aptian-Albian reefal limestones, at the end of the arc building. In the continent-based segment, the arc magmas are exclusively differentiated calc-alkaline suites depleted in HREE and Y, formed of predominantly siliceous lavas and pyroclastic rocks, emitted in a sub-aerial or shallow marine environment.

Thus, taking into account this above mentioned model, the Cretaceous volcanic series, accreted to the margins of cratonal America, in Colombia, Venezuela, Greater Antilles and Mexico, could be related to the same west-south-west dipping subduction of oceanic basins, fringing the North and South American continental cratons and connected directly with the inter-American Tethys. While the subduction was proceeding, this magmatic arc drifted towards the North and South American cratons and finally, collided with the continental margins at different periods during the Cretaceous.     

Solomon Sea circulation and water mass modifications: response at ENSO timescales

The South Pacific low latitude western boundary currents (LLWBCs) carry waters of subtropical origin through the Solomon Sea before joining the equatorial Pacific. Changes in their properties or transport are assumed to impact El Niño Southern Oscillation (ENSO) dynamics. At ENSO timescales, the LLWBCs transport tends to counterbalance the interior geostrophic one. When transiting through the complex geography of the Solomon Sea, the main LLWBC, the New Guinea Coastal Undercurrent, cannot follow a unique simple route to the equator. Instead, its routes and water mass properties are influenced by the circulation occurring in the Solomon Sea. In this study, the response of the Solomon Sea circulation to ENSO is investigated based on a numerical simulation. The transport anomalies entering the Solomon Sea from the south are confined to the top 250 m of the water column, where they represent 7.5 Sv (based on ENSO composites) for a mean transport of 10 Sv. The induced circulation anomalies in the Solomon Sea are not symmetric between the two ENSO states because of (1) a bathymetric control at Vitiaz Strait, which plays a stronger role during El Niño, and (2) an additional inflow through Solomon Strait during La Niña events. In terms of temperature and salinity, modifications are particularly notable for the thermocline water during El Niño conditions, with cooler and fresher waters compared to the climatological mean. The surface water at Vitiaz Strait and the upper thermocline water at Solomon Strait, feeding respectively the equatorial Pacific warm pool and the Equatorial Undercurrent, particularly affect the heat and salt fluxes. These fluxes can change by up to a factor of 2 between extreme El Niño and La Niña conditions.     

Morphology,structure and kinematics of a rainfall controlled slow‐moving Andean landslide,Peru

The large slow‐moving landslide of Maca is located in the upper Colca valley (southern Peru), a region characterized by a well pronounced rainy period, and intense and recurrent sustained seismicity. The landslide, developed in deep lacustrine deposits, has recently accelerated, threatening the Maca village. This work aims at understanding the rupture mechanism and the causes of the recent landslide reactivation/acceleration. We present a multidisciplinary characterization of the Maca landslide that includes: (i) geological and morphological mapping in the field; (ii) remote sensing analysis using an historical aerial photograph of 1955 and the Pléiades satellite images (2013); (iii) global positioning system (GPS) including time‐series of surveys over 13 years, and continuous measurements over 14 months; (iv) a geophysical campaign with deep electrical resistivity tomography profiles acquired across the landslide mass. Our study shows that this 60 Mm³ landslide, which can be classified as a clay/silt compound landslide, moved by 15 m between 2001 and 2014 with a large inter‐annual velocity variation (up to a factor of 500) depending on the rainfall intensity. We suggest that these dramatic changes in velocity are the result of the combination of a threshold mechanism and the short intense rainy season in Peru. This study reveals three main driving factors acting at different timescales: (i) over several decades, the river course has significantly changed, causing the Maca landslide reactivation in the 1980s due to the erosion of its toe; (ii) at the year scale, a minimum amount of rainfall is required to trigger the motion and this amount controls the landslide velocity; (iii) transient changes in slide velocity may occur anytime due to earthquakes. This study particularly highlights the non‐linear behaviour of the motion with rainfall. Copyright © 2016 John Wiley & Sons, Ltd.     

Investigation and sensitivity analysis of a mechanistic phytoplankton model implemented in a new modular numerical tool (Eco3M) dedicated to biogeochemical modelling

A new class of phytoplankton models with a mechanistic basis has been presented in a companion paper (Baklouti, M., Diaz, F., Pinazo, C., Faure, V., Queguiner, B., 2006. Investigation of mechanistic formulations depicting phytoplankton dynamics for models of marine pelagic ecosystems. Progress in Oceanography). It is the default class of models implemented in our new numerical tool Eco3M, which is dedicated to Ecological, Mechanistic and Modular Modelling. A brief overview of its main features is given in Section 2 of the present paper. In the next sections, a particular phytoplankton model among the aforementioned class has been tested with special emphasis on the mechanistic photosynthesis component relating the photosynthetic rate to the proportion of open photosystems II. The present study encompasses several essential steps that are inherent to any modelling, including model reduction, model sensitivity analysis and comparison of model outputs with experiments. The global sensitivity analysis of the plankton model for one-at-a-time parameter perturbations revealed a restricted set of parameters having major influence on the model outputs. Sensitivity tests involving simultaneous parameter perturbations within the range actually encountered in the literature provided a confidence interval for the outputs. Chemostat experiments performed on nitrate-limited diatoms grown under low (LL) and high-light (HL) conditions have been used for comparison with model outputs. The good fit between measured data and model outputs using the same parameter values in both the LL and HL cases demonstrates the ability of our model to represent the main features of phytoplankton dynamics including photoacclimation. Finally, Eco3M is ultimately intended to include explicit bacterial and zooplankton compartments, as well as to be coupled with ocean circulation models, but the intrinsic behavior of the phytoplankton model has been investigated first, independently of physical forcing.     

Erosive flood events on the surface of Mars: Application to Mangala and Athabasca Valles

We address key factors involved in determining water flow conditions in outflow channels on Mars, including the temperature of the sub-surface water being released and the environmental conditions of low temperature, low atmospheric pressure, and low acceleration due to gravity. We suggest how some of the assumptions made in previous work may be improved. Our model considers the thermodynamic effects of simultaneous evaporation and freezing of water, and fluid dynamical processes including changes in flow rheology caused by assimilation of cold rock and ice eroded at the channel bed, and ice crystal growth due to water freezing. We model how far initially turbulent water could flow in a channel before it erodes and entrains enough material to become laminar, and subsequently ceases to erode the bed. An ice raft will begin to form on the flood while transition occurs between turbulent and laminar flow. Estimates are given for water transit times, ～17–19 h, initial water depths, 50–62 m, and average flow speeds, 5–12 m s^?1, in the Mangala and Athabasca Valles. We show that these two outflow channels, and by implication others like them, could plausibly have been formed in single water release events. Resulting mean erosion rates are approximately 0.7 mm s^?1, a factor of three greater than previous estimates based on combinations of estimates of flood duration and required water volumes. This is explained by the consideration of the effects of eroded ice and the physics of thermal erosion in the present study.     

Assessment of the ecological quality status of soft-bottoms in Reunion Island (tropical Southwest Indian Ocean) using AZTI marine biotic indices

The ability of the two synthetic marine biotic indices, AMBI and M-AMBI, to account for changes in the ecological quality of coastal soft bottoms of Reunion Island according to disturbances was assessed from macrobenthic samples collected in five sectors between 1994 and 2004. Samples were collected under non-perturbed conditions and at two sites subjected to heavy organic enrichment. Both indices are based on a classification of macrofauna into ecological groups (EG), and their transfer to tropical waters required some adaptations. These indices proved efficient in detecting a degradation of habitat quality. Their use resulted in the classification of all sites sampled between 1996 and 1998 as "good" or "high". M-AMBI nevertheless tended to result in the attribution of a slightly worse ecological quality status than AMBI. Together with an update of the EG species list for the Indian Ocean area, our results support the extension of both indices for the assessment of tropical soft bottoms.     

Modeling of in-situ crystallization processes in the Permian mafic layered intrusion of Mont Collon (Dent Blanche nappe, western Alps)

The Mont Collon mafic complex is one of the best preserved examples of the Early Permian magmatism in the Central Alps, related to the intra-continental collapse of the Variscan belt. It mostly consists (> 95 vol.%) of ol + hy-normative plagioclase-wehrlites, olivine- and cpx-gabbros with cumulitic structures, crosscut by acid dikes. Pegmatitic gabbros, troctolites and anorthosites outcrop locally. A well-preserved cumulative sequence is exposed in the Dents de Bertol area (center of intrusion). P–T calculations indicate that this layered magma chamber emplaced at mid-crustal levels at about 0.5 GPa and 1100 °C. The Mont Collon cumulitic rocks record little magmatic differentiation, as illustrated by the restricted range of clinopyroxene mg-number (Mg#_cpx = 83–89). Whole-rock incompatible trace-element contents (e.g. Nb, Zr, Ba) vary largely and without correlation with major-element composition. These features are characteristic of an in-situ crystallization process with variable amounts of interstitial liquid L trapped between the cumulus mineral phases. LA-ICPMS measurements show that trace-element distribution in the latter is homogeneous, pointing to subsolidus re-equilibration between crystals and interstitial melts. A quantitative modeling based on Langmuir's in-situ crystallization equation successfully duplicated the REE concentrations in cumulitic minerals of all rock facies of the intrusion. The calculated amounts of interstitial liquid L vary between 0 and 35% for degrees of differentiation F of 0 to 20%, relative to the least evolved facies of the intrusion. L values are well correlated with the modal proportions of interstitial amphibole and whole-rock incompatible trace-element concentrations (e.g. Zr, Nb) of the tested samples. However, the in-situ crystallization model reaches its limitations with rock containing high modal content of REE-bearing minerals (i.e. zircon), such as pegmatitic gabbros. Dikes of anorthositic composition, locally crosscutting the layered lithologies, evidence that the Mont Collon rocks evolved in open system with mixing of intercumulus liquids of different origins and possibly contrasting compositions. The proposed model is not able to resolve these complex open systems, but migrating liquids could be partly responsible for the observed dispersion of points in some correlation diagrams. Absence of significant differentiation with recurrent lithologies in the cumulitic pile of Dents de Bertol points to an efficiently convective magma chamber, with possible periodic replenishment.     

Consequences of explosive eruptions on small Solar System bodies: the case of the missing basalts on the aubrite parent body

Aubrite meteorites are composed of constituents which are almost certainly of igneous origin. If they were generated by the melting and fractionation of enstatite chondrite-like parental material, as seems very likely, then plagioclase-rich, basaltic complements to the aubrites should have formed. However, such materials are not known as individual meteorites, and the compositions of two plagioclase-silica clasts and one albite-silica-(diopside-anorthite) clast (probably an impact melt) in the Norton County aubrite breccia suggest that they are not the putative enstatite-plagioclase basalts. We propose a new mechanism that explains the absence of these materials, showing that the expansion of even very small amounts of volatiles present in a melt approaching the surface of a small, low-gravity body will be enough to disrupt the melt into a spray of droplets moving faster than the local escape velocity. This explosive volcanic process of melt removal requires larger melt volatile contents on larger bodies, and data on the solubility of volatiles in basaltic melts suggest that the process was limited to bodies smaller than about 100 km in radius.