modeling mountain building and the seismic cycle in the himalaya of nepal a host of information is now available regarding the geological and thermal structure as well as deformation rate across the himalaya of central nepal these data are reconciled in a two dimensional mechanical model that incorporates the rheological layering of the crust which depends on the local temperature and surface processes over geological timescale 5 ma the 20 mm yr estimated shortening rate across the range is accommodated by localized thrust faulting along the main himalayan thrust fault mht the mht reaches the surface along the foothills where it is called the main frontal thrust fault mft the mht flattens beneath the lesser himalaya and forms a midcrustal ramp at the front of the higher himalaya consistent with the river incision and the anticlinal structure of the lesser himalaya farther northward the mht roots into a subhorizontal shear zone that coincides with a midcrustal seismic reflector aseismic slip along this shear zone is accommodated in the interseismic period by elastic straining of the upper crust increasing the coulomb stress beneath the front of the higher himalaya where most of the microseismic activity clusters negligible deformation of the hanging wall requires a low apparent friction coefficient less than 0 3 on the flat portion of the mht on the ramp might be as high as 0 6 sensitivity tests show that a rather compliant quartz rich rheology and a high radioactive heat production in the upper crust of 2 5 w m 3 is required erosion affects the thermal structure and interplays with crustal deformation a dynamic equilibrium is obtained in which erosion balances tectonic uplift maintaining steady state thermal structure topography and deformation field using a linear diffusion model of erosion we constrain the value of the mass diffusivity coefficient to 0 5 1 6 l0 4 m 2 yr this study demonstrates that the data are internally consistent and compatible with current understanding of the mechanics of crustal deformation and highlight the role of viscous flow in the lower crust and of surface erosion in orogeny processes on the long term as well as during interseismic period

anomalies on geologic maps from multispectral and textural classification the bleida mining district morocco the variogram analysis of different reference image facies on a spot xs scene enables the definition of a discriminant texture index the values of this index are extracted locally in a given neighborhood this textural information is matched with the multispectral information in order to realize a combined classification the result is subsequently compared to a digitized geological map in order to extract zones of anomaly regarded as the differences between the combined classification and the geological map a majority of the anomalies require further interpretation while certain others are interesting to mineral prospecting

source parameters and time dependent slip distributions of slow slip events on the cascadia subduction zone from 1998 to 2008 1 we invert for the time dependent slip history of slow slip events on the cascadia subduction zone using gps data from 1998 to 2008 the 16 slip transients have sufficient station coverage to solve for the slip distribution on the plate interface gps time series are inverted for fault slip using the extended network inversion filter limited station coverage south of portland 45 5 n latitude restricts our analysis to events on the northern half of the subduction zone slip is resolved at the base of the seismogenic zone and the slip distributions suggest a potential segment boundary near seattle 47 6 n that correlates roughly with geologic and tectonic boundaries events that initiate to the north and south tend to overlap at about this latitude we compile statistics on source parameters such as propagation rate recurrence interval and stress drop which can be used to constrain proposed models of the source mechanics over a 10 year period total strain release from slow slip events is nonuniform along strike with the greatest cumulative slip 27 cm centered beneath port angeles 48 1 n this slip patch also exhibits the most regular recurrence of mw 6 events relative to other locations along strike the spatial extent of the slip patch beneath port angeles correlates with the along strike bend of the cascadia subduction zone in northwestern washington suggesting that plate geometry plays an important role in controlling the along strike characteristics of slow slip

variety and spatial heterogeneity of tectonic tremor worldwide 1 deep tectonic tremor occurs at various sites worldwide and the source characteristics are heterogeneous even at small scales in this study tremor sources were determined using data sets of seismic waveforms from various locations some of which are not well recognized as being a site of tremor activity the regions of interest are subduction zones at nankai and kyushu japan cascadia in western north america mexico southern chile and new zealand tremor locations are consistent with the geometry of subducting plates and tremor depth tends to be shallower where younger plate is subducting tremor duration shows a negative correlation with tidal sensitivity as reported previously in western shikoku in the nankai subduction zone the duration is also related to the width of the tremor zone in the dip direction of the subducting plate linear structures resembling striations are observed in hypocenter distributions in western shikoku and in the probability of short term migration direction these spatial heterogeneities in tremor activity may be controlled by the heterogeneous frictional properties of the downgoing slab due to long term plate subduction thereby providing information on the maturity of the plate interface

