neotectonic faulting and forearc sliver motion along the atirro rio sucio fault system costa rica central america the atirro rio sucio fault system forms a major northwest trending strike slip fault zone in east central costa rica we examined the kinematics and temporal evolution of this fault system through geomorphic structural and seismologic analysis this 150 km long strike slip fault zone traverses the northern flank of the paleovolcanic cordillera de talamanca and extends northwestward into the active cordillera volcanica central historical seismicity includes frequent minor swarms and occasional moderate magnitude m 5 0 6 5 damaging earthquakes field geomorphic evidence fault kinematic data and earthquake focal mechanisms are consistent in showing dextral slip along the mapped traces of northwest striking faults continuity with other transcurrent faults in northwest costa rica indicates that the atirro rio sucio fault system may form the southeastern end of a regional network of northwest trending dextral faults that accommodate margin parallel displacement of the central american forearc sliver the atirro rio sucio fault system originates within the central costa rica deformed belt inboard of the indenting cocos ridge we infer that ridge collision drives lateral escape of crustal fragments northwestward along an array of dextral central costa rica deformed belt faults including the major structures of the atirro rio sucio fault system this zone of arc parallel extrusion thus represents the root of the central american forearc sliver consistent with recent geodynamic models we propose that northwestward sliver escape along the atirro rio sucio faults is driven by rigid indentation of the aseismic cocos ridge into southern costa rica

the evolution of the western ligurian flysch units and the role of mud diapirism in ancient accretionary prisms maritime alps northwestern italy the western ligurian flysch units represent an upper cretaceous to paleogene accretionary prism overthrust onto the paleo european continental margin during the collisional stages of the alpine orogeny eocene their precollisional evolution has been reconsidered according to the data collected through extensive field mapping the study of lowgrade metamorphism and by considering more strictly the processes documented in modern convergent zones

deep magma body beneath the summit and rift zones of kilauea volcano hawaii a magnitude 7 2 earthquake in 1975 caused the south flank of kilauea volcano hawaii to move seaward in response to slippage along a deep fault since then a large part of the volcano s edifice has been adjusting to this perturbation the summit of kilauea extended at a rate of 0 26 meter per year until 1983 the south flank uplifted more than 0 5 meter and the axes of both the volcano s rift zones extended and subsided the summit continues to subside these ground surface motions have been remarkably steady and much more widespread than those caused by either recurrent inflation and deflation of the summit magma chamber or the episodic propagation of dikes into the rift zones kilauea s magmatic system is therefore probably deeper and more extensive than previously thought the summit and both rift zones may be underlain by a thick near vertical dike like magma system at a depth of 3 to 9 kilometers

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

relations among subduction parameters clues to the dynamics of the subduction process are found in the many measurable parameters of modern subduction zones based on a critical appraisal of the geophysical and geological literature 26 parameters are estimated for each of 39 modern subduction zones to isolate causal relationships among these parameters multivariate analysis is applied to this data set this analysis yields empirical quantitative relations that predict strain regime and strike slip faulting in the overriding plate maximum earthquake magnitude benioff zone length slab dip arc trench gap and maximum trench depth excellent correlation is found between length of the benioff zone and the product of convergence rate and age of the downgoing slab this relationship is consistent with the conductive heating model of molnar et al 1979 if the model is modified in one respect the rate of heating of the slab is not constant it is substantially slower during passage of the slab beneath the accretionary prism and overriding plate the structural style in the overriding plate is determined by its stress state though the stress state of overriding plates cannot be quantified one can classify each individual subduction zone into one of seven semiquantitative strain classes that form a continuum from strongly extensional class 1 back arc spreading to strongly compressional class 7 active folding and thrusting this analysis indicates that strain class is probably determined by a linear combination of convergence rate slab age and shallow slab dip interplate coupling controlled by convergence rate and slab age is an important control on strain regime and the primary control on earthquake magnitude arc parallel strike slip faulting is a common feature of convergent margins forming a forearc sliver between the strike slip fault and trench optimum conditions for the development of forearc slivers are oblique convergence a compressional environment and a continental overriding plate the primary factor controlling presence of strike slip faulting is coupling strongly oblique convergence is not required the rate of strike slip faulting is affected by both convergence obliquity and convergence rate maximum trench depth is a response to flexure of the underthrusting plate the amount of flexural deflection at the trench depends on the vertical component of slab pull force which is very sensitive to slab age and shallow slab dip shallow slab dip conforms to the cross sectional shape of the overriding plate which is controlled by width of the accretionary prism and duration of subduction deep slab dip is affected by the mantle trajectory established at shallow depth but may be modified by mantle flow much of the structural complexity of convergent margins is probably attributable to terrane juxtaposition associated with temporal changes in both forearc strike slip faulting and strain regime empirical equations relating subduction parameters can provide both a focus for future theoretical studies and a conceptual and kinematic link between plate tectonics and the geology of subduction zones

statistical tests of additional plate boundaries from plate motion inversions we have investigated the application of the f ratio test a standard statistical technique to the results of relative plate motion inversions the method tests whether the improvement in fit of the model to the data resulting from the addition of another plate to the model is greater than that expected purely by chance this approach appears to be useful in determining whether additional plate boundaries are justified we confirm previous results favoring separate north american and south american plates with a boundary located between 30 n and the equator using chase s global relative motion data we show that in addition to separate west african and somalian plates separate west indian and australian plates with a best fitting boundary between 70 e and 90 e can be resolved these results are generally consistent with the observation that the indian plate s internal deformation extends somewhat westward of the ninetyeast ridge the relative motion pole is similar to minster and jordan s and predicts the nw se compression observed in earthquake mechanisms near the ninetyeast ridge

a re evaluation of the stratigraphy and evolution of the kohistan arc sequence pakistan himalaya implications for magmatic and tectonic arc building processes new field mapping and structural data combined with published geochemical data from the kohistan arc in the nw himalaya enable a re evaluation of the arc stratigraphy key lithological units and their relationships are more clearly defined permitting the construction of a revised magmatic tectonic history for the arc the oldest units are transitional oceanic type basalts which form the basement to the subduction related sequence arc type gabbroic sheets and plutons intrude the oceanic basalts together these form the kamila amphibolite belt metasediments and basaltic lavas were deposited within an extensional basin onto the kamila amphibolite belt basement this sequence exposed across the arc forms a distinct stratigraphic unit which is formally defined here as the jaglot group sediment charged turbidity currents transported material into the basin whilst submarine eruptions contributed the basaltic component this period of extension culminated in the eruption of high mg boninites of the chalt volcanic group which overlie the rocks of the jaglot group the earliest granitoids of the kohistan batholith predate suturing and intrude the jaglot and chalt sequences at c 100 ma kohistan sutured to asia suturing being accompanied by thickening of the arc with the development of major intra arc shear zones and a penetrative regionally developed steep cleavage at c 85 ma intra arc rifting permitted the emplacement into the arc of the voluminous gabbronorites of the chilas complex which clearly intrudes the kamila amphibolite belt to the south and the jaglot group to the north the chilas complex has been regarded as part of the pre suturing juvenile arc sequence field evidence summarized here show this to be not so heat advection associated with emplacement of the complex caused amphibolite facies regional metamorphism melting of the lower arc and plutonism some of the resultant granitoid plutons were unroofed and eroded during a compressional phase at between 80 and 55 ma before emplacement of further plutons and extrusion of basaltic through to rhyolitic volcanic rocks at between 55 and 40 ma at least three phases of extension and rifting each separated by short lived phases of compression charac terized arc evolution much of the magmatism is controlled by extensional tectonics within the overriding plate of the kind commonly associated with a retreating subduction zone

growth of the lava dome and extrusion rates at soufriere hills volcano montserrat west indies 2005 2008 the third episode of lava dome growth at soufriere hills volcano began 1 august 2005 and ended 20 april 2007 volumes of the dome and talus produced were measured using a photo based method with a calibrated camera for increased accuracy the total dense rock equivalent dre volume of extruded andesite magma 306 51 mm3 was similar within error to that produced in the earlier episodes but the average extrusion rate was 5 6 0 9 m3s 1 dre higher than the previous episodes extrusion rates varied in a pulsatory manner from 0 5 m3s 1 to 20 m3s 1 on 18 may 2006 the lava dome had reached a volume of 85 mm3 dre and it was removed in its entirety during a massive dome collapse on 20 may 2006 extrusion began again almost immediately and built a dome of 170 mm3 dre with a summit height 1047 m above sea level by 4 april 2007 there were few moderate sized dome collapses 1 10 mm3 during this extrusive episode in contrast to the first episode of dome growth in 1995 8 when they were numerous the first and third episodes of dome growth showed a similar pattern of low 0 5 m3s 1 but increasing magma flux during the early stages with steady high flux after extrusion of 25 mm3

identification of structural controls in an active lava dome with high resolution dems volcan de colima mexico monitoring the topography of active lava domes is critical for detecting changes that may trigger or influence collapse or explosive activity internal dome structure and conditions are more difficult to elucidate but also play vital roles here we describe the exposure following an explosion of significant scarps in the active dome at volcan de colima mexico that are interpreted as evidence of brittle failure planes and a complex internal dome morphology in the first use of automated 3d computer vision reconstruction techniques structure from motion and multi view stereo sfm mvs on an active volcanic dome we derive high resolution surface models from oblique and archive photographs taken with a consumer camera the resulting 3d models were geo referenced using features identified in a web sourced orthoimage no ground based measurements were required in december 2010 the dome 2 14 106 m3 had a flat upper surface reflecting an overall ductile emplacement regime between then and may 2011 a period of low explosivity was accompanied by a small volume loss 0 4 105 m3 and arcuate steps appeared in the dome surface suggesting the presence of localized planes of weakness the complex array of summit scarps was exposed following a significant explosion in june 2011 and is interpreted to be the surface expression of fault planes in the dome the 1 m resolution dems indicated that the region of greatest volume loss was not coincident with the assumed location of the conduit and that heterogeneity within the dome may have been important during the june explosion

time averaged discharge rate of subaerial lava at kīlauea volcano hawai i measured from tandem x interferometry implications for magma supply and storage during 2011 2013 differencing digital elevation models dems derived from terrasar add on for digital elevation measurements tandem x synthetic aperture radar imagery provides a measurement of elevation change over time on the east rift zone ezr of k lauea volcano hawai i the effusion of lava causes changes in topography when these elevation changes are summed over the area of an active lava flow it is possible to quantify the volume of lava emplaced at the surface during the time spanned by the tandem x data a parameter that can be difficult to measure across the entirety of an 100 km2 lava flow field using ground based techniques or optical remote sensing data based on the differences between multiple tandem x derived dems collected days to weeks apart the mean dense rock equivalent time averaged discharge rate of lava at k lauea between mid 2011 and mid 2013 was approximately 2 m3 s which is about half the long term average rate over the course of k lauea s 1983 present erz eruption this result implies that there was an increase in the proportion of lava stored versus erupted a decrease in the rate of magma supply to the volcano or some combination of both during this time period in addition to constraining the time averaged discharge rate of lava and the rates of magma supply and storage topographic change maps derived from space based tandem x data provide insights into the four dimensional evolution of k lauea s erz lava flow field tandem x data are a valuable complement to other space air and ground based observations of eruptive activity at k lauea and offer great promise at locations around the world for aiding with monitoring not just volcanic eruptions but any hazardous activity that results in surface change including landslides floods earthquakes and other natural and anthropogenic processes

