numerical modeling of the emplacement of socompa rock avalanche chile 1 the 7 5 ka socompa sector collapse emplaced 25 km3 of fragmented rock as a thin but widespread 500 km2 avalanche deposit followed by late stage sliding of 11 km3 as toreva blocks most of the avalanche mass was emplaced dry although saturation of a basal shear layer cannot be excluded modeling was carried out using the depth averaged granular flow equations in order to provide information on the flow behavior of this well preserved long run out avalanche results were constrained using structures preserved on the surface of the deposit as well as by deposit outline and run up a proxy for velocity models assuming constant dynamic friction fail to produce realistic results because the low basal friction angles 1 to 3 5 necessary to generate observed run out permit neither adequate deposition on slopes nor preservation of significant morphology on the deposit surface a reasonable fit is obtained however if the avalanche is assumed simply to experience a constant retarding stress of 50 100 kpa during flow this permits long run out as well as deposition on slopes and preservation of realistic depositional morphology in particular the model explains a prominent topographic escarpment on the deposit surface as the frozen front of a huge wave of debris reflected off surrounding hills the result that socompa avalanche experienced a small approximately constant retarding stress during emplacement is consistent with a previously published analysis of avalanche data

formation of recurring slope lineae on mars by rarefied gas triggered granular flows active dark flows known as recurring slope lineae have been observed on the warmest slopes of equatorial mars the morphology composition and seasonality of the lineae suggest a role of liquid water in their formation however internal and atmospheric sources of water appear to be insufficient to sustain the observed slope activity experimental evidence suggests that under the low atmospheric pressure at the surface of mars gas can flow upwards through porous martian soil due to thermal creep under surface regions heated by the sun and disturb small particles here we present numerical simulations to demonstrate that such a dry process involving the pumping of rarefied gas in the martian soil due to temperature contrasts can explain the formation of the recurring slope lineae in our simulations solar irradiation followed by shadow significantly reduces the angle of repose due to the resulting temporary temperature gradients over shaded terrain and leads to flow at intermediate slope angles the simulated flow locations are consistent with observed recurring slope lineae that initiate in rough and bouldered terrains with local shadows over the soil we suggest that this dry avalanche process can explain the formation of the recurring slope lineae on mars without requiring liquid water or co2 frost activity transient streaks on martian slopes have been attributed to liquid water simulations show that a dry avalanche process involving the flow of gas in the martian soil due to temperature contrasts can instead explain these recurring features

effects of near fault seismic loadings on run out of large scale landslide a case study this study presents the run out analysis of the daguangbao landslides subjected to near fault multi direction earthquake forces using discontinuous deformation analysis dda the daguangbao landslide is the largest landslide induced by the 2008 wenchuan earthquake in order to investigate the effects of near fault seismic force on landslide run out kinematic behavior of sliding mass is simulated by a dynamic discrete numerical analysis method called dda in this simulation based on the shape of failure surface and the feature of slope geology the whole slope is divided into three parts base block upper sliding mass and lower sliding mass then two sliding masses are divided into the smaller discrete deformable blocks based on pre existing discontinuities size effect of the huge landslide is also considered baseline corrected real horizontal and vertical ground motions are taken as volume force acting to the base block the results show that seismic force has a significant influence on the landslide progression sliding distance and shape of post failure results of the horizontal and vertical situation are in good agreement with those obtained from post earthquake field investigation remote sensing image and description from the survivors

interpreting rockfall activity on an outcrop talus slope system in the southern japanese alps using an integrated survey approach abstract understanding relationships between meteorological conditions and rockfall activity is important for the management of rockfall risks as well as for the estimation of future rockfall activity under climate change seasonal activity of rockfall on talus slopes depends on the actual triggers of rockfall which in turn may vary as well between outcrops because of different rockfall sources and the influence of topography in this study we interpret spatial and temporal variations in rockfall activity and related sediment transfer processes in an outcrop talus slope system in central japan through monitoring using rock temperature and moisture sensors sediment traps and time lapse cameras in addition tree ring analysis has been added to understand possible changes of rockfall activity over longer time periods results suggest that the frequency of grouped gravel transport i e soil creep was higher than that of individual transport i e rockfall and dry ravel on the talus slope whereas rockfall likely is the predominant process from the outcrops observations from sediment traps showed that coarse gravel is characterized by higher bounce heights and that it travelled longer distances than fine gravel rockfall from the outcrops was most active during dormancy of trees which locally lasts from autumn to early spring and can be related to freeze thaw cycles in contrast rockfall and other sediment transfer processes were monitored on the talus slope all year round secondary movement of gravels on the talus slope such as soil creep during rainfall events likely triggers rockfall during times when the outcrops are not producing rockfall i e from late spring to summer consequently our results suggest that the spatio temporal rockfall activity varies substantially between outcrops and talus slopes and that differences in process behavior should be taken into account in future disaster mitigation measures

shedding dynamics and mass exchange by dry granular waves flowing over erodible beds abstract a continuous exchange of particles between an erodible substrate and the granular flow above it occurs during almost all geophysical events involving granular material such as snow avalanches debris flows and pyroclastic flows the balance between eroded and deposited material can drastically influence the runout distance and duration of the flow in certain conditions a perfect balance between erosion and deposition may occur leading to the steady propagation of material in which the flow maintains its shape and velocity throughout it is shown experimentally how the erosion deposition process in dense flows of sand 160 200 m on an erodible bed of the same material produces steadily propagating avalanches that deposit subtle levees at their lateral extent moreover it is shown in this paper by using two colours of the same sand that although the avalanche is propagating at constant velocity and maintaining a constant shape the grains that are initially released are deposited along the flow path and that the avalanche will eventually be composed entirely of particles that are eroded from the bed different steady travelling wave regimes are obtained depending on the slope angle thickness of the erodible layer and the amount of material released outside of the range of parameters where steady travelling waves form the avalanches loose mass and decay if the initial amount of material released is too small or if the initial release is too large they re adjust to a steadily propagating regime by shedding material and breaking into smaller avalanches at its rear side numerical simulations are performed using a shallow water like avalanche model together with a friction law that captures the erosion deposition process in flowing to static regimes and a transport equation for the interface between layers of the two colours the characteristic behaviours observed in the experiments are qualitatively reproduced specifically the complex processes such as the exchange of particles leading to a change in colour of the avalanche and the formation of lateral levees are captured by the model finally a comparison is made with deposits in lunar craters which are interpreted as closely analogous to the deposits formed in our laboratory experiments

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 m3 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

flow of variably fluidized granular masses across three dimensional terrain 2 numerical predictions and experimental tests numerical solutions of the equations describing flow of variably fluidized coulomb mixtures predict key features of dry granular avalanches and water saturated debris flows measured in physical experiments these features include time dependent speeds depths and widths of flows as well as the geometry of resulting deposits three dimensional 3 d boundary surfaces strongly influence flow dynamics because transverse shearing and cross stream momentum transport occur where topography obstructs or redirects motion consequent energy dissipation can cause local deceleration and deposition even on steep slopes velocities of surge fronts and other discontinuities that develop as flows cross 3 d terrain are predicted accurately by using a riemann solution algorithm the algorithm employs a gravity wave speed that accounts for different intensities of lateral stress transfer in regions of extending and compressing flow and in regions with different degrees of fluidization field observations and experiments indicate that flows in which fluid plays a significant role typically have high friction margins with weaker interiors partly fluidized by pore pressure interaction of the strong perimeter and weak interior produces relatively steep sided flat topped deposits to simulate these effects we compute pore pressure distributions using an advection diffusion model with enhanced diffusivity near flow margins although challenges remain in evaluating pore pressure distributions in diverse geophysical flows riemann solutions of the depth averaged 3 d coulomb mixture equations provide a powerful tool for interpreting and predicting flow behavior they provide a means of modeling debris flows rock avalanches pyroclastic flows and related phenomena without invoking and calibrating rheological parameters that have questionable physical significance