fault slip and seismic moment of the 1700 cascadia earthquake inferred from japanese tsunami descriptions 1 the 1700 cascadia earthquake attained moment magnitude 9 according to new estimates based on effects of its tsunami in japan computed coseismic seafloor deformation for hypothetical ruptures in cascadia and tsunami modeling in the pacific ocean reports of damage and flooding show that the 1700 cascadia tsunami reached 1 5 m heights at seven shoreline sites in japan three sets of estimated heights express uncertainty about location and depth of reported flooding landward decline in tsunami heights from shorelines and post 1700 land level changes we compare each set with tsunami heights computed from six cascadia sources each source is vertical seafloor displacement calculated with a three dimensional elastic dislocation model for three sources the rupture extends the 1100 km length of the subduction zone and differs in width and shallow dip for the other sources ruptures of ordinary width extend 360 670 km to compute tsunami waveforms we use a linear long wave approximation with a finite difference method and we employ modern bathymetry with nearshore grid spacing as small as 0 4 km the various combinations of japanese tsunami heights and cascadia sources give seismic moment of 1 9 1022 n m equivalent to moment magnitude 8 7 9 2 this range excludes several unquantified uncertainties the most likely earthquake of moment magnitude 9 0 has 19 m of coseismic slip on an offshore full slip zone 1100 km long with linearly decreasing slip on a downdip partial slip zone the shorter rupture models require up to 40 m offshore slip and predict land level changes inconsistent with coastal paleoseismological evidence

imaging the source region of cascadia tremor and intermediate depth earthquakes the subduction of hydrated oceanic crust releases volatiles that weaken the plate boundary interface trigger earthquakes and regulate transient phenomena such as episodic tremor and slip ets it is not clear how dehydration can separately induce earthquakes within the subducting plate and ets partly because few data exist on their relationship to subduction zone structures we present results of a seismic experiment in the washington cascades united states that images a region producing both earthquake types migration of scattered teleseis mic waves provides images of low velocity subducting crust at depths 40 45 km with sharp boundaries above and below it the sharp upper boundary indicates a layer of weak sediment or an overpressured fault zone that terminates abruptly downdip at 40 45 km depth regular earthquakes are at the top of the mantle within the downgoing plate everywhere the plate is 95 km deep but ets only exists where the sharp upper boundary occurs the ets location supports models of slow slip that require near lithostatic fluid pressure whereas regular earthquakes nucleate closer to the origin of metamorphic dehydration very low shear stresses on the plate boundary may limit seismicity to ets and similar phenomena

investigation of land ice ocean interaction with a fully coupled ice ocean model 1 model description and behavior 1 antarctic ice shelves interact closely with the ocean cavities beneath them with ice shelf geometry influencing ocean cavity circulation and heat from the ocean driving changes in the ice shelves as well as the grounded ice streams that feed them we present a new coupled model of an ice stream ice shelf ocean system that is used to study this interaction the model is capable of representing a moving grounding line and dynamically responding ocean circulation within the ice shelf cavity idealized experiments designed to investigate the response of the coupled system to instantaneous increases in ocean temperature show ice ocean system responses on multiple timescales melt rates and ice shelf basal slopes near the grounding line adjust in 1 2 years and downstream advection of the resulting ice shelf thinning takes place on decadal timescales retreat of the grounding line and adjustment of grounded ice takes place on a much longer timescale and the system takes several centuries to reach a new steady state during this slow retreat and in the absence of either an upward or downward sloping bed or long term trends in ocean heat content the ice shelf and melt rates maintain a characteristic pattern relative to the grounding line

spatial correlation of interseismic coupling and coseismic rupture extent of the 2011 mw 9 0 tohoku oki earthquake 1 imaging the extent to which the rupture areas of great earthquakes coincide with regions of pre seismic interplate coupling is central to understanding patterns of strain accumulation and release through the earthquake cycle both geodetic and seismic estimates of the coseismic rupture extent for the march 11 2011 mw 8 9 9 0 earthquake tohoku oki earthquake may be spatially correlated 0 26 0 05 to 0 82 0 05 with a region estimated to be partially to fully coupled in the interseismic period preceding the earthquake though there is substantial variation in the estimated distribution and magnitude of coseismic slip the 400 km long region estimated to have slipped 4 m corresponds to an area of the subduction zone interface that was coupled at 30 of long term plate convergence rate with peak slip near a region coupled 80 the northern termination of rupture is collocated with a region of relatively low 20 interseismic coupling near the epicenter of the 1994 mw 7 6 sanriku oki earthquake and near a region of potential long term low coupling or ongoing slow slip slip on the subduction interface beneath the coastline 40 50 km depth is best constrained by the land based gps data and least constrained on the shallowest portion of the plate interface due to the 230 km distance between geodetic observations and the japan trench