ni in chrome pyrope garnets a new geothermometer proton microprobe analyses of the minerals in garnet peridotite xenoliths from kimberlites show that the partitioning of ni between chrome pyrope garnet and olivine is strongly temperature t dependent the range of ni contents in olivines is small relative to that in the analyzed garnets a geothermometer therefore can be derived based only on the ni content of garnet this allows estimation of t for single cr pyrope grains such as the inclusions in diamonds if these can be assumed to have equilibrated with olivine

the multi scale influence of topography on lava flow morphology predicting lava flow pathways is important for understanding effusive eruptions and for volcanic hazard assessment one particular challenge is understanding the interplay between flow pathways and substrate topography that is often rough on a variety of scales 1 m to 10 s km to study this problem we develop a lava flow modeling framework that combines spectral analysis of substrate roughness with a new lava flow model multiflow the mutliflow model includes a multiple flow direction routing algorithm in conjunction with a thresholding function that limits the extent of the flow comparison of multiflow results with recent flows on mt etna mauna loa and kilauea suggests that multiflow accurately predicts the influence of topography on lava flow paths and morphology spectral filtering of pre flow digital elevation models allow us to quantify which topographic wavelengths are most important for influencing lava flow pathways filtering constrains the height of pre existing objects that flow overtop as a measure of flow thickness and aids in identifying divergence of flow pathways from pre eruptive topography that would result from modification of the substrate by the flow low pass filter cutoffs in the range of 10 100 m significantly improve the fit between modeled flow and real flow paths suggesting a sensitivity of lava flow paths to such mesoscale topographic wavelengths the three case studies are generalized with a parameter investigation of lava flow branching on synthetic red noise topography we demonstrate that the tendency of a flow to branch or not can be predicted on the basis of spectral characteristics of the underlying surface we end by examining spectral signatures associated with volcanic and erosional processes on the island of hawai i that may aid in forecasting lava flow pathways and deconvolving competing topographic construction and erosion in volcanic landscapes on longer timescales

edifice strength and magma transfer modulation at piton de la fournaise volcano 1 from 2003 to 2007 eruptive activity at piton de la fournaise was shown to follow cycles comprising many summit proximal eruptions and finishing by a distal eruption gps measurements evidenced striking asymmetric deformation between its western and eastern flanks horizontal displacements recorded during interdistal periods showed a characteristic amplitude at the top of the eastern flank displacements recorded at the base of the summit cone showed a bimodal distribution with low amplitudes during interdistal periods and large ones during distal eruptions to account for displacement asymmetry characteristic amplitude and large flank displacement we modeled the volcanic edifice using a drucker prager elastoplastic rheology friction angles of 15 and 30 were needed to model the displacements respectively during distal eruptions and interdistal periods this change shows that strain weakening occurred during distal events large plastic displacement that occurred in the eastern flank during distal eruptions relaxed the horizontal elastic stress accumulated during interdistal periods it triggered summit deflation horizontal magma transfer and distal flank eruption and reset the eruptive cycle our elastoplastic models also show that simple source geometries may induce large eastern flank displacements that would be explained by a complex geometry in a linear elastic edifice magma supply is often thought to control volcano s eruptive activity with surface deformation reflecting changes in magma supply rate the volcano s response being linear our results bring some evidences that on piton de la fournaise time space discretization of magma transfer may be the result of the edifice s nonlinear response rather than changes in magma supply

measuring large topographic change with insar lava thicknesses extrusion rate and subsidence rate at santiaguito volcano guatemala lava flows can produce changes in topography of the order of 10s 100s of metres a knowledge of the resulting volume change provides evidence about the dynamics of an eruption using differential insar phase delays it is possible to estimate height differences between the current topography and a digital elevation model dem this does not require a pre event sar image so it does not rely on interferometric phase remaining coherent during eruption and emplacement synthetic tests predict that we can estimate lava thicknesses of as little as 9 m given a minimum of five interferograms with suitably large orbital baseline separations in the case of continuous motion such as lava flow subsidence we invert interferometric phase simultaneously for topographic change and displacement we apply this to santiaguito volcano guatemala and measure increases in lava thickness of up to 140 m between 2000 and 2009 largely associated with activity between 2000 and 2005 we find a mean extrusion rate of 0 437 0 06 m

overriding plate deformation and variability of fore arc deformation during subduction insight from geodynamic models and application to the calabria subduction zone in nature subducting slabs and overriding plate segments bordering subduction zones are generally embedded within larger plates such large plates can impose far field boundary conditions that influence the style of subduction and overriding plate deformation here we present dynamic laboratory models of progressive subduction in three dimensional space in which the far field boundary conditions at the trailing edges of the subducting plate sp and overriding plate op are varied four configurations are presented free both plates free sp fixed op fixed and sp op fixed we investigate their impact on the kinematics and dynamics of subduction particularly focusing on overriding plate deformation the results indicate that the variation in far field boundary conditions has an influence on the slab geometry subduction partitioning and trench migration partitioning our models also indicate that in natural narrow subduction zones assuming a homogeneous overriding plate the formation of back arc basins e g tyrrhenian sea aegean sea and scotia sea is generally expected to occur at a comparable location 250 700 km from the trench irrespective of the boundary condition in addition our models indicate that the style of fore arc deformation shortening or extension is influenced by the mobility of the overriding plate through controlling the force normal to the subduction zone interface trench suction our geodynamic model that uses the sp op fixed setup is comparable to the calabria subduction zone with respect to subduction kinematics slab geometry trench curvature and accretionary configuration furthermore the model can explain back arc and fore arc extension at the calabria subduction zone since the latest middle miocene as a consequence of subduction of the narrow calabrian slab and the immobility of the subducting african plate and overriding eurasian plate this setting induced strong trench suction driving fore arc extension and forced subduction to be accommodated almost entirely by slab rollback not trenchward subducting plate motion while trench retreat was accommodated almost entirely by back arc and fore arc extension not trenchward overriding plate motion comparable to our sp op fixed model

lava flow volume and morphology from digitised contour maps a case study at mount etna sicily abstract the volume and morphology of a lava flow field can be measured by mapping and comparing the topography before and after lava emplacement contour maps are a widely available source of topographic data containing a record of the geomorphological changes at volcanoes due to lava emplacement this paper explores the use of 1 25 000 scale digitised contour maps for mapping two lava flows of diverse eruption history and morphology at mount etna sicily we find that subtracting dems created from these maps gives a good representation of volume and shape for thicker lava flows 10 m mean thickness in the case of the 1983 lava flow however thinner

using tandem x to measure pyroclastic flow source location thickness and volume application to the 3rd june 2018 eruption of fuego volcano guatemala abstract the estimation of the volume of volcanic flows during an ongoing eruption is challenging but this information is crucial for improving risk assessment and for forecasting future events although previous studies have shown the ability of tandem x satellite data to derive the thickness and the volume of lava flow fields during effusive eruptions the method has not been explored yet for pyroclastic flows using bi static interferometry we produce tandem x dem on fuego volcano guatemala to measure the significant topographic changes caused by the 3rd june 2018 eruption which destroyed the town of san miguel los lotes we estimate the volume of the pyroclastic density currents pdcs to be 15 1 4 2 106 m3 the deposits are likely to be the source of lahars during future rainy seasons we identify the main channel of deposition positive elevation changes and the source region of pyroclastic material areas of significant substrate erosion and vegetation destruction negative elevation changes our results show that the june 3rd 2018 pyroclastic flow was predominantly composed of material which had gravitationally collapsed from a location close to the vent the eroded material increased the volume of the flow bulking and likely caused the run out distance of the 2018 pdc to be larger than previous eruptions 1999 2017 this study highlights the potential of remote sensing techniques for actively monitoring topography changes in inaccessible locations and to rapidly derive deposit volumes

evolution of dike opening during the march 2011 kamoamoa fissure eruption kīlauea volcano hawai i 1 the 5 9 march 2011 kamoamoa fissure eruption along the east rift zone of k lauea volcano hawai i followed months of pronounced inflation at k lauea summit we examine dike opening during and after the eruption using a comprehensive interferometric synthetic aperture radar insar data set in combination with continuous gps data we solve for distributed dike displacements using a whole k lauea model with dilating rift zones and possibly a deep decollement modeled surface dike opening increased from nearly 1 5 m to over 2 8 m from the first day to the end of the eruption in agreement with field observations of surface fracturing surface dike opening ceased following the eruption but subsurface opening in the dike continued into may 2011 dike volumes increased from 15 to 16 to 21 million cubic meters mcm after the first day eruption end and 2 months following respectively dike shape is distinctive with a main limb plunging from the surface to 2 3 km depth in the up rift direction toward k lauea s summit and a lesser projection extending in the down rift direction toward pu u at 2 km depth volume losses beneath k lauea summit 1 7 mcm and pu u 5 6 mcm crater relative to dike plus erupted volume 18 3 mcm yield a dike to source volume ratio of 2 5 that is in the range expected for compressible magma without requiring additional sources inflation of k lauea s summit in the months before the march 2011 eruption suggests that the kamoamoa eruption resulted from overpressure of the volcano s magmatic system

volumetric change quantification of the 2010 merapi eruption using tandem x insar abstract accurate and up to date information on the time varying topography of active volcanoes is crucial for assessing material flows through these volcanic systems and along their surfaces as dome building volcanoes can change on meter scale within short time intervals in active phases conventionally used repeat pass sar interferometry insar often fails due to decorrelation effects while other methods to measure topography like aerial photogrammetry or lidar depend on the availability of the local infrastructure and on good weather conditions here we show the potential of the spaceborne earth observation mission tandem x to monitor topographic changes at active dome building volcanoes in this innovative satellite mission two nearly identical satellites fly in a close formation taking radar images of the earth s surface at the same time in bistatic mode this opens the possibility to generate digital elevation models dems of volcanoes analyzing a series of repeated dems enables assessment of volume changes caused by lava dome activity we apply this method to merapi in indonesia by analyzing the volume changes during the hazardous 2010 eruption during this eruption the old lava dome collapsed and a new lava dome was extruded which collapsed again within twelve days we show three dems derived from bistatic tandem x insar in descending orbit taken before and after the eruption we found a mean volume change of 19 10 6 m 3 in the summit area the two estimates differ by 0 4 10 6 m 3 which is 2 of the average dome volume change the tandem x based results are about 25 lower than reference data which can be explained by the shadow masks that were introduced to exclude geometrically decorrelated areas resulting in conservative volume change estimates

