mechanics of low stress forearcs nankai and cascadia the nankai and cascadia subduction zones have low stress forearcs where the margin normal horizontal compressive stress is similar to or less than the vertical stress except in the accretionary prisms in the absence of active back arc spreading the forearc margin normal stress is controlled mainly by the total shear force along the subduction fault or the plate coupling force and the gravitational force the plate coupling force creates compression and the gravitational force in the presence of margin topography creates lateral tension we give a simple expression that relates the forearc stress to topography and average fault strength with the latter represented by an effective coefficient of friction we also use a finite element model of two converging plates in frictional contact to model the forearc stresses the analyses indicate that a low stress forearc requires not only the subduction fault to be very weak but also the mechanically coupled fault area to be narrow in the downdip direction and shallow in depth for nankai and cascadia where the coupled fault areas are narrow and shallow because of the young and hot subducting plates is about 0 05 or less and in any case less than 0 09 the low fault strength is probably not a unique property of low stress forearcs since the same low values combined with wider and deeper coupled areas can be shown to create large margin normal compression in addition there does not seem to be any relation between subduction earthquake magnitudes and plate coupling forces great earthquakes occur at subduction zones having low stress forearcs such as nankai and cascadia

estimates of stress directions by inversion of earthquake fault plane solutions in sicily summary the tectonic stress orientation is estimated in the lithosphere of northern sicily the southern tyrrhenian sea and southern calabria and in the wadati benioff zone below the tyrrhenian by inversion of fault plane solutions of earthquakes covering a magnitude range from 2 5 to 7 1 focal mechanisms of 97 earthquakes are taken from the literature after a critical evaluation of their data quality an average misfit of f 13 indicates that the set of all shallow 50 km earthquakes is generated by a heterogeneous stress field for three subsets based on regional and magnitude separation f was small enough 2 8 lfs5 9 to support the assumption of a homogeneous stress direction for an additional subset with f 7 4 such a condition is close to being fulfilled even though some heterogeneity appears to be present the number of earthquakes in these subsets ranged from nine to 22 and the uncertainties of the principal stress directions were generally of the order of 20 at the 90 per cent confidence level the earthquakes with a425 define a regional stress field with the greatest principal stress ul dipping at a shallow angle to the south in north eastern sicily and south western calabria the stress field estimated by earthquakes is extensional with o3 in a direction of wnw and a near vertical ol in agreement with the graben tectonics mapped geologically in this area in western sicily the c1 direction is oriented wnw but this result is judged less reliable than the others based on the broader confidence limits of the solution and the average misfit of 7 4 the earthquakes in the wadati benioff zone define ol dipping at about 70 to the nw subparallel to the zone with o2 horizontal and striking parallel to the zone

the active portion of the campi flegrei caldera structure imaged by 3 d inversion of gravity data 1 we present an improved density model and a new structural map of the neapolitan yellow tuff caldera the active portion of the nested campi flegrei caldera the model was built using a new 3 d inversion of the available high precision gravity data and a new digital terrain and marine model the inversion procedure based on a variable depth lumped assembling of the subsurface gravity distribution via cell aggregation gives better defined insights into the internal caldera architecture that well agree with the available geological geophysical and geochemical data the adopted 3 d gravity method is highly efficient for characterizing the shallow caldera structure down to 3 km depth and defining features related to regional or volcano tectonic lineaments and dynamics in particular the resulting density distribution highlights a pronounced density low in correspondence of the central portion of the caldera with a detail not available till now the joint interpretation of the available data suggests a subsurface structural setting that supports a piecemeal collapse of the caldera and allows the identification of its headwall positive gravity anomalies localize dense intrusions presently covered by late volcanic deposits along the caldera marginal faults and the main structural lineaments both bordering the resurgent block and cutting the caldera floor these results allow us to both refine the current geological structural framework and propose a new structural map that highlights the caldera boundary and its internal setting this map is useful to interpret the phenomena occurring during unrest and to improve both short term and long term volcanic hazards assessment

conceptual and numerical models of ring fault formation abstract most ring faults of collapse calderas are primarily shear fractures the initiation and development of which depends on the state of stress in the host rock the state of stress in a volcano is controlled by the loading conditions such as magma chamber geometry and pressure but also by the mechanical properties of its rock units and structures such as existing contacts faults and joints ring fault formation is thus essentially a problem in rock physics field observations show that some ring faults are dip slip whereas others are partly faults shear fractures and partly ring dykes extension fractures although slip on existing ring faults is much more common in basaltic edifices shield volcanoes than in true composite volcanoes in both types of volcanoes most caldera unrest periods do not result in ring fault slip here i present new conceptual and numerical models of caldera formation in volcanoes with shallow spherical circular or oblate ellipsoidal sill like magma chambers in the layered models the host rock above the chamber is composed of 30 comparatively thin layers with stiffnesses young s moduli alternating between 1 gpa and 100 gpa the chamber itself is located in a single thick layer the crustal segment hosting the chamber is either 20 km or 40 km wide but has a constant thickness of 20 km the loading conditions considered are 1 a crustal segment subject to 5 mpa tension 2 crustal segment subject to excess magmatic pressure of 10 mpa at the bottom doming of the volcanic field containing the chamber 3 a combination of tension and doming and 4 chamber subject to underpressure negative excess pressure of 5 mpa the main results are as follows 1 excess pressure and underpressure in a chamber normally favour dyke injection rather than ring fault formation 2 for doming or tension a spherical magma chamber favours dyke injection except when the layer hosting the chamber is very soft 10 gpa or one with recent dyke injections in which case the surface stress field favours ring fault formation 3 for a sill like chamber in a 20 km wide crustal segment a ring fault can be generated by either tension or tension and doming for a 40 km wide segment doming alone is sufficient to generate a ring fault 4 since individual layers in a volcano may develop different local stresses stress field homogenisation through all the layers between the chamber and the surface is a necessary condition for ring fault formation 5 because the mechanical properties of the layers that constitute basaltic edifices are more uniform than those that constitute true composite volcanoes it follows that stress field homogenisation and thus ring fault formation or slip is more commonly reached in basaltic edifices than in composite volcanoes 6 both for basaltic edifices and composite volcanoes the stress fields most likely to initiate ring faults are those generated around sill like chambers subject to tension doming or both

aftershock seismicity of the 27 february 2010 mw 8 8 maule earthquake rupture zone on 27 february 2010 the mw 8 8 maule earthquake in central chile ruptured a seismic gap where significant strain had accumulated since 1835 shortly after the mainshock a dense network of temporary seismic stations was installed along the whole rupture zone in order to capture the aftershock activity here we present the aftershock distribution and first motion polarity focal mechanisms based on automatic detection algorithms and picking engines by processing the seismic data between 15 march and 30 september 2010 from stations from iris ipgp gfz and university of liverpool we determined 20 205 hypocentres with magnitudes mw between 1 and 5 5 seismic activity occurs in six groups 1 normal faulting outer rise events 2 a shallow group of plate interface seismicity apparent at 25 35 km depth and 50 120 km distance to the trench with some variations between profiles along strike the aftershocks occur largely within the zone of coseismic slip but extend 50 km further north and with predominantly shallowly dipping thrust mechanisms along dip the events are either within the zone of coseismic slip or downdip from it depending on the coseismic slip model used 3 a third band of seismicity is observed further downdip at 40 50 km depth and further inland at 150 160 km trench perpendicular distance with mostly shallow dipping 28 thrust focal mechanisms indicating rupture of the plate interface significantly downdip of the coseismic rupture and presumably above the intersection of the continental moho with the plate interface 4 a deep group of intermediate depth events between 80 and 120 km depth is present north of 36 s within the maule segment a large portion of events during the inter seismic phase originated from this depth range 5 the magmatic arc exhibits a small amount of crustal seismicity but does not appear to show significantly enhanced activity after the mw 8 8 maule 2010 earthquake 6 pronounced crustal aftershock activity with mainly normal faulting mechanisms is found in the region of pichilemu 34 5 s these crustal events occur in a 30 km wide region with sharp inclined boundaries and oriented oblique to the trench the best located events describe a plane dipping to the southwest consistent with one of the focal planes of the large normal faulting aftershock mw 6 9 on 11 march 2010

the july august 2000 seismic swarms at campi flegrei volcanic complex italy two swarms of microearthquakes ml 2 2 occurred on july 2 7 and august 22 2000 at campi flegrei accompained by a ground uplift episode 4 cm which interrupted on early march 2000 the descending trend started on 1985 spectral analysis indicates a direct involvement of magmatic hydrothermal fluids in the source process of the july swarm while the august events are typical of shear failure similar to most of the earthquakes that occurred during the last 1982 1984 bradyseismic crisis precise 3 d relative location applied to similar earthquakes allows for the recognition of two parallel alignments trending ne sw at depths of 1 7 and 3 2 km this trend is consistent with the direction of the main focal plane obtained from fault plane solutions and evidences tensile failure in close proximity to the zone of maximum uplift as depicted by geodetic measurements a fault weakening mechanism triggered by increasing pore pressure is invoked as the cause of these earthquakes

a scattered wave image of subduction beneath the transverse ranges over 5600 short period recordings of teleseismic events were used to create deterministic maps of p wave scatterers in the upper mantle beneath southern california between depths of 50 and 200 kilometers the southern flank of the slab subducting beneath the transverse ranges is marked by strong scattering the marked scattering indicates that the edge of the slab is a sharp thermal boundary such a boundary could be produced by slab shearing or small scale convection in the surrounding mantle the northern limb of the slab is not a strong scatterer consistent with thicker lithosphere north of the transverse ranges

hydrothermal circulation in subducting crust reduces subduction zone temperatures most thermal models of subduction zones assume no advection of heat by fluid flow because slow flow through underthrusting sediment the decollement and wedge likely transports only a minor amount of heat we model coupled fluid and heat transport in a subduction zone and show that hydrothermal circulation in subducting basaltic basement rocks can have a great influence on subduction zone temperatures fractured basaltic basement has permeability several orders of magnitude higher than a typical decollement allowing fluid circulation to redistribute and extract heat from a subduction zone we simulate systems with upper basaltic basement permeability ranging from 10 13 to 10 10 m2 in addition we incorporate the effect of permeability reduction within the basaltic basement as it is subducted the models with fluid transport show suppressed temperatures along the subducting slab relative to models with no fluid transport with continuous sediment cover heat is extracted from under the margin wedge to the trench in models where faulted ocean crust exposes high permeability basement to the ocean floor cooling from ocean bottom water results in highly suppressed heat flow relative to conductive models hydrothermally cooled ocean crust also acts to slow thermally controlled diagenetic reaction progress within subducting sediment

age spreading rates and spreading asymmetry of the world s ocean crust we present four companion digital models of the age age uncertainty spreading rates and spreading asymmetries of the world s ocean basins as geographic and mercator grids with 2 arc min resolution the grids include data from all the major ocean basins as well as detailed reconstructions of back arc basins the age spreading rate and asymmetry at each grid node are determined by linear interpolation between adjacent seafloor isochrons in the direction of spreading ages for ocean floor between the oldest identified magnetic anomalies and continental crust are interpolated by geological estimates of the ages of passive continental margin segments the age uncertainties for grid cells coinciding with marine magnetic anomaly identifications observed or rotated to their conjugate ridge flanks are based on the difference between gridded age and observed age the uncertainties are also a function of the distance of a given grid cell to the nearest age observation and the proximity to fracture zones or other age discontinuities asymmetries in crustal accretion appear to be frequently related to asthenospheric flow from mantle plumes to spreading ridges resulting in ridge jumps toward hot spots we also use the new age grid to compute global residual basement depth grids from the difference between observed oceanic basement depth and predicted depth using three alternative age depth relationships the new set of grids helps to investigate prominent negative depth anomalies which may be alternatively related to subducted slab material descending in the mantle or to asthenospheric flow a combination of our digital grids and the associated relative and absolute plate motion model with seismic tomography and mantle convection model outputs represents a valuable set of tools to investigate geodynamic problems

abrupt thermal transition reveals hydrothermal boundary and role of seamounts within the cocos plate 1 new thermal data from 18 24 ma lithosphere on the cocos plate delineate contrasting subsurface thermal conditions in adjacent sections of crust heat flow through seafloor created at the east pacific rise is generally suppressed by 70 relative to conductive lithospheric cooling models whereas heat flow through adjacent similarly aged lithosphere generated at the cocos nazca spreading center is consistent with these models the transition between thermal regimes is remarkably abrupt only 2 5 km wide indicating a shallow hydrothermal origin the transition is more closely associated with differences in the distribution of basement outcrops than with tectonic boundaries demonstrating the importance of the former in extracting heat from the lithosphere on a regional basis

