flowgo a kinematic thermo rheological model for lava flowing in a channel we present a kinematic self adaptive numerical model to describe the down flow thermal and rheological evolution of channel contained lava as our control volume of lava advances down a channel it cools and crystallizes an increasingly thick and extensive surface crust grows and its heat budget and rheology evolve by estimating down flow heat and velocity loss our model calculates the point at which the control volume becomes stationary giving the maximum distance lava flowing in the channel can extend modeled effusion rates velocities widths surface crust parameters heat budget cooling temperature crystallinity viscosity and yield strength all compare well with field data collected during eruptions at mauna loa k lauea and etna modeled lengths of 25 27 2 5 5 7 and 0 59 0 83 km compare with measured lengths of 25 27 4 and 0 75 km for the three flows respectively over proximal flow portions we calculate cooling crystallization viscosity and yield strength of 1 10 c km 1 0 001 0 01 volume fraction km 1 103 104 pa s and 10 3 102 pa respectively at the flow front cooling crystallization viscosity and yield strength increase to 100 c km 1 0 1 volume fraction km 1 106 107 pa s and 103 104 pa respectively all of which combine to cause the lava to stop flowing our model presents a means of a analyzing lava flow thermo rheological relationships b identifying important factors in determining how far a channel fed flow can extend c assessing lava flow hazard and d reconstructing flow regimes at prehistoric unobserved or remote flows

long submarine lava flows observations and results from numerical modeling long 100 km lava flows are relatively common on mars and venus and have been identified on the moon but they are rarely documented on earth however although 75 of the earth s surface is covered by water only a small percentage of the ocean floor has been investigated at a resolution sufficient to unequivocally identify the boundaries of long submarine lava flows even so basaltic lava flows as long as 110 km have been identified on the deep 1500 m seafloor near hawaii and the east pacific rise ambient conditions on the deep ocean floor may favor the development of long lava flows for the following reasons first high pressures 15 mpa keep volatiles dissolved in basaltic lavas preventing viscosity increases associated with exsolution and vesiculation second seawater rapidly quenches the surface of submarine basalt flows so that an insulating glass layer 1 5 cm thick encases submarine flows within seconds after their emplacement this glass rind effectively insulates the molten flow interior from additional heat loss making submarine basalt flows behave as well insulated subaerial tube fed flows thus for identical basalt flows emplaced on the deep seafloor and subaerially a submarine flow could advance farther before stopping results of numerical modeling indicate that thin 1 m submarine basalt flows behave similarly to identical subaerial flows but thicker submarine flows may advance significantly farther than their subaerial counterparts

models of seismic anisotropy in the deep continental lithosphere seismological observations sks wave polarizations systematic p residual variations azimuthal dependence of pn and surface wave velocities or a dispersion of surface waves are not consistent with isotropic if laterally heterogeneous upper mantle structure therefore an anisotropy should be considered as an a priori aspect of future large scale studies of mantle structure most studies of anisotropy however have assumed horizontal or vertical axes of symmetry but such orientations cannot explain bipolar patterns of spatial variations of p residuals which we have observed at many seismological stations on the basis of the petrophysical properties of real upper mantle rocks we consider anisotropy formed either by hexagonal or by orthorhombic aggregates composed of olivine orthopyroxene and clinopyroxene rotations of the aggregates about vertical and horizontal axes allow us to find the three dimensional orientations of symmetry axes that fit combinations of both p and s seismological observations in central europe and in western north america the orientations with plunging symmetry axes velocity extremes seem to be consistent across large spatially uniform tectonic units and change abruptly at important suture zones

thermal erosion by laminar lava flows thermal erosion by hot lavas during their laminar flow over cold ground is investigated both theoretically and experimentally an analysis of the steady forced convective heat transfer by a laminar channel flow at large peclet numbers is presented which is used to determine the final steady state erosion velocity as well as the thicknesses and timescales of the associated thermal boundary layers in both the lava and ground the initial transient period involving the growth and remelting of a basal chill layer is also quantified laboratory experiments are described in which hot molten wax flowed over and thermally eroded an underlying bed of solidified wax the experimental observations of chill formation and thermal erosion are found to be in excellent agreement with the theoretical predictions the theory is applied to predict rates of thermal erosion consistent with those observed in recent flows of basaltic lava in hawaii and of carbonatite lava in tanzania

improved remote sensing estimates of lava flow cooling a case study of the 1991 1993 mount etna eruption infrared radiance data acquired by the landsat thematic mapper tm provide information regarding the thermal state of active lavas which can be used to estimate flow surface temperatures and thus lava cooling saturation of one or both of tm s short wave infrared swir wavebands has meant however that previous studies of lava flow cooling have approximated what is almost certainly an n thermal component flow surface temperature distribution where n 1 with a single pixel integrated temperature pit we show that the high sensitivity of tm s swir wavebands to the small amount of high temperature e g 1000 c material exposed in hot fractures within the flows cooler e g 200 c surface crust results in the calculation of pits that are higher than the temperature of this crust which is usually the dominant radiative component this causes heat loss from such pixels to be overestimated we present adaptations that allow more realistic heat loss estimates to be produced using tm data by assuming that a second high temperature component contributes to the radiance measured from saturated pixels and by defining an appropriate range for the temperature of the lava crust by analyzing a suite of five tm images acquired during the 1991 1993 eruption of mount etna we present more realistic estimates of how this flow cooled as well as results that show how the flow regime of the lava evolved over time the methods we present will continue to be of use as although the next generation of high spatial resolution sensors such as the landsat 7 enhanced thematic mapper plus etm and the terra advanced space borne thermal emission and reflection radiometer aster will be less susceptible to saturation than tm the problem will still persist

subsidence in the paricutin lava field causes and implications for interpretation of deformation fields at volcanoes abstract assessment of volcanic hazards includes interpretation of ground deformation signal which at polygenetic volcanoes often results from the superposition of deformation due to pressure changes in the magmatic system and due to surficial processes such as cooling of emplaced lava the deformation signal associated with emplaced lava is sometimes considered negligible if fields are decades old but if the lava thickness is great deformation may still be occurring possibly leading to misinterpretation of the observed deformation here i evaluate the 2007 2011 ground motion of the 1943 1952 lava field of the paricutin monogenetic cinder cone mexico interferometric synthetic aperture radar insar time series reveal patchy subsidence restricted to the lava field and following linear rates up to 5 5 cm year there is a clear correlation between subsidence rates and topography suggesting a causal relationship with deposits or lava thickness i estimate these thicknesses in the subsiding areas using pre and post eruption topographic maps and show that they reach up to 200 m a numerical model for lava flow cooling was developed considering radiation and convection from the surface conductive transfer inside the flow and to the ground and vesiculation and latent heat generation at the top and bottom of the flow the model shows that compaction induced by cooling of the thick deposits emplaced 60 years ago explains the observed subsidence when conductive transfer to the ground is considered these results demonstrate that thick deposits can keep deforming significantly even decades after their emplacement emphasizing the importance of considering cooling processes when interpreting deformation fields at polygenetic volcanoes producing massive lava fields

fracture and surface crust development in a holocene pahoehoe lava flow on the island of tenerife canaries abstract an almost horizontal pahoehoe surface in a holocene plagioclase basanite lava on tenerife displays three scales of fracture within the surface crust an early formed set of large scale fractures divides up the surface into an orthogonal set of rectangular slabs with dimensions of several metres and depths of 10 12 cm the shortest slab dimension is parallel to the flow emplacement direction inferred from a strong surface lineation the slabs are domed with the centre an average of 9 6 cm with range 4 19 6 cm above the edges of the slabs profiles of the slabs normal and perpendicular to the margins and through the crest indicate that they can be described by a power law in which the deflection of the slab h is related to the distance from the crest x with an exponent between 2 and 3 analysis of joints within the slabs indicates two smaller scale networks an intermediate scale joint network bounds blocks with rectilinear to polygonal shapes in plan view and has a characteristic mean spacing of 24 2 cm range 10 5 48 cm the major fractures in this set are normal and parallel to the slab margins a smaller scale joint network bounds polygonal equant blocks in plan view and has characteristic spacing of 6 4 cm range 3 7 10 5 cm a model of cooling from the pahoehoe surface is used to constrain the growth of the crust and timing of fracture development the large scale slabs are attributed to localised accumulation of gas beneath the growing crust causing buoyant forces the tensile stresses caused by uplift are sufficient to form the large scale fractures after 2 or 3 h of cooling the intermediate scale fracture network is attributed to the flexure of the slab crust the smaller scale polygonal joint network is related to the build up of isotropic tensile stresses in the cooling slab crust due to thermal contraction with fracture development being promoted by the flexure of the slabs an analysis of the slab deformation indicates that lava crust is weak the weakness is explained by division of the crust into three zones an outer zone with small scale joints that cause negligible strength a middle zone of elastic behaviour in which stress can build up and a lower zone of plastic deformation the crustal slabs display profiles similar to that expected in a bending elastic plate the deformation of the 10 cm thick crust can be explained if the elastic zone was about 2 cm thick this result agrees with an independent calculation of elastic zone thickness based on the position of the brittle ductile transition being located at the 600 c isotherm at a depth of about 2 cm when the crustal slabs were rifted apart

