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.

Previous work has shown that Mo isotopes measurably fractionate between metal and silicate liquids, even at temperatures appropriate for core formation. However, the effect of variations in the structural environment of Mo in the silicate liquid, especially as a function of valence state, on Mo isotope fractionation remained poorly explored. We have investigated the role of valence state in metal-silicate experiments in a gas-controlled furnace at 1400 degrees C and at oxygen fugacities between 10(-12.7) and 10(-9.9), i.e. between three and 0.2 log units below the iron-wustite buffer. Two sets of experiments were performed, both with a silicate liquid in the CaO-Al2O3-SiO2 system. One set used molybdenum metal wire loops as the metal source, the other liquid gold alloyed with 2.5 wt % Mo contained in silica glass tubes. X-ray absorption near-edge spectroscopy analysis indicates that Mo6+/Sigma Mo in the silicate glasses varies between 0.24 and 0.77 at oxygen fugacities of 10-(12.0) and 10(-9.9) in the wire loop experiments and between 0.15 and 0.48 at 10(-11.4) and 10(-9.9) in the experiments with Au-Mo alloys. Double spiked analysis of Mo isotope compositions furthermore shows that Mo isotope fractionation between metal and silicate is a linear function of Mo6+/Sigma Mo in the silicate glasses, with a difference of 0.51 parts per thousand in Mo-98/Mo-95 between purely Mo4+-bearing and purely Mo6+-bearing silicate liquid. The former is octahedrally and the latter tetrahedrally coordinated. Our study implies that previous experimental work contained a mixture of Mo4+ and Mo6+ species in the silicate liquid. Our refined parameterisation for Mo isotope fractionation between metal and silicate can be described as

The Alaska Volcano Observatory (AVO) monitors volcanoes in Alaska and issues notifications and warnings of volcanic unrest and eruption. We evaluate the timeliness and accuracy of eruption forecasts for 53 eruptions at 20 volcanoes, beginning with Mount Redoubt's 1989-1990 eruption. Successful forecasts are defined as those where AVO issued a formal warning before eruption onset. These warning notifications are now part of AVO's Aviation Color Code and Volcanic Alert Level. This analysis considers only the start of an eruption, although many eruptions have multiple phases of activity. For the 21 eruptions at volcanoes with functioning local seismic networks, AVO has high forecasting success at volcanoes with: >15 years repose intervals and magmatic eruptions (4 out of 4, 100%); or larger eruptions (Volcanic Explosivity Index (VEI) 3 or greater; 6 out of 10, 60%). Therefore, AVO successfully forecast all four monitored, longer-repose period, VEI 3+ eruptions: Redoubt 1989-1990 and 2009, Spurr 1992, and Augustine 2005-2006. For volcanoes with functioning seismic monitoring networks, success rates are lower for: volcanoes with shorter repose periods (3 out of 16, 19%); more mafic compositions (3 out of 18, 17%); or smaller eruption size (VEI 2 or less, 1 out of 11, 9%). These eruptions (Okmok, Pavlof, Veniaminof, and Shishaldin) often lack detectable precursory signals. For 32 eruptions at volcanoes without functioning local seismic networks, the forecasting success rate is much lower (2, 6%; Kasatochi 2008 and Shishaldin 2014). For remote volcanoes where the main hazard is to aviation, rapid detection is a goal in the absence of in situ monitoring. Eruption detection has improved in recent years, shown by a decrease in the time between eruption onset and notification. Even limited seismic monitoring can detect precursory activity at volcanoes with certain characteristics (intermediate composition, longer repose times, larger eruptions), but difficulty persists in detecting subtle precursory activity at frequently active volcanoes with more mafic compositions. This suggests that volcano-specific characteristics should be considered when designing monitoring programs and evaluating forecasting success. More proximally-located sensors and data types are likely needed to forecast eruptive activity at frequently-active, moremafic volcanoes that generally produce smaller eruptions.