Granular flow experiments on the interaction with stationary runout path materials and comparison to rock avalanche events

The central focus of this work is to study the processes acting well below the surface of a moving rock or debris avalanche during travel over stationary substrate material. Small‐scale physical models at a linear scale of 1:10⁴ used coal as avalanche analogue material and different granular material simulating sedimentary substrates varying in frictional resistance, thickness and relative basal boundary roughness, as well as inerodible, non‐deformable runout path conditions. Substrate materials with the least frictional resistance showed the greatest response to granular flow overriding, becoming entirely mobilized beneath and ahead of the moving mass and producing the longest runout observed with a unique deposit profile shape. With a smooth substrate basal contact, failure occurred along this plane and avalanche and substrate became coupled during runout. With a rough base, however, temporary force chains of grain contacts in the substrate prevailed longer, imparted their resistance to motion/shear into the granular flow, and the flow rear section consequently halted earlier than when moving over substrates with a weak base. Reducing substrate thickness diminished the effect of basal contact roughness on granular flow runout and deposit length. Inerodible, non‐deformable substrate conditions caused changes in granular flow behaviour from essentially en masse sliding on low‐friction surfaces to increasing granular agitation over rougher paths. Copyright © 2012 John Wiley & Sons, Ltd.     

Insights into rock avalanche emplacement processes from detailed morpho‐lithological studies of the Tschirgant deposit (Tyrol,Austria)

Large, rapid rockslope failures generate deposits with complex morphologies due to a number of causal and influencing factors. To investigate these, we conducted a detailed case study at the carbonate Tschirgant deposit (Tyrol, Austria). It preserved evidence of simultaneous rock sliding (very large, coherent hummocks) and rock avalanche spreading (smaller, more scattered hummocks and ridges). Motion indicators, such as longitudinal ridges furthermore pinpoint the transition between linear sliding and radial spreading. The lithological distribution in the Tschirgant deposit shows that it retained source stratigraphy despite being split into two accumulation lobes by a high bedrock ridge. Furthermore, lithology had a very strong control on the final deposit morphology in that the different lithologic units form individual deposit surfaces. River erosion has created fortuitous outcrops that reveal the basal rock avalanche contact. The underlying valley‐fill sediments (substrates) have been intricately involved in shaping the rock avalanche morphology and, where entrained, highlight internal rock avalanche deformation features. This study shows that intrinsic dynamic properties of granular media (e.g. tendency for longitudinal alignments), emplacement mode, lithology (and source predisposition), runout path topography, and substrates form the quintet of causal factors that shape rock avalanche morphology. Copyright © 2015 John Wiley & Sons, Ltd.     

Paleomagnetic estimate of the emplacement temperature of the long-runout Nevado de Colima volcanic debris avalanche deposit, Mexico

Stoopes and Sheridan have mapped a volcanic debris avalanche of Nevado de Colima which has an exceptionally long runout (120 km) and low fall-height to length ratio (H/L = 0.04). We present paleomagnetic results from this volcanic debris avalanche deposit which provide evidence that this avalanche was emplaced at elevated temperatures. The majority of samples, collected from lithic clasts in the volcanic debris avalanche deposit, exhibit two-component remanent magnetizations with a low-temperature component (25–350°C) which is well grouped about the geomagnetic field direction at Colima and a high-temperature component (350–580°C) which is randomly oriented. Although the temperature of the deposit most likely varied with distance from the volcanic source and the thickness of the deposit, our results suggest an emplacement temperature of approximately 350°C at intermediate distances (18–26 km) from the source. In order for the rock clasts (20–40 cm diameter) to be heated to these temperatures, the avalanche was most likely the results of a magmatic, Bezymianny-type eruption. The mixing of hot, juvenile gases with the clasts provides an explanation for the high degree of fluidization of this material, as evidenced by the long runout of this avalanche deposit.     

Discrete Element Modeling of Tangjiagou Two-Branch Rock Avalanche Triggered by the 2013 Lushan MW6.6 Earthquake, China

Two branches of Tangjiagou rock avalanche were triggered by Lushan earthquake in Sichuan Province, China on April 20^th, 2013. The rock avalanche has transported about 1 500 000 m³ of sandstone from the source area. Based on discrete element modeling, this study simulates the deformation, failure and movement process of the rock avalanche. Under seismic loading, the mechanism and process of deformation, failure, and runout of the two branches are similar. In detail, the stress concentration occur firstly on the top of the mountain ridge, and accordingly, the tensile deformation appears. With the increase of seismic loading, the strain concentration zone extends in the forward and backward directions along the slipping surface, forming a locking segment. As a result, the slipping surface penetrates and the slide mass begin to slide down with high speed. Finally, the avalanche accumulates in the downstream and forms a small barrier lake. Modeling shows that a number of rocks on the surface exhibit patterns of horizontal throwing and vertical jumping under strong ground shaking. We suggest that the movement of the rock avalanche is a complicated process with multiple stages, including formation of the two branches, high-speed sliding, transformation into debris flows, further movement and collision, accumulation, and the final steady state. Topographic amplification effects are also revealed based on acceleration and velocity of special monitoring points. The horizontal and vertical runout distances of the surface materials are much greater than those of the internal materials. Besides, the sliding duration is also longer than that of the internal rock mass.     

Impact of an 0.2 km3 Rock Avalanche on Lake Eibsee (Bavarian Alps,Germany) – Part I: Reconstruction of the paleolake and Effects of the Impact

Rock avalanches destroy and reshape landscapes in only a few minutes and are among the most hazardous processes on Earth. The surface morphology of rock avalanche deposits and the interaction with the underlying material are crucial for runout properties and reach. Water within the travel path is displaced, producing large impact waves and reducing friction, leading to long runouts. We hypothesize that the 0.2 km³ Holocene Eibsee rock avalanche from Mount Zugspitze in the Bavarian Alps overran and destroyed Paleolake Eibsee and left a unique sedimentological legacy of processes active during the landslide. We captured 9.5 km of electrical resistivity tomography (ERT) profiles across the rock avalanche deposits, with up to 120 m penetration depth and more than 34 000 datum points. The ERT profiles reveal up to ~50 m thick landslide debris, locally covering up to ~30 m of rock debris with entrained fine-grained sediments on top of isolated remnants of decametre-wide paleolake sediments. The ERT profiles allow us to infer processes involved in the interaction of the rock avalanche with bedrock, lake sediments, and morainal sediments, including shearing, bulging, and bulldozing. Complementary data from drilling, a gravel pit exposure, laboratory tests, and geomorphic features were used for ERT calibration. Sediments overrun by the rock avalanche show water-escape structures. Based on all of these datasets, we reconstructed both position and size of the paleolake prior to the catastrophic event. Our reconstruction of the event contributes to process an understanding of the rock avalanche and future modelling and hazard assessment. Here we show how integrated geomorphic, geophysical, and sedimentological approaches can provide detailed insights into the impact of a rock avalanche on a lake. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Catastrophic rock avalanche 3600 years BP from El Capitan,Yosemite Valley,California

Large rock slope failures from near‐vertical cliffs are an important geomorphic process driving the evolution of mountainous landscapes, particularly glacially steepened cliffs. The morphology and age of a 2·19 × 10⁶ m³ rock avalanche deposit beneath El Capitan in Yosemite Valley indicates a massive prehistoric failure of a large expanse of the southeast face. Geologic mapping of the deposit and the cliff face constrains the rock avalanche source to an area near the summit of ～8·5 × 10⁴ m². The rock mass free fell ～650 m, reaching a maximum velocity of 100 m s^?1, impacted the talus slope and spread across the valley floor, extending 670 m from the base of the cliff. Cosmogenic beryllium‐10 exposure ages from boulders in the deposit yield a mean age of 3·6 ± 0·2 ka. The ～13 kyr time lag between deglaciation and failure suggests that the rock avalanche did not occur as a direct result of glacial debuttressing. The ～3·6 ka age for the rock avalanche does coincide with estimated late Holocene rupture of the Owens Valley fault and/or White Mountain fault between 3·3 and 3·8 ka. The coincidence of ages, combined with the fact that the most recent (AD 1872) Owens Valley fault rupture triggered numerous large rock falls in Yosemite Valley, suggest that a large magnitude earthquake (≥M7.0) centered in the south‐eastern Sierra Nevada may have triggered the rock avalanche. If correct, the extreme hazard posed by rock avalanches in Yosemite Valley remains present and depends on local earthquake recurrence intervals. Published in 2010 by John Wiley & Sons, Ltd.     

Surface exposure dating of the Veliki vrh rock avalanche in Slovenia associated with the 1348 earthquake

Over 30 samples from bedrock and boulders from the Veliki vrh rock avalanche have been collected for surface exposure dating. The limestone rocks have been radiochemically treated to isolate and determine long-lived ³⁶Cl by accelerator mass spectrometry. It could be shown that the Veliki vrh rock avalanche from the Košuta Mountain (Slovenia) event can be very likely linked to one of the major historical earthquakes in Europe happening on the 25^th of January 1348. Taken into account independently determined denudation rates, inherited ³⁶Cl originating from pre-exposure at shallow depths (20–55 m) could be calculated. The high amount of inherited ³⁶Cl, i.e. 17–46% of the total ³⁶Cl, makes this site not suitable for a precise determination of the ³⁶Cl production rate as it was originally anticipated. Veliki vrh is a “classic” rock avalanche of high velocity. The slope failed in the upper part with a translational slide predominantly along the bedding planes, whereas dynamic fragmentation is the cause for further crushing of the material and the long runout.     

Morphology and emplacement of an unusual debris-avalanche deposit at Jocotitlán volcano,Central Mexico

A pre-historic collapse of the northeastern flank of Jocotitlán Volcano (3950 m), located in the central part of the Trans Mexican Volcanic Belt, produced a debris-avalanche deposit characterized by surficial hummocks of exceptional size and conical shape. The avalanche covered an area of 80 km², had an apparent coefficient of friction (H/L)_of 0.11, a maximum runout distance of 12 km, and an estimated volume of 2.8 km³. The most remarkable features of the Jocotitlán debris avalanche deposit are: the several steep (29–32°) conical proximal hummocks (up to 165 m high), large tansverse ridges (up to 205 m high and 2.7 km long) situated at the base of the volcano, and the steep 15–50 m thick terminal scarp. Proximal conical hummocks and parallel ridges that can be visually fitted back to their pre-collapse position on the mountain resulted from a sliding mode of emplacement. Steep primary slopes developed as a result of the accumulation of coarse angular clasts at the angle of repose around core clasts that are decameters in size. Distal hummocks are commonly smaller, less conical, and clustered with more diffuse outlines. Field evidence indicates that the leading distal edge of the avalanche spilled around certain topographic barriers and that the distal moving mass had a yield strength prior to stopping. In the NE sector, the avalanche was suddenly confined by topographically higher lacustrine and volcaniclastic deposits which as a result were intensely thrust-faulted, folded, and impacted by large clasts that separated from the avalanche front. Post-emplacement loading also induced normal faulting of these soft, locally water-rich sediments. The regional tectonic pattern, N-NE direction of flank failure, and the presence of a major normal fault which intersects the volcano and is parallel to the orientation of the Acambay graben located 10 km to the N suggest a genetic relationship between the extensional tectonic stress regime and triggering of catastrophic slope failure. The presence of a 3-m-thick sequence of pumice and obsidian-rich pyroclastic surge and fall tephra directly overlying the debris-avalanche deposit indicates that magma must have been present within the edifice just prior to the catastrophic flank failure. The breached crater left by the avalanche has mostly been filled by dacitic domes and lava flows. The youngest pryroclastic surge deposits on the upper flanks of the volcano have an historical C¹⁴ age of 680±80 yearsBp (Ad 1270±80). Thus Jocotitlán volcano, formerly believed to be extinct, should be considered potentially active. Because of its close proximity to Mexico-City (60 km), the most populous city in the world, reactivation could engender severe hazards.     

Multi‐method (14C, 36Cl, 234U/230Th) age bracketing of the Tschirgant rock avalanche (Eastern Alps): implications for absolute dating of catastrophic mass‐wasting

Correct and precise age determination of prehistorical catastrophic rock‐slope failures prerequisites any hypotheses relating this type of mass wasting to past climatic regimes or palaeo‐seismic records. Despite good exposure, easy accessibility and a long tradition of absolute dating, the age of the 230 million m³ carbonate‐lithic Tschirgant rock avalanche event of the Eastern Alps (Austria) still is relatively poorly constrained. We herein review the age of mass‐wasting based on a total of 17 absolute ages produced with three different methods (¹⁴C, ³⁶Cl, ²³⁴U/²³⁰Th). Chlorine‐36 (³⁶Cl) cosmogenic surface exposure dating of five boulders of the rock avalanche deposit indicates a mean event age of 3.06 ± 0.62 ka. Uranium‐234/thorium‐230 (²³⁴U/²³⁰Th) dating of soda‐straw stalactites formed in microcaves beneath boulders indicate mean precipitation ages of three individual soda straws at 3.20 ± 0.26 ka, 3.04 ± 0.10 ka and 2.81 ± 0.15 ka; notwithstanding potential internal errors, these ages provide an ‘older‐than’ (ante quam) proxy for mass‐wasting. Based on radiocarbon ages (nine sites) only, it was previously suggested that the present rock avalanche deposit represents two successive failures (3.75 ± 0.19 ka bp , 3.15 ± 0.19 ka bp ). There is, however, no evidence for two events neither in surface outcrops nor in LiDAR derived imagery and drill logs. The temporal distribution of all absolute ages (¹⁴C, ³⁶Cl, ²³⁴U/²³⁰Th) also does not necessarily indicate two successive events but suggest that a single catastrophic mass‐wasting took place between 3.4 and 2.4 ka bp . Taking into account the maximum age boundary given by reinterpreted radiocarbon datings and the minimum U/Th‐ages of calcite precipitations within the rock avalanche deposits, a most probable event age of 3.01 ± 0.10 ka bp can be proposed. Our results underscore the difficulty to accurately date catastrophic rock slope failures, but also the potential to increase the accuracy of age determination by combining methods. Copyright © 2016 John Wiley & Sons, Ltd.     

Runout of the Socompa volcanic debris avalanche, Chile: a mechanical explanation for low basal shear resistance

We propose a mechanical explanation for the low basal shear resistance (about 50 kPa) previously used to simulate successfully the complex, well-documented deposit morphology and lithological distribution produced by emplacement of the 25 km³ Socompa volcanic debris avalanche deposit, Chile. Stratigraphic evidence for intense basal comminution indicates the occurrence of dynamic rock fragmentation in the basal region of this large granular mass flow, and we show that such fragmentation generates a basal shear stress, retarding motion of the avalanche, that is a function of the flow thickness and intact rock strength. The topography of the Socompa deposit is realistically simulated using this fragmentation-derived resistance function. Basal fragmentation is also compatible with the evidence from the deposit that reflection of the avalanche from topography caused a secondary wave that interacted with the primary flow.     

Landslide Amplification by Liquefaction of Runout‐Path Material after the 2008 Wenchuan (M 8·0) Earthquake,China

Here, we propose that an earthquake can trigger the failure of a landslide mass while simultaneously triggering liquefaction of runout‐path materials before the arrival of the landslide mass, thus greatly increasing the size and mobility of an overriding landslide. During the 2008 Wenchuan earthquake, about 60 000 landslides were triggered, directly resulting in about 20 000 casualties. While these landslides mainly originated from steep slopes, some landslides with high mobility formed in colluvial valley deposits. Among these, the most catastrophic was the Xiejiadian landslide in Pengzhou city, which traveled hundreds of meters before coming to rest. Through field investigation and laboratory testing, we conclude that this landslide primarily formed from colluvial deposits in the valley and secondarily from failure of slopes in granitic rock located uphill. Much of the granitic slope failure was deposited in the upper part of the travel path (near the slide head); the remainder was dispersed throughout the main landslide deposit. Superposition of deposits at the landslide toe indicates that landslide debris derived from colluvial soil was deposited first. The deposits at the landslide toe displayed flow characteristics, such as fine materials comprising basal layers and large boulders covering the deposit surface. We hypothesize that the main part of the landslide resulted from seismogenic liquefaction of valley colluvium, rather than from liquefaction potentially caused by undrained loading from the granitic slope failures impacting the colluvium. To examine the likelihood that seismogenic liquefaction occurred, we took samples from different areas of the landslide deposit and performed undrained cyclic shear tests on them in the laboratory. The results showed that the sandy soils that comprise most of the deposit are highly liquefiable under seismic loading. Therefore, we conclude that liquefaction of the colluvium in the valley during the earthquake was the main reason for this rapid (~46 m/s) long‐runout (1·7 km) landslide. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling the spatio-temporal repartition of right-truncated data: an application to avalanche runout altitudes in Hautes-Savoie

In this paper, we propose a novel approach for generating avalanche hazard maps based on the spatial dependence of avalanche runout altitudes. The right-truncated data are described with a Bayesian hierarchical model in which the spatio-temporal process is assumed to be the sum of independent spatial and temporal terms. Topography is roughly taken into account according to valley altitude and path exposition, and the spatial dependence is modelled with a Matérn covariance function. An application is performed to the Haute-Savoie region, French Alps. A spatial dependence in runout altitudes is identified, and an effective range of about 10 km is inferred. The temporal trend extracted highlights the increase of avalanche runout altitudes from 1955, attributed to both anthropogenic factors and climate warming. In a cross validation scheme, spatial predictions are provided on undocumented paths using kriging equations. All in all, although our model is unable to take into account small topographic features, it is a first-ever approach that produces very encouraging results. It could be enhanced in future work by incorporating a numerical physically-based code into the modelling.     

A New Method for the Estimation of Avalanche Distance Exceeded Probabilities

A crucial point in any methodology for avalanche hazard assessment is the evaluation of avalanche distance exceeded probability, i.e., the annual probability that any assigned location along a given path is reached or exceeded by an avalanche. Typically this problem is faced by estimating the snow volume in the starting zone that is likely to accumulate an average every T years by statistical analysis of snowfall record, and then using this volume as input to an appropriately calibrated avalanche dynamics model to determine the runout distancesfor this design event. This methodology identifies the areas that canbe affected by an avalanche for the considered value of the return period (i.e. the average interval of time for a certain event to repeat itself), ¯T. However, it does not allow us to evaluate the actual avalanche encounter probability for any given point in the runout zone. In the present work this probability is computed by numerical integration of the expression P(x) = ∫₀ ^∞ P^*(V)f(V) dV, where f is the probabilitydensity function (PDF) of the avalanche release volume V, and P^* is the probability of the point x being reached or passed by an avalanche if the release volume is V; this latter probability is calculated by avalanche dynamics simulations. The procedure is implemented using a one-dimensional hydraulic-continuum avalanche dynamic model, calibrated on data from different Italian Alpine ranges, and is applied to a real world hazard mapping problem.     

Geophysical investigation of the Sandalp rock avalanche deposits

In the study of rock avalanche phenomena, numerical modelling makes use of back analyses of the rock avalanche propagation for calibration of the modelling assumptions and parameters. The back analyses require knowledge of the run-out area boundaries and the thickness distribution of the deposit. Geophysical methods can be applied to retrieve the thickness distribution, but, due to strong heterogeneities and logistic problems they are seldom applied. The aim of this work is to assess the potential of integrated geophysical methods to recognise and characterise a deposit created by two rock avalanches which occurred in the Sandalp valley (Switzerland) in 1996. The topography of the site before and after the rock avalanche is known and can be used as a benchmark. Resistivity tomography, seismic P-wave tomography, and active and passive surface wave analysis have been applied on several profiles deployed both on the rock avalanche deposit and in the surrounding area. Innovative approaches for surface wave analysis based on laterally constrained inversion and multimodal inversion have been applied to the data. A comparison of the results of the geophysical investigations with the topographic benchmark has shown the capability of the geophysical methods to locate the bottom of the deposit in the areas where the contrast with the host sediments properties is significant. In these areas, the deposit has higher resistivities and lower seismic velocities than the underlying materials. In the areas where the deposit is thicker and richer in fine-grained materials the geophysical parameters are not able to discriminate between the rock avalanche deposit and the underlying sediments. As a secondary task, the geophysical methods also allowed the bedrock pattern to be outlined.     

Bayesian optimal design of an avalanche dam using a multivariate numerical avalanche model

For snow avalanches, passive defense structures are generally designed by considering high return period events. However, defining a return period turns out to be tricky as soon as different variables are simultaneously considered. This problem can be overcome by maximizing the expected economic benefit of the defense structure, but purely stochastic approaches are not possible for paths with a complex geometry in the runout zone. Therefore, in this paper, we include a multivariate numerical avalanche propagation model within a Bayesian decisional framework. The influence of a vertical dam on an avalanche flow is quantified in terms of local energy dissipation with a simple semi-empirical relation. Costs corresponding to dam construction and the damage to a building situated in the runout zone are roughly evaluated for each dam height–hazard value pair, with damage intensity depending on avalanche velocity. Special attention is given to the poor local information to be taken into account for the decision. Using a case study from the French avalanche database, the Bayesian optimal dam height is shown to be more pessimistic than the classical optimal height because of the increasing effect of parameter uncertainty. It also appears that the lack of local information is especially critical for a building exposed to the most extreme events only. The residual hazard after dam construction is analyzed and the sensitivity to the different modelling assumptions is evaluated. Finally, possible further developments of the approach are discussed.     

使用宽频带地震记录研究2017年6月24日茂县新磨滑坡的动力学过程

北京时间2017年6月24日5时39分左右,四川省茂县叠溪镇新磨村发生大型岩质滑坡.体积约4.3×10⁶ m³的巨型岩体从山顶脱落,顺坡滑行约2.6 km后破碎沉积;碎屑物掩埋了整个新磨村,造成了巨大的人员伤亡和财产损失.本文使用来自滑坡周围的10个地震台站的宽频带观测资料的长周期信号反演了这次滑坡的受力时间函数;同时使用逐步细化的格点搜索方法得到了滑坡的位置,与其真实位置一致;根据反演的受力时间函数计算了滑坡过程中滑体的运动学参数,得到的滑体运动轨迹与实际路径吻合.综合分析地震信号、受力时间函数和运动学参数表明,本次滑坡主运动的持续时间约为79 s;脱落岩体在5∶38∶50.2启动后持续加速,在5∶39∶37.2达到速度峰值,约为52.1 m·s^-1;这段时间内岩体没有明显的破碎;之后,岩体开始铲刮并裹挟古滑坡造成的碎屑沉积物,自身也开始破碎解体,总体开始减速运动,直到5∶40∶9.2主运动停止;此后,小规模的碎屑散落又持续了约10 s的时间.     

A gravitational spreading origin for the Socompa debris avalanche

Socompa Volcano arguably provides the world's best-exposed example of a sector collapse-derived debris avalanche deposit. New observations lead us to re-interpret the origin of the sector collapse. We show that it was triggered by failure of active thrust-anticlines in sediments and ignimbrites underlying the volcano. The thrust-anticlines were a result of gravitational spreading of substrata under the volcano load. About 80% of the resulting avalanche deposit is composed of substrata formerly residing under the volcano and in the anticlines. The collapse scar can be traced up to 5 km from the edifice, truncating two spreading-related anticlines, which collapsed in the event. Outcrops near the volcano preserve evidence of edifice material being carried along on top of mobilised substrata. On the north side of the scar, the avalanche motion was initially at right angles to the failure edge. Structural relations indicate that immediately prior to collapse the substrata disintegrated, became effectively liquidised, and were ejected from beneath the edifice. Catastrophic mobilisation of substrata probably resulted from breakdown of ignimbrite clasts and cement. It may have occurred through progressive rock fracture by high shear strain during spreading. Material ejected from under Socompa formed a layer on which volcanic edifice debris was transported. This interpretation of events explains the puzzling observation that avalanche units with the lowest gravitational potential energy moved the furthest. It can also account for avalanche motion normal to the collapse scar walls. Ignimbrites and other rock types probably capable of similar behaviour underlie many other volcanoes. Identification of spreading at other sites could therefore be a first step towards assessment of the potential for this style of catastrophic sector collapse.     

The flow dynamics of an extremely large volume pyroclastic flow,the 2.08-Ma Cerro Galán Ignimbrite,NW Argentina,and comparison with other flow types

The 2.08-Ma Cerro Galán Ignimbrite (CGI) represents a >630-km³ dense rock equivalent (VEI 8) eruption from the long-lived Cerro Galán magma system (∼6 Ma). It is a crystal-rich (35–60%), pumice (<10% generally) and lithic-poor (<5% generally) rhyodacitic ignimbrite, lacking a preceding plinian fallout deposit. The CGI is preserved up to 80 km from the structural margins of the caldera, but almost certainly was deposited up to 100 km from the caldera in some places. Only one emplacement unit is preserved in proximal to medial settings and in most distal settings, suggesting constant flow conditions, but where the pyroclastic flow moved into a palaeotopography of substantial valleys and ridges, it interacted with valley walls, resulting in flow instabilities that generated multiple depositional units, often separated by pyroclastic surge deposits. The CGI preserves a widespread sub-horizontal fabric, defined by aligned elongate pumice and lithic clasts, and minerals (e.g. biotite). A sub-horizontal anisotropy of magnetic susceptibility fabric is defined by minute magnetic minerals in all localities where it has been analysed. The CGI is poor in both vent-derived (‘accessory’) lithics and locally derived lithics from the ground surface (‘accidental’) lithics. Locally derived lithics are small (<20 cm) and were not transported far from source points. All data suggest that the pyroclastic flow system producing the CGI was characterised throughout by high sedimentation rates, resulting from high particle concentration and suppressed turbulence at the depositional boundary layer, despite being a low aspect ratio ignimbrite. Based on these features, we question whether high velocity and momentum are necessary to account for extensive flow mobility. It is proposed that the CGI was deposited by a pyroclastic flow system that developed a substantial, high particle concentration granular under-flow, which flowed with suppressed turbulence. High particle concentration and fine-ash content hindered gas loss and maintained flow mobility. In order to explain the contemporaneous maintenance of high particle concentration, high sedimentation rate at the depositional boundary layer and a high level of mobility, it is also proposed that the flow(s) was continuously supplied at a high mass feeding rate. It is also proposed that internal gas pressure within the flow, directed downwards onto the substrate over which the flow was passing, reduced the friction between the flow and the substrate and also enhanced its mobility. The pervasive sub-horizontal fabric of aligned pumice, lithic and even biotite crystals indicates a consistent horizontal shear force existed during transport and deposition in the basal granular flow, consistent with the existence of a laminar, shearing, granular flow regime during the final stages of transport and deposition.     

Pyroclastic flows from the 1991 eruption of Unzen volcano,Japan

Pyroclastic flows from Unzen were generated by gravitational collapse of the growing lava dome. As soon as the parental lobe failed at the edge of the dome, spontaneous shattering of lava occurred and induced a gravity flow of blocks and finer debris. The flows had a overhanging, tongue-like head and cone- or rollershaped vortices expanding outward and upward. Most of the flows traveled from 1 to 3 km, but some flows reached more than 4 km, burning houses and killing people in the evacuated zone of Kita-kamikoba on the eastern foot of the volcano. The velocities of the flows ranged from 15 to 25 m/s on the gentle middle flank. Observations of the flows and their deposits suggest that they consisted of a dense basal avalanche and an overlying turbulent ash cloud. The basal avalanche swept down a topographic low and formed to tongue-like lobe having well-defined levees; it is presumed to have moved as a non-Newtonian fluid. The measured velocities and runout distances of the flows can be matched to a Bingham model for the basal avalanche by the addition of turbulent resistance. The rheologic model parameters for the 29 May flow are as follows: the density is 1300 kg/m³, the yield strength is 850 Pa, the viscosity is 90 Pa s, and the thickness of the avalanche is 2 m. The ash cloud is interpreted as a turbulent mixing layer above the basal avalanche. The buoyant portions of the cloud produced ash-fall deposits, whereas the dense portions moved as a surge separated from the parental avalanche. The ash-cloud surges formed a wide devastated zone covered by very thin debris. The initial velocities of the 3 June surges, when they detached from avalanches, are determined by the runout distance and the angle of the energy-line slope. A comparison between the estimated velocities of the 3 June avalanches and the surges indicates that the surges that extended steep slopes along the avalanche path, detached directly from the turbulent heads of the avalanches. The over-running surge that reached Kita-Kamikoba had an estimated velocity higher than that of the avalanche; this farther-travelled surge is presumed to have been generated by collapse of a rising ash-cloud plume.