space geodesy constrains ice age terminal deglaciation the global ice 6g_c vm5a model a new model of the last deglaciation event of the late quaternary ice age is here described and denoted as ice 6g c vm5a it differs from previously published models in this sequence in that it has been explicitly refined by applying all of the available global positioning system gps measurements of vertical motion of the crust that may be brought to bear to constrain the thickness of local ice cover as well as the timing of its removal additional space geodetic constraints have also been applied to specify the reference frame within which the gps data are described the focus of the paper is upon the three main regions of last glacial maximum ice cover namely north america northwestern europe eurasia and antarctica although greenland and the british isles will also be included if peripherally in the discussion in each of the three major regions the model predictions of the time rate of change of the gravitational field are also compared to that being measured by the gravity recovery and climate experiment satellites as an independent means of verifying the improvement of the model achieved by applying the gps constraints several aspects of the global characteristics of this new model are also discussed including the nature of relative sea level history predictions at far field locations in particular the caribbean island of barbados from which especially high quality records of postglacial sea level change are available but which records were not employed in the development of the model although ice 6g c vm5a is a significant improvement insofar as the most recently available gps observations are concerned comparison of model predictions with such far field relative sea level histories enables us to identify a series of additional improvements that should follow from a further stage of model iteration

slip rate and tremor genesis in cascadia at many plate boundaries conditions in the transition zone between seismogenic and stable slip produce slow earthquakes in the cascadia subduction zone these events are consistently observed as slow aseismic slip on the plate interface accompanied by persistent tectonic tremor however not all slow slip at other plate boundaries coincides spatially and temporally with tremor leaving the physics of tremor genesis poorly understood here we analyze seismic geodetic and strainmeter data in cascadia to observe for the first time a large tremor generating slow earthquake change from tremor genic to silent and back again the tremor falls silent at reduced slip speeds when the migrating slip front pauses as it loads the stronger adjacent fault segment to failure the finding suggests that rheology and slip speed regulated stressing rate control tremor genesis and the same section of fault can slip both with and without detectable tremor limiting tremor s use as a proxy for slip

cascadia subducting plate fluids channelled to fore arc mantle corner ets and silica deposition in this study we first summarize the constraints that on the cascadia subduction thrust there is a 70 km gap downdip between the megathrust seismogenic zone and the episodic tremor and slip ets that lies further landward there is not a continuous transition from unstable to conditionally stable sliding seismic rupture occurs mainly offshore for this hot subduction zone ets lies onshore we then suggest what does control the downdip position of ets we conclude that fluids from dehydration of the downgoing plate focused to rise above the fore arc mantle corner are responsible for ets there is a remarkable correspondence between the position of ets and this corner along the whole margin hydrated mineral assemblages in the subducting oceanic crust and uppermost mantle are dehydrated with downdip increasing temperature and seismic tomography data indicate that these fluids have strongly serpentinized the overlying fore arc mantle laboratory data show that such fore arc mantle serpentinite has low permeability and likely blocks vertical expulsion and restricts flow updip within the underlying permeable oceanic crust and subduction shear zone at the fore arc mantle corner these fluids are released upward into the more permeable overlying fore arc crust an indication of this fluid flux comes from low poisson s ratios and vp vs found above the corner that may be explained by a concentration of quartz which has exceptionally low poisson s ratio the rising fluids should be silica saturated and precipitate quartz with decreasing temperature and pressure as they rise above the corner

periodic slow earthquakes from the cascadia subduction zone continuous geodetic measurements from convergent margins have shown that deep transient creep events can release large amounts of strain energy without detectable seismic shaking and they are thus known as slow or silent earthquakes because subduction zones generate the largest earthquakes