a model for the analysis of rapid landslide motion across three dimensional terrain a new numerical model for the dynamic analysis of rapid flow slides debris flows and avalanches has been developed the model is an extension of an earlier algorithm and is implemented using a nu

forensic investigations of the cima salti landslide northern italy using runout simulations abstract over the last decades a movement has begun to reclassify deposits previously misidentified as having origins other than landsliding the reverse problem of incorrectly assigning deposits a mass movement origin however has been addressed less in the landslide community the cima salti landslide in the lake garda region of northern italy is a cautionary tale of assuming source areas and volumes it was traditionally thought to have dammed tenno lake in the middle ages and to have a volume of 20 30 mm3 we show through geological field data and simple runout simulations with the codes dan3d and shaltop which produced comparable results that the volume of the landslide was likely significantly overestimated in the past and that it most likely did not dam tenno lake as has been assumed we propose that a smaller volume landslide 2 5 mm3 was deposited on stagnant ice melting in situ in the lateglacial period a relatively minor event in the complex history of the magnone valley this interpretation emphasizes the importance of careful field investigations and assumption validation at cima salti and in a broader context it also shows the unique capacity of landslide simulation to guide field observation and discriminate mass emplacement processes

numerical modeling of self channeling granular flows and of their levee channel deposits when not laterally confined in valleys pyroclastic flows create their own channel along the slope by selecting a given flowing width furthermore the lobe shaped deposits display a very specific morphology with high parallel lateral levees a numerical model based on saint venant equations and the empirical variable friction coefficient proposed by pouliquen and forterre 2002 is used to simulate unconfined granular flow over an inclined plane with a constant supply numerical simulations successfully reproduce the self channeling of the granular lobe and the levee channel morphology in the deposits without having to take into account mixture concepts or polydispersity numerical simulations suggest that the quasi static shoulders bordering the flow are created behind the front of the granular material by the rotation of the velocity field due to the balance between gravity the two dimensional pressure gradient and friction for a simplified hydrostatic model competition between the decreasing friction coefficient and increasing surface gradient as the thickness decreases seems to play a key role in the dynamics of unconfined flows the description of the other disregarded components of the stress tensor would be expected to change the balance of forces the front s shape appears to be constant during propagation the width of the flowing channel and the velocity of the material within it are almost steady and uniform numerical results suggest that measurement of the width and thickness of the central channel morphology in deposits in the field provides an estimate of the velocity and thickness during emplacement citation mangeney a f bouchut n thomas j p vilotte and m o bristeau 2007 numerical modeling of self channeling granular flows and of their levee channel deposits

dynamics of pyroclastic density currents conditions that promote substrate erosion and self channelization mount st helens washington usa abstract the may 18th 1980 eruption of mount st helens msh produced multiple pyroclastic density currents pdcs burying the area north of the volcano under 10s of meters of deposits detailed measurements of recently exposed strata from these pdcs provide substantial insight into the dynamics of concentrated currents including inferences on particle particle interactions current mobility due to sedimentation fluidization and internal pore pressure particle support mechanisms the influence of surface roughness and the conditions that promote substrate erosion and self channelization four primary flow units are identified along the extensive drainage system north of the volcano each flow unit has intricate vertical and lateral facies changes and complex cross cutting relationships away from source each flow unit is an accumulation from an unsteady but locally sustained pdc or an amalgamation of several pdc pulses the pdcs associated with units i and ii likely occurred during the pre climactic waxing phase of the eruption these currents flowed around and filled in the hummocky topography leaving the massive to diffusely stratified deposits of units i and ii the deposits of both units i and ii are generally more massive in low lying areas and more stratified in areas of high surface roughness suggesting that surface roughness enhanced basal shear stress within the flow boundary units iii and iv are associated with the climactic phase of the eruption which produced the most voluminous and wide spread pdcs both flow units are characteristically massive and enriched in vent derived lithic blocks these currents flowed over and around the debris avalanche deposits as evidenced by the erosion of blocks from the hummocks unit iii is massive poorly sorted and shows little to no evidence of elutriation or segregation of lithics and pumice suggesting a highly concentrated current where size density segregation was suppressed unit iv shows similar depositional features but typically has a basal lithic rich region is variably fines depleted and contains pumice lobes suggesting density segregation in a less concentrated current relative to unit iii deep erosive channels cut by the unit iii current and thick lithic levee deposits within unit iv occur in an area where debris avalanche relief is limited suggesting self channelization developed as a function of internal flow dynamics an increase in the proportion and size of lithic blocks is found 1 downstream of debris avalanche hummocks suggesting the pdcs were energetic enough to locally entrain accidental lithics from the hummocks and transport them tens of meters downstream and 2 within large channels cut by later pdcs into earlier pdc deposits suggesting self channelization of the flows increased the carrying capacity of the subsequent channelized currents finally the combination of thick massive deposits with a high percentage of fine ash within unit iii and in the medial distal depositional regions of units ii iv suggests that the pdcs developed and maintained a high internal pore pressure during transport and deposition the most important results include our ability to understand the role of internal pore pressure on current mobility the influence of self channelization on carrying capacity of the currents and the influence of surface roughness on substrate erosion these observations have critical consequences for understanding the flow dynamics and hazard potential of pdcs

formation of levees and en echelon shear planes during snow avalanche run out snow avalanches often form levees and en echelon shear planes in the run out zone we describe the formation of these depositional structures using a simple model that accounts for the role of granular fluctuations in avalanche motion a mathematical feature of this model is the existence of a bifurcation saddle point describing how granular fluctuations control the avalanche velocity in the run out zone the saddle point discriminates between a flowing and stopping regime and defines the physical boundary between the flow and non flow regions of the avalanche i e the location of shear planes in the avalanche deposits the formation of a shear plane depends on the interplay between terrain slope and avalanche mass flux which varies from avalanche head to tail levees can form immediately at the avalanche front or for steep slopes and low mass fluxes at the avalanche tail at ravine and gully shoulders the mass flux is restricted thus initiating levee formation we find that the levee lines are parallel to the flow direction when the mass flux is constant en echelon shear lines occur when the mass flux is decreasing we test the model using several case studies where we have accurate laser scans of avalanche deposits our results suggest that avalanche flow parameters can be determined from simple levee measurements or conversely formation of levees and flow fingers can be predicted once the parameters governing the granular fluctuations are known

