dike intrusion associated with the 2000 eruption of miyakejima volcano japan the model for the 2000 dike intrusion event between kozushima and miyakejima volcano japan was reinvestigated after the sudden earthquake swarm in miyakejima volcano a dike intrusion of large volume was detected by the nationwide gps network geonet the displacements detected with gps stations over an area with a radius of about 200 km shows a distribution that is consistent with the dike source being located near miyakejima volcano

propagation of rifting along the arabia somalia plate boundary into afar it is generally accepted that the aden ridge has propagated westward from 58 e to the western tip of the gulf of aden tadjoura at the edge of afar here we use new tectonic and geochronological data to examine the geometry and kinematics of deformation related to the penetration of that ridge on dry land in the republic of djibouti we show that it veers northward forming a narrow zone of dense faulting along the northeastern edge of the afar depression the zone includes two volcanic rifts asal ghoubbet and manda inakir connected to one another and to the submarine part of the ridge by transfer zones both rifts are composite divided into two or three disconnected parallel nw se striking subrifts all of which appear to have propagated northwestward in asal ghoubbet as in manda inakir the subrifts appear to have formed in succession through north directed jumps from subrifts more farther south at present the northernmost subrifts manda and dirko koma of the manda inakir rift form the current tip of the northward propagating arabia somalia plate boundary in afar we account for most observations by a mechanical model similar to that previously inferred for the gulf of aden in which propagation is governed by the intensity and direction of the minimum horizontal principal stress 3 we interpret the northward propagation on land almost orthogonal to that in the gulf to be related to necking of the central afar lithosphere where it is thinnest such necking may be a consequence of differential magmatic thickening greater in the center of the afar depression where the ethiopian hot spot enhanced profuse basaltic effusion and underplating than along the edges of the depression the model explains why the aden ridge foregoes its wsw propagation direction constant from 58 e to asal ghoubbet at a smaller scale individual rifts and subrifts keep opening perpendicular to the arabia somalia or danakil somalia motion vector and propagate northwestward concurrently such lithospheric cracks are forced to jump northward such that the plate boundary remains inside the regional n s necking zone changes of obliquity between the directions of overall and local propagation may account for different segmentation patterns a small angle promoting long en echelon subrifts and a high angle smaller nested subrifts within subrifts the propagation mechanism is thus similar whether in oceanic or continental lithosphere the principal change being the overall propagation path here governed by thickness changes rather than by the geometry in map view as previously inferred for the rest of the aden ridge finally because the same mechanism has led rifting along the red sea to propagate southward and jump to the western edge of afar the arabia somalia and arabia nubia plate boundaries tips have missed each other and keep overlapping further leading to strain transfer by large scale bookshelf faulting

a sea floor spreading event captured by seismometers two thirds of earth s surface is formed at mid ocean ridges yet sea floor spreading events are poorly understood because they occur far beneath the ocean surface at 9 50 n on the east pacific rise ocean bottom seismometers recently recorded the microearthquake character of a mid ocean ridge eruption including precursory activity a gradual ramp up in activity rates since seismic monitoring began at this site in october 2003 suggests that eruptions may be forecast in the fast spreading environment the pattern culminates in an intense but brief 6 hour inferred diking event on 22 january 2006 followed by rapid tapering to markedly decreased levels of seismicity

segmented lateral dyke growth in a rifting event at barðarbunga volcanic system iceland crust at many divergent plate boundaries forms primarily by the injection of vertical sheet like dykes some tens of kilometres long1 previous models of rifting events indicate either lateral dyke growth away from a feeding source with propagation rates decreasing as the dyke lengthens2 3 4 or magma flowing vertically into dykes from an underlying source5 6 with the role of topography on the evolution of lateral dykes not clear here we show how a recent segmented dyke intrusion in the bar arbunga volcanic system grew laterally for more than 45 kilometres at a variable rate with topography influencing the direction of propagation barriers at the ends of each segment were overcome by the build up of pressure in the dyke end then a new segment formed and dyke lengthening temporarily peaked the dyke evolution which occurred primarily over 14 days was revealed by propagating seismicity ground deformation mapped by global positioning system gps interferometric analysis of satellite radar images insar and graben formation the strike of the dyke segments varies from an initially radial direction away from the bar arbunga caldera towards alignment with that expected from regional stress at the distal end a model minimizing the combined strain and gravitational potential energy explains the propagation path dyke opening and seismicity focused at the most distal segment at any given time and were simultaneous with magma source deflation and slow collapse at the bar arbunga caldera accompanied by a series of magnitude m 5 earthquakes dyke growth was slowed down by an effusive fissure eruption near the end of the dyke lateral dyke growth with segment barrier breaking by pressure build up in the dyke distal end explains how focused upwelling of magma under central volcanoes is effectively redistributed over long distances to create new upper crust at divergent plate boundaries

piezometric and hydrogeochemical characterization of groundwater circulation in complex karst aquifers a case study the mancha real pegalajar aquifer southern spain a small karst aquifer of great structural complexity has been subjected to significant resource withdrawal over recent decades this exploitation aroused social conflict due to the effect it has had on emblematic springs this research has analysed piezometric data collected over the course of 12 years and the spatial hydrochemical data supplied by the main water points associated with it the spatial and temporal evolution of the main chemical species in the groundwater and the hydrogeochemical processes affecting them have been studied modelling them with the programme phreeqc these data suggest a complicated model of hydrogeological function with sectors storing water at different depths and connected to each other locally as determined by the geological structure

the volcano seismic crisis in afar ethiopia starting september 2005 we describe the seismo volcanic crisis that occurred in northern afar in late 2005 which involved 15 earthquakes greater than m5 and a small explosive silicic eruption from a vent called da ure at 12 651 n 40 519 n close to dabbahu volcano the purpose is to pull together the different sources of information into a coherent preliminary interpretation of what happened the main geophysical data are the locations of the largest earthquakes and a radar interferogram of unusually high quality that reveals injection of a 60 km long dyke with surface deformation expressed as normal faulting subsidence occurred around the dabbahu volcanic edifice most of the dyke is likely to have been basaltic rather than silicic although the eruption was silicic the volume of the subsidence represents at most 25 of the magma injected into the dyke the silicic eruption was possibly triggered by interaction of basaltic magma with a shallow silicic reservoir at about the same time as the eruption and dyking episode some activity appears also to have taken place at the lava lake at erta ale volcano some 150 km to the north of the eruption site we evaluate the possibility that there may be some link by calculating stresses associated with opening of the fissure and looking at the activity of the lava lake as revealed by the thermal anomaly seen by weather satellites

controlling factors on earthquake swarms associated with magmatic intrusions constraints from iceland abstract the currently two most important surveying tools in volcano monitoring are observations of seismicity and surface deformation magma migration within the earth s crust is frequently associated with seismic activity often occurring as distinct earthquake swarms understanding the exact nature and significance of these swarms is important when evaluating volcanic crisis situations there seems however to be no general correlation between the amount of seismic energy release and the rate and volume of magma on the move which may complicate an immediate risk evaluation it has previously been shown that stressing rate appears as the controlling factor on the occurrence of seismicity during intrusion however in this paper we emphasize that this is only true under specific circumstances other factors influence the evolution and resulting seismic energy release and each intrusion scenario needs to be evaluated separately three icelandic cases where seismic swarms have been confirmed through deformation measurements to be related to magmatic movements were selected we show how the relationship between volume change and resulting seismicity varies greatly between cases and suggest that in a general approach the most important factor governing the level of seismic energy release is the background stress state tectonic setting regional stresses and tectonic history are therefore all of great importance when evaluating a magmatically induced seismic crisis the stressing rate only plays the dominant role when either a stress in the intruded area is close to the yield strength before intrusion or possibly b when the intrusion is extraordinarily voluminous and rapid

the june july 1993 seismo acoustic event at coaxial segment juan de fuca ridge evidence for a lateral dike injection the coaxial segment juan de fuca ridge was the site of an intense swarm of earthquakes that began at 21 43 gmt on june 26 1993 julian day 177 the swarm started near 46 15 n and migrated northward over the next 40 hours to 46 36 n where the majority of 676 events occurred during the following 3 weeks of activity the earthquakes propagated nne at a velocity of 0 3 0 1 m s 1 from the southern to northern swarm sites the activity went undetected by land based seismic networks along the oregon and washington coasts suggesting that the earthquakes were all m 4 0 the character of this earthquake swarm is very similar to dike injections observed at krafla and kilauea volcanoes the earthquake activity and subsequent migration probably represent a lateral dike injection into faults and fissures comprising the coaxial segment the reservoir acting as the dike s source likely resides beneath or to the south of the initial southern swarm of earthquakes the large magma supply at axial volcano cannot be ruled out as the source for the coaxial dike even though axial volcano exhibited no earthquake activity related to the coaxial swarm the t wave earthquake swarm reported here is the first deep ocean observation of volcanic seismicity associated with what most likely was a mid ocean ridge dike injection event

stress transfer between thirteen successive dyke intrusions in ethiopia large earthquakes are known to trigger subsequent earthquakes in nearby regions but similar triggering has not been confirmed for volcanic hazards analysis of the progressive deformation associated with 13 volcanic dykes intruded in ethiopia between 2005 and 2009 indicates that magmatic intrusions can help to trigger subsequent intrusions

a reinterpretation of seismicity associated with the january 1983 dike intrusion at kilauea volcano hawaii in january 1983 a dike intrusion fissure eruption generated a swarm of 375 magnitude 1 to 3 earthquakes along a 16 km segment of kilauea s middle east rift zone we searched the hawaiian volcano observatory catalog for multiplets of similar events from this region from 1980 through 1985 and obtained precise relative locations by waveform cross correlation over 150 of the intrusion earthquakes could be grouped into 14 multiplets of five or more events with sufficient similarity for accurate relocation some multiplets were active for only a few minutes during the downrift migration phase of the seismic swarm consistent with generation near the propagating dike tip while others were active for several days the two multiplets nearest the origin of the seismic swarm include events from the preceding days and months most multiplets span only 50 to 100 m following relocation are located at about 3 to 4 km depth and appear to deepen downrift the catalog depths of those earthquakes in multiplets and those not in multiplets are similar suggesting that most of the recorded seismicity may have come from a very limited depth interval despite the fact that the dike breached the surface by analogy with a mechanical model used to explain a similar clustering of background seismicity in the upper east rift in 1991 we infer that the earthquakes are generated in regions of high stress concentration immediately above kilauea s deforming deep rift body this conclusion is consistent with the depth of the top of the deep rift body inferred from geodetic data and with numerical calculations suggesting that a significant ambient differential stress is required for dikes to produce earthquakes larger than magnitude 1

