Syn-eruptive , soft-sediment deformation of dilute pyroclastic density current deposits : triggers from granular shear , dynamic pore pressure , ballistic impacts and shock waves

Introduction Conclusions References

and occur after deposition but penecontemporaneous to the eruption.Isolated, cmscale, overturned beds with vortex forms have been interpreted to be the signature of shear instabilities occurring at the boundary of two granular media.They may represent the frozen record of granular, pseudo Kelvin-Helmholtz instabilities.Their recognition can be a diagnostic for flows with a granular basal boundary layer.The occurrence of degassing pipes together with basal intrusive dikes suggest fluidization during flow stages, and can facilitate the development of Kelvin-Helmholtz structures.The occurrence at the base of flow units of injection dikes in some outcrops compared with suction-driven local uplifts in others indicates the role of dynamic pore pressure.Variations of the latter are possibly related to local changes between depletive and accumulative dynamics of flows.Ballistic impacts can trigger unconventional sags producing local displacement or liquefaction.Based on the deformation depth, these can yield precise insights into depositional unit boundaries.Such impact structures may also be at the origin of some of the steep truncation planes visible at the base of the so-called "chute and pool" structures.Finally, the passage of shock waves emanating from the vent may be preserved in the form of trains of isolated, fine-grained overturned beds which may disturb the surface bedding without occurrence of a sedimentation phase in the vicinity of a vent.

Introduction
The dynamics of pyroclastic density currents (PDCs) remain poorly understood.This is despite the fact that they are one of the most efficient transport means on the flanks of volcanoes exhibiting explosive eruptions, thereby yielding a major risk potential for life, environment and infrastructures.Analogue and numerical modeling approaches are well-suited to investigate targeted hypothesized processes, but the question of which process to model can only be answered through real PDC data.Cross-bedded, dilute PDC deposits can contain intriguing overturned and deformed patterns attributed to soft-sediment deformation (SSD).The understanding of these structures can yield insight into the syn-and post-depositional processes surrounding the bed interface: i.e. the basal boundary layer (BBL), the bed state, and conditions extant in the emplacement environment.In particular, syn-depositional SSD structures provide constrain on the shearing and dynamic pore pressure at the BBL that controls the sedimentation of PDCs, whereas syn-eruptive SSD records information on the eruptive dynamics and depositional units.PDCs are largely emplaced subaerially under metastable conditions favoring SSD.Thus a variety of specific SSD triggers may occur during an eruption and PDC deposits represent excellent targets for studies of SSD.Introduction

Conclusions References
Tables Figures

Back Close
Full
A variety of triggers can be involved, predominantly related to seismogenic fluidization and/or liquefaction (special issue of Sedimentary Geology 235, Mohindra and Bagati, 1996;Owen, 1996), but also to tsunami waves (Alsop and Marco, 2012), storms (Chen and Lee, 2013) or volcanic base surges (Crowe and Fisher, 1973).

Conclusions References
Tables Figures

Back Close
Full tached from the layer to which they initially belong (Mills, 1983).Dikes are elongated and oriented quaquaversal to the layer they originate from.
Overturned laminae/beds: a few laminations or layers that show a coherent overturning recumbent towards the parent flow direction, laterally confined in otherwise undisturbed bedding.They can occur in sets of downstream repetitive but isolated patterns.They are distinguished from overturned stratification, which is an overturning of a stratal package as a whole (Allen and Banks, 1972;Røe and Hermansen, 2006;Bridge and Demicco, 2008, p. 357-358).
Vortex bedding: similar as overturned laminae, but with a vortex shape (Rowley et al., 2011)."Vorticity" is preferred to "rotation", since any simple shear deformation includes a rotational component.
The interpretation of the trigger mechanism(s) for SSD is not always straightforward and can include a combination of different effects.Here, distinction is made between the deformation, the agent of deformation, and the trigger.The deformation tensor in rock mechanics can be written as a sum of components of stretching, pure shear and simple shear (rotation).Identification of the relationships with the surroundings permits the interpretation of the physical agents responsible for the deformation as well as possible triggers.
Of interest here is the distinction between: (1) syn-sedimentary BBL (flow) shearing and dynamic pore pressure effects, (2) intra-deposit movements, and (3) post depositional mass movements.Bioturbation and biochemical effects are not dealt with here.
BBL shearing includes the effects of the flow drag during or directly after sedimentation.It can be enhanced by the sediment state and the nature of the BBL.Intra-deposit movements lead to sediment fabric rearrangement and deformation.These are often related to the expelling of trapped fluids during or after sedimentation, in situ releases, Introduction

Conclusions References
Tables Figures

Back Close
Full or compaction and loading.Mass movements are understood here as slumping, i.e. a short-scale, rather coherent sediment re-mobilization, the limit of which is taken to be debris flows.At the origin of the deformation can occur a trigger, a phenomenon that is not directly described in terms of the forces producing the deformation, but is causally responsible for their generation (e.g., ground-shaking facilitating fluidization of sediment, favoring fluid movements and producing pseudonodules).

