The dynamics and thermodynamics of large ash flows

 Ash flow deposits, containing up to 1000 km³ of material, have been produced by some of the largest volcanic eruptions known. Ash flows propagate several tens of kilometres from their source vents, produce extensive blankets of ash and are able to surmount topographic barriers hundreds of metres high. We present and test a new model of the motion of such flows as they propagate over a near horizontal surface from a collapsing fountain above a volcanic vent. The model predicts that for a given eruption rate, either a slow (10–100 m/s) and deep (1000–3000 m) subcritical flow or a fast (100–200 m/s) and shallow (500–1000 m) supercritical flow may develop. Subcritical ash flows propagate with a nearly constant volume flux, whereas supercritical flows entrain air and become progressively more voluminous. The run-out distance of such ash flows is controlled largely by the mass of air mixed into the collapsing fountain, the degree of fragmentation and the associated rate of loss of material into an underlying concentrated depositional system, and the mass eruption rate. However, in supercritical flows, the continued entrainment of air exerts a further important control on the flow evolution. Model predictions show that the run-out distance decreases with the mass of air entrained into the flow. Also, the mass of ash which may ascend from the flow into a buoyant coignimbrite cloud increases as more air is entrained into the flow. As a result, supercritical ash flows typically have shorter runout distances and more ash is elutriated into the associated coignimbrite eruption columns. We also show that one-dimensional, channellized ash flows typically propagate further than their radially spreading counterparts. As a Plinian eruption proceeds, the erupted mass flux often increases, leading to column collapse and the formation of pumiceous ash flows. Near the critical conditions for eruption column collapse, the flows are shed from high fountains which entrain large quantities of air per unit mass. Our model suggests that this will lead to relatively short ash flows with much of the erupted material being elutriated into the coignimbrite column. However, if the mass flux subseqently increases, then less air per unit mass is entrained into the collapsing fountain, and progressively larger flows, which propagate further from the vent, will develop. Our model is consistent with observations of a number of pyroclastic flow deposits, including the 1912 eruption of Katmai and the 1991 eruption of Pinatubo. The model suggests that many extensive flow sheets were emplaced from eruptions with mass fluxes of 10⁹–10¹⁰ kg/s over periods of 10³–10⁵ s, and that some indicators of flow "mobility" may need to be reinterpreted. Furthermore, in accordance with observations, the model predicts that the coignimbrite eruption columns produced from such ash flows rose between 20 and 40 km. Received: 25 August 1995 / Accepted: 3 April 1996     

Potential cauchy-poisson waves generated by submarine eruptions of kick 'em Jenny volcano

Kick 'em Jenny is the only known currently active submarine volcano in the Lesser Antilles. The volcano has erupted at least 10 times since first being discovered in 1939 and the summit has shoaled from a depth of 232 m in 1962 to its present-day depth of 150 m. Kick 'em Jenny is located in a province of explosive volcanism, has a known history of explosive eruptions and erupts magma of an explosive type. Future eruptions are likely to become increasingly more violent as the effect of the overlying water pressure becomes less. A preliminary study (Smith and Shepherd, 1993) suggests that Kick 'em Jenny is a prime candidate for tsunamigenic eruptions on a potentially hazardous scale, possibly affecting the whole of the eastern Caribbean region.The classic approach to problems of water waves generated by sudden disturbances of the free surface makes use of the Cauchy-Poisson-Lamb theory. A large number of theoretical developments to this theory have been made for specific forms of surface disturbance. A development by Unoki and Nakano (1953a, b) considers both two- and three-dimensional Cauchy-Poisson waves generated by finite initial elevations and impulses applied to a free surface of infinitely deep water. Unoki and Nakano's results compared well to the wave systems recorded following submarine eruptions of the Myojinsho Reef volcano in 1952–53.Given the similarity of the two situations, Unoki and Nakano's theory is applied to Kick 'em Jenny to provide estimates of potential Cauchy-Poisson wave heights throughout the eastern Caribbean for a range of eruption magnitudes. The results show that, although the waves generated are unlikely to pose much of a threat to the eastern Caribbean as a whole, they should be considered a hazard to the islands immediately adjacent to the volcano including Grenada, the Grenadines, and St Vincent.     

The Cauvery Shear Zone, Southern Granulite Terrain, India: A crustal-scale flower structure

The Cauvery Shear Zone (CSZ) is a crustal-scale shear system within the Southern Granulite Terrain along the southern margin of the Archaean Dharwar craton. Structural interpretation of satellite data and field observations reveal four major shear zones within the CSZ system. They show dextral shear kinematics synchronous with a major Neoproterozoic tectono-metamorphic event (D2) associated with intracrustal melting and migmatisation. The disposition, geometry and contemporaneity of shear fabrics of the CSZ system are modelled in terms of a crustal-scale flower structure akin to transpressional and collisional orogens. In the light of recent seismic evidence for a displaced Moho structure and a mid- to lower-crustal low velocity zone, the flower structure across the CSZ may extend to mantle depths.     

The Tornquist Zone, a north east inclining lithospheric transition at the south western margin of the Baltic Shield: Revealed through a nonlinear teleseismic tomographic inversion

From July 1996 to August 1997 the TOR project operated 130 seismographs in North Germany, Denmark and South Sweden, with the aim of collecting signals from local, regional and teleseismic earthquakes. This data set is particularly interesting since the seismic antenna crosses the most significant geological boundary in Europe, the Tornquist Zone, which in the northern part is the border between the Baltic Shield and the younger European lithosphere. Previous studies have shown significant physical changes in the crust and upper mantle across this transition zone, including two independent teleseismic tomographic studies of the TOR data set. But these two studies disagree on the orientation of the slope of the transition. Both studies used an iterative linearized inversion method. We will in this work Preprint submitted to Elsevier Science 27 July 2005 present an inversion based on Bayesian statistics, where the solution space is examined in order to study a very large number of tomographic solutions and to examine the solution uniqueness and uncertainty. The method is applied to measurements of 3345 relative teleseismic P-phase travel times from 48 teleseismic earthquakes with good azimuthal coverage with respect to the great circle arc of the TOR array. We find the lithospheric transition to be a north east inclination of around 30° to 45° off vertical.     

Geodynamic significance of granitoid magmatism in the southeast Anatolian orogen: geochemical and geochronogical evidence from Göksun–Afşin (Kahramanmaraş, Turkey) region

In southeast Anatolia, there are number of tectonomagmatic units in the Kahramanmaraş–Malatya–Elazığ region that are important in understanding the geological evolution of the southeast Anatolian orogenic belt during the Late Cretaceous. These are (a) metamorphic massifs, (b) ophiolites, (c) ophiolite-related metamorphics and (d) granitoids. The granitoids (i.e. Göksun–Afşin in Kahramanmaraş, Doğanşehir in Malatya and Baskil in Elazığ) intrude all the former units in a NE–SW trending direction. The granitoid in Göksun–Afşin (Kahramanmaraş) region is mainly composed of granodioritic and granitic in composition. The granodiorite contains a number of amphibole-bearing mafic microgranular enclaves of different sizes, whereas the granite is intruded by numerous aplitic dikes. The granitoid rocks have typical calcalkaline geochemical features. The REE- and Ocean ridge granite-normalized multi-element patterns and tectonomagmatic discrimination diagrams, as well as biotite geochemistry suggest that the granitoids were formed in a volcanic arc setting. The K–Ar geochronology of the granitoid rocks yielded ages ranging from 85.76±3.17 to 77.49±1.91 Ma. The field, geochemical and geochronological data suggest the following Late Cretaceous tectonomagmatic scenario for southeast Anatolia. The ophiolites were formed in a suprasubduction zone tectonic setting whereas the ophiolite-related metamorphic rocks formed either during the initiation of intraoceanic subduction or late-thrusting (∼90 Ma). These units were then overthrust by the Malatya–Keban platform during the progressive elimination of the southern Neotethys. Thrusting of the Malatya–Keban platform over the ophiolites and related metamorphic rocks was followed by the intrusion of the granitoids (88–85 Ma) along the Tauride active continental margin in the southern Neotethys.     

Post-mid-Cretaceous eastern North Sea evolution inferred from 3D thermo-mechanical modelling

The Late Cretaceous–Cenozoic evolution of the eastern North Sea region is investigated by 3D thermo-mechanical modelling. The model quantifies the integrated effects on basin evolution of large-scale lithospheric processes, rheology, strength heterogeneities, tectonics, eustasy, sedimentation and erosion.

The evolution of the area is influenced by a number of factors: (1) thermal subsidence centred in the central North Sea providing accommodation space for thick sediment deposits; (2) 250-m eustatic fall from the Late Cretaceous to present, which causes exhumation of the North Sea Basin margins; (3) varying sediment supply; (4) isostatic adjustments following erosion and sedimentation; (5) Late Cretaceous–early Cenozoic Alpine compressional phases causing tectonic inversion of the Sorgenfrei–Tornquist Zone (STZ) and other weak zones.

The stress field and the lateral variations in lithospheric strength control lithospheric deformation under compression. The lithosphere is relatively weak in areas where Moho is deep and the upper mantle warm and weak. In these areas the lithosphere is thickened during compression producing surface uplift and erosion (e.g., at the Ringkøbing–Fyn High and in the southern part of Sweden). Observed late Cretaceous–early Cenozoic shallow water depths at the Ringkøbing–Fyn High as well as Cenozoic surface uplift in southern Sweden (the South Swedish Dome (SSD)) are explained by this mechanism.

The STZ is a prominent crustal structural weakness zone. Under compression, this zone is inverted and its surface uplifted and eroded. Contemporaneously, marginal depositional troughs develop. Post-compressional relaxation causes a regional uplift of this zone.

The model predicts sediment distributions and paleo-water depths in accordance with observations. Sediment truncation and exhumation at the North Sea Basin margins are explained by fall in global sea level, isostatic adjustments to exhumation, and uplift of the inverted STZ. This underlines the importance of the mechanisms dealt with in this paper for the evolution of intra-cratonic sedimentary basins.