monitoring volcano slope instability with synthetic aperture radar a review and new data from pacaya guatemala and stromboli italy volcanoes abstract volcano slope instability manifests in many forms ranging from steady state to punctuated movement or shallow erosion to deep seated spreading the interplay of gravity magmatic or hydrothermal fluids and often active tectonics on a volcanic edifice results in complex spatial and temporal variations in deformation kinematics while this makes the recognition assessment and monitoring of volcano slope instability challenging advancements in synthetic aperture radar sar technology have significantly accelerated our knowledge of instability phenomena and our ability to assess their hazards this review discusses the applications and challenges of sar imagery to various slope instabilities at volcanoes around the world sar amplitude images are powerful tools for mapping areas of geomorphological changes these snapshots can be combined using interferometric sar insar to create multi temporal deformation maps that provide unique information on the evolution of slope failures space borne insar has become an economic way to detect changes at volcanoes at very high resolution ground based insar gbinsar can produce frequent sar images on the order of seconds to minutes propelling insar from monitoring to surveillance and early warning applications however interpreting insar derived ground deformation signals related to volcano slope instability is still challenging deformation from magma rise or variations in hydrothermal systems may be inseparable from persistent deep seated flank motion similarly shallow or localized ground motion may occur in overlap with thermal contraction or subsidence of newly emplaced lava or tephra deposits if triggered to failure landslides can vary over several magnitudes from small volume rock falls that pose only a localized hazard to large volume debris avalanches capable of travelling tens of kilometers away from the source these collapses can have cascading hazards if accompanied by directed blasts or tsunamis emphasizing the need for continued advancement of our understanding of these events and their consequences to highlight the utility of sar for measuring and monitoring mass movement two case studies of stromboli italy and pacaya guatemala volcanoes are discussed in detail where recent instability events persistent volcanic activity and ground truth constraints have resulted in excellent case histories in applying sar imagery to understand these potentially hazardous slope instabilities

chemostratigraphic correlation of the middle and upper proterozoic between the yanshan and shennongjia basins abstract this paper is concerned with the distribution of major elements transition elements ree sulphur isotopes and organic carbon in parts of the proterozoic of china lithophile and ore forming elements such as ba sr mn k v fe and cu are good chemical markers and are used for correlation in the yanshan basin the ba content is particularly high in red beds of several sections ba fe k mn v ca mg and organic carbon have provided useful data for correlation between the yanshan and shennongjia basins these elements reflect in varying degree widespread important geological events we propose several chemostratigraphic stages with divisions approximately at 1950 ma 1600 ma 1200 ma and 1000 ma these divisions correspond approximately to major time boundaries in the proterozoic of china 34s values for sulphur isotopes in proterozoic strata in china are between 13 and 27 ree patterns normalized to chondritic abundances for middle upper proterozoic rocks of china are similar to post archaean sedimentary rocks of australia they also show eu depletion and eu eu ratios between 0 6 and 0 8

the largest local tsunami in 20th century hawaii abstract the 1975 kalapana tsunami had waves that ran up to 14 6 m high it deposited a discontinuous basalt boulder and carbonate sand veneer on the halape apua point coast of the island of hawaii these deposits fine up to 320 m inland and 10 4 masl the deposit is an unusual sedimentary association in which sand sized material is a minor component both boulders and sand fine inland over the same distance not in a continuum and they terminate at approximately the same point boulders were transported and deposited by only one wave the second whereas the sand deposit was in part reworked by subsequent smaller waves boulder transport by one wave in association with a relative paucity of finer sediment reduced internal clast interference and helped to impart a marked fining inland sequence this fining inland sequence is at odds with work that proposes a chaotic or random nature of boulder emplacement by tsunamis where there is abundant sediment clast interference during entrainment will impede the fining inland of larger material in cases where it is difficult to determine the nature of inundation or the number of tsunamis a chaotic pattern may simply reflect deposition by storm tsunami or a combination of the two

active thrusting and folding along the northern tien shan and late cenozoic rotation of the tarim relative to dzungaria and kazakhstan we have studied geometries and rates of late cenozoic thrust faulting and folding along the northern piedmont of the tien shan mountain belt west of urumqi where the m 8 3 manas earthquake occurred on december 23 1906 the northern range of the tien shan rising above 5000 m overthrusts a flexural foredeep filled with up to 11 000 m of sediment of the dzungarian basement our fieldwork reveals that the active thrust reaches the surface 30 km north of the range front within a 200 km long zone of neogene quaternary anticlines fault scarps are clearest across inset terraces within narrow valleys incised through the anticlines by large rivers flowing down from the range in all the valleys the scarps offset vertically the highest terrace surface by the same amount 10 2 0 7 m inferring an early holocene age 10 2 kyr for this terrace which is continuous with the largest recent fans of the piedmont yields a rate of vertical throw of 1 0 0 3mm yr on the main active thrust at the surface a quantitative morphological analysis of the degradation of terrace edges that are offset by the thrust corroborates such a rate and yields a mass diffusivity of 5 5 2 5 m 2 kyr a rather fresh surface scarp 0 8 0 15 m high that is unlikely to result from shallow earthquakes with 6 m 7 in the last 230 years is visible at the extremities of the main fold zone we associate this scarp with the 1906 manas earthquake and infer that a structure comprising a deep basement ramp under the range gently dipping flats in the foreland and shallow ramps responsible for the formation of the active fault propagation anticlines could have been activated by that earthquake if so the return period of a 1906 type event would be 850 380 years the small size of the scarp for an earthquake of this magnitude suggests that a large fraction of the slip at depth 2 3 is taken up by incremental folding near the surface comparable earthquakes might activate flat detachments and ramp anticlines at a distance from the front of other rising quaternary ranges such as the san gabriel mountains in california or the mont blanc aar massifs in the alps we estimate the finite cenozoic shortening of the folded dzungarian sediments to be of the order of 30 km and the cenozoic shortening rate to have been 3 1 5 mm yr assuming comparable shortening along the tarim piedmont and minor additional active thrusting within the mountain belt we infer the rate of shortening across the tien shan to be at least 6 3 mm yr at the longitude of manas 85 5 e a total shortening of 125 30 km is estimated from crustal thickening assuming local airy isostatic equilibrium under the same assumption serial n s sections imply that cenozoic shortening across the belt increases westwards to 203 50 km at the longitude of kashgar 76 e as reflected by the westward increase of the width of the belt this strain gradient implies a clockwise rotation of tarim relative to dzungaria and kazakhstan of 7 2 5 around a pole located near the eastern extremity of the tien shan west of hami 96 e 43 5 n comparable to that revealed by paleomagnetism between tarim and dzungaria 8 6 8 7 a 6 mm yr rate of shortening at the longitude of manas would imply a rate of rotation of 0 45 m y and would be consistent with a shortening rate of 12 mm yr north of kashgar taking such values to be representative of late cenozoic rates would place the onset of reactivation of the tien shan by the india asia collision in the early to middle miocene 16 22 9 m y in accord with the existence of particularly thick late neogene and quaternary deposits such reactivation would thus have started much later than the collision roughly at the time of the great mid miocene changes in tectonic regimes denudation and sedimentation rates observed in southeast asia the himalayas and the bay of bengal and of the correlative rapid change in seawater sr isotopic ratio 20 to 15 ma like these other changes the rise of the tien shan might be a distant consequence of the end of indochina s escape

