Fluid oscillations in a laboratory geyser with a bubble trap

Author Posting. © The Author(s), 2018. This is the author's version of the work. It is posted here under a nonexclusive, irrevocable, paid-up, worldwide license granted to WHOI. It is made available for personal use, not for redistribution. The definitive version was published in Rudolph, M. L., Sohn, R. A., & Lev, E.. Fluid oscillations in a laboratory geyser with a bubble trap. Journal of Volcanology and Geothermal Research, 368, (2018):100-110. doi:10.1016/j.jvolgeores.2018.11.003.Geysers are rare geologic features that episodically erupt water and steam. While it is understood that the eruptions are triggered by the conversion of thermal to kinetic energy during decompression of hot uids, geysers commonly exhibit a range of dynamic behaviors in-between and during eruptions that have yet to be adequately explained. In-situ measurements of temperature and pressure as well as remote geophysical techniques have revealed oscillatory behavior across a range of timescales, ranging from eruption cycles to impulsive bubble collapse events. Many geysers, including Old faithful in Yellowstone National Park, USA, are believed to have o set subsurface reservoirs (referred to as a `bubble trap') that can trap and accumulate noncondensable gas or steam entering the system. The impact of a bubble trap on the dynamic behaviors of the system, however, has not been fully established. We constructed a laboratory bubble trap and performed a series of experiments to study how uids oscillate back and forth between the eruption conduit and laterally-offseet reservoir in-between eruptions. We present a new theoretical model based on Hamiltonian mechanics that successfully predicts the oscillation frequencies observed in our experiments based on the conduit system geometry, the amount of gas that has accumulated in the bubble trap, and the amount of liquid water in the system. We demonstrate that when scaled to Old Faithful Geyser, this mechanism is capable of producing oscillations at the observed frequencies.The authors thank Paul Fucile and Glenn Macdonald for engineering support in designing and constructing the laboratory analog geyser rig. Funding for the laboratory geyser was provided by the US National Science Foundation grant EAR-1516361. EL was funded through a RISE award from Columbia University.2019-11-1

Periodic outgassing as a result of unsteady convection in Ray Lava Lake, Mount Erebus, Antarctica

Persistently active lava lakes show continuous outgassing and open convection over years to decades. Ray Lake, the lava lake at Mount Erebus, Ross Island, Antarctica, maintains long-term, near steady-state behavior in temperature, heat flux, gas flux, lake level, and composition. This activity is superposed by periodic small pulses of gas and hot magma every 5-18 minutes and disrupted by sporadic Strombolian eruptions. The periodic pulses have been attributed to a variety of potential processes including unstable bidirectional flow in the conduit feeding the lake. In contrast to hypotheses invoking a conduit source for the observed periodicity, we test the hypothesis that the behavior could be the result of dynamics within the lake itself, independent of periodic influx from the conduit. We perform numerical simulations of convection in Ray Lake driven by both constant and periodic inflow of gas-rich magma from the conduit to identify whether the two cases have different observational signatures at the surface. Our simulations show dripping diapirs or pulsing plumes leading to observable surface behavior with periodicities in the range of 5-20 minutes. We conclude that a convective speed faster than the inflow speed can result in periodic behavior without requiring periodicity in conduit dynamics. This finding suggests that the surface behavior of lava lakes might be less indicative of volcanic conduit processes in persistently outgassing volcanoes than previously thought, and that dynamics within the lava lake itself may modify or overprint patterns emerging from the conduit

Rayleigh–Taylor instabilities with anisotropic lithospheric viscosity

Rocks often develop fabric when subject to deformation, and this fabric causes anisotropy of physical properties such as viscosity and seismic velocities. We employ 2-D analytical solutions and numerical flow models to investigate the effect of anisotropic viscosity (AV) on the development of Rayleigh–Taylor instabilities, a process strongly connected to lithospheric instabilities. Our results demonstrate a dramatic effect of AV on the development of instabilities—their timing, location, and, most notably, their wavelength are strongly affected by the initial fabric. Specifically, we find a significant increase in the wavelength of instability in the presence of AV which favours horizontal shear. We also find that an interplay between regions with different initial fabric gives rise to striking irregularities in the downwellings. Our study shows that for investigations of lithospheric instabilities, and likely of other mantle processes, the approximation of isotropic viscosity may not be adequate, and that AV should be included