an efficient and stable hydrodynamic model with novel source term discretization schemes for overland flow and flood simulations numerical models solving the full 2 d shallow water equations swes have been increasingly used to simulate overland flows and better understand the transient flow dynamics of flash floods in a catchment however there still exist key challenges that have not yet been resolved for the development of fully dynamic overland flow models related to 1 the difficulty of maintaining numerical stability and accuracy in the limit of disappearing water depth and 2 inaccurate estimation of velocities and discharges on slopes as a result of strong nonlinearity of friction terms this paper aims to tackle these key research challenges and present a new numerical scheme for accurately and efficiently modeling large scale transient overland flows over complex terrains the proposed scheme features a novel surface reconstruction method srm to correctly compute slope source terms and maintain numerical stability at small water depth and a new implicit discretization method to handle the highly nonlinear friction terms the resulting shallow water overland flow model is first validated against analytical and experimental test cases and then applied to simulate a hypothetic rainfall event in the 42 km2haltwhistle burn uk

flow of variably fluidized granular masses across three dimensional terrain 1 coulomb mixture theory rock avalanches debris flows and related phenomena consist of grain fluid mixtures that move across three dimensional terrain in all these phenomena the same basic forces govern motion but differing mixture compositions initial conditions and boundary conditions yield varied dynamics and deposits to predict motion of diverse grain fluid masses from initiation to deposition we develop a depth averaged three dimensional mathematical model that accounts explicitly for solid and fluid phase forces and interactions model input consists of initial conditions path topography basal and internal friction angles of solid grains viscosity of pore fluid mixture density and a mixture diffusivity that controls pore pressure dissipation because these properties are constrained by independent measurements the model requires little or no calibration and yields readily testable predictions in the limit of vanishing coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods and in the limit of vanishing fluid stress they describe one phase granular avalanches analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations which can simulate interaction of high friction surge fronts with more fluid debris that follows special numerical methods described in the companion paper are necessary to solve the full equations but exact analytical solutions of simplified equations provide critical insight an analytical solution for translational motion of a coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically a solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris solutions for static equilibrium of source areas describe conditions of incipient slope instability and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits simplified equations and solutions may apply in additional situations identified by a scaling analysis assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full scale flows in which fluid effects are significant therefore large geophysical flows can exhibit dynamics not evident at laboratory scales

martian gullies a comprehensive review of observations mechanisms and insights from earth analogues upon their discovery in 2000 martian gullies were hailed as the first proof of recent i e less than a few million years flowing liquid water on the surface of a dry desert planet many processes have been proposed to have formed martian gullies ranging from liquid water seepage from aquifers melting of snow ice and frost to dry granular flows potentially lubricated by co2 terrestrial analogues have played a pivotal role in the conception and validation of gully formation mechanisms comparison with the terrestrial landscape argues for gully formation by liquid water debris flows originating from surface melting however limited knowledge of sediment transport by sublimation is a critical factor in impeding progress on the co2 sublimation hypothesis we propose avenues towards resolving the debate a laboratory simulations targeting variables that can be measured from orbit b applications of landscape evolution models c incorporation of the concept of sediment connectivity d using 3d fluid dynamic models to link deposit morphology and flow rheology and e a more intense exchange of techniques between terrestrial and planetary geomorphology including quantitative and temporal approaches finally we emphasize that the present may not accurately represent the past and that martian gullies likely formed by a combination of processes

sinuous gullies on mars frequency distribution and implications for flow properties recent gullies on mars are suspected to be the result of liquid water bearing flows a formation from wet flows has been challenged by studies invoking granular dry flows our study focuses on the sinuous shapes observed for some of the recent martian gullies sinuous gullies are found in locations and slopes of 10 15 similar to straight gullies and they are therefore related to the same formation processes numerical simulations of granular flows are performed here by introducing topographic variations such as obstacles roughness or slope changes that could possibly generate flow sinuosity none of these simulations was able to reproduce sinuous shapes on a slope lower than 18 with friction angles typical of dry granular material the only way to simulate sinuous shapes is to create small amplitude periodic variations of the topography of the deposit an origin not supported by current martian imagery given the presence of sinuosity in natural terrestrial debris flows we have concluded that sinuous martian gullies are better reproduced by liquid water bearing debris flows sinuous shapes in leveed flows are used to derive mechanical parameters from several martian gullies using photoclinometry values in yield strength of 100 2200 pa velocities of 1 1 3 3 m s 1 and viscosities from 40 to 1040 pa s are found which are all within the range of values for terrestrial debris flows with various proportions of liquid water 20 40

mobility and topographic effects for large valles marineris landslides on mars 1 recent experiments on dry granular flows over horizontal plane bare some similarities with large martian landslides observed in valles marineris vm however martian normalized runout are twice as large as those that observed in dry granular flow experiments numerical simulations on theoretical 2d and real 3d topographies reconstructed from remote sensing data show that slope effects significantly reduce the shift between experimental results and martian observation however topography effects are not strong enough to explain the high mobility of martian landslides as a result other physical and or geological processes should play a key role into the dynamics of martian landslides a new mobility is defined that makes it possible to characterize the dynamics of the flow regardless of the geometry of the released mass and of the underlying topography

erosion and mobility in granular collapse over sloping beds 1 we describe laboratory experiments of granular material flowing over an inclined plane covered by an erodible bed designed to mimic erosion processes of natural flows travelling over deposits built up by earlier events two controlling parameters are the inclination of the plane and the thickness of the erodible layer we show that erosion processes can increase the flow mobility i e runout over slopes with inclination close to the repose angle of the grains r by up to 40 even for very thin erodible beds erosion efficiency is shown to strongly depend on the slope of the topography entrainment begins to affect the flow at inclination angles exceeding a critical angle c r 2 runout distance increases almost linearly as a function of the thickness of the erodible bed suggesting that erosion is mainly supply dependent two regimes are observed during granular collapse a first spreading phase with high velocity followed by a slow thin flow provided either the slope or the thickness of the erodible bed is high enough surprisingly erosion affects the flow mostly during the deceleration phase and the slow regime the avalanche excavates the erodible layer immediately at the flow front waves are observed behind the front that help to remove grains from the erodible bed steep frontal surges are seen at high inclination angles over both rigid or erodible bed finally simple scaling laws are proposed making it possible to obtain a first estimate of the deposit and emplacement time of a granular collapse over a rigid or erodible inclined bed

pressures at the base of dry flows of angular rock fragments as a function of grain size and flow volume experimental results geophysical granular flows such as pyroclastic flows and rock avalanches kill people and damage properties worldwide the pressures exerted at their base affect the retarding forces that act on them and for this reason affect also their mobility that is important to foresee when assessing natural hazards in mountain regions here we present the results of experiments obtained by measuring with a load cell the basal pressures exerted by dry and cohesionless granular flows that descend a curved chute in the laboratory the interaction between these flows and the chute surface on which they travel is dominated by collisions of particles and or clusters of particles a dimensional analysis suggests that the energy dissipation of these flows increases as grain size increases and as flow volume decreases all the other features equal therefore the smaller the grain size and the larger the volume the larger is expected to be flow mobility although the longer travel distances of the centre of mass of finer grain size flows are easily discernible in our experiments the effect of volume is probably hidden by additional phenomena such as the deposition first of the frontal portion of longer flows on the less steep more distal part of the slope that prevents the rear portion and the centre of mass of the flows to travel further downhill

