Steep schlieren and associated enclaves in the Vinalhaven granite,Maine: possible indicators for granite rheology

Schlieren that form above the solidus are potentially diagnostic of the rheology of crystallizing granitic magmas and could provide insights into magma chamber structure. Here we focus on steep schlieren associated with comagmatic enclaves in the Vinalhaven granite. Steep schlieren extend 1–3 m upward from the sides of enclaves and appear to have steep cylindrical shapes in 3D. Analyses of schlieren widths and the sizes of associated enclaves suggest that granitic crystal mush in which they occur had a plastic rheology probably characterized by a yield strength. The enclaves, now enclosed in coarse-grained granite, must have existed at higher levels in a crystal-poor part of a magma chamber and settled downward until reaching material with a yield strength exceeding the buoyancy of the enclaves. In addition to constraining the local rheology of the granite, their relative positions may indicate vertical rheological variations and possibly the chamber floor.     

Quantitative field constraints on the dynamics of silicic magma chamber rejuvenation and overturn

The catastrophic eruption of large-volume, crystal-rich silicic magmas is often proposed to be a consequence of reheating, melting and overturn of partially molten, buoyant silicic material following repeated injection of dense, hot mafic magma. To test this “rejuvenation hypothesis”, we analyze at high spatial resolution 33 examples of deformed interfaces between intrusive mafic and silicic layers in two plutons of the Coastal Maine Magmatic Province, USA. These deformed interfaces are thought to record the buoyant overturn of silicic crystal mush layers, apparently in response to the injection and cooling of hot, dense mafic magmas. We use spectral analysis and scaling theory along with petrologic and textural data to identify, characterize, and understand periodic deformations from the scale of individual crystals (≈1 cm) to the thicknesses of mafic and silicic layers. Deformations at the largest scale lengths (>100 m) are at wavelengths comparable to, or greater than, silicic layer thicknesses and support a conjecture that mafic recharge can cause large-scale Rayleigh–Taylor-type overturning of silicic mushy layers. By contrast, the smallest scales of individual crystals probably record effects related to production and buoyancy-driven rise of melt from the tops of silicic mushes in contact with overlying hot basalt, whereas intermediate scales are explained by compaction. Our results constrain the evolution of a thermal rejuvenation event and potentially identify a condition for a large-scale overturn of the magma chamber that may lead to eruption. This work provides the first quantitative field-based constraints on some of the key physical processes inherent to the rejuvenation hypothesis.     

The volcanic history of Volcán Alcedo,Galápagos Archipelago: A case study of rhyolitic oceanic volcanism

Volcán Alcedo is one of the seven western Galápagos shields and is the only active Galápagos volcano known to have erupted rhyolite as well as basalt. The volcano stands 4 km above the sea floor and has a subaerial volume of 200 km³, nearly all of which is basalt. As Volcán Alcedo grew, it built an elongate domal shield, which was partly truncated during repeated caldera-collapse and partial-filling episodes. An outward-dipping sequence of basalt flows at least 250 m thick forms the steepest (to 33°) flanks of the volcano and is not tilted; thus a constructional origin for the steep upper flanks is favored. About 1 km³ of rhyolite erupted late in the volcano's history from at least three vents and in 2–5 episodes. The most explosive of these produced a tephra blanket that covers the eastern half of the volcano. Homogeneous rhyolitic pumice is overlain by dacite-rhyolite commingled pumice, with no stratigraphic break. The tephra is notable for its low density and coarse grain size. The calculated height of the eruption plume is 23–30 km, and the intensity is estimated to have been 1.2x10⁸ kg/s. Rhyolitic lavas vented from the floor of the caldera and from fissures along the rim overlie the tephra of the plinian phase. The age of the rhyolitic eruptions is 120 ka, on the basis of K-Ar ages. Between ten and 20 basaltic lava flows are younger than the rhyolites. Recent faulting resulted in a moat around part of the caldera floor. Alcedo most resently erupted sometime between 1946 and 1960 from its southern flank. Alcedo maintains an active, transient hydrothermal system. Acoustic and seismic activity in 1991 is attributed to the disruption of the hydrothermal system by a regional-scale earthquake.     

Influence of supercontinents on deep mantle flow

The assembly of supercontinents should impact mantle flow fields significantly, affecting the distribution of subduction, upwelling plumes, lower mantle chemical heterogeneities, and thus plausibly contributing to voluminous volcanism that is often associated with their breakup. Alternative explanations for this volcanism include insulation by the continent and thus elevated subcontinental mantle temperatures. Here we model the thermal and dynamic impact of supercontinents on Earth-like mobile-lid convecting systems. We confirm that insulating supercontinents (over 3000 km extent) can impact mantle temperatures, but show the scale of temperature anomaly is significantly less for systems with strongly temperature-dependent viscosities and mobile continents. Additionally, for continents over 8000 km, mantle temperatures are modulated by the development of small-scale convecting systems under the continent, which arise due to inefficient lateral convection of heat at these scales. We demonstrate a statistically robust association between rising plumes supercontinental interiors for a variety of continental configurations, driven largely by the tendency of subducting slabs to lock onto continental margins. The distribution of slabs also affects the spatial positioning of deep mantle thermochemical anomalies, which demonstrate stable configurations in either the sub-supercontinent or intraoceanic domains. We find externally forced rifting scenarios unable to generate significant melt rates, and thus the ultimate cause of supercontinent breakup related volcanism is probably related to dynamic continental rifting in response to mantle reconfiguration events.     