Seismic anisotropy in Eastern Tibet from shear wave splitting reveals changes in lithospheric deformation

Knowledge about seismic anisotropy can provide important insight into the deformation of the crust and upper mantle beneath tectonically active regions. Here we focus on the southeastern part of the Tibetan plateau, in Sichuan and Yunnan provinces, SW China. We measured shear wave splitting of core-refracted phases (SKS and SKKS) at a temporary array of 25 IRIS-PASSCAL stations. We calculated splitting parameters using a multi-channel and a single-record cross-correlation method. Multiple layers of anisotropy cannot be ruled out but are not required by the data. A Fresnel zone analysis suggests that the shallow mantle (between 60 and 160 km depth) is the most likely source of anisotropy. The polarization directions reveal a pronounced transition from primarily north–south in the north (Sichuan) to mostly east–west orientations in the south (Yunnan). In the southern part of the study region, that is, south of ~26°N, the fast polarization directions do not correlate well with known surface features and geodetic estimates of the crustal displacement fields. Whereas GPS campaigns provide evidence suggesting north–south crustal flow across the Red River Fault, the pattern of anisotropy argues against such flow in the upper mantle. These observations support models that allow differential movement of upper crust relative to lithospheric mantle. In the northern part of the study region the relationships are more ambiguous and coherent deformation of the crust and mantle lithosphere cannot be excluded. The interpretation of the shear wave splitting results is non-unique, but we suggest that the observed N–S transition reflects a fundamental change in deformation regime across our study region. It may be related to lateral variations in lithospheric rheology, or may mark a transition from the direct impact of the continental collision to dominance of the far-field strain field associated with regional subduction processes. Understanding the nature of the lateral change in deformation regime may prove critical for our understanding the geotectonic evolution of (eastern) Tibet

The effect of environmental conditions on the physiological response during a stand-up paddle surfing session

Stand Up Paddleboard (SUP) surfing entails riding breaking waves and maneuvering the board on the wave face in a similar manner to traditional surfing. Despite some scientific investigations on SUP, little is known about SUP surfing. The aim of this study was to investigate the physiological response during SUP surfing sessions and to determine how various environmental conditions can influence this response. Heart rate (HR) of an experienced male SUP surfer aged 43 was recorded for 14.9 h during ten surfing sessions and synced with on board video footage to enable the examination of the effect of different surfing modes and weather conditions on exercise intensity. Results indicated that the SUP surfer’s HR was above 70% of HRmax during 85% of each session, with the greatest heart rates found during falls off the board (~85% HRmax) and while paddling back to the peak (~83% HRmax). Total time surfing a wave was less than 5%, with the majority of time spent paddling back into position. Wind speed positively correlated with HR (r = 0.75, p < 0.05) and wave height negatively correlated with wave caching frequency (r = 0.73, p < 0.05). The results highlight the aerobic fitness for SUP surfing, where wave riding, paddling back to the peak, and falls appear to be associated with the greatest cardiovascular demand and demonstrate that environmental conditions can have an effect on the physiological response during SUP surfing sessions

Stable isochronal synchronization of mutually coupled chaotic lasers

The dynamics of two mutually coupled chaotic diode lasers are investigated experimentally and numerically. By adding self feedback to each laser, stable isochronal synchronization is established. This stability, which can be achieved for symmetric operation, is essential for constructing an optical public-channel cryptographic system. The experimental results on diode lasers are well described by rate equations of coupled single mode lasers

Insights on lava–ice/snow interactions from large-scale basaltic melt experiments

Quantitative measurements of interactions between lava and ice/snow are critical for improving our knowledge of glaciovolcanic hazards and our ability to use glaciovolcanic deposits for paleoclimate reconstructions. However, such measurements are rare because the eruptions tend to be dangerous and not easily accessible. To address these difficulties, we conducted a series of pilot experiments designed to allow close observation, measurements, and textural documentation of interactions between basaltic melt and ice. Here we report the results of the first experiments, which comprised controlled pours of as much as 300 kg of basaltic melt on top of ice. Our experiments provide new insights on (1) estimates for rates of heat transfer through boundary layers and for ice melting; (2) controls on rates of lava advance over ice/snow; (3) formation of lava bubbles (i.e., Limu o Pele) by steam from vaporization of underlying ice or water; and (4) the role of within-ice discontinuities to facilitate lava migration beneath and within ice. The results of our experiments confirm field observations about the rates at which lava can melt snow/ice, the efficacy with which a boundary layer can slow melting rates, and morphologies and textures indicative of direct lava-ice interaction. They also demonstrate that ingestion of external water by lava can create surface bubbles (i.e., Limu) and large gas cavities. We propose that boundary layer steam can slow heat transfer from lava to ice, and present evidence for rapid isotopic exchange between water vapor and melt. We also suggest new criteria for identifying ice-contact features in terrestrial and martian lava flows

Investigating lava flow rheology using video analysis and numerical flow models

Lava rheology is a major control on lava flow behavior and a critical parameter in flow simulations, but is very difficult to measure at field conditions or correctly extrapolate from the lab scale. We present a new methodology for investigating lava rheology through a combination of controlled experiments, image analysis and numerical forward modeling. Our experimental setup, part of the Syracuse University Lava Project (http://lavaproject.syr.edu) includes a large furnace capable of melting up to 450 kg of basalt, at temperatures well above the basalt liquidus. The lava is poured onto either a tilted bed of sand or a steel channel to produce meter-long flows. This experimental setup is probably the only facility that allows such large scale controlled lava flows made of natural basaltic material. We document the motion of the lava using a high-resolution video camera placed directly above the flows, and the temperature using infrared probes and cameras. After collecting the footage, we analyze the images for lava deformation and compare with numerical forward-models to constrain the rheological parameters and laws which best describe the flowing lava. For the video analysis, we employ the technique of Differential Optical Flow, which uses the time-variations of the spatial gradients of the image intensity to estimate velocity between consecutive frames. An important benefit for using optical flow, compared with other velocimetry methods, is that it outputs a spatially coherent flow field rather than point measurements. We demonstrate that the optical flow results agree with other measures of the flow velocity, and estimate the error due to noise and time-variability to be under 30% of the measured velocity. Our forward-models are calculated by solving the Stokes flow equations on an unstructured finite-element mesh defined using the geometry of the observed flow itself. We explore a range of rheological parameters, including the lava's apparent viscosity, the power-law exponent m and the thermal activation energy. Our measurements of apparent viscosity agree well with predictions of the composition-based Shaw (1972) and GRD model (Giordano, Russell and Dingwell, 2008). We find that for the high-temperature portion of the flow a weakly shear-thinning or Newtonian rheology (m > 0.7) with an effective activation energy of B = 5500 J gives the best fit to the data. Our methodology is the first time that high-resolution optical flow analysis of flowing lava is combined with numerical flow models to constrain rheology. The methodology we present here can be used in field conditions to obtain in-situ information on lava rheology, without physical interaction with the flow and without being limited to point-wise, low strain-rate, local measurements currently available through the use of rotational viscometers in the field

Benchmarking computational fluid dynamics models of lava flow simulation for hazard assessment, forecasting, and risk management

Abstract Numerical simulations of lava flow emplacement are valuable for assessing lava flow hazards, forecasting active flows, designing flow mitigation measures, interpreting past eruptions, and understanding the controls on lava flow behavior. Existing lava flow models vary in simplifying assumptions, physics, dimensionality, and the degree to which they have been validated against analytical solutions, experiments, and natural observations. In order to assess existing models and guide the development of new codes, we conduct a benchmarking study of computational fluid dynamics (CFD) models for lava flow emplacement, including VolcFlow, OpenFOAM, FLOW-3D, COMSOL, and MOLASSES. We model viscous, cooling, and solidifying flows over horizontal planes, sloping surfaces, and into topographic obstacles. We compare model results to physical observations made during well-controlled analogue and molten basalt experiments, and to analytical theory when available. Overall, the models accurately simulate viscous flow with some variability in flow thickness where flows intersect obstacles. OpenFOAM, COMSOL, and FLOW-3D can each reproduce experimental measurements of cooling viscous flows, and OpenFOAM and FLOW-3D simulations with temperature-dependent rheology match results from molten basalt experiments. We assess the goodness-of-fit of the simulation results and the computational cost. Our results guide the selection of numerical simulation codes for different applications, including inferring emplacement conditions of past lava flows, modeling the temporal evolution of ongoing flows during eruption, and probabilistic assessment of lava flow hazard prior to eruption. Finally, we outline potential experiments and desired key observational data from future flows that would extend existing benchmarking data sets