post emplacement cooling and contraction of lava flows insar observations and a thermal model for lava fields at hekla volcano iceland lava flows contract as they cool causing progressive subsidence of the flow surface here we study this process by measuring and modelling the deformation of emplaced lava flows and the surrounding substrate the temporal trend of vertical lava movements was investigated using interferometric analysis of synthetic aperture radar insar images from the 1991 and 2000 hekla eruptions covering periods of 23 and 12 years respectively data from six tracks from three satellites including both ascending and descending passes were used to create 99 interferograms from which trends of accumulated subsidence and subsidence velocities were derived subsidence rates are similar for both lava flows and decay approximately exponentially from about 20 mm year five years after emplacement to about 2 mm year 15 years after emplacement a one dimensional semi analytical model was fitted to the observed subsidence rates with subsidence due to phase change calculated analytically and subsidence due to thermal contraction calculated numerically using dilatometic constraints obtained experimentally the initial thicknesses of the 1991 and 2000 lava fields d1991 and d2000 scaled thermal expansivity and thermal diffusivity are the crucial parameters influencing lava subsidence and subsidence rate inversion for these parameters reveal linear correlations between them best fitting results of inversions for d1991 range from 10 m to 27 m for d2000 from 10 m to 30 m 9 24 10 6k 1 and 1 7 10 7m2s 1

the occurrence vesiculation and solidification of dense blue glassy pahoehoe dense blue glassy pahoehoe occurs on the lower kalapana flow field of hawaii as limited breakouts at the base of tumuli and pressure rises or in localized depressions dense blue glassy pahoehoe differs from the more abundant spongy pahoehoe in that it has a distinctive silvery blue color a poorly vesiculated glassy outer selvage and lower extrusion temperature 10 c cooler the color of blue glassy pahoehoe may be attributed to the concentration of iron and magnesium at the surface whereas the golden color of the spongy pahoehoe is possibly the result of fluorine enrichment the outer margin of the dense pahoehoe consists of three layers developed in response to an increased cooling rate toward the surface the outer glass selvage is less than 1 cm thick containing minor amounts of chromite olivine plagioclase and augite inside the glassy rind is a 5 mm thick transition zone characterized by increasing numbers and sizes of quenched plagioclase crystals and augite clusters the interior of dense blue glassy pahoehoe is largely crystalline crystal size distribution csd theory demonstrates that the formation and textures of plagioclase pyroxene and oxides are related to the cooling rates growth rates g and nucleation rates j present in each layer of the dense pahoehoe in contrast to the outer selvage the crystalline interior typically contains numerous large coalesced vesicles as compared to the evenly vesiculated spongy pahoehoe the enigmatic contrast between the poorly vesicular outer glassy selvage and the vesicular interior has been attributed to both lower and higher vapor contents of the dense pahoehoe the higher density lower extrusion temperatures degassed outer glassy selvage well developed crystalline interior and fluorine deficient surface of dense blue glassy as compared to spongy pahoehoe indicate that blue glassy pahoehoe was stagnant long enough to exsolve considerable vapor in the lava tube system before extrusion

spongy pahoehoe in hawaii a study of vesicle distribution patterns in basalt and their significance a well defined pahoehoe lava type that is very common medially and distally in hawaii is characterized by a high concentration and fairly uniform distribution of spherical or near spherical vesicles measurements of vesicle sizes and concentrations have been made on ten of these spongy pahoehoe lava flow units the vesicles increase in size toward the middle accompanied by a moderate increase in lava porosity the close approach to bilateral symmetry on either side of the horizontal median plane and the common occurrence of a median gas blister shows that no significant upward movement of vesicles occurred suggesting that the lava possessed a yield strength and was more or less static olivine phenocrysts when present are however concentrated in the lower half of the same flow units showing that the lava previously lacked a significant yield strength the vesicles are regarded as early inherited from the vent but the size characteristics of the vesicle population are a late formed feature vesicles grew in static lava mostly by coalescence and it is postulated that coalescence was promoted by the presence of abundant diktytaxitic voids which punctured the walls of contiguous vesicles zones in which the vesicle concentration is lower and the vesicles are larger and strongly deformed interrupt the symmetry of some spongy pahoehoe units and gas blisters higher than the median plane occur in many examples these zones are interpreted to result from late stage shearing and point to a mechanism by which vesicles may be eliminated from a lava

ambient effects on basalt and rhyolite lavas under venusian subaerial and subaqueous conditions both subaerial and subaqueous environments have been used as analog settings for venus volcanism to assess the merits of this the effects of ambient conditions on the physical properties of lava on venus the seafloor and land on earth are evaluated rhyolites on venus and on the surface of earth solidify before basalts do because of their lower eruption temperatures rhyolite crust is thinner than basalt crust at times less than about an hour especially on venus at later times rhyolite crust is thicker because of its lower latent heat relative to basalt the high pressure on the seafloor and venus inhibits the exsolution of volatiles in lavas vesicularity and bulk density are proportional so that lavas of the same composition should be more dense on the seafloor and less dense on land because viscosity depends partly upon the fraction of unvesiculated water in a melt basalts with the same initial volatile abundance will be least viscous on the seafloor and most viscous on land assuming the same preemptive h2o contents molten rhyolites on venus will have viscosities 10 that of rhyolites on land despite lower expected viscosities underwater flows are more buoyant and should have heights like subaerial and venusian lavas of the same composition and extrusive history in cases where the influence of crust is insignificant a volume of rhyolite will have a higher aspect ratio than the same volume of basalt no matter what the environment if flow rheology is dominated by the presence of strong crust aspect ratios differ little among environments or between compositions these analyses support a rhyolitic interpretation for the composition of venusian festooned flows and a basaltic interpretation for the composition of venusian steep sided domes although ambient effects are significant extrusion rate and eruption history must also be considered to explain analogous volcanic landforms on earth and venus

numerical modeling of lava flow cooling applied to the 1997 okmok eruption approach and analysis 1 throughout february and march of 1997 okmok volcano in the eastern aleutian islands of alaska erupted a 6 km long lava flow of basaltic a a within its caldera a numerical model for lava flow cooling was developed and applied to the flow to better understand the nature of its cooling radiation and convection from the surface as well as conduction to the ground were used to transport the flow s heat to its surroundings in the model internally a conduction only approach moved heat from the interior outward vesiculation latent heat generation and thermal conductivity changes with temperature are among the other factors that were dynamically accounted for results indicate that ambient temperature fluctuations on the scale of days to weeks must be taken into account to create an accurate short term prediction of lava surface temperature daily data of rainfall and ambient temperature as opposed to yearly averages greatly increased the accuracy of the model furthermore convective cooling of the lava surface was observed to be a dominant heat loss process during the first 200 days indicating the convective heat transfer coefficient is a prime determinant of the accuracy of the model for predicting surface temperatures over a longer cooling period 2 years thermal conductivity and porosity proved to be among the dominating factors for heat loss because of the limiting role of conductive heat flow in the interior the model s flexibility allows application to flows other than the 1997 okmok eruption