interseismic coupling and slow slip events on the cascadia megathrust in this study we model geodetic strain accumulation along the cascadia subduction zone between 2007 0 and 2017 632 using position time series from 352 continuous gps stations first we use the secular linear motion to determine interseismic locking along the megathrust we determine two end member models assuming that the megathrust is either a priori locked or creeping which differ essentially along the trench where the inversion is poorly constrained by the data in either case significant locking of the megathrust updip of the coastline is needed the downdip limit of the locked portion lies 20 80 km updip from the coast assuming a locked a priori but very close to the coast for a creeping a priori second we use a variational bayesian independent component analysis vbica decomposition to model geodetic strain time variations an approach which is effective to separate the geodetic strain signal due to non tectonic and tectonic sources the slow slip events sses kinematics is retrieved by linearly inverting for slip on the megathrust the independent components related to these transient phenomena the procedure allows the detection and modelling of 64 sses which spatially and temporally match with the tremors activity sees and tremors occur well inland from the coastline and follow closely the estimated location of the mantle wedge corner the transition zone between the locked portion of the megathrust and the zone of tremors is creeping rather steadily at the long term slip rate and probably buffers the effect of sses on the megathrust seismogenic portion

2010 maule earthquake slip correlates with pre seismic locking of andean subduction zone analysis of seismic data and global positioning system gps observations from the region around the site of the recent magnitude 8 8 maule earthquake in chile have been combined with modelling of the andean subduction zone to build up a detailed image of the event the earthquake occurred at a seismic gap part of the fault that has not experienced a major earthquake since 1835 and it seems to have released most of the stress accumulated since that time the observations suggest that coseismic slip heterogeneity at the scale of single asperities and therefore the seismic potential of future great earthquakes might be anticipated by geodetic observations these authors show evidence for a high coherence between the slip distribution inferred from the 2010 maule earthquake chile and the patchwork of interseismic locking distribution derived from global positioning system observations during the previous decade their work suggests that coseismic slip heterogeneity at the scale of single asperities should indicate the seismic potential of future earthquakes which thus might be anticipated by geodetic observations the magnitude 8 8 maule chile earthquake of 27 february 2010 ruptured a segment of the andean subduction zone megathrust that has been suspected to be of high seismic potential1 2 3 4 5 6 it is the largest earthquake to rupture a mature seismic gap in a subduction zone that has been monitored with a dense space geodetic network before the event this provides an image of the pre seismically locked state of the plate interface of unprecedentedly high resolution allowing for an assessment of the spatial correlation of interseismic locking with coseismic slip pre seismic locking might be used to anticipate future ruptures in many seismic gaps6 7 8 9 10 11 12 given the fundamental assumption that locking and slip are similar this hypothesis however could not be tested without the occurrence of the first gap filling earthquake here we show evidence that the 2010 maule earthquake slip distribution correlates closely with the patchwork of interseismic locking distribution as derived by inversion of global positioning system gps observations during the previous decade the earthquake nucleated in a region of high locking gradient and released most of the stresses accumulated in the area since the last major event in 1835 two regions of high seismic slip asperities appeared to be nearly fully locked before the earthquake between these asperities the rupture bridged a zone that was creeping interseismically with consistently low coseismic slip the rupture stopped in areas that were highly locked before the earthquake but where pre stress had been significantly reduced by overlapping twentieth century earthquakes our work suggests that coseismic slip heterogeneity at the scale of single asperities should indicate the seismic potential of future great earthquakes which thus might be anticipated by geodetic observations

funginite resinite associations in coal abstract petrographic examination of a series of indian tertiary age coal seams some with 38 resinite mineral free basis showed that several had associations of funginite in resinite there are three primary mechanisms for incorporation of fungus into resinite 1 as an accidental inclusion for example as a fungal infection on an insect encapsulated in resin 2 infestation of fungus into tree whether into a wound or other point of entry with encapsulation by resin as a defense mechanism by the tree alternatively fungal hyphae could also be growing from another portion of the tree that is already infected and 3 resinicolous fungus using the resin as a food source in the coal samples examined while incidental fungal inclusion on invading insects cannot be eliminated we have seen no direct evidence for insects both encapsulation of invading fungus by resin and resinicolous relationships are among the associations observed

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