the restless resurgent campi flegrei nested caldera italy constraints on its evolution and configuration abstract we have collected new stratigraphic data in the campi flegrei caldera through surface geological mapping and study of cores drilled to various depths we integrate these data with geological geochronological and geophysical constraints available in the literature to reconstruct the volcanic and deformational history of the caldera the caldera is a resurgent nested structure resulting from two main collapses related to the campanian ignimbrite 37 ka and the neapolitan yellow tuff 12 ka eruptions while the structure is affected by a broad subsidence ongoing local resurgence and unrest occur inside the younger nested caldera structure although most of the caldera floor consists of multiple structural blocks due to long term deformations magmas were able to erupt only in one sector of the resurgent part the caldera has shown signs of unrest during the past 25 years with short term deformations that have generated a net uplift of 3 5 m the short term deformations are related to the neapolitan yellow tuff caldera resurgence and interpreted as the net result of brittle and ductile components the long term deformations likely represent the sum of the permanent brittle component of each short term deformational event

structure and evolution of the bay of pozzuoli italy using marine seismic reflection data implications for collapse of the campi flegrei caldera digital marine seismic reflection data acquired in 1973 in the bay of pozzuoli and recently reprocessed were used to study the volcanological evolution of the marine sector of campi flegrei caldera during the last 37 ka in order to gain more information interpretation also involved estimation of the pseudo velocity and the pseudo density from the resistivity logs of two onshore deep exploration wells the main results are 1 discovery of ancient pre 18 ka and post 37 ka submarine and mainly effusive volcanic activity along coeval emission centers located at the edges of campi flegrei caldera 2 confirmation that the caldera collapse in the marine sector of campi flegrei seems strongly controlled by regional ne sw and nw se structural discontinuities 3 the finding of at least two episodes of collapse in the bay and 4 identification of a post 18 ka volcanic deflation phase that has caused about 150 200 m of subsidence in the central sector of the bay of pozzuoli in the last 18 ka

seismicity and geometry of the south chilean subduction zone 41 5 s 43 5 s implications for controlling parameters in 2005 an amphibious seismic network was deployed on the chilean forearc between 41 75 s and 43 25 s 364 local events were observed in a 11 month period a subset of the p and s arrival times were inverted for hypocentral coordinates 1 d velocity structure and station delays main seismic activity occurred predominantly in a belt parallel to the coast of chiloe island in a depth range of 12 30 km presumably related to the plate interface the 30 inclination of the shallow part of the wadati benioff zone is similar to observations further north indicating that oceanic plate age is not controlling the subduction angle of the shallower part for the chilean subduction zone the down dip termination of abundant intermediate depth seismicity at approximately 70 km depth seems to be related to the young age and high temperature of the oceanic plate crustal seismicity is associated with the liquine ofqui fault zone and active volcanoes

p sv conversions at a shallow boundary beneath campi flegrei caldera italy evidence for the magma chamber seismograms from an active seismic experiment carried out at campi flegrei caldera near naples italy show a large amplitude 5v polarized shear wave following by less than 1 5 s p waves reflected at wide angle from a deep crustal interface early arriving sv polarized waves with the same delay to direct p waves are also observed in seismograms from a regional 280 km deep magnitude 5 1 earthquake such short delays of s to p waves are consistent with a p sv conversion on transmission occurring at a shallow boundary beneath the receivers the large amplitude of the converted sv phase along with that the p waves are near vertical requires a boundary separating a very low rigidity layer from the upper caldera fill the converted phases are interpreted as a seismic marker of a magma chamber the top of this magma chamber is located slightly deeper than the deepest earthquakes observed during the 1982 1984 unrest of campi flegrei

remarkably consistent thermal state of the south central chile subduction zone from 36 s to 45 s delineating rupture areas of subduction zone earthquakes is necessary for understanding the controls on seismic and aseismic slip for the largest recorded earthquake the 1960 chile event with moment magnitude 9 5 the rupture area is only loosely defined due to limitations in the global seismic network at the time the rupture extends 900 km along strike coastal deformation is consistent with either a constant rupture width of 180 200 km along the entire length or a narrower 115 km rupture in the southern half a southward narrowing of the seismogenic zone has been hypothesized to result from warming of the subduction zone to the south where the subducting plate is younger we present results of thermal models at 36 s 38 s 43 s and 45 s to examine potential along strike changes in thermal state models most consistent with observed surface heat flux include fluid circulation in the oceanic crust that advects heat to the ocean this ventilated hydrothermal circulation preferentially cools transects with young subducting lithosphere frictional heating preferentially warms transects with older subducting lithosphere the combined effects of frictional heating and hydrothermal circulation increase decollement temperatures in the 36 s and 38 s transects by up to 155 c and decrease temperatures in the 45 s transect by up to 150 c in our preferred models decollement temperatures 200 km landward of the trench in all four transects are 350 400 c this is consistent with a constant 200 km wide seismogenic zone for the 1960 mw 9 5 rupture with decreasing slip magnitude in the southern half of the rupture

deep crustal heterogeneity along and around the san andreas fault system in central california and its relation to the segmentation the three dimensional distribution of scatterers in the crust along and around the san andreas fault system in central california is estimated using an inversion analysis of coda envelopes from local earthquakes i analyzed 3801 wave traces from 157 events recorded at 140 stations of the northern california seismic network the resulting scatterer distribution shows a correlation with the san gregorio san andreas hayward and calaveras faults these faults seem to be almost vertical from the surface to 15 km depth some of the other scatterers are estimated to be at shallow depths 0 5 km below the diablo range and these may be interpreted as being generated by topographic roughness the depth distribution of scatterers shows relatively stronger scattering in the lower crust at 15 25 km depth especially between the san andreas fault and the hayward calaveras faults this suggests a subhorizontal detachment structure connecting these two faults in the lower crust several clusters of scatterers are located along the san andreas fault at intervals of 20 30 km from south of san francisco to the intersection with the calaveras fault this part of the san andreas fault appears to consist of partially locked segments also 20 30 km long which rupture during m6 7 events and segment boundaries characterized by stronger scattering and stationary microseismicity the segment boundaries delineated by the present analysis correspond with those estimated from the slip distribution of the great 1906 san francisco earthquake and from the fault geometry as reported by the working group on california earthquake probabilities 1990 although the segment boundaries along the san andreas fault in and around the san francisco bay area are still uncertain

the chile ridge a tectonic framework a new chile ridge tectonic framework is developed based on satellite altimetry data shipboard geophysical data and primarily 38 500 km of magnetic data gathered on a joint u s chilean aeromagnetic survey eighteen active transforms with fossil fracture zones fzs including two complex systems the chile fz and valdivia fz systems have been mapped between the northern end of the antarctic nazca plate boundary chile ridge at 35 s and the chile margin triple junction at 47 s chile ridge spreading rates from 23 ma to present have been determined and show slowdowns in spreading rates that correspond to times of ridge trench collisions the valdivia fz system previously mapped as two fzs with an uncharted seismically active region between them is now recognized to be a multiple offset fz system composed of six fzs separated by short ridge segments 22 to 27 km long at chron 5a 12 ma the chile ridge propagated from the valdivia fz system northward into the nazca plate through crust formed 5 myr earlier at the pacific nazca ridge evidence for this propagation event includes the friday and crusoe troughs located at discontinuities in the magnetic anomaly sequence and interpreted as pseudofaults this propagation event led to the formation of the friday microplate which resulted in the transferal of crust from the nazca plate to the antarctic plate and in a 500 km northward stepwise migration of the pacific antarctic nazca triple junction rift propagation microplate formation microplate extinction and stepwise triple junction migration are found to occur during large scale plate motion changes and plate boundary changes in the southeast pacific

magma chamber of the campi flegrei supervolcano at the time of eruption of the campanian ignimbrite a supereruption that occurred in the campi flegrei area italy ca 39 ka had regional and global scale environmental impacts and deposited the campanian ignimbrite ci we attempt to shed light on critical aspects of the eruption depth of magma chamber intensive pre eruptive magma conditions and the large volume magma plumbing system on the basis of information derived from analyzing melt inclusion mi data to achieve these aims we provide new measurements of homogenization temperatures and values of dissolved h2o within phenocryst hosted mis from pumices erupted during different phases of the ci eruption the mi data indicate that a relatively homogeneous overheated trachytic magma resided within a relatively deep magma chamber dissolved water contents in mis indicate that prior to the eruption the magma chamber underwent radical changes related to differential upward movement of magma decompression of the rising trachytic magma caused a decrease in water solubility and crystallization and trachytic bodies were emplaced at very shallow depths the proposed eruptive model links portions of the main magma chamber and apophyses with specific eruptive units

heat flow along the nantroseize transect results from iodp expeditions 315 and 316 offshore the kii peninsula japan mal conductivity varies between about 1 0 w m 1 k 1 near the seafloor to 1 6 w m 1 k 1 at a depth of 1 km thermal conductivity generally increases with depth and correlates with variations in porosity and lithology thermal gradients along the transect are characterized by 48 sediment temperature measurements from 6 sites thermal corrections for the effects of bathymetric relief and sedimentation improve the confidence with which the heat flow values can be interpreted heat flow generally decreases with landward distance from the deformation front and varies from 70 mw m 2 just landward of the deformation front to 54 mw m 2 at sites characterizing the outer fore arc high and to 57 mw m 2 at the kumano basin site iodp heat flow measurements are significantly lower than nearby seafloor heat flow measurements this difference is most likely due to variations in bottom water temperature that have a large effect on values of seafloor heat flow thus the heat flow of the nankai accretionary prism is lower than previously thought we present thermal models of subduction along this transect and explore the impact of the initial geotherm conductive plate cooling based on the age of subducting seafloor 20 myr under predicts the observed heat flow we find a good fit to the data using a geotherm appropriate for 10 myr seafloor the extra heat is interpreted in terms of back arc thermal environments

large heat and fluid fluxes driven through mid plate outcrops on ocean crust geophysical data for the cocos plate sea floor suggest that basement outcrops along mid ocean ridge flanks can discharge very large quantities of heat and fluid this is indicative of high crustal permeability at the regional scale

deep lithospheric structures along the southern central chile margin from wide angle p wave modelling p crustal and upper mantle structures of the subduction zone in south central chile between 42 degrees s and 46 degrees s are determined from seismic wide angle reflection and refraction data using the seismic ray tracing method to calculate minimum parameter models three profiles along differently aged segments of the subducting nazca plate were analysed in order to study subduction zone structure dependencies related to the age that is thermal state of the incoming plate the age of the oceanic crust at the trench ranges from 3 ma on the southernmost profile immediately north of the chile triple junction to 6 5 ma old about 100 km to the north and to 14 5 ma old another 200 km further north off the island of chiloe remarkable similarities appear in the structures of both the incoming as well as the overriding plate the oceanic nazca plate is around 5 km thick with a slightly increasing thickness northward reflecting temperature changes at the time of crustal generation the trench basin is about 2 km thick except in the south where the chile ridge is close to the deformation front and only a small 800 m thick trench infill could develop in south central chile typically three quarters 1 5 km of the trench sediments subduct below the decollement in the subduction channel to the north and south of the study area only about one quarter to one third of the sediments subducts the rest is accreted above similarities in the overriding plate are the width of the active accretionary prism 35 50 km and a strong lateral crustal velocity gradient zone about 75 80 km landward from the deformation front where landward upper crustal velocities of over 5 0 5 4 km s su 1 su decrease seaward to around 4 5 km s su 1 su within about 10 km which possibly represents a palaeo backstop this zone is also accompanied by strong intraplate seismicity differences in the subduction zone structures exist in the outer rise region where the northern profile exhibits a clear bulge of uplifted oceanic lithosphere prior to subduction whereas the younger structures have a less developed outer rise this plate bending is accompanied by strongly reduced rock velocities on the northern profile due to fracturing and possible hydration of the crust and upper mantle the southern profiles do not exhibit such a strong alteration of the lithosphere although this effect may be counteracted by plate cooling effects which are reflected in increasing rock velocities away from the spreading centre overall there appears little influence of incoming plate age on the subduction zone structure which may explain why the m w 9 5 great chile earthquake from 1960 ruptured through all these differing age segments the rupture area however appears to coincide with a relatively thick subduction channel