thermal and dynamical regimes of single and two phase magmatic flow in dikes finite element calculations of magma flow in dikelike channels with length to width ratios of 1000 1 or more have been used to investigate the coupling between thermal and dynamical regimes due to temperature dependent viscosity and dissipation steady state solutions with realistic thermal and dynamical parameter values have been obtained the models show that the onset of solidification on the boundaries of a basaltic or andesitic dike as predicted by idealized laminar flow models can be prevented or significantly delayed by a small amount of transverse flow induced by rising bubbles boundary roughness or turbulence this effect will reduce the critical initial widths of dikes estimated by bruce and huppert 1989 in the absence of transverse flow the bulk temperature of the magma may actually increase slightly with distance along the dike as a result of viscous dissipation even while solidification is occurring on the walls of the dike with converging or necking dike walls it is found that boundary temperatures fall to a minimum and then increase with distance along a dike even if viscous heating is neglected for viscous heating to offset significant rates of heat loss 2 kw m2 in a plane parallel 1 m wide basaltic dike an average flow velocity of 2 7 m s driven by a pressure gradient of 1 7 mpa km is required a 7 void fraction caused by exsolution of volatiles or chamber gas in the magma will produce this pressure gradient the ease of producing such a gradient by reducing the density of the magmatic column with addition of a gas phase makes it likely that flows of basaltic magma could be maintained in dikes tens of kilometers long furthermore a gas phase may be important for the propagation of the fracture prior to the initial injection of magma during dike emplacement rapid transport of more silicic magmas through dikes is inhibited by the requirement of large driving pressure gradients exceeding several hundred megapascals per kilometer however the pressure gradient can be substantially reduced if the magma is a heterogenous mixture of a predominantly silicic component with a more mafic component pipe flow experiments involving molten polymers that exhibit dynamical similarity to magmas strongly suggest that unmixing occurs when a two component magma mixture in which the components have different viscosities rises within a dike within a few dike widths of the inlet the less viscous mafic component encapsulates the more viscous silicic part effectively lubricating the passage of the more viscous component compositional variations in the obsidian dome volcano support the occurrence of this pressure reduction or self lubrication mechanism such a self lubrication process often may be necessary to permit very viscous magmas to reach earth s surface if so chemical and lithologic zoning would be anticipated as a common feature of near surface intrusions or lava flows that are characterized by a high silica and or high crystal content such as at deadman dome in long valley california

solidification dynamics in channeled viscoplastic lava flows 1 the influence of a viscoplastic lava rheology on the dynamics of channeled lava flows is analyzed using analogue experiments the experiments used slurries of polyethylene glycol and kaolin which flowed with a constant flux down an inclined channel under water three sets of complementary experiments are presented isothermal cooling and solidifying flows which quantified the effects of the viscoplastic rheology on shear internal convection and surface crust formation the isothermal and cooling experiments showed the formation of unyielded central plug regions which were not broken up by the convective overturning in the solidifying experiments flows fell into one of three regimes a tube regime in which crust covered the entire flow surface a shear controlled regime with a mobile raft of crust in the channel center and open shear zones near the walls and a plug controlled regime where the width of the central crust was determined by the width of the central plug region the crust coverage is parameterized in terms of two dimensionless parameters the ratio wp of central plug region width to channel width and a parameter which characterizes the relative importance of the strain and solidification rates finally the dynamics of typical lava flows on mt etna and the 1984 mauna loa lava flow are examined we show that our parameterization agrees with lava flow crust widths observed in the field and find that even small yield strengths have a major effect on crust coverage

conclusion recommendations and findings of the red seed working group red seed stands for risk evaluation detection and simulation during effusive eruption disasters and combines stakeholders from the remote sensing modelling and response communities with experience in tracking volcanic effusive events the group first met during a three day long workshop held in clermont ferrand france between 28 and 30 may 2013 during each day presentations were given reviewing the state of the art in terms of a volcano hot spot detection and parameterization b operational satellite based hot spot detection systems c lava flow modelling and d response protocols during effusive crises at the end of each presentation set the four groups retreated to discuss and report on requirements for a truly integrated and operational response that satisfactorily combines remote sensors modellers and responders during an effusive crisis the results of collating the final reports and follow up discussions that have been on going since the workshop are given here we can reduce our discussions to four main findings 1 hot spot detection tools are operational and capable of providing effusive eruption onset notice within 15 min 2 spectral radiance metrics can also be provided with high degrees of confidence however if we are to achieve a truly global system more local receiving stations need to be installed with hot spot detection and data processing modules running on site and in real time 3 models are operational but need real time input of reliable time averaged discharge rate data and regular updates of digital elevation models if they are to be effective the latter can be provided by the radar photogrammetry community 4 information needs to be provided in an agreed and standard format following an ensemble approach and using models that have been validated and recognized as trustworthy by the responding authorities all of this requires a sophisticated and centralized data collection distribution and reporting hub that is based on a philosophy of joint ownership and mutual trust while the next chapter carries out an exercise to explore the viability of the last point the detailed recommendations behind these findings are detailed here

the deep lithosphere in the alps a model inferred from p residuals abstract the tomographic studies of the alps and their surroundings based on the relative p residuals indicated two lithosphere roots attaining depths over 200 km separated by an upwelling of the lithosphere asthenosphere transition beneath the central part of this mountain range the directional dependence of relative p velocities indicates an anisotropic wave propagation in the subcrustal lithosphere beneath the alps and their surroundings the presented model of the deep lithosphere structure suggests several phases of the tectonic development of this region collision of a tip of the adriatic plate with europe in the eastern alps and dextral strike slip movement of the adriatic ending with underthrusting of europe beneath the western alps and counterclockwise rotation the latter contributes to a steep plunge of both plates beneath the eastern alps and to stress accumulation in the friuli region a model with delamination of the rigid upper crust and its rotation resolves some of the inconsistencies between the near surface and deep tectonics and also explains the distribution of seismicity

calculated effect of vesicles on the thermal properties of cooling basaltic lava flows abstract this paper investigates the effects of porosity in the form of vesicles and bubbles on the transfer of heat within basaltic lava this study was undertaken to provide input for realistic cooling models of basaltic lava flows and to help explain some recent field measurements by jones 1992 1993 these field observations indicate that the surfaces of vesicular pahoehoe flows cool significantly more rapidly than that of dense flows it has been suggested that thermal radiation across vesicles is responsible for this enhanced cooling rate it is shown here that for vesicles in typical pahoehoe flows radiation across vesicles may enhance the effective thermal conductivity and thermal diffusivity but only at high temperatures t 800 c and high vesicularities 40 it is also found that convection of the gas phase within bubbles will not occur unless the cavities are larger than about 1 cm furthermore previous work has demonstrated that porosity greatly reduces the thermal conductivity and thermal inertia of cold lava it is shown here that this should also be true at high temperatures examining only radiation across vesicles severely understates the effect of porosity on the cooling of the surface of lava flows thermal inertia developed to describe diurnal surface temperature variations is the most appropriate thermal parameter in this case thermal inertia is a measure of how quickly surface temperature responds to heating or cooling with low thermal inertia materials responding more quickly it is shown here that porosity greatly reduces thermal inertia at all temperatures and that this provides a more general explanation for the field observation

thermal inertia mapping of mars from 60 s to 60 n twenty micrometer brightness temperatures are used to derive the thermal inertia for 81 of the martian surface between latitudes 60 these data were acquired by the two viking infrared thermal mappers in 1977 and 1978 following the two global dust storms of 1977 the spatial resolution used is 2 in latitude by 2 in longitude and the total range in derived inertia is 1 to 15 10 3 cal cm 2 sec 12 k 1 the distribution of thermal inertia is strongly bimodal with all values of thermal inertia less than 4 10 3 cal cm 2 sec 12 k 1 being associated with three disjoint bright regions mostly in the northern hemisphere sufficient dust is raised in global storms to provide fine material adequate to produce these low inertia areas but the specific deposition mechanism has not been defined at the low resolution used no complete exposures of clean rock were found there is some tendency for darker material to be associated with higher thermal inertia although the trend is far from one to one the distribution of high and low inertia areas is sufficiently nonrandom to produce a variation in whole disk brightness temperature with central meridian longitude this variation and the change in surface kinetic temperature associated with dust storms are factors in establishing the whole disk brightness temperature at radio and infrared wavelengths and will be important for those who use mars as a calibration source

p type and s type pahoehoe a study of vesicle distribution patterns in hawaiian lava flows abstract two main types of hawaiian pahoehoe lava p type pipe vesicle bearing and s type spongy are characterized their origin is inferred from features of the vesicle sizes and abundances particularly in their chilled margins which record frozen in the vesicle population at the time and place of solidification s type margins are distinctly more vesicular than p type porosity and vesicle size and range in s type margins are highest at the vent vesicle size reaches a minimum in medial areas and then increases while porosity decreases onto the flatter ground of the coastal terrace the vesicles are inferred to be part of the complement of bubbles that was present at the vent modified by loss and coalescence during travel in lava tubes the minimum porosity occurs on the coastal terrace and is the result of degassing during a longer residence time in lava tubes there the systematically lower porosity of p type margins than s type is attributed to a still longer residence allowing more time for bubble coalescence and loss p type pahoehoe is common only on flatter ground such as the coastal terrace it is found where tumuli and lava rises are widespread and where extensive networks of lava expansion clefts occur through which lava moving below the surface crust can lose gas without significant cooling additionally p type pahoehoe shows clear evidence that a strong inward movement and growth of bubbles took place from the outermost selvage into an inner selvage concentration zone a mechanism that might cause this inward movement is surface tension operating against deforming forces in the steep viscosity gradient of flow margins causing down gradient bubble translation inward movement does not usually occur in the upper crust of thick 2 m flow units because of the absence of strong deforming forces in the shallow velocity gradient there many of the larger bubbles that entered the inner selvage rose buoyantly to form pipe vesicles p type pahoehoe is characterized by sampled profiles across 20 flow units continued gas loss through surface clefts caused a general overall reduction in vesicle content within each flow unit to less than that in the upper crust and the pronounced depletion of vesicles from its lower half is attributed to scavenging by ascending pipe bubbles in a small proportion of flow units overall gas loss was small and a median gas blister fed by pipe bubbles developed instead