plate tectonics and orogenic research after 25 years a tethyan perspective abstract orogeny the process by which the earth s prominent mountain ranges are constructed has been a central topic of interest in the earth sciences since at least the end of the 18th century the recognition that strains and displacements of very considerable magnitude occur along all of the three dimensions within an orogenic belt during its evolution has grown gradually during the last two centuries emphasis on primary vertical movements dominated the ideas on the nature of orogeny during the first half of the 19th century whereas compression and consequent uplift across mountain belts were believed to be the main cause for their origin during the subsequent one hundred years or so that were spent under the dominance of the fixist contraction theory mobilist tectonicians realised that continental drift in places also required motion along the trend of orogenic belts but this view did not gain general acceptance recognition of significant strike slip motion parallel or subparallel with mountain ranges evolved independently and mostly within the fixist camp by the 1960 s presence of important motions both along and across mountain belts had become common knowledge but no theoretical basis existed to account for them all

laser ablation inductively coupled plasma mass spectrometry a new technique for the determination of trace and ultra trace elements in silicates abstract this paper describes recent work applying a laser ablation system coupled to an inductively coupled plasma mass spectrometer la icp ms for the direct analysis of solid geological materials this work demonstrates the potential of la icp ms for the determination of a wide range of petrogenetically important trace and ultra trace elements including for example ree hf ta nb th u following a routine method of sample preparation powdered geological materials have been prepared as both pressed powder disks and fused glasses both common methods of sample preparation for x ray fluorescence xrf analysis the solid materials were sampled by ablation using a pulsed nd yag laser operating at 1064 nm analyses can be produced at approximately 10 samples per hour this instrumental method has limits of detection at or close to those in chondritic meteorites and gives linear calibrations over four orders of magnitude the accuracy of the technique has been evaluated using reference materials to calibrate the instrument and treating geological survey of japan basalts jb 1a jb 2 and jb 3 as unknowns detection limits are better than routine xrf analysis and compare favourably with instrumental neutron activation analysis laser ablation overcomes the problems of sample dissolution employed in standard wet chemical techniques whilst the fused glasses provide homogeneous solid samples the fused glass technique has been applied to a wide range of reference materials from ultra basic rocks through basalts and andesites to granites as well as syenite mica schist and black shale for all of the elements commonly used to generate multi element discrimination diagrams the data obtained define straight line calibrations this method is therefore capable of analysing the complete range of silicate compositions normally encountered with a single calibration i e there is no apparent matrix effect

high resolution tandem x dem an accurate method to estimate lava flow volumes at nyamulagira volcano d r congo nyamulagira and nyiragongo are two of the most active volcanoes in africa but their eruptive histories are poorly known assessing lava flow volumes in the region remains difficult as field surveys are often impossible and available digital elevation models dems do not have adequate spatial or temporal resolutions we therefore use terrasar x add on for digital elevation measurement tandem x interferometry to produce a series of 0 15 arc sec 5 m dems from between 2011 and 2012 over these volcanoes tandem x dems have an absolute vertical accuracy of 1 6 m resulting from the comparison of elevation with gps measurements acquired around nyiragongo the difference between tandem x derived dems from before and after the 2011 2012 eruption of nyamulagira provides an accurate thickness map of the lava flow emplaced during that activity values range from 3 m along the margins to 35 m in the middle with a mean of 12 7 m the erupted volume is 305 2 36 0 106 m3 height errors on thickness depend on the land covered by the flow and range from 0 4 m in old lavas to 5 5 m in dense vegetation we also reevaluate the volume of historical eruptions at nyamulagira since 2001 from the difference between tandem x and srtm 1 arc sec dems and compare them to previous work planimetric methods used in literature are consistent with our results for short duration eruptions but largely underestimate the volume of the long lived 2011 2012 eruption our new estimates of erupted volumes suggest that the mean eruption rate and the magma supply rate were relatively constant at nyamulagira during 2001 2012 respectively 23 1 m3 s 1 and 0 9 m3 s 1

variation of cenozoic extension and volcanism across the southern sierra madre occidental volcanic province mexico the middle to late cenozoic tectonic magmatic evolution of the sierra madre occidental volcanic province south of the tropic of cancer is summarized and analyzed for the first time based on new geologic and structural work and published information in the eastern part of the study region mesa central physiographic province silicic volcanism occurred in a short lived episode culminating at ca 30 ma and was followed by crustal scale extension between 30 and 27 ma in the western part of the study area sierra madre occidental physiographic province a voluminous episode of ignimbrite volcanism at 24 21 ma was succeeded by east west extension that produced regularly spaced grabens affecting only the upper crust in the westernmost part of the study region an andesitic to rhyolitic arc formed between 17 and 12 ma was affected by crustal scale north northwest trending extensional faulting leading to the formation of the gulf of california in the mesa central the maximum extension was oriented approximately east west and amounted to 20 in the eastern sierra madre occidental physiographic province extension was only 8 and oriented approximately east west we observe that trenchward shifting of the climax of subduction volcanism and extension occurred during late oligocene early miocene and late miocene time comparison with the offshore tectonics indicates that the first two tectonic magmatic pulses coincide with periods of fast spreading at the pacific farallon boundary south of the shirley fracture zone we propose that increases in the spreading rate are related to periods of high subduction rate which in turn correspond to episodes of retreating subduction a retreating slab may have generated a flux of hotter asthenospheric material into the mantle wedge producing widespread melting at the base of the crust as well as intraarc extension in the overriding plate boundary conditions i e plate tectonics ultimately determined timing magnitude and orientation of extension whereas volcanic and tectonic styles are controlled by the internal structure of crustal blocks and by the gravitational and thermal effects of magmatism

photogrammetric monitoring of lava dome growth during the 2009 eruption of redoubt volcano abstract the 2009 eruption of redoubt volcano alaska began with a phreatic explosion on 15 march followed by a series of at least 19 explosive events and growth and destruction of at least two and likely three lava domes between 22 march and 4 april on 4 april explosive activity gave way to continuous lava effusion within the summit crater we present an analysis of post 4 april lava dome growth using an oblique photogrammetry approach that provides a safe rapid and accurate means of measuring dome growth photogrammetric analyses of oblique digital images acquired during helicopter observation flights and fixed wing volcanic gas surveys produced a series of digital elevation models dems of the lava dome from 16 april to 23 september the dems were used to calculate estimates of volume and time averaged extrusion rates and to quantify morphological changes during dome growth effusion rates ranged from a maximum of 35 m 3 s 1 during the initial two weeks to a low of 2 2 m 3 s 1 in early summer 2009 the average effusion rate from april to july was 9 5 m 3 s 1 early rapid dome growth was characterized by extrusion of blocky lava that spread laterally within the summit crater in mid to late april the volume of the dome had reached 36 10 6 m 3 roughly half of the total volume and dome growth within the summit crater began to be limited by confining crater walls to the south east and west once the dome reached the steep north sloping gorge that breaches the crater growth decreased to the south but the dome continued to inflate and extend northward down the gorge effusion slowed during 16 april 1 may but in early may the rate increased again this rate increase was accompanied by a transition to exogenous dome growth from mid may to july the effusion rate consistently declined the decrease is consistent with observations of reduced seismicity gas emission and thermal anomalies as well as declining rates of geodetic deflation or inflation these trends suggest dome growth ceased by july 2009 the volume of the dome at the end of the 2009 eruption was about 72 10 6 m 3 more than twice the estimated volume of the largest dome extruded during the 1989 1990 eruption in total the 2009 dome extends over 400 m down the glacial gorge on the north end of the crater with a total length of 1 km width of 500 m and an average thickness of 200 m

flat subduction dynamics and deformation of the south american plate insights from analog modeling received 14 june 2007 revised 13 january 2008 accepted 12 march 2008 published 21 june 2008 1 we present lithospheric scale analog models investigating how the absolute plates motion and subduction of buoyant oceanic plateaus can affect both the kinematics and the geometry of subduction possibly resulting in the appearance of flat slab segments and how it changes the overriding plate tectonic regime experiments suggest that flat subductions only occur if a large amount of a buoyant slab segment is forced into subduction by kinematic boundary conditions part of the buoyant plateau being incorporated in the steep part of the slab to balance the negative buoyancy of the dense oceanic slab slab flattening is a long term process 10 ma which requires the subduction of hundreds of kilometers of buoyant plateau the overriding plate shortening rate increases if the oceanic plateau is large enough to decrease the slab pull effect slab flattening increases the interplate friction force and results in migration of the shortening zone within the interior of the overriding plate the increase of the overriding plate topography close to the trench results from 1 the buoyancy of the plate subducting at trench and 2 the overriding plate shortening experiments are compared to the south american active margin where two major horizontal slab segments had formed since the pliocene along the south american subduction zone flat slab segments below peru and central chile nw argentina appeared at 7 ma following the beginning of buoyant slab segments subduction in northern ecuador and northern chile the process of slab flattening resulting from the carnegie and iquique ridges subductions respectively seems to be active but not completed the formation of flat slab segments below south america from the pliocene may explain the deceleration of the nazca plate trenchward velocity citation espurt n f funiciello j martinod b guillaume v regard c faccenna and s brusset 2008 flat subduction dynamics and deformation of the south american plate insights

the roughness of natural terrain a planetary and remote sensing perspective we examine the various methods and parameters in common use for quantifying and reporting surface topographic roughness it is shown that scale dependent roughness parameters are almost always required though not widely used we suggest a method of standardizing the parameters that are computed and reported so that topographic data gathered by different workers using different field techniques can be directly and easily intercompared we illustrate the proposed method by analyzing topographic data from 60 different surfaces gathered by five different groups and examine the information for common features we briefly discuss the implications of our analysis for studies of planetary surface roughness lander safety and radar remote sensing modeling and analysis

extensional compressional tectonics associated with the aegean arc comparison with the andean cordillera of south peru north bolivia extensional compressional tectonics are analysed along the convergent boundaries of the aegean arc and of the andean cordillera of south peru and north bolivia this paper deals only with deformations of miocene pliocene and pleistocene ages the prevailing regime in the aegean basin is extensional yet two compressional tectonic phases of u miocene l pliocene and of l pleistocene age have been recorded they are also known in many places all around the mediterranean thus forces that induce compressional tectonics in the aegean seem to have their origin outside the region they seem to be linked with changes in the conditions of convergence change in rate between the arabo african and eurasian plates beneath the andean cordillera of south peru north bolivia the subduction zone plunges at an angle of 25 30 as it does under the aegean arc there as in the aegean sea tectonics are predominantly extensional and several short lived compressional phases break up this extensional regime the latter may also be induced by changes of the conditions of convergence here between the nazca and south american plates beneath central peru north of the nazca ridge the subduction zone plunges at a small angle about 10 there neotectonic deformations are different in the recent they are essentially compressional or expressed by strike slip faulting geological data from the aegean and andes seem to demonstrate that the compressional deformation in the crust above subduction zones is strongly controlled by the stress conditions along convergent boundaries