a new depth averaged model for flow like landslides over complex terrains with curvatures and steep slopes abstract flow like landslides are one of the most catastrophic types of natural hazards due to their high velocity and long travel distance they travel like fluid after initiation and mainly fall into the flow movement type in the updated varnes classification hungr et al 2014 in recent years depth averaged models have been widely reported to predict the velocity and run out distance of flow like landslides however most of the existing depth averaged models present different shortcomings for application to real world simulations this paper presents a novel depth averaged model featured with a set of new governing equations derived in a mathematically rigorous way based on the shallow flow assumption and mohr coulomb rheology particularly the new mathematical formulation takes into account the effects of vertical acceleration and curvature effects caused by complex terrain topographies the model is derived on a global cartesian coordinate system so that it is easy to apply in real world applications a godunov type finite volume method is implemented to numerically solve these new governing equations together with a novel scheme proposed to discretise the friction source terms the hydrostatic reconstruction approach is implemented and improved in the context of the new governing equations providing well balanced and non negative numerical solutions for mass flows over complex domain topographies the model is validated against several test cases including a field scale flow like landslide satisfactory results are obtained demonstrating the model s improved simulation capability and potential for wider applications

debris flows over sand dunes on mars evidence for liquid water 1 this study focuses on the formation and physical properties of the gullies observed over large martian dunes especially those of the russell crater 54 s 347 w geomorphic features like sinuosities and connections of the channels show that gullies over dunes involve flows with a significant proportion of liquid the occurrence of levees implies that these flows are debris flows with a yield strength characteristics of bingham plastic materials we apply terrestrial methods to estimate viscosity and velocity of these flows from levee size and sinuosities we obtain average velocities in the range of 1 to 7 m s 1 and apparent viscosities of 2 8 to 46 000 pa s with an average at 740 pa s compared with the 0 001 pa s of pure water these viscosities and velocities are in the range of terrestrial debris flows with a proportion of 10 to 40 of h2o these properties are typical of water holding debris flows but not of pure water surface runoff or co2 driven flows the debris flows over dunes are oriented on south facing slopes like other recent gullies meltwater from ground ice formed during a recent period of high obliquity is the more likely explanation for the formation of such flows over dunes

modeling of the steady state temperature field in lava flow levees the rationale of lava flow deviation is to prevent major damage and among the possible techniques the opening of the flow levees has often been demonstrated to be suitable and reliable the best way to open the levees in the right point in order to obtain the required effect is to produce an explosion in situ and it is then necessary to map with the highest precision the temperature field inside the levees in order to design a safe and successful intervention the levees are formed by lava flows due to their non newtonian rheology where the shear stress is lower than the yield stress the levees then cool and solidify due to heat loss into the atmosphere in this work we present analytical solutions of the steady state heat conduction problem in a levee using the method of conformal mapping for simple geometrical shapes of the levee cross section triangular or square numerical solutions are obtained with a finite element code for more complex realistic geometries

emplacement of pyroclastic flows during the 1998 1999 eruption of volcan de colima mexico after three years of quiescence volcan de colima reawakened with increasing seismic and rock fall activity that reached its peak on november 20 1998 when a new lava dome forced its way to the volcano s summit the new lava rapidly reached the s sw edge of the summit area beginning the generation of merapi type pyroclastic flows that traveled down la lumbre and the el cordoban western and eastern ravines reaching distances of 3 4 5 and 3 km respectively on december 1 1998 the lava flow split into three fronts that in early 1999 had reached 2 8 3 1 and 2 5 km in length advancing down the el cordoban ravines the lava flow fronts disaggregated into blocks forming pyroclastic flows one of the best examples occurred on december 10 1998 as the lava flow ceased moving in early 1999 activity became more explosive strong blasts were recorded on february 10 may 10 and july 17 1999 the last event developed a 10 km high eruptive column from which a pyroclastic flow developed from the base traveling 3 3 km sw from the summit into the san antonio montegrande ravines regardless of the mechanism of pyroclastic flow generation each flow immediately segregated into a basal avalanche that moved as a granular flow and an upper ash cloud in which particles were sustained in turbulent suspension when the basal avalanche lost velocity and eventually stopped the upper ash cloud continued to move independently as a dilute pyroclastic flow that produced a massive pyroclastic flow deposit and an upper dune bedded surge deposit the dilute pyroclastic flow scorched and toppled maguey plants and trees and sandblasted vegetation in the direction of the flow at the end of the dilute pyroclastic flow path the suspended particles lifted off in a cloud from which a terminal ash fall was deposited the basal avalanche emplaced block and ash flow deposits up to 8 m thick that filled the main ravines and consisted of several flow units each flow unit was massive monolithologic matrix supported and had a clast supported steep front ca 1 5 to 2 m thick composed of boulders up to 1 7 m in diameter the juvenile lithic clasts had an average density of 1800 kg m3 the dilute pyroclastic flow emplaced overbank deposits found on valley margins or beyond the tip of block and ash flow deposits they consist from bottom to top of a massive medium to coarse sand size flow layer 2 4 cm thick a dune bedded surge layer 2 10 cm thick and a massive silt size layer 0 5 cm thick the total estimated volume of the pyroclastic flow deposits produced during the 1998 1999 eruption is 24 105 m3

fine grained linings of leveed channels facilitate runout of granular flows abstract catastrophic dense granular flows such as occur in rock avalanches debris flows and pyroclastic flows move as fully shearing mixtures that have approximately 60 vol solids and tend to segregate to form coarse grained fronts and leveed channels levees restrict spreading of unconfined flows and form as coarse particles that become concentrated in the top of the flow are transported to the front and then advect to the sides in the flow head channels from which most material has drained away down slope are commonly lined with fine grained deposit widely thought to remain from the tail of the waning flow we show how segregation in experimental dense flows of carborundum or sand 300 425 m mixed with spherical fine ballotini 150 250 m on rough slopes of 27 29 produces fine grained channel linings that are deposited with the levees into which they grade laterally maximum runout distance is attained with mixtures containing 30 40 sand just sufficient to segregate and form levees that are adequately robust to restrict the spreading attributable to the low friction fines resin impregnation and serial sectioning of deliberately arrested experimental flows shows how fines lined levees form from the flow head the flows create their own stable conduit entirely from the front which in a geophysical context can play an important mechanistic role in facilitating runout the flow self organization ensures that low friction fines at the base of the segregated channel flow shear over fine grained substrate in the channel thus reducing frictional energy losses we propose that in pyroclastic flows and debris flows which have considerable mobility attributable to pore fluid pressures such fine grained flow contact zones form similarly and not only reduce frictional energy losses but also reduce flow substrate permeability so as to enhance pore fluid pressure retention thus the granular flow self organization that produces fine grained channel linings can be an important factor in facilitating long runout of catastrophic geophysical flows on the low slopes few degrees of depositional fans and aprons around mountains and volcanoes

influence of the scar geometry on landslide dynamics and deposits application to martian landslides received 24 january 2011 revised 23 june 2011 accepted 30 june 2011 published 4 october 2011 1 landslides dynamics prediction remains difficult in spite of a considerable number of studies the runout distance is widely used in analysis of landslide dynamics and in the calibration of the rheological parameters involved in numerical modeling however the unknown impact of the significant uncertainty in the shape of the initial released mass on the runoutdistanceandontheoverallshapeofthedepositraisesquestionsabouttherelevanceof these approaches the impact of the initial scar geometry on flow and distribution of the deposits is studied here using satellite data and numerical modeling of theoretical landslides and martian landslides informed by geomorphological analysis by varying the initial scar geometry from spoon shaped to steep wall geometry our results show that the runout distance is a very robust parameter that is only slightly affected by the change in the geometry of the initial scar on the contrary the lateral extent of the deposit is shown to be controlled by the scar geometry providing unique insights into the initial landsliding conditionsonmarsandmakesitpossibletoaccuratelyrecoverthevolumeinitiallyinvolved an essential ingredient for volume balance calculation a feedback analysis of valles marineris landslides can be drawn showing good agreement between numerical results and geomorphological analysis the geometry of the initial scar inferred from numerical modeling is strongly correlated with the regional tectonic history in valles marineris area