PDCs and their possible SSD triggers
As particulate density currents, the depositional processes of PDCs are fundamental in their dynamics, since particles are both the agent of excess density driving momentum and the resulting sediment.Extreme and varied flow-substrate BBL processes may occur.The classification of Branney and Kokelaar (2002) emphasizes BBL processes and theorizes a classification into 4 types: granular-, fluid-escape-, falloutand traction-dominated BBLs.Douillet et al. (2014) discussed different types of crossbedding aggradation patterns as an upper or lower limit of the saltation threshold (the minimum shearing required to put grains in motion by wind), thus supposing a tractional BBL scheme.Alternatively, emplacement can be envisioned as a series of pulses with high basal concentration and no relationship to saltation, regardless of averaged concentration (Doronzo and Dellino, 2014), or stepwise en-masse deposition (Sulpizio and Dellino, 2008).The understanding of the nature and significance of BBL processes for PDCs may be further augmented by the study of syn-depositional SSD.

Conclusions References
Tables Figures

Back Close
Full sults in metastable deposits prone to further re-arrangement (Smith and Kokelaar, 2013).Moreover, the variations from very fine to very coarse beds typical of pyroclastic deposits as well as common inverse grading make them susceptible to SSD after deposition (Gibert et al., 2011).
In addition to their metastable nature, the eruptive environment itself is subject to a variety of triggers.Seismic activity associated with eruption further destabilizes freshly emplaced pyroclasts.Syn-PDC processes can be recorded in SSD (Crowe and Fisher, 1973), and the likely formation of traction carpets and granular BBL can produce granular shear instabilities (Rowley, 2010;Rowley et al., 2011;Smith and Kokelaar, 2013)."Flame-like" structures are often reported (McDonough et al., 1984;Valentine et al., 1989;Brand and White, 2007;Brand and Clarke, 2009) and when interpreted as sheared structures, can serve to reconstruct palaeoflow directions (Giannetti and Luongo, 1994;Brown et al., 2008).Fluid escape SSD (dikes, pipes, plumes, pillars) can occur by escape of water accompanying phreatomagmatic eruptions (Nocita, 1988), degassing of fresh pyroclasts (Gernon et al., 2008(Gernon et al., , 2009)), burning underlying vegetation, or be due to thermal expansion (Branney and Kokelaar, 2002, p. 61-66, and references therein).Interestingly, the high deposition rates combined with possible fluidized state of the flow can trap gases in the deposits that subsequently escape as degassing pipes within seconds after deposition (Komorowski et al., 2013).These can occur as fines-depleted pipes, few cms in length and diameter (Pistolesi et al., 2011;Smith and Kokelaar, 2013), or large dm to m scale depressions at the surface of deposits (Charbonnier and Gertisser, 2008).The high deposition rates also trigger simple load casts (Mattsson and Tripoli, 2011).Blocks ejected ballistically during an eruptive event deform the fresh deposits by landing (Gençalioğlu-Kuşcu et al., 2007;Jordan et al., 2013).Post eruptive processes are also common on steep sided volcanic edifices, with freshly deposited material likely to be unstable and slump (Fiske and Tobisch, 1978;Voight et al., 1983;Branney and Kokelaar, 1994;Ward and Day, 2006) as well as inherent contraction and compaction fractures following emplacement (Whelley et al., 2012).Introduction

Conclusions References
Tables Figures

Back Close
Full

Granular shear instabilities
Observations of syn-flow shear structures bring further insights into the BBL processes of PDCs.Simple shear is often invoked for the formation of overturned stratification (e.g., Allen and Banks, 1972;Mills, 1983;Røe and Hermansen, 2006).For such structures, the flow transmits and imposes part of its shear stress to the ground and thus translates the uppermost beds.In the other hand, sheat instabilities can be produced at the boundary between two fluids to form recurrent, vortex-shaped, Kelvin-Helmholtz instabilities.Valentine et al. (1989) suggested that flame-like SSD structures could be related to Kelvin-Helmholtz instabilities "between the bedload fluid and the overlying surge".Several analogue experimental studies with granular flows over grain beds have evidenced isolated but recurrent wave-like instabilities at the bed-flow interface (Goldfarb et al., 2002;Mangeney et al., 2010;Rowley, 2010;Roche et al., 2013;Farin et al., 2014).Goldfarb et al. (2002) have produced trains of wave instabilities with the shape of overturned laminae and noted that those were "likely produced by shearing differences" and "lacked any kind of vorticity".However, a rotational component must be present to produce the observed shark fin patterns.Rowley (2010) and Rowley et al. (2011) have imaged trains of shear-instabilities with well-developed vortex bedding convincingly interpreted as granular Kelvin-Helmholz instabilities.They further demonstrate the periodicity of these structures and document field examples.The wavy nature of those instabilities was further demonstrated in Farin et al. (2014), which also noted that the wavelength and amplitude are greatest for slopes close to the repose angle (highest speed).Roche et al. (2013) provided videos of the instabilities and an explanation for the fluid-like behavior of these instabilities.They suggested as a mechanism that negative dynamic pore pressures fluidize fine-grained beds and deform them as a whole rather than as individual grains.Other experimental work with granular flows has evidenced intriguing inter-penetration of beds over sinusoidal surfaces (Caicedo-Carvajal et al., 2006), longitudinal vortices in the flow direction (Forterre and Pouliquen, 2001), or Taylor vortices (Conway et al., 2004).Introduction