risk evaluation detection and simulation during effusive eruption disasters lava ingress into a vulnerable population will be difficult to control so that evacuation will be necessary for communities in the path of the active lava followed by post event population infrastructural societal and community replacement and or relocation there is a pressing need to set up a response chain that bridges scientists and responders during an effusive crisis to allow near real time delivery of globally standard products for a timely and adequate humanitarian response in this chain the scientific research groups investigating lava remote sensing and modelling need to provide products that are both useful to and trusted by the crisis response community requirements for these products include a formats that can be immediately integrated into a crisis management procedure and b in an agreed and stable standard a review of current capability reveals that we are at a point where the community can provide such a response as is the aim of the red seed risk evaluation detection and simulation during effusive eruption disasters working group this book is the first production of this group and is intended not only as a directory of current capabilities and operational service providers but also as a statement of intent and need while providing a simulation designed to demonstrate how a truly pan disciplinary response to an effusive crisis could work

observations and initial modeling of lava so2 interactions at prometheus io we present observations and initial modeling of the lava so2 interactions at the flow fronts in the prometheus region of io recent high resolution observations of prometheus reveal a compound flow field with many active flow lobes many of the flow lobes are associated with bright streaks of what is interpreted to be volatilized and recondensed so2 radiating away from the hot lava lower resolution color data show diffuse blue to violet areas also near the active flow front perhaps from active venting of so2 not clearly visible in any of the images is a single source vent for the active plume while the size of the proposed vent is probably near the limit of the resolution we expected to see radial or concentric albedo patterns or other evidence for gas and entrained particles above the flow field the lack of an obvious plume vent earlier suggestions that the prometheus type plumes may originate from the advancing flow lobes and the high resolution images showing evidence for large scale volatilization of the so2 rich substrate at prometheus encouraged us to develop a model to quantify the heat transfer between a basaltic lava flow and a substrate of so2 snow we calculate that the vaporization rate of so2 snow is 2 5 10 6 m s 1 per unit area using an estimated 5 m2 s 1 lava coverage rate from change detection images we show that the gas production rate of so2 at the flow fronts is enough to produce a resurfacing rate of 0 24 cm yr 1 at the annulus of prometheus this is much less than other estimates of resurfacing by the prometheus plume while not easily explaining the main prometheus plume our model readily accounts for the bright streaks

pitted cones and domes on mars observations in acidalia planitia and cydonia mensae using moc themis and tes data 1 domes and cones with summit pits located in acidalia planitia and cydonia mensae were studied using moc and themis images and a tes derived thermal inertia map north of 40 5 n latitude the features have a dome like morphology and south of that latitude the morphology is more cone like layering is apparent in the summit craters of fresher looking southern cones and asymmetric aprons were observed in some instances some of the northern domes also display layering in their summit craters but asymmetric aprons were not observed the northern domes can also display multiple summit pits or no summit pits at all and can occur in association with higher albedo pancake features the northern domes are higher in albedo but have apparent thermal inertias that are lower than the surrounding plains the apparent thermal inertia values of the southern cones range from values comparable to the surrounding plains to slightly lower from the tes thermal inertia map we infer that the thermal inertia values of the pitted cones are between those of basaltic fine dust and sand while those of the surrounding plains are closer to that of basaltic sand while a unique interpretation of the origin of the pitted cones is not possible with the available data we do not find compelling evidence to suggest an origin related to either basaltic volcanism or ground ice instead an origin for these features through some combination of mud volcanism and evaporite deposition around geysers and or springs is most consistent with the observations

relative dating of hawaiian lava flows using multispectral thermal infrared images a new tool for geologic mapping of young volcanic terranes the weathering of hawaiian basalts in arid and semiarid environments is accompanied by change in their thermal infrared emittance spectra the spectral differences can be measured and mapped with multispectral imaging systems the differences appear to be related to the degree of development preservation and alteration of glassy crusts the oxidation of iron and the accretion of silica rich surface veneers because the measurements are quantitative and in image format they are useful for estimating relative ages in geologic mapping of lava flows in hawaii this technique is most diagnostic for distinguishing among sparsely vegetated flows less than 1 5 ka in age

refraction seismic investigations of the northern massif central france abstract new refraction seismic data have been obtained in the northern french massif central seven shots were recorded simultaneously on two parallel nnw sse i e perpendicular to the surface structures oriented profiles allowing for inline and fan interpretations the results of a 2d interpretation reveal a large scale updoming of the moho of about 2 km compared to average western europe corresponding to the cenozoic uplift of the massif central local crustal thinning of up to 20 is restricted to the cenozoic grabens as seen at the surface beneath the volcanic area of cantal one of the largest cenozoic european shield volcanoes which was active in miocene to recent times significant crustal thickening of 2 3 km is observed velocities in the upper mantle are generally in the range of 7 9 8 1 km s beneath the limagne graben and surrounding undisturbed areas underneath the volcanic area of cezallier however velocity is reduced to about 7 7 km s corresponding to the area of crustal thickening therefore the low velocity in this zone is explained by the remainders of cooling magma chambers and heating by ascending hot plume material the cenozoic tectonic evolution of the massif central has affected crust and upper mantle in different ways oligocene rifting corresponds well to local crustal thinning but no correlation is detected to velocities of the upper mantle in contrast the volanic event occurring since miocene times is related to upper mantle velocity reduction and crustal thickening due to underplating

a theory for plastic deformation and textural evolution of olivine polycrystals seismic anisotropy in the upper mantle is due primarily to preferred orientation of olivine crystals induced by progressive deformation in order to understand better the origin of seismic anisotropy we present a simple theory for plastic deformation and textural evolution of olivine polycrystals each crystal in the aggregate is assumed to deform by intracrystalline slip on three major slip systems whose hardnesses and stress exponents are known from experiments the evolving grain orientation distribution in the aggregate is calculated by minimizing the difference between the local crystal deformation and the global aggregate deformation subject to the constraint of global strain compatibility the axial compression texture predicted by our model agrees well to first order with that determined experimentally by nicolas et al 1973 although there are significant discrepancies in the details our results suggest that the orientation texture developed during progressive plane strain deformation is a nearly unique function of the finite strain such that the crystallographic axes 100 010 and 001 are concentrated around the finite strain axes a c and b respectively a b c this result may allow the state of finite strain at depth to be estimated from observations of seismic anisotropy

radiative temperature measurements at kupaianaha lava lake kilauea volcano hawaii field spectroradiometer data in the wavelength range of 0 4 2 5 m and spectral resolution of 1 5 nm have been used to compute the radiative temperature of the surface of kupaianaha lava lake kilauea volcano hawaii two sets of observations a total of 120 spectra were made on october 12 1987 and january 23 1988 when the lava lake was in a period of active overturning the area of the surface for which temperatures were measured was 0 23 0 55 m2 two numerical models of two and three components have been used to match the measured radiant flux ratios and to describe the surface of the lava pond in terms of radiant area and temperature three stages of activity on the lake surface are identified stage 1 characterized by magma fountaining and overturning events exhibited the hottest crustal temperatures 180 572 c and the largest fractional hot areas 10 3 stage 1 average flux densities were 2 2 104 w m2 the highest recorded for the three stages of activity on either day the largest radiative area of fresh magma was 29 at 1100 c while cooling from magmatic temperatures to newly formed crust at 790 c took place in a matter of seconds stage 2 marked by rifting events between plates of crust exhibited crustal temperatures between 100 and 340 c with fractional hot areas at least an order of magnitude lower than those found for stage 1 average flux densities calculated for three examples of stage 2 activity were 5 3 103 w m2 stage 3 which was quiescent periods when the lake was covered by a thick crust dominated the activity of the lake both temporally and spatially over 90 of the time the characteristic crustal temperature of stage 3 was 80 345 c with most solutions near 200 300 c and fractional hot areas of 10 5 of the viewing area average flux densities for stage 3 were 4 9 103 w m2 for many stage 3 examples a two component model was sufficient to describe the spectral data however for almost all of the stage 1 and 2 examples and the remainder of the stage 3 examples a three component model was required these determinations of lava temperature and radiant area have relevance for satellite and airborne measurements of the thermal characteristics of active volcanoes and indicate that temporal variability of the thermal output of lava lakes occurs on the time scale of seconds to minutes