gravity driven granular free surface flow around a circular cylinder snow avalanches and other hazardous geophysical granular flows such as debris flows lahars and pyroclastic flows often impact on obstacles as they flow down a slope generating rapid changes in the flow height and velocity in their vicinity it is important to understand how a granular material flows around such obstacles to improve the design of deflecting and catching dams and to correctly interpret field observations in this paper small scale experiments and numerical simulations are used to investigate the supercritical gravity driven free surface flow of a granular avalanche around a circular cylinder our experiments show that a very sharp bow shock wave and a stagnation point are generated in front of the cylinder the shock standoff distance is accurately reproduced by shock capturing numerical simulations and is approximately equal to the reciprocal of the froude number consistent with previous approximate results for shallow water flows as the grains move around the cylinder the flow expands and the pressure gradients rapidly accelerate the particles up to supercritical speeds again the internal pressure is not strong enough to immediately push the grains into the space behind the cylinder and instead a grain free region or granular vacuum forms on the lee side for moderate upstream froude numbers and slope inclinations the granular vacuum closes up rapidly to form a triangular region but on steeper slopes both experiments and numerical simulations show that the pinch off distance moves far downstream

constraining landslide characteristics with bayesian inversion of field and seismic data using a fully nonlinear bayesian approach based on forward modelling of granular flow we invert for landslide parameters volume release geometry and rheology from different kinds of observations synthetic tests show that the runout distance and the deposit area by themselves do not constrain landslide parameters better constraints on landslide parameters are obtained from the thickness distribution of the landslide deposits as well as from the force history applied by the landslide to the ground which contains information on the landslide dynamics therefore inverting force histories calculated from seismic broad band records is an important alternative to inverting thickness distributions of landslide deposits which are usually difficult to obtain we test the method on the 1997 boxing day debris avalanche on montserrat island which involved 40 50 mm3 the bayesian inversion and granular flow model provide good estimates for volume release geometry and effective friction coefficient this study thus underlines the value of broad band seismic records as observations to monitor landslides and validation for their numerical flow models

generation of block and ash flows during the 1990 1995 eruption of unzen volcano japan processes generating block and ash flows by gravitational dome collapse merapi type pyroclastic flow were observed in detail during the 1990 1995 eruption of unzen volcano japan two different types were identified by analysis of video records and observations during helicopter flights most of the block and ash flows erupted during the 1991 1993 exogenous dome growth stage initially involved crack propagation due to cooling and flowage of the dome lava lobes the mass around the crack became unstable locally decreasing in tensile strength finally a slab separated from the lobe front fragmented progressively from the base to the top within a few seconds and became a block and ash flow rock falls immediately followed in response to local instability of the lobe front clasts in these rock falls fragmented and merged with the preceding flow in contrast block and ash flows during the endogenous dome growth stage in 1994 were generated due to local bulge of the dome unstable lava blocks collapsed and subsequently fragmented to produce block and ash flows

the indication of martian gully formation processes by slope area analysis the formation process of recent gullies on mars is currently under debate this study aims to discriminate between the proposed formation processes pure water flow debris flow and dry mass wasting through the application of geomorphological indices commonly used in terrestrial geomorphology high resolution digital elevation models dems of earth and mars were used to evaluate the drainage characteristics of small slope sections data from earth were used to validate the hillslope debris flow and alluvial process domains previously found for large fluvial catchments on earth and these domains were applied to gullied and ungullied slopes on mars in accordance with other studies our results indicate that debris flow is one of the main processes forming the martian gullies that were being examined the source of the water is predominantly distributed surface melting not an underground aquifer evidence is also presented indicating that other processes may have shaped martian crater slopes such as ice assisted creep and solifluction in agreement with the proposed recent martian glacial and periglacial climate our results suggest that within impact craters different processes are acting on differently oriented slopes but further work is needed to investigate the potential link between these observations and changes in martian climate

solution precipitation creep and fluid flow in halite a case study of zechstein z1 rocksalt from neuhof salt mine germany zechstein z1 rocksalt from the fulda basin from the immediate vicinity of the hessen potash bed is folded into tight to isoclinal folds which are cut by an undeformed 1 cm thick coarse grained halite vein microstructures were investigated in etched gamma irradiated thin sections from both the wall rock and the vein the lack of synsedimentary dissolution structures and the widespread occurrence of plate shaped and hopper grains in the wall rock suggests that the sedimentary environment was perennial lake deformation microstructures are in good agreement with solution precipitation creep process and salt flow under very low differential stress strength contrast between anhydrite rich and anhydrite poor layers caused the small scale folding in the halite beds the vein is completely sealed and composed mainly of euhedral to subhedral halite grains which often overgrow the wall rock grains those microstructures together with the presence of occasional fluid inclusion bands suggest that the crystals grew into a solution filled open space based on considerations on the maximum value of in situ differential stress the dilatancy criteria the amount of released fluids from the potash bed during metamorphism and the volume change it is proposed that the crack was generated by hydrofracturing of the rocksalt due to the presence of the salt metamorphic fluid at near lithostatic pressure

effects of pore pressure in pyroclastic flows numerical simulation and experimental validation pyroclastic flows are mixtures of volcanic gases and particles that can be very hazardous owing to their fluid like behavior one possible mechanism to explain this behavior is the reduction of particles friction due to the internal gas pore pressure to verify this hypothesis we present a numerical model of a granular flow with high initial pore pressure that decreases with time as the gas particle mixture propagates first we validate the model by reproducing laboratory experiments then the numerical code is applied to pyroclastic flows of lascar volcano 1993 eruption chile the simulation reproduces the runout and the morphological features of the deposits with lateral levees a central channel and a lobate front our results support the hypothesis of the role of gas pore pressure in pyroclastic flows and explain both the fluid like behavior of the flows and the formation of lateral levees

the controls and consequences of substrate entrainment by pyroclastic density currents at mount st helens washington usa evidence in the deposits from the may 18 1980 eruption at mount st helens demonstrates that pyroclastic density currents pdcs produced during the afternoon of the eruption became intermittently erosive using detailed componentry and granulometry we constrain the sources for lithic blocks in the deposits and identify deposits from pdcs that became locally erosive the componentry of the lithics in the fall deposits is used as a proxy for vent erosion and assumed to represent the starting componentry for pdcs prior to entrainment from any other source we find little evidence in the pdc deposits nearest to the base of the volcano for entrainment from the steep flanks however significant evidence indicates that pdcs eroded into the debris avalanche hummocks suggesting that entrainment is favored as pdcs interact with highly irregular topography evidence for locally entrained material downstream from debris avalanche hummocks decreases with height in the outcrop suggesting that less entrainment occurs as local relief decreases and upstream topography is buried the prevalence of lithofacies containing locally entrained material at the base of unit contacts and only 10s of meters downstream from debris avalanche hummocks suggests that the majority of entrainment occurs at or near the head of the current occasionally entrained material is located high above unit contacts and deposited well after the initial head of the current is inferred to have passed indicating that entrainment can occur during periods of non deposition either from the semi sustained body of the current or from a pulsating current additionally self channelization of pdcs either by levee deposition or scouring into earlier pdc deposits occurs independently of interaction with topographic obstacles and can affect carrying capacity and runout distance while we begin to explore the mechanisms and effects of erosion on current dynamics additional laboratory and numerical studies are necessary to fully understand these processes