magma supply and storage at kilauea volcano hawaii 1956 1983 abstract shallow crustal magma reservoirs beneath the summit of kilauea volcano and within its rift zones are linked in such a way that the magma supply to each can be estimated from the rate of ground deformation at the volcano s summit our model builds on the well documented pattern of summit inflation as magma accumulates in a shallow summit reservoir followed by deflation as magma is discharged to the surface or into the rift zones magma supply to the summit reservoir is thus proportional to summit uplift and supply to the rift zones is proportional to summit subsidence the average proportionality constant is 0 33 10 6 m 3 rad this model yields minimum supply estimates because it does not account for magma which escapes detection by moving passively through the summit reservoir or directly into the rift zones calculations suggest that magma was supplied to kilauea during july 1956 april 1983 at a minimum average rate of 7 2 10 6 m3 month roughly 35 of the net supply was extruded the rest remains stored within the volcano s east rift zone 55 and southwest rift zone 10 periods of relatively rapid supply were associated with the large kapoho eruption in 1960 and the sustained mauna ulu eruptions in 1969 1971 and 1972 1974 bursts of harmonic tremor from the mantle beneath kilauea were also unusually energetic during 1968 1975 suggesting a close link between kilauea s deep magma supply region and shallow storage reservoirs it remains unclear whether pulses in magma supply from depth give rise to corresponding increases in shallow supply or if instead unloading of a delicately balanced magma transport system during large eruptions or intrusions triggers more rapid ascent from a relatively constant mantle source

lava flow mapping and volume calculations for the 2012 2013 tolbachik kamchatka fissure eruption using bistatic tandem x insar the bistatic acquisition mode of the german tandem x radar satellite mission provides a reliable source for measuring morphological changes associated with volcanic activity we present the use of this system to measure key lava flow parameters including thickness volume runout and flow extent by using two tandem x data pairs to generate digital elevation models dems prior to and immediately following the 2012 2013 eruption of tolbachik volcano kamchatka morphometric parameters and areal distribution of the new lava flow field are determined using a cell by cell elevation difference between the two dems a total flow volume of 0 53 0 07 km3 a mean flow thickness of 14 5 m and a modal thickness of 7 8 m are calculated we use these calculated flow parameters as input to a volume limited lava flow emplacement model model simulations are able to reproduce the sw portion of the 2012 2013 tolbachik lava flow using a 75 m shuttle radar topography mission srtm dem and the 15 m tandem x derived dem with goodness of fit measures of 56 3 and 59 6 respectively based on the jaccard similarity coefficient the flow simulation done using srtm data underestimates the observed 14 4 km flow runout by over 3 km while the simulation with tandem x data overestimates flow runout by about 1 5 km performance of the lava flow modeling algorithm is highly dependent on the modal lava thickness highlighting the importance of using tandem x dems to provide precise lava flow measurements in order to constrain input parameters for numerical modeling of lava flows

petroleum systems of the simeulue fore arc basin offshore sumatra indonesia forearc basins result from plate convergence these basins are situated offshore between an outer arc high and the mainland historically these regions have not been considered important petroleum provinces partly because low heat flow may limit significant thermal hydrocarbon generation the simeulue forearc basin extends between simeulue island and northern sumatra it is a frontier shallow shelf area with few wells and no wells in the basin center therefore it is studied using geophysical data and geologic surface samples multichannel seismic data show bright spots above potential hydrocarbon reservoirs in carbonate buildups amplitude versus offset analyses indicate the presence of gas and surface geochemical prospecting suggests thermal hydrocarbon generation heat flow in the simeulue basin ranges between 37 and 74 mw m2 as deduced from one dimensional petroleum system modeling and bottom simulating reflector depths two possible source rocks eocene and lower middle miocene were assigned for three dimensional petroleum system modeling in the simeulue basin because of a similar pre miocene geologic evolution of the present day back arc and the fore arc it can be assumed that the back arc source rocks also occur in the fore arc modeling based on two heat flow scenarios 40 and 60 mw m2 reveals that oil and gas generation is possible within and below the main depocenters of the central and southern simeulue basin this study shows that deep burial 6 km 3 7 mi of source rocks can compensate for low heat flow therefore forearc basins may be more prolific for hydrocarbons than previously considered and each forearc basin should be studied carefully to evaluate its hydrocarbon potential

causes of compression and extension behind trenches abstract three factors have hampered attempts to find the causes of compression and extension within overriding plates of subduction zones 1 1 adjacent subduction zones are often similar in many subduction parameters 2 2 proposed causal variables are themselves correlated and 3 3 physical models are limited by uncertainties in physical quantities and complexity of the subduction process nevertheless modern subduction zones provide a rich source of observations that help to distinguish among models we consider ten variables or combinations of variables which may control strain regime behind trenches convergence rate vc slab age as shallow slab dip ds absolute motion voa asthenospheric flow slab pull force as vc as voa vc ds voa and vc ds as ten observations are used to test these models a key prerequisite to these tests is the semiquantitative classification of the strain regime of each modern subduction zone into one of seven strain classes from active back arc spreading class 1 through active folding and foreland thrusting class 7 statistical analysis indicates that the modern variance in strain classes is not adequately accounted for by any single or dual variable model but the combination of vc ds and either asorvoa can account for 77 of this variance r 0 88 rapid changes in strain class both with time and along a subduction zone are similarly more consistent with the 3 variable models the role of vc is also supported by maximum compression orientations parallel to vc the role of ds is also supported by the tendency for oceanic overriding plates to be more extensional than continental overriding plates two factors suggest that the vc ds as model is more likely than vc ds voa 1 1 the dominance of one or more of the variables vc as andds in other subduction relations slab length earthquake moment trench relative depth deep dip and arc trench gap and 2 2 the correlation between strain class and maximum earthquake moment which is affected by vc and as however no confident discrimination between the two 3 variable models is yet possible

the tectonic expression slab pull at continental convergent boundaries examination of five thrust belt systems developed at continental subduction boundaries suggests that they comprise two distinct groups that display pronounced and systematic differences in structural style topographic elevation denudation metamorphism postcollisional convergence and foredeep basin geometry and facies the distinctive geological features developed within each thrust belt group appear to be causally linked to the relative rates of subduction and convergence via the magnitude of horizontal compressional stress transmitted across the subduction boundary at subduction boundaries where the rate of overall plate convergence is less than the rate of subduction termed here retreating subduction boundaries the transmission of horizontal compressive stress across the plate boundary is small and regional deformation of the overriding plate is by horizontal extension the tectonic expression of these retreating subduction boundaries includes topographically low mountains little erosion or denudation low grade to no metamorphism little to no involvement of crystalline basement in shortening little to no postcollisional convergence anomalously deep foredeep basins and a protracted history of flysch deposition within the adjacent foredeep basin analysis of deflection and gravity data across three retreating subduction boundaries apennine carpathian and hellenic systems shows that subduction is driven by gravitational forces acting on dense subducted slabs at depths between about 40 and 80 km carpathians 50 and 150 km apennines and 50 and 250 km hellenides the total mass anomalies represented by the slabs are approximately 3 1012 6 1012 and 12 1012 n m respectively the slabs are partially supported by flexural stresses transmitted through the subducted lithosphere to the foreland and partially supported by dynamic viscous stresses in the asthenosphere at subduction boundaries where the rate of overall plate convergence is greater than the rate of subduction termed here advancing subduction boundaries the transmission of horizontal compressive stress across the plate boundary is large and regional deformation of the overriding plate is by horizontal shortening the tectonic expression of these advancing subduction boundaries includes topographically high mountains antithetic thrust belts large amounts of erosion and denudation exposure of high grade metamorphic rocks at the surface extensive deformation of crystalline basement to midcrustal depths protracted postcollisional convergence tens of millions of years and a protracted history of molasse deposition within the adjacent foredeep basins analysis of gravity and deflection data across two advancing subduction boundaries developed within the continental lithosphere western to eastern and southern alps and himalayas shows that the thrust sheets have been translated for great distances over the foreland lithosphere relative to the point at which the subduction forces are applied thus obscuring any flexural and gravity signals from the subducted slab however it appears that far field stresses presumably related to global plate motions drive most of the convergent motion across these subduction boundaries the concept that orogenic belts formed above retreating subduction boundaries have recognizable tectonic signatures that differ from those of orogenic belts formed above advancing subduction boundaries suggests that it may be possible to interpret the plate boundary settings in which ancient orogenic belts evolved appendix b is available with entire article on microfiche order from the american geophysical union 2000 floridaavenue n w washington d c 20009 document t92 004 2 50 payment must accompany order

satellite laser ranging and geological constraints on plate motion we have analyzed satellite laser ranging slr observed baseline rates of change and compared them with rates determined from sea floor spreading rates and directions and earthquake solutions with the number of years of observation now over six for many of the baselines the inaccuracy of determining baseline rates of change has diminished so that in some cases it is less than a few mm per year thus a direct comparison between baseline rates of change and rates of change established using geophysical information which we call geological rates is now feasible in most cases there is good agreement between the rates determined from slr and geological rates but in some cases there appear to be discrepancies these discrepancies involve many of the data for which one end of the baseline is either quincy california huahine french polynesia or simosato japan we have devised a method for looking at the discrepancies for these slr observatories which allows us to calculate the motion not modelled by the geologic information the results are discussed in terms of what is known about plate margins and other information

controls of subduction geometry location of magmatic arcs and tectonics of arc and back arc regions most variation in geometry and angle of inclination of subducted oceanic lithosphere is caused by four interdependent factors combinations of 1 rapid absolute upper plate motion toward the trench and active overriding of the subducted plate 2 rapid relative plate convergence and 3 subduction of intraplate island seamount chains aseismic ridges and oceanic plateaus anomalously low density oceanic lithosphere cause low angle subduction under conditions of low angle subduction the upper surface of the subducted plate is in contact with the base of the overlying plate the wedge of low density asthenosphere is replaced by subducted lithosphere and the width of the arc trench gap either is significantly increased or a magmatic arc is not developed within the overlying plate the fourth factor is age of the subducting lithosphere subduction of young lithosphere produces two opposing tendencies 1 low angle subduction and increased arc trench distance owing to its low density and 2 decreased arc trench distance owing to its higher temperature two factors of secondary importance contribute to variation in subduction zone geometry and arc trench distance accretion of sediment in trenches depresses the upper portion of the subducting oceanic plate and causes the trench axis to migrate seaward prolonged subduction thickens the upper plate depresses the isotherms in the subducted plate and may create a broader arc both factors increase the arc trench gap the four primary factors also control development of other tectonic elements such as regional subsidence for example the amazon basin and a portion of the cretaceous interior seaway of western united states intra arc extension for example the basin and range province foreland fold and thrust belts and laramide style tectonics