the morphologies of volcanic landforms at central elysium planitia evidence for recent and fluid lavas on mars abstract this paper focuses on physical parameters flow rates and rheological properties of lava flows observed in the central elysium planitia cep region of mars the flows are modeled as newtonian fluids using the jeffrey s equation and the concept of graetz number or alternatively as bingham fluids in addition to these approaches a theoretical model of the shape of shield volcanoes based on the solution for the porous flow of an unconfined aquifer is applied to 5 shields providing independent quantifications of rheological variations between the shields this analysis indicates that of the five volcanoes studied two are partially buried by lava postdating their formation a result which has been confirmed independently in one case by high resolution images our observations reveal that two types of lava flows may be found in the cep region the first group is composed of large lava flows with viscosities around 2 5 10 5 pa s or yield strengths ranging from 100 to 500 pa the second group includes small lava flows of the shield volcanoes and large leveed lava channels on the plains with viscosities below 103 pa s or yield strengths less than 200 pa when compared with other volcanic regions on mars investigated with similar approaches these latter values are at present the lowest inferred for martian lava flows several hypotheses for the formation of these lavas are discussed in the context of cep given that low viscosity can be the result of 1 high temperature 2 low crystal content 3 low si abundance of the liquid phase and or 4 the presence of dissolved volatiles two scenarios are considered in the first one it is demonstrated that low viscosity lavas of low silica content can be produced in the context proposed by schumacher and breuer schumacher s breuer d 2007 geophys res lett 34 l12202 for recent volcanism however geochemical maps derived from grs measurements do not provide support for anomalously low silica concentrations in this region in the second scenario a water rich magma is proposed although arguments in favor of a water rich mantle source below the cep are not available at the present time

the 1960 chile earthquake inversion for slip distribution from surface deformation summary a total of 166 observations of sea level change 130 measurements of elevation difference and 16 determinations of horizontal strain provide an excellent view of the quasi static source process of the great 1960 chilean earthquake these surface deformation data were employed in classical uniform slip fault models as well as more recently developed models that allow spatial variability of slip the best uniform slip planar usp model is 850km long 130km wide and dips 20 seventeen metres of fault displacement contributed to a usp moment of 9 4 times 1022 n m the variable slip planar vsp model concentrates slip on a 900 km long 150 km wide band parallel to the coast several peaks of slip with dimensions of 50 100 km appear in this band and are thought to represent major subduction zone asperities important fractures of the oceanic lithosphere bound the 1960 rupture and are offered as a potential source of fault segmentation within the chilean subduction zone the vsp moment for 1960 earthquake totals 9 5 times 1022 n m about one fifth of the value estimated for the foreshock mainshock sequence from seismic methods except for areas out to sea geodetic resolution on the fault is fairly uniform thus it is unlikely that slip missed by the network could increase the vsp moment much beyond 1 8 times 1023 n m several patches of moment isolated from the main body at 80 110 km depth are found down dip in the vsp model and are presumably indicative of aseismic slip one patch at the northern end of the rupture is probably associated with the initiation phase of the mainshock although the time sequence of the relationship is unknown tide gauge records suggest that another patch between 40 and 43 s responsible for the observed strain and uplifts inland at those latitudes is not of coseismic origin but derives from in place post seismic creep over several years apparently great 1960 type events are not typical members of the 128 yr earthquake cycle in south central chile the nazca south america boundary here is characterized by a variable rupture mode in which major asperities are completely broken by great earthquakes only once in four or five earthquake cycles the more frequent large earthquakes that geographically overlap the great events fill in between the locked zones

the age of the neapolitan yellow tuff caldera forming eruption campi flegrei caldera italy assessed by 40ar 39ar dating method abstract the neapolitan yellow tuff nyt is the product of the largest known trachytic phreatoplinian eruption it covered an area larger than 1000 km2 with an estimated volume of about 40 km3 of erupted magma during the course of the eruption a caldera collapsed within the previously formed campanian ignimbrite caldera the resulting nested structure strongly influenced the following volcanic activity in the campi flegrei caldera as previous dating of the nyt does not converge toward a unique result a new set of 40ar 39ar age determinations has been carried out to better constrain the age of the eruption two variants of the 40ar 39ar dating method were applied to determine the age of the nyt eruption 1 single crystal total fusion sctf on an individual phenocryst of feldspar and 2 laser incremental heating lih on bulk aliquots of feldspar phenocrysts the results of the sctf analyses show that the overall sample weighted mean age derived from the conventional age calculation is 15 6 0 8 ka a weighted mean of the isochron age is 15 3 1 2 ka 2 and has been assumed as the best indicator of age to be derived from the sctf analyses the lih analyses results show that plateau ages vary from 15 4 0 5 to 14 5 0 5 ka the overall weighted mean age of the isochron results is 14 9 0 4 ka 2 this result has been assumed as the reference age for the nyt eruption and agrees with the sctf age the new age obtained for the nyt deposits is of great relevance for the understanding of the evolution and the present state of the campi flegrei caldera and collocates the nyt in a crucial stratigraphical position to date the climatic oscillations that occurred between the late glacial and the holocene

three dimensional model of the feeder conduit of the 1983 eruption of mt etna volcano from ground deformation measurements a levelling network and a horizontal distance measuring network down mt etna s south flank were measured before during and after the 1983 eruption we analyse the large movements recorded using theoretical models of surface displacement to determine the shape and position of the feeder conduit that supplied magma to the eruption the results indicate that a dyke the top surface extending nearly horizontally connects the eruption site to a point near the southeast crater the top of the dyke lies at an altitude of 2400 m to 2500 m and dips between 75 and 90 to the west more tentative evidence indicates that the top of the active dyke dropped by about 150 m between the 25th and the 95th day of the eruption and that a small branching dyke 1 km long was emplaced to the east of the main dyke but did not attain the surface the implications of these observation are discussed

the southern andes between 36 and 40 s latitude seismicity and average seismic velocities abstract the project issa 2000 integrated seismological experiment in the southern andes consists of a temporary seismological network and a seismic refraction profile a network of 62 seismological stations was deployed across the southern andes at 38 s three hundred thirty three local seismic events were observed in a 3 month period p and s arrival times of a subset of high quality data were inverted simultaneously for 1 d velocity structure hypocentral coordinates and station delays seismic refraction data along a transect at 39 s provide further constraints on the crustal structure low crustal velocities beneath the forearc may be either due to subducted trench sediments or serpentinized mantle material of the continental lithosphere the continental moho is not clearly observed in this region average velocities of the crust beneath the arc are higher than those beneath the forearc crustal thickness is about 40 km crustal seismicity concentrates in the forearc region along the bio bio and gastre fault zones the area between these two prominent fault zones seems to be nearly devoid of crustal seismicity but shows highest uplift and topography in the forearc region benioff seismicity is observed down to 150 km depth resulting in the first accurate image of the benioff zone in the southern andes a maximum of seismicity at 60 km depth may be caused by dehydration embrittlement

tectonic control on sediment accretion and subduction off south central chile implications for coseismic rupture processes of the 1960 and 2010 megathrust earthquakes based on a compilation of published and new seismic refraction and multichannel seismic reflection data along the south central chile margin 33 46 s we study the processes of sediment accretion and subduction and their implications on megathrust seismicity in terms of the frontal accretionary prism fap size the marine south central chile fore arc can be divided in two main segments 1 the maule segment south of the juan fernandez ridge and north of the mocha block characterized by a relative large fap 20 40 km wide and 2 the chiloe segment south of the mocha block and north of the nazca antarctic south america plates junction characterized by a small fap 10 km wide in addition the maule and chiloe segments correlate with a thin 1 km thick and thick 1 5 km thick subduction channel respectively the mocha block lies between 37 5 and 40 s and is configured by the chile trench mocha and valdivia fracture zones this region separates young 0 25 ma oceanic lithosphere in the south from old 30 35 ma oceanic lithosphere in the north and it represents a fundamental tectonic boundary separating two different styles of sediment accretion and subduction respectively a process responsible for this segmentation could be related to differences in initial angles of subduction which in turn depend on the amplitude of the down deflected oceanic lithosphere under trench sediment loading on the other hand a small fap along the chiloe segment is coincident with the rupture area of the trans pacific tsunamigenic 1960 earthquake mw 9 5 while a relatively large fap along the maule segment is coincident with the rupture area of the 2010 earthquake mw 8 8 differences in earthquake and tsunami magnitudes between these events can be explained in terms of the fap size along the chiloe and maule segments that control the location of the updip limit of the seismogenic zone the rupture area of the 1960 event also correlates with a thick subduction channel chiloe segment that may provide enough smoothness at the subduction interface allowing long lateral earthquake rupture propagation

moho topography and subducting oceanic slab of the chilean continental margin in the maximum slip segment of the 1960 mw 9 5 valdivia chile earthquake from p receiver functions the maximum slip area of the 1960 mw 9 5 valdivia chile earthquake has been investigated down to 100 km depth by a receiver function survey along two trench perpendicular profiles at 39 3 and 39 5 s the data are compared with seismicity and local earthquake tomography we find a slab of ca 35 dip at depths similar to what is observed to the north and south of the study area the slab and the overlying mantle wedge show a local decrease in seismic velocities possibly caused by hydration or underplating of sediments the presence of hydrated or sedimentary rocks would reduce plate coupling in this region in agreement with recent gps measurements the moho is found at depths of 30 to 40 km deepening from the coast to the andes the moho topography is close to isostatic equilibrium under the volcanic centre of villarrica quetrupillan and lanin the moho conversion occurs about 10 km shallower than expected for a standard isostatic model this can be explained by a density reduction in the lower crust due to magmatic processes or by non compensated stress

activating and reactivating pairs of nested collapses during caldera forming eruptions campi flegrei italy 1 a structural model for campi flegrei caldera cfc italy one of earth s most hazardous calderas is provided data from the north outer rim show that this formed by downsag soon after the neapolitan yellow tuff emplacement 15 ka campanian ignimbrite s cfc s other major eruption 39 ka displacement at depth also suggests the activity of a buried caldera fault therefore during the campanian ignimbrite eruption the north rim was activated as caldera bounding fault and during the neapolitan yellow tuff eruption as downsag previously published data reveal an inner collapsed structure to the south formed during the campanian ignimbrite eruption and reactivated during the neapolitan yellow tuff eruption so cfc consists of 2 nested depressions both active during the campanian ignimbrite and neapolitan yellow tuff eruptions but with different amounts of collapse cfc is a major example of how 2 nested calderas may 1 form during the same collapse episode 2 be reactivated together during caldera forming eruptions

space geodetic estimation of the nazca south america euler vector abstract gps data from four sites in the nazca plate easter island galapagos robinson crusoe and san felix islands and from five sites in the stable core of the south american plate enabled us to estimate the euler vector of the nazca plate with respect to south america the observed velocities of easter island 6 6 cm yr at 102 3 galapagos 5 1 cm yr at 90 0 robinson crusoe 6 6 cm yr at 80 1 and san felix 6 0 cm yr at 82 1 are significantly slower than the global plate model nuvel 1a velocites for those four sites our estimated euler pole is located at 48 8 n 91 7 w with a rate of 0 59 m y the nuvel 1a and earlier space geodetic studies give rotation rates that are 20 faster

geomorphology and structural development of the nested summit crater of lascar volcano studied with terrestrial laser scanner data and analogue modelling abstract many volcano summits host craters that are partially overlapping the formation of such nested craters has been commonly interpreted as vent migration here we present an additional mechanism that may explain the geometry of nested craters at volcanoes lascar volcano the most active volcano of the central volcanic zone in the chilean andes hosts ene wsw trending summit craters that are partially overlapping nested details on the evolution and interaction between the different craters remain unclear to create a robust dataset terrestrial laser scanner tls data were collected at the summit of lascar in 2013 the resulting topographic data set consisting of more than 15 million data points with centimetre sampling allows visualising almost the complete eastern edifice of the volcano s summit from the tls data a digital elevation model dem and a slope map were generated allowing us to create a lineament map and quantify the observed morphological and structural features to further improve our understanding of the processes responsible for the formation of the craters and geomorphology we designed sandbox analogue models results suggest that one of the craters is a parasite crater formed as a consequence of ongoing activity in the adjacent crater our data suggest that the nested craters have all been modified since the last major eruption in 1993 by near surface effects associated with cooling compaction and gravitational sliding of the crater floor infill as the active crater deepens the adjacent inactive crater extends and partially slumps towards the active one understanding the structural development of these nested craters is relevant for assessing potential future eruption sites thus making lascar a dynamic target for a detailed morphology study these findings may similarly be applied to other volcanoes where nested craters have developed