a sph depth integrated model for popocatepetl 2001 lahar mexico sensitivity analysis and runout simulation abstract lahars are debris flows of volcanic origin which can endanger or even destroy communities located near the flanks of volcanoes lahars are not always triggered by eruptions they can form during heavy rainfall or after hydrothermal alteration and volcanic edifice collapse decision makers need lahar hazard maps to devise hazard prevention measures that will prevent casualties so lahar modelling is an important tool for assessing flow behavior and determining inundation areas the depth integrated numerical model used in this study is derived from the velocity pressure of the biot zienkiewicz model and was discretized using the smoothed particle hydrodynamics sph method to simulate a lahar that occurred at popocatepetl volcano in 2001 in order to investigate the convergence of the model we used a range of different sph mesh resolutions once the optimum mesh resolution was bounded we analyzed the model s sensitivity to the initial lahar volume the density of the geomaterial and the rheological parameter of the bingham fluid the results show that the sph depth integrated model produced a highly accurate simulation of the distribution and velocity of the 2001 lahar the study also shows the effects of sph mesh resolution and the relevant influence of rheological parameters

simulation of tsaoling landslide taiwan based on saint venant equations over general topography the 1999 chi chi earthquake triggered the catastrophic tsaoling landslide in taiwan the geomorphological change measured from the data of the 1989 and 2000 aerial photos reveals that the scar and deposit volumes are about 0 126 km3 and 0 15 km3 respectively the debris material ran over a distance of 1 6 km with 500 m descent in elevation in this paper we use the continuum model of hydraulic flow shaltop2d based on the equations of bouchut and westdickenberg to simulate numerically the landslide dynamics when the mass is moving the flow is assumed hydrostatic with a basal coulomb friction the best fit is obtained using the basal friction angle equal to 6 the only parameter of the simulation uniformly applied in the calculation domain the landslide front reaches the chinshui river valley the foothill of the slope within only about 25 s after initiation and the motion settles in about 113 s the maximum speed is estimated 75 m s the spreading of the deposit agrees well with the field measurement

suitability of simple rheological laws for the numerical simulation of dense pyroclastic flows and long runout volcanic avalanches 1 the rheology of volcanic rock avalanches and dense pyroclastic flows is complex and it is difficult at present to constrain the physics of their processes the problem lies in defining the most suitable parameters for simulating the behavior of these natural flows existing models are often based on the coulomb rheology sometimes with a velocity dependent stress e g voellmy but other laws have also been used here i explore the characteristics of flows and their deposits obtained on simplified topographies by varying source conditions and rheology the coulomb rheology irrespective of whether there is a velocity dependent stress forms cone shaped deposits that do not resemble those of natural long runout events a purely viscous or a purely turbulent flow can achieve realistic velocities and thicknesses but cannot form a deposit on slopes the plastic rheology with e g bingham or without a velocity dependent stress is more suitable for the simulation of dense pyroclastic flows and long runout volcanic avalanches with this rheology numerical flows form by pulses which are often observed during natural flow emplacement the flows exhibit realistic velocities and deposits of realistic thicknesses the plastic rheology is also able to generate the frontal lobes and lateral levees which are commonly observed in the field with the plastic rheology levee formation occurs at the flow front due to a divergence of the driving stresses at the edges once formed the levees then channel the remaining flow mass the results should help future modelers of volcanic flows with their choice of which mechanical law corresponds best to the event they are studying

can the collapse of a fly ash heap develop into an air fluidized flow reanalysis of the jupille accident 1961 a fly ash heap collapse occurred in jupille liege belgium in 1961 the subsequent flow of fly ash reached a surprisingly long runout and had catastrophic consequences its unprecedented degree of fluidization attracted scientific attention as drillings and direct observations revealed no water saturated zone at the base of the deposits scientists assumed an air fluidization mechanism which appeared consistent with the properties of the material in this paper the air fluidization assumption is tested based on two dimensional numerical simulations the numerical model has been developed so as to focus on the most prominent processes governing the flow with parameters constrained by their physical interpretation results are compared to accurate field observations and are presented for different stages in the model enhancement so as to provide a base for a discussion of the relative influence of pore pressure dissipation and pore pressure generation these results show that the apparently high diffusion coefficient that characterizes the dissipation of air pore pressures is in fact sufficiently low for an important degree of fluidization to be maintained during a flow of hundreds of meters

erosion transport and segregation of pumice and lithic clasts in pyroclastic flows inferred from ignimbrite at lascar volcano chile abstract investigations have been made on the distribution of pumice and lithic clasts in the lithic rich soncor ignimbrite 26 5 ka and the 1993 pumice flow deposits of lascar volcano chile the soncor ignimbrite shows three main lithofacies which grade into one another coarse lithic breccias range from matrix poor stratified varieties irregular shaped sheets and elongate hummocks in proximal environments to breccia lenses with pumiceous ignimbrite matrix massive lithic rich facies comprise the bulk of the ignimbrite pumice rich facies are bimodal with abundant large pumice clasts often with reverse grading rare lithic clasts and occur distally and on high ground adjacent to deep proximal valleys in the 1993 pyroclastic flow deposits lithic rich facies are deposited on slopes up to 14 whereas pumice rich facies are deposited only on slopes

field observations of pristine block and ash flow deposits emplaced april 16 17 1991 at volcan de colima mexico abstract for the first time since the effusive activity of 1981 1982 volcan de colima extruded scoriaceous andesitic lava on march 1 1991 this followed a period of increasing seismicity summit dome deformation and fumarolic activity that began in february of 1991 growth of the new lava dome proceeded until april 16 at 16 00 hrs local time gmt 6 hours when a series of gravitational failures of the upper edifice and the new lava dome flow triggered the emplacement of successive small volume non glowing block and ash flows to a distance of 4 km within the upper reaches of barranca el cordoban south southwest of the present summit aerial photos taken at the time of the event clearly show the lack of explosions at the summit dome efficient elutriation of the finest silt and clay sized fraction from the main body of the avalanching material produced large reddish yellowish and grayish ash clouds which rose to an altitude of 1 5 2 5 km for several hours and were carried southeast to northeast by predominant winds vegetation in the path of the debris was totally stripped sand blasted and bent in the sense of flowage but was not burnt no fatalities were reported we present a summary of our field observations of the pristine deposits which we studied prior the rainy season on april 20 27 between 2 200 and 2 500 m elevation

new views of granular mass flows concentrated grain fluid mixtures in rock avalanches debris flows and pyroclastic flows do not behave as simple materials with fixed rheologies instead rheology evolves as mixture agitation grain concentration and fluid pressure change during flow initiation transit and deposition throughout a flow however normal forces on planes parallel to the free upper surface approximately balance the weight of the superincumbent mixture and the coulomb friction rule describes bulk intergranular shear stresses on such planes pore fluid pressure can temporarily or locally enhance mixture mobility by reducing coulomb friction and transferring shear stress to the fluid phase initial conditions boundary conditions and grain comminution and sorting can influence pore fluid pressures and cause variations in flow dynamics and deposits