a quantitative look at the demise of a basaltic vent the death of kupaianaha kilauea volcano hawai i the kupaianaha vent the source of the 48th episode of the 1983 to present pu u o o kupaianaha eruption erupted nearly continuously from july 1986 until february 1992 this investigation documents the geophysical and geologic monitoring of the final 10 months of activity at the kupaianaha vent detailed very low frequency vlf electromagnetic profiles across the single lava tube transporting lava from the vent were used to determine the cross sectional area of the molten lava within the tube combined with measurements of lava velocity these data provide an estimate of the lava output of kupaianaha in addition lava temperatures calculated from analysis of quenched glass and bulk rock chemistry were obtained for samples taken from the tube at the same site the combined data set shows the lava flux from kupaianaha vent declining linearly from 250 000 m3 day in april 1991 to 54 000 m3 day by november 1991 during that time surface breakouts of lava from weak points along the tube occurred progressively closer to the vent consistent with declining efficiency in lava transport there were no significant changes in lava temperature or in bulk mgo content during this period another eruptive episode the 49th began uprift of kupaianaha on 8 november 1991 and erupted lava concurrently with kupaianaha for 18 days lava flux from kupaianaha decreased in response to this new episode but the response was delayed by approximately 1 day after 14 november 1991 lava velocities were no longer measurable in the tube because the lava stream beneath the skylight had crusted over however the vlf derived electrical conductances documented the decreasing flux of molten lava through the tube kupaianaha remained active but output continued to decrease until early february 1992 when the last active surface flows were seen in november 1991 we used the linearly decreasing effusion rate to accurately predict the date for the death of the kupaianaha vent the linear nature of the decline in lava tube conductance and the delayed and slow response of the waha ula tube conductances to the 49th eruptive episode led us to speculate that a the kupaianaha vent shut down because of a decrease in driving pressure and not because of a freeze up of the vent and b that pu u o o episode 49 and kupaianaha were fed nearly vertically from a source deep within the rift zone

analysis of the 2001 lava flow eruption of mt etna from three dimensional mapping 1 the 2001 etna eruption was characterized by a complex temporal evolution with the opening of seven eruptive fissures each feeding different lava flows this work describes a method adopted to obtain the three dimensional geometry of the whole lava flow field and for the reconstruction based on topographic data of the temporal evolution of the largest lava flow emitted from a vent located at 2100 m a s l preeruption and posteruption digital elevation models dem were extracted from vector contour maps comparison of the two dems and analysis of posteruption orthophotos allowed us to estimate flow area thickness and bulk volume additionally the two dimensional temporal evolution of the 2100 flow was precisely reconstructed by means of maps compiled during the eruption these data together with estimates of flow thickness allowed us to evaluate emitted lava volumes and in turn the average volumetric flow rates the analysis performed in this paper provided a total lava bulk volume of 40 1 106 m3 for the whole lava flow field most of which emitted from the 2100 vent 21 4 106 m3 the derived effusion rate trend shows an initial period of waxing flow followed by a longer period of waning flow this is in agreement not only with the few available effusion rate measurements performed during the eruption but also with the theoretical model of wadge 1981 for the temporal variation in discharge during the tapping of a pressurized source

pyroclastic density current volume estimation after the 2010 merapi volcano eruption using x band sar abstract pyroclastic density current deposits remobilized by water during periods of heavy rainfall trigger lahars volcanic mudflows that affect inhabited areas at considerable distance from volcanoes even years after an eruption here we present an innovative approach to detect and estimate the thickness and volume of pyroclastic density current pdc deposits as well as erosional versus depositional environments we use sar interferometry to compare an airborne digital surface model dsm acquired in 2004 to a post eruption 2010 dsm created using cosmo skymed satellite data to estimate the volume of 2010 merapi eruption pdc deposits along the gendol river kali gendol kg results show pdc thicknesses of up to 75 m in canyons and a volume of about 40 106 m3 mainly along kg and at distances of up to 16 km from the volcano summit this volume estimate corresponds mainly to the 2010 pyroclastic deposits along the kg material that is potentially available to produce lahars our volume estimate is approximately twice that estimated by field studies a difference we consider acceptable given the uncertainties involved in both satellite and field based methods our technique can be used to rapidly evaluate volumes of pdc deposits at active volcanoes in remote settings and where continuous activity may prevent field observations

publisher s noteoverview of the precursors and dynamics of the 2012 13 basaltic fissure eruption of tolbachik volcano kamchatka russia

tracking lava flow emplacement on the east rift zone of kīlauea hawai i with synthetic aperture radar coherence 1 lava flow mapping is both an essential component of volcano monitoring and a valuable tool for investigating lava flow behavior although maps are traditionally created through field surveys remote sensing allows an extraordinary view of active lava flows while avoiding the difficulties of mapping on location synthetic aperture radar sar imagery in particular can detect changes in a flow field by comparing two images collected at different times with sar coherence new lava flows radically alter the scattering properties of the surface making the radar signal decorrelated in sar coherence images we describe a new technique sar coherence mapping scm to map lava flows automatically from coherence images independent of look angle or satellite path we use this approach to map lava flow emplacement during the pu u kupaianaha eruption at k lauea hawai i the resulting flow maps correspond well with field mapping and better resolve the internal structure of surface flows as well as the locations of active flow paths however the scm technique is only moderately successful at mapping flows that enter vegetation which is also often decorrelated between successive sar images along with measurements of planform morphology we are able to show that the length of time a flow stays decorrelated after initial emplacement is linearly related to the flow thickness finally we use interferograms obtained after flow surfaces become correlated to show that persistent decorrelation is caused by post emplacement flow subsidence

magnetic lineations in marginal basins of the western pacific this paper separates the small oceanic basins around the western perimeter of the pacific ocean into marginal basins that have formed through back arc extension and those that apparently have not and reviews our knowledge of magnetic lineation patterns observed in possible and probable back arc basins magnetic lineations in these basins resemble lineations commonly associated with the world s mid oceanic spreading systems indicating that similar processes of crustal accretion occur in both tectonic environments in some back arc basins of the southwestern pacific magnetic lineation and other evidence suggest that back arc basins can evolve through the interaction and growth of multi plate systems because of the small time and space characteristics of back arc basins compared with the world s major oceanic spreading systems tectonic conditions favourable for the generation of back arc basins are either relaxed rapidly or easily interrupted models proposed to account for back arc basin development include a local models where back arc extension is mechanically driven by the downgoing slab and b global plate kinematics models where conditions favourable for back arc extension are governed by the motion of the overriding plate relative to the trench axis

a new method for measuring deformation on volcanoes and other natural terrains using insar persistent scatterers 1 we present here a new insar persistent scatterer ps method for analyzing episodic crustal deformation in non urban environments with application to volcanic settings our method for identifying ps pixels in a series of interferograms is based primarily on phase characteristics and finds low amplitude pixels with phase stability that are not identified by the existing amplitude based algorithm our method also uses the spatial correlation of the phases rather than a well defined phase history so that we can observe temporally variable processes e g volcanic deformation the algorithm involves removing the residual topographic component of flattened interferogram phase for each ps then unwrapping the ps phases both spatially and temporally our method finds scatterers with stable phase characteristics independent of amplitudes associated with man made objects and is applicable to areas where conventional insar fails due to complete decorrelation of the majority of scatterers yet a few stable scatterers are present

high resolution dem generation from spaceborne and terrestrial remote sensing data for improved volcano hazard assessment a case study at nevado del ruiz colombia abstract volcanoes with rugged terrain remain a challenging target for generating high resolution digital elevation models dems especially in tropical areas with frequent cloud cover using nevado del ruiz volcano as an example we combined dems from the tandem x tdx satellite mission terrestrial radar interferometry tri and structure from motion sfm to generate a new dem with 10 m spatial resolution this is the first study combining satellite radar ground based radar photography and freely available global dems to generate a high resolution dem without data gaps tdx data from ascending and descending orbits were combined to generate the base dem instead of using a raster format to fuse dems generated from different data sets with different resolutions we developed a methodology based on 3 d point clouds 1 re georeference the 5 m tri and 1 m sfm dems to the 10 m tdx dem using the iterative closest point icp algorithm to minimize the horizontal and vertical discrepancy between dems then 2 merge the multiple point clouds to generate a final dem without data gaps using an adaptive algorithm that uses two search distances to smooth the transition at the edges of different data sets we assess the new 10 m dem by comparing simulated inundation zones obtained with two volcano flow models laharz for lahars and volcflow for pyroclastic flows and find significant differences with respect to the 30 m srtm dem our laharz simulation over the new dem shows a longer lahar run out distance for pyroclastic flows the volcflow simulation over the new dem produces highly channelized flows over the steep portions of a river channel and gives a larger extent of thicker deposits compared to those obtained with the 30 m srtm dem quantitative and qualitative geomorphic analysis suggests that up to date dems with high spatial resolution 10 m or even better need to be generated to improve volcano hazard assessment for active volcanoes

terrane motion by strike slip faulting of forearc slivers forearc slivers bounded by a trench and an active strike slip fault occur in about 50 of modern subduction zones analysis of earthquake slip vectors indicates that modern sliver terranes typically migrate at rates of 1 2 cm yr tertiary transport of some forearc slivers by 1000 km or more is therefore expected as a consequence of normal subduction active arc parallel strike slip faulting occurs whenever interplate coupling is strong and convergence is somewhat oblique strike slip rate increases with greater convergence obliquity

seismic structure in southern peru evidence for a smooth contortion between flat and normal subduction of the nazca plate rapid changes in slab geometry are typically associated with fragmentation of the subducted plate however continuous curvature of the slab is also possible the transition from flat to normal subduction in southern peru is one such geometrical change the morphology of the subducted nazca plate along this transition is explored using intraslab earthquakes recorded by temporary regional seismic arrays observations of a gradual increase in slab dip coupled with a lack of any gaps or vertical offsets in the intraslab seismicity suggest warping of the slab concentrations of focal mechanisms at orientations which are indicative of slab bending are also observed along the change in slab geometry the presence of a thin ultra slow velocity layer usl atop the horizontal nazca slab is identified and located the lateral extent of this usl is coincident with the margin of the projected linear continuation of the subducting nazca ridge implying a causal relationship wherein increased hydration of the ridge results in the formation of the usl downdip waveform modelling of the 2 d structure in southern peru using a finite difference algorithm provides constraints on the velocity and geometry of the slab s seismic structure and confirms the absence of any tears in the slab the seismicity and structural evidence suggests smooth contortion of the nazca plate along the transition from flat to normal subduction the slab is estimated to have experienced 10 per cent strain in the along strike direction across this transition