impact of megathrust geometry on inversion of coseismic slip from geodetic data application to the 1960 chile earthquake 1 we analyze the role of megathrust geometry on slip estimation using the 1960 chile earthquake mw 9 5 as an example a variable slip distribution for this earthquake has been derived by barrientos and ward 1990 applying an elastic dislocation model with a planar fault geometry their model shows slip patches at 80 110 km depth isolated from the seismogenic zone interpreted as aseismic slip we invert the same geodetic data set using a finite element model fem with precise geometry derived from geophysical data isoparametric fem is implemented to constrain the slip distribution of curve shaped elements slip resolved by our precise geometry model is limited to the shallow region of the plate interface suggesting that the deep patches of moment were most likely an artifact of the planar geometry our study emphasizes the importance of fault geometry on slip estimation of large earthquakes

effects of frictional behavior and geometry of subduction fault on coseismic seafloor deformation abstract the seismogenic zone of subduction faults appears to have an updip limit seaward of which the fault exhibits velocity strengthening behavior we use a two dimensional finite element model including a frictional subduction fault to explore how coseismic strengthening of the updip segment affects seafloor deformation for a stress drop of a few mpa along the seismogenic zone strengthening of the updip segment by a comparable amount can prevent the rupture from breaking the trench the resultant seafloor uplift is much larger than that predicted by a model of equal seismic moment or maximum fault slip in which the seismogenic zone extends to the trench with a curved fault geometry although a lower degree of coseismic strengthening of the updip segment leads to greater slip in the shallowest part of the fault it produces smaller seafloor uplift contrary to a popular belief given the paucity of direct observations of the coseismic behavior of the updip segment the models yield important information for understanding rupture mechanics of subduction faults and tsunami generation we also illustrate how the results may be used to guide fault slip parameterization in analytical dislocation models

on the geometrical form of volcanoes abstract many volcanic edifices have a remarkably symmetric geometrical form an example is mount fuji in japan we model this form assuming that the surface of the volcano is a surface of uniform hydraulic potential that an erupting magma will follow the path of minimum resistance to the surface in order to model the resistance to fluid flow we assume the volcanic edifice is a uniform porous medium the vertical flow of magma is also resisted by the gravitational body force if the volcano becomes too tall flank eruptions will widen it if the volcano becomes too wide summit eruptions will increase its elevation using the dupuit approximation for an unconfined aquifer it is shown that the percolation equation is applicable as magma reaches the surface it is assumed to extend the solid porous matrix a similarity solution is obtained to this moving boundary problem the solution predicts a uniform shape for all volcanoes this shape is shown to be in excellent agreement with the geometrical form of mount fuji

interpretation of gravity and magnetic data in the phlegraean fields geothermal area naples italy abstract available gravity and magnetic data of the phlegraean fields geothermal area naples italy have been interpreted and the obtained structural models discussed in the light of the other available geological volcanological and geophysical data on the basis of the results of a previous seismic reflection survey in the gulf of naples and in the pozzuoli bay which delineated a basement characterized by a seismic velocity of 4 6 km s it has been possible to evaluate the gravity anomaly connected with the morphology of this horizon 2 7 g cm3 the residual anomaly map obtained after subtraction of the regional long wavelength components relative to mantle and deep crustal structures and the computed components relative to the above mentioned seismic basement shows up as a circular low with an amplitude of 10 mgal centred in the pozzuoli bay this gravity low has been interpreted as due to the occurrence in the centre of pozzuoli bay of light 0 2 g cm3 material with a maximum thickness of about 2 km however a contribution to the anomaly due to a narrow magmatic body intruded in the basement as suggested by volcanological and ground deformation data cannot be excluded the aeromagnetic map of the phlegraean fields is characterized by three main anomalies which have been fitted by superficial tridimensional parallelepipedic bodies schematically representing lava flows and domes their anomalies have been subsequently subtracted from the observed field obtaining as a residual a large anomaly centred in the southwestern area of the pozzuoli bay it has been interpreted as being due to a lowmagnetized body which taking into account the thermal state of the area should represent that part of the pyroclastic sequence which has lost part of its magnetization by thermo chemical alteration

frictional slip of granite at hydrothermal conditions sliding on faults in much of the continental crust likely occurs at hydrothermal conditions i e at elevated temperature and elevated pressure of aqueous pore fluids yet there have been few relevant laboratory studies to measure the strength sliding behavior and friction constitutive properties of faults at hydrothermal conditions we slid laboratory granite faults containing a layer of granite powder simulated gouge velocity stepping experiments were performed at temperatures of 23 to 600 c pore fluid pressures ph2o of 0 dry and 100 mpa wet effective normal stress of 400 mpa and sliding velocities v of 0 01 to 1 m s 0 32 to 32 m yr conditions were similar to those in earlier tests on dry granite to 845 c by lockner et al 1986 the mechanical results define two regimes the first regime includes dry granite up to at least 845 and wet granite below 250 c in this regime the coefficient of friction is high 0 7 to 0 8 and depends only modestly on temperature slip rate and ph2o the second regime includes wet granite above 350 c in this regime friction decreases considerably with increasing temperature temperature weakening and with decreasing slip rate velocity strengthening these regimes correspond well to those identified in sliding tests on ultrafine quartz we infer that one or more fluid assisted deformation mechanisms are activated in the second hydrothermal regime and operate concurrently with cataclastic flow slip in the first cool and or dry regime is characterized by pervasive shearing and particle size reduction slip in the second hot and wet regime is localized primarily onto narrow shear bands adjacent to the gouge rock interfaces weakness of these boundary shears may result either from an abundance of phyllosilicates preferentially aligned for easy dislocation glide or from a dependence of strength on gouge particle size major features of the granite data set can be fit reasonably well by a rate and temperature dependent three regime friction constitutive model chester this issue we extrapolate the experimental data and model fit in order to estimate steady state shear strength versus depth along natural slipping faults for sliding rates as low as 31 mm yr we do this for two end member cases in the first case pore pressure is assumed hydrostatic at all depths shallow crustal strength in this case is similar to that calculated in previous work from room temperature friction data while at depths below about 9 13 km depending on slip rate strength becomes less sensitive to depth but sensitive to slip rate in the second case pore pressure is assumed to be near lithostatic at depths below 5 km strength is low at all depths in this case 20 mpa in agreement with observations of weak faults such as the san andreas the predicted depth of transition from velocity weakening to velocity strengthening lies at about 13 km depth for a slip rate of 31 mm yr in rough agreement with the seismic aseismic transition depth observed on mature continental faults these results highlight the importance of fluid assisted deformation processes active in faults at depth and the need for laboratory studies on the roles of additional factors such as fluid chemistry large displacements higher concentrations of phyllosilicates and time dependent fault healing

lava tubes terraces and megatumuli on the 1614 24 pahoehoe lava flow field mount etna sicily the 1614 1624 lava flow of mt etna was formed during a long duration flank eruption involving predominantly pahoehoe flows which produced unusual surface features including mega tumuli here defined and terraces detailed mapping of the flow units surface features and associated tubes reveals a complex sequence of emplacement for the field the stair stepped terraces appear to have been formed as a consequence of self damming of tube fed flows which developed perched ponds of lava surges of lava through tubes elevated sections of crusted lava at the distal ends of the flow to generate tumuli some as high as 130 m as a consequence of pressure via hydrostatic head conditions within the tube although pahoehoe lavas and the related features described here are atypical of mt etna they may reflect styles of eruption and lava emplacement found on volcanoes elsewhere

the reflection seismic survey of project tipteq the inventory of the chilean subduction zone at 38 2 s summary we describe results of an active source seismology experiment across the chilean subduction zone at 38 2 s the seismic sections clearly show the subducted nazca plate with varying reflectivity below the coast the plate interface occurs at 25 km depth as the sharp lower boundary of a 2 5 km thick highly reflective region which we interpret as the subduction channel that is a zone of subducted material with a velocity gradient with respect to the upper and lower plate further downdip along the seismogenic coupling zone the reflectivity decreases in the area of the presumed 1960 valdivia hypocentre the plate interface itself can be traced further down to depths of 50 60 km below the central valley we observe strong reflectivity at the plate interface as well as in the continental mantle wedge the sections also show a segmented forearc crust in the overriding south american plate major features in the accretionary wedge such as the lanalhue fault zone can be identified at the eastern end of the profile a bright west dipping reflector lies perpendicular to the plate interface and may be linked to the volcanic arc

a geophysical model for the origin of volcano vent clusters in a colorado plateau volcanic field variation in spatial density of quaternary volcanic vents and the occurrence of vent clusters correlate with boundaries in proterozoic crust in the springerville volcanic field svf arizona usa inverse modeling using 538 gravity measurements shows that vent clusters correlate with gradients in the gravity field due to lateral variation in crustal density these lateral discontinuities in the crustal density can be explained by boundaries in the north american crust formed during proterozoic accretion spatial density of volcanic vents is low in regions of high density proterozoic crust high in areas of relatively low density proterozoic crust and is greatest adjacent to crustal boundaries vent alignments parallel these boundaries we have developed 2d and 3d numerical models of magma ascent through the crust to simulate long term average magma migration that led to the development of vent clusters in the svf assuming that a viscous fluid flow through a porous media is statistically equivalent to magma migration averaged over geological time in the full field scale the location and flux from the uniform magma source region are boundary conditions of the model changes in model diffusivity associated with changes in the bulk properties of the lithosphere can simulate preferential magma migration paths and alter estimated magma flux at the surface implying that large scale crustal structures such as inherited tectonic block boundaries influence magma ascent and clustering of volcanic vents probabilistic models of volcanic hazard for distributed volcanic fields can be improved by identifying crustal structures and assessing their impact on volcano distribution with the use of numerical models

shape rheology and emplacement times of small martian shield volcanoes abstract this paper focuses on the shape of 31 shield volcanoes of small to intermediate size a few tens of kilometers in diameter in 5 regions of mars tempe terra syria planum pavonis mons arsia mons and central elysium planitia a model for the shape of these shield volcanoes is applied based upon the concept of porous flow of an unconfined aquifer proposed by lacey et al lacey a ockendon j r turcotte d l 1981 on the geometrical form of volcanoes earth planet sci lett 54 1 139 143 and developed by turcotte and schubert turcotte d schubert g 2002 geodynamics cambridge university press second edition the agreement between the topographic profiles and the theoretical shape suggests that the model is appropriate and can be used to derive the product of flow rate and viscosity q estimates of q are found to be homogeneous within a given volcanic region but show large differences from one region to another values generally decreasing with age from oldest to youngest intrinsically the method cannot separate the flow rate from the viscosity but independent observations on individual lava flows in syria planum and central elysium planitia suggest that this evolution may be explained by decreasing magma viscosity with age independent constraints on the viscosity are used to interpret q variations between volcanic regions except for tempe terra and within error bars relative effusion rates are essentially independent of volcano volumes suggesting that the size of a shield volcano is principally controlled by the duration of volcanic activity the specificity of central elysium planitia is emphasized with highest flow rates associated with lowest lava viscosity however the lower viscosities at central elysium planitia cannot totally explain low q values considering the equation of magma flow through a dike from a magma chamber at depth to the shield we suggest that lower volatile contents and or deeper magma chambers may contribute to the decrease of q values with time finally the emplacement times of the 31 shield volcanoes were calculated using a calibration based upon terrestrial volcanoes mauna loa eldborgir and skjalbrei ur emplacement times of martian volcanoes calculated in this way range from a few hundred thousand years tempe terra to a few tens of thousands of years central elysium planitia despite the simplicity of the model the quality of the fits and the internal consistency of the results are reasonable encouraging development of more advanced models linking shield volcano shape to rheological properties of lavas and their characteristic time scales of emplacement