Magma dynamics, crystallization, and chemical differentiation of the 1959 Kilauea Iki lava lake, Hawaii, revisited

Using constraints from an extensive database of geological and geochemical observations along with results from fluid mechanical studies of convection in magma chambers, we identify the main physical processes at work during the solidification of the 1959 Kilauea Iki lava lakes. In turn, we investigate their quantitative influence on the crystallization and chemical differentiation of the magma, and on the development of the internal structure of the lava lake. In contrast to previous studies, vigorous stirring in the magma, driven predominately by the descent of dense crystal-laden thermal plumes from the roof solidification front and the ascent of buoyant compositional plumes due to the in situ growth of olivine crystals at the floor, is predicted to have been an inevitable consequence of very strong cooling at the roof and floor. The flow is expected to have caused extensive but imperfect mixing over most of the cooling history of the magma, producing minor compositional stratification at the roof and thermal stratification at the floor. The efficient stirring of the large roof cooling is expected to have resulted in significant internal nucleation of olivine crystals, which ultimately settled to the floor. Additional forcing due to either crystal sedimentation or the ascent of gas bubbles is not expected to have increased significantly the amount of mixing. In addition to convection in the magma, circulation driven by the convection of buoyant interstitial melt in highly permeable crystal-melt mushes forming the roof and the floor of the lava lake is envisaged to have produced a net upward flow of evolved magma from the floor during solidification. In the floor zone, mush convection may have caused the formation of axisymmetric chimneys through which evolved magma drained from deep within the floor into the overlying magma and potentially the roof. We hypothesize that the highly evolved, pipe-like ‘vertical olivine-rich bodies’ (VORBs) [Bull. Volcanol. 43 (1980) 675] observed in the floor zone, of the lake are fossil chimneys. In the roof zone, buoyant residual liquid both produced at the roof solidification front and gained from the floor as a result of incomplete convective mixing is envisaged to have percolated or ‘leaked‘ into the overlying highly-permeable cumulate, displacing less buoyant interstitial melt downward. The results from Rayleigh fractionation-type models formulated using boundary conditions based on a quantitative understanding of the convection in the magma indicate that most of the incompatible element variation over the height of the lake can be explained as a consequence of a combination of crystal settling and the extensive but imperfect convective mixing of buoyant residual liquid released from the floor solidification front. The remaining chemical variation is understood in terms of the additional influences of mush convection in the roof and floor on the vertical distribution of incompatible elements. Although cooling was concentrated at the roof of the lake, the floor zone is found to be thicker than the roof zone, implying that it grew more quickly. The large growth rate of the floor is explained as a consequence of a combination of the substantial sedimentation of olivine crystals and more rapid in situ crystallization due to both a higher liquidus temperature and enhanced cooling resulting from imperfect thermal and chemical mixing.     

Osmo-power. Theory and performance of an osmo-power pilot plant

A theoretical and experimental study of the production of useful energy by the natural process of osmosis is presented. Using the results of the study a conceptual design of an osomotic pilot plant is performed. The power produced by a 1.6 MW km⁻² plant has a competitive cost with that produced by both fossil power plants and nuclear power plants.     

What controls the disequilibrium state of gravel-bed rivers?

A river at equilibrium is described by a statistically-stationary mean bed elevation profile that arises in response to steady supplies of relief, water and sediment. Outside of the profile shape, how is the equilibrium state of a river most reliably identified and rigorously defined? Motivated by a proposed link between equilibrium and physical processes, we use scaling theory to develop the dimensionless channel response number ξ=KU_b/U_p. ξ is a metric for the local disequilibrium state of gravel-bed mountain streams, which reflects a balance between the rate of topographic adjustment U_b, and the rate of bed sediment texture adjustment U_p. The coefficient K can take one of two forms depending on choice of length scale for topographic adjustment. We hypothesize that equilibrium occurs where and when ξ≈O(1), and consequently, disequilibrium is the more general state captured by conditions of ξ≉O(1). The rates U_b and U_p are controlled by the mechanics of sediment deposition and entrainment at the local scale of the channel width. The extent to which either process regulates disequilibrium depends on the bed strength, which is set by the time-varying grain size distribution and packing. We use flume experiments to understand ξ and find that in the limit ξ>>1, the time-varying response of an experimental channel depends sensitively on the spatially-averaged bed shear stress ratio τ/τ_ref. When τ/τ_ref≈1.5, U_b was the dominant control on disequilibrium. However, when τ/τ_ref≈2.0, U_p contributed more significantly to disequilibrium. These results suggest that after an upstream supply perturbation, the equilibrium timescale is governed by U_p, which we show is consistent with expectations from linear damping theory. Our experimental test of ξ is promising, but inconclusive with respect to our hypothesis. This uncertainty can be readily addressed with numerical or additional physical experiments. © 2019 John Wiley & Sons, Ltd.