Stratigraphy, chemistry, and eruptive dynamics of the 12.4 ka plinian eruption of Apoyeque volcano, Managua, Nicaragua

Apoyeque volcano, located 9 km northwest of Managua city, erupted explosively at 12.4 ka. The Plinian eruption deposited a widespread pumice fall deposit known as the Upper Apoyeque Tephra (UAq). The UAq is massive, reversely graded, and consists of white juvenile pumice (~78 vol.%), a variety of cognate lithics and accidental altered lithics. The whole-rock pumice composition is rhyodacitic (SiO₂?=?66.9–68.5 wt.%) with a mineral paragenesis of plagioclase, orthopyroxene, clinopyroxene, amphibole, titanomagnetite, and ilmenite in a rhyolitic glass groundmass (SiO₂?=?74.4?±?0.6 wt.%). The deposit’s dispersal axis is to the south, with the deposit covering a minimum area of 877 km² within the 50 cm isopach and has a total volume of 3 km³ (dense rock equivalent, 1.15 km³). The eruption column reached a maximum height of ca.28 km. The eruption ejected a total mass of 3?×?10¹² kg at an average rate of 2?×?10⁸ kg/s, and based on available models, we infer duration of almost 4 h. Petrographic and geochemical characteristics suggest that the eruption was triggered by magma mixing.     

Precursory activity of the 161 ka Kos Plateau Tuff eruption, Aegean Sea (Greece)

The Kos Plateau Tuff (KPT) eruption of 161 ka was the largest explosive Quaternary eruption in the eastern Mediterranean. We have discovered an uplifted beach deposit of abraded pumice cobbles, directly overlain by the KPT. The pumice cobbles resemble pumice from the KPT in petrography and composition and differ from Plio-Pleistocene rhyolites on the nearby Kefalos Peninsula. The pumice contains enclaves of basaltic andesite showing chilled lobate margins, suggesting co-existence of two magmas. The deposit provides evidence that the precursory phase of the KPT eruption produced pumice rafts, and defines the paleoshoreline for the KPT, which elsewhere was deposited on land. The beach deposit has been uplifted about 120 m since the KPT eruption, whereas the present marine area south of Kos has subsided several hundred metres, as a result of regional neotectonics. The basaltic andesite is more primitive than other mafic rocks known from the Kos–Nisyros volcanic centre and contains phenocrysts of Fo₈₉ olivine, bytownite, enstatite and diopside. Groundmass amphibole suggests availability of water in the final stages of magma evolution. Geochemical and mineralogical variation in the mafic products of the KPT eruption indicate that fractionation of basaltic magma in a base-of-crust magma chamber was followed by mixing with rhyolitic magma during eruption. Low eruption rates during the precursory activity may have minimised the extent of mixing and preserved the end-member magma types.     

Deep-bedded ultramafic diatremes in the Missouri River Breaks volcanic field,Montana, USA: 1 km of syn-eruptive subsidence

The ultramafic Eocene Missouri River Breaks volcanic field (MRBVF, Montana, USA) includes over 50 diatremes emplaced in a mostly soft substrate. The current erosion level is 1.3–1.5 km below the pre-eruptive surface, exposing the deep part of the diatreme structures and some dikes. Five representative diatremes are described here; they are 200-375 m across and have sub-vertical walls. Their infill consists mostly of 55-90 % bedded pyroclastic rocks (fine tuffs to coarse lapilli tuffs) with concave-upward bedding, and 45–10 % non-bedded pyroclastic rocks (medium lapilli tuffs to tuff breccias). The latter zones form steep columns 15–135 m in horizontal dimension, which cross-cut the bedded pyroclastic rocks. Megablocks of the host sedimentary formations are also present in the diatremes, some being found 1 km or more below their sources. The diatreme infill contains abundant lithic clasts and ash-sized particles, indicating efficient fragmentation of magma and country rocks. The spherical to sub-spherical juvenile clasts are non-vesicular. They are accompanied by minor accretionary lapilli and armored lapilli. The deposits of dilute pyroclastic density currents are locally observed. Our main interpretations are as follows: (1) the observations strongly support phreatomagmatic explosions as the energy source for fragmentation and diatreme excavation; (2) the bedded pyroclastic rocks were deposited on the crater floor, and subsided by 1.0–1.3 km to their current location, with subsidence taking place mostly during the eruption; (3) the observed non-bedded pyroclastic columns were created by debris jets that punched through the bedded pyroclastic material; the debris jets did not empty the mature diatreme, occupying only a fraction of its width, and some debris jets probably did not reach the crater floor; (4) the mature diatreme was nearly always filled and buttressed by pyroclastic debris at depth – there was never a 1.3–1.5-km-deep empty hole with sub-vertical walls, otherwise the soft substrate would have collapsed inward, which it only did near the surface, to create the megablocks. We infer that syn-eruptive subsidence shifted down bedded pyroclastic material and shallow sedimentary megablocks by 0.8–1.1 km or more, after which limited post-eruptive subsidence occurred. This makes the MRBVF diatremes an extreme end-member case of syn-eruptive subsidence in the spectrum of possibilities for maar-diatreme volcanoes worldwide.     