seismogenic zone temperatures and heat flow anomalies in the to nankai margin segment based on temperature data from iodp expedition 333 and thermal model the location of the seismogenic zone across the to nankai margin segment has been widely investigated using geodetic tsunami seismologic and heat flow data as well as thermal models seafloor heat flow measurements show large scatter raising questions about the thermal state of the incoming and overriding plates in this study new temperature and conductivity measurements recorded in the philippine sea plate psp and in the accretionary prism during iodp expedition 333 are integrated with recent data on the margin structure and seismicity to improve constraints on the thermal structure of the incoming plate and the subduction zone as well as the thermal conditions within the seismogenic zone iodp expedition 333 measurements provide heat flow values that are lower than the average of nearby seafloor measurements thermal modeling for the psp suggests that hydrothermal warming is insignificant in the trench thus an observed widespread positive thermal anomaly in the shikoku basin that peaks at the kashinozaki knoll more likely results from vigorous thermal convection in the upper mantle in back arc context the modeled thermal structure of the margin from the trench to the mantle shows a significantly colder interplate contact than in previous studies the seismogenic zone is likely to be 30 170 km from the deformation front corresponding to the 100 340 c temperature range the 1944 co seismic slip zone extends farther southward where temperatures are about 60 c slow slip earthquakes are located in the downdip transition zone between 340 and 440 c across the intersection of the moho and the megathrust fault this improved thermal structure from the trench to the mantle results mainly from a better estimate of oceanic plate age slab dip and sedimentation rate variations in the trench this last parameter cools the interplate contact until a depth of 35 km shifting the thermally defined seismogenic zone landward by 20 km and significantly diminishing its influence near the corner flow area

limited overlap between the seismic gap and coseismic slip of the great 2010 chile earthquake the mw 8 8 chile earthquake on 27 february 2010 occurred in a broad region where tectonic loading has been accumulating since the great 1835 earthquake a joint inversion of tsunami and geodetic data suggests that the 2010 earthquake rupture only partly overlapped with the zone of preseismic locking

lithospheric studies based on array analysis of p coda and microseisms abstract a simple coherence measure the semblance coefficient has been used to sample the lithosphere beneath the norsar array in terms of its ability to produce coherent seismic energy the data set is comprised of two samples of teleseismic p coda waves south of honshu and hindu kush events and two samples of microseismic noise data with duration of about one hour each concerning the constitution of the earth structure there are two possibilities for coherent energy to be present in noise like wavefields the first is tied to scattering effects and the second may be brought about by seismic emission effects the emissive component represents the seismic response of the medium to variations in the background stress field which may arise from various sources passing seismic waves body or surface free oscillations of the earth tides an abnormal level of heat flow or geodynamic processes in general the presence of a well known source of seismic signal hunderfossen dam in the vicinity of norsar area provided an opportunity to test the data processing algorithm results obtained for two independent microseismic data samples point to the same location in close vicinity to the dam for p coda data the signal energy from the upper mantle areas as indicated by the maximum of semblance measure values at the depth around 116 km reaches about 3 8 and 5 of the total p coda energy for honshu and hindu kush events respectively for microseismic noise data this value is about 1 at a depth around 100 km comparison with the results of holographic studies supports the assumption that in the case of microseismic data the detected signal energy is of emissive nature while for p coda data we need more coda samples to make a sound judgement about the nature scattering or emission of the detected signals

the updip and downdip limits to great subduction earthquakes thermal and structural models of cascadia south alaska sw japan and chile this study examines thermal and structural controls of the updip and downdip rupture limits of great subduction thrust earthquakes data on past great earthquake seismic limits have been compiled for four continental subduction zones cascadia sw japan south alaska and chile these limits have then been compared to the predictions of several models for what constrains great earthquake rupture temperatures on the subduction thrusts have been estimated by finite element numerical models the landward limits of the observed updip aseismic zones correspond to the position where the thrust temperature reaches about 100 c that is depths of about 2 to 10 km for the subduction zones studied this temperature agrees with the dehydration of stable sliding smectite clays to illite chlorite the temperatures in this region are controlled mainly by the thickness of sediment on the incoming crust and by the crustal age and thus heat flow the downdip limits correspond to the depth on the thrust where either 1 the temperature reaches about 350 c which corresponds to thermally activated stable sliding behavior for crustal rocks with a transition to 450 c or 2 about 40 km depth if 350 c is reached at greater depth depths of about 40 km approximately correspond to the intersection of the thrust with the continental forearc moho and this downdip limit may be a consequence of stable sliding serpentinite or talc and other hydrated forearc mantle rocks the primary temperature controls on the downdip region are the age of the subducting oceanic plate and the thrust dip profile secondary control comes from the thickness of incoming sediment the convergence rate and the radioactive heat generation in the overlying forearc the 100 c updip limit occurs near the trench for young subducting plates with a thick sediment section such as cascadia 6 8 ma and up to 80 km landward for older oceanic crust such as south alaska 50 ma the 350 c downdip thermal limit is applicable for young oceanic plates e g cascadia and nankai whereas the forearc mantle limit applies for older plates e g south alaska and chile except near the chile rise for the margins studied that have experienced great earthquakes there is generally good agreement between the postulated thermal and moho limits and the rupture or seismogenic zone as defined by the distribution of aftershocks and by waveform tsunami and dislocation modeling the downdip limit of the interseismic locked zone from dislocation modeling is also in agreement with these limits

global analysis of the effect of fluid flow on subduction zone temperatures 1 knowledge of the controls on temperature distributions at subduction zones is critical for understanding a wide range of seismic metamorphic and magmatic processes here we present the results of 220 thermal model simulations covering the majority of known subduction zone convergence rates incoming plate ages and slab dips we quantify the thermal effects of fluid circulation in the subducting crust by comparing results with and without advective heat transfer in the oceanic crustal aquifer we find that hydrothermal cooling of a subduction zone is maximized when the subducting slab is young slowly converging steeply dipping and the crustal aquifer is ventilated near the trench incoming plate age is one of the primary controls on the effectiveness of advective heat transfer in the aquifer and the greatest temperature effects occur with an incoming plate 50 ma the thermal effects of fluid circulation decrease dramatically with increasing age of the incoming plate temperatures in the cascadia nankai southern chile colombia ecuador mexico and solomon islands subduction zones are likely strongly affected by fluid circulation for these systems only thermal models of cascadia and nankai have included fluid flow in subducting crust

crustal thickness and mantle wedge structure from receiver functions in the chilean maule region at 35 s a temporary passive seismic network of 21 broad band stations was deployed in central chile between 35 s to 36 s the network recorded data prior to the magnitude mw 8 8 2010 maule earthquake at a latitude of the major slip and surface deformation the experiment was conducted to survey crustal and mantle structures and to assess the state of hydration of the mantle wedge we present results of a teleseismic p receiver function study supporting a continental moho at approx 38 km depth phase conversion at this boundary could be observed continuously from the intersection of the subducting slab with the continental moho towards the andes the character of receiver functions indicated little evidence for infiltration of water from the subducting plate into the overlying mantle wedge suggesting that only a small amount of water is released from the subducting plate aftershocks of the maule earthquake and post seismic slip reached depths of 50 km and hence slip spread down dip of the continental moho in the post seismic phase co seismic rupture however occurred updip of the continental moho spare aftershock seismicity is observed at the intersection of the continental moho with the subducting slab

downdip landward limit of cascadia great earthquake rupture 1 this paper examines the constraints to the downdip landward limit of rupture for the cascadia great earthquakes off western north america this limit is a primary control for ground motion hazard at near coastal cities the studies also provide information on the physical controls of subduction thrust rupture globally the constraints are 1 locked transition zones from geodetic deformation gps repeated leveling tide gauges 2 rupture zone from paleoseismic coastal marsh subsidence paleogeodesy 3 temperature on the thrust for the seismic aseismic transition 4 change in thrust seismic reflection character downdip from thin seismic to thick ductile 5 fore arc mantle corner aseismic serpentinite and talc overlying the thrust 6 updip limit of episodic tremor and slip ets slow slip 7 rupture area associations with shelf slope basins 8 depth limit for small events on the thrust and 9 landward limit of earthquakes on the nootka transform fault zone the most reliable constraints for the limit of large rupture displacement 10 m are generally just offshore in agreement with thermal control for this hot subduction zone but well offshore central oregon and near the coast of northern washington the limit for 1 2 m rupture that can still provide strong shaking is less well estimated 25 50 km farther landward the fore arc mantle corner and the updip extent of ets slow slip are significantly landward from the other constraints surprisingly there is a downdip gap between the best other estimates for the great earthquake rupture zone and the ets slow slip in this gap plate convergence may occur as continuous slow creep

understanding caldera structure and development an overview of analogue models compared to natural calderas abstract understanding the structure and development of calderas is crucial for predicting their behaviour during periods of unrest and to plan geothermal and ore exploitation geological data including that from analysis of deeply eroded examples allow the overall surface setting of calderas to be defined whereas deep drillings and geophysical investigations provide insights on their subsurface structure collation of this information from calderas worldwide has resulted in the recent literature in five main caldera types downsag piston funnel piecemeal trapdoor being viewed as end members despite its importance such a classification does not adequately examine a the structure of calderas particularly the nature of the caldera s bounding faults and b how this is achieved including the genetic relationships among the five caldera types various sets of analogue models specifically devoted to study caldera architecture and development have been recently performed under different conditions apparatus materials scaling parameters stress conditions the first part of this study reviews these experiments which induce collapse as a result of underpressure or overpressure within the chamber analogue the experiments simulating overpressure display consistent results but the experimental depressions require an exceptional amount of doming seldom observed in nature to form therefore these experiments are not appropriate to understand the structure and formation of most natural calderas the experiments simulating underpressure reveal a consistent scenario for caldera structure and development regardless of their different boundary conditions these show that complete collapse proceeds through four main stages proportional to the amount of subsidence progressively characterized by 1 downsag 2 reverse ring fault 3 peripheral downsag 4 peripheral normal ring fault the second part of this study verifies the possibility that these latter calderas constitute a suitable analogue to nature and consists of a comprehensive comparison of the underpressure experiments to natural calderas this shows that all the experimental structures as well as their progressive development are commonly observed at natural calderas highlighting a consistency between models and nature as the shallow structure of experimental calderas corresponds to a precise architecture at depth it provides a unique key to infer the deeper structure of natural calderas recognizing diagnostic surface features within a caldera will thus allow it to be categorized within a precise structural and evolutionary context the general relationship between the evolutionary stage of a caldera and its d s diameter subsidence ratio allows such a quantification with stage 1 calderas characterized by d s 40 stage 2 by 18 the four stages adequately explain the architecture and development of the established caldera end members along a continuum where one or more end members downsag piston funnel piecemeal trapdoor may correspond to a specific stage while such a continuum is controlled by progressive subsidence specific collapse geometries will result from secondary contributory factors roof aspect ratio collapse symmetry pre existing faults these considerations allow proposing an original classification of calderas incorporating their structural and genetic features

episodic development of a convergent margin inferred from structures and processes along the southern chile margin seismic reflection data acquired in the vicinity of isla mocha across the southern coast of chile image structures formed along the continental margin and reveal an episodic history of accretion nonaccretion and possibly erosion structures formed at the toe of the continental slope suggest frontal accretion of to 1 km of trench fill seismic images also reveal that a small accretionary wedge 20 30 km wide abuts the truncated continental metamorphic basement that extends seaward from beneath the shelf the small size of the accretionary wedge on three profiles examined here is not consistent with a long history of accretion with the current deformational style as current rates of frontal accretion could have accumulated all of the existing accretionary wedge in less than 1 2 m y this is a small fraction of convergence history along this margin and the current accretionary mode has not been consistently maintained in the past the isla mocha region is located between the temperate climate of central chile and the glacial climate of southern chile and climatic conditions in this region have likely fluctuated sufficiently to cause significant variation in trench sediment supply accretionary and nonaccretionary or erosional episodes are probably linked to temporal variations in trench sediment thickness as suggested by observations along the chile margin currently thick trench sediment correlates with accretion along the southern chile margin and thin trench sediment correlates with nonaccretion tectonic erosion as near the chile ridge and from the juan fernandez ridge to northern chile the isla mocha region also lies 900 1000 km north of the chile triple junction and the chile ridge lies approximately 2000 km to the west and has not yet collided and affected the margin near isla mocha this part of the precollision zone provides an excellent reference to examine the effects of chile ridge collision in the development of the chile margin the most apparent effect of subduction of the buoyant young crust of the chile ridge is a shallow trench that is nearly devoid of sediment consequently the triple junction is undergoing nonaccretion or erosion and the accretionary complex near the triple junction remains smaller than to the north or south because the current phase of rapid accretion elsewhere in the trench has bypassed the triple junction region the interplay of subduction zone processes such as trench sedimentation and ridge collision has resulted in an episodic development of the margin and produced a discontinuous record of convergence history within the accretionary wedge