a global study of pn anisotropy beneath continents we apply a two station method of velocity determination to p n arrival times reported in the international seismological centre bulletin during the period 1964 1993 the azimuthal variation of these velocities provides evidence of p n anisotropy as large as 9 6 anisotropic regions are identified by examination of the azimuthal variation of the p n velocities within circular caps of radius 1 5 and 3 spaced at 1 intervals across the earth this procedure leads to 255 geographically distributed estimates of anisotropy the length scale over which the anisotropy varies indicates that uppermost mantle anisotropy beneath continents is controlled by regional orogenic strains and is inconsistent with a general interpretation in terms of plate scale mantle flow in addition the correlation between p n anisotropy and surficial tectonic features indicates that the anisotropy is primarily the result of the most recent deformation event we compare our p n anisotropy results fast directions and magnitudes with those of sks splitting studies to attempt to constrain the depth extent of the anisotropy this comparison suggests that in some areas where the fast directions agree sks splitting observations may arise from a layer beginning just below the moho whereas in others it must have a deeper origin in these areas the relative magnitudes of the measured anisotropy are consistent with a thin 100 150 km lithospheric layer of anisotropy

features of lava lake filling and draining and their implications for eruption dynamics lava lakes experience filling circulation and often drainage depending upon the style of activity and location of the vent features formed by these processes have proved difficult to document due to dangerous conditions during the eruption inaccessibility and destruction of features during lake drainage k lauea iki lava lake k lauea hawai i preserves many such features because lava ponded in a pre existing crater adjacent to the vent and eventually filled to the level of and interacted with the vent and lava fountains during repeated episodes a cyclic pattern of lake filling to above vent level followed by draining back to vent level preserved features associated with both filling and draining field investigations permit us to describe the characteristic features associated with lava lakes on length scales ranging from centimeters to hundreds of meters in a fashion analogous to descriptions of lava flows multiple vertical rinds of lava coating the lake walls formed during filling as the lake deepened and lava solidified against vertical faces drainage of the lake resulted in uneven formation of roughly horizontal lava shelves on the lakeward edge of the vertical rinds the shelves correlate with stable staggered lake stands shelves either formed as broken relict slabs of lake crust that solidified in contact with the wall or by accumulation accretion and widening at the lake surface in a dynamic lateral flow regime thin upper lava shelves reflect an initially dynamic environment in which rapid lake lowering was replaced by slower and more staggered drainage with the formation of thicker more laterally continuous shelves at all lava lakes experiencing stages of filling and draining these processes may occur and result in the formation of similar sets of features

a komatiite analog to potential ultramafic materials on io ultramafic volcanism is one model which fits currently available data for some eruptions on jupiter s moon io assuming that such activity does occur it is possible to apply komatiitic analogs and terrestrial ultramafic flow models to io we have used such analogs and models along with galileo data to investigate the nature of ultramafic materials that are assumed to occur at high temperature hotspots on io the unusual komatiites of the commondale greenstone belt south africa are consistent with available galileo data on the temperatures and composition of potential ionian ultramafic materials these komatiites had high sio2 50 wt and mgo 31 wt contents which are inferred to result in high liquidus temperatures 1610 c low dynamic viscosities 0 2 pa s and low densities 2680 kg m and in the crystallization of orthopyroxene phenocrysts not found in other komatiites we used the commondale composition and the model of williams et al 1998 to investigate potential ultramafic flow behavior on io compared to their terrestrial counterparts ionian ultramafic lava flows would have less potential for turbulent flow lower maximum emplacement distances and lower thermal erosion rates because of the lower ionian gravity and their more silicic composition shallow 10 m thermal erosion channels seem less likely on io unless 1 lavas are superheated 2 lavas are more ultramafic 31 mgo 3 eruption rates eruption durations or flow volumes are high or 4 substrates are unconsolidated partly consolidated or volatile rich and thus more erodable if volatiles are present in the lava as vesicles then increasing lava vesicularity results in decreasing maximum emplacement distances decreasing thermal erosion rates and decreasing erosion channel depths for given eruption durations

lateral variations of pn wave velocity in northwestern italy the pn arrival times recorded from seismic networks operating throughout northwestern italy and surrounding regions were inverted to map the structural variations of the uppermost mantle over the area and to estimate the crustal static delays at each station by means of careful data selection a quality data set was obtained removing statistical outliers poorly recorded events and scantily sampled stations moreover synthetic data were used to evaluate the resolution power of the available data set and the adopted iterative inversion technique the agreement between synthetic and calculated models is more satisfactory where the path coverage of the rays is quite complete a low velocity zone was found beneath the western and northwestern side of the alps it may be related to an increase in depth of the moho as supported by other geophysical data the high pn velocities found in the eastern side of the western alps indicate the presence of high velocity and density lower crust rocks of the so called ivrea body the high pn velocity underlying the ligurian sea could be related to a high velocity structure existing in the upper mantle the shape of the velocity anomalies matches not only the major tectonic features but also the bouguer gravity anomalies of the area

numerical modeling of lava regolith heat transfer on the moon and implications for the preservation of implanted volatiles 1 we have performed a series of numerical simulations of heat transfer between lunar lava flows and the underlying regolith with the goal of determining the depths in the substrate beneath which implanted extralunar volatiles would survive outgassing by the downward propagating heat pulse exogenous materials of interest include solar wind and solar flare particles and the cosmogenic products of galactic cosmic ray gcr particles emplaced in the lunar regolith early in solar system history extraction and analysis during future lunar missions would yield information about the evolution of the sun and inner solar system environment particles implanted in regolith deposits may be protected from gardening and saturation if buried by a lava flow but must be sufficiently deep in the regolith to survive the consequent heating our simulations include detailed treatments of lava and regolith thermophysical properties thermal conductivity specific heat capacity which vary widely over the temperature range relevant to lava flows in the lunar environment 200 to 1500 k simulations adopting temperature dependent properties together with a treatment of latent heat of lava crystallization indicate that implanted volatiles would be fully preserved at depths greater than 20 cm beneath 1 m thick lava flows 60 deeper than predicted by simulations employing constant properties these results highlight the necessity of appropriate prescription of material properties consideration of the range of lunar lava flow thicknesses allows us to determine the range of depths from which pristine samples may be recovered

thermal infrared spectral character of hawaiian basaltic glasses multispectral thermal infrared remote sensing of hawaiian basalts younger than a few hundred years old can be used to discriminate basalt lava flows that are difficult to map with visible and near infrared images laboratory thermal infrared reflectance spectra of the exposed surfaces of mauna loa and kilauea basalts show systematic changes with age fresh glass collected from active lava flows shows a broad reflectance feature between 8 and 12 m indicative of a strong degree of disorder after a few weeks of exposure to the environment the spectra of the top surfaces of samples begin to exhibit spectral features suggestive of ordering into silicate chainlike and sheetlike units with progressive aging features of apparent sheetlike structures become the preferred mode as soon as 3 years after emplacement a silica rich coating can begin to develop on the surfaces adding another spectral feature to the laboratory reflectance measurements the thermal infrared multispectral scanner tims with six channels between 8 and 12 m has been used to take advantage of these spectral differences in basalts to map flows in hawaii remote sensing instruments with increased spectral resolution in the thermal infrared would aid in further discrimination of basalts with more subtle differences in age quenching and weathering histories

the role of porosity in thermal inertia variations on basaltic lavas abstract thermal inertia is inversely proportional to porosity for hawaiian basalts extreme porosities 80 are required if the observed low thermal inertias on the martian shield volcanoes are the result of pristine lava flow surface properties such volcanic surfaces are anticipated to have a short lifetime in the martian environment and an aeolian origin appears to be the most likely interpretation of the thermal measurements on mars

mapping the distribution of vesicular textures on silicic lavas using the thermal infrared multispectral scanner the thermal infrared multispectral scanner tims is a potentially powerful tool for mapping detailed chemical variations in silicic lava flows which in turn could expand knowledge of pre eruption compositional gradients and mixing processes in silicic magma bodies however thermal infrared image data may be greatly influenced by the range of vesicular textures which occur on silicic flows to investigate the effect of vesicularity on tims imagery independent of chemical variations we studied little glass mountain at the medicine lake volcano of northern california a large rhyolitic flow of uniform composition but textural heterogeneity the imagery was recalibrated so that the digital number values for a lake in the scene matched a calculated ideal spectrum for water tims spectra for the lava show useful differences in coarsely and finely vesicular pumice data particularly in tims bands 3 and 4 images generated by ratioing these bands accurately map out those areas known from field studies to be coarsely vesicular pumice these texture related emissivity variations are probably due to the larger vesicles being relatively deeper and separated by smaller septa leaving less smooth glass available to give the characteristic emission of the lava in studies of inaccessible lava flows as on mars areas of coarsely vesicular pumice must be identified and avoided before chemical variations can be interpreted remotely determined distributions of vesicular and glassy textures can also be related to the volatile contents and potential hazards associated with the emplacement of silicic lava flows on earth