inversion of complex body waves we have developed a method that inverts seismic body waves to determine the mechanism and rupture pattern of earthquakes the rupture pattern is represented as a sequence of subevents distributed on the fault plane this method is an extension of our earlier method in which the subevent mechanisms were fixed in the new method the subevent mechanisms are determined from the data and are allowed to vary during the sequence when subevent mechanisms are allowed to vary however the inversion often becomes unstable because of the complex trade offs between the mechanism the timing and the location of the subevents many different subevent sequences can explain the same data equally well and it is important to determine the range of allowable solutions some constraints must be imposed on the solution to stabilize the inversion we have developed a procedure to explore the range of allowable solutions and appropriate constraints in this procedure a network of grid points is constructed on the i plane where and i are respectively the onset time and the distance from the epicenter of a subevent the best fit subevent is determined at all grid points then the correlation is computed between the synthetic waveform for each subevent and the observed waveform the correlation as a function of and i and the best fit mechanisms computed at each i grid point depict the character of allowable solutions and facilitate a decision on the appropriate constraints to be imposed on the solution the method is illustrated using the data for the 1976 guatemala earthquake

constraining the dip of shallow shallowly dipping thrust events using long period love wave radiation patterns applications to the 25 october 2010 mentawai indonesia and 4 may 2018 hawaii island earthquakes constraining precise faulting geometry for shallow shallowly dipping thrust earthquakes is a common challenge plate boundary megathrust faults near the trench and decollement faults beneath volcanic islands and nappes may dip only a few degrees long period point source moment tensor waveform inversions provide limited resolution of shallow fault dip angle the possibility of splay faulting requires precise dip estimation high sensitivity to dip is provided by long period love wave amplitude azimuthal radiation patterns which undergo rapid change from four lobed to two lobed as dip decreases from 10 to 0 modeling variability in love wave nodal amplitudes allows the dip to be determined to within a few degrees this is demonstrated for the 25 october 2010 mentawai m ww 7 8 tsunami earthquake which ruptured the 2 0 5 0 dipping megathrust beneath the shallow sedimentary wedge offshore of indonesia and the 4 may 2018 hawaii island m ww 6 9 thrust earthquake which ruptured the 2 5 7 5 dipping decollement under the island flank

a geologically constrained monte carlo approach to modeling exposure ages from profiles of cosmogenic nuclides an example from lees ferry arizona we present a user friendly and versatile monte carlo simulator for modeling profiles of in situ terrestrial cosmogenic nuclides tcns our program available online at http geochronology earthsciences dal ca downloads models html permits the incorporation of site specific geologic knowledge to calculate most probable values for exposure age erosion rate and inherited nuclide concentration while providing a rigorous treatment of their uncertainties the simulator is demonstrated with 10be data from a fluvial terrace at lees ferry arizona interpreted constraints on erosion based on local soil properties and terrace morphology yield a most probable exposure age and inheritance of 83 9 14 1 19 1 ka and 9 49 2 52 1 21 104 atoms g 1 respectively 2 without the ability to apply some constraint to either erosion rate or age shallow depth profiles of any cosmogenic nuclide except for nuclides produced via thermal and epithermal neutron capture e g 36cl cannot be optimized to resolve either parameter contrasting simulations of 10be data from both sand and pebble sized clasts within the same deposit indicate grain size can significantly affect the ability to model ages with tcn depth profiles and when possible sand not pebbles should be used for depth profile exposure dating

an integrated crustal velocity field for the central andes 1 we present an integrated velocity field for the central andes derived from gps observations collected between january 1993 and march 2001 that eliminates the velocity bias between the south america nazca plate project snapp and central andes gps project cap velocity fields published by norabuena et al 1998 and bevis et al 1999 the reference frame is realized by minimizing the motion of eight continuous gps stations and one rover gps station located in the stable core of the south american plate the rms horizontal motion of these stations is just 1 1 mm yr the amplitude of these residual motions is roughly compatible with expected levels of measurement error in our new solution five of the six snapp stations located just outside the orogenic belt are effectively stationary and the velocities for adjacent cap and snapp stations now agree at a level consistent with their formal uncertainties