calderas and caldera structures a review abstract calderas are important features in all volcanic environments and are commonly the sites of geothermal activity and mineralisation yet it is only in the last 25 years that a thorough three dimensional study of calderas has been carried out utilising studies of eroded calderas geophysical analysis of their structures and analogue modelling of caldera formation as more data has become available on calderas their individuality has become apparent a distinction between caldera caldera complex cauldron and ring structure is necessary and new definitions are given in this paper descriptions of calderas based on dominant composition of eruptives basaltic peralkaline andesitic dacitic rhyolitic can be used and characteristics of each broad group are given styles of eruption may be effusive or explosive with the former dominant in basaltic calderas and the latter dominant in andesitic dacitic rhyolitic and peralkaline calderas four end member collapse styles occur plate or piston piecemeal trapdoor and downsag but many calderas have multiple styles features of so called funnel and chaotic calderas proposed in the literature can be explained by other collapse styles and the terms are considered unnecessary ground deformation comprises subsidence or collapse essential characteristics of a caldera and uplifting doming and fracturing due to tumescence and or resurgence frequent but not essential collapse may occur on pre existing structures such as regional faults or on faults created during the caldera formation and the shape of the collapse area will be influenced by depth size and shape of the magma chamber the final morphology of a caldera will depend on how the caldera floor breaks up whether collapse takes place in one event or multiple events whether vertical movement is spasmodic or continuous throughout the eruptive sequence and whether blocks subside uniformly or chaotically at one or more collapse centres a meaningful description of any caldera should therefore include number of collapse events presence or absence of resurgence caldera floor coherency caldera floor collapse geometry and dominant composition of eruptives

adaptive finite element modelling of two dimensional magnetotelluric fields in general anisotropic media summary an adaptive unstructured mesh finite element fe procedure is presented for improving the quality of numerical solutions to the magnetotelluric forward problem in a general 2 d anisotropic conductivity structure we implement a self adaptive goal oriented grid refinement algorithm in which a finite element analysis is performed on a sequence of refined meshes the mesh refinement process is guided by a dual error estimate weighting to bias refinement towards elements that affect the solution at the em receiver locations we validate the finite element code against a layered 1 d model with a sea water layer further we compare the fe results with those obtained by a finite difference fd scheme for both a block seamountain and a sea bottom hill model both fe and fd schemes show very good agreement for the block seamountain model for the sea bottom hill model however only on the flat seafloor segments both the fe and fd solutions fit very well but on the seafloor slope fd results are oscillating due to a simplistic staircase approximation of the bathymetric undulations the fd scheme for 2 d anisotropic conductors developed primarily for the modelling of magnetotelluric data on a flat earth surface is thus not an adequate tool for dealing with structures with sloping bathymetry and topography whereas the fe method with adapting mesh can easily handle such structures at almost any level of complexity

alteration of the subducting oceanic lithosphere at the southern central chile trench outer rise hydrothermal circulation and brittle faulting processes affecting the oceanic lithosphere are usually confined to the upper crust for oceanic lithosphere created at intermediate to fast spreading rates lower crust and mantle rocks are therefore relatively dry and undeformed however recent studies at subduction zones suggest that hydration of the oceanic plate is most vigorous at the trench outer rise where extensional bending related faulting affects the hydrogeology of the oceanic crust and mantle to understand the degree of hydration we studied the seismic velocity structure of the incoming nazca plate offshore of southern central chile 43 s here the deep sea trench is heavily filled with up to 2 km of sediments seismic refraction and wide angle data complemented by seismic reflection imaging of sediments are used to derive a two dimensional velocity model using joint refraction and reflection traveltime tomography the seismic profile runs perpendicular to the spreading ridge and trench axes the velocity model derived from the tomography inversion consists of a 5 3 km thick oceanic crust and shows p wave velocities typical for mature fast spreading crust in the seaward section of the profile with uppermost mantle velocities as fast as 8 3 km s approaching the chile trench seismic velocities are significantly reduced however suggesting that the structures of both the oceanic crust and uppermost mantle have been altered possibly due to a certain degree of fracturing and hydration the decrease of the velocities roughly starts at the outer rise 120 km from the deformation front and continues into the trench even though the trench is filled with sediment basement outcrops in the outer rise frequently pierce the sedimentary blanket anomalously low heat flow values near outcropping basement highs indicate an efficient inflow of cold seawater into the oceanic crust hydration and crustal cracks activated by extensional bending related faulting are suggested to govern the reduced velocities in the vicinity of the trench considering typical flow distances of 50 km water might be redistributed over most of the trench outer rise area where trapped in faults seawater may migrate down to mantle depth causing up to 9 of serpentinization in at least the uppermost 2 km of the mantle between the outer rise and the trench axis

a porous flow model of magma migration within mt etna the influence of extended sources and permeability anisotropy abstract the isotropic point source model of porous flow proposed by bonafede and boschi is employed to model magma flow within the edifice of mt etna in probabilistic terms the flow pattern is found to be nearly horizontal suggesting magma propagation through lateral dykes departing radially from the central conduits the point source model is able to provide a detailed description of eruption site distribution for a few sections of mt etna but it fails to explain simultaneously the eruption distribution along the mutually orthogonal n s and e w sections the solution is then generalized to deal with extended sources and with anisotropic media employed respectively to simulate a hypothetical n s trending fracture system at the base of the volcano or the influence of the regional stress field on the directional probability of dyke opening the results show that both models are able to explain the presence on mt etna of the n s trending rift in particular the anisotropic model provides the best fits to the north south west and east topographic sections several features of the eruption distribution on mt etna remain unexplained however by a porous flow model and apparently require the intervention of specific structural effects for instance a detailed reproduction of the three lobe distribution of eruptive sites observed in recent historic times seems to require permeability heterogeneities the effect of local topographic loading or the presence of a local axial stress field related to the inflation of mt etna

importance of gravitational spreading in the tectonic and volcanic evolution of mount etna the interaction between gravity and thermal effects largely determines the structural and magmatic evolution of volcanic constructs at scales spanning several orders of magnitude from small cones to the oceanic crust1 3 although gravitational spreading is a direct consequence of this interaction and a fundamental process in volcano growth it is rarely recognized as such here we describe a striking example of this process at etna volcano in italy the volcano and its clay rich substratum are slowly spreading towards the east and south driven by gravity spreading produces extensional structures in the summit region and contpressional structures at the base of the volcano eastward movement of the volcanic edifice over a stationary magma supply may be the cause of an apparent westward migration of volcanic activity as gravitational spreading seems to control the location and magnitude of shallow seismicity and flank eruptions an appreciation of its effects could become an essential element of future volcanic hazard evaluation our model proposed here for etna may also be relevant to a reinterpretation of the geological history of a number of other volcanoes

case for very low coupling stress on the cascadia ssubduction fault a fundamental problem in plate tectonics is the shear strength of major plate boundary faults this translates to the question whether the generally observed small earthquake stress drops of 3 10 mpa on major faults release most of the accumulated stress or only a small fraction of it there is strong evidence that the san andreas fault a major transform plate boundary is weak 20 mpa shear resistance it is not yet clear whether subduction thrust faults are also weak we present two types of evidence from the northern cascadia subduction zone that indicate very low coupling shear stress on that plate interface and hence very low strength of the subduction thrust fault comparable to that estimated for the san andreas fault first the well defined surface heat flow and heat generation allow negligible frictional heating on the plate interface the average shear stress on the fault must thus be very low over a time scale of a few million years second focal mechanism solutions for small crustal earthquakes in the southern vancouver island area indicate that the horizontal stress in the direction of plate convergence has a similar magnitude to the vertical stress this inferred stress state requires the present tectonic stress coupled across the subduction thrust fault to be very low one explanation for the weakness of the fault is the presence of near lithostatic pore fluid pressure in the region of the fault zone for which there is independent evidence the conclusion of a weak subduction thrust fault does not conflict with geodetic observations of contemporary surface deformation which indicate that the fault is currently locked accumulating strain energy toward a future great earthquake the surface deformation responds to the small 20 mpa temporal changes of the stress field associated with the subduction earthquake cycle this transient stress is superimposed on the larger background regional stress field in which the maximum compression is parallel to the margin the weakness of the cascadia subduction thrust fault and the unusual stress state of the forearc region have important implications for earthquake hazards for example a subduction earthquake may induce large strike slip earthquakes in the forearc that affect a large area

crustal motion in the zone of the 1960 chile earthquake detangling earthquake cycle deformation and forearc sliver translation temporary deformation in great earthquake cycles and permanent shear deformation associated with oblique plate convergence both provide critical clues for understanding geodynamics and earthquake hazard at subduction zones in the region affected by the mw 9 5 great chile earthquake of 1960 we have obtained gps observations that provide information on both types of deformation our velocity solutions for the first time span the entire latitudinal range of the 1960 earthquake the new observations revealed a pattern of opposing roughly arc normal motion of coastal and inland sites consistent with what was reported earlier for the northern part of this region this finding supports the model of prolonged postseismic deformation as a result of viscoelastic stress relaxation in the mantle the new observations also provide the first geodetic evidence for the dextral motion of an intravolcanic arc fault system and the consequent northward translation of a forearc sliver the sliver motion can be modeled using a rate of 6 5 mm a accommodating about 75 of the margin parallel component of nazca south america relative plate motion with the rate diminishing to the north furthermore the new gps observations show a southward decrease in margin normal velocities of the coastal area we prefer explaining the southward decrease in terms of changes in the width or frictional properties of the megathrust seismogenic zone because of the much younger age of the subducting plate and warmer thermal regime in the south the currently locked portion of the plate interface may be narrower using a three dimensional viscoelastic finite element model of postseismic and interseismic deformation following the 1960 earthquake we demonstrate that this explanation although not unique is consistent with the gps observations to the first order

three dimensional density model of the nazca plate and the andean continental margin 1 we forward modeled the bouguer anomaly in a region encompassing the pacific ocean 85 w and the andean margin 60 w between northern peru 5 s and patagonia 45 s the three dimensional density model that reproduces the gravity field is a continental scale representation of density structure to 410 km depth that characterizes the mantle and crust of the oceanic nazca plate subducted slab and continental margin with a minimum number of bodies we predefined the density of each body after studying the dependency of density on composition of crustal and mantle materials and pressure temperature conditions appropriate for the andean setting a database of independent geophysical information constrains the geometry of the top of the subducted slab locally the moho of the oceanic and continental crusts and indirectly the lithosphere asthenosphere boundary underneath the continental plate other boundaries notably the intracrustal density discontinuity separating upper from lower crust below the continent were not constrained and their geometry is the result of fitting the observed and calculated bouguer anomaly during forward modeling this contribution presents the model to the andean geoscientific community and contains some tools like a sensitivity analysis that helps potential users of the model to interpret its results we describe and discuss some of these results in order to illustrate the application of the model to the study of a wide range of phenomena e g modification of oceanic plate structure by hot spots shape of the subducted slab thermal structure of the continental lithosphere compensation mechanism and formation of orogenic relieve causes of andean segmentation

water input and water release from the subducting nazca plate along southern central chile 33 s 46 s the age of the subducting nazca plate off chile increases northward from 0 ma at the chile triple junction 46 s to 37 ma at the latitude of valparaiso 32 s age related variations in the thermal state of the subducting plate impact on a the water influx to the subduction zone as well as on b the volumes of water that are released under the continental fore arc or alternatively carried beyond the arc southern central chile is an ideal setting to study this effect because other factors for the subduction zone water budget appear constant we determine the water influx by calculating the crustal water uptake and by modeling the upper mantle serpentinization at the outer rise of the chile trench the water release under fore arc and arc is determined by coupling fem thermal models of the subducting plate with stability fields of water releasing mineral reactions for upper and lower crust and hydrated mantle results show that both the influx of water stored in and the outflux of water released from upper crust lower crust and mantle vary drastically over segment boundaries in particular the oldest and coldest segments carry roughly twice as much water into the subduction zone as the youngest and hottest segments but their release flux to the fore arc is only about one fourth of the latter this high variability over a subduction zone of 1500 km length shows that it is insufficient to consider subduction zones as uniform entities in global estimates of subduction zone fluxes