Conclusions References
Tables Figures

Back Close
Full 2 Geological settings and data The SSD structures presented here belong to different volcanoes and both magmatic and phreatomagmatic eruptions of various intensities and depositional environments.
As pointed by Mills (1983), SSD should be studied within their environment, and thus a brief context is introduced.Several types of SSD are identified with orders of magnitude between their dimensions as well as between the grain size of layers involved.Description of all discussed SSD structures is presented in Table 1.

Ubehebe crater (California, USA)
Ubehebe tuff ring is part of the Holocene/Pleistocene Ubehebe Craters complex and may have erupted between 0.8-2.1 ka (Sasnett et al., 2012).They erupted onto ancient lake sediments, at least partially phreatomagmatically.The arid climate does not explain the phreatomagmatic activity and interaction with a shallow water table is preferred (Sasnett et al., 2012).Crowe and Fisher (1973) reported SSD structures such as: contorted beds without preferred orientation, flame structures oriented with the flow direction and disrupted layers of thin tuff curled and pulled apart.They mapped the orientation of ballistic impact sags, mention post-eruption slumping on the Northwestern and Southeastern parts of the crater, and noted that SSD occurs within pre-existing channels filled with massive deposits but is absent in cross-bedded dominated overbanks.Here, a variety of SSD structures are documented from the Southern flank: folds-and-faults, curled layers, ballistic impact sags, a diapiric flame-like horizon and vortex features (Fig. 1).

Tungurahua (Ecuador)
The 17 August 2006 PDCs (Kelfoun et al., 2009;Hall et al., 2013;Bernard et al., 2014) are not linked to phreatomagmatic processes but rather to accumulation and subsequent destabilization of pyroclasts near the crater.The overbank sediments containing Introduction

Conclusions References
Tables Figures

Back Close
Full the SSD structures have been interpreted to have formed from dilute PDCs originating from dense PDCs by flow stripping (Douillet et al., 2013b).SSD was identified in a lacquer peel within well-developed mm-scale ash lamination (Fig. 2) located on the lee side (approx.20 cm from the crest) of an aggrading, transverse dune bedform that indicated very high sedimentation rates (Douillet et al., 2013a) approx.6 km from the vent.Two clusters of small-scale overturned and recumbent laminae occur at different height in the same horizons.The upper structure exhibits a single, well developed overturned laminaset (Fig. 2b), whereas the lower one is a cluster of several recumbent anticlines followed by relatively massive material with diffuse oversteepened bedding in the upstream direction (Fig. 2c and d).

Soufrière Hills (Montserrat)
The 11 February 2010 partial dome-collapse event of Soufrière Hills (Montserrat) produced a series of 6 block and ash flows, 5 of them occurring within 15 min and was the largest event since the 1995 awakening (Wadge et al., 2014;Stinton et al., 2014).Numerous degassing pipes were observed in block and ash flow deposits as well as massive ash units (Stinton et al., 2014).Other post-depositional structures are described by Stinton et al. (2014) as "rootless phreatic explosion craters", textiti.e.structures related to hot blocks turning water into steam explosively.They can have diameters between 1 and 30 m, consist of "decimetre-sized blocks in a coarse ash-rich matrix derived from the underlying primary PDC deposits" and have a contact to underlying cross-bedded units or down to the pre-collapse surface.Here, SSD structures are documented from the Belham river valley less than 6 km from the vent (Fig. 3).According to Stinton et al. (2014), only three PDCs flowed in this drainage (stage 3-H, 4-K and 4-6), Wadge et al. (2014) also mentioning PDCs in this zone for the 11 February 2010 collapse.Basal, small-scale dikes and pseudonodules intrude in underlying diffusely cross-stratified ash from a massive lapilli and ash lens, whereas fines-poor, small-scale pipes are found in the otherwise ash-rich, massive, overlying layer (Fig. 3a and b).The top of the latter has a contact with a series of 3 vortex and undulating forms (Fig. 3b).These deposits Introduction

Conclusions References
Tables Figures

Back Close
Full are found in the thalweg of the river valley, which may have contained some water.Another outcrop exhibits a large scale circular depression (ca. 3 m diam.) with ca. 10 cm deflation at the surface of the deposits (Fig. 3c).