reconciling seismic and geodetic models of the 1989 kilauea south flank earthquake 1 the inconsistency in depth between geodetic and seismic models of the 1989 kilauea earthquake has long been puzzling previous attempts to incorporate elastic heterogeneity have deepened geodetic models substantially bringing them closer to the seismic depth however recent studies that have included heterogeneity for other faults indicate that the effect is not so great we show here that elastic heterogeneity has a relatively minor effect on the geodetic model depth for this earthquake also however by combining three different sets of geodetic data we are able to get a more accurate estimate of the depth which does in fact coincide with the hypocentral depth at the 95 confidence level when we consider that static elastic parameter values are commonly less than dynamic values the agreement is even better furthermore the depth is consistent with the earthquake having occurred at the interface between the volcanic pile and the pre volcanic seafloor

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

editorial flank dynamics sector collapses lahars and rockfalls analysis monitoring and modelling of small to large scale volcanic slope instability slope dynamics in volcanic environments comprise a wide spectrum of phenomena from large lateral collapse to shallow debris remobilization which may represent a major threat for human communities and infrastructures many volcanos built up from the ocean floor and large portions of the volcano edifice are submerged in these settings only the edifice s summit can be investigated by terrestrial remote sensing and in situ approaches growth and destruction including tectonics and gravitational phenomena affect entire volcano flanks and are not limited to the physical boundary of the sea level but could comprise their subaqueous parts

vertical deformation measurements on the submerged south flank of kīlauea volcano hawai i reveal seafloor motion associated with volcanic collapse 1 a four year seafloor geodetic study was conducted to measure vertical deformation of the submerged south flank of k lauea volcano on the active hilina slump the hilina slump is a site of significant deformation major earthquakes with ground cracking and associated tsunami a new technique was developed to measure vertical deformation on the seafloor using pressure sensors in campaign style surveys the data revealed the midslope bench of the offshore slump is uplifting at 9 0 2 4 cm a but transitions to no significant deformation on the outer bench and slump toe elastic half space inverse modeling of these data show that the observed deformation can be fit by seaward slip of 28 1 7 3 cm a extending from the east rift zone 27 0 0 5 km on a gently island dipping decollement fault plane at a depth of approximately 7 km modeling suggests that the outer bench is not currently deforming because the majority of k lauea is submerged these offshore measurements are necessary to constrain the seaward extent of fault slip and the motion of the outer bench

rapid deformation of the south flank of kilauea volcano hawaii the south flank of kilauea volcano has experienced two large magnitude m 7 2 and m 6 1 earthquakes in the past two decades global positioning system measurements conducted between 1990 and 1993 reveal seaward displacements of kilauea9s central south flank at rates of up to about 10 centimeters per year in contrast the northern side of the volcano and the distal ends of the south flank did not displace significantly the observations can be explained by slip on a low angle fault beneath the south flank combined with dilation deep within kilauea9s rift system both at rates of at least 15 centimeters per year

rapid deformation of kilauea volcano global positioning system measurements between 1990 and 1996 campaign global positioning system gps measurements from 1990 to 1996 are used to calculate surface displacement rates on kilauea volcano hawaii the gps data show that the south flank of the volcano which has generated several large earthquakes in the past 3 decades is displacing at up to 8 cm yr to the south southeast the summit and rift zones are subsiding with maximum subsidence rates of 8 cm yr observed a few kilometers south of the summit caldera elastic dislocation modeling of the gps data suggests that the active sources of deformation include deep rift opening along the upper east and east rift zone fault slip along a subhorizontal fault near the base of the volcano and deflation near the summit caldera a nonlinear optimization algorithm was used to explore the parameter space and to find the best fitting source geometry there is a broad range of model geometries that fit the data reasonably well however certain models can be ruled out including those that have shallow rift opening or shallow fault slip some offshore aseismic slip on a fault plane that dips between 258 north northwest and 88 south southeast is required best fitting slip and rift opening rates are 23 28 cm yr although rates as low as 10 cm yr are permitted by the data

long period source characteristics of the 1975 kalapana hawaii earthquake the 1975 kalapana hawaii earthquake occurred under the highly mobile south flank of kilauea volcano it has been interpreted variously as a normal faulting earthquake a thrust faulting earthquake and a landslide primary evidence for the landslide model has been the failure of previous faulting models to explain the observed love wave radiation pattern and the tsunami amplitudes generated by the kalapana event here we present a reanalysis of the long period digital seismic data for this event centroid moment tensor analysis shows that the seismic radiation pattern can be explained well by thrust faulting on a plane dipping shallowly landward the seismic moment of 3 8 10 20 n m m w 7 7 that we determine is approximately twice as large as earlier estimates the geometry and seismic moment of the focal mechanism determined here are consistent with the observed tsunami amplitudes inversion of long period body wave waveforms shows that the earthquake source duration 72 sec is unusually long for an earthquake of this size

slope failure and volcanic spreading along the submarine south flank of kilauea volcano hawaii 1 new multichannel reflection data and high resolution bathymetry over the submarine slopes of kilauea volcano provide evidence for current and prior landsliding suggesting a dynamic interplay among slope failure regrowth and volcanic spreading disrupted strata along the upper reaches of kilauea s flank denote a coherent slump correlated with the active hilina slump the slump comprises mostly slope sediments underlain by a detachment 3 5 km deep extension and subsidence along the upper flank is compensated by uplift and folding of the slump toe which surfaces about midway down the submarine flank uplift of strata forming papa u seamount and offset of surface features along the western boundary of kilauea indicate that the slump has been displaced 3 km in a south southeast direction this trajectory matches coseismic and continuous ground displacements for the hilina slump block on land and contrasts with the southeast vergence of the rest of the creeping south flank to the northeast slope sediments are thinned and disrupted within a recessed region of the central flank demonstrating catastrophic slope failure in the recent past debris from the collapsed flank was shed into the moat in front of kilauea building an extensive apron seaward sliding of kilauea s flank offscraped these deposits to build an extensive frontal bench a broad basin formed behind the bench and above the embayed flank uplift and back tilting of young basin fill indicate recent and possibly ongoing bench growth the hilina slump now impinges upon the frontal bench this buttress may tend to reduce the likelihood of future catastrophic detachment