coseismic strengthening of the shallow portion of the subduction fault and its effects on wedge taper 1 according to the dynamic coulomb wedge model coseismic strengthening of the shallowest part of the subduction interface can cause permanent deformation of the overlying outer wedge so that the cumulative effects of great interplate earthquakes control the geometry of the wedge in this work we use a numerical model which is a hybrid of the frictional contact model and the classical crack model to study how stress is coseismically transferred from the seismogenic zone to the strengthening updip zone to cause wedge deformation in this static model the critical strengthening required to prevent the rupture from breaking the trench depends on the force drop of the seismogenic zone defined as the product of the average shear stress drop and the area of the seismogenic zone in a simple model of uniform material properties with a few megapascals average stress drop over a seismogenic zone of 120 km downdip width the critical strengthening for a 30 km wide updip zone is an increase in the effective friction coefficient by about 0 05 using the coulomb wedge theory we demonstrate that this level of strengthening can readily push the overlying wedge into a critical state of failure with much greater strengthening the rupture is able to extend into the updip zone only slightly causing localized wedge compression in the area of slip termination we examined wedge geometry of 22 subduction zones in the light of the model results we found that the surface slope of these wedges is generally too high to be explained using the classical wedge theory but can be explained using the dynamic coulomb wedge model including coseismic strengthening of the shallow portion of the megathrust

long distance fluid and heat transport in the oceanic crust entering the nankai subduction zone nantroseize transect abstract i examine the potential causes of anomalous seafloor heat flux on the oceanic plate in the nankai trough seismogenic zone experiment nantroseize transect offshore southern japan the most prominent anomaly is a 50 mw m 2 change in heat flux between integrated ocean drilling program sites c0011 and c0012 over a distance of 30 km wide transition at the seafloor the observed surface heat flux pattern is indicative of hydrothermal circulation in the basement aquifer and advection of heat from the subducted crust into the aquifer on the incoming plate for a 600 m thick aquifer the permeability is likely 7 10 11 m 2 and hydrothermal circulation transports at least 300 times more heat than conduction alone the heat flux from the subduction zone seaward to the incoming plate is consistent with hydrothermal circulation in the subducting crust persisting to 100 km landward of the deformation front vigorous fluid circulation in the basaltic basement is consistent with both the seafloor thermal anomalies and geochemical anomalies near the sediment basement interface

thermal regime of the costa rican convergent margin 2 thermal models of the shallow middle america subduction zone offshore costa rica at the costa rica margin along the middle america trench along strike variations in heat flow are well mapped these variations can be understood in terms of either ventilated fluid flow where exposed basement allows fluids to freely advect heat between the crustal aquifer and ocean or insulated fluid flow where continuous sediment cover restricts heat advection to within the crustal aquifer we model fluid flow within the subducting aquifer using nusselt number approximations coupled with finite element models of subduction and explore its effect on temperatures along the subduction thrust the sensitivity of these models to the initial thermal state of the plate and styles of fluid flow either ventilated or insulated is explored heat flow measurements on cool crust accreted at the east pacific rise are consistent with ventilated hydrothermal cooling that continues with subduction these models yield much cooler temperatures than predicted from simulations initialized with conductive predictions and without hydrothermal circulation heat flow transects on warm crust accreted at the cocos nazca spreading center are consistent with models of insulated hydrothermal circulation that advects heat updip within the subducting crustal aquifer near the trench these models are warmer than conductive predictions and cooler than conductive predictions downdip of the trench comparisons between microseismicity and modeled isotherms suggest that the updip limit of microseismicity occurs at temperatures warmer than 100 c and that the downdip extent of microseismicity is bounded by the intersection of the subduction thrust with the base of the overriding crust

a wider seismogenic zone at cascadia due to fluid circulation in subducting oceanic crust in the cascadia subduction zone the extent of the seismogenic portion of the plate interface is poorly resolved by seismicity due to the lack of a large megathrust event during the instrumental record therefore fault zone temperatures 150 to 350 c are used to estimate the limits of the seismogenic zone previous thermal models for the cascadia margin estimated that 350 c on the plate boundary occurs 40 70 km offshore in contrast models of interseismic deformation have been interpreted to indicate a seismogenic zone extending landward of the coastline to the updip edge of a region of episodic tremor and slip ets we examined cascadia subduction zone temperatures with thermal models that include the effects of fluid circulation in an ocean crust aquifer fluid circulation cools the subduction zone and widens the thermally defined seismogenic zone by shifting the 350 c isotherm at the plate boundary 30 55 km landward relative to results from simulations without fluid flow our thermal models indicate a 60 80 km wide transition zone between 350 c on the fault and the updip edge of ets under british columbia canada washington and oregon united states ets occurs at 410 550 c the location of the basalt to eclogite transition in the subducting crust provides an important constraint on the thermal models because hydrothermal circulation in the ocean crust aquifer produces only small surface heat flux anomalies on this margin with thick sediment in the trench and on the incoming plate

a porous flow model for the geometrical form of volcanoes critical comments a critical evaluation is presented of the assumptions on which the mathematical model for the geometrical form of a volcano arising from the flow of magma in a porous medium of lacey et al 1981 is based the lack of evidence for an equipotential surface or its equivalent in volcanoes prior to eruption is pointed out and the preference of volcanic eruptions for low ground is attributed to the local stress field produced by topographic loading rather than a rising magma table other difficulties with the model involve the neglect of the surface flow of lava under gravity away from the vent and the use of the dupuit approximation for unconfined flow and the assumption of essentially horizontal magma flow comparisons of model predictions with the shapes of actual volcanoes reveal the model not to fit lava shield volcanoes for which the cone represents the solidification of small lava flows and to provide a poor fit to composite central volcanoes

thermal models of the mexico subduction zone implications for the megathrust seismogenic zone 1 it has been proposed that the seismogenic zone of subduction thrust faults is primarily controlled by temperature or rock composition changes we have developed numerical models of the thermal structure of the mexico subduction zone to examine the factors that affect the temperature of the subduction thrust fault although the oceanic plates subducting beneath mexico are young the top of the oceanic plate at the trench is cool because of the lack of a thick cover of insulating sediments marine heat flow observations suggest that hydrothermal circulation may further cool the oceanic plate this results in a cool subduction thrust fault where the brittle part of the fault extends to depths of over 40 km at these depths even slight frictional heating may have significant effects on temperature along the thrust fault particularly for regions with a high convergence rate and shallow plate dip with the addition of a small amount of frictional heating the temperatures of the deep 30 40 km thrust fault are increased by over 200 c as the observed downdip limit of rupture in recent well constrained megathrust earthquakes is confined to depths above the intersection of the thrust fault and the continental moho a temperature of 350 c may control the downdip extent of the seismogenic zone thus in order to be consistent with the observed shallow rupture areas it is necessary to include a small amount of frictional heating corresponding to an average shear stress of 15 mpa

caldera structure amount of collapse and erupted volumes the case of bolsena caldera italy calderas are common on volcanoes but their structure is seldom visible the 19 km wide bolsena caldera italy formed between 0 6 and 0 2 ma the largely preserved structural rim and subsurface data make bolsena ideal to investigate caldera structure in relation to the subsidence and erupted volume in this paper we use remote sensing field analysis and available subsurface data at the surface the caldera passes from a downsag south rim to a narrow and densely faulted area north rim with outer normal and inner reverse faults the caldera structure on the widely faulted east rim appears to be scale dependent with a staircase like fault zone larger scale horst and graben like structures intermediate scale and domino like structures smaller scale subsurface data indicate asymmetric collapse with a northward increase in subsidence ranging from diffuse to the south to focused to the north deformation at the surface the collapse rate constant between ca 490 and 175 ka was more than magma output between ca 330 and 130 ka highlighting significant 200 m and prolonged 200 ka posteruptive subsidence as the nearby latera caldera west rim of bolsena was mostly active between ca 265 and 160 ka much of the subsidence at bolsena may be related to this activity suggesting a common magmatic reservoir the subsidence related structural variations along the caldera rim and the significant posteruptive subsidence found at bolsena have not been found in other calderas

regional heat flow variations across the sedimented juan de fuca ridge eastern flank constraints on lithospheric cooling and lateral hydrothermal heat transport seafloor heat flow has been estimated continuously along a transect on the sedimented eastern flank of the juan de fuca ridge to provide constraints on the scale and rate of lateral heat and fluid transport in the upper igneous crust and on the total flux from the young lithosphere beneath the profile extends from 20 km east of the ridge axis where turbidite sediments lap onto the oceanic crust at an age of less than 1 ma and ends 110 km east of the axis over 3 6 ma crust estimated heat flow is derived from the combination of the depth to igneous basement determined from seismic reflection profiles and basement temperatures measured in nine ocean drilling program boreholes along the transect the strategy is based on the knowledge that in this young area upper basement temperatures are locally homogenized by vigorous hydrothermal circulation and hence vary smoothly between drilling sites at the young end of the transect near the region of extensive outcrop heat flow is only 15 of that predicted for lithosphere of this age heat flow increases systematically to the east reaching the full predicted level about 20 km from the point of sediment basement onlap the anomaly can be accounted for by mixing and or eastward flow of cool water in the upper igneous crust beneath the sediment at a darcian rate of the order of 1 m yr a geochemical signature of fresh seawater in basement is observed many tens of kilometers farther to the east although the thermal signature is no longer present along the remainder of the transect the heat flow conforms very well to that predicted by simple boundary layer cooling theory once the thermal effects of rapid sedimentation on this ridge flank are accounted for

the campanian plain and phlegrean fields structural setting from potential field data abstract a boundary analysis of the gravity and magnetic fields of the phlegrean fields volcanic area and of the surrounding campanian plain reveals a complex structural setting the campanian plain results well defined from sharp density boundaries a set of major e w lineaments occurs within the plain we derive a framework of the structural control to the volcanological evolution of campanian plain which appears strongly influenced by a set of ne sw lineaments bordering the acerra depression the density and magnetization boundaries of the phlegrean fields are significantly consistent both indicating for the collapsed phlegrean caldera an area considerably less extended than that based on previous geologically defined models many other structures are evidenced in the phlegrean area often with agreement between density and magnetization boundaries inside the caldera a strong consistence between the structures we identified and recent seismicity has been noted

the bay of naples southern italy constraints on the volcanic structures inferred from a dense seismic survey 1 the bays of naples and pozzuoli host volcanic areas that have been active within recent history in addition the vicinities areas around the vesuvius volcano and the campi flegrei caldera are extremely densely inhabited over the past decade this situation has prompted several active seismic experiments focused on the vesuvius volcano more recently the serapis active seismic survey covered the entire bay of naples with particular attention to pozzuoli bay the processed serapis p wave travel time collection extended with a previously acquired data set has allowed the computation of the most comprehensive three dimensional p wave velocity distribution for the bay of naples which includes a small scale high resolution model of the campi flegrei caldera the joint interpretation of the velocity distribution together with both the available gravity measurements and models and the geothermal drilling information has allowed us to locate and map lithological units at depths we propose a model of the campi flegrei caldera in relation to the underlying mesozoic limestone unit in this model no magma chamber exists in the upper 6 km beneath the bay of naples the feeding system of the campi flegrei caldera is directly controlled by the tectonic structures affecting the limestone unit beneath the bay of naples two ne sw normal faults run beneath the vesuvius volcano and campi flegrei which are related to the pre pleistocene apenninic tectonics

slab1 0 a three dimensional model of global subduction zone geometries 1 we describe and present a new model of global subduction zone geometries called slab1 0 an extension of previous efforts to constrain the two dimensional non planar geometry of subduction zones around the focus of large earthquakes slab1 0 describes the detailed non planar three dimensional geometry of approximately 85 of subduction zones worldwide while the model focuses on the detailed form of each slab from their trenches through the seismogenic zone where it combines data sets from active source and passive seismology it also continues to the limits of their seismic extent in the upper mid mantle providing a uniform approach to the definition of the entire seismically active slab geometry examples are shown for two well constrained global locations models for many other regions are available and can be freely downloaded in several formats from our new slab1 0 website http on doi gov d9arbs we describe improvements in our two dimensional geometry constraint inversion including the use of average active source seismic data profiles in the shallow trench regions where data are otherwise lacking derived from the interpolation between other active source seismic data along strike in the same subduction zone we include several analyses of the uncertainty and robustness of our three dimensional interpolation methods in addition we use the filtered subduction related earthquake data sets compiled to build slab1 0 in a reassessment of previous analyses of the deep limit of the thrust interface seismogenic zone for all subduction zones included in our global model thus far concluding that the width of these seismogenic zones is on average 30 larger than previous studies have suggested