trace metal fluxes to the ocean the importance of high standing oceanic islands 1 trace metal analyses of river channel sediments from taiwan and new zealand provide new estimates for global budgets of trace metal inputs to the world s ocean mean trace metal values from high standing island river sediments are less chemically weathered than observed in sediments collected from large continental rivers and more similar to upper continental crust ucc values our mean values for cu mn co and pb are similar to the ucc zn ni and cr are more highly enriched compared to ucc values our new estimates of global trace metal budgets combining previous estimates with these new data show that geochemical fluxes of metals into the oceans are lower than previous estimates ranging from 9 lower for co to 24 lower for cu

use of clumped isotope thermometry to constrain the crystallization temperature of diagenetic calcite we describe an approach to estimating the crystallization temperatures of diagenetic calcites using clumped isotope thermometry a paleothermometer based on the 13 c 18 o bond enrichment in carbonates application of this thermometer to calcified gastropod shells and calcite cements in an early eocene limestone from the colorado plateau reveals a record of calcite precipitation and replacement at temperatures varying from 14 to 123 c the early eocene host sediments were never deeply buried but they experienced a significant thermal pulse associated with the emplacement of a late miocene basalt flow the combination of independent constraints on thermal history with clumped isotope thermometry petrographic including cathodoluminescence observations and oxygen isotopic data provides an improved basis for estimation of the temperature and timing of diagenetic events and fluid sources the petrography and calcite 18 o values taken alone suggest that the aragonite to calcite transformation of gastropod shell material occurred simultaneously with early formation of cements and lithification of the matrix in the same sample however addition of clumped isotope thermometry demonstrates that this phase transformation of shell material occurred at temperatures of 94 123 c in a highly rock buffered microenvironment i e with the isotopic composition of fluid buffered by coexisting carbonate millions of years after lithification of the matrix and formation of initial low temperature 14 19 c calcite cements within shell body cavities clumped isotope temperatures in excess of reasonable earth surface conditions recorded by later formed cements demand that cement growth occurred in association with the lava emplacement our results illustrate the potential for clumped isotope thermometry to constrain conditions of diagenesis and guide interpretations that would not be possible on the basis of conventional stable isotopic and petrographic data alone and demonstrate how petrographic characterization of clumped isotope thermometry samples can benefit paleoclimate studies

measurements of the cooling at the base of pahoehoe flows temperature measurements were made in august and october of 1994 at the base of small active pahoehoe flow lobes on kilauea volcano hawai i surprisingly after about 2 minutes of initial rapid cooling the basal temperatures increased significantly before resuming a slower cooling trend the recorded temperatures are also warmer than those predicted by simple conductive cooling models preliminary calculations suggest that these observations can be explained by the insulating effect of microfractures in the underlying flow and the kinetics of crystallization the initial cooling may be rapid enough to super cool the lava such that the onset of crystallization is delayed when crystallization does begin it is initially rapid releasing sufficient latent heat to re heat the lava as the system equilibrates the flow resumes a slow cooling

io s loki patera modeling of three brightening events in 2013 2016 abstract loki patera is one of the most dramatically time variable volcanic features on io exhibiting episodic brightening events every 1 3 years that may produce over 15 of io s global heat flow we observed three such brightening events with adaptive optics imaging at the keck and gemini n telescopes over the course of 70 nights of observation in 2013 2016 the high cadence and multi wavelength nature of the observations provides constraints on models for activity at loki patera the matson et al 2006 model for loki patera as an overturning basaltic magma sea is adapted to fit the observations of all three events in particular we adjust the details of the overturn progression and modify the lava thermal properties to include dependencies on temperature and porosity to improve the fit to the data the preferred models find overturn front propagation velocities of 1 2 1 7 km day corresponding to resurfacing rates of 1500 2200 m2 s the time intervals of 440 540 days between successive events are roughly consistent with the 540 day period calculated by rathbun et al 2002 for events prior to 2001 the best coverage was obtained for the 2016 brightening model fits to this event require a lava bulk thermal conductivity of 0 55 0 75 w m k with the best fit obtained for a value of 0 7 w m k and an average porosity that decreases during cooling for all three events the overturn front appears to propagate around the patera in the clockwise direction opposite to what has been inferred for past brightening events there is evidence that the overturn may be more complex than a single propagating wave perhaps involving multiple simultaneous resurfacing waves as well as portions of the patera that are active even after the nominal bright phase has ended the measured intensities are anomalously low when loki patera is viewed at high emission angles suggestive of topographic shadowing due to a raised area or the edge of the depression in which the magma sea resides

emplacement and inflation of pahoehoe sheet flows observations and measurements of active lava flows on kilauea volcano hawaii inflated pahoehoe sheet flows have a distinctive horizontal upper surface which can be several hundred meters across and are bounded by steep monoclinal uplifts the inflated sheet flows we studied ranged from 1 to 5 m in thickness but initially propagated as thin sheets of fluid pahoehoe lava generally 20 30 cm thick individual lobes originated at outbreaks from the inflated front of a prior sheet flow lobe and initially moved rapidly away from their source velocities slowed greatly within hours due to radial spreading and to depletion of lava stored within the source flow as the outward flow velocity decreases cooling promotes rapid crustal growth at first the crust behaves plastically as pahoehoe toes form after the crust attains a thickness of 2 5 cm it behaves more rigidly and develops enough strength to retain incoming lava thus increasing the hydrostatic head at the flow front the increased hydrostatic pressure is distributed evenly through the liquid lava core of the flow resulting in uniform uplift of the entire sheet flow lobe initial uplift rates are rapid flows thicken to 1 m in 1 2 hours but rates decline sharply as crustal thickness increases and as outbreaks occur from the margins of the inflating lobe one flow reached a final thickness of nearly 4 m after 350 hr inflation data define power law curves whereas crustal cooling follows square root of time relationships the combination of data can be used to construct simple models of inflated sheet flows as the flow advances preferred pathways develop in the older portions of the liquid cored flow these pathways can evolve into lavatube systems within a few weeks formation of lava tubes results in highly efficient delivery of lava at velocities of several kilometers per hour to a flow front that may be moving 1 2 orders of magnitude slower if advance of the sheet flow is terminated the tube remains filled with lava that crystallizes in situ rather than draining to form the cave like lava tubes commonly associated with pahoehoe flows inflated sheet flows from kilauea and mauna loa are morphologically similar to some thick icelandic and submarine sheet flows suggesting a similar mechanism of emplacement the planar sheet like geometry of flood basalt flows may also result from inflation of sequentially emplaced flow lobes rather than nearly instantaneous emplacement as literal floods of lava

limited role for thermal erosion by turbulent lava in proximal athabasca valles mars the athabasca valles flood lava is among the most recent 50 ma and best preserved effusive lava flows on mars and was probably emplaced turbulently the williams et al 2005 model of thermal erosion by lava has been applied to what we term proximal athabasca the 75 km long upstream portion of athabasca valles for emplacement volumes of 5000 and 7500 km3 and average flow thicknesses of 20 and 30 m the duration of the eruption varies between 11 and 37 days the erosion of the lava flow substrate is investigated for three eruption temperatures 1270 c 1260 c and 1250 c and volatile contents equivalent to 0 65 vol bubbles the largest erosion depths of 3 8 7 5 m are at the lava source for 20 m thick and bubble free flows that erupted at their liquidus temperature 1270 c a substrate containing 25 vol ice leads to maximum erosion a lava temperature 20 c below liquidus reduces erosion depths by a factor of 2 2 if flow viscosity increases with increasing bubble content in the lava the presence of 30 50 vol bubbles leads to erosion depths lower than those relative to bubble free lava by a factor of 2 4 the presence of 25 vol ice in the substrate increases erosion depths by a factor of 1 3 nevertheless modeled erosion depths consistent with the emplacement volume and flow duration constraints are far less than the depth of the channel 35 100 m we conclude that thermal erosion does not appear to have had a major role in excavating athabasca valles