mechanical models of trishear like folds previous workers have formulated velocity descriptions of the trishear kinematic model of fault propagation folds which are inherently non unique we present two complete mechanical models of fault related folding and assess the validity of the assumptions used in the assignment of velocities in the trishear description and to eliminate the untenable situation of an infinite number of possible solutions for velocity fields the mechanical model of forced folding forced fold based on viscous folding theory is used to derive the velocity fields in an anisotropic sedimentary cover overlying faulted and displaced rigid basement blocks the solution of displacements around a stress free fault in an elastic body is used to model fault arrest folds that form around a fault imbedded in a deformed medium the results indicate that the velocity fields assumed in the trishear model more closely resemble the velocities derived in the mechanical forced fold model than the mechanical model of fault arrest folding the forced fold model produces a triangular region of concentrated deformation similar to the trishear region assumed in the kinematic models while the deformation produced by the fault arrest model is not concentrated within a triangular zone

the use of earthquake rate changes as a stress meter at kilauea volcano stress changes in the earth s crust are generally estimated from model calculations that use near surface deformation as an observational constraint but the widespread correlation of changes of earthquake activity with stress has led to suggestions that stress changes might be calculated from earthquake occurrence rates obtained from seismicity catalogues although this possibility has considerable appeal because seismicity data are routinely collected and have good spatial and temporal resolution the method has not yet proven successful owing to the non linearity of earthquake rate changes with respect to both stress and time here however we present two methods for inverting earthquake rate data to infer stress changes using a formulation for the stress and time dependence of earthquake rates application of these methods at kilauea volcano in hawaii yields good agreement with independent estimates indicating that earthquake rates can provide a practical remote sensing stress meter

retrieving real time co seismic displacements using gps glonass a preliminary report from the september 2015 mw 8 3 illapel earthquake in chile

holocene shortening rates of an andean front thrust southern precordillera argentina abstract a significant part of the quaternary shortening between the chilean trench and the relative stable interior of the south american plate at the pampean flat slab 27 33 s has been accommodated at the eastern foothills of the andes and mainly within a narrow neotectonic belt along the eastern side of the argentine precordillera tectonic geodesy results point out that this area is being shortened at a 2 4 mm a rate whereas shortening rates estimated over longer time periods 1 20 ma suggest values ranging from 1 mm a to 16 mm a geomorphic and geologic evidence indicate that the east directed las higueras thrust system is one of the main structures that has accommodated quaternary deformation at this section of the andean orogenic front 32 05 32 35 s an outcrop exhibiting the thrust propagation into fluvial sediments allows the holocene shortening rates at the northern end of this structure to be estimated based on retrodeformation of holocene strata and radiocarbon dating of two charcoal bearing beds estimated shortening rates yielded mean values of 1 90 0 28 mm a for the last 4495 143 cal yr bp and 1 53 0 26 for the last 8245 48 cal yr bp these results pose some uncertainties due to the incompleteness of the exposed deformation at the hanging wall however they correspond to a key timescale which helps to bridge the gap between rates derived from the short term gps data and the long term permanent deformation rates obtained through geologic studies although the estimated rates suggest that slip on the thrust could have accelerated during the last 4 ka more data are necessary to reliably address this key issue

seismic and aseismic moment budget and implication for the seismic potential of the parkfield segment of the san andreas fault this study explores methods to assess the seismic potential of a fault based on geodetic measurements geological information of fault slip rate and seismicity data the methods are applied to the parkfield section along the san andreas fault saf at the transition zone between the saf creeping segment in the north and the locked section of cholame to the south where m w 6 earthquakes occurred every 24 5 yrs on average since the 1857 m w 7 7 fort tejon earthquake we compare the moment released by the known earthquakes and associated postseismic deformation with the moment deficit accumulated during the interseismic period derived from geodetic measurement of interseismic strain we find that the recurring m w 6 earthquakes are insufficient to balance the slip budget we discuss and evaluate various possible scenarios which might account for the residual moment deficit and implications of the possible magnitude and return period of m w 6 earthquakes on that fault segment the most likely explanation is that this fault segment hosts m w 6 5 7 5 earthquakes with a return period of 140 300 yrs such events could happen as independent earthquakes in conjunction with ruptures of the carrizo plain segment of the saf we show how the results from our analysis can be formally incorporated in probabilistic seismic hazard assessment assuming various magnitude frequency distribution and renewal time models

kinematic earthquake source inversion and tsunami runup prediction with regional geophysical data rapid near source earthquake source modeling relying only on strong motion data is limited by instrumental offsets and magnitude saturation adversely affecting subsequent tsunami prediction seismogeodetic displacement and velocity waveforms estimated from an optimal combination of high rate gps and strong motion data overcome these limitations supplementing land based data with offshore wave measurements by seafloor pressure sensors and gps equipped buoys can further improve the image of the earthquake source and prediction of tsunami extent inundation and runup we present a kinematic source model obtained from a retrospective real time analysis of a heterogeneous data set for the 2011 mw9 0 tohoku oki japan earthquake our model is consistent with conceptual models of subduction zones exhibiting depth dependent behavior that is quantified through frequency domain analysis of slip rate functions the stress drop distribution is found to be significantly more correlated with aftershock locations and mechanism types when off shore data are included the kinematic model parameters are then used as initial conditions in a fully nonlinear tsunami propagation analysis notably we include the horizontal advection of steeply sloping bathymetric features comparison with post event on land survey measurements demonstrates that the tsunami s inundation and runup are predicted with considerable accuracy only limited in scale by the resolution of available topography and bathymetry we conclude that it is possible to produce credible and rapid kinematic source models and tsunami predictions within minutes of earthquake onset time for near source coastal regions most susceptible to loss of life and damage to critical infrastructure regardless of earthquake magnitude