Laacher See (Germany)
Laacher See was the location of a large eruption commonly attributed to phreatomagmatic explosions around 11 800 yr B.P. (Schmincke et al., 1973).Dune bedforms crossstratification made of coarse lapilli to fine ash intercalated with lapilli to volcanic dust fall horizons occur over tens of km 2 .Three isolated SSD structures are found around the "Wingertsbergwand" area, several km southward from the inferred vents (Fig. 4).
A composite SSD structure several m long and ca. 1 m thick occurs as a lateral series of tilted blocks that evolve into folds-and-faults beds in the (approximate) downstream direction, accommodating a local compression (Fig. 4a-f).It is abruptly confined in depth by the lower ash layer and underlying beds show no sign of deformation.A few tens of m distant, a structure of similar dimensions characterized by oversteepened lamination and downward oriented, massive, lapilli pseudonodules resembles a "chute and pool" structure (Fig. 4g-i).A further structure approx.150 m away, has a convex symmetrical form (ca. 10 cm vertical displacement) in an initially planar fine-grained bed.It is intercalated above a massive ash bed and below coarse-ash to lapilli, subplanar, diffuse bedsets (Fig. 4j and k).Flow direction inferred from overlying cross-beds is roughly oriented from left to right but may be sub-parallel to the outcrop wall.The bed is partly missing on the right from the deformation.Similar ash layers pinch out above the convex shape and may represent an overlap of the same unit.

Tower Hill (Victoria, Australia)
Tower Hill maar (ca.35 000 yr B.P., Sherwood et al., 2004;Prata and Cas, 2012) exhibits intriguing trains of oversteepened laminations contained within a single bedset (Fig. 5).They outcrop in the upper part of the southern rim (CRB quarry), parallel to the Introduction

Conclusions References
Tables Figures

Back Close
Full crater wall.Underlying beds fine up from massive coarse ash and lapilli by increasing occurrence of thin, sub-planar, ash beds forming a diffusely stratified lapilli-ash facies.This grades into the fine-grained ripple beds with topping SSD and the reverse sequence occurs above.This sequence suggests a fall phase progressively influenced by pseudo base-surge (in the sense of Waters and Fisher, 1971) with increasingly efficient fragmentation related to phreatomagmatic explosions at the fine-grained SSD bedsets (optimally efficient water : magma ratio in Prata, 2012).The flow direction inferred from the underlying ripple bedding is oriented roughly parallel to the lateral extension of the outcrop (Prata, 2012).The SSD consists of isolated, oversteepened laminations with coherent orientation.They are recurrent with wavelength of ca.50 cm and over hundreds of m.

Purrumbete Lake (Victoria, Australia)
The deposits forming the Purrumbete maar (ca.20 000 yr B.P.) are characterized by three temporally separated eruption phases and vent locations, with relatively dry as well as wet phreatomagmatic conditions (Jordan et al., 2013).Ballistic bombs with impact sags are widespread in these deposits, suggesting wet deposits (Jordan et al., 2013).The SSD documented here outcrops with two faces at right angles.Perpendicular to the crater, folds-and-faults structures increase in size, faulting and recumbence outward from the vent (Fig. 6a, b and e), but parallel to the rim, only chaotic flamelike structures are visible (Fig. 6c and d).The overlying deposits are planar laminated ash with individual laminae followed over several m.They lie conformably on the SSD horizon and are related to fallout.

Discussion and interpretation
As a general observation, many of the examples documented have a fine-grained underlying or basal layer (Purrumbete, Laacher See, Ubehebe, Merapi).Fine-grained

SED Introduction Conclusions References
Tables Figures

Back Close
Full layers are likely to have distinct rheological properties which yield the onset of deformation (Mills, 1983).Ash is no exemption, at both small (Gibert et al., 2011) and regional scales (Wiemer, 2014, and references therein).Fine deposits are likely to have low permeability, impacting either on the transfer and diffusion of dynamic pore pressure from subsequent flows (Roche et al., 2013), or on fluid escapes from enclosed layers (Peltier et al., 2012).

Slumps
For the Ubehebe thrusting folds-and-faults (Fig. 1a, top), observations point toward a gravitational slump: (1) deformed beds confined between undeformed strata, (2) large number of beds involved, (3) overturn orientation toward the thalweg of a channel, (4) axial planes of folds dipping upslope with folds' strike normal to microfaults, (5) imbrications and overlap (stacking) of deformed layers.The base of the folds-andfaults (Fig. 1a, base, and 1c) shows components of rotation, mixing and layer pull apart, indicating the floor thrust of slumping.Given the coherent state of the beds involved combined with their great variations in grain size and the thickness of SSD, a subaerial slump of wet sediment is favored.This would enhance cohesion and explain the gravity-sliding by overweight due to water.The interpretation is further supported by the characteristics of the nearby ballistic impact sags.Overlying strata are not deformed and so the slump has occurred during or between eruptive phases.

Granular Rayleigh-Taylor instabilities
The diapiric flame-like structures at Ubehebe (Fig. 1d) are remarkably similar to structures produced by granular Rayleigh-Taylor instabilities between a granular medium and air (Niebling et al., 2010).Usual interpretations of such features point toward liquefaction, loading, or water escape structures (e.g., Owen, 1996), a situation prone to Introduction

Conclusions References
Tables Figures

