Deep submarine pyroclastic eruptions: theory and predicted landforms and deposits

Submarine pyroclastic eruptions at depths greater than a few hundred meters are generally considered to be rare or absent because the pressure of the overlying water column is sufficient to suppress juvenile gas exsolution so that magmatic disruption and pyroclastic activity do not occur. Consideration of detailed models of the ascent and eruption of magma in a range of sea floor environments shows, however, that significant pyroclastic activity can occur even at depths in excess of 3000 m. In order to document and illustrate the full range of submarine eruption styles, we model several possible scenarios for the ascent and eruption of magma feeding submarine eruptions: (1) no gas exsolution; (2) gas exsolution but no magma disruption; (3) gas exsolution, magma disruption, and hawaiian-style fountaining; (4) volatile content builds up in the magma reservoir leading to hawaiian eruptions resulting from foam collapse; (5) magma volatile content insufficient to cause fragmentation normally but low rise speed results in strombolian activity; and (6) volatile content builds up in the top of a dike leading to vulcanian eruptions. We also examine the role of bulk-interaction steam explosivity and contact-surface steam explosivity as processes contributing to volcaniclastic formation in these environments. We concur with most earlier workers that for magma compositions typical of spreading centers and their vicinities, the most likely circumstance is the quiet effusion of magma with minor gas exsolution, and the production of somewhat vesicular pillow lavas or sheet flows, depending on effusion rate. The amounts by which magma would overshoot the vent in these types of eruptions would be insufficient to cause any magma disruption. The most likely mechanism of production of pyroclastic deposits in this environment is strombolian activity, due to the localized concentration of volatiles in magma that has a low rise rate; magmatic gas collects by bubble coalescence, and ascends in large isolated bubbles which disrupt the magma surface in the vent, producing localized blocks, bombs, and pyroclastic deposits. Another possible mode of occurrence of pyroclastic deposits results from vulcanian eruptions; these deposits, being characterized by the dominance of angular blocks of country rocks deposited in the vicinity of a crater, should be easily distinguishable from strombolian and hawaiian eruptions. However, we stress that a special case of the hawaiian eruption style is likely to occur in the submarine environment if magmatic gas buildup occurs in a magma reservoir by the upward drift of gas bubbles. In this case, a layer of foam will build up at the top of the reservoir in a sufficient concentration to exceed the volatile content necessary for disruption and hawaiian-style activity; the deposits and landforms are predicted to be somewhat different from those of a typical primary magmatic volatile-induced hawaiian eruption. Specifically, typical pyroclast sizes might be smaller; fountain heights may exceed those expected for the purely magmatic hawaiian case; cooling of descending pyroclasts would be more efficient, leading to different types of proximal deposits; and runout distances for density flows would be greater, potentially leading to submarine pyroclastic deposits surrounding vents out to distances of tens of meters to a kilometer. In addition, flows emerging after the evacuation of the foam layer would tend to be very depleted in volatiles, and thus extremely poor in vesicles relative to typical flows associated with hawaiian-style eruptions in the primary magmatic gas case. We examine several cases of reported submarine volcaniclastic deposits found at depths as great as 3000 m and conclude that submarine hawaiian and strombolian eruptions are much more common than previously suspected at mid-ocean ridges. Furthermore, the latter stages of development of volcanic edifices (seamounts) formed in submarine environments are excellent candidates for a wide range of submarine pyroclastic activity due not just to the effects of decreasing water depth, but also to: (1) the presence of a summit magma reservoir, which favors the buildup of magmatic foams (enhancing hawaiian-style activity) and episodic dike emplacement (which favors strombolian-style eruptions); and (2) the common occurrence of alkalic basalts, the CO₂ contents of which favor submarine explosive eruptions at depths greater than tholeiitic basalts. These models and predictions can be tested with future sampling and analysis programs and we provide a checklist of key observations to help distinguish among the eruption styles.     

Exploration and discovery in Yellowstone Lake: results from high-resolution sonar imaging, seismic reflection profiling, and submersible studies

‘No portion of the American continent is perhaps so rich in wonders as the Yellow Stone’ (F.V. Hayden, September 2, 1874)Discoveries from multi-beam sonar mapping and seismic reflection surveys of the northern, central, and West Thumb basins of Yellowstone Lake provide new insight into the extent of post-collapse volcanism and active hydrothermal processes occurring in a large lake environment above a large magma chamber. Yellowstone Lake has an irregular bottom covered with dozens of features directly related to hydrothermal, tectonic, volcanic, and sedimentary processes. Detailed bathymetric, seismic reflection, and magnetic evidence reveals that rhyolitic lava flows underlie much of Yellowstone Lake and exert fundamental control on lake bathymetry and localization of hydrothermal activity. Many previously unknown features have been identified and include over 250 hydrothermal vents, several very large (>500 m diameter) hydrothermal explosion craters, many small hydrothermal vent craters (1–200 m diameter), domed lacustrine sediments related to hydrothermal activity, elongate fissures cutting post-glacial sediments, siliceous hydrothermal spire structures, sublacustrine landslide deposits, submerged former shorelines, and a recently active graben. Sampling and observations with a submersible remotely operated vehicle confirm and extend our understanding of the identified features. Faults, fissures, hydrothermally inflated domal structures, hydrothermal explosion craters, and sublacustrine landslides constitute potentially significant geologic hazards. Toxic elements derived from hydrothermal processes also may significantly affect the Yellowstone ecosystem.     

An Mw based earthquake Catalogue for central, northern and northwestern Europe using a hierarchy of magnitude conversions

Data from 25 local catalogues and 30special studies of earthquakes in central,northern and northwestern Europe have beenincorporated into a Databank. The dataprocessing includes discriminating eventtypes, eliminating fake events and dupletsand converting different magnitudes andintensities to M_w if this is not givenby the original source. The magnitudeconversion is a key task of the study andimplies establishment of regressionequations where no local relations exist.The Catalogue contains tectonic events fromthe Databank within the area44°N–72°N,25°W–32°E and the time period1300–1993. The lower magnitude level forthe Catalogue entries is setat M_w == 3.50. The area covered by thedifferent catalogues are associated withpolygons. Within each polygon only datafrom one or a small number of the localcatalogues, supplemented by data fromspecial studies, enter the Catalogue. Ifthere are two or more such catalogues orstudies providing a solution for an event,a priority algorithm selects one entry forthe Catalogue. Then M_w is calculatedfrom one of the magnitude types, or frommacroseismic data, given by the selectedentry according to another priority scheme.The origin time, location, M_w magnitude and reference are specified for eachentry of the Catalogue. So is theepicentral intensity, I₀, if providedby the original source. Following thesecriteria, a total of about 5,000earthquakes constitute the Catalogue.Although originally derived for the purposeof seismic hazard calculation within GSHAP,the Catalogue provides a data base for manytypes of seismicity and seismic hazardstudies.     

On the circulation in the East Frisian Wadden Sea: numerical modeling and data analysis

 In this paper we use a combination of numerical modeling and data analysis to gain a better understanding of the major characteristics of the circulation in the East Frisian Wadden Sea. In particular, we concentrate on the asymmetry of the tidal wave and its modulation in the coastal area, which results in a complex pattern of responses to the sea-level forcing from the North Sea. The numerical simulations are based on the 3-D primitive equation General Estuarine Transport Model (GETM) with a horizontal resolution of 200 m and terrain-following vertical coordinates. The model is forced at its open boundaries with sea-level data from an operational model for the German Bight (German Hydrographic Office). The validation data for our model simulations include time series of tidal gauge data and surface currents measured at a pile in the back-barrier basin of the Island Langeoog, as well as several ADCP transects in the Accumer Ee tidal inlet. Circulation and turbulence characteristics are investigated for typical situations driven by spring and neap tides, and the analysis is focused on dominating temporal and spatial patterns. By investigating the response of five back-barrier basins with rather different morphologies to external forcing, an attempt is made to elucidate the dominating physical balances controlling the circulation in the individual sub-basins. It is demonstrated that the friction at the seabed tends to slow down the tidal signal in the shallow water. This leads to the establishment of flood dominance in the shallow sea north of the barrier islands. South of the islands, where the water volume of the channels at low tide is smaller than the tidal prism, the asymmetry of the tidal signal is shifted towards ebb dominance, a feature which is particularly pronounced at spring tide. At the northern open boundary, the tidal wave propagating from west to east generates a sea-level difference of ∼1 m along the boundary, and thereby triggers vigorous alongshore currents. The frictional control in the model is located in the inlets, as well as along the northern boundary. The correlation between velocity and turbulent kinetic energy tends to the establishment of a net southward transport, giving theoretical support to the observed accumulation of sediments on the intertidal flats. Weak turbulence along the northern shores of the barrier islands and the small magnitude of the residual currents there promote accumulation of suspended matter in these areas, although wave action will generally counteract this effect. Received: 29 May 2002 / Accepted: 26 September 2002 Responsible Editor: Jean-Marie Beckers Acknowledgements We are indebted to S. Dick for providing the data from the operational model of BSH and to B. Flemming for the useful discussions. The topography data and Fig. 1 have been prepared in cooperation with F. Meyer. Figure 2 has been prepared by G. Brink-Spalink. We also thank for the comments from an anonymous reviewer which helped to improve our paper.     

Earth-mass planets in exoplanetary systems

No exoplanets with masses anywhere near as small as that of the Earth have yet been found. Barrie Jones considers whether they could, nevertheless, exist.     

Characterization of Earthquake Strong Ground Motion

— Some underwater landslides are triggered by strong ground motions caused by earthquakes. This paper reviews current concepts and trends in the characterization of strong ground motion. Improved empirical ground motion models have been derived from a strong motion data set that has grown markedly over the past decade. However, these empirical models have a large degree of uncertainty because the magnitude-distance-soil category parameterization of these models often oversimplifies reality. This reflects the fact that other conditions that are known to have an important influence on strong ground motions, such as near-fault rupture directivity effects, crustal waveguide effects, and basin response effects, are not treated as parameters of these simple models. Numerical ground motion models based on seismological theory that include these additional effects have been developed and extensively validated against recorded ground motions, and used to estimate the ground motions of past earthquakes and predict the ground motions of future scenario earthquakes.     

Calcrete profile development in Quaternary alluvial sequences,southeast Spain: implications for using calcretes as a basis for landform chronologies

A detailed study of the morphology and micro‐morphology of Quaternary alluvial calcrete profiles from the Sorbas Basin shows that calcretes may be morphologically simple or complex. The ‘simple’ profiles reflect pedogenesis occurring after alluvial terrace formation and consist of a single pedogenic horizon near the land surface. The ‘complex’ profiles reflect the occurrence of multiple calcrete events during terrace sediment aggradation and further periods of pedogenesis after terrace formation. These ‘complex’ calcrete profiles are consequently described as composite profiles. The exact morphology of the composite profiles depends upon: (1) the number of calcrete‐forming events occurring during terrace sediment aggradation; (2) the amount of sediment accretion that occurs between each period of calcrete formation; and (3) the degree of pedogenesis after terrace formation. Simple calcrete profiles are most useful in establishing landform chronologies because they represent a single phase of pedogenesis after terrace formation. Composite profiles are more problematic. Pedogenic calcretes that form within them may inherit carbonate from calcrete horizons occurring lower down in the terrace sediments. In addition erosion may lead to the exhumation of older calcretes within the terrace sediment. Calcrete ‘inheritance’ may make pedogenic horizons appear more mature than they actually are and produce horizons containing carbonate embracing a range of ages. Calcrete exhumation exposes calcrete horizons whose morphology and radiometric ages are wholly unrelated to terrace surface age. Composite profiles are, therefore, only suitable for chronological studies if the pedogenic horizon capping the terrace sequence can be clearly distinguished from earlier calcrete‐forming events. Thus, a detailed morphological/micro‐morphological study is required before any chronological study is undertaken. This is the only way to establish whether particular calcrete profiles are suitable for dating purposes. Copyright © 2003 John Wiley & Sons, Ltd.     

Testing a model for predicting the timing and location of shallow landslide initiation in soil‐mantled landscapes

The growing availability of digital topographic data and the increased reliability of precipitation forecasts invite modelling efforts to predict the timing and location of shallow landslides in hilly and mountainous areas in order to reduce risk to an ever‐expanding human population. Here, we exploit a rare data set to develop and test such a model. In a 1·7 km² catchment a near‐annual aerial photographic coverage records just three single storm events over a 45 year period that produced multiple landslides. Such data enable us to test model performance by running the entire rainfall time series and determine whether just those three storms are correctly detected. To do this, we link a dynamic and spatially distributed shallow subsurface runoff model (similar to TOPMODEL) to an in?nite slope model to predict the spatial distribution of shallow landsliding. The spatial distribution of soil depth, a strong control on local landsliding, is predicted from a process‐based model. Because of its common availability, daily rainfall data were used to drive the model. Topographic data were derived from digitized 1 : 24 000 US Geological Survey contour maps. Analysis of the landslides shows that 97 occurred in 1955, 37 in 1982 and ?ve in 1998, although the heaviest rainfall was in 1982. Furthermore, intensity–duration analysis of available daily and hourly rainfall from the closest raingauges does not discriminate those three storms from others that did not generate failures. We explore the question of whether a mechanistic modelling approach is better able to identify landslide‐producing storms. Landslide and soil production parameters were ?xed from studies elsewhere. Four hydrologic parameters characterizing the saturated hydraulic conductivity of the soil and underlying bedrock and its decline with depth were ?rst calibrated on the 1955 landslide record. Success was characterized as the most number of actual landslides predicted with the least amount of total area predicted to be unstable. Because landslide area was consistently overpredicted, a threshold catchment area of predicted slope instability was used to de?ne whether a rainstorm was a signi?cant landslide producer. Many combinations of the four hydrological parameters performed equally well for the 1955 event, but only one combination successfully identi?ed the 1982 storm as the only landslide‐producing storm during the period 1980–86. Application of this parameter combination to the entire 45 year record successfully identi?ed the three events, but also predicted that two other landslide‐producing events should have occurred. This performance is signi?cantly better than the empirical intensity–duration threshold approach, but requires considerable calibration effort. Overprediction of instability, both for storms that produced landslides and for non‐producing storms, appears to arise from at least four causes: (1) coarse rainfall data time scale and inability to document short rainfall bursts and predict pressure wave response; (2) absence of local rainfall data; (3) legacy effect of previous landslides; and (4) inaccurate topographic and soil property data. Greater resolution of spatial and rainfall data, as well as topographic data, coupled with systematic documentation of landslides to create time series to test models, should lead to signi?cant improvements in shallow landslides forecasting. Copyright © 2003 John Wiley & Sons, Ltd.