University of Oxford

The Tierra Blanca Joven (TBJ) eruption from Ilopango volcano deposited thick ash over much of El Salvador when it was inhabited by the Maya, and rendered all areas within at least 80 km of the volcano uninhabitable for years to decades after the eruption. Nonetheless, the more widespread environmental and climatic impacts of this large eruption are not well known because the eruption magnitude and date are not well constrained. In this multifaceted study we have resolved the date of the eruption to 431 ± 2 CE by identifying the ash layer in a well-dated, high-resolution Greenland ice-core record that is >7,000 km from Ilopango; and calculated that between 37 and 82 km3 of magma was dispersed from an eruption coignimbrite column that rose to ∼45 km by modeling the deposit thickness using state-of-the-art tephra dispersal methods. Sulfate records from an array of ice cores suggest stratospheric injection of 14 ± 2 Tg S associated with the TBJ eruption, exceeding those of the historic eruption of Pinatubo in 1991. Based on these estimates it is likely that the TBJ eruption produced a cooling of around 0.5 °C for a few years after the eruption. The modeled dispersal and higher sulfate concentrations recorded in Antarctic ice cores imply that the cooling would have been more pronounced in the Southern Hemisphere. The new date confirms the eruption occurred within the Early Classic phase when Maya expanded across Central America.

Measuring the low bromine abundances in Earth's materials remains an important challenge in order to constrain the geodynamical cycle of this element. Suitable standard materials are therefore required to establish reliable analytical methods to quantify Br abundances. In this study we characterise 21 Br-doped glasses synthesized from natural volcanic rocks of mafic to silicic compositions, in order to produce a new set of standards for Br analyses using various techniques. The nominal Br contents (amounts of Br loaded in the experimental samples) of 15 of 21 glasses were confirmed within 20% by instrumental neutron activation analysis (INAA). Using this newset of standards, we compare three micro-analytical approaches to measure Br contents in silicate glasses: synchrotron X-ray fluorescence (SR-XRF), laser ablation-inductively coupled plasma mass spectrometry (LA-ICP-MS), and secondary ion mass spectrometry (SIMS). With SR-XRF, the Br contents of the standard glasses were determined with the highest accuracy (b10% for Br ≥ 10 ppm; N25% for Br ≤ 5 ppm), and high precision (b10% for Br contents N10 ppm; 20–30% for Br ≤ 10 ppm). The detection limit was estimated to be b1 ppm Br. All those factors combined with a high spatial resolution (5 × 5 μm for the presented measurements), means that SR-XRF iswell suited to determine the lowBr abundance in natural volcanic glasses (crystal-hosted melt inclusions or matrix glasses of crystallized samples). At its current stage of development, the LA-ICP-MS method allows the measurement of hundreds to thousands ppm Br in silicate glasses with a precision and accuracy generally within 20%. The Br detection limit of this method has not been estimated but its low spatial resolution (90 μm) currently prevents its use to characterise natural volcanic glasses, however it is fully appropriate to analyse super liquidus or sparsely phyric, Br-rich experimental charges. Our study shows that SIMS appears to be a promising technique to measure the lowBr contents of natural volcanic glasses. Its spatial resolution is relatively good (~15 μm) and, similarly to SR-XRF, the detection limit is estimated to be ≤ 1 ppm. Using our new set of standards, the Br contents of two MPI-DING reference glasses containing ≤ 1.2 ppm of Br were reproduced with precision b5% and accuracy b20%. Moreover, SIMS presents the advantage of being a more accessible instrument than SR-XRF and data processing is more straightforward.

We investigate numerically the passage of spontaneous, dynamic in-plane shear ruptures from initiation to their final rupture speed, using very fine grids. By carrying out more than 120 simulations, we identify two different mechanisms controlling supershear transition. For relatively weaker faults, the rupture speed always passes smoothly and continuously through the range of speeds between the Rayleigh and shear wave speeds (the formerly considered forbidden zone of rupture speeds). This, however, occurs in a very short time, before the ruptures reach the compressional wave speed. The very short time spent in this range of speeds may explain why a jump over these speeds was seen in some earlier numerical and experimental studies and confirms that this speed range is an unstable range, as predicted analytically for steady state, singular cracks. On the other hand, for relatively stronger faults, we find that a daughter rupture is initiated by the main (mother) rupture, ahead of it. The mother rupture continues to propagate at sub-Rayleigh speed and eventually merges with the daughter rupture, whose speed jumps over the Rayleigh to shear wave speed range. We find that this daughter rupture is essentially a “pseudorupture,” in that the two sides of the fault are already separated, but the rupture has negligible slip and slip velocity. After the mother rupture merges with it, the slip, the slip velocity, and the rupture speed become dominated by those of the mother rupture. The results are independent of grid sizes and of methods used to nucleate the initial rupture.

Research on global ice-volume changes during Pleistocene glacial cycles is hindered by a lack of detailed sea-level records for time intervals older than the last interglacial. Here we present the first robustly dated, continuous and highly resolved records of Red Sea sea level and rates of sea-level change over the last 500,000 years, based on tight synchronization to an Asian monsoon record. We observe maximum ‘natural’ (pre-anthropogenic forcing) sea-level rise rates below 2m per century following periods with up to twice present-day ice volumes, and substantially higher rise rates for greater ice volumes. We also find that maximum sea-level rise rates were attained within 2 kyr of the onset of deglaciations, for 85% of such events. Finally, multivariate regressions of orbital parameters, sea-level and monsoon records suggest that major meltwater pulses account for millennial-scale variability and insolation-lagged responses in Asian monsoon records.

We use velocity measurements from a network of continuous GPS sites spanning the Apennines of peninsular Italy to test the hypothesis that the active deformation of the region is explained by variations in gravitational potential energy of the lithosphere. The simple geometry of the mountain chain allows us to treat the deformation as two-dimensional, neglecting gradients of velocity along the strike of the chain. Under this assumption, the integral of gravitational potential energy per unit area of the lithosphere (GPE) in the direction perpendicular to the chain is related by a simple expression to the velocity in the same direction. We show that the observed velocities match this expression with an RMS misfit of 0.5 mm/yr. This agreement suggests that deformation of the Apennines reflects a balance, within the mountain chain itself, between lateral variations in GPE and the stresses required to deform the lithosphere. Forces arising from processes external to the belt are not required to explain the observations.

The recent availability of small, cheap ultraviolet spectrometers has facilitated the rapid deployment of automated networks of scanning instruments at several volcanoes, measuring volcanic SO2 gas flux at high frequency. These networks open up a range of other applications, including tomographic reconstruction of the gas distribution which is of potential use for both risk mitigation, particularly to air traffic, and environmental impact modelling. Here we present a methodology for visualising reconstructed plumes using virtual globes, such as GoogleEarth, which allows animations of the evolution of the gas plume to be displayed and easily shared on a common platform. We detail the process used to convert tomographically reconstructed cross-sections into animated gas plume models, describe how this process is automated and present results from the scanning network around Mt.Etna, Sicily. We achieved an average rate of one frame every12 min, providing a good visual representation of the plume which can be examined from all angles. Increating these models, an approximation to turbulent diffusion in the atmosphere was required. To this end we derived the value of the turbulent diffusion coefficient for quiescent conditions near Etna to be around 200–500 m2s-1.