the quantum event of oceanic crustal accretion impacts of diking at mid ocean ridges seafloor diking eruptive events represent the irreducible quantum events of upper oceanic crustal accretion they record events by which a large portion of the oceanic crust has formed through geological history since 1993 the u s navy s real time sound surveillance system has allowed location of ongoing acoustic signatures of dike emplacement and basalt eruptions at ridge crests in the northeast pacific these diking eruptive events trigger a sequence of related rapidly evolving physical chemical and biological processes magmatic volatiles released during these events may provide nutrients for communities of subseafloor microorganisms some of which thrive in high temperature anaerobic environments many of the organisms identified from these systems are archaea if microorganisms can thrive in the water saturated pores and cracks within deep volcanically active portions of our planet other hydrothermally active planets may harbor similar life forms

satellite radar interferometry 1993 1999 suggests deep accumulation of magma near the crust mantle boundary at the krafla volcanic system iceland 1 deep magma accumulation near the crust mantle boundary 21 km depth at the krafla volcanic system is suggested from insar observations a best fit model derived from four interferograms covering 1993 1999 comprises an opening dike representing plate spreading and post rifting deformation and two mogi sources a mogi source deflating at a rate of 0 3 106 m3 yr coincides with the shallow krafla magma chamber while a deeper inflating mogi source further north at 21 km depth inflates at a rate of 26 106 m3 yr the inflating source is at or near the crust mantle boundary as identified by seismic studies and is interpreted as accumulating magma

modes of faulting at mid ocean ridges abyssal hill bounding faults that pervade the oceanic crust are the most common tectonic feature on the surface of the earth the recognition that these faults form at plate spreading centres came with the plate tectonic revolution recent observations reveal a large range of fault sizes and orientations numerical models of plate separation dyke intrusion and faulting require at least two distinct mechanisms of fault formation at ridges to explain these observations plate unbending with distance from the top of an axial high reproduces the observed dip directions and offsets of faults formed at fast spreading centres conversely plate stretching with differing amounts of constant rate magmatic dyke intrusion can explain the great variety of fault offset seen at slow spreading ridges very large offset normal faults only form when about half the plate separation at a ridge is accommodated by dyke intrusion

why did arabia separate from africa insights from 3 d laboratory experiments abstract we have performed 3 d scaled lithospheric experiments to investigate the role of the gravitational force exerted by a subducting slab on the deformation of the subducting plate itself experiments have been constructed using a dense silicone putty plate to simulate a thin viscous lithosphere floating in the middle of a large box filled with glucose syrup simulating the upper mantle we examine three different plate configurations i subduction of a uniform oceanic plate ii subduction of an oceanic continental plate system and iii subduction of a more complex oceanic continental system simulating the asymmetric africa eurasia system each model has been performed with and without the presence of a circular weak zone inside the subducting plate to test the near surface weakening effect of a plume activity our results show that a subducting plate can deform in its interior only if the force distribution varies laterally along the subduction zone i e by the asymmetrical entrance of continental material along the trench in particular extensional deformation of the plate occurs when a portion of the subduction zone is locked by the collisional process the results of this study can be used to analyze the formation of the arabian plate we found that intraplate stresses similar to those that generated the africa arabia break up can be related to the neogene evolution of the northern convergent margin of the african plate where a lateral change from collision mediterranean and bitlis to active subduction makran has been described second intraplate stress and strain localization are favored by the presence of a weakness zone such as the one generated by the afar plume producing a pattern of extensional deformation belts resembling the red sea gulf of aden rift system

brittle creep damage as the seismic signature of dyke propagations within basaltic volcanoes contemporary to nine dyke intrusions on piton de la fournaise etna and miyakejima volcanoes we recover stationary seismicity rate and energy release over time whether the dyke reaches the surface or not this generic seismicity pattern for the dyke propagation of low viscosity magma argues for the fluid driven crack propagation to be a scale independent stationary process this prevents any prediction of the time to eruption during the dyke propagation phase using seismicity rate alone the seismic signature of the volcano deformation triggered by dyke injections corresponds to brittle creep damage in a strain driven setting whether mechanical properties of host rock structure or geometrical effects influence this generic stationary response is not resolved by the seismic data because a few if any aftershocks are resolved contemporary to dyke intrusion the seismicity is purely driven by the dyke dynamics that is a proxy for the dyke volumetric growth

evidence of a recent magma dike intrusion at the slow spreading lucky strike segment mid atlantic ridge 1 mid ocean ridge volcanic activity is the fundamental process for creation of ocean crust yet the dynamics of magma emplacement along the slow spreading mid atlantic ridge mar are largely unknown we present acoustical seismological and biological evidence of a magmatic dike intrusion at the lucky strike segment the first detected from the deeper sections 1500 m of the mar the dike caused the largest teleseismic earthquake swarm recorded at lucky strike in 20 years of seismic monitoring and one of the largest ever recorded on the northern mar hydrophone records indicate that the rate of earthquake activity decays in a nontectonic manner and that the onset of the swarm was accompanied by 30 min of broadband 3 hz intrusion tremor suggesting a volcanic origin two submersible investigations of high temperature vents located at the summit of lucky strike seamount 3 months and 1 year after the swarm showed a significant increase in microbial activity and diffuse venting this magmatic episode may represent one form of volcanism along the mar where highly focused pockets of magma are intruded sporadically into the shallow ocean crust beneath long lived discrete volcanic structures recharging preexisting seafloor hydrothermal vents and ecosystems index terms 3035 marine geology and geophysics midocean ridge processes 7280 seismology volcano seismology 8419 8149 tectonophysics planetary tectonics 5475 4259 oceanography general ocean acoustics 9325 information related to geographic region atlantic ocean keywords mid atlantic ridge earthquake hydroacoustic

transient deformation in the asal ghoubbet rift djibouti since the 1978 diking event is deformation controlled by magma supply rates the asal ghoubbet rift ag rift in djibouti lies in the subaerial continuation of the aden ridge system thereby constituting a unique location to study rifting processes and mechanisms involved in continental breakup and oceanic spreading continually upgraded and expanded geodetic technology has been used to record the 1978 asal rifting event and postdiking deformation in light of recent results obtained for the manda hararo dabbahu rifting event 2005 2010 we propose that the horizontal and vertical geodetic data can be modeled with a double source involving a dike like inflation component aligned along the rift axis and a spherical pressure source located at midsegment below the fieale caldera by revisiting the codiking data we propose that the reservoir below fieale could have fed at least partially the 1978 injection and the contemporaneous ardoukoba eruption and potentially induced local subsidence due to magma draining out of the central reservoir as an alternative to previously proposed viscoelastic relaxation models we reinterpret postdiking observations using a purely elastic rheology we determine the relative contribution of a midsegment reservoir inflation and a dike like opening component together with their respective time evolutions our results suggest that interactions between steadily accumulating tectonic strain and temporal variations in melt supply to the shallow magma plumbing system below the ag rift may entirely explain the geodetic observations and that viscoelastic deformation processes played a minor role in the 30 years following the 1978 rifting event

seafloor spreading event in western gulf of aden during the november 2010 march 2011 period captured by regional seismic networks evidence for diking events and interactions with a nascent transform zone in november 2010 intense seismic activity including 29 events with a magnitude above 5 0 started in the western part of the gulf of aden where the structure of the oceanic spreading ridge is characterized by a series of n115 trending slow spreading segments set within an ew trending rift using signals recorded by permanent and temporary networks in djibouti and yemen we located 1122 earthquakes with a magnitude ranging from 2 1 to 5 6 from 01 november 2010 to 31 march 2011 by looking in detail at the space time distribution of the overall seismicity and both the frequency and the moment tensor of large earthquakes we reexamine the chronology of this episode in addition we also interpret the origin of the activity using high resolution bathymetric data as well as from observations of sea floor cable damage caused by high temperatures and lava flows the analysis allows us to identify distinct active areas firstly we interpret that this episode is mainly related to a diking event along a specific ridge segment located at e044 in light of previous diking episodes in nearby subaerial rift segments for which field constraints and both seismic and geodetic data exist we interpret the space time evolution of the seismicity of the first few days migration of earthquakes suggests initial magma ascent below the segment center this is followed by a southeastward dike propagation below the rift immediately followed by a northwestward dike propagation below the rift ending below the northern ridge wall the cumulative seismic moment associated with this sequence reaches 9 1 1017 nm and taking into account a very low seismic versus geodetic moment we estimate an horizontal opening of 0 58 to 2 9 m the seismic activity that followed occurred through several bursts of earthquakes aligned along the segment axis which are interpreted as short dike intrusions implying fast replenishment of the crustal magma reservoir feeding the dikes over the whole period the opening is estimated to 1 76 8 8 m across the segment a striking feature of this episode is that the seismicity remained confined within one individual segment whereas the adjacent en echelon segments were totally quiescent suggesting that the magma supply system of one segment is disconnected from those of the neighboring segments second we identify activity induced by the first intrusion with epicenters aligned along a n035 e trending 30 km long at the northwestern end of the active opening segment this group encompasses more than 7 earthquakes with magnitude larger than 5 0 and with strike slip focal mechanisms consistent with the faults identified in the bathymetry and the structural pattern of the area we propose that a transform fault is currently in formation which indicates an early stage of the ridge segmentation at the locus of the trend change of the spreading ridge which also corresponds to the boundary between a clear oceanic lithosphere and the zone of transform between continental and oceanic crust

time dependent dike propagation from joint inversion of seismicity and deformation data dike intrusions both deform earth s surface and induce propagating earthquake swarms we develop methods to utilize both deformation and seismicity from brittle volcano tectonic earthquakes to image time dependent dike propagation dieterich s 1994 seismicity rate theory is used to relate dike induced stress changes to seismicity rate and is combined with elastic green s functions relating dike opening to deformation different space time patterns of seismicity develop if earthquakes occur at the same depth as the dike compared to above below the dike in the former seismicity initiates near the dike s leading edges but shuts off as the dike tips pass and seismogenic volumes fall into stress shadows in the latter seismicity continues at a decaying rate after the tips pass we focus on lateral propagation and develop a nonlinear inversion method that estimates dike length and pressure as a function of time the method is applied to the 2007 father s day intrusion in kilauea volcano seismicity is concentrated at 3 km depth comparable to geodetic estimates of dike depth and decays rapidly in time with lateral propagation only and a vertical dike tip line it is difficult to fit both gps data and the rapid down rift jump in seismicity suggesting significant vertical propagation for the events to have occurred below the dike requires a very short aftershock decay time hence unreasonably high background stressing rate the rapid decay is better explained if the dike extends somewhat below the seismicity we suggest that joint inversion is useful for studying the diking process and may allow for improved short term eruption forecasts