structural analysis and thermal remote sensing of the los humeros volcanic complex implications for volcano structure and geothermal exploration abstract the los humeros volcanic complex lhvc is an important geothermal target in the trans mexican volcanic belt understanding the structure of the lhvc and its influence on the occurrence of thermal anomalies and hydrothermal fluids is important to get insights into the interplay between the volcano tectonic setting and the characteristics of the geothermal resources in the area in this study we present a structural analysis of the lhvc focused on quaternary tectonic and volcano tectonic features including the areal distribution of monogenetic volcanic centers morphostructural analysis and structural field mapping revealed the geometry kinematics and dynamics of the structural features in the study area also thermal infrared remote sensing analysis has been applied to the lhvc for the first time to map the main endogenous thermal anomalies these data are integrated with newly proposed unconformity bounded stratigraphic units to evaluate the implications for the structural behavior of the caldera complex and geothermal field the lhvc is characterized by a multistage formation with at least two major episodes of caldera collapse los humeros caldera 460 ka and los potreros caldera 100 ka the study suggests that the geometry of the first collapse recalls a trap door structure and impinges on a thick volcanic succession 10 5 1 55 ma now hosting the geothermal reservoir the main ring faults of the two calderas are buried and sealed by the widespread post calderas volcanic products and for this reason they probably do not have enough permeability to be the main conveyers of the hydrothermal fluid circulation an active previously unrecognized fault system of volcano tectonic origin has been identified inside the los potreros caldera this fault system is the main geothermal target probably originated by active resurgence of the caldera floor the active fault system defines three distinct structural sectors in the caldera floor where the occurrence of hydrothermal fluids is controlled by fault induced secondary permeability the resurgence of the caldera floor could be induced by an inferred magmatic intrusion representing the heat source of the geothermal system and feeding the simultaneous monogenetic volcanic activity around the deforming area the operation of the geothermal field and the plans for further exploration should focus on both the active resurgence fault system and the new endogenous thermal anomalies mapped outside the known boundaries of the geothermal field

different styles of deformation in the calabrian arc southern italy implications for a seismotectonic zoning abstract the most recent seismic profiles in southern italy show the existence of a sudden thinning of the crust at the boundary between the chain and the tyrrhenian margin the abrupt change in thickness can be followed along a zone which has the same geometry in plan as the calabrian arc by considering these data and the gravimetric anomalies it is possible to associate this crustal anomaly with a deep seated shear zone which determines a crustal shortening of about 40 60 km which can only be a consequence of post tortonian tectonics the surface equivalents of the deep seated shear zone can be recognized from north to south in the alignment of the intra apenninic basins vallo di diano san arcangelo potenza in the crati mesima graben and in the mount kumeta alcantara fault zone the distribution of the seismicity and its connection with the surface structures shows that the largest earthquakes occur along the deep shear zone in particular shocks with the highest magnitude and the longest recurrence intervals are located in those areas where the deep shear zone is at an angle of about 90 with the direction of maximum shortening crati mesima graben the fault zones nearly parallel to the regional compression axis e g mount kumeta alcantara fault zone are characterized by earthquakes of lower magnitude taking into account the neotectonic evolution of the regional structures as well as the orientation of the stress field and its connection with the deep seated shear zone it is possible to distinguish the following seismotectonic zones upper crati mesima graben transverse throughs mount pollino mount raparo fault zone external ionian area san arcangelo basin zone of the external flysch sicani mountains mount kumeta alcantara fault zone caltanissetta basin and iblean plateau bradanic trough murge ridge the definition of the geometry of the shear zone at depth is one of the most important but still unresolved problems

history of the exploitation of thermo mineral resources in campi flegrei and ischia italy abstract thermo mineral springs are widely spread over the volcanic areas of campi flegrei and ischia and on the western side of the city of naples their use has since long been known by local inhabitants because they were one of the main curative therapies in the past the fortune and use of the thermal springs have been affected by the volcanic eruptions that occurred in different times in these areas the exploitation of geothermal springs in campi flegrei suffered a decisive decline after the eruption of monte nuovo in 1538 conversely after that eruption in ischia island there was a progressive increase of the use of thermal spring therapy the use of geothermal resources for electric generation has not been implemented in campi flegrei because of the difficulty in using hypersaline fluids whereas in ischia local authorities preferred only the development of touristic activity based on geothermal springs

3d scattering image of the campi flegrei caldera southern italy new hints on the position of the old caldera rim abstract a scattering image of the campi flegrei caldera has been obtained using coda wave envelopes from local microearthquake data recorded at a portable 12 station digital network this was operated during the paroxysmic phase of the ground uplift bradyseism that occurred in the years 1983 1984 the nishigami technique has been used to obtain the 3d distribution of the strong scatterers in an earth volume with a horizontal extension of 20 km 20 km and a vertical thickness of 6 km this earth volume was segmented in blocks of 2 km 2 km 2 km for the outer section and in blocks of 1 km 1 km 1 km for the inner section the inversion was performed using a hybrid jackknife svd method that has allowed qualitative control of the stability and robustness of the solution the results show that the scatterers revealed by the coda envelopes which were band pass filtered at a center frequency of 18 hz match the inner border of the caldera rim which is well outlined on land by geological observations the scatterer position also delineates this border well beneath sea level these results are in close agreement with recently obtained seismic velocity tomography the scatterers imaged at the 6 and 12 hz band pass frequencies are located around the zone of maximum uplift beneath solfatara crater where most of the present residual volcanic activity is concentrated

interaction between the chile ridge and chile trench geophysical and geothermal evidence geophysical and geothermal data are examined from the three southernmost sections of the chile ridge starting at 44 s and continuing south to the triple junction of the nazca antarctic and south america plates at 47 s these sections represent three progressively younger stages in a ridge trench collision event corresponding to 3 m y before the collision 1 m y before the collision and culminating in an ongoing collision at the triple junction magnetic and bathymetric data across the chile rise indicate that there is little change in the configuration of the spreading center as the ridge approaches and then collides with the trench in the collision zone a normal looking rift valley can be traced for 40 km before it disappears beneath the toe of the landward trench slope there is no evidence for the complex pattern of ridge jumps and spreading center rotations that occurred when the pacific farallon spreading center collided with north america in contrast the overriding plate is greatly affected by the ridge collision the landward trench slope steepens and narrows as the collision zone is approached at the southern end of the collision zone a ridge associated with the taitao fracture zone is colliding with the trench and may be in the process of being obducted onto the landward trench slope geothermal measurements were made along three transects of the margin corresponding to the time 3 m y before the collision during the collision and 6 m y after the collision the heat flow measurements in the collision zone document a large pulse of heat associated with the subduction of the ridge values as high as 345 mw m2 were recorded on the lower trench slope a prominent bottom simulating reflector bsr observed over a wide area of the landward trench slope north of the triple junction and more locally south of the triple junction is used to expand the grid of heat flow observations excellent agreement is found between measured heat flow and estimates of heat flow based on the depth to the bsr the heat flow measurements compare favorably with a theoretical model assuming conductive heat flow we estimate that the accretionary prism has been substantially removed in the collision zone and conclude that the landward trench slope is undergoing an episode of rapid tectonic erosion periodic ridge collisions in the past may account for the apparent truncation of the andean forearc region

modes and times of caldera resurgence the 10 ka evolution of ischia caldera italy from high precision archaeomagnetic dating abstract ischia is a well exposed and densely populated late quaternary caldera in the campanian magmatic province of italy ischia caldera experienced an average uplift rate of 3 3 cm year in the last ca 30 ka and is still actively resurging during the last 10 ka coeval with the resurgence a volcanic field of alkali trachytic to trachyandesitic lava domes lava flows tuff and scoria rings and pumice cones developed mainly on the eastern sector of the caldera along both resurgence related faults and regional nnw and ne striking faults in order to improve both our understanding of the recent volcanic history and the evaluation of future risks on ischia island a high precision archaeomagnetic dating method was applied to the products of 12 volcanic centres with probable age 14 c ages archaeomagnetic and other geochronological data as well as stratigraphic constraints show that during the studied time interval ischia volcanism occurred in five periods separated by phases of quiescence and coeval with earthquake and landslide events this fact suggests a pulsating mode of uplifting and deformations of the ischia resurgence during the two oldest periods of activity 7200 6800 bc and 4100 2300 bc resurgence probably produced a dome shaped structure location and geometry of vents suggest the occurrence of magma uprise along the fractures produced by bending of the overburden crustal block most of magma was emplaced as intrusions at the interior of the resurgent block whereas volcanism was represented by very viscous differentiated and crystallized lavas that emplaced as domes and high aspect ratio flows the resurgent dome caused recurrent lateral collapses that removed about 2 5 km 3 of rocks during the three youngest periods of activity 1800 1000 bc 650 bc 355 ad and 1302 ad resurgence affected a fault bounded asymmetric block this resulted from both a hydrostatic rebound of the crustal block after removal of material involved into huge slope instability triggering an increase in uplift rate and b new influxes of less evolved magma batches into the shallow reservoir that in turn favoured the intense volcanism of the last 4000 years

observations of high frequency scattered waves using dense arrays at teide volcano abstract a study of the seismic coda of seven small earthquakes recorded on the teide volcano canary islands spain was carried out using a temporary small aperture 12 station seismic array the purpose was to measure backazimuth apparent velocity and the type of waves that compose the coda in the frequency range 4 to 6 hz we used the zero lag cross correlation zlc method to obtain the components of the wave vector and three component analysis techniques based on the covariance matrix of the signal in the time domain to investigate the polarization properties of the signals the results show that a great part of the coda signals in the analyzed frequency range are almost uncorrelated while a low number of isolated correlated arrivals show an apparent slowness between 2 10 3 and 2 5 10 4 sec m and an almost random distribution of backazimuths the correlated arrivals have been interpreted as generated by strong scatterers probably related to the presence in the study area of surface topography irregularities the volcanic cone and the caldera rim the wave type varies from pure s type waves to mixed surface waves with some rayleigh components

the campi flegrei caldera italy formation and evolution in interplay with sea level variations since the campanian ignimbrite eruption at 39 ka abstract to date the origin of the campi flegrei caldera is still under debate and may be related to 1 a single caldera collapse associated with the neapolitan yellow tuff nyt eruption 2 two subsequent caldera collapses associated with the nyt and the preceding campanian ignimbrite ci eruptions forming a nested caldera complex or 3 not related to a caldera collapse after all here we study the submerged portion of the caldera which has favored a marine depositional setting and thus represents an ideal location for the reconstruction of its formation history utilizing multichannel seismic data volcanic deposits and edifices were seismically distinguished from sedimentary successions and the stratigraphy could be refined and extended back to the campanian ignimbrite eruption at 39 ka high resolution multichannel reflection seismic data revealed the existence of a nested caldera complex formed during the ci eruption at 39 ka and the more recent nyt eruption at 15 ka a ring fault bounding an inner caldera collapse structure was clearly imaged it appears that this inner ring fault was initially activated during the ci caldera collapse and later reactivated during the nyt caldera collapse with different amounts of subsidence the nyt caldera probably formed during an asymmetrical collapse with a maximum subsidence of 75 m in the offshore portion the vertical displacement related to the ci caldera collapse may be significantly larger the submerged caldera depression accommodates post eruption sediments within this high resolution archive two major unconformities developed at 8 6 ka and 5 ka when resurgence related uplift exceeded the rate of sea level rise concurrent with the emersion of the la starza terrace a previously unknown post collapse submarine volcanic mound located between nisida island and nisida bank probably formed between 4 8 and 3 7 ka also the penta palummo bank appears to be constructed of at least two monogenetic volcanic edifices the penta palummo volcano formed at 100 ka and a younger mushroom shaped unit formed between 39 and 15 ka the main outcome of this study is a conceptual evolutionary model providing novel insights on the formation of the campi flegrei nested caldera in the course of two large scale eruptions ci and nyt with associated caldera collapses along mutual i e reactivated faults and subsequent caldera resurgence as both the genesis and subsurface architecture of the campi flegrei caldera are still hotly debated topics in literature our discoveries can be regarded as a substantial advancement in the understanding of the campi flegrei volcanic area