three dimensional seismic velocity structure of mauna loa and kilauea volcanoes in hawaii from local seismic tomography we present a new three dimensional seismic velocity model of the crustal and upper mantle structure for mauna loa and kilauea volcanoes in hawaii our model is derived from the first arrival times of the compressional and shear waves from about 53 000 events on and near the island of hawaii between 1992 and 2009 recorded by the hawaiian volcano observatory stations the vp model generally agrees with previous studies showing high velocity anomalies near the calderas and rift zones and low velocity anomalies in the fault systems the most significant difference from previous models is in v p vs structure the high vp and high v p vs anomalies below mauna loa caldera are interpreted as mafic magmatic cumulates the observed low vp and high v p vs bodies in the kaoiki seismic zone between 5 and 15 km depth are attributed to the underlying volcaniclastic sediments the high vp and moderate to low vp vs anomalies beneath kilauea caldera can be explained by a combination of different mafic compositions likely to be olivine rich gabbro and dunite the systematically low vp and low vp vs bodies in the southeast flank of kilauea may be caused by the presence of volatiles another difference between this study and previous ones is the improved vp model resolution in deeper layers owing to the inclusion of events with large epicentral distances the new velocity model is used to relocate the seismicity of mauna loa and kilauea for improved absolute locations and ultimately to develop a high precision earthquake catalog using waveform cross correlation data 2014 american geophysical union all rights reserved

coastal lithofacies and biofacies associated with syndepositional dolomitization and silicification draken formation upper riphean svalbard the draken formation 120 250 m of northeast spitsbergen svalbard forms part of a thick upper proterozoic carbonate platform succession it consists predominantly of intraformational dolomitic conglomerates with excellent textural preservation six main lithofacies were recognized in the field quartz sandstones stromatolitic mats conglomerates with silicified intraclasts dolostone conglomerates with desiccated mudrocks oolitic pisolitic grainstones and fenestral dolostones a series of five main gradational biofacies were recognized from silicified and rare calcified microfossils biofacies 1 represents low energy subtidal benthos erect filaments and plankton acritarchs and vase shaped microfossils whereas biofacies 2 to 5 are microbial mat assemblages with filamentous mat builders and associated dwellers and washed in plankton ranging from basal intertidal to high intertidal supratidal colour values a measure of the lightness of the colour shade of sawn rock samples were quantified using a munsell chart and exhibit a pronounced variation means of major groups varying from 4 0 to 5 95 across the spectrum of subtidal to supratidal sediments as inferred from other criteria the lightening in progressively more exposed sediments is related to lowering of organic carbon contents probably mainly by oxidation six types of early cement have been recognized calcite microspar type 1 is common as a subtidal cement in many proterozoic formations whereas types 2 subtidal isopachous fringes 3 subtidal hardground dolomicrite and 4 intertidal meniscus dolomicrite are very similar to phanerozoic examples except for their dolomitic mineralogy types 5 and 6 are complex and variable dolomite growths associated with expansion and replacive phenomena they characterize the fenestral lithofacies and compare with modern supratidal cements consideration of diagenetic fabrics and truncation textures of intraclasts indicates that leaching dolomitization silicification were all significant syndepositional processes altering the original metastable carbonates the data set provides evidence for a spectrum of peritidal environments including ooid shoals protected subtidal tidal sandflats and protected carbonate mudflats different sections show a preponderance of particular facies the coastal lithofacies continuum was completely dolomitized unlike offshore to ooid shoal facies of adjacent formations dolomitization thus bears a relationship to depositional bathymetry although hydrodynamics clearly have a role the potential importance of whiting precipitation in raising mg ca in marginal marine environments is also stressed

mechanical models of the 1975 kalapana hawaii earthquake and tsunami abstract our objective is to produce a mechanically realistic model for the 1975 kalapana event that explains the overall deformation and geological evolution of the south flank of kilauea and reproduces most known earthquake and tsunami observations to do this we present a new structural interpretation of geological data from kilauea along with modeling of the tsunami using recent seismic analyses in so doing we hypothesize an offshore north facing normal fault that we call the kalapana fault because of its limited onshore expression near kalapana we argue that several different interpretations of seismic data are simultaneously true each coinciding with a specific geological structure we perform a direct numerical simulation of the 1975 kalapana tsunami and we report results for the near field and the far field simultaneously our interpretation of geological structures provides sufficient constraints to derive three tsunami sources the kalapana fault a slump and a thrust fault these three tsunami sources are mechanically coupled according to our model of kilauea evolution although their contributions to tsunami generation remain relatively distinct slump displacement predicts a small shear stress beneath kilauea volcano we find that our geological interpretation and tsunami sources are sufficiently robust to reproduce almost all tsunami observations repeated earthquakes like the 1975 kalapana event will produce extension and subsidence in the upper flank of kilauea as well as compression and uplift of the toe region accompanied by limited slump movement these deformations appear to stabilize kilauea and make catastrophic failure of the volcano flank less likely at this time