model space exploration for determining landslide source history from long period seismic data the seismic signals generated by two large volcanic debris avalanches montserrat lesser antilles 1997 and mount st helens usa 1980 and a large rock ice avalanche mount steller usa 2005 have been analyzed for the two debris avalanches given the times and locations of such landslides their signals were recorded by only a few seismic stations moreover these signals cover only a very narrow frequency band and include considerable noise the mount steller on the contrary was precisely recorded for each event the source mechanism i e point force has been determined by waveform inversion using at most two broadband seismic stations the resulting force is very difficult to interpret in terms of landslide characteristics a monte carlo inversion was therefore performed by imposing a simple force model associated with the landslide based on the schematic view of an accelerating decelerating mass traveling down the slope the best parameter set of the force model was then found by minimizing misfits and maximizing correlations between data and synthetic signals this model appears to contain the minimum degree of complexity required to well reproduce the seismic data we detail here the method for the montserrat debris avalanche and then present it s validation on the well studied mount st helens debris avalanche and the well recorded mount steller rock ice avalanche the horizontal and vertical components of the resulting force have different source time functions the best force model compares well with the force obtained by waveform inversion finally this simple force model was interpreted using analytical and empirical relations derived from the sliding block model granular flow model and landslide studies this made it possible to estimate the order of magnitude of the mass flow duration and direction initial topography slope mean velocity and travel distance of the avalanches for these three avalanches the calculated characteristics are consistent with former studies

high resolution bathymetry reveals contrasting landslide activity shaping the walls of the mid atlantic ridge axial valley 1 axial valleys are found along most slow spreading mid ocean ridges and are one of the most prominent topographic features on earth in this paper we present the first deep tow swath bathymetry for the axial valley walls of the mid atlantic ridge these data allow us to analyze axial valley wall morphology with a very high resolution 0 5 to 1 m compared to 50 m for shipboard multibeam bathymetry revealing the role played by landslides slow spreading ridge axial valleys also commonly expose mantle derived serpentinized peridotites in the footwalls of large offset normal faults detachments in our map of the ashadze area lat 13 n ultramafic outcrops have an average slope of 18 and behave as sliding deformable rock masses with little fragmentation by contrast the basaltic seafloor in the krasnov area lat 16 38 n has an average slope of 32 and the erosion of the steep basaltic rock faces leads to extensive fragmentation forming debris with morphologies consistent with noncohesive granular flow comparison with laboratory experiments suggests that the repose angle for this basaltic debris is 25 we discuss the interplay between the normal faults that bound the axial valley and the observed mass wasting processes we propose that along axial valley walls where serpentinized peridotites are exposed by detachment faults mass wasting results in average slopes 20 even in places where the emergence angle of the detachment is larger

the dependence of caribbean rainfall on the interaction of the tropical atlantic and pacific oceans abstract seasonally stratified analyses of rainfall anomalies over the intra americas sea and surrounding land areas and of onset and end dates of the central american rainy season show that the variability of the tropical atlantic sea surface temperature anomaly ssta is more strongly associated with rainfall over the caribbean and central america than is tropical eastern pacific ssta seasonal differences include the importance of antisymmetric configurations of tropical atlantic ssta in the dry season but not in the rainy season both oceans are related to rainfall but the strength of the rainfall response appears to depend on how ssta in the tropical atlantic and eastern pacific combine the strongest response occurs when the tropical atlantic is in the configuration of a meridional dipole antisymmetric across the itcz and the eastern tropical pacific is of opposite sign to the tropical north atlantic when the tropical north atlantic and tropical pacific are of the same sign the rainfall respons

evaluating titan2d mass flow model using the 1963 little tahoma peak avalanches mount rainier washington the titan2d geophysical mass flow model is evaluated by comparing its simulation results and those obtained from another flow model flow3d with published data on the 1963 little tahoma peak avalanches on mount rainier washington the avalanches totaling approximately 10 10 6 m 3 of broken lava blocks and other debris traveled 6 8 km horizontally and fell 1 8 km vertically h l 0 246 velocities calculated from runup range from 24 to 42 m s and may have been as high as 130 m s while the avalanches passed over emmons glacier titan2d is a code for an incompressible coulomb continuum it is a depth averaged dshallow watert granular flow model the conservation equations for mass and momentum are solved with a coulomb type friction term at the basal interface the governing equations are solved on multiple processors using a parallel adaptive mesh godunov scheme adaptive gridding dynamically concentrates computing power in regions of special interest mesh refinement and coarsening key on the perimeter of the moving avalanche the model flow initiates as a pile defined as an ellipsoid by a height z and an elliptical base defined by radii in the x and y planes flow parameters are the internal friction angle and bed friction angle results from the model are similar in terms of velocity history lateral spreading location of runup areas and final distribution of the little tahoma peak deposit the avalanches passed over the emmons glacier along their upper flow paths but lower in the valley they traversed stream gravels and glacial outwash deposits this presents difficulty in assigning an appropriate bed friction angle for the entire deposit incorporation of variable bed friction angles into the model using gis will help to resolve this issue d 2004 elsevier b v all rights reserved

numerical modeling of landquakes the thurwieser landslide that occurred in italy in 2004 and the seismic waves it generated have been simulated and compared to the seismic signal recorded a few tens of kilometers from the landslide source i e landquake the main features of the low frequency seismic signal are reproduced by the simulation topography effects on the flowing mass have a major impact on the generated seismic signal whereas they weakly affect low frequency wave propagation simulation of the seismic signal makes it possible to discriminate between possible alternative scenarios for flow dynamics and to provide first estimates of the rheological parameters during the flow as landquakes are continuously recorded by seismic networks our results provide a new way to collect data on the dynamics and rheology of natural flows

constraints on eolian grain flow dynamics through laboratory experiments on sand slopes abstract to illuminate both the dynamics of eolian grain flows and the controls on their geometry we conducted experiments on an inclined platform on which grain flows were initiated by broadcasting sand near the top of the platform in a pattern consistent with wind blown deposition the growing depositional bumps eventually over steepened and failed generating steep headward and laterally migrating scarps mass was fed to the grain flow in this way for approximately 1 4 of the duration of the flow after which the flows were solitary waves spreading both downslope and cross slope the flows abandoned small levees at their widest points and eventually stopped commonly before reaching the base of the platform our grain flows were typically 2 m long by 20 cm wide and slightly less than 1 cm thick the surface velocities at the toe were 15 cm s increased to 20 cm s at the nose and then steadily decreased to 10 cm s towards the tail of the flow the longitudinal velocity gradients created considerable longitudinal strain rates the toe was in compression with mean strain rates of order 0 1 s while the tail was in extension with mean strain rates of order 0 1 s averaged over many seconds the velocities were quite steady over shorter periods however distinct velocity fluctuations dominated the record with variations up to 100 and s andard deviations 10 20 of the mean there were both temporal and spatial correlations in velocity histories of surface tracers we interpret spatial correlations to reflect fixed topographic perturbations at the self defined bottom boundary or ramps and temporal correlations to indicate coherent regions within the flow or rafts ramps were more common than rafts the longitudinal length scales of these features may be up to many times the thickness of the flow and endure up to two seconds based on the average velocity and the magnitude of the velocity response to perturbations at the bed the slightly rate dependent coulomb rheological model proposed by savage and hutter 1991 best describes the general behavior of these flows grain flows stop when the toe either thins below a critical thickness or reaches the toe of the slope the thickness of the nose of a grain flow is controlled by a combination of gradients in longitudinal volume flux and volume losses to the levees at the edge of the flow the nose eventually thins sufficiently to stop the initial spatial distribution of velocities is set by the volume flux history of the evolving scarp within the depositional bump this sets the spatial distribution of strain within the flow which in turn controls the thickness of the nose assuming that the sizes and geometries of depositional bumps are not a function of dune size we expect a minimum flow thickness to occur on dunes whose lee face length equals the natural runout length of a typical grain flow for dunes taller than this grain flow thickness should increase because of coalescence of flow tongues at some mid point on the lee face loading the mid face region with coalescing grain flow tongues that subsequently fail as a single large flow