propagation anomalies in northwestern italy by inversion of teleseismic residuals travel time residuals of teleseismic p waves were analysed in order to elucidate the crust upper mantle structure in northwestern italy and the western alpine arc using digital data obtained from both fixed seismograph networks operating in nw italy notably liguria piedmont and temporary arrays with the aid of cross correlation techniques reliable travel time residuals were calculated which were then inverted to obtain models of propagation anomalies the reliability of the inversion procedure was tested using synthetic data the model thus obtained appears to be stable and shows strong lateral heterogeneities at a litho asthenospheric level in particular it confirms the high velocity contrast caused by the ivrea body in the shallower layers and the presence of alpine roots reaching down to at least 200 km a statistical analysis performed on the propagation times of rays crossing the resulting four layered model reveals rms below 0 1 s

mapping high pn velocity beneath the colorado plateau constrains uplift models the massive international seismological centre data set of the past 20 years and the two station method are used to determine pn velocities in the mantle lid beneath the colorado plateau in this method the event is located at distances where pn is the first arrival 2 16 and the path is in or very near the azimuth of the two station pair and crosses the plateau this technique to a large extent minimizes the hypocenter mislocation effect and possible errors due to variations in the crustal structure near the source since only the difference in travel times at the two stations is used however this technique has a few underlying assumptions and possible sources of errors such as the quality of the pn data base and station delays caused by varying crustal structure that require an extremely careful application of the method a detailed study of the source of errors and a methodology of selection of the data are presented application of this method to the colorado plateau using all possible two station pairs from 53 stations located within or along the margin of the plateau yields an average high pn velocity of 8 12 0 09 km s this value is considerably larger than the average value of 7 83 km s based on available but very limited seismic refraction profiles but is remarkably similar to the average value of 8 1 km s for the relatively stable midcontinent region our new pn velocity for the colorado plateau eliminates the paradox in the literature that emphasizes the rather close similarity between average pn velocities beneath the colorado plateau and the basin and range province while their tectonic and magmatic cenozoic history is dramatically different previous models for the structure and evolution of the plateau have used the low pn velocity as an important constraint on density and thermal state of the lithosphere hence such models should be reexamined on the basis of this new uppermost mantle pn velocity determination there are two main models that have been proposed to explain the 2 km uplift of the colorado plateau one is based on a combination of thermal thinning of the lithosphere and crustal thickening and the other involves a combination of the delamination of the subducted subhorizontal farallon oceanic plate from the overriding north american plate and crustal thickening we show that the delamination model is more readily consistent not only with our velocity determination and the elevation of the plateau but also with varied geological observations reported in the literature that concern the cenozoic evolution of western north america

thermal mechanical modeling of cooling history and fracture development in inflationary basalt lava flows abstract thermal mechanical analyses of isotherms in low volume basalt flows having a range of aspect ratios agree with inferred isotherm patterns deduced from cooling fracture patterns in field examples on the eastern snake river plain idaho and highlight the caveats of analytical models of sheet flow cooling when considering low volume flows our field observations show that low volume lava flows have low aspect ratios width divided by thickness typically

constraints on cooling and degassing of pumice during plinian volcanic eruptions based on model calculations during explosive volcanic eruptions pumice clasts are transported into the atmosphere and their thermal and degassing histories determine how much volatiles are lost syneruptively this process is studied by combining a steady state eruption column model with a model for the cooling and degassing of pumice in the model we investigate the impact of various parameters e g mass eruption rate biot number eruption temperature and geometry on the cooling and degassing of pumice and find that the biot number and the eruption temperature are the most influential during typical risetimes of pumices inside eruption columns 200 300 s those smaller than 0 5 cm in diameter are found to lose little of their volatiles syneruptively and remain in thermal equilibrium with the plume in the case of larger pumices the ratio of the cooling timescale r p 2 where rp is the pumice radius and is the thermal diffusivity to the degassing timescale w 2 d 0 where w is the vesicle wall thickness and d 0 is the species diffusivity controls degassing if this ratio is larger than 50 degassing of the pumice will be nearly complete with a few percent of the original volatiles left in the outer rind of the pumice for a ratio of less than 0 1 no volatiles escape the pumice typical values for thermal and species diffusivity as well as vesicle wall thickness indicate that in the case of degassing cl and h 2 o the transition from 0 1 to 50 will be in the pumice size range of 1 10 cm here volatile concentration is a strong function of position inside the pumice suggesting that values for volatile contents measured on picked matrix glass from crushed samples may not be the best way to estimate the volatile content of matrix glass a number which is frequently used to estimate the total volatile input of volcanic eruptions this effect is even more pronounced when looking at hydrogen isotopic fractionation with larger pumices tending toward d 80

spatter stability constraining accumulation rates and temperature conditions with experimental bomb morphology we have developed the first experimental methodology to create a volcanic spatter pile using molten basalt this method permits reproduction of thermal conditions that yield the wide variety of spatter morphologies observed in nature the morphology of the clasts is most strongly controlled by the time the clast spends above the glass transition temperature which is in turn affected by the rate of accumulation and cooling of the deposit also spatter piles that remain hotter over longer durations experience increased fusion between clasts less void space between clasts and generally larger aspect ratios our experimental method successfully replicated natural microcrystal textures rheology and clast size work is still therefore required to achieve realistic vesicle distribution and deposit void space based on presented experimental work we estimate emplacement conditions of southern idaho spatter vents to have been 850 900 c with eruption temperatures closer to 1000 1100 c the rapid decrease from eruption temperature to effective emplacement temperature is the result of clast flight as well as equilibrating with the cooler surrounding material the morphology of the natural clasts matches experiments that have accumulation rates of 2 5 4 5 m h which also is consistent with the few measurements made at active eruptions finally we provide a constraint on the temperatures and accumulation rates that can lead to the construction of fused spatter features as well as provide the steps for future experiments to investigate other aspects such as compression impact and larger sizes of spatter formation by adapting our methodology

a relation between lava discharge rate thermal insulation and flow area set using lidar data 1 a simple linear relation can be used to link time averaged discharge rate tadr and lava flow area a the relation applies to given insulation conditions as described by the characteristic flow surface temperature te and will vary from case to case depending on rheological and topographic influences on flow spreading most flows have insulation conditions that change through time modifying the relationship between tadr and area as insulation conditions evolve using lidar data we can define tadr the flow area that the discharge feeds and te allowing generation of a case specific relation to convert satellite data derived flow areas to tadr for etna s 2006 lava flow field we obtain a relation whereby tadr 5 6 10 6 a for well insulated conditions te 100 c and tadr 1 5 10 4 a for poorly insulated conditions te 600 c

a reassessment of the emplacement and erosional potential of turbulent low viscosity lavas on the moon we have reevaluated the role of thermal erosion by low viscosity lunar lavas as a mechanism for the formation of the lunar sinuous rules we have adapted the model of williams et al 1998 and used the compositions of an apollo 12 basalt and a terrestrial komatiitic basalt to investigate the compositional and environmental effects on the flow behavior of low viscosity lavas on the moon and the earth our model predicts that lunar lavas could have erupted as turbulent flows that were capable of flowing hundreds of kilometers on a sufficiently flat unobstructed substrate these results are consistent with previous studies modeling of lavas over a substrate of the same composition shows that thermal erosion rates would have been low 10 cm d 1 as a result long duration eruptions approximately months to years would have been required to incise deep tens to hundreds of meters channels partial melting and mechanical removal of the substrate a mechanism suggested by hulme 1973 to enhance erosion only slightly increases thermal erosion rates other factors such as higher flow rates or lava superheating could have produced deep rules by thermal erosion during shorter duration eruptions a superheated lunar lava not only would have had a higher erosion rate 40 cm d 1 but also would have remained uncrusted for tens of kilometers which is consistent with the open channel morphology of most sinuous rules for lunar lavas with large volatile i e vesicle contents the presence of vesicles would have tended to increase viscosity at low strain rates resulting in shorter turbulent flow distances lower thermal erosion rates and thus shallower erosion channel depths for given eruption durations

viscoelastic thermal stress in cooling basalt flows polygonal fractures commonly form in basaltic lava flows during cooling of the molten lava a transient thermomechanical model is constructed to apply the temperature history in a cooling basalt flow as a boundary condition in a mechanical model of viscoelastic stress in cooling basalt nearly all past modeling of stress in cooling basalt has treated the material as purely elastic below a well defined glass transition temperature an important distinction of this model is that the viscosity of the material is allowed to be temperature dependent and the effects of viscoelastic relaxation are included on the basis of this research faster cooling and higher strain rates lead to markedly higher effective glass transition temperatures consistent with previous observations by researchers the total thermal strain in the system is constant yet the percentage of strain that is not relaxed by viscous flow varies with position relative to the cooling surface this elastic strain partition is the portion of the total strain that creates stress in the material and it is directly dependent on cooling rate thus a genetic link between cooling rate and magnitude of stress is identified