the compositional evolution of meso and neoproterozoic carbonate rocks abstract petrochemical data on meso and neoproterozoic riphean and vendian carbonates from the urals and siberia suggest that upper riphean carbonate rocks differ from both older and younger equivalents in their lower dolomite limestone and higher strontium abundance it is suggested that the directional increase in the ca content of the earth s sedimentary shell results from the hydrothermal interaction of seawater with oceanic basalts manifest in a progressive reduction of dolomite formation throughout geological history superimposed cyclical variation in dolomite limestone may reflect secular changes in the rate of co 2 introduction into the atmosphere in particular the relative low dolomite limestone ratio seen in upper riphean carbonates may reflect low p co 2

flexure and seismicity beneath the south flank of kilauea volcano and tectonic implications we present a refined model for the tectonic behavior of kilauea volcano s south flank including flexure calculations and microearthquake data supportive of the model in our model the south flank moves seaward over the downwardly flexed ocean crust in a manner analogous to an accretionary prism in a subduction zone the flank is driven seaward by the gravitational stresses inherent in its shape this driving force is augmented by high density material in kilauea s rift zones elastic flexure calculations predict a configuration for the downwardly flexed crust which agrees with previous seismic refraction and gravity modeling this configuration in turn is consistent with the hypothesis that many earthquakes occur along the volcanic pile ocean crust interface focal mechanisms for these events predominantly indicate southeastward directed overthrusting we also propose a relationship between rift zone intrusion and south flank displacement including a suggestion that the decollement configuration controls the relative activity of kilauea s east and southwest rift zones

crustal motion in the southern andes 26 36 s do the andes behave like a microplate 1 a new global positioning system gps derived velocity field for the andes mountains 26 36 s allows analysis of instantaneous partitioning between elastic and anelastic deformation at the orogen s opposing sides adding an andes microplate to the traditional description of nazca south america plate convergence provides the kinematic framework for nearly complete explanation of the observed velocity field the results suggest the oceanic nazca boundary is fully locked while the continental backarc boundary creeps continuously at 4 5 mm yr the excellent fit of model to data 1 7 mm yr rms velocity misfit and the relative aseismicity of the upper crust in the interior andean region in comparison with its boundaries supports the notion that the mountains are not currently accruing significant permanent strains additionally the model implies permanent deformation is not accumulating throughout the backarc contractional wedge but rather that the deformation is accommodated only within a narrow deformational zone in the backarc

mechanism of the 1975 kalapana hawaii earthquake inferred from tsunami data we investigated the source mechanism of the 1975 kalapana hawaii earthquake m s 7 2 by modeling the tsunamis observed at three tide gauge stations hilo kahului and honolulu we computed synthetic tsunamis for various fault models the arrival times and the amplitudes of the synthetic tsunamis computed for ando s fault model fault length 40 km fault width 20 km strike n70 e dip 20 se rake 90 fault depth 10 km and slip 5 6 m are 10 min earlier and 5 times smaller than those of the observed respectively we tested fault models with different dip angles and depths models with a northwest dip direction yield larger tsunami amplitudes than those with a southeast dip direction models with shallower fault depths produce later first arrivals than deeper models we also considered the effects of the hilina fault system but its contribution to tsunami excitation is insignificant this suggests that another mechanism is required to explain the tsunamis one plausible model is a propagating slump model with a 1 m subsidence along the coast and a l m uplift offshore this model can explain the arrival times and the amplitudes of the observed tsunamis satisfactorily an alternative model is a wider fault model that dips 10 nw with its fault plane extending 25 km offshore well beyond the aftershock area of the kalapana earthquake these two models produce a similar uplift pattern offshore and kinematically are indistinguishable as far as tsunami excitation is concerned the total volume of displaced water is estimated to be 2 5 km 3 from the comparison of slump model and the single force model suggested earlier from seismological data we prefer a combination of faulting and large scale slumping on the south flank of kilauea volcano as the most appropriate model for the 1975 kalapana earthquake two basic mechanisms have been presented for explaining the deformation of the south flank of kilauea 1 pressure and density variation along the rift zone caused by magma injection and 2 gravitational instability due to the steep topography of the south flank of kilauea in either mechanism large displacements on the south flank are involved that are responsible for the observed large tsunamis

earthquakes triggered by silent slip events on kīlauea volcano hawaii silent earthquakes or slow slip events happen so slowly that they cannot be detected by conventional seismographs but now the availability of gps data means that they can be followed in detail segall et al report on a series of small earthquakes on k lauea volcano hawaii the earthquakes at depths of 7 to 8 km appear to be triggered by otherwise silent slow slips though small triggered earthquakes of this type can be used to help quantify the increased hazard due to slow slip events slow slip events or silent earthquakes have recently been discovered in a number of subduction zones including the nankai trough1 2 3 in japan cascadia4 5 and guerrero6 in mexico but the depths of these events have been difficult to determine from surface deformation measurements although it is assumed that these silent earthquakes are located along the plate megathrust this has not been proved slow slip in some subduction zones is associated with non volcanic tremor7 8 but tremor is difficult to locate and may be distributed over a broad depth range9 except for some events on the san andreas fault10 slow slip events have not yet been associated with high frequency earthquakes which are easily located here we report on swarms of high frequency earthquakes that accompany otherwise silent slips on k lauea volcano hawaii for the most energetic event in january 2005 the slow slip began before the increase in seismicity the temporal evolution of earthquakes is well explained by increased stressing caused by slow slip implying that the earthquakes are triggered the earthquakes located at depths of 7 8 km constrain the slow slip to be at comparable depths because they must fall in zones of positive coulomb stress change triggered earthquakes accompanying slow slip events elsewhere might go undetected if background seismicity rates are low detection of such events would help constrain the depth of slow slip and could lead to a method for quantifying the increased hazard during slow slip events because triggered events have the potential to grow into destructive earthquakes