evidence that coupling to magma chambers controls the volume history and velocity of laterally propagating intrusions 1 in this paper i present a simple analytical solution for the unsolved problem of a propagating dike coupled to a magma chamber as recognized by segall et al 2001 the flow of the magma from a chamber to a dike is described by an ordinary differential equation for the unknown pressures of chamber and dike the most intuitive assumption on magma exchange is that the volume gained by the intrusion equals the volume lost by the chamber this constraint however implies a chamber so large that its pressure does not drop during the feeding process and an incompressible magma that does not expand or vesiculate as it intrudes models assuming constant driving pressure fail to explain observations here mass conservation constrains pressure as magma flows from the chamber to the dike and the compressibilities of magma and sources control volume changes these assumptions allow me to decouple the equations and solve the system analytically the model predicts that chamber and dike volume change exponentially with time as v t v 1 exp t the asymptotic volume v and the timescale are found to be functions of rock magma chamber and dike parameters and of the initial pressure conditions assuming that the intrusion is shaped as an elongating cuboid its velocity is found to change as v v0exp t where v0 is the initial dike velocity the model is successfully tested against the best observations available for lateral dike propagation events the volume history of chamber and dike during the 2000 miyakejima intrusion japan and dike velocity during the 1978 krafla event iceland and during some intrusions following the 2005 event in afar ethiopia this paper confirms and extends the results of a previous study by rivalta and segall 2008 who found the final ratio between dike volume and the volume withdrawn from the chamber to be rv 1 4 m 3 1 where is the host rock rigidity and m is the magma compressibility here i demonstrate that the formula for rv holds at any time during the intrusion not just at the end this model confirms that some magma chambers behave as stiff magma tanks able to inflate large dikes as balloons and demonstrates that high volume ratios are unlikely if the chamber has a pure sill shape fitting observations with the model allows us to estimate several dike chamber magma and rock parameters

capturing magma intrusion and faulting processes during continental rupture seismicity of the dabbahu afar rift continental rupture models emphasize the role of faults in extensional strain accommodation extension by dyke intrusion is commonly overlooked a major rifting episode that began in 2005 september in the afar depression of ethiopia provides an opportunity to examine strain accommodation in a zone of incipient plate rupture earthquakes recorded on a temporary seismic array 2005 october to 2006 april direct observation of fault patterns and geodetic data document ongoing strain and continued dyke intrusion along the 60 km long dabbahu rift segment defined in earlier remote sensing studies epicentral locations lie along a 3 km wide 50 km long swath that curves into the se flank of dabbahu volcano a second strand continues to the north toward gab ho volcano considering the 8 m of opening in the september crisis we interpret the depth distribution of microseismicity as the dyke intrusion zone the dykes rise from 10 km to the near surface along the 60 km long length of the tectono magmatic segment focal mechanisms indicate slip along nnw striking normal faults perpendicular to the arabia nubia plate opening vector the seismicity insar continuous gps and structural patterns all suggest that magma injection from lower or subcrustal magma reservoirs continued at least 3 months after the main episode persistent earthquake swarms at two sites on dabbahu volcano coincide with areas of deformation identified in the insar data 1 an elliptical northwestward dipping zone of seismicity and subsidence interpreted as a magma conduit and 2 a more diffuse 8 km radius zone of shallow seismicity 2 km above a shadow zone interpreted as a magma chamber between 2 5 and 6 km subsurface insar and continuous gps data show uplift above a shallow source in zone 2 and uplift above the largely aseismic gab ho volcano the patterns of seismicity provide a 3 d perspective of magma feeding systems maintaining the along axis segmentation of this incipient seafloor spreading segment

constraints on dike propagation from continuous gps measurements the january 1997 east rift zone eruption on kilauea volcano hawaii occurred within a network of continuous global positioning system gps receivers the gps measurements reveal the temporal history of deformation during dike intrusion beginning 8 hours prior to the onset of the eruption the dike volume as a function of time estimated from the gps data using elastic green s functions for a homogeneous half space shows that only two thirds of the final dike volume accumulated prior to the eruption and the rate of volume change decreased with time these observations are inconsistent with simple models of dike propagation which predict accelerating dike volume up to the time of the eruption and little or no change thereafter deflationary tilt changes at kilauea summit mirror the inferred dike volume history suggesting that the rate of dike propagation is limited by flow of magma into the dike a simple lumped parameter model of a coupled dike magma chamber system shows that the tendency for a dike to end in an eruption rather than intrusion is favored by high initial dike pressures compressional stress states large compressible magma reservoirs and highly conductive conduits linking the dike and source reservoirs comparison of model predictions to the observed dike volume history the ratio of erupted to intruded magma and the deflationary history of the summit magma chamber suggest that most of the magma supplied to the growing dike came from sources near to the eruption through highly conductive conduits interpretation is complicated by the presence of multiple source reservoirs magma vesiculation and cooling as well as spatial variations in dike normal stress reinflation of the summit magma chamber following the eruption was measured by gps and accompanied a rise in the level of the pu u o o lava lake for a spheroidal chamber these data imply a summit magma chamber volume of 20 km3 consistent with recent estimates from seismic tomography continuous deformation measurements can be used to image the spatiotemporal evolution of propagating dikes and to reveal quantitative information about the volcanic plumbing systems

mechanical models for dikes a third school of thought abstract geological and geophysical data from continental volcanic centers and giant radial swarms and from oceanic shield volcanoes and rift zones indicate that dikes propagate laterally for distances that can be 10 to over 100 times their height what traps dikes within the shallow lithosphere and promotes these highly eccentric shapes gravity induced stress gradients in the surrounding rock and pressure gradients in the magma are the primary loading mechanisms pressure gradients due to magma flow are secondary to insignificant because the flow direction is dominantly horizontal this configuration of vertical blade shaped dikes with horizontal dike propagation and magma flow is fundamentally different from the two dike model configurations described in a recent review paper as two schools of thought for mechanical models of dikes in school i a dike is disconnected from its source and ascends under the influence of buoyancy in school ii a dike is connected to a magma reservoir and is driven upward by magma flux from the source we review the geological and geophysical data supporting the vertical dike horizontal flow propagation configuration and suggest the abundance and veracity of these data in many different geological settings and the modeling results that address this physical process warrant adding this as a third school of thought a new analytical solution for the boundary value problem of a homogeneous isotropic and linear elastic solid with a vertical fluid filled crack is used to investigate the effects of gravitationally induced stress and pressure gradients on the aperture distribution dike tip stress intensity and stable height model results indicate that in a homogeneous crust dikes can achieve stable heights greater than a kilometer only if the host rock fracture toughness k ic 100 mpa m 1 2 however density stratification of the crust is an effective mechanism for trapping kilometer scale dikes even if the host rock is very weak k ic 0 this analysis may explain why vertical dikes propagate laterally for great distances but reside within a narrow range of depths in the crust

strain transfer between disconnected propagating rifts in afar we showed before that both the aden and red sea plate boundaries are currently rifting and propagating along two distinct paths into afar through the opening of a series of disconnected propagating rifts here we use new geochronological tectonic and paleomagnetic data that we acquired mostly in the southeastern part of afar to examine the geometry kinematics and time space evolution of faulting related to strain transfer processes it appears that transfer of strain is accommodated by a bookshelf faulting mechanism wherever rifts or plate boundaries happen to overlap without being connected this mechanism implies the rotation about a vertical axis of small rigid blocks along rift parallel faults that are shown to slip with a left lateral component which is as important as their normal component of slip rates of 2 3 mm yr by contrast where rifts do not overlap either a classic transform fault maskali or an oblique transfer zone mak arrasou kinematically connects them the length of the aden red sea overlap has increased in the last 0 9 myr as the aden plate boundary propagated northward into afar as a consequence the first order blocks that we identify within the overlap did not all rotate during the same time span nor by the same amounts similarly the major faults that bound them did not necessarily initiate and grow as their neighboring faults did despite these variations in strain distribution and kinematics the overlap kept accommodating a constant amount of strain 7 to 15 of the extension amount imposed by plate driving forces which remained distributed on a limited number seven or eight of major faults each one having slipped at constant rates 3 and 2 mm yr for vertical and lateral rates respectively the fault propagation rates and the block rotation rates that we either measure or deduce are so fast 30 130 mm yr and 15 38 myr respectively that they imply that strain transfer processes are transient as has been shown to be the case for the processes of tearing rift propagation and strain jumps in afar

local earthquake magnitude scale and seismicity rate for the ethiopian rift a calibrated local earthquake magnitude scale is essential for quantitative analyses of seismicity in ethiopia effective monitoring of earthquakes and resulting assessment of seismic hazard are especially important as regions with seismic and volcanic activity coincide with regions of economic significance and population growth we have developed a local magnitude m l scale for the northern main ethiopian rift mer using earthquake data collected during 2001 2003 on 122 three component broadband seismic stations waveform data from 2139 local earthquakes were corrected for instrument response and convolved with the nominal wood anderson torsion seismograph response appropriate for the original definition of local magnitude the hypocentral distances considered are 5 to 800 km with the best represented range from 5 to 150 km a total of 30 908 maximum zero to peak amplitudes awa were incorporated into a direct linear inversion for individual earthquake local magnitudes m l 244 station factors c and 2 linear distance dependent factors n k in the distance correction term log a o of the equation for local magnitude m l log awa log a o c the resulting distance correction is given by log a o 1 196997log r 17 0 001066 r 17 2 which implies that ground motion attenuation is relatively high consistent with ongoing magma intrusion and the presence of shallow magma reservoirs beneath the mer station corrections significantly reduce m l residuals and range between 0 42 m l units the catalog of earthquakes is complete above m l 2 1 and the annual cumulative seismicity rate follows the relation log n 4 5 1 13 m l our results are critical for accurate routine quantitative analysis of past current and future seismicity in ethiopia

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

september 2005 manda hararo dabbahu rifting event afar ethiopia constraints provided by geodetic data faults the volume of the 2005 dike 1 5 2 0 km 3 is not balanced by sufficient volume loss at dabbahu and gabho volcanoes 0 42 and 0 12 km 3 respectively taking into account the deflation of a suspected deep midsegment magma chamber simultaneously to dike intrusion produces a smoother opening distribution along the southern segment above the dike faults slipped by an average 3 m yielding an estimated geodetic moment of 3 5 10 19 nm one order of magnitude larger than the cumulative seismic moment released during the earthquake swarm between dabbahu and ado ale volcanic complexes significant opening occurred on the western side of the dike the anomalous location of the dike at this latitude offset to the east of the axial depression may explain this phenomenon a two stage intrusion scenario is proposed whereby rifting in the northern manda hararo rift was triggered by magma upwelling in the dabbahu area at the northern extremity of the magmatic segment although vigorous dike injection occurred during the september 2005 event the tectonic stress deficit since the previous rifting episode was not fully released leading to further intrusions in 2006 2009