two dimensional anisotropic tomography of lithosphere beneath france using regional arrival times this paper presents an anisotropic tomography study beneath france travel times of pn waves were inverted in order to simultaneously retrieve lateral variations of both seismic velocity and seismic anisotropy within the subcrustal lithosphere the parameterization scheme uses a classical approach the subcrustal lithosphere is modeled using a two dimensional horizontal grid of nonoverlapping cells 0 5 0 5 each cell being characterized by an average slowness the isotropic term of slowness and an azimuthal anisotropy term the velocity and anisotropy perturbation amplitudes do not exceed 2 in the paris basin and in the aquitaine basin more significant perturbations can be observed in the tectonically active regions the images of the modeled subcrustal lithosphere emphasize the not recently reactivated areas of the hercynian range by high velocities 2 to 5 from the armorican massif to the west part of the massif central and from the northeast massif central to the vosges mountains and black forest mountains regions these areas are also characterized by an anisotropy of 2 5 with the fast directions parallel to the hercynian structures low velocities 2 to 4 are observed in the neogene volcanic areas of the massif central which confirm the results obtained from a short wavelength teleseismic tomography these low velocities are most probably caused by high temperatures another low velocity region 3 to 4 roughly follows the alpine belt and is interpreted as the image of the crustal root of the range in the inner part of the alps fast pn directions lie perpendicular to the belt with an anisotropy amplitude of 3 while in the outer part they are parallel to the range with amplitudes ranging from 3 to 5 to the east and from 2 to 3 to the west the anisotropy feature of the pyrenean range follows a simple pattern one observes high amplitudes up to 5 and a fast direction which is parallel with the strike of the range the observations show seismic anisotropy for most of the tectonically active zones across france recently formed structural units and old hercynian structures as well and the amplitudes reach the same order as the amplitudes of the velocity perturbations

numerical experiments on the dynamics of channelised lava flows at mount cameroon volcano with the flowgo thermo rheological model abstract as for many other effusive volcanoes mount cameroon mc is a volcano for which only limited information exists on lava flow properties and emplacement dynamics for recent eruptions this study provides new quantitative constraints for rheological and dynamic characteristics of lava flow effusion for the 1982 and 2000 eruptive events used to calibrate the flowgo thermo rheological model for these lava flows the flowgo 1 d physical model is used to simulate down flow evolution of the geometry and rheology of channel contained cooling limited lava flows morphometric data from historical lava flows were acquired from the field e g channel geometry levee and background slope in order to estimate lava yield strength velocity and effusion rate lava density and viscosity were also estimated from compositional data and laboratory methods to account for uncertainty in the input rheological and geometrical data three end member scenarios were used to bracket the potential range in lava channel initial dimension initial lava temperature and phenocryst content for each of these scenarios two crustal growth models were used one assuming strong insulation due to lava flow surface crusting the other a much lower crusting rate twelve numerical simulations were made per flow and the results were compared against the channel geometry microlite content yield strength and viscosity estimates from field and laboratory investigations best fit models where obtained for both the 1982 and 2000 lava flows using a low rate of surface crusting high initial temperature and low phenocryst content model predicted lengths were within 5 of the actual lengths modelled mean effusion rates for the 1982 52 64 m 3 s 1 and 2000 10 m 3 s 1 flows closely matched field data derived estimations 26 68 and 9 5 m 3 s 1 respectively flowgo model results are highly sensitive to initial channel dimensions phenocryst content and the flowgo model is unable to reproduce the observed rapid near vent increase in microlite content near vent rapid cooling induced by lava degassing is yet to be fully documented and integrated into numerical lava flow models such as flowgo but obtained results already provide valuable insights for lava flow hazard assessment at mount cameroon

emplacement history and inflation evidence of a long basaltic lava flow located in southern payenia volcanic province argentina abstract the el corcovo lava flow from the huanul shield volcano in the southern mendoza province central western argentina traveled a distance of 70 km and covered a minimum area of 415 km 2 the flow emplacement was controlled both by extrinsic e g topography and intrinsic e g lava supply rate lava physicochemical characteristics factors the distal portion of the lava flow reached the colorado river valley in la pampa province where it spread and then was confined by earlier river channels cross sections through the flow surveyed at several localities show two vesicular layers surrounding a dense central section where vesicles are absent or clustered in sheet shaped and cylindrical shaped structures lavas of the el corcovo flow are alkaline basalts with low values of viscosity the morphological and structural characteristics of the flow and the presence of landforms associated with lava accumulation are the evidence of inflation this process involved the formation of a tabular sheet flow up to 4 m of thick with a large areal extent in the proximal sectors while at terminal sectors frontal lobes reached inflation values up to 10 m the numerous swelling structures present at these portions of the flow suggest the movement of lava in lava tubes we propose that this aspect and the low viscosity of the lava allowed the flow travel to a great distance on a gentle slope relief

the lithosphere in central europe seismological and petrological aspects abstract the lithosphere thickness in the variscan belt of central europe varies between about 60 and 150 km with typical values of 100 120 km our estimates derived from directionally independent representative p wave residuals are in good agreement with magnetotelluric determinations of a layer with increased conductivity in the upper mantle the large scale anisotropies of the subcrustal lithosphere beneath four seismological stations determined from spatial variations of relative p residuals vary between 6 5 and 15 2 for p velocities the s wave anisotropies determined from sks polarizations vary between 2 2 and 6 7 these values are in reasonable agreement with the anisotropies of peridotites determined in laboratory systematic spatial variations of the directional terms of relative residuals in dependence on azimuths and incidence angles suggest the existence of large dipping anisotropic structures in the subcrustal lithosphere the residual patterns at most stations in the saxothuringicum rhenohercynicum and in the massif central imply northwesterly orientated dips of the anisotropic structures while stations in the moldanubicum the alpine foredeep and most of the alps north of the insubric line suggest southeasterly orientated dips in our interpretation the dipping anisotropic structures may represent paleosubductions which retain olivine preferred orientations originating from an ancient oceanic lithosphere the variscides of central europe may thus represent a collision zone characterized by two systems of paleosubductions divergent relative to the suture between the moldanubicum and the saxothuringicum

a deterministic methodology for prediction of fracture distribution in basaltic multiflows the fracture distribution in basalt flows is a direct result of thermal processes thus basalt flows present a unique opportunity to characterize a nearly perfect deterministic system with its fundamental physical parameters fracture distribution data collected on cliff exposures of basalt flows near the idaho national engineering and environmental laboratory ineel are combined with calculations of cooling rate and temperature distribution from a finite element model to construct a predictive methodology for fracture spacing the methodology is based on an empirical power law relationship between inverse cooling rate and fracture spacing the methodology may be applied to unexposed basalt flows of approximately elliptical cross section whose thickness and width are constrained only by geophysical or borehole data if sufficient fracture data on nearby exposed flows are available the methodology aids waste remediation efforts at sites involving contaminant transport through fractured basalt such as the ineel and the hanford site in washington as well as involving transport and fluid flow through volcanic or intrusive rocks where thermal processes are responsible for fracturing

the emplacement of pahoehoe toes field observations and comparison to laboratory simulations we observed active pahoehoe lobes erupted on kilauea during may june 1996 and found a range of emplacement styles associated with variations in local effusion rate flow velocity and strain rate these emplacement styles were documented and quantified for comparison with earlier laboratory experiments at the lowest effusion rates velocities and strain rates smooth surfaced lobes were emplaced via swelling where new crust formed along an incandescent lip at the front of the lobe and the rest of the lobe was covered with a dark crust at higher effusion rates strain rates and velocities lobes were emplaced through tearing or cracking tearing was characterized by ripping of the ductile crust near the initial breakout point and most of the lobe surface was incandescent during its emplacement this mechanism was observed to generate both smooth surfaced lobes and when the lava encountered an obstacle folded lobes cracking lobes were similar to those emplaced via tearing but involved breaking of a thicker brittle crust at the initial breakout of the lobe and therefore required somewhat higher flow rates than did tearing cracking lobes typically formed ropy folds in the center of the lobe and smooth margins at the highest effusion rates strain rates and flow velocities the lava formed open channels with distinct levees the final lobe morphologies were compared to results from laboratory simulations which were designed to infer effusion rate from final flow morphology to quantitatively test the laboratory results on the scale of individual natural pahoehoe lobes there is general agreement between results from laboratory simulations and natural lavas on the scale of individual pahoehoe lobes but there are disparities between laboratory flows and lava flows on the scale of an entire pahoehoe lava flow field