Reconstructing the eruption magnitude and energy budgets for the pre-historic eruption of the monogenetic ∼5 ka Mt. Gambier Volcanic Complex, south-eastern Australia

Understanding explosive volcanic eruptions, especially phreatomagmatic eruptions, their intensities and energy budgets is of major importance when it comes to risk and hazard studies. With only a few historic occurrences of phreatomagmatic activity, a large amount of our understanding comes from the study of pre-historic volcanic centres, which causes issues when it comes to preservation and vegetation. In this research, we show that using 3D geometrical modelling it is possible to obtain volume estimates for different deposits of a pre-historic, complex, monogenetic centre, the Mt. Gambier Volcanic Complex, south-eastern Australia. Using these volumes, we further explore the energy budgets and the magnitude of this eruption (VEI 4), including dispersal patterns (eruption columns varying between 5 and 10 km, dispersed towards north-east to south), to further our understanding of intraplate, monogenetic eruptions involving phreatomagmatic activity. We also compare which thermodynamic model fits best in the creation of the maar crater of Mt. Gambier: the major-explosion-dominated model or the incremental growth model. In this case, the formation of most of the craters can best be explained by the latter model.     

Volcanology and petrology of Mathews Tuya,northern British Columbia,Canada: glaciovolcanic constraints on interpretations of the 0.730 Ma Cordilleran paleoclimate

Petrological, volcanological and geochronological data collected at Mathews Tuya together provide constraints on paleoclimate conditions during formation of the edifice. The basaltic tuya was produced via Pleistocene glaciovolcanism in northern British Columbia, Canada, and is located within the Tuya volcanic field (59.195°N/130.434°W), which is part of the northern Cordilleran volcanic province (NCVP). The edifice comprises a variety of lithofacies, including columnar-jointed lava, pillow lava, massive dikes, and volcaniclastic rocks. Collectively these deposits record the transition from an explosive subaqueous to an effusive subaerial eruption environment dominated by Pleistocene ice. As is typical for tuyas, the volcaniclastic facies record multiple fragmentation processes including explosive, quench and mechanical fragmentation. All samples from Mathews Tuya are olivine-plagioclase porphyritic alkali olivine basalts. They are mineralogically and geochemically similar to nearby glaciovolcanic centers from the southeastern part of the Tuya volcanic field (e.g., Ash Mountain, South Tuya, Tuya Butte) as well as the dominant NCVP rock type. Crystallization scenarios calculated with MELTS account for variations between whole rock and glass compositions via low pressure fractionation. The presence of olivine microphenocrysts and the absence of pyroxene phenocrysts constrain initial crystallization pressures to less than 0.6 GPa. The eruption of Mathews Tuya occurred between 0.718 ± 0.054 Ma and 0.742 ± 0.081 Ma based on ⁴⁰Ar/³⁹Ar geochronology (weighted mean age of 0.730 Ma). The age determinations provide the first firm documentation for large (>700 m thick), pre-Fraser/Wisconsin glaciers in north-central British Columbia ~0.730 Ma, and correlate in age with glaciovolcanic deposits in Russia (e.g., Komatsu et al. Geomorph 88: 352-366, 2007) and with marine isotopic evidence for large global ice volumes ~0.730 Ma.     

Improved estimates of mean sea level changes in the German Bight over the last 166 years

In this paper, mean sea level changes in the German Bight, the south-eastern part of the North Sea, are analysed. Records from 13 tide gauges covering the entire German North Sea coastline and the period from 1843 to 2008 have been used to derive high quality relative mean sea level time series. Changes in mean sea level are assessed using non-linear smoothing techniques and linear trend estimations for different time spans. Time series from individual tide gauges are analysed and then ‘virtual station’ time series are constructed (by combining the individual records) which are representative of the German Bight and the southern and eastern regions of the Bight. An accelerated sea level rise is detected for a period at the end of the nineteenth century and for another one covering the last decades. The results show that there are regional differences in sea level changes along the coastline. Higher rates of relative sea level rise are detected for the eastern part of the German Bight in comparison to the southern part. This is most likely due to different rates of vertical land movement. In addition, different temporal behaviour of sea level change is found in the German Bight compared to wider regional and global changes, highlighting the urgent need to derive reliable regional sea level projections for coastal planning strategies.     

Study of the earthquake of the January 23, 1880, in San Cristóbal,Cuba and the Guane fault