evidence for a large earthquake and tsunami 100 400 years ago on western vancouver island british columbia abstract a peaty marsh soil is sharply overlain by a sand sheet and intertidal mud at tidal marshes near tofino and ucluelet vancouver island british columbia foraminifera and vascular plant fossils show that the buried soil was submerged suddenly and was covered quickly by sand radiocarbon ages place this event between 100 and 400 yr ago the coastal subsidence suggested by the submergence occurred in an area of net late holocene emergence perhaps during the most recent great earthquake on the northern part of the cascadia subduction zone the sand sheet overlying the peaty soil records the tsunami triggered by this earthquake similar stratigraphic sequences of about the same age have been reported from estuaries along the outer coasts of washington and northern oregon suggesting that hundreds of kilometers of the cascadia subduction zone may have ruptured during one or a series of plate boundary earthquakes less than 400 yr ago

juan de fuca slab geometry and its relation to wadati benioff zone seismicity 1 a new model of the subducted juan de fuca plate beneath western north america allows first order correlations between the occurrence of wadati benioff zone earthquakes and slab geometry temperature and hydration state the geo referenced 3d model constructed from weighted control points integrates depth information from earthquake locations and regional seismic velocity studies we use the model to separate earthquakes that occur in the cascadia forearc from those that occur within the underlying juan de fuca plate and thereby reveal previously obscured details regarding the spatial distribution of earthquakes seismicity within the slab is most prevalent where the slab is warped beneath northwestern california and western washington suggesting that slab flexure in addition to expected metamorphic dehydration processes promotes earthquake occurrence within the subducted oceanic plate earthquake patterns beneath western vancouver island are consistent with slab dehydration processes conversely the lack of slab earthquakes beneath western oregon is consistent with an anhydrous slab double differenced relocated seismicity resolves a double seismic zone within the slab beneath northwestern california that strongly constrains the location of the plate interface and delineates a cluster of seismicity 10 km above the surface that includes the 1992 m7 1 mendocino earthquake we infer that this earthquake ruptured a surface within the cascadia accretionary margin above the juan de fuca plate we further speculate that this earthquake is associated with a detached fragment of former farallon plate other subsurface tectonic elements within the forearc may have the potential to generate similar damaging earthquakes

a silent slip event on the deeper cascadia subduction interface continuous global positioning system sites in southwestern british columbia canada and northwestern washington state usa have been moving landward as a result of the locked state of the cascadia subduction fault offshore in the summer of 1999 a cluster of seven sites briefly reversed their direction of motion no seismicity was associated with this event the sudden displacements are best explained by 2 centimeters of aseismic slip over a 50 kilometer by 300 kilometer area on the subduction interface downdip from the seismogenic zone a rupture equivalent to an earthquake of moment magnitude 6 7 this provides evidence that slip of the hotter plastic part of the subduction interface and hence stress loading of the megathrust earthquake zone can occur in discrete pulses

space time correlation of slip and tremor during the 2009 cascadia slow slip event 1 episodic tremor and slip ets involving transient deformations accompanied by emergent low frequency tremor occurs in subduction zones around the world ets events increase the shear stress on locked megathrusts and may potentially trigger damaging earthquakes despite the clear association of tremor and slip the physical relationship between them is unresolved tremor appears to result from slip on small asperities on the plate interface due to either creep on the surrounding fault or stress increases ahead of the propagating slow slip front previous studies of migrating slow slip events have not had sufficient spatial and temporal resolution to differentiate between these two models to address this we invert gps data from the august 2009 ets event in central cascadia for the space time evolution of fault slip rate we find a correlation in both space and time between tremor epicenters and the independently determined position of high fault slip rate this supports the first hypothesis that tremor asperities are loaded directly by slow slip rather than by stress increases ahead of the slip front and provides new insights into the mechanics of ets

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

slow slip evolves into megathrust earthquakes in 2d numerical simulations 1 slow slip events sse in many subduction zones incrementally stress the adjacent locked megathrust suggesting that they could potentially either trigger or evolve into damaging earthquakes we explore this with 2d quasi dynamic simulations with rate state friction dilatancy and coupled 1d pore fluid and heat transport steady state weakening friction allows transient slip to nucleate but is inhibited by dilatant strengthening and destabilized by thermal pressurization sse spontaneously nucleate in low effective stress velocity weakening lesvw regions if the dimension of the lesvw is relatively small the sse are trapped at its updip end imparting a strong stress concentration in the locked zone after several centuries sse penetrate into the region of higher effective stress where thermal pressurization eventually leads to dynamic rupture for larger lesvw regions sse tend to increase in length with time ultimately higher slip speeds enhance thermal weakening leading to dynamic instability within the sse zone in both cases the onset of the ultimate sse is essentially indistinguishable from preceding events