the opening of the woodlark basin subduction of the woodlark spreading system and the evolution of northern melanesia since mid pliocene time abstract magnetic anomaly and seismological data define segments of active seafloor spreading and associated magnetic lineations trending ene in the woodlark basin the total opening rate has been approximately 6 cm yr for the last 1 m y spreading rates diminish by over 10 from east to west along the woodlark spreading system implying a pole of current opening 15 20 to the west commencement of seafloor spreading in the basin has apparently been time transgressive beginning prior to 3 5 m y in the east and at successively later times to the west earthquake focal mechanisms and geological evidence suggest that the land areas bounding the western end of the woodlark basin are undergoing tensional deformation we believe that eventually the woodlark basin plate boundary will propagate westward through the d entrecasteaux islands into the papuan peninsula hitherto unreported shallow seismicity associated with the northern margin of the ne trending section of the woodlark rise probably reflects partial decoupling of the woodlark and solomon basins possibly due to mechanical difficulties in subducting the young woodlark lithosphere analysis of the relative motions between the solomon indo australian and pacific plates shows that the woodlark spreading system has been subducted at high rates 10 cm yr beneath the solomon islands during the opening of the woodlark basin several tectonic and geological features limited to the region of interaction of the woodlark basin with the solomon trench and arc may be symptomatic of ridge subduction these features include high heat flow in the solomon trench which shoals to 4 km low levels of seismicity and only shallow hypocenters and voluminous eruptions of high k 2 o tio 2 olivine basalts and basaltic andesites extremely close to the trench axis this close association in space and time of an unusual volcanic suite with ridge subduction implies a strong dependence of the petrogenesis on the tectonic regime a combination of this study of the woodlark basin and the previous study of the bismarck basin taylor 1979 provides a reconstruction of the positions of the continents ocean basins and island chains in northern melanesia for mid pliocene time in accepting the existence of a solomon plate we can explain the trench like structure off the trobriand margin of new guinea the occurrence of late cenozoic calc alkaline volcanism along the papuan peninsula and the presence of intermediate depth seismicity beneath the north papuan peninsula the rapid changes in relative motions along or across the new ireland solomons chain over the past 3 5 m y may explain the spatial and temporal changes in igneous activity observed on these islands

topographic and stochastic influences on pahoehoe lava lobe emplacement a detailed understanding of p hoehoe emplacement is necessary for developing accurate models of flow field development assessing hazards and interpreting the significance of lava morphology on earth and other planetary surfaces active p hoehoe lobes on k lauea volcano hawai i were examined on 21 26 february 2006 using oblique time series stereo photogrammetry and differential global positioning system measurements during this time the local discharge rate for peripheral lava lobes was generally constant at 0 0061 0 0019 m3 s but the areal coverage rate of the lobes exhibited a periodic increase every 4 13 0 64 min this periodicity is attributed to the time required for the pressure within the liquid lava core to exceed the cooling induced strength of its margins the p hoehoe flow advanced through a series of down slope and cross slope breakouts which began as 0 2 m thick units i e toes that coalesced and inflated to become approximately meter thick lobes the lobes were thickest above the lowest points of the initial topography and above shallow to reverse facing slopes defined relative to the local flow direction the flow path was typically controlled by high standing topography with the zone directly adjacent to the final lobe margin having an average relief that was a few centimeters higher than the lava inundated region this suggests that toe scale topography can at least temporarily exert strong controls on p hoehoe flow paths by impeding the lateral spreading of the lobe observed cycles of enhanced areal spreading and inflated lobe morphology are also explored using a model that considers the statistical likelihood of sequential breakouts from active flow margins and the effects of topographic barriers

rapid low cost photogrammetry to monitor volcanic eruptions an example from mount st helens washington usa we describe a low cost application of digital photogrammetry using commercially available photogrammetric software and oblique photographs taken with an off the shelf digital camera to create sequential digital elevation models dems of a lava dome that grew during the 2004 2008 eruption of mount st helens msh volcano renewed activity at msh provided an opportunity to devise and test this method because it could be validated against other observations of this well monitored volcano the datasets consist of oblique aerial photographs snapshots taken from a helicopter using a digital single lens reflex camera twelve sets of overlapping digital images of the dome taken during 2004 2007 were used to produce dems and to calculate lava dome volumes and extrusion rates analyses of the digital images were carried out using photogrammetric software to produce three dimensional coordinates of points identified in multiple photos the evolving morphology of the dome was modeled by comparing successive dems results were validated by comparison to volume measurements derived from traditional vertical photogrammetric surveys by the us geological survey cascades volcano observatory our technique was significantly less expensive and required less time than traditional vertical photogrammetric techniques yet it consistently yielded volume estimates within 5 of the traditional method this technique provides an inexpensive rapid assessment tool for tracking lava dome growth or other topographic changes at restless volcanoes

high resolution digital elevation model from tri stereo pleiades 1 satellite imagery for lava flow volume estimates at fogo volcano resolving changes in topography through time using accurate high resolution digital elevation models dems is key to understanding active volcanic processes for the first time in a volcanic environment we utilize very high resolution tri stereo optical imagery acquired by the pleiades 1 satellite constellation and generate a 1 m resolution dem of fogo volcano cape verde the most active volcano in the eastern atlantic region point cloud density is increased by a factor of 6 5 compared to conventional stereo imagery and the number of 1 m2 pixels with no height measurements is reduced by 43 we use the dem to quantify topographic changes associated with the 2014 2015 eruption at fogo height differences between the posteruptive pleiades 1 dem and the preeruptive topography from tandem x give a lava flow volume of 45 83 0 02 106 m3 emplaced over an area of 4 8 km2 at a mean rate of 6 8 m3 s 1

volumetric characteristics of lava flows from interferometric radar and multispectral satellite data the 1995 fernandina and 1998 cerro azul eruptions in the western galapagos we have used a suite of remotely sensed data numerical lava flow modeling and field observations to determine quantitative characteristics of the 1995 fernandina and 1998 cerro azul eruptions in the western galapagos islands flank lava flow areas volumes instantaneous effusion rates and average effusion rates were all determined for these two eruptions for which only limited syn eruptive field observations are available using data from spot topsar ers 1 and ers 2 we determined that the 1995 fernandina flow covers a subaerial area of 6 5 106 m2 and has a subaerial dense rock equivalent dre volume of 42 106 m3 field observations atsr satellite data and the flowgo numerical model allow us to determine that the effusion rate declined exponentially from a high of 60 200 m3 s 1 during the first few hours to 5 m3 s 1 prior to ceasing after 73 days with a mean effusion rate of 4 16 m3 s 1 integrating the atsr derived exponentially declining effusion rate over the eruption duration produces a total subaerial submarine dre volume of between 27 and 100 106 m3 the range in values being due to differing assumptions about heat loss characteristics only values in the higher part of this range are consistent with the independently derived subaerial volume using spot topsar ers 1 and ers 2 data we determine that the 1998 cerro azul flow is 16 km long covers 16 km2 and has a dre volume of 54 106 m3 flowgo produces at vent velocity and effusion rate values of 11 m s 1 and 600 m3 s 1 respectively the velocity value agrees well with the 12 m s 1 estimated in the field the mean effusion rate total dre volume duration was 7 47 m3 s 1 dike dimensions fissure lengths and pressure gradients along the conduit based on magma chamber depth estimates of 3 5 km produce mean effusion rates for the two eruptions that range over nearly four orders of magnitude the range being due to uncertainty in the magma viscosity dike dimensions and pressure gradient between magma chamber and vent although somewhat consistent with mean effusion rates from other techniques their wide range makes them less useful the exponentially declining effusion rates during both eruptions are consistent with release of elastic strain being the driving mechanism of the eruptions our results provide independent input parameters for previously published theoretical relationships between magma chamber pressurization and eruption rates that constrain chamber volumes and increases in volume prior to eruption as well as time constants of exponential decay during the eruption the results and theoretical relationships combine to indicate that at both volcanoes probably 25 30 of the volumetric increase in the magma chamber erupted as lava onto the surface in both eruptions the lava flow volumes are less than 1 of the magma chamber volume

slope stability models for rainfall induced lahars during long lasting eruptions abstract in this study we analyse the spatial distribution of potential lahar sources during long lasting eruptions using two slope stability models shalstab and trigrs the analysis is based on observed deposit grain size and thickness rainfall data and slope angle of the lahar events in the area surrounding the cordon caulle volcano chile that occurred during the 2011 eruption the main phase of the eruption 4 7 june was characterized by eruptive plumes from 7 to 11 km high which dispersed most of the tephra eastward toward argentina with a total estimated volume of about 1 km3 tephra fallout blanketed the region from ese to ene of the volcano with the thickness of the tephra fallout layers between approximately 1 m 15 km from vent and 0 06 m 240 km from the vent on 10 june 2011 a major lahar occurred close to the argentina chile border that reached the national road 231 28 km from the vent three field campaigns were undertaken to collect samples and data from tephra fallout deposits and triggering mechanisms in the lahar source area model input parameters were obtained from geotechnical tests and field measurements several metrics were used for model evaluation and best fit to the data were obtained for simulations considering non cohesive ash layers for shalstab and a cohesion of 0 5 kpa for trigrs both models are sensitive to the physical properties of the tephra fallout deposit and the hydraulic and materiel strength properties of the study area they both also show good agreement with field data but provide different information trigrs provides an estimate of the timing based on a storm event and location of a potentially unstable area while shalstab simulations result in landslide susceptibility classes based on critical rainfall value these outcomes provide fundamental insights into lahar triggering during long lasting volcanic eruptions and are crucial to the compilation of lahar hazard maps and emergency management plans in the south andes volcanic region

the 2007 8 volcanic eruption on jebel at tair island red sea observed by satellite radar and optical images we use high resolution optical images and interferometric synthetic aperture radar insar data to study the september 2007 january 2008 jebel at tair eruption comparison of pre and post eruption optical images reveals several fresh ground fissures a new scoria cone near the summit and that 5 9 0 1 km2 of new lava covered about half of the island decorrelation in the insar images indicates that lava flowed both to the western and to the northeastern part of the island after the start of the eruption while later lavas were mainly deposited near the summit and onto the north flank of the volcano from the insar data we also estimate that the average thickness of the lava flows is 3 8 m resulting in a bulk volume of around 2 2 107 m3 we observe no volcano wide pre or post eruption uplift which suggests that the magma source may be deep the co eruption interferograms on the other hand reveal local and rather complex deformation we use these observations to constrain a tensile dislocation model that represents the dike intrusion that fed the eruption the model results show that the orientation of the dike is perpendicular to the red sea rift implying that the local stresses within the volcanic edifice are decoupled from the regional stress field

lava effusion rate definition and measurement a review measurement of effusion rate is a primary objective for studies that model lava flow and magma system dynamics as well as for monitoring efforts during on going eruptions however its exact definition remains a source of confusion and problems occur when comparing volume flux values that are averaged over different time periods or spatial scales or measured using different approaches thus our aims are to 1 define effusion rate terminology and 2 assess the various measurement methods and their results we first distinguish between instantaneous effusion rate and time averaged discharge rate eruption rate is next defined as the total volume of lava emplaced since the beginning of the eruption divided by the time since the eruption began the ultimate extension of this is mean output rate this being the final volume of erupted lava divided by total eruption duration whether these values are total values i e the flux feeding all flow units across the entire flow field or local i e the flux feeding a single active unit within a flow field across which many units are active also needs to be specified no approach is without its problems and all can have large error up to 50 however good agreement between diverse approaches shows that reliable estimates can be made if each approach is applied carefully and takes into account the caveats we detail here there are three important factors to consider and state when measuring giving or using an effusion rate first the time period over which the value was averaged second whether the measurement applies to the entire active flow field or a single lava flow within that field and third the measurement technique and its accompanying assumptions