inverting geodetic time series with a principal component analysis based inversion method 1 the global positioning system gps system now makes it possible to monitor deformation of the earth s surface along plate boundaries with unprecedented accuracy in theory the spatiotemporal evolution of slip on the plate boundary at depth associated with either seismic or aseismic slip can be inferred from these measurements through some inversion procedure based on the theory of dislocations in an elastic half space we describe and test a principal component analysis based inversion method pcaim an inversion strategy that relies on principal component analysis of the surface displacement time series we prove that the fault slip history can be recovered from the inversion of each principal component because pcaim does not require externally imposed temporal filtering it can deal with any kind of time variation of fault slip we test the approach by applying the technique to synthetic geodetic time series to show that a complicated slip history combining coseismic postseismic and nonstationary interseismic slip can be retrieved from this approach pcaim produces slip models comparable to those obtained from standard inversion techniques with less computational complexity we also compare an afterslip model derived from the pcaim inversion of postseismic displacements following the 2005 8 6 nias earthquake with another solution obtained from the extended network inversion filter enif we introduce several extensions of the algorithm to allow statistically rigorous integration of multiple data sources e g both gps and interferometric synthetic aperture radar time series over multiple timescales pcaim can be generalized to any linear inversion algorithm

from geodetic imaging of seismic and aseismic fault slip to dynamic modeling of the seismic cycle understanding the partitioning of seismic and aseismic fault slip is central to seismotectonics as it ultimately determines the seismic potential of faults thanks to advances in tectonic geodesy it is now possible to develop kinematic models of the spatiotemporal evolution of slip over the seismic cycle and to determine the budget of seismic and aseismic slip studies of subduction zones and continental faults have shown that aseismic creep is common and sometimes prevalent within the seismogenic depth range interseismic coupling is generally observed to be spatially heterogeneous defining locked patches of stress accumulation to be released in future earthquakes or aseismic transients surrounded by creeping areas clay rich tectonites high temperature and elevated pore fluid pressure seem to be key factors promoting aseismic creep the generally logarithmic time evolution of afterslip is a distinctive feature of creeping faults that suggests a logarithmic dependency of fault friction on slip rate

an increment of volcano collapse kinematics of the 1975 kalapana hawaii earthquake abstract earthquakes on kilauea volcano can be even more hazardous than volcanic eruptions the most recent large earthquake was the 1975 kalapana quake which generated a locally damaging tsunami over 8 m of horizontal surface displacement and 3 m of subsidence at the coast the seismic magnitude estimates for this earthquake range from ms 7 2 to mw 7 7 the complexity of the seismic data is matched by the complex surface deformation that was observed with a combination of leveling tilt and electronic distance measurements edm the geodetic data shows evidence of collapse of the caldera due to magma withdrawal localized subsidence over the rift zone normal faulting associated with shallow slumping along the hilina pali and displacement on the basal detachment fault we model the summit caldera rift zones and basal detachment fault as dislocations in an elastic half space to better quantify the sources of deformation related to this devastating event using inversion techniques that allow us to let the data constrain the geometry and magnitude of these sources we find 0 04 km3 of magma withdrawal at the summit between 3 5 m of opening along the rift zone and 7 1 m of slip along the basal detachment at 8 3 km depth models that allow for finer spatial resolution of slip on the detachment show that the largest slip occurred west of the seismically recorded earthquake hypocenter along the region of the fault below the coast and that the majority of the fault slip occurred south of the region of microseismicity our best fitting model has a geodetic moment of 4 1 1020 nm mw 7 7 which is consistent with tsunami models and recent analysis of long period seismic data the residual displacements of sites located in the hanging wall of the hilina pali slump system which were not used in the inversion of the other sources suggest that triggered shallow slumping contributed several meters of horizontal displacement to the coseismic displacement at these sites

magmatically triggered slow slip at kilauea volcano hawaii we demonstrate that a recent dike intrusion probably triggered a slow fault slip event sse on kilauea volcano9s mobile south flank our analysis combined models of advanced land observing satellite interferometric dike intrusion displacement maps with continuous global positioning system gps displacement vectors to show that deformation nearly identical to four previous sses at kilauea occurred at far field sites shortly after the intrusion we model stress changes because of both secular deformation and the intrusion and find that both would increase the coulomb failure stress on possible sse slip surfaces by roughly the same amount these results in concert with the observation that none of the previous sses at kilauea was directly preceded by intrusions but rather occurred during times of normal background deformation suggest that both extrinsic intrusion triggering and intrinsic secular fault creep fault processes can lead to sses

tectonic controls of phanerozoic sedimentary rock cycling the cyclic nature of the phanerozoic sedimentary rock distribution carbon sulphur coupling and material transfer among sedimentary reservoirs appears to be controlled by tectonic factors the distribution of preserved sedimentary mass in terms of rock mass remaining u geologic age shows a minimum c 300 350 ma ago which separates two subcycles of erosion and deposition of sedimentary rocks the older subcycle was interrupted because of the major continental collisions of the devonian and late carboniferous these collisions resulted in a reduction of outcrop areas of rocks of the older cycle relative to their masses leading to a decline in the probability of destruction and an increase in half life of these older sediments a strong correlation exists between the long term cyclicity in the phanerozoic global sea level curve and the distribution of carbon and sulphur among their major exogenic reservoirs this correlation is related to two principal tectonic modes of the phanerozoic oscillatory and submergent it is postulated that the submergent mode of active plate convergence obduction and subduction of sediments large ridge volume and high sea level gave rise to low erosion and sedimentation rates less restricted environments of carbonate deposition and relatively high atmospheric co 2 levels high temperatures resulting from an increased rate of production of co 2 from diagenetic and metamorphic reactions at subduction zones as sea level rose carbon was transferred from the sedimentary reservoir of reduced organic carbon to that of oxidized inorganic carbon in limestones whereas sulphur moved from the oxidized sulphate reservoir to the reduced sulphide reservoir as sea level fell reservoir transfers were opposite to those above culminating in the oscillatory mode of generally elevated continental interiors these reservoir transfers are consistent with secular changes observed in the distribution of 13 c and 34 s in sedimentary materials during the phanerozoic petrographic examination of phanerozoic oolite formations shows that ooids with preserved calcitic relict textures are characteristic of pre carboniferous carbonate rocks whereas ooids with relict textures indicative of initial aragonite mineralogy are dominant in rocks of younger age these changes in ooid mineralogy may be interpreted as reflecting changes in co 2 levels of the ocean atmosphere system consistent with the above tectonic considerations atmospheric co 2 levels were higher prior to carboniferous time favouring formation of calcitic ooids and skeletal parts after the carboniferous co 2 levels fell and aragonite and mg calcites of greater than 8 mol mg increased in abundance as precipitates

fault frictional parameters and material properties revealed by slow slip events at kilauea volcano hawai i 1 we categorize slow slip events at kilauea volcano into two distinct families based on gps measurements of the surface displacement patterns an event correlation filter confirms that eastern and western families are statistically distinguishable with the western family notably self similar the western family exhibits quasi periodicity with regular repeat times while eastern family events are aperiodic or have complicated periodicity if the decollement is the source fault for both families of events it must have varying frictional properties at the 10 km scale of separation the temporal slip and spatial scaling behavior are consistent with a simplistic rate and state dependent frictional formalism provided that the characteristic slip distance for state evolution dc is of the order of millimeters rather than the 10 100 m typically found in lab studies and the shear rigidity is around 2 gpa consistent with fault gouge material

pressurized magma reservoir within the east rift zone of kīlauea volcano hawai i evidence for relaxed stress changes from the 1975 kalapana earthquake we use 2000 2012 insar data from multiple satellites to investigate magma storage in k lauea s east rift zone erz the study period includes a surge in magma supply rate and intrusion eruptions in 2007 and 2011 the kupaianaha area inflated by 5 cm prior to the 2007 intrusion and the n pau crater area by 10 cm following the 2011 intrusion for the n pau crater area elastic modeling suggests an inflation source at 5 2 km depth or more below sea level the reservoir is located in the deeper section of the rift zone for which secular magma intrusion was inferred for the period following the 1975 mw7 7 decollement earthquake reservoir pressurization suggests that in this section of the erz extensional stress changes due to the earthquake have largely been compensated for and that this section is approaching its pre 1975 state reservoir pressurization also puts the molten core model into question for this section of k lauea s rift zone

inversion of strong ground motion and teleseismic waveform data for the fault rupture history of the 1979 imperial valley california earthquake a least squares point by point inversion of strong ground motion and teleseismic body waves is used to infer the fault rupture history of the 1979 imperial valley california earthquake the imperial fault is represented by a plane embedded in a half space where the elastic properties vary with depth the inversion yields both the spatial and temporal variations in dislocation on the fault plane for both right lateral strike slip and normal dip slip components of motion inversions are run for different fault dips and for both constant and variable rupture velocity models effects of different data sets are also investigated inversions are compared which use the strong ground motions alone the teleseismic body waves alone and simultaneously the strong ground motion and teleseismic records the inversions are stabilized by adding both smoothing and positivity constraints the moment is estimated to be 5 0 10 25 dyne cm and the fault dip 90 5 dislocation in the hypocentral region south of the united states mexican border is relatively small and almost dies out near the border dislocation then increases sharply north of the border to a maximum of about 2 m under interstate 8 dipslip motion is minor compared to strike slip motion and is concentrated in the sediments the best fitting constant rupture velocity is 80 per cent of the local shear wave velocity however there is a suggestion that the rupture front accelerated from the hypocenter northward the 1979 imperial valley earthquake can be characterized as a magnitude 5 earthquake at the hypocenter which then grew into or triggered a magnitude 6 earthquake north of the border

rupture in the 4 may 2018 m w 6 9 earthquake seaward of the kilauea east rift zone fissure eruption in hawaii

dating fluvial terraces with and profiles application to the wind river wyoming abstract fluvial strath terraces provide a record of river incision and the timing of climatic perturbations to the fluvial system dating depositional surfaces like terraces that are older than the range of 14 c however is difficult we employ a cosmogenic radionuclide crn profile technique that addresses a major problem of crn dating on such surfaces nuclide inheritance by measuring 10 be and 26 al profiles we constrain the exposure age and the mean crn inheritance for the deposit the crn profile also yields a self check on the assumptions underlying the method we report our attempts to date terraces along the wind river wy like many sequences of western north american fluvial terraces these are inferred to reflect oscillation between glacial and interglacial conditions in the headwaters previous dating of some of these terraces and the associated terraces and glacial deposits makes this a unique location to compare dating methods dates from five sites along the bull lake glacial correlative terrace wr 3 are 118 125 ka which agrees with dates on bull lake age moraines and independent age estimates on the terrace and is consistent with the model of terrace glacial relationship crn inheritance is significant and highly variable requiring it be considered despite the additional sampling complexity assuming all inheritance in wr 3 deposits arises during exhumation in the headwaters we obtain minimum mean rates of exhumation of 13 130 m my for the source rocks alternatively assuming the crns are inherited during clast transport the time of fluvial transport from source to terrace is 10 ka it increases downstream and is lower for sand than cobbles the crn ages for older terraces wr 7 300 ka and wr 15 510 ka are lower by 50 than previous estimates based on tephrochronology the most plausible explanation is eolian deflation of a once thicker loess cover on the terrace surfaces mean thicknesses of loess of 0 5 1 5 m are required to reconcile these concentrations of crn with the previous estimates of age difficulty in dating the older terraces emphasizes that geologic caution independent estimates of age and multiple sample sites should still be part of dating depositional surfaces with crns even when employing the inheritance correction technique