numerical modelling of quaternary deformation and post rifting displacement in the asal ghoubbet rift djibouti africa over the last three decades a host of information on rifting process relating to the geological and thermal structure long time scale deformation quaternary and holocene and rifting cycle displacement across the asal ghoubbet rift has been made available these data are interpreted with a two dimensional thermo mechanical model that incorporates rheological layering of the lithosphere dyke inflation and faulting active fault locations and geometry are mainly controlled by both thermal structure and magma intrusion into the crust the distributed slip throughout the inner rift is related to the closeness of magma chamber leading to additional stress into the upper thinned crust assuming a constant arabia somalia motion of 11 mm year the variation of subsidence rate between the last 100 and 9 ka is associated with a decrease of the average injection rate from 10 to 5 mm year these values about equal to the regional opening rate suggest that both volcanism and tectonic play an equivalent role in the rifting process our modelled sequence of events gives one possible explanation for both vertical and horizontal displacements observed since the 1978 seismovolcanic crisis although part of the post rifting deformation could be due to viscous relaxation the high opening rate in the first years after the event and the abrupt velocity change in 1984 1986 argue for a large dyke inflation of 12 cm year ending in 1985 the asymmetric and constant pattern of the gps velocity since 1991 suggests that present post rifting deformation is mainly controlled by fault creep and regional stretching this study demonstrates the internal consistency of the data set highlights the role of magmatism in the mechanics of crustal stretching and reveals a complex post rifting process including magma injection fault creep and regional stretching d 2005 elsevier b v all rights reserved

recent rift related volcanism in afar ethiopia abstract rift zones are the most common magmatic environment on earth however opportunities to observe active rifting are rare and consequently the volcanological characteristics of rift systems are not well understood an ongoing phase of magmatic rifting along a section of the red sea system in afar ethiopia presents an exceptional opportunity to constrain relationships between volcanism and crustal growth here by integrating analyses of satellite images i e modis omi aster and ali with field observations we characterise two recent august 2007 and june 2009 basaltic fissure eruptions in afar and evaluate the role and significance of volcanism in the rifting process both events were brief 36 72 h and erupted 4 4 18 106 m3 of lava from a fissure system 4 6 5 km in length data from the spaceborne ozone monitoring instrument omi suggests total so2 emissions for each eruption of 26 5 kt 2007 and 34 7 kt 2009 consistent with complete degassing of the erupted magma volumes using geodetic models for the intrusive activity in afar we estimate the partitioning of magma between intrusive and extrusive components up to july 2009 to be 180 1 comparing the first order volcanic characteristics and the intrusive extrusive volume balance for the afar volcanism with data from the 1975 1984 krafla rifting cycle iceland suggests that the volcanic flux in afar will rise significantly over the next few years as the stresses are increasingly relieved by dyking and subsequent dykes are able to propagate more easily to the surface as a consequence basaltic fissure eruptions in this section of the afar rift will become of increasing large magnitude as the rifting event matures over the next 5 10 yr using available models of magmatic rifting we forecast the likely size and location of future eruptions in afar

dike induced faulting in rift zones of iceland and afar geodetic data and field observations demonstrate that the emplacement of dikes in volcanic rift zones frequently generates normal faulting and graben subsidence at the earth9s surface elastic modeling of the vertical ground surface displacements above dikes and faults indicates that the extent of graben subsidence can be achieved only if fault slip extends virtually to or beyond the dike plane at depth a mechanical model that includes dikes and frictional faults shows that dike opening tends to compress and lock faults located to either side of the dike therefore slip extending into or beyond the dike cavity must occur either 1 on faults that intersect the dike near its top above the zone of dike induced compression or 2 on faults that slip ahead of the dike as it propagates laterally data from iceland indicate that slip occurred on deep faults that presumably slipped in advance of the laterally propagating dike

evolution of a lateral dike intrusion revealed by relatively relocated dike induced earthquakes the 2014 15 barðarbunga holuhraun rifting event iceland abstract understanding dikes is vital as they serve both as bodies that build the crust and as conduits that feed eruptions and must be monitored to evaluate volcanic hazard during the 2014 15 bar arbunga rifting event iceland intense seismicity accompanied the intrusion of a 50 km lateral dike which culminated in a 6 month long eruption we here present relocations of earthquakes induced by the lateral dike intrusion using cross correlated sub sample relative travel times the 100 m spatial resolution achieved reveals the complexity of the dike propagation pathway and dynamics jerky segmented and allows us to address the precise relationship between the dike and seismicity with direct implications for hazard monitoring the spatio temporal characteristics of the induced seismicity can be directly linked in the first instance to propagation of the tip and opening of the dike and following this after dike opening indicate a relationship with magma pressure changes i e dike inflation deflation followed by a general post opening decay seismicity occurs only at the base of the dike where dike imposed stresses combined with the background tectonic stress from regional extension over 200 yr since last rifting are sufficient to induce failure of pre existing weaknesses in the crust while the greatest opening is at shallower depths emplacement oblique to the spreading ridge resulted in left lateral shear motion along the distal dike section studied here and a prevalence of left lateral shear failure fault plane strikes are predominately independent of the orientation of lineations delineated by the hypocenters indicating that they are controlled by the underlying host rock fabric this high resolution study provides unprecedented opportunity for comparison with both geodetic and field frozen dike observations and development and consolidation of analytical and analogue models with implications for rifting processes and real time monitoring of magma intrusion

length and timescales of rift faulting and magma intrusion the afar rifting cycle from 2005 to present although fault and magmatic processes have achieved plate spreading at mid ocean ridges throughout earth s history discrete rifting episodes have rarely been observed this paper synthesizes ongoing seismic structural space based geodetic and petrologic studies from the subaerial red sea rift in ethiopia where a major rifting episode commenced in september 2005 our aims are to determine the length and timescales of magmatism and faulting the partitioning of strain between faulting and magmatism and their implications for the maintenance of along axis segmentation most of the magma for the initial and subsequent 12 intrusions was sourced from the center of the dabbahu manda hararo rift segment strain is accommodated primarily by axial dike intrusions fed from mid segment magma chamber s these findings show that episodic approximate century interval rapid opening of discrete rift segments is the primary mechanism of plate boundary deformation the scale 65 km 8 km and intensity of crustal

january 2006 seafloor spreading event at 9 50 n east pacific rise ridge dike intrusion and transform fault interactions from regional hydroacoustic data 1 an array of autonomous underwater hydrophones is used to investigate regional seismicity associated with the 22 january 2006 seafloor spreading event on the northern east pacific rise near 9 50 n significant earthquake activity was observed beginning 3 weeks prior to the eruption where a total of 255 earthquakes were detected within the vicinity of the 9 50 n area this was followed by a series of 252 events on 22 january and a rapid decline to background seismicity levels during the subsequent 3 days because of their small magnitudes accurate locations could be derived for only 20 of these events 18 of which occurred during a 1 h period on 22 january these earthquakes cluster near 9 45 n and 9 55 n at the distal ends of the young lava flows identified posteruption where the activity displays a distinct spatial temporal pattern alternating from the north to the south and then back to the north this implies either rapid bilateral propagation along the rift or the near simultaneous injection of melt vertically from the axial magma lens short duration t wave risetimes are consistent with the eruption of lavas in the vicinity of 9 50 n on 22 january 2006 eruptions on 12 and 15 16 january also may be inferred from the risetime data however the locations of these smaller magnitude events cannot be determined accurately roughly 15 h after the last earthquakes were located adjacent to the eruption site a sequence of 16 earthquakes began to the north northeast at a distance of 25 40 km from the 9 50 n site these events are located in vicinity of the clipperton transform and its western inside corner an area from which the regional hydrophone network routinely detects seismicity coulomb stress modeling indicates that a dike intrusion spanning the known eruptive zone to the south 9 46 9 56 n would act to promote normal faulting or a combination of normal faulting and transform slip within this region with stress changes on the order of 1 10 kpa

surface deformation and strike slip faulting controlled by dyking and host rock lithology a compendium from the krafla rift iceland abstract in the present work we have obtained a detailed complete view of surface deformation at a volcano tectonic structure belonging to the krafla fissure swarm ne iceland through the integration of field mapping and the aerial structure from motion technique asfm as well as through orthomosaics and digital surface models dsms the documented 2 5 km long deformation zone is marked by a range of features such as topographic bulging parallel extension fractures and narrow grabens affected by local floor uplift which we have explained as the effect of the northward shallow propagation of a dyke from the krafla magma chamber our lithostratigraphic field survey of the area at a 1 5000 scale indicates that major changes in the pattern of surface deformation take place across the contacts between deposits with different rheological properties the transition from very stiff lavas to softer hyaloclastites produces changes from extensional fracturing to normal faulting though the width of the deformation zone does not change we have also documented a number of ne sw left lateral and nnw sse right lateral strike slip fault zones that are rotated in clockwise and anticlockwise sense to the main nne ssw graben trend and extend outwards as much as about 17 m the graben zone terminates to the north by splitting into two strike slip fault zones these were likely generated at the front of the dyke tip during its propagation and were later on bypassed by the advancement of the dyke as the dyke entered the soft hyaloclastites propagation was slowed down and the dyke became thicker

widening of the krafla fissure swarm during the 1975 1981 volcano tectonic episode frequent distance measurements across the krafla fissure swarm north iceland recorded the extension accompanying the sequence of rifting events which started in december 1975 and lasted for 6 years an 80 to 90 km long section of the fissure swarm extended during this sequence of rifting and volcanic events maximum widening of about 8 m occurred 10 to 12 km north of leirhnjukur which is located above the center of the krafla magma reservoir from that location the amount of widening decreased north wards and is estimated to exceed 2 m where the seismicity indicated the northern termination of the present rifting off the north coast about 70 km north of leirhnjukur the amount of widening also decreased southwards and approached zero at 15 to 20 km south of leirhnjukur the ground deformation associated with these rifting events can be summarised as a narrow strip 1 to 2 km wide along the fissure swarm is heavily fractured with numerous open cracks parallel to the fissure swarm this fractured strip has subsided 2 to 3 m relative to its flanks the flanks of the fractured zone have been uplifted relative to regions farther away the uplift is not well constrained but tilt observations at several locations indicate about 1 m uplift the flanks of the rift zone have contracted perpendicular to the fissure swarm the maximum contractional strain exceeds 300 mm per km the amount of areal expansion windening of the fissure swarm times the length of the fractured zone associated with these rifting events is estimated to be about 0 30 km2 for individual events the area of expansion has been roungly proportional to the volume of subsidence above the krafla magma reservoir if the width of a new dike is equal widening of surface fissures the ratio of the subsidence volume to the area of expansion for the best observed events indicated a height of a new dike system as 2 4 to 2 8 km this ratio is significantly less for events of large lava production but even during these events the majority of magma leaving the krafla reservoir was apparently emplaced in subsurface fissures

the afar depression transition between continental rifting and sea floor spreading abstract detailed geophysical and geological studies were undertaken in the 1970 s in the afar depression ethiopia in order to study the processes of continental break up occurring in this unique setting the afar is located between the ethiopian and somali plateaus and the red sea we have re evaluated the available seismic refraction and gravity data using modern interpretation techniques and incorporated new information derived from magnetic gravity seismicity and seismic refraction data the results of the re evaluation indicate that the afar depression is a transition zone between the continental rifts of kenya and the present ongoing sea floor spreading of the red sea and the gulf of aden in the depression itself the crust is greatly attenuated and is underlain by a low velocity 7 4 7 5 km s high temperature upper mantle material however it is still a continental crust which thins to the northeast and east towards the red sea and the gulf of aden it is not a newly formed oceanic crust as has been suggested in the past the trends of the areas of crustal thinning are indicated by the gravity anomalies the trend of the ethiopian rift reaches south afar where the trend of crustal thinning changes and becomes parallel to the red sea trend it is offset en echelon from the spreading area of the red sea which is marked by an alignment of gravity maxima the gulf of aden trend continues west towards the ethiopian rift the gravity trends coincide with areas of maximum crustal thinning and may thus suggest a possible triple junction in development

integrated brachiopod based bioevents and sequence stratigraphic framework for a late ordovician subpolar platform eastern anti atlas morocco the upper ordovician katian hirnantian brachiopods of tafilalt eastern anti atlas are locally abundant diverse and well preserved providing a near continuous record of faunal change on a high latitude siliciclastic dominated platform a chronostratigraphic framework based on brachiopod distribution and preservation in shell accumulation events and integrated with sequence stratigraphy has been generated for the katian interval which has allowed correlation with the chitinozoan based chronostratigraphic and sequence stratigraphic framework erected for the central anti atlas in tafilalt two katian transgressive regressive composite depositional sequences c 60 and 170 m thick and related to third order fluctuations in sea level were unaffected by hirnantian glaciogenic erosion they were deposited on a mixed platform with a bryonoderm association dominated by brachiopods bryozoans and echinoderms brachiopods developed in high energy inner shelf areas whereas bryozoans mainly trepostomates and fenestrates and pelmatozoans cystoids and crinoids dominated in low energy outer shelf areas brachiopod accumulations mark distinct event surfaces such as lag and event concentrations hydraulic simple and composite concentrations related to transgressive surfaces and hiatal condensed concentrations marking maximum flooding surfaces the taphonomic condensation displayed by the hirnantian alnif member which onlaps the erosive base of glaciogenic tunnel channels is explained as reworking and resedimentation of allochthonous robust biogenic hard parts sourced from the underlying katian ktaoua group

acoustic detection of a seafloor spreading episode on the juan de fuca ridge using military hydrophone arrays until recently no practical method has been available to continuously monitor seismicity of seafloor spreading centers the availability of the u s navy s sound surveillance system sosus for environmental research has allowed the continuous monitoring of low level seismicity of the juan de fuca ridge in the northeast pacific on june 22 1993 noaa installed a prototype system at u s naval facility whidbey island to allow real time acoustic monitoring of the juan de fuca ridge using sosus on june 26 2145 gmt a burst of low level seismic activity with accompanying harmonic tremor was observed and subsequently located near 46 15 n 129 53 w on the spreading axis of the juan de fuca ridge over the following 2 days the activity migrated to the nne along the spreading axis with the final locus of activity near 46 31 5 n 129 35 w the nature of the seismicity was interpreted to represent a lateral dike injection with the possibility of eruption on the seafloor based on this interpretation a response effort was initiated by u s and canadian research vessels and both warm water plumes and fresh lavas were subsequently identified at the site

propagation of rifting along the arabia somalia plate boundary the gulfs of aden and tadjoura the localization and propagation of rifting between arabia and somalia are investigated by assessing the deformation geometry and kinematics at different scales between the eastern gulf of aden and the gulf of tadjoura using bathymetric magnetic seismological and structural evidence large scale southwestward propagation of the aden ridge markedly oblique to the arabia somalia relative motion vector began about 30 myr ago between the error and sharbithat ridges it was an episodic process with stages of rapid propagation mostly at rates 10 cm yr interrupted by million year pauses on transverse discontinuities coinciding with rheological boundaries between different crustal provinces of the arabia somalia plate the longest pause was at the shukra el sheik discontinuity 45 e where the ridge tip stalled for 13 myr between 17 and 4 ma west of that discontinuity rifting and spreading took place at an azimuth n25 10 e and rate 1 2 0 3 cm yr different from those of the global arabia somalia motion vector n39 1 73 cm yr implying an additional component of movement n65 10 e 0 7 0 2 cm yr due to rotation of the danakil microplate at shukra el sheik the typical oceanic ridge gives way to a narrow wsw trending axial trough resembling a large fissure across a shallow shelf this trough is composed of about eight rift segments which result from normal faulting and fissuring along n110 n130 e trends all the segments step to the left southwestward mostly through oblique transfer zones with en echelon normal faults only two segments show clear significant overlap there is one clear transform the maskali fault between the obock and tadjoura segments the latter segment which encroaches onland is composed of two parallel subrifts iboli ambabbo that propagated northwestward and formed in succession the most recent southwestern subrift ambabbo represents the current tip of the aden ridge we propose a mechanical model in which the large scale propagation of the ridge followed a wsw trending zone of maximum tensile stress while the small scale propagation of its nw trending segments was dictated by the orientation of that stress oblique propagation was a consequence of passive lithospheric necking of the arabia somalia plate along its narrow section in map view between socotra and the kink of the red sea ethiopian rift above the afar plume individual ridge segments oriented roughly perpendicular to plate motion like lithospheric cracks were forced to jump southward because of confinement within the necking zone self sustaining plate scale necking may explain why the aden ridge did not connect with the red sea through bab el mandeb but continued straight into afar

evidence from the ad 2000 izu islands earthquake swarm that stressing rate governs seismicity magma intrusions and eruptions commonly produce abrupt changes in seismicity far from magma conduits that cannot be associated with the diffusion of pore fluids or heat such swarm seismicity also migrates with time and often exhibits a dog bone shaped distribution the largest earthquakes in swarms produce aftershocks that obey an omori type exponential temporal decay but the duration of the aftershock sequences is drastically reduced relative to normal earthquake activity here we use one of the most energetic swarms ever recorded to study the dependence of these properties on the stress imparted by a magma intrusion a 1 000 fold increase in seismicity rate and a 1 000 fold decrease in aftershock duration occurred during the two month long dyke intrusion we find that the seismicity rate is proportional to the calculated stressing rate and that the duration of aftershock sequences is inversely proportional to the stressing rate this behaviour is in accord with a laboratory based rate state constitutive law suggesting an explanation for the occurrence of earthquake swarms any sustained increase in stressing rate whether due to an intrusion extrusion or creep event should produce such seismological behaviour

comparison of dike intrusions in an incipient seafloor spreading segment in afar ethiopia seismicity perspectives 1 oceanic crust is accreted through the emplacement of dikes at spreading ridges but the role of dike intrusion in plate boundary deformation during continental rupture remains poorly understood between 2005 and 2009 the 70 km long dabbahu manda hararo rift segment in ethiopia has experienced 14 large volume dike intrusions 9 of which were recorded on temporary seismic arrays a detailed comparison of the seismic characteristics of the seismically monitored dikes is presented with implications for dike intrusion processes and magmatic plumbing systems all of the migrating swarms of earthquakes started from a 80 of energy is released during the propagation phase with minimal seismic energy release after the dike propagation ceases we interpret that faulting and graben formation above the dikes occurs hours after the passage of the dike tip coincident with the onset of low frequency earthquakes dike lengths show no systematic reduction in length with time suggesting that topographic loading and stress barriers influence dike length as well as changes in tectonic stress the propagation velocities of all the dikes follow a decaying exponential northward propagating dikes had faster average velocities than those that propagated southward suggesting preconditioning by the 2005 megadike or ongoing heating from a subcrustal magma source north of the midsegment

uppermost mantle pn velocity model for the afar region ethiopia an insight into rifting processes the afar depression ethiopia offers unique opportunities to study the transition from continental rifting to oceanic spreading because the process is occurring onland using traveltime tomography and data from a temporary seismic deployment we describe the first regional study of uppermost mantle p wave velocities vpn we find two separate low vpn zones as low as 7 2 km s 1 beneath regions of localized thinned crust in northern afar indicating the existence of high temperatures and potentially partial melt the zones are beneath and off axis from contemporary crustal magma intrusions in active magmatic segments the dabbahu manda hararo and erta ale segments this suggests that these intrusions can be fed by off axis delivery of melt in the uppermost mantle and that discrete areas of mantle upwelling and partial melting thought to characterize segmentation of the uppermost mantle at seafloor spreading centres are initiated during the final stages of break up

refining fault slip rates using multiple displaced terrace risers an example from the honey lake fault ne california usa abstract faulted terrace risers are semi planar features commonly used to constrain quaternary slip rates along strike slip faults these landforms are difficult to date directly and therefore their ages are commonly bracketed by age estimates of the adjacent upper and lower terrace surfaces however substantial differences in the ages of the upper and lower terrace surfaces a factor of 2 4 difference observed globally produce large uncertainties in the slip rate estimate in this investigation we explore how the full range of displacements and bounding ages from multiple faulted terrace risers can be combined to yield a more accurate fault slip rate we use 0 25 m cell size digital terrain models derived from airborne lidar data to analyze three sites where terrace risers are offset right laterally by the honey lake fault in ne california usa we use ages for locally extensive subhorizontal surfaces to bracket the time of riser formation an upper surface is the bed of abandoned lake lahontan having an age of 15 8 0 6 ka and a lower surface is a fluvial terrace abandoned at 4 7 0 1 ka we estimate lateral offsets of the risers ranging between 6 6 and 28 3 m median values a greater than fourfold difference in values the amount of offset corresponds to the riser s position relative to modern stream meanders the smallest offset is in a meander cutbank position whereas the larger offsets are in straight channel or meander point bar positions taken in isolation the individual terrace riser offsets yield slip rates ranging from 0 3 to 7 1 mm a however when the offset values are collectively assessed in a probabilistic framework we find that a uniform linear slip rate of 1 6 mm a 1 4 1 9 mm a at 95 confidence can satisfy the data within their respective uncertainties this investigation demonstrates that integrating observations of multiple offset elements crest midpoint and base from numerous faulted and dated terrace risers at closely spaced sites can refine slip rate estimates on strike slip faults

faroe iceland ridge experiment 2 crustal structure of the krafla central volcano the seismic velocity structure of the krafla central volcano is characterized by large variations in compressional velocity a 40 km wide high velocity dome extends from the lower crust 11 14 km depth beneath the volcano narrowing upward a magma chamber sits at its top near 3 km depth it is defined by both 0 2 0 3 s compressional wave delays and shear wave shadowing to be 2 3 km n s 8 10 km e w and 0 7 1 8 km thick the near surface structure uppermost 2 5 km of the krafla caldera is approximately flat lying with only minor lateral heterogeneities the crust beneath the magma chamber has low shear wave attenuation and anomalously high compressional and shear wave velocities shear waves reflected from a 19 km deep moho are clearly visible for some paths through the crustal volume below the magma chamber even though the more shallow diving s waves are severely attenuated the midcrust beneath the shallow magma chamber cannot contain partial melt or even be at near solidus temperatures the krafla central volcano plays a major role in crustal genesis along the plate boundary the high velocity dome in our view represents crust generated in and around the magma chamber which has subsequently been advected to greater depths

a constant influx model for dike propagation implications for magma reservoir dynamics 1 most observations of seismicity rate during dike propagation on basaltic volcanoes show 1 rate stationarity despite possible variations of the dike tip velocity 2 frequent lack of clear and monotonic hypocenter migration following dike propagation and 3 event occurrences located backward with respect to the dike tip position on these bases the origin of the seismicity contemporary to dike intrusion within basaltic volcanoes cannot be solely related to the crack tip propagation seismicity rather appears to be the response of the edifice itself to the volumetric deformation induced by the magma intruding the solid matrix the volume change induced into the volcano edifice over time by the intruding magma is equal to the magma flux injected into the dike from the reservoir the consequence of this is that the stationary seismicity rate observed during the intrusion is a proxy for the magma flux withdrawn from the reservoir we consider a two phase dike propagation model including a first vertical propagation followed by a lateral migration along a lithological discontinuity we explore 1 under which geophysical conditions the vertical dike is fed at constant flow rate of magma and 2 dike propagation patterns implications entailed by constant volumetric flux on the piton de la fournaise volcano case study suggest a minimum size for the magma reservoir of about 1 km3 and a maximum value for the initial magma reservoir overpressure of about 2 2 mpa considering similar magma inflow rates during vertical and lateral dike propagation phases we reproduce independent estimates of propagation velocities rise times and injected volumes when applying the model to the august 2003 piton de la fournaise eruption

characteristics of tectonomagmatic earthquake swarms at the southwest indian ridge between 16 e and 25 e summary the ultraslow spreading southwest indian ridge swir is a prominent end member of the global mid ocean ridge system it spreads with a full rate of 14 16 mm y 1 and shows several segments of various obliquities the western swir consists of the oblique and orthogonal supersegments lying at an epicentral distance of 21 to the vna2 seismic array operated by the german neumayer station in east antarctica the array monitors backazimuth apparent velocity and signal to noise ratio of arriving waves and provides a data set of seismicity from the western swir over several years compared to the global seismological network its detection threshold for earthquakes occurring at the western swir is more than 0 5 mb lower enabling a more comprehensive study of mid ocean ridge processes than the teleseismic earthquake catalogues we identified a total number of 743 earthquakes occurring at the western part of the swir and calculated the body wave magnitudes mb from p wave amplitude picks on the vna2 broad band sensor obtaining a magnitude range from mb 3 18 to mb 5 34 in the years of 2001 2004 2005 and 2008 significantly increased event rates indicated four earthquake swarms with up to 164 events lasting for several days all swarms had strong events registered in the international seismological centre catalogue the relocalization of these events confirmed that all swarms occurred in the same region on the orthogonal supersegment we analysed event and moment release rate histories b values and aftershock decay rates modified omori law finding that the swarms of 2001 2004 and 2005 have similarities in the temporal distribution of seismic moment and event numbers the swarm of 2008 is smaller with high magnitude events at the swarm s onset which represent shear failure on normal faults the application of the modified omori law and the b value show that the earthquakes of the swarms do not follow the classical main shock aftershock pattern of purely tectonic earthquake sequences at the orthogonal supersegment a continuous positive magnetic anomaly along the rift axis a negative mantle bouguer anomaly basalts at the seafloor and potentially volcanic edifices indicate robust magmatic crustal accretion the high resolution bathymetry shows ubiquitous rift parallel ridges with steep flanks towards the rift axis indicating high angle normal faults the high magnitude earthquakes detected teleseismically during the swarms are generated there we interpret that the swarms are caused by magmatic accretion episodes at a suggested volcanic centre of the orthogonal supersegment and that possible magma injection activates the steeply dipping fault planes

the afar depression interpretation of the 1960 2000 earthquakes summary we studied the seismic activity of the afar depression ad and adjacent regions during the period 1960 2000 we define seven distinct seismogenic regions using geological tectonic and seismological data based on the frequency magnitude relationships we obtain b values of about 1 for the different regions the pattern of the distribution of the location of epicentres fits with the known active fault zone in the ad and the axial volcanic ridges the bab el mandab area and the danakil aysha a blocks are less active for 125 intermediate to strong earthquakes the seismic moment and source parameters were calculated the results of the fault plane solutions for the afar depression indicate mainly strike slip and normal sense of movement originating from fault planes striking nw se these results indicate a clockwise block rotation described previously as a bookshelf model in central ad there are a few rightlateral faults east of massawa with e w striking fault planes at the southern red sea north of the danakil block the mixed focal mechanisms with axial plane striking nw se comprise several reverse faulting strike slip motion and normal faulting right lateral movement was also calculated for a cluster of seismic events between the manda hararo and alyata volcanic ridges along nw se striking faults along the n s striking faults in the escarpment at the western afar margins there are two distinct clusters of epicentres the strong earthquakes at the southern cluster exhibit normal or strike slip motions the intermediate to small earthquakes in the northern cluster exhibit reverse and strike slip motions mainly normal faults were calculated along ne sw striking faults of the ethiopian east african rift estimates of the seismic efficiency suggest that the maximal values are about 50 per cent or less implying that most of the motion is taken aseismically

shallow axial magma chamber at the slow spreading erta ale ridge thermal models predict that spreading velocity at mid ocean rifts should influence the geometry of the underlying magma chamber insar data from the slow spreading ethiopian rift identify a shallow elongated magma chamber a feature usually associated with fast spreading rifts beneath the erta ale segment implying that spreading velocity may not be so important after all

temporal and spatial history of the 1999 2000 endeavour segment seismic series juan de fuca ridge 1 two large clusters of earthquake activity in june of 1999 and january of 2000 have dominated recent seismicity along the endeavour segment of the juan de fuca ridge the impacts of the june 1999 sequence on the hydrologic system which include changes in vent temperature and chemistry within the main endeavour vent field have been well documented previously analysis of seismic and hydroacoustic data indicates that both sequences exhibit a swarm like behavior characterized by the absence of a dominant main shock event the epicentral locations of events within the two swarms overlap spatially with centroid positions near 47 49 and 47 46 n latitude during the june 1999 swarm the initial activity spans the along axis region where a shallow axial magma chamber reflector was later imaged the epicenters then migrate 12 km to the south at a rate of 0 3 m s consistent with lateral dike propagation a distinct subcluster of events also occurred in the vicinity of surveyor volcano on the overlapping portion of the cobb segment given its distance from the main swarm this activity may represent a triggered response to dynamic shaking the january 2000 swarm has a more limited along axis extent relative to the june 1999 swarm with no indication of lateral migration much of this activity is concentrated in a region predicted to have undergone extension due to dike propagation in 1999 although it contains fewer total events and is of shorter duration relative to the june 1999 swarm the january 2000 activity exhibits a higher peak rate of seismicity and greater mean event magnitude as in situ temperature monitoring was not in place during january 2000 and vent fluids were not sampled until june 2000 the impacts of this swarm on the hydrothermal system are unknown the southernmost tip of the endeavour segment also is found to be a region of repeating swarm activity although morphologic evidence indicates the cobb segment has been propagating northward recently this seismic activity suggests that the western limb of the endeavour cobb overlap zone remains active

structure and evolution of the eastern gulf of aden conjugate margins from seismic reflection data the gulf of aden is a young and narrow oceanic basin formed in oligo miocene time between the rifted margins of the arabian and somalian plates its mean orientation n75 e strikes obliquely 50 to the n25 e opening direction the western conjugate margins are masked by oligo miocene lavas from the afar plume this paper concerns the eastern margins where the 19 35 ma breakup structures are well exposed onshore and within the sediment starved marine shelf those passive margins about 200 km distant are non volcanic offshore during the encens sheba cruise we gathered swath bathymetry single channel seismic reflection gravity and magnetism data in order to compare the structure of the two conjugate margins and to reconstruct the evolution of the thinned continental crust from rifting to the onset of oceanic spreading between the alula fartak and socotra major fracture zones two accommodation zones trending n25 e separate the margins into three n110 e trending segments the margins are asymmetric offshore the northern margin is narrower and steeper than the southern one including the onshore domain the southern rifted margin is about twice the breadth of the northern one we relate this asymmetry to inherited jurassic cretaceous rifts the rifting obliquity also influenced the syn rift structural pattern responsible for the normal faults trending from n70 e to n110 e the n110 e fault pattern could be explained by the decrease of the influence of rift obliquity towards the central rift and or by structural inheritance the transition between the thinned continental crust and the oceanic crust is characterized by a 40 km wide zone our data suggest that its basement is made up of thinned continental crust along the southern margin and of thinned continental crust or exhumed mantle more or less intruded by magmatic rocks along the northern margin abstract from author

three dimensional seismic structure of the mid atlantic ridge 35 n evidence for focused melt supply and lower crustal dike injection 1 we gathered seismic refraction and wide angle reflection data from several active source experiments that occurred along the mid atlantic ridge near 35 n and constructed three dimensional anisotropic tomographic images of the crust and upper mantle velocity structure and crustal thickness the tomographic images reveal anomalously thick crust 8 9 km and a low velocity bull s eye from 4 to 10 km depth beneath the center of the ridge segment the velocity anomaly is indicative of high temperatures and a small amount of melt up to 5 and likely represents the current magma plumbing system for melts ascending from the mantle in addition at the segment center seismic anisotropy in the lower crust indicates that the crust is composed of partially molten dikes that are surrounded by regions of hot rock with little or no melt fraction our results indicate that mantle melts are focused at mantle depths to the segment center and that melt is delivered to the crust via dikes in the lower crust our results also indicate that the segment ends are colder receive a reduced magma supply and undergo significantly greater tectonic stretching than the segment center

spatial and temporal patterns of deformation at the tendaho geothermal prospect ethiopia abstract observations of ground deformation in east africa have been fundamental for unveiling the tectonics of continental rifting assessing the seismic and volcanic hazard to development and identifying geothermal resources here we investigate the active natural and anthropogenic processes in the tendaho graben afar using interferometric synthetic aperture radar insar collected by the envisat satellite in 2004 2010 we used the poly interferometric rate and time series estimation rate method to calculate displacement in satellite line of sight and a least square inversion to decompose the line of sight displacement into vertical and rift perpendicular components we observe two zones of deformation a 20 km wide circular region of subsidence located 10 km northeast of the town of semera with a maximum displacement rate of 5 cm yr and elongated zone 50 km of subsidence in the area of the geothermal prospect maximum rate of 4 cm yr the temporal characteristics of subsidence varies between these zones with an increase in subsidence rate observed in the circular region in august 2008 we used a bayesian inversion to find the best fitting source models and compared this to locations of seismicity and other geophysical observations the pattern of deformation is consistent with a combination of magmatic and geothermal processes but there does not appear to be a direct link to a sequence of dyke intrusions during 2005 2010 at manda hararo graben 60 km away but dynamic stress changes or deep crustal flow could account for the observations

imaging composite dike propagation etna 2002 case 1 late on the night of 26 october 2002 a dike intrusion started suddenly at mount etna producing intense explosive activity and lava effusion on the southern flank five to six hours afterward a long field of eruptive fractures propagated radially along the northeastern flank of the volcano producing marked variations at the permanent tilt network the dike propagation velocity was inferred by the associated seismicity we modeled the temporal evolution of the continuously recorded tilt data both during the vertical dike propagation on the high south flank on 26 october and during the radial propagation along the northeast flank between 27 and 28 october the reproduction of the recorded tilt signal allowed us to describe the geometry and characteristics of the two dikes in greater detail than the previous static inversion we deduced that the eruption was characterized by an unusual composite mechanism clearly showing a transition from a nearly pure opening mode displacement to a mechanism characterized by an equally strong normal dip slip component and a smaller left lateral strike slip component in this study we demonstrate the interaction between the final segment of the dike and a preexisting structure that was reactivated in response to the intrusion we show that tilt and its modeling represent a powerful tool to verify and constrain dike intrusions in detail

precursory swarms of long period events at redoubt volcano 1989 1990 alaska their origin and use as a forecasting tool abstract during the eruption of redoubt volcano from december 1989 through april 1990 the alaska volcano observatory issued advance warnings of several tephra eruptions based on changes in seismic activity related to the occurrence of precursory swarms of long period lp seismic events dominant period of about 0 5 s the initial eruption on december 14 occurred after 23 years of quiescence and was heralded by a 23 hour swarm of lp events that ended abruptly with the eruption after a series of vent clearing explosions over the next few days dome growth began on december 21 another swarm with lp events similar to those of the first began on the 26th and ended in a major tephra eruption on january 2 eruptions continued over the next two weeks and then ceased until february 15 when a large eruption initiated a long phase of repetitive dome building and dome destroying episodes that continued into april warnings were issued before the major events on december 14 and january 2 but as the eruptive sequence continued after january 2 the energy of the swarms decreased and forecasting became more difficult a significant but less intense swarm preceded the february 15 eruption which was not forecast this eruption destroyed the only seismograph on the volcanic edifice and stymied forecasting until march 4 when the first of three new stations was installed within 3 km of the active vent from march 4 to the end of the sequence on april 21 there were eight eruptions six of which were preceded by detectable swarms of lp events although weak these swarms provided the basis for warnings issued before the eruptions on march 23 and april 6 the initial swarm on december 13 had the following features 1 short duration 23 hours 2 a rapidly accelerating rate of seismic energy release over the first 18 hours of the swarm followed by a decline of activity during the 5 hours preceding the eruption 3 a magnitude range from 0 4 to 1 6 4 nearly identical lp signatures with a dominant period near 0 5 s 5 dilatational first motions everywhere and 6 a stationary source location at a depth of 1 4 km beneath the crater this occurrence of long period events suggests a model involving the interaction of magma with groundwater in which magmatic gases steam and water drive a fixed conduit at a stationary point throughout the swarm the initiation of that sequence of events is analogous to the failure of a pressure relief valve connecting a lower supercharged magma dominated reservoir to a shallow hydrothermal system a three dimensional model of a vibrating fluid filled crack recently developed by chouet is found to be compatible with the seismic data and yields the following parameters for the lp source crack length 280 380 m crack width 140 190 m crack thickness 0 05 0 20 m crack stiffness 100 200 sound speed of fluid 0 8 1 3 km s compressional wave speed of rock 5 1 km s density ratio of fluid to rock 0 4 and ratio of bulk modulus of fluid to rigidity of rock 0 03 0 07 the fluid filled crack is excited intermittently by an impulsive pressure drop that varies in magnitude within the range of 0 4 to 40 bar such disturbance appears to be consistent with a triggering mechanism associated with choked flow conditions in the crack

correction to transient rift opening in response to multiple dike injections in the manda hararo rift afar ethiopia imaged by time dependent elastic inversion of interferometric synthetic aperture radar data 1 interferometric synthetic aperture radar insar data spanning the time intervals separating thirteen dike intrusions in the manda hararo dabbahu rift afar ethiopia from 2005 to 2009 show that transient deformation occurs in the inter diking period this deformation can be explained by the presence of seven inflating or deflating pressure sources by combining the data acquired on four different insar tracks through time dependent elastic models we are able to track these deformation modes with a time resolution smaller than 1 month sustained deflation of a deep magma reservoir at dabbahu in the 6 months following the main rifting event of 2005 and slow decelerating post eruptive re inflation of two shallow magma reservoirs below dabbahu and gabho volcanoes are monitored a deflation signal of deep origin on the neighboring rift system is also detected possibly caused by outflow of material from a preexisting reservoir into the deep plate boundary in contrast rapidly evolving deformation is observed at the center of the manda hararo rift segment transient deformation events are monitored in the weeks months following the diking events with pulses of localized rift opening after the dike intrusions followed by an exponential like decay of opening rate this signal may be associated with the replenishment of the central magma reservoir involved in feeding the 2005 2009 dikes alternatively the predominantly horizontal mode of deformation suggests an interaction between the response of the lithosphere to tectonic strain accumulation and the process of hydraulic connectivity within the central magma plumbing system

geodetic constraints on present day motion of the arabian plate implications for red sea and gulf of aden rifting 1 five years of continuously recording gps observations in the kingdom of saudi arabia together with new continuous and survey mode gps observations broadly distributed across the arabian peninsula provide the basis for substantially improved estimates of present day motion and internal deformation of the arabian plate we derive the following relative geodetic euler vectors latitude n longitude e rate myr counterclockwise for arabia nubia 31 7 0 2 24 6 0 3 0 37 0 01 arabia somalia 22 0 0 5 26 2 0 5 0 40 0 01 arabia india 18 0 3 8 87 6 3 3 0 07 0 01 arabia sinai 35 7 0 8 17 1 5 0 0 15 0 04 and arabia eurasia 27 5 0 1 17 6 0 3 0 404 0 004 we use these euler vectors to estimate present day stability of the arabian plate the rate and direction of extension across the red sea and gulf of aden and slip rates along the southern dead sea fault south of the lebanon restraining bend 4 5 4 7 0 2 mm yr left lateral 0 8 1 1 0 3 mm yr extension and the owens fracture zone 3 2 2 5 0 5 mm yr right lateral increasing from north to south 1 2 mm yr extension on a broad scale the arabian plate has no resolvable internal deformation weighted root mean square of residual motions for arabia equals 0 6 mm yr although there is marginally significant evidence for n s shortening in the palmyride mountains syria at 1 5 mm yr we show that present day arabia plate motion with respect to eurasia is consistent within uncertainties i e 10 with plate tectonic estimates since the early miocene when arabia separated from nubia we estimate the time of red sea and gulf of aden rifting from present day arabia motion plate tectonic evidence for a 70 increase in arabia nubia relative motion at 13 ma and the width of the red sea and gulf of aden and find that rifting initiated roughly simultaneously 2 2 myr along the strike of the red sea from the gulf of suez to the afar triple junction as well as along the west gulf of aden at 24 2 2 ma based on the present kinematics we hypothesize that the negative buoyancy of the subducted ocean lithosphere beneath the makran and the zagros fold thrust belt is the principle driver of arabia eurasia convergence and that resisting forces associated with arabia eurasia continental collision have had little impact on plate motion

strike slip faulting during the 2014 barðarbunga holuhraun dike intrusion central iceland over a 13 day period magma propagated laterally from the subglacial bar arbunga volcano in the northern rift zone iceland it created 30 000 earthquakes at 5 7 km depth along a 48 km path before erupting on 29 august 2014 the seismicity which tracked the dike propagation advanced in short bursts at 0 3 4 7 km h separated by pauses of up to 81 h during each surge forward seismicity behind the dike tip dropped moment tensor solutions from the leading edge show exclusively left lateral strike slip faulting subparallel to the advancing dike tip releasing accumulated strain deficit in the brittle layer of the rift zone behind the leading edge both left and right lateral strike slip earthquakes are observed the lack of non double couple earthquakes implies that the dike opening was aseismic

reply to comment on co2 variability in mid ocean ridge basalts from syn emplacement degassing constraints on eruption dynamics by soule et al earth planet sci lett 2012 327 328 39 49

on the small scale fractal geometrical structure of a living coral reef barrier topographical complexity of coral reefs is of primary importance for a number of hydrodynamical and ecological processes the present study is based on a series of high resolution seabottom elevation measurements along the maupiti barrier reef french polynesia several statistical metrics and spectral analysis are used to characterize the spatial evolution of the coral geometrical structure from the reef crest to the backreef a consistent fractal like power law exists in the spectral density of bottom elevation for length scales between 0 1 and 7m while at larger scale the reef structure shows a different pattern such a fine characterization of the reef geometrical structure provides key elements to reconstruct the reef history to improve the representation of reef roughness in hydrodynamical models and to monitor the evolution of coral reef systems in the global change context

effects of the crustal stress on the growth of dikes conditions of intrusion and extrusion of magma dike growth is analyzed on the basis of the fluid dynamical behavior of incompressible newtonian magma in a spatially uniform fissure of temporally variable width the length and width of the dike vary with the magma pressure at the dike entrance connected to the magma chamber and with the mean stress working normal to the dike walls the mean stress consists of the regional crustal stress and the mean elastic stress which is approximately evaluated using an analytic solution of the crack with a constant aperture the length and width of the dike are determined as a function of time by simple ordinary differential equations the dike can grow when the initial pressure of the entering magma is above the lithostatic pressure and high enough to suppress compression of the crustal stress the assumed size of the initial crack affects the growth only in the duration needed for an initial adjustment stage after this stage the width and length first grow exponentially and then the growth slows as magma is released from the chamber the growth in the later stage is classified into two types when the initial magma pressure is high enough the dike elongates unlimitedly while the width finally decreases this extrusive growth is likely to lead to eruptions the other intrusive growth in which the dike finally stops its growth with a finite length and width occurs when the crustal stress component normal to the fissure is of sufficiently strong extension buoyancy of magma is not very effective in the earlier growth stage but it shifts the marginal condition separating the extrusive and intrusive growths the extrusive and intrusive events that recently happened in izu oshima volcano and izu tobu volcano group japan are interpreted as the examples of these two growth types

sequence of rifting in afar manda hararo rift ethiopia 2005 2009 time space evolution and interactions between dikes from interferometric synthetic aperture radar and static stress change modeling thirteen dike intrusions in the manda hararo rift afar ethiopia from september 2005 to june 2009 studied using an extensive interferometric synthetic aperture radar insar data set provide insight into the mechanics of a major active rift kinematic inversions of insar data reveal that dikes opened by 0 8 3 5 m at an average 5 km depth with volumes of 0 04 0 2 km 3 with up to 12 m opening and a volume greater than 1 km 3 for the september 2005 megadike dikes have their source in a midsegment magma reservoir which induces a local shallowing of the brittle ductile boundary presumably due to thermal weakening of the lithosphere the smaller dikes in 2006 2009 were emplaced in regions of minimum opening of the september 2005 megadike above the central magma reservoir in contrast the most voluminous dike intrusions in 2006 2009 occurred near the locus of the peak of maximum opening of the september 2005 megadike 10 km north of the magma source this may suggest that tension on the plate boundary was highest there both prior to 2005 and possibly also after 2005 evolution and distribution of normal stress on the plate boundary throughout the rifting episode may indicate that tension near to the magma reservoir is lower than toward segment ends average relief of normal stresses of tectonic origin coeval to dike intrusions is comparable with shear stress drops for earthquakes presumably because dikes in the manda hararo rift are intruded at low magma pressure and high tectonic stress time space evolution and interactions between dikes from interferometric synthetic aperture radar and static stress change modeling

low frequency hybrid earthquakes near a magma chamber in afar quantifying path effects areas of active volcanism contain elaborate velocity structures that complicate interpretations of earthquake source mechanisms we examine the spectral characteristics of 805 earthquakes that immediately followed a large volume basaltic dike intrusion and associated silicic flank eruption of dabbahu volcano in the afar depression as recorded on near source seismometers we use these results to quantify the contribution of scattering and attenuation to the observed spectra of low frequency hybrid and volcano tectonic earthquake clusters from beneath dabbahu volcano and around the dike zone we find strong variations in the signal amplitude and frequency content of earthquakes recorded at stations separated by as little as 2 km caused by preferential attenuation of high frequencies depending on the vantage point these observations imply that there are large impedance contrasts near the cooling solidifying and recently intruded dike we estimate the intrinsic absorption attenuation coefficient qi and inverse scattering length g0 averaged over a 300 sq km area beneath dabbahu our results are consistent with the highest attenuation coefficients from studies of volcanic provinces in italy qi 1 0 02 g0 0 1 km 1 for a signal at 2 hz the magnitude of these two parameters indicates there are large impedance contrasts present in the area due to the recent intrusion of magma and associated fracturing

magmatic cycles pace tectonic and morphological expression of rifting afar depression ethiopia the existence of narrow axial volcanic zones of mid oceanic ridges testifies of the underlying concentration of both melt distribution and tectonic strain as a result of repeated diking and faulting axial volcanic zones therefore represent a spectacular topographic expression of plate divergence however the submarine location of oceanic ridges makes it difficult to constrain the interplay between tectonic and magmatic processes in time and space in this study we use the dabbahu manda hararo dmh magmatic rift segment afar ethiopia to provide quantitative constraints on the response of tectonic processes to variations in magma supply at divergent plate boundaries the dmh magmatic rift segment is considered an analogue of an oceanic ridge exhibiting a fault pattern extension rate and topographic relief comparable to intermediate to slow spreading ridges here we focus on the northern and central parts of dmh rift where we present quantitative slip rates for the past 40 kyr for major and minor normal fault scarps in the vicinity of a recent september 2005 dike intrusion the data obtained show that the axial valley topography has been created by enhanced slip rates that occurred during periods of limited volcanism suggestive of reduced magmatic activity probably in association with changes in strain distribution in the crust our results indicate that the development of the axial valley topography has been regulated by the lifetimes of the magma reservoirs and their spatial distribution along the segment and thus to the magmatic cycles of replenishment differentiation 100 kyr our findings are also consistent with magma induced deformation in magma rich rift segments the record of two tectonic events of metric vertical amplitude on the fault that accommodated the most part of surface displacement during the 2005 dike intrusion suggests that the latter type of intrusion occurs roughly every 10 kyr in the northern part of the dmh segment

earthquake classification location and error analysis in a volcanic environment implications for the magmatic system of the 1989 1990 eruptions at redoubt volcano alaska abstract determination of the precise locations of seismic events associated with the 1989 1990 eruptions of redoubt volcano posed a number of problems including poorly known crustal velocities a sparse station distribution and an abundance of events with emergent phase onsets in addition the high relief of the volcano could not be incorporated into the hypoellipse earthquake location algorithm this algorithm was modified to allow hypocenters to be located above the elevation of the seismic stations the velocity model was calibrated on the basis of a posteruptive seismic survey in which four chemical explosions were recorded by eight stations of the permanent network supplemented with 20 temporary seismographs deployed on and around the volcanic edifice the model consists of a stack of homogeneous horizontal layers setting the top of the model at the summit allows events to be located anywhere within the volcanic edifice detailed analysis of hypocentral errors shows that the long period lp events constituting the vigorous 23 hour swarm that preceded the initial eruption on december 14 could have originated from a point 1 4 km below the crater floor a similar analysis of lp events in the swarm preceding the major eruption on january 2 shows they also could have originated from a point the location of which is shifted 0 8 km northwest and 0 7 km deeper than the source of the initial swarm we suggest this shift in lp activity reflects a northward jump in the pathway for magmatic gases caused by the sealing of the initial pathway by magma extrusion during the last half of december volcano tectonic vt earthquakes did not occur until after the initial 23 hour long swarm they began slowly just below the lp source and their rate of occurrence increased after the eruption of 01 52 ast on december 15 when they shifted to depths of 6 to 10 km after january 2 the vt activity migrated gradually northward this migration suggests northward propagating withdrawal of magma from a plexus of dikes and or sills located in the 6 to 10 km depth range precise relocations of selected events prior to january 2 clearly resolve a narrow steeply dipping pencil shaped concentration of activity in the depth range of 1 7 km which illuminates the conduit along which magma was transported to the surface a third event type named hybrid which blends the characteristics of both vt and lp events originates just below the lp source and may reflect brittle failure along a zone intersecting a fluid filled crack the distribution of hybrid events is elongated 0 2 0 4 km in an east west direction this distribution may offer constraints on the orientation and size of the fluid filled crack inferred to be the source of the lp events

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

topographic controls on dike injection in volcanic rift zones abstract dike emplacement in volcanic rift zones is often associated with the injection of blade like dikes which propagate long distances parallel to the rift but frequently remain trapped at depth and erupt only near the tip of the dike over geologic time this style of dike injection implies that a greater percentage of extension is accommodated by magma accretion at depth than near the surface in this study we investigate the evolution of faulting topography and stress state in volcanic rift zones using a kinematic model for dike injection in an extending 2 d elastic viscoplastic layer we show that the intrusion of blade like dikes focuses deformation at the rift axis leading to the formation of an axial rift valley however flexure associated with the development of the rift topography generates compression at the base of the plate if the magnitude of these deviatoric compressive stresses exceeds the deviatoric tensile stress associated with far field extension further dike injection will be inhibited in general this transition from tensile to compressive deviatoric stresses occurs when the rate of accretion in the lower crust is greater than 50 60 of the far field extension rate these results indicate that over geologic timescales the injection of blade like dikes is a self limiting process in which dike generated faulting and topography result in an efficient feedback mechanism that controls the time averaged distribution of magma accretion within the crust

crustal structure and magmato tectonic processes in an active rift asal ghoubbet afar east africa 2 insights from the 23 year recording of seismicity since the last rifting event 1 the asal ghoubbet ag rift has sustained a major volcano tectonic rifting episode in 1978 and has been subsequently monitored with continuous geodetic and seismological surveys it is therefore an ideal place to study the transient magmato tectonic processes that operate after a rifting episode we examine the space time evolution of 2500 md 2 8 earthquakes recorded in the rift from 1979 to 2001 we focus on the relationships between this seismic activity and both the three dimensional structure of the rift and its postrifting behavior depicted from geodesy the results highlight the major role of the central magmatic system fieale shark bay on the structure seismic activity and overall behavior of the rift from 1978 to 1986 the rift opens at a fast rate yet mainly aseismically the opening is magmatically driven and accommodated since 1986 when the opening rate abruptly decreased the seismicity is concentrated in the central part of the rift and reveals pulses of activity of the central volcanic system these pulses result from that magmatic zone undergoing alternating stretching and inflating episodes thus while the plate driven induced stresses have been rebuilding in the rift since 1987 the rift opening is still essentially accommodated in the axial magmatic zone the ag rift has thus sustained a postrifting unsteady opening over more than 23 years following its stretching episode that transient opening has essentially occurred aseismically and most tectonic faults remain relaxed and locked