extension in the northern mariana inner forearc seismic reflection seamarc ii acoustic imagery and bathymetric data across the northern mariana island arc reveal dense high angle normal faulting of the arc and inner forearc faults with throws up to 780 m cut a thick sequence of stratified but otherwise undeformed sediments the major faults display a sinuous or zigzag trace that is commonly associated with rifting this pattern can be indicative of a single event of triaxial strain that requires two contemporaneous conjugate fault pairs in orthorhombic symmetry interpretation and statistical analyses of the data show activity since the last sedimentary deposition and concurrent fault development of two distinct fault populations one associated with the backarc basin s east boundary fault and the other with extension in the forearc extensional faults along the eastern side of the backarc basin near the active arc at 22 n where active rifting is occurring have approximately the same trend as the east boundary fault in the survey area we estimate the direction of extensional minimum principal strain for the east boundary fault at 22 n to be 39 1 orthogonal to the extensional strain direction measured in the inner forearc the change in strain is distinct and abrupt occurring across a 20 km wide region along the arc massif within this narrow region of transition the intermediate and minimum principal strain directions are nearly equal in magnitude with orientations midway between the dominant forearc and backarc basin principal strain directions statistical analyses show that the strain from either backarc or forearc regime is not transmitted across this narrow region the two opposing strain regimes coexist virtually juxtaposed the coexistence of these two orthogonal extensional fault regimes within the overriding plate between the mariana backarc basin and inner forearc may be explained by assuming that a weak transitional region exists along the east boundary fault zone where stresses can be more easily relieved this allows active basin wide extension in the backarc basin while the forearc data are consistent with models that require radial forearc deformation resulting from increasing arc trench curvature of the northern half of the mariana island arc

simulations of measures to control lava flows we evaluated the quantitative effects of artificial barriers water cooling and guiding channels on lava flow using the lava simulation program lavasim lava flow is basically subject to the topography around the path effusive rate and viscosity to prevent damage due to lava flow we conducted experiments in controlling the flow direction velocity and temperature the simulation demonstrated that artificial barriers can successfully change the direction of a lava flow and is more effective when placed nearly parallel to the flow direction at a point where the topography is not very steep while a barrier placed perpendicular to the flow direction can only stop the flux temporarily ultimately allowing the solidified crust to accumulate and causing the following mass to go over the barrier the water cooling trial was also effective in controlling the direction and temperature although the amount of water was as much order as 105 m3 the guiding channels successfully control the direction and inundated area but only in local areas

topographic mapping from interferometric synthetic aperture radar observations interferometric synthetic aperture radar observations provide a means for obtaining high resolution topographic terrain maps from data acquired simultaneously at two slightly displaced antennas calculation of the three dimensional coordinates of all the points in a radar image can be made from the combination of along track slant range and interferometer fringe measurements thus the result compensates for the layover exhibited by conventional radar maps and removes the false targets generated by multiple signal paths to a given object in the scene we have derived a topographic map of a portion of the san francisco bay area utilizing data that were recorded by a radar system mounted on a nasa cv990 aircraft and processed by a general purpose digital computer this map displays the height above sea level of a region approximately 11 km by 10 km in size sampled on an 11 m pixel grid uncertainties in the estimated height result from imprecise knowledge of the observational geometry radar bandwidth limitations and finite signal to noise ratios for our system each factor contributed approximately 3 4 m to a relative rms error of about 6 m additional global errors of the order of 10 m can result from inaccurate monitoring of aircraft attitude the statistical variation of the height measurements from the ocean portion of the map which we can presume to be quite flat is 2 10 m rms which is consistent with the theoretical estimate comparison of our results with u s geological survey contour maps indicates a high degree of correlation between the two sets of altitude data

on the mechanism of crustal block rotations in the central andes abstract the results of 62 paleomagnetic studies of jurassic and younger rock units from the andes of argentina bolivia chile and peru are examined for indications of the timing and mechanism of crustal block rotations rotations and their 95 confidence intervals are recalculated in pole space using a new south american reference apw path as observed by many investigators the arica deflection the change in topographic and tectonic trends located near arica northernmost chile divides the data into a counterclockwise rotated set north of the deflection and a clockwise rotated set to the south rotations range from 75 counterclockwise to nearly 50 clockwise and show substantial differences between neighboring crustal blocks south of arica there is a weak indication of a decrease in amount of rotation with distance from the plate margin but no such correlation is observed north of arica rotation increases with age both north and south of arica for rock units of 90 ma for older rock units no correlation between age and amount of rotation exists unlike the western margin of north america no significant margin parallel displacements accompany the rotations several models for rotation are examined it is concluded that a combination of late cenozoic oroclinal bending and post mid cretaceous local block rotations in response to shear driven sub crustal ductile flow can explain most of the paleomagnetic and geologic observations

landscape evolution on the eastern part of lombok indonesia related to the 1257 ce eruption of the samalas volcano abstract one of the most powerful eruptions of the holocene known as the samalas eruption had a vei of 7 and took place on lombok eastern indonesia in 1257 ce thick tephra fall covered the entire island and pyroclastic density currents pdcs up to 50 m thick buried almost half of lombok until now there has been no detailed study of the geomorphological evolution of lombok related to this eruption especially for the eastern part of lombok the main goal of this paper is therefore to reconstruct the paleo topography on the eastern part of lombok before the eruption of the samalas volcano in 1257 ce and to analyse the subsequent landscape evolution following this eruption data were collected from more than 1300 points in order to gather geomorphological geological and geochemical information for the eastern part of lombok over the last millennium data were obtained from wells which are widespread in the studied area natural outcrops mainly cliffs along the shoreline and the rivers and human made outcrops in quarries as well as from the two dimensional resistivity profiling dipole dipole array carried out on eastern lombok in 2016 and 2017 we estimate that 4435 5 5 106 m3 of pumice rich pdc over a surface of 171 km2 buried the eastern part of lombok in the aftermath of the 1257 ce eruption and as a result of this rapid pumice rich pdc deposition there has been some aggradation of the shoreline at some points resulting in most of the pre 1257 coral reef being buried after several intense rainfall events the deposited volcanic material was probably eroded during the first decades following the 1257 ce eruption by mountain streams producing lahars we estimate that the remaining volume of pumice rich pdc deposits following this erosion is 625 5 5 106 m3 over a surface of 89 km2 which is equivalent to 14 of the initial volume

development of the 1990 kalapana flow field kilauea volcano hawaii the 1990 kalapana flow field is a complex patchwork of tube fed pahoehoe flows erupted from the kupaianaha vent at a low effusion rate approximately 3 5 m3 s these flows accumulated over an 11 month period on the coastal plain of kilauea volcano where the pre eruption slope angle was less than 2 the composite field thickened by the addition of new flows to its surface as well as by inflation of these flows and flows emplaced earlier two major flow types were identified during the development of the flow field large primary flows and smaller breakouts that extruded from inflated primary flows primary flows advanced more quickly and covered new land at a much higher rate than breakouts the cumulative area covered by breakouts exceeded that of primary flows although breakouts frequently covered areas already buried by recent flows lava tubes established within primary flows were longer lived than those formed within breakouts and were often reoccupied by lava after a brief hiatus in supply tubes within breakouts were never reoccupied once the supply was interrupted during intervals of steady supply from the vent the daily areal coverage by lava in kalapana was constant whereas the forward advance of the flows was sporadic this implies that planimetric area rather than flow length provides the best indicator of effusion rate for pahoehoe flow fields that form on lowangle slopes

observations on basaltic lava streams in tubes from kilauea volcano island of hawai i from 1986 to 1997 the pu u o o kupaianaha eruption of kilauea produced a vast pahoehoe flow field fed by lava tubes that extended 10 12 km from vents on the volcano s east rift zone to the ocean within a kilometer of the vent tubes were as much as 20 m high and 10 25 m wide on steep slopes 4 10 a little farther away from the vent some tubes formed by roofing over of lava channels lava streams were typically 1 2 m deep flowing within a tube that here was typically 5 m high and 3 m wide on the coastal plain 1 tubes within inflated sheet flows were completely filled typically 1 2 m high and several tens of meters wide tubes develop as a flow s crust grows on the top bottom and sides of the tubes restricting the size of the fluid core the tubes start out with nearly elliptical cross sectional shapes many times wider than high broad flat sheet flows evolve into elongate tumuli with an axial crack as the flanks of the original flow were progressively buried by breakouts temperature measurements and the presence of stalactites in active tubes confirmed that the tube walls were above the solidus and subject to melting sometimes the tubes began downcutting progressive downcutting was frequently observed through skylights a rate of 10 cm d was measured at one skylight for nearly 2 months

radar interferogram filtering for geophysical applications the use of sar interferometry is often impeded by decorrelation from thermal noise temporal change and baseline geometry power spectra of interferograms are typically the sum of a narrow band component combined with broad band noise we describe a new adaptive filtering algorithm that dramatically lowers phase noise improving both measurement accuracy and phase unwrapping while demonstrating graceful degradation in regions of pure noise the performance of the filter is demonstrated with sar data from the ers satellites over the jakobshavns glacier of greenland

modeling the october 2005 lahars at panabaj guatemala an extreme rainfall event in october of 2005 triggered two deadly lahars on the flanks of toliman volcano guatemala that caused many fatalities in the village of panabaj we mapped the deposits of these lahars then developed computer simulations of the lahars using the geologic data and compared simulated area inundated by the flows to mapped area inundated computer simulation of the two lahars was dramatically improved after calibration with geological data specifically detailed field measurements of flow inundation area flow thickness flow direction and velocity estimates collected after lahar emplacement were used to calibrate the rheological input parameters for the models including deposit volume yield strength sediment and water concentrations and manning roughness coefficients simulations of the two lahars with volumes of 240 200 55 400 and 126 000 29 000 m3 using the flo 2d computer program produced models of lahar runout within 3 of measured runouts and produced reasonable estimates of flow thickness and velocity along the lengths of the simulated flows we compare areas inundated using the jaccard fit model sensitivity and model precision metrics all related to bayes theorem these metrics show that false negatives areas inundated by the observed lahar where not simulated and false positives areas not inundated by the observed lahar where inundation was simulated are reduced using a model calibrated by rheology the metrics offer a procedure for tuning model performance that will enhance model accuracy and make numerical models a more robust tool for natural hazard reduction

quantitative trace element analyses of silicate reference materials and a stainless steel using the proton microprobe abstract the accuracy of and detection limits for trace element analysis by the proton microprobe were determined using a stainless steel sample gk 19 and glass samples prepared from whole rock reference powders ja 1 andesite jb 1a basalt jf 1 feldspar and jgb 1 gabbro in the steel sample the measured concentrations of the major transition metals cr mn fe co ni and cu agreed with the recommended values within an average uncertainty of 3 thus the analytical software for the proton microprobe showed an excellent correction for the secondary fluorescence effect and executed appropriate theoretical calculations for accurate quantification the measured concentrations of all the other trace elements also agreed with recommended values within an average uncertainty of 7 for the glass samples the measured contents of trace ni in jb 1a and jgb 1 agreed with the recommended values within uncertainties of 14 and 3 but the contents of trace mn and cu in all the samples were an average of 30 and 40 lower than the recommended values respectively the discrepancy in mn and cu contents is considered to be due to an effect of overlap of the low energy side tail of the intense fek peak for the mnk peak and a loss of cu during preparation of the glasses respectively the concentrations of all the trace elements except cu agreed with the recommended values within 10 average for ja 1 jb 1a and jf 1 and 15 for jgb 1 in addition almost all trace elements were detected at levels of only a few micrograms per gram

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

emissions of trace elements during the 2012 2013 effusive eruption of tolbachik volcano kamchatka enrichment factors partition coefficients and aerosol contribution abstract gases and aerosols from the 2012 13 effusive eruption of tolbachik basaltic volcano kamchatka were sampled in february and may 2013 from a lava tube window located 300 m from the eruptive crater temperature at the sampling point was 1060 1070 c the chemical and isotopic compositions of the sampled gases 92 4 h 2 o 3 5 co 2 2 3 so 2 on average d from 25 0 to 38 6 correspond to a typical volcanic arc gas without dilution by meteoric or hydrothermal water halogen contents in the gases 1 37 hcl 0 5 hf were higher than average arc values the total amount of analyzed metallic and metalloid trace elements in the gas exceeded 665 ppm six most abundant trace elements k 250 ppm na 220 ppm si 74 ppm br 48 ppm cu 21 ppm and fe 12 ppm accounted for 95 of the total content of trace elements in the gas the gases contained 24 ppb re 12 ppb ag 4 9 ppb au and 0 45 ppb pt refractory rock forming elements mg al ca and some other elements such as ba and th were transported mainly in the form of silicate microspheres and altered rock particles the concentrations of metals in the eruptive tolbachik gases are higher than the corresponding concentrations in high temperature fumaroles worldwide although the mutual ratios of the elements are approximately the same the gas magma partition coefficients of eleven elements exceed unity including the non metals f s cl br as se and te and the rare metals cd re tl and bi despite the relatively low concentrations of trace elements in the volcanic gases at the highest temperatures superficial magma degassing provides information on the sources and sinks of metals

magma compressibility and the missing source for some dike intrusions 1 dike intrusions are often accompanied by localized deflation interpreted as depressurizing magma chambers feeding the dike in some cases the inferred volume decrease is a factor of 4 or 5 less than the volume increase of the dike here we explore whether this discrepancy can be explained by compressibility of the magma combined with the fact that cracks are much more compliant than equidimensional magma chambers if pressure changes are small the magma compressibility m is constant and the dike ends up in hydrostatic equilibrium with an ellipsoidal magma chamber at the same depth the ratio rv of the volume of the crack to the volume lost by the chamber is rv 1 4 m 3 1 where is the host rock rigidity for gas poor magmas m 0 6 2 10 10 pa 1 and 3 25 gpa we find 1 2 rv 7 7 large changes in magma compressibility due to gas exsolution increase rv

stress fields of the overriding plate at convergent margins and beneath active volcanic arcs tectonic stress fields in the overriding plate at convergent plate margins are complex and vary on local to regional scales volcanic arcs are a common element of overriding plates stress fields in the volcanic arc region are related to deformation generated by subduction and to magma generation and ascent processes analysis of moment tensors of shallow and intermediate depth earthquakes in volcanic arcs indicates that the seismic strain field in the arc region of many convergent margins is subhorizontal extension oriented nearly perpendicular to the arc a process capable of generating such a globally consistent strain field is induced asthenospheric corner flow below the arc region

topographic analyses of k lauea volcano hawai i from interferometric airborne radar we analyze digital topographic data collected in september 1993 over a 500 km2 portion of k lauea volcano hawai i by the c band 5 6 cm wavelength topographic synthetic aperture radar topsar airborne interferometric radar field surveys covering an 1 km2 area of the summit caldera and the distal end of an 8 m thick a flow indicate that the 10 m spatial resolution topsar data have a vertical accuracy of 1 2 m over a variety of volcanic surfaces after conversion to a common datum topsar data agree favorably with a digital elevation model dem produced by the u s geological survey usgs with the important exception of the region of the ongoing eruption which postdates the usgs dem this dem comparison gives us confidence that subtracting the usgs data from topsar data will produce a reasonable estimate of the erupted volume as of september 1993 this subtraction produces dense rock equivalent dre volumes of 392 439 and 90 106 m3 for the pu u k pa ianah and episode 50 53 stages of the eruption respectively these are 124 89 and 94 of the volumes calculated by staff of the hawaiian volcano observatory hvo but do not include lava of k pa ianah and episodes 50 53 that flowed into the ocean and are thus invisible to topsar accounting for this lava increases the topsar volumes to 124 159 and 129 of the hvo volumes including the 2 m uncertainty derived from the field surveys produces topsar derived volumes for the eruption as a whole that range between 81 and 125 of the usgs derived values the vesicularity and ocean corrected topsar volumes yield volumetric eruption rates of 4 5 4 5 and 2 7 m3 s for the three stages of the eruption which compare with hvo derived values of 3 6 2 8 and 2 1 m3 s respectively our analysis shows that care must be taken when vertically registering the topsar and usgs dems to a common datum because c band topsar penetrates only partially into thick forest and therefore produces a dem within the tree canopy whereas the usgs dem is adjusted for vegetation

on the evolution of marginal basins this paper presents a critique of the models proposed to explain the formation of back arc basins those marginal basins located behind active or inactive trench systems and whose origin is inferred to be subduction related several marginal basins are not generically back arc but rather are related to the initiation of subduction continental rifting or plate boundary readjustments models of back arc basin formation must be consistent with the observed characteristics of back arc basins which include the following 1 the crustal structure and magnetic lineation patterns of back arc and major ocean basins are similar implying that the processes of crustal accretion in both tectonic environments are also similar 2 arc trench systems without currently spreading back arc basins are in the majority 3 although petrographically similar to and within the compositional range of midocean ridge basalts morb back arc basin basalts babb show consistent geochemical differences from n type morb and are in many respects transitional toward island arc tholeiites iat 4 distinct babb and iat mantle sources may be in extremely close proximity during the initial rifting of an island arc 5 rifting of the volcanic arc may occur on either side of the line of active volcanoes 50 km and may vary from one side to the other along strike 6 a remnant arc is not always developed 7 back arc spreading may be initiated above subducted lithosphere but with time the center of spreading migrates away from the arc and no longer overlies a benioff zone and 8 there is no simple correlation between the timing of global plate reorganizations and the formation of back arc basins the three main classes of models proposed to explain back arc basin formation are mantle diapirism induced aesthenospheric convection and global plate kinematics because the driving mechanism for the first two classes of models is provided by the local subducting slab these models fail to explain the temporal and spatial distribution of back arc basins the third class of models proposes that back arc basins should form whenever global plate interactions require divergence between the overriding plate and the trench line the latter being observed to move only in a seaward direction with respect to a hot spot reference frame however the tectonic setting of several back arc basins suggests that they represent more than just a passive response to kinematic boundary conditions for example some dynamic back arc force seems required to move the new hebrides tonga and south sandwich forearc plates seaward neither the local nor the global models adequately specify the necessary and or sufficient conditions for back arc basin formation a better understanding of the processes associated with the initiation of rifting and of subduction is required arc rifting and back arc spreading are both secondary phenomena to subduction but they are also primary extensional phenomena fundamentally analogous to continental rifting and the initial spreading of major ocean basins

hydrothermal manganese deposits in island arcs and related to subduction processes a possible model for genesis abstract submarine hydrothermal manganese deposits have recently been recovered from a number of island arc settings such as the tonga kermadec ridge ogasawa arc trench bonin arc mariana arc and bismarck archipelago it has also been shown that large volumes of fluid are released during dewatering of subducting sediments the amount of water subducted globally 1 5 10 15 g yr 1 is about two orders of magnitude less than that advected at ocean spreading centres hydrothermal manganese deposits may be precipitated from these waters the distribution of terrestrial manganese deposits around the circum pacific belt and elsewhere suggests that some of these deposits may be of hydrothermal origin related to subduction processes at present most of these terrestrial deposits are inadequately described to permit the establishment of the precise tectonic setting at the time of deposition of the manganese deposits and therefore an unambiguous assessment of this hypothesis

image based measurement of flux variation in distal regions of active lava flows 1 understanding the processes involved with the advance of lava flows is critical for improving hazard assessments at many volcanoes here we describe the application of computer vision and oblique photogrammetric techniques to visible and thermal images of active a flows in order to investigate distal flow processes at mount etna sicily photogrammetric surveys were carried out to produce repeated topographic data sets for calculation of volumetric lava flux at the flow fronts velocity profiles from a distal channel were obtained by rectification of a thermal image sequence and are used to investigate the rheological properties of the lava significant variations of the magma flux were observed and pulses of increased flux arrived within the flow front region on timescales of several hours the pulses are believed to be the distal result of more frequent flux changes observed in the vent region hence they reflect the importance of flow processes which are believed to cause the coalescence of flux pulses along the channel system as well as short period variations in effusion rate in considering advance processes for the individual flow fronts it must be assumed that they were fed by a highly unsteady flux which was volumetrically at least an order of magnitude lower than that observed near the vent

intra arc sedimentation in a low lying marginal arc eocene clarno formation central oregon abstract the clarno formation consists of andesitic volcaniclastic rocks interstratified with clayey paludal sediments and lava flows and cut locally by irregular hypabyssal stocks dikes and sills lateral lithofacies variations are pronounced and intrusive and extrusive volcanic rocks appear haphazardly emplaced throughout the formation a range of sedimentary environments is represented in the volcaniclastic sequences including near vent flow and breccia accumulations bouldery high gradient braided streams and relatively low gradient sandy tuff braidplains associated with paludal deposits we infer that the coarse grained volcaniclastic rocks of the clarno formation accumulated largely in volcanic flank and apron settings the formation s stratigraphy reflects an abundance of small volcanoes erupting into sedimentary lowlands and only a few scattered remnants of large central vent volcanoes are known the absence of systematic variation across a large outcrop belt argues against derivation of the succession from a line of volcanoes beyond the reaches of the present outcrop we infer that the arc was composed of small to medium size volcanoes arranged haphazardly across a broad area the sedimentary succession most probably accumulated in a series of intra arc depressions formed by diffuse block rotation driven by oblique subduction during the eocene