the emplacement dynamics of pumice lobes ascertained from morphology and granulometry examples from the 1993 deposits at lascar volcano chile abstract the work presented here focuses on lobe shapes and clast populations within lobate termini of the 1993 pumice flow deposits at lascar volcano chile a new method to analyze a coarse tail grain size population with field photographs is presented using this method 33 000 0 5 cm clasts from the pumice lobes of the 1993 pumice flow deposits were measured at 36 sites and the resultant grain size distributions were then related to lobe morphology lobe margins i e levees clefts and snouts were found to contain significantly larger pumice clasts and be more poorly sorted than lobe central channels i e locations away from the margins previous laboratory experiments suggest lobe margins form by the floatation and deflection of larger clasts to the margins of an advancing flow lobe results here indicate that the same sorting process efficiently segregates clasts into two flow regimes 1 a mobile central channel depleted in coarse clasts and 2 friction dominated margins enriched in clasts 15 cm the lobe margins 60 enriched in larger particles with matrix

transferability of a calibrated numerical model of rock avalanche run out application to 20 rock avalanches on the nuussuaq peninsula west greenland long run out rock avalanches are one of the most hazardous geomorphic processes and risk assessments of the potential threat they pose are often reliant on numerical modelling of their potential run out distance the development of such models requires a thorough understanding of past flow behaviour inferred from deposits emplaced by previous events despite this few records exist of multiple rock avalanches that occurred in conditions sufficiently consistent to develop a set of more generalised and hence transferrable rules we conduct field and imagery based mapping and use numerical modelling to investigate the emplacement of 20 adjacent rock avalanches on the southern flanks of the nuussuaq peninsula west greenland the rock avalanches run out towards the vaigat strait and are sourced from a range of coastal mountains of relatively uniform geology we calibrate a three dimensional continuum dynamic flow code volcflow with data from a modern well constrained event that occurred at paatuut ad 2000 the best fit model assumes a constant retarding stress with a collisional stress coefficient simulating run out to within 0 3 of that observed this calibration was then used to model the emplacement of deposits from five other neighbouring rock avalanches before simulating the general characteristics of a further 14 rock avalanche deposits on simplified topography our findings illustrate that a single calibration of volcflow can account for the observed deposit morphology of a uniquely large collection of rock avalanche deposits emplaced by a series of events spanning a large volume range although the prevailing approach of tuning models to a specific case may be useful for detailed back analysis of that event we show that more generally applied models even using a single pair of rheological parameters can be used to model potential rock avalanches of varied volumes in a region and therefore to assess the risks that they pose

influence of slope angle on pore pressure generation and kinematics of pyroclastic flows insights from laboratory experiments the influence of slope angle on pore pressure generation and kinematics of fines rich pyroclastic flows was investigated through laboratory experiments granular flows were generated by the release of a column of fine glass beads d 0 08 mm in an inclined channel 0 30 the granular column could be fluidized while the channel base was either smooth or made rough by glued beads of 3 mm diameter pore pressure measurements reveal that the degree of autofluidization caused by air escaping from the substrate interstices into which flow particles settled was high at all slope angles flow runout increase due to autofluidization however was reduced at slope angle higher than 12 because of the occurrence of a strong deceleration phase that limited the flow duration this is probably caused by the combination of flow head thinning at increased slope angle and settling of particles into the substrate interstices until the flow ran out of mass analysis of high speed videos suggests that ingestion of ambient air at the flow front did not occur even on steep slopes of 30 experiments at inclinations close to 25 or slightly higher 30 than the repose angle of the granular material 28 5 revealed the formation of a thin basal deposit that was then eroded as the flow thickness and velocity gradually decreased our study suggests that air escape from substrate interstices in nature can be a significant external cause of pore pressure generation that favors low energy dissipation and long runout distances of pyroclastic flows on moderate topographies

channeled flow analytic solutions laboratory experiments and applications to lava flows although channeled lava flows are common in basaltic volcanism relationships between channel morphology eruption and emplacement parameters and lava properties are not well understood several models have commonly been used to constrain these relations but they have not been well tested on natural or simulated lava flows over a wide range of parameter space here we test the accuracy and assumptions of a moderately simple analytic rectangular channel solution by comparing the behavior of well controlled laboratory polyethylene glycol peg channeled flows to the analytic solution for isothermal steady newtonian flow in a rectangular channel with constant dimensions this analytic solution agrees well with laboratory measurements volumetric effusion rates q m3 s 1 calculated from the analytical model using measured peg flows as input yield ratios of qcalculated qpumped of 0 2 to 3 6 and flow rates calculated from a best fit surface velocity profile to measured velocities give more accurate ratios of 0 8 to 1 2 we find a very weak dependence of solution accuracy on slope attributable to flow front effects within the laboratory flows we subsequently apply the solution to several subaerial and submarine terrestrial flows as well as extraterrestrial channeled flows over a wide range of flow parameters viscosity ranges and flow rates obtained using measured channel dimensions and assumed lava properties are plausible interestingly the resulting extraterrestrial estimates of viscosities and flow rates tend to fall closer to known terrestrial measurements and estimates of channel flow than to previous planetary estimates we therefore suggest that the analytic newtonian rectangular channel flow model is a more appropriate physical model for many channeled terrestrial and planetary flows than the newtonian infinite sheet flow and approximation to bingham channel flow widely used previously

relation between dry granular flow regimes and morphology of deposits formation of levees in pyroclastic deposits experiments on dry granular matter flowing down an inclined plane are performed in order to study the dynamics of dense pyroclastic flows the plane is rough and always wider than the flow focusing this study on the case of laterally unconfined free boundary flows we found that several flow regimes exist depending on the input flux and on the inclination of the plane each flow regime corresponds to a particular morphology of the associated deposit in one of these regimes the flow reaches a steady state and the deposit exhibits a levee channel morphology similar to those observed on small pyroclastic flow deposits the levees result from the combination between lateral static zones on each border of the flow and the drainage of the central part of the flow after the supply stops particle segregation features are created during the flow corresponding to those observed on the deposits of pyroclastic flows moreover the measurements of the deposit morphology thickness of the channel height of the levees width of the deposit are quantitatively related to parameters of the dynamics of the flow flux velocity height of the flow leading to a way of studying the flow dynamics from only measurements of the deposit some attempts to make extensions to natural cases are discussed through experiments introducing the polydispersity of the particle sizes and the particle segregation process