All available data on the January 23, 1880, earthquake near San Cristobal, Western Cuba, are compiled and presented here. The earthquake reached a maximum intensity of eight degrees (MSK) and caused three fatalities. It was accompanied by 65 aftershocks and was felt as far away as the Florida Keys. Twentieth century specialists has associated this event, in its day the strongest recorded (Ms = 6.2) in the region, with the Pinar fault. The Pinar fault is well expressed topographically as the boundary between the Guaniguanico Range in the north and an alluvial plain to the south. Most of the major damage caused by the earthquake was located on the alluvial plain, which in consequence has been considered the epicenter area. In the study presented here, the data compiled from the first reports of Father Benito Vines Martorell, S.J., and Pedro Salteraín y Legarra, indicate that the seismic structure was located in the alluvial plain, and that it was the Guane fault, and not the Pinar fault, that was responsible for the earthquake. The Guane fault, found below the alluvial sediments, extends NE-SW for over 110 km. Its eastern extreme, near San José de las Lajas (La Habana), is linked to another active fault which represents a seismoactive knot responsible for the earthquake of March 9, 1995 (I = 5 degrees, MSK). Seismic events of the Western Cuban region are related to the transpressive interaction of the North American and Caribbean Plates, damped by oceanic structures.     

Gravity Spectra from the Density Distribution of Earth’s Uppermost 435 km

The Earth masses reside in a near-hydrostatic equilibrium, while the deviations are, for example, manifested in the geoid, which is nowadays well determined by satellite gravimetry. Recent progress in estimating the density distribution of the Earth allows us to examine individual Earth layers and to directly see how the sum approaches the observed anomalous gravitational field. This study evaluates contributions from the crust and the upper mantle taken from the LITHO1.0 model and quantifies the gravitational spectra of the density structure to the depth of 435 km. This is done without isostatic adjustments to see what can be revealed with models like LITHO1.0 alone. At the resolution of 290 km (spherical harmonic degree 70), the crustal contribution starts to dominate over the upper mantle and at about 150 km (degree 130) the upper mantle contribution is nearly negligible. At the spatial resolution \(<150\,\hbox {km},\) the spectra behavior is driven by the crust, the mantle lid and the asthenosphere. The LITHO1.0 model was furthermore referenced by adding deeper Earth layers from ak135, and the gravity signal of the merged model was then compared with the observed satellite-only model GOCO05s. The largest differences are found over the tectonothermal cold and old (such as cratonic), and over warm and young areas (such as oceanic ridges). The misfit encountered comes from the mantle lid where a velocity–density relation helped to reduce the RMS error by 40%. Global residuals are also provided in terms of the gravitational gradients as they provide better spatial localization than gravity, and there is strong observational support from ESA’s satellite gradiometry mission GOCE down to the spatial resolution of 80–90 km.     

Relocation of the Ms≥2.0 Earthquakes in the Northern Tianshan Region, Xinjiang, Using the Double-Difference Earthquake Relocation Algorithm

We applied the double-difference earthquake relocation algorithm to 1348 earthquakes with M_S≥2.0 that occurred in the northern Tianshan region, Xinjiang, from April 1988 to June 2003, using a total of 28701 P- and S-wave arrival times recorded by 32 seismic stations in Xinjiang. Aiming to obtain most of these M_S≥2.0 earthquakes relocations, and considering the requirements of the DD method and the condition of data, we added the travel time data of another 437 earthquakes with 1.5≤M_S<2.0. Finally, we obtained the relocation results for 1253 earthquakes with M_S≥2.0, which account for 93% of all the 1348 earthquakes with M_S≥2.0 and includes all the M_S≥3.0 earthquakes. The reason for not relocating the 95 earthquakes with 2.0≤M_S<3.0 is analyzed in the paper. After relocation, the RMS residual decreased from 0.83s to 0.14s, the average error is 0.993 km in E-W direction, 1.10 km in N-S direction, and 1.33 km in vertical direction. The hypocenter depths are more convergent than before and distributed from 5 km to 35 km, with 94% being from 5km to 35 km, 68.2% from 10 km to 25 km. The average hypocenter depth is 19 km.     

Relocation of the Ms≥2.0 Earthquakes in the Northern Tianshan Region, Xinjiang, Using the Double-Difference Earthquake Relocation Algorithm

We applied the double-difference earthquake relocation algorithm to 1348 earthquakes with M_S≥2.0 that occurred in the northern Tianshan region, Xinjiang, from April 1988 to June 2003, using a total of 28701 P- and S-wave arrival times recorded by 32 seismic stations in Xinjiang. Aiming to obtain most of these M_S≥2.0 earthquakes relocations, and considering the requirements of the DD method and the condition of data, we added the travel time data of another 437 earthquakes with 1.5≤M_S<2.0. Finally, we obtained the relocation results for 1253 earthquakes with M_S≥2.0, which account for 93% of all the 1348 earthquakes with M_S≥2.0 and includes all the M_S≥3.0 earthquakes. The reason for not relocating the 95 earthquakes with 2.0≤M_S<3.0 is analyzed in the paper. After relocation, the RMS residual decreased from 0.83s to 0.14s, the average error is 0.993 km in E-W direction, 1.10 km in N-S direction, and 1.33 km in vertical direction. The hypocenter depths are more convergent than before and distributed from 5 km to 35 km, with 94% being from 5km to 35 km, 68.2% from 10 km to 25 km. The average hypocenter depth is 19 km.