Volcanic eruptions are generally forecast based on strong increases in monitoring parameters such as seismicity or gas emissions above a relatively low background level (e.g., Voight, 1988; Sparks, 2003). Because of this, forecasting individual explosions during an ongoing eruption, or at persistently restless volcanoes, is difficult as seismicity, gas emissions, and other indicators of unrest are already in a heightened state. Therefore, identification of short-term precursors to individual explosions at volcanoes already in heightened states of unrest, and an understanding of explosion trigger mechanisms, is important for the reduction of volcanic risk worldwide. Seismic and visual observations at Telica Volcano, Nicaragua, demonstrate that a) episodes of seismic quiescence reliably preceded explosions during an eruption in May 2011 and b) the duration of precursory quiescence and the energy released in the ensuing explosion were strongly correlated. Precursory seismic quiescence is interpreted as the result of sealing of shallow gas pathways, leading to pressure accumulation and eventual catastrophic failure of the system, culminating in an explosion. Longer periods of sealing and pressurization lead to greater energy release in the ensuing explosion. Near-real-time observations of seismic quiescence at restless or erupting volcanoes can thus be useful for both timely eruption warnings and for forecasting the energy of impending explosions. (C) 2016 Elsevier B.V. All rights reserved.

Pressure has an essential role in the production(1) and control(2,3) of superconductivity in iron-based superconductors. Substitution of a large cation by a smaller rare-earth ion to simulate the pressure effect has raised the superconducting transition temperature T-c to a record high of 55 K in these materials(4,5). In the same way as T-c exhibits a bell-shaped curve of dependence on chemical doping, pressure-tuned T-c typically drops monotonically after passing the optimal pressure(1-3). Here we report that in the superconducting iron chalcogenides, a second superconducting phase suddenly re-emerges above 11.5 GPa, after the T-c drops from the first maximum of 32 K at 1 GPa. The T-c of the re-emerging superconducting phase is considerably higher than the first maximum, reaching 48.0-48.7 K for Tl0.6Rb0.4Fe1.67Se2, K0.8Fe1.7Se2 and K0.8Fe1.78Se2.

A joint analysis of geodetic and seismic datasets from Kilauea Volcano during a period of magmatic unrest in 2006 demonstrates the effectiveness of this combination for testing and constraining models of magma dynamics for a complex, multi-source system. At the end of 2003, Kilauea's summit began a four-year-long period of inflation due to a surge in magma supply to the volcano. In 2006, for the first time since 1982, Kilauea's Southwest Rift Zone (SWRZ) also experienced inflation. To investigate the characteristics of active magma sources and the nature of their interactions with faults in the SWRZ during 2006, we integrate, through Coulomb stress modeling, contemporary geodetic data from InSAR and GPS with a new catalogue of double-couple fault-plane solutions for volcano-tectonic earthquakes. We define two periods of inflation during 2006 based on the rate of deformation measured in daily GPS data, spanning February to 15 March 2006 (Period 1) and 16 March to 30 September 2006 (Period 2). InSAR data for these two periods are inverted to determine the position, change in size, and shape of inflation sources in each period. Our new models are consistent with microseismic activity from each period. They suggest that, during Period 1, deformation in the SWRZ can be explained by pressurization of magma in a spherical reservoir beneath the south caldera, and that, during Period 2, magma was also aseismically intruded farther to the southwest into the SWRZ along a sub-horizontal plane. Our Coulomb stress analysis shows that the microseismicity recorded in the SWRZ is induced by overpressurization of the south caldera reservoir, and not by magma intrusion into the SWRZ. This study highlights the importance of a joint analysis of independent geophysical datasets to fully constrain the nature of magma accumulation. (C) 2016 Elsevier B.V. All rights reserved.

A primary consequence of plate tectonics is that basaltic oceanic crust subducts with lithospheric slabs into the mantle. Seismological studies extend this process to the lower mantle, and geochemical observations indicate return of oceanic crust to the upper mantle in plumes. There has been no direct petrologic evidence, however, of the return of subducted oceanic crustal components from the lower mantle. We analyzed superdeep diamonds from Juina-5 kimberlite, Brazil, which host inclusions with compositions comprising the entire phase assemblage expected to crystallize from basalt under lower-mantle conditions. The inclusion mineralogies require exhumation from the lower to upper mantle. Because the diamond hosts have carbon isotope signatures consistent with surface-derived carbon, we conclude that the deep carbon cycle extends into the lower mantle.

We report a finding of a pressure-induced quantum critical transition in K0.8FexSe2 (x = 1.7 and 1.78) superconductors through in situ high-pressure electrical transport and x-ray diffraction measurements in diamond anvil cells. Transitions from metallic Fermi liquid behavior to non-Fermi liquid behavior and from antiferromagnetism to paramagnetism are found in the pressure range of 9.2-10.3 GPa, in which superconductivity tends to disappear. The change around the quantum critical point from the coexisting antiferromagnetism state and the Fermi liquid behavior to the paramagnetism state and the non-Fermi liquid behavior in the iron-selenide superconductors demonstrates a unique mechanism for their quantum critical transition.

Harmonic tremor is a common feature of volcanic, hydrothermal, and ice sheet seismicity and is thus an important proxy for monitoring changes in these systems. However, no automated methods for detecting harmonic tremor currently exist. Because harmonic tremor shares characteristics with speech and music, digital signal processing techniques for analyzing these signals can be adapted. I develop a novel pitch-detection-based algorithm to automatically identify occurrences of harmonic tremor and characterize their frequency content. The algorithm is applied to seismic data from Popocatepetl Volcano, Mexico, and benchmarked against a monthlong manually detected catalog of harmonic tremor events. During a period of heightened eruptive activity from December 2014 to May 2015, the algorithm detects 1465 min of harmonic tremor, which generally precede periods of heightened explosive activity. These results demonstrate the algorithm's ability to accurately characterize harmonic tremor while highlighting the need for additional work to understand its causes and implications at restless volcanoes.

We present a method to optimize absorption cells for precise wavelength calibration in the near-infrared. We apply it to design and optimize methane isotopologue cells for precision radial velocity measurements in the K band. We also describe the construction and installation of two such cells for the CSHELL spectrograph at NASA's IRTF. We have obtained their high-resolution laboratory spectra, which we can then use in precision radial velocity measurements and which can also have other applications. In terms of obtainable RV precision, methane should outperform other proposed cells, such as the ammonia cell ((NH3)-N-14) recently demonstrated on CRIRES/VLT. The laboratory spectra of the ammonia and methane cells show strong absorption features in the H band that could also be exploited for precision Doppler measurements. We present spectra and preliminary radial velocity measurements obtained during our first-light run. These initial results show that a precision down to 20-30 m s(-1) can be obtained using a wavelength interval of only 5 nm in the K band and S/N similar to 150. This supports the prediction that a precision down to a few meters per second can be achieved on late-M dwarfs using the new generation of NIR spectrographs, thus enabling the detection of terrestrial planets in their habitable zones. Doppler measurements in the NIR can also be used to mitigate the radial velocity jitter due to stellar activity, enabling more efficient surveys on young active stars.

During the run-up to eruption, volcanoes often show geophysically detectable signs of unrest. However, there are long-standing challenges in interpreting the signals and evaluating the likelihood of eruption, especially during the early stages of volcanic unrest. Considerable insight can be gained from combined geochemical and geophysical studies. Here we take such an approach to better understand the beginning of eruption run-up, viewed through the lens of the 1999 sub-Plinian basaltic eruption of Shishaldin volcano, Alaska. The eruption is of interest due to its lack of observed deformation and its apparent long run-up time (9 months), following a deep long-period earthquake swarm. We evaluate the nature and timing of recharge by examining the composition of 138 olivine macrocrysts and 53 olivine-hosted melt inclusions and through shear-wave splitting analysis of regional earthquakes. Magma mixing is recorded in three crystal populations: a dominant population of evolved olivines (Fo(60-69)) that are mostly reversely zoned, an intermediate population (Fo(69-76)) with mixed zonation, and a small population of normally zoned more primitive olivines (Fo(76-80)). Mixing-to-eruption timescales are obtained through modeling of Fe-Mg interdiffusion in 78 olivines. The large number of resultant timescales provides a thorough record of mixing, demonstrating at least three mixing events: a minor event similar to 11 months prior to eruption, overlapping within uncertainty with the onset of deep long-period seismicity; a major event similar to 50 days before eruption, coincident with a large (M5.2) shallow earthquake; and a final event about a week prior to eruption. Shear-wave splitting analysis shows a change in the orientation of the local stress field about a month after the deep long-period swarm and around the time of the M5.2 event. Earthquake depths and vapor saturation pressures of Raman-reconstructed melt inclusions indicate that the recharge magma originated from depths of at least 20 km, and that mixing with a shallow magma or olivine cumulates occurred in or just below the edifice (<3 km depth). Deformation was likely outside the spatial and temporal resolution of the satellite measurements. Prior to eruption magma was stored over a large range of depths (similar to 0-2.5 km below the summit), suggesting a shallow, vertical reservoir that could provide another explanation for the lack of detectable deformation. The earliest sign of unrest (deep long-period seismicity) coincides temporally with magmatic activity (magma mixing and a change in the local stress state), possibly indicating the beginning of eruption run-up. The more immediate run-up began with the major recharge event similar to 50 days prior to eruption, after which the signs of unrest became continuous. This timescale is long compared to the seismic run-up to other basaltic eruptions (typically hours to days). Other volcanoes classified as open-system, based on their lack of precursory deformation, also tend to have relatively long run-up durations, which may be related to the time required to fill the shallow reservoir with magmas sourced from greater depth. (C) 2018 Elsevier B.V. All rights reserved.

High pressure x-ray diffraction measurements on perovskite RMnO3 (R = Dy, Ho and Lu) reveal that varying structural changes occur for different R ions. Large lattice changes (orthorhombic strain) occur in DyMnO3 and HoMnO3 while the Jahn-Teller (JT) distortion remains stable. Conversely, in the small R-ion system LuMnO3, Mn-O bond distortions are observed between 4 and 8 GPa with a broad minimum in the JT distortion. High pressure infrared measurements indicate that a phonon near 390 cm(-1) corresponding to the complex motion of the Mn and O ions changes anomalously for LuMnO3. It softens in the 4-8 GPa region, which is consistent with the structural change in Mn-O bonds and then hardens at higher pressures. By contrast, the phonons continuously harden with increasing pressure for DyMnO3 and HoMnO3. Density functional theory methods show that E-phase LuMnO3 is the most stable phase up to the 10 GPa pressure examined. Simulations indicate that the distinct structural change under pressure in LuMnO3 can possibly be used to optimize the electric polarization by pressure/strain.