Back Close
Full occur during the eruption with high sedimentation rate and wet deposits.These explanations correspond to forms of granular Rayleigh-Taylor instabilities (Selker, 1993;Vinningland et al., 2007Vinningland et al., , 2010)).Some of the flame-like structures are overturned toward the flow direction, which may indicate shearing and syn-PDC SSD.
The curled and pulled apart, coarse-grained, isolated, flat "pseudonodules" from Ubehebe (Fig. 1b) are interpreted as detached load casts.These form in the presence of an inverse density gradient resulting from changes of porosity driven by the grain size distribution of successive layers (Mills, 1983;Bridge and Demicco, 2008, p. 353-354).As such, these also share the configuration for granular Rayleigh-Taylor instabilities.A shock (seismicity or impact) may trigger detachment, but is not necessary, and those structures may be post-eruptive.Their localized nature is taken to rule out remote triggers such as seismicity and no subsequent impact is visible above the structures.Further dynamic considerations coupled with the pseudo-wavelength of the structures and interface characteristics may resolve the question of their similarity with Rayleigh-Taylor instabilities (see Selker, 1993, and Appendix).

Granular shear and pseudo-Kelvin-Helmholz instabilities
At Tungurahua, the imbrication of overturned laminae with confinement within an otherwise undeformed bedset suggests syn-depositional processes (Fig. 2).SSD cannot be correlated with any impact sag.The orientation parallel to the flow direction suggests the influence of the latter.The vortex-shaped SSD structures are interpreted as granular shear instabilities related to Kelvin-Helmholtz vortices, based on reports and interpretations from analogue experiments (Rowley et al., 2011;Farin et al., 2014).If a pure wind BBL had moved the sediments, they would have begun to saltate as individual grains rather than deform as a whole (Douillet et al., 2014), and since the deposits were dry, no water can have triggered cohesion.Roche et al. (2013) explains the formation of wave instabilities at the interface between a fine-grained erodible bed and

SED Introduction Conclusions References
Tables Figures

Back Close
Full granular flow as linked with movements as a whole related to fluidization.This suggests that the observed features are indicative of a granular BBL and possible occurrence of traction carpets on the lee of the dune bedform.Although cross-stratification is generally interpreted as indicative of low particle concentration at the BBL, experiments by Leclair and Arnott (2005) have shown that laminations can be produced at more than 35 % particle concentration, a concentration at which a granular BBL can occur.The scale of the structures being similar to experimental results, the granular BBL is interpreted to be of the same order of thickness and velocities (few cm thick and few m s The downflow evolution of SSD at Tungurahua (Fig. 2c and d) brings further support to the discussion of Rowley et al. (2011).Indeed, they suggested that pseudo Kelvin-Helmholtz vortices may only be anecdotal in sedimentary records, since they intrinsically mix the deposits and create graded massive units.In the outcrop, well defined and thin lamination is visible downstream of the deformation front highlighted by vortices.In contrast, upstream from the front, stratification is comparatively thick and massive, with diffuse oversteepened laminations contained within the layers (Fig. 2d).This indicates the ploughing effect of the downstream migrating vortices that tend to mix and homogenize the initial bedding, as predicted by Rowley et al. (2011).
The recumbent and vortex structures at Ubehebe (Fig. 1d) have an overturning orientation with flow and a vortex shape.They only differ from Tungurahua by their occurrence in otherwise massive deposits.This may be an effect of successive ploughing by Kelvin-Helmholtz vortices or simply result from massive deposition.A vortex form is also observed at Soufrière Hills (Fig. 3b, top).In this case, the vortex is followed downstream by a gentle undulation and a steep step.Although the second and third structures have not a vortex shape, they are interpreted as proto, granular Kelvin-Helmholtz instabilities at different development stages, and the downstream repetition of deformation is taken as sign of the wavy nature of the instability.
Interestingly, sheared structure with a vortex-like structure are also present on the stoss and crest of dune bedforms covered by aggrading bedsets at Roccamonfina

Conclusions References
Tables Figures

Back Close
Full volcano (Italy, Fig. 5 in Giannetti and Luongo, 1994).If all these structures represent granular Kelvin-Helmholtz instabilities, they could share similar dynamics to their fluid analogue and quantitative information could be derived (Rowley et al., 2011, developed in Appendix).From theoretical considerations, BBL velocities of more than 2.5 m s −1 for 1 % relative particle concentration are necessary for instabilities to develop (Appendix, Fig. 8).This number rapidly drops for higher flow concentrations, and shear instabilities thus plausibly develop for basal granular BBL few cm in thickness.

Influence of dynamic pore pressure
SSDs from Soufrière Hills seem to originate from the dark mLA lensoidal layer that connects to the small basal dikes and intrusions and to the overlying pseudonodules and pipes (Fig. 3a and b).Komorowski et al. (2013) interpret small degassing pipes in the deposits of the Merapi 2010 block and ash flows as related to rapidly deposited and fluidized flows.Here, the dark mLA layer is interpreted as fluidized and overpressurized in dynamic pore pressure during flow in order to explain the basal dikes and intrusions as injection features.Basal dikes in subglacial deposits are indeed usually interpreted as indicating overpressure of the flows and injections (e.g., Douillet et al., 2012).The associated mLA layer would have held part of the overpressure through rapid sedimentation, and subsequently released the gas during deflation and compaction after burial by the overlying layer.This could further have destabilized the overlying beds and eased the formation of shear instabilities found at the upper interface of the mA layer above the pipes.Alternatively, the influence of bed water turned into steam cannot be ruled out in the river thalweg.The large scale depletion of the surface (Fig. The convex deformation of the fine-grained planar strata at Laacher See (Fig. 4j and  k) lacks any recumbent component, is isolated and no impact is visible.A localized decrease of the dynamic pore pressure of subsequent flows may have been transferred to the ground and slightly deformed the fine bed by suction.The deformed layer would have acted as a comparatively impermeable seal, containing the underpressure above it and lifting up.Negative dynamic pore pressures are indeed recorded from experimental granular flows and PDCs and trigger remobilization (Roche et al., 2013;Farin et al., 2014;Bernard et al., 2014).The dune bedform found just above the SSD may also have had some influence and produced a slight overweight.Aeolian dunes produce SSD on underlying beds, yet rather as load structures and in static settings (Chan and Bruhn, 2014).
The influence of dynamic pore pressure in destabilizing the sediment bed may be of importance in all syn-flow SSD scenarios.Both over-or under-pressurization can occur, and may connect with a flow's inhomogeneity.An overpressure is usually associated with flow front, and is directly followed by underpressure.Overpressurization may also be linked with accumulative flows, whereas underpressure relate to depletive evolution.

Impact records
At Laacher See, the tilted blocks (or domino/bookshelf structures) and thrust folds-andfaults packages locally share characteristics with slump folds (Fig. 4a-f).However, the very limited throw, absence of significant slope or possible slump trigger and the nested nature appear to exclude this interpretation.Given the great thickness of the disturbed beds, a surface instability, simple shearing or granular Kelvin-Helmholtz instability, also seems unlikely.dominated by tilted blocks indicates extension (Fig. 4a and d).
(2) The central part with the concave shape of the upper beds together with upward-penetrating flame-like beds suggest vertical compression (Fig. 4b and e).
(3) The folds-and-faults and local decollement in the downstream part record lateral displacement away from the central part ("escape zone", Fig. 4c and f).The source of the SSD can thus be localized above the central part, in the vertical compression zone, and have forced local displacement to the right.In light of this, the SSD is interpreted as the print of a large block bouncing on the bed and transmitting a deformation oriented with its trajectory.This is further supported by the presence of large blocks (> 3 m diam.) in nearby areas in deposits otherwise dominated by ash and lapilli.Noteworthy, the abrupt confinement of the deformation in depth indicates a higher state of compaction of the undeformed beds, and thus their belonging to an older event separated by sufficiently long time for compaction.The basal ash layer would represent an initial fall event belonging to the deformed unit.Thus impact sags may also be used to trace genetic units.
The diagnosis is easier at Ubehebe (Fig. 1e), where impacting blocks are nested in deformed beds and just above pseudonodules and dikes.Thorough observation indicates that the coarse and massive layer escaped into the enclosing fine-grained beds: it is the most disturbed and exhibits pseudonodules and small dikes with respect to both the over and underlying layers, which still contain stratification.The isotropic nature of the leakage with apparent absence of preferential escape directions supports a liquefaction mechanism.To account for the coarse-grained nature, water saturation is inferred, in agreement with the other Ubehebe SSD structures.A grain-flow triggered by an impact-induced liquefaction of the porous and water saturated coarse-ash enclosed in impermeable fine-grained layers has likely produced the nodules and dikes.
The pseudo "chute-and-pool" from Laacher See (Fig. 4g-i) shares similarities with both impact structures.The central part exhibits a depression with concave beds indicating compression.The right part is disturbed by massive lapilli material with downward-oriented pillows/pseudonodules (mL and PN in Fig. 4h and i).These are related to a liquefied grain flow of porous and water saturated lapilli beds.The pil-Introduction

Conclusions References
Tables Figures

Back Close
Full lows/pseudonodules are underlined on the right by a ca. 10 cm thick, massive, ashdominated contour with a diffuse front to the undisturbed cross-stratified bedsets (mA and dxstA in Fig. 4i).The massive fore front is interpreted as representing the final escape of water that was less coupled with sediments.The liquefaction event is related to a large block impact that could have bounced and compacted the concave central depression.The extreme right of the structure containing stoss stratification dipping at more than the repose angle (upper part of Fig. 4h) may have been oversteepened by the rearrangement of the underlying sediment, a process readily evidenced by Nocita (1988).The coarse lag breccia on top of the central depression may either indicate that the impacting block stayed in place and acted upon the depositional dynamics, resuspended fines during impact, or be a simple infill of the topography.

A trigger for "chute and pool" structures
The two impact SSDs from Laacher See share remarkable similarities with the basal oversteepened truncations observed in structures generally interpreted as "chute and pools" (types I to IV of Schmincke et al., 1973).If the disturbed beds had been slightly more destabilized and permitted entrainment, the same configuration would be observed.Such impact SSDs would explain the oversteepened truncations and be at the origin of some of the "chute and pools" structures (see also Nocita, 1988).This would also explain the observation by Schmincke et al. (1973) that "chute and pools" occur in rather proximal parts, since ballistic blocks are likely to land closer to the crater than the total distance travelled by a PDC.This interpretation does not contradict the subsequent hydraulic jump dynamics of the structures, but the jump would be a consequence of the bed morphology rather than the other way round as usually suggested.
A hydraulic jump would however not be necessary and simple morphological blocking of the bedload equally well explains the "chute and pools" depositional patterns (basal blocking and stop-and-go models in Douillet et al., 2013a;Martínez et al., 2007, resp.).
The answer likely lies upstream from these structures, at the proximal truncation limit.

SED Introduction Conclusions References
Tables Figures

Back Close
Full

Deformation driven by shock waves
At Tower Hill (Fig. 5), the regularity of patterns, high degree of preservation and absence of slope appear to discredit slumping.The lateral persistence indicates a largescale effect and discredit shear instabilities.Indeed, a flow with thin granular BBL forming pseudo Kelvin-Helmholtz instabilities is unlikely to stay in this state over several hundred m.Moreover, either lateral flow velocities were slow enough for fine ash and volcanic dust with ripple lamination to deposit, or the ground was covered with a stretch of water.During the phreatomagmatic phase with efficient fragmentation associated with the fine beds, shock waves may have been produced by the explosions (e.g., Scolamacchia and Schouwenaars, 2009).These could propagate close to the rim, quaquaversal to the southern vent, and destabilize the fine-grained bedsets by transmitting their orientation to the ground.Valentine et al. (1989) suggested shock waves as a possible trigger for overturned flame-like structures.They noted that "when a shock passes over a granular deposit, bed particles experience a lift force due to the change in velocity across the shock" and "the bed immediately behind the shock has been observed in experiments to take on a wavelike configuration" citing the convincing experiments by Borisov et al. (1967).Recent shock experiments by Wayne et al. (2013) developed recumbent vortex-like shapes on dust beds and further support the interpretation (see also Fedorov, 2004).The Purrumbete structure (Fig. 6) has a preferential direction away from crater: (1) all beds are overturned outside of the crater, (2) the deformation, vorticity degree and thickness of beds involved increase away from crater, whereas (3) the crater-parallel face is chaotic.Microfaults suggest cohesion, and there is neither evidence of traction nor of granular flow in the overlying planar deposits related to fallout, thus granular shear is excluded.An envisaged interpretation is that these beds are involved in a small-scale slump.However, overlying beds lie conformably on top of the deformed strata, and are thus emplaced after deformation, implying that a very small amount of material would have slumped, unlikely to be sufficient to yield a consequent gravita-

SED Introduction Conclusions References
Tables Figures

Back Close
Full tional potential for deformation.Moreover, the bed is only gently sloping (< 5 • ).The overturn, vorticity direction as well as the imbrication fabric at thrust faults would suggest an outward oriented slump, but no scar is visible at the deformation onset (left part).Rather, the evolution of intensity of deformation, absence of scar in the proximal side, and imbrication suggest that deformation could origin from the distal part, with some force pushing the sediment toward the crater.Similarly as for the Tower Hill maar, these structures can be the record of shock waves that destabilized the uppermost deposits and transferred them a tilt.Here again, the vortex-like entrainment evidenced in shock experiments supports the interpretation (Borisov et al., 1967;Wayne et al., 2013), and the proximity from the vent make shock wave influence probable.The passage of a shock wave is likely asedimentary (not associated with deposition), and the conformity of overlying fall beds supports a trigger mechanism without sedimentation, the signature being uniquely present as deformation.As SSD triggered by shock waves, these can share similarities with Richtmyer-Meshkov instabilities (Brouillette, 2002).

Conclusions
The exercise presented here has demonstrated the richness of information contained in SSD structures from the dilute PDC environment.SSD contains a record of syn-flow and syn-eruptive processes combined with post-depositional bed-state information.Syn-flow processes were evidenced through granular, pseudo Kelvin-Helmholtz instabilities as well as evidences of suction and injection related to dynamic pore pressure of the flows.These observations feed the understanding of BBL processes of PDCs.Basal intrusions support the interpretation of fluidized flows with dynamic pore overpressure.Vortex-shaped laminae may be a valid indicator of granular-based flows or traction carpets.The suction vs. injection at the base of flows can relate to depletive resp.accumulative phases of a flow.

SED Introduction Conclusions References
Tables Figures

Back Close
Full Several syn-eruptive processes are recorded by SSD.Ballistic impacts may take more evolved forms than simple sags.They yield information on the bed state such as the compaction degree and water content, which helps to delimit eruptive units and environmental context.Some forms of impact may be at the origin of the so-called "chute and pools" bedforms.Slumps similarly inform on water content and eruptive units.The understanding of prevailing eruption type (wet vs. dry eruptions) may thus benefit from thorough analysis of SSD.Finally, we suggest that shock waves may leave a signature in the sediments by destabilization and overturning of the surface beds close to the vent without any direct deposits.
SSD from PDCs are of interest in the context of sedimentary research since they record subaerial, syn-and post-flow SSD structures, emphasizing that water is not a prerequisite for SSD.Moreover, PDC deposits can be unstable and have large permeability contrasts that facilitate SSD formation.Finally, the recognition of structures similar to instabilities occurring at fluid boundaries (Kelvin-Helmholz, Rayleigh-Taylor) further emphasizes the similarities between fluids and granular mixtures.SSD seems widespread in deposits of dilute PDCs, especially from phreatomagmatic eruptions, and should be addressed more attention.The variety of possible triggers, especially in the context of explosive volcanic eruptions, calls for further field and experimental work.

Appendix A: Instabilities between two fluids A1 Granular Kelvin-Helmholtz instabilities at a bed-flow interface
Given structures interpreted as granular Kelvin-Helmholtz instabilities, a theoretical resolution similar to the fluid instability can be expressed (Rowley et al., 2011).Any fluid dynamics analysis is based on the integration of "infinitesimal fluid elements", a notion comparable to grains in a granular mixture.The fluid-dynamics analytical method just justify in itself its applicability to granular mediums.

SED Introduction Conclusions References
Tables Figures

Back Close
Full The problem is taken in 2-D with reference frame (e x -flow parallel direction-, e z -upward direction parallel to g, the gravity acceleration-).Consider two homogenous mediums F 1 and F 2 , F 2 lying above F 1 and the interface an infinite horizontal plane.Suppose the fluids of densities ρ 1 and ρ 2 , incompressible (Dρ 1,2 /Dt = 0), inviscid (ν = 0), with constant horizontal velocity u 1,2 (z) = U 1,2 e x and irrotational.The surface disturbance (ξ) can be written in the form (see Drazin, 2002;Douillet, 2014, Chap.2): with k the wave number.Linearization of the problem posed by the boundary conditions has solution (see Drazin, 2002;Douillet, 2014, Chap. II.2): Assumptions can be made for the case of an instability between a granular flow and deposit.The deposit does not move (U 1 = 0) and the flow density is a portion of the deposit density (ρ 2 = xρ 1 with 0 ≤ x ≤ 1).Thus Eq. (A2) simplifies into: In order that a wave occurs, Eq. (A3) must have an imaginary component (the angular velocity w = Im(s)).The second term in s must be real for an exponential decay or increase to develop, and thus, be an instability.Thus the term under the square root must be positive and a condition for a bed-flow instability is (see also Rowley et al., 2011):

Conclusions References
Tables Figures

Back Close
Full This condition is granted for large wavenumber (k), i.e. short waves, high particle concentrations (x), or large flow velocities (Fig. 8).Further, the phase velocity of an instability (c = w/k) can be derived: Under the assumptions, the wave velocity is thus entirely characterized by the concentration difference between the bed and flow (x) and the velocity of the latter (U 2 ), and the wavelength of the instability (λ = 2π/k) does not appear explicitly.

A2 Granular Rayleigh-Taylor instabilities
A Rayleigh-Taylor instability is a surface instability between two resting fluids of different densities.Thus Eq. ( A2) can be equally used with U 1,2 = 0.For the case of the curled and pulled appart structures at Ubehebe (Fig. 1b), the upper coarse grained layer was sinking in the massive fine-grained layer underneath, thus ρ 2 = xρ 1 with x ≥ 1, and Eq.(A2) simplifies into: The field observation is the lengthscale of the curled layers (λ = 2π/k = ca.15-30 cm).
The missing variables are a timescale for the growth of the instability and the density ratio.Estimating one permits to quantify the other.Introduction

Conclusions References
Tables Figures

Conclusions References
Tables Figures

Back Close
Full Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | 3c) may relate to similar deflation of fluidized pockets, although simple re-arrangement of the grains underneath or any depletion could lead to similar surface expressions.The surface mainly consists of coarse particles and small deflation cracks developed, thus the structure may relate to fines' elutriation.A relation with the "phreatic rootless explosions" in Stinton et al. (2014) is not favored.Discussion Paper | Discussion Paper | Discussion Paper | Yet, there is a preferential orientation directed roughly outward from the postulated crater.The lateral evolution of the SSD leads another insight.(1) The upstream part Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Brouillette, M.: The Richtmyer-Meshkov instability, Annu.Rev.
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Figure 1 .Figure 2 .Figure 3 .Figure 4 .Figure 6 .Figure 7 .Figure 8 .
Figure 1.SSD structures from Ubehebe craters.(a) folds-and-faults structure (F&F) related to slump, with interpretation of the outlined lower part in (c).(b) Curled and pulled apart coarsegrained layers interpreted as detached pseudonodules.(d) Diapiric flame-like structures in upper part, recumbent overturned bed in the middle right, vortex beds in the lower part.(e) Interpenetrating coarse bed with pseudonodules (PN) and dikes at the base of ballistic impact sags.

Table 1 .
Main characteristics of SSD structures presented in figures and discussed in text.