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

assessment of teleseismically determined source parameters for the april 25 2015 mw 7 9 gorkha nepal earthquake and the may 12 2015 mw 7 2 aftershock on april 25 2015 a major m w 7 9 thrust earthquake ruptured the deeper portion of the seismogenic plate boundary beneath nepal along which india is underthrusting eurasia an m w 7 2 aftershock on may 12 2015 extended the eastern down dip edge of the rupture these destructive events caused about 9000 fatalities and 23 000 injuries the overall rupture zone is about 170 km long and 40 80 km wide this region of the plate boundary previously experienced a large earthquake in 1833 and in 1934 a larger m s 8 0 event located to the east ruptured all the way to the surface the main himalayan thrust mht on which slip occurred in 2015 has a very low dip angle of 6 and the depth of the mainshock slip distribution is very shallow extending from 7 to 18 km the shallow dip and depth present challenges for resolving faulting characteristics using teleseismic data we analyze global teleseismic signals for the mainshock and aftershock to estimate source parameters evaluating the stability of various procedures used for remotely characterizing kinematics of such shallow faulting back projection and finite fault slip inversion are used to assess the spatio temporal rupture history and evidence for frequency dependent radiation along dip slip zone width constraints from near field geodetic observations are imposed on the preferred models to overcome some limitations of purely teleseismic methods radiated energy stress drop and moment rate functions are determined for both events

primary surface ruptures of the great himalayan earthquakes in 1934 and 1255 it is unclear where plate boundary thrusts generate giant rather than great earthquakes along the himalayas the source sizes and recurrence times of large seismic events are particularly uncertain since no surface signatures were found for those that shook the range in the twentieth century here we challenge the consensus that these events remained blind and did not rupture the surface we use geomorphological mapping of fluvial deposits palaeo seismological logging of river cut cliffs and trench walls and modelling of calibrated 14c ages to show that the mw 8 2 bihar nepal earthquake on 15 january 1934 did break the surface traces of the rupture are clear along at least 150 km of the main frontal thrust fault in nepal between 85 50 and 87 20 e furthermore we date collapse wedges in the sir valley and find that the 7 june ad 1255 earthquake an event that devastated kathmandu and mortally wounded the nepalese king abhaya malla also ruptured the surface along this stretch of the mega thrust thus in the past 1 000 years two great earthquakes 679 years apart rather than one giant eleventh century ad event contributed to the frontal uplift of young river terraces in eastern nepal the rare surface expression of these earthquakes implies that surface ruptures of other reputedly blind great himalayan events might exist the recurrence times of great himalayan earthquakes are difficult to assess because they rarely rupture the surface field mapping and 14c dating of offset fluvial deposits are used to identify two great himalayan quakes that ruptured the surface along the main plate boundary fault in ad 1255 and 1934

duplex in the main himalayan thrust illuminated by aftershocks of the 2015 m w 7 8 gorkha earthquake in april 2015 the lower locked portion of the main himalayan thrust ruptured beneath nepal causing the disastrous mw 7 8 gorkha earthquake elucidating the enigmatic geometry of this plate boundary fault is important for understanding the nucleation and arrest of large earthquake ruptures as well as the seismic hazard topography and tectonics of the himalaya here we interpret the geometry of the main himalayan thrust from the spatial distribution and rupture patterns of a dynamic sequence of aftershocks following the gorkha earthquake which were recorded by a rapidly deployed dense seismic network we find that the thrust comprises two north dipping subhorizontal planes that are connected by a system of bounded imbricate thrust faults this structure is known as a duplex we propose that this duplex acts as an impediment to plate convergence and accommodates tectonic stress along its complex system of faults such a prominent structure to the main himalayan thrust is consistent with surface geological studies but challenges geophysically derived conventional models with simpler geometries the main himalayan thrust comprises two fault planes connected by imbricated faults a structure that impedes convergence according to an analysis of the distribution and orientation of aftershocks of the 2015 gorkha earthquake in nepal

the 5 september 2012 nicoya costa rica mw 7 6 earthquake rupture process from joint inversion of high rate gps strong motion and teleseismic p wave data and its relationship to adjacent plate boundary interface properties 1 on 5 september 2012 a large thrust earthquake mw 7 6 ruptured a densely instrumented seismic gap on the shallow dipping plate boundary beneath the nicoya peninsula costa rica ground motion recordings directly above the rupture zone provide a unique opportunity to study the detailed source process of a large shallow megathrust earthquake using very nearby land observations hypocenter relocation using local seismic network data indicates that the event initiated with small emergent seismic waves from a hypocenter 10 km offshore 13 km deep on the megathrust a joint finite fault inversion using high rate gps strong motion ground velocity recordings gps static offsets and teleseismic p waves reveals that the primary slip zone slip 1 m is located beneath the peninsula the rupture propagated downdip from the hypocenter with a rupture velocity of 3 0 km s the primary slip zone extends 70 km along strike and 30 km along dip with an average slip of 2 m the associated static stress drop is 3 mpa the seismic moment is 3 5 1020 nm giving mw 7 6 the coseismic large slip patch directly overlaps an onshore interseismic locked region indicated by geodetic observations and extends downdip to the intersection with the upper plate moho at deeper depths below the upper plate moho seismic tremor and low frequency earthquakes have been observed most tremor locates in adjacent areas of the megathrust that have little coseismic slip a region of prior slow slip deformation to the southeast also has no significant coseismic slip or aftershocks an offshore locked patch indicated by geodetic observations does not appear to have experienced coseismic slip and aftershocks do not overlap this region allowing the potential for a comparable size rupture offshore in the future

imaging the moho and the main himalayan thrust in western nepal with receiver functions the crustal structure of western nepal is studied for the first time by performing receiver function analysis on teleseismic waveforms recorded at 16 seismic stations the moho geometry is imaged as it deepens from 40 km depth beneath the foothills and the lesser himalaya to 58 km depth beneath the higher himalayan range a midcrustal low velocity zone is detected at 15 km depth along 55 km horizontal distance and is interpreted as the signature of fluids expelled from rocks descending in the footwall of the main himalayan thrust our new image allows structural comparison of the moho and of the main himalayan thrust geometry along strike of the himalayas and documents long wavelength lateral variations the general crustal architecture observed on our images resembles that of central nepal therefore western nepal is also expected to be able to host large m w 8 megathrust earthquakes as the 1505 ce event plain language summary we investigate the earth s crust in western nepal we use 16 stations to detect waves from faraway earthquakes and using these we make an image of the structure of the crust we find that the crust is getting thicker beneath the himalaya and that fluids are locally present in the upper part of the crust as this image is similar to the one we know in central nepal where major earthquakes happened it is possible that major earthquakes will also happen in western nepal in the future without knowing when

convergence rate across the nepal himalaya and interseismic coupling on the main himalayan thrust implications for seismic hazard we document geodetic strain across the nepal himalaya using gps times series from 30 stations in nepal and southern tibet in addition to previously published campaign gps points and leveling data and determine the pattern of interseismic coupling on the main himalayan thrust fault mht the noise on the daily gps positions is modeled as a combination of white and colored noise in order to infer secular velocities at the stations with consistent uncertainties we then locate the pole of rotation of the indian plate in the itrf 2005 reference frame at longitude 1 34 3 31 latitude 51 4 0 3 with an angular velocity of 0 5029 0 0072 myr the pattern of coupling on the mht is computed on a fault dipping 10 to the north and whose strike roughly follows the arcuate shape of the himalaya the model indicates that the mht is locked from the surface to a distance of approximately 100 km down dip corresponding to a depth of 15 to 20 km in map view the transition zone between the locked portion of the mht and the portion which is creeping at the long term slip rate seems to be at the most a few tens of kilometers wide and coincides with the belt of midcrustal microseismicity underneath the himalaya according to a previous study based on thermokinematic modeling of thermochronological and thermobarometric data this transition seems to happen in a zone where the temperature reaches 350 c the convergence between india and south tibet proceeds at a rate of 17 8 0 5 mm yr in central and eastern nepal and 20 5 1 mm yr in western nepal the moment deficit due to locking of the mht in the interseismic period accrues at a rate of 6 6 0 4 10 19 nm yr on the mht underneath nepal for comparison the moment released by the seismicity over the past 500 years including 14 m w 7 earthquakes with moment magnitudes up to 8 5 amounts to only 0 9 10 19 nm yr indicating a large deficit of seismic slip over that period or very infrequent large slow slip events no large slow slip event has been observed however over the 20 years covered by geodetic measurements in the nepal himalaya we discuss the magnitude and return period of m 8 earthquakes required to balance the long term slip budget on the mht

interseismic strain accumulation on the himalayan crustal ramp nepal the departement of mines and geology has been monitoring the seismicity of the central himalayas of nepal since 1985 intense microseismicity and frequent medium size earthquakes ml 4 tend to cluster beneath the topographic front of the higher himalaya this 10 20km deep seismicity also correlates with a zone of localized uplift that has been evidenced from geodetic data both microseismic and geodetic data indicate strain accumulation on a mid crustal ramp that had been previously inferred from geological and geophysical evidence this ramp connects a flat decollement under the lesser and sub himalaya with a deeper decollement under the higher himalaya and probably acts as a geometric asperity where strain and stress build up during the interseismic period the large himalayan earthquakes could nucleate there and probably activate the whole flat and ramp system up to the blind thrusts of the sub himalaya

interseismic coupling and asperity distribution along the kamchatka subduction zone 1 gps measurements of interseismic horizontal surface velocities reveal the degree of kinematic coupling of the plate boundary thrust along the kamchatka subduction zone from about 51 degrees to 57 degrees n latitude inversions for the distribution of aseismic slip rate along the similar to 15 degrees nw dipping underthrust suggest a nonslipping plate interface in southern kamchatka above similar to 50 km depth along the segment that ruptured in the m w 9 1952 earthquake north of similar to 53 degrees n the subduction interface experiences significant aseismic slip consistent with the lower seismic moment release in m 8 5 earthquakes along this portion of the subduction zone the gps velocities are consistent with a boundary element forward model in which historic earthquake rupture zones are represented as locked asperities surrounded by a zero shear stress subduction interface loaded by plate convergence models in which the complete rupture zones of historic earthquakes are considered locked greatly overpredict the degree of kinematic coupling reducing the area of the locked model asperities to the central 25 area of historic rupture zones fits the data well suggesting that large earthquakes involve small fully locked core asperities surrounded by conditionally stable portions of the plate interface areas of low aseismic slip rate appear to be roughly correlated with areas of low isostatic gravity anomalies over offshore forearc basins while less coupled portions of the kamchatka subduction zone coincide with high gravity anomalies offshore of two peninsulas possibly related to the subduction of the emperor meji seamount chain and the kruzenstern fracture zone

petroleum geology of pennsylvanian and lower permian strata tucumcari basin east central new mexico abstract the tucumcari basin of east central new mexico is a structural depression that existed as a depositional basin from strawn middle pennsylvanian until late wolfcampian early permian time depth to precambrian ranges from 6500 ft to more than 9000 ft high angle faults form the north east and west edges of the basin no major structural discontinuities separate the basin from a shallow shelf to the south faults cut pennsylvanian and wolfcampain strata but generally do not offset post wolfcampian strata faults and regional structure control facies and thickness patterns within the pennsylvanian and wolfcampian coarse arkosic sands in the pennsylvanian and wolfcampian are good reservoirs those sands were deposited in northern and western parts of the basin and were derived from highlands of precambrian granitic rocks that formed the northern and western margins of the basin high energy limestones are possible reservoirs in the southern part of the basin porous wolfcampian dolostones cover the frio uplift on the east side of the basin pennsylvanian and wolfcampian marine shales and micritic limestones are source rocks within the tucumcari basin post wolfcampian strata are thermally immature two presently noncommercial pools of oil and gas the latigo ranch pool and the t 4 ranch pool have been discoveredmore in strawn sands in the northern part of the tucumcari basin oil generated in the upper paleozoic has migrated vertically into the triassic two oil accumulations the santa rosa tar sands and the newkirk oil pool are in triassic sandstones and have combined reserves of 153 million bbls of oil less

the mechanism of partial rupture of a locked megathrust the role of fault morphology assessment of seismic hazard relies on estimates of how large an area of a tectonic fault could potentially rupture in a single earthquake vital information for these forecasts includes which areas of a fault are locked and how the fault is segmented much research has focused on exploring downdip limits to fault rupture from chemical and thermal boundaries and along strike barriers from subducted structural features yet we regularly see only partial rupture of fully locked fault patches that could have ruptured as a whole in a larger earthquake here we draw insight into this conundrum from the 25 april 2015 m w 7 8 gorkha nepal earthquake we invert geodetic data with a structural model of the main himalayan thrust in the region of kathmandu nepal showing that this event was generated by rupture of a decollement bounded on all sides by more steeply dipping ramps the morphological bounds explain why the event ruptured only a small piece of a large fully locked seismic gap we then use dynamic earthquake cycle modeling on the same fault geometry to reveal that such events are predicted by the physics depending on the earthquake history and the details of rupture dynamics however great earthquakes that rupture the entire seismogenic zone are also possible these insights from nepal should be applicable to understanding bounds on earthquake size on megathrusts worldwide

large paleoearthquake timing and displacement near damak in eastern nepal on the himalayan frontal thrust an excavation across the himalayan frontal thrust near damak in eastern nepal shows displacement on a fault plane dipping 22 has produced vertical separation across a scarp equal to 5 5 m stratigraphic structural geometrical and radiocarbon observations are interpreted to indicate the displacement is the result of a single earthquake of 11 3 3 5 m of dip slip displacement that occurred 1146 1256 ad empirical scaling laws indicate that thrust earthquakes characterized by average displacements of this size may produce rupture lengths of 450 800 km and moment magnitudes mw of 8 6 to 9 sufficient strain has accumulated along this portion of the himalayan arc during the roughly 800 years since the 1146 1256 ad earthquake to produce another earthquake displacement of similar size

measuring ground displacements from sar amplitude images application to the landers earthquake ers sar amplitude images are utilized to map ground displacements from a sub pixel correlation method it yields a ground two dimensional displacement field with independent measurements about every 128m in azimuth and 250m in range the accuracy depends on the characteristics of the images for the landers test case the 1 uncertainty is 0 8m in range and 0 4m in azimuth we show that this measurement provides a map of major surface fault ruptures accurate to better than 1km and information on coseismic deformation comparable to the 92 gps measurements available although less accurate this technique is more robust than sar interferometry and provides complementary information since interferograms are only sensitive to the displacement in range

ductile shearing to brittle thrusting along the nepal himalaya linking miocene channel flow and critical wedge tectonics to 25th april 2015 gorkha earthquake the 25th april 2015 magnitude 7 8 gorkha earthquake in nepal ruptured the main himalayan thrust mht for 140 km east west and 50 km across strike the earthquake nucleated at a depth of 15 18 km approximating to the brittle ductile transition and propagated east along the mht but did not rupture to the surface leaving half of the fault extent still locked beneath the siwalik hills coseismic slip shows that motion is confined to the ramp flat geometry of the mht and there was no out of sequence movement along the main central thrust mct below 20 km depth the mht is a creeping aseismic ductile shear zone cumulated deformation over geological time has exhumed the deeper part of the himalayan orogen which is now exposed in the greater himalaya revealing a tectonic history quite different from presently active tectonics there early miocene structures including the mct are almost entirely ductile with deformation occurring at temperatures higher than 400 c and were active between 22 16 ma kyanite and sillimanite grade gneisses and migmatites approximately 5 20 km thick in the core of the greater himalayan sequence ghs together with leucogranite intrusions along the top of the ghs were extruded southward between 22 15 ma concomitant with ages of partial melting thermobarometric constraints show that ductile extrusion of the ghs during the miocene occurred at muscovite dehydration temperatures 650 775 c and thus brittle thrusting and critical taper models for ghs deformation are unrealistic as partial melting and channel flow ceased at 15 ma brittle thrusting and underplating associated with duplex formation occurred along the lesser himalaya passively uplifting the ghs

towards inferring earthquake patterns from geodetic observations of interseismic coupling ultimately seismotectonic studies seek to provide ways of assessing the timing magnitude and spatial extent of future earthquakes ample observations document the spatial variability in interseismic coupling defined as a degree of locking of a fault during the period of stress build up between seismic events fully or nearly locked fault patches are often surrounded by aseismically creeping areas however it is unclear how these observations could help assess future earthquakes here we simulate spontaneous seismic and aseismic fault slip with a fully dynamic numerical model our simulations establish the dependence of earthquake rupture patterns and interseismic coupling on spatial variations of fault friction we consider the long term evolution of slip on a model fault where two seismogenic locked segments are separated by an aseismically slipping patch where rupture is impeded we find that the probability for a large earthquake to break through the rupture impeding patch is correlated with the interseismic coupling averaged over this patch in addition the probability that an earthquake breaks through the rupture impeding patch and interseismic coupling are both related to fault friction properties through a single non dimensional parameter our study opens the possibility of predicting seismic rupture patterns that a fault system can produce on the basis of observations of its interseismic coupling and suggests that regions of low interseismic coupling may reveal permanent barriers to large earthquakes

one dimensional reference velocity model and precise locations of earthquake hypocenters in the kumaon garhwal himalaya we report an optimum 1d velocity model of the kumaon garhwal himalaya based on analysis of 5000 seismograms generated at 50 digital broadband seismic stations in the region during april 2005 and june 2008 these include data from 385 local earthquakes with azimuth gaps of 180 each of which had a minimum of seven p and five s phase readings this velocity model was then used to map the locations of 1150 earthquakes 1 m 5 that occurred in the region a majority 75 of these earthquakes occurred in the upper 20 km of the crust largely above the main himalayan thrust and forms an 10 km wide band along the surface trace of the munsiari thrust however we also find another parallel band of earthquakes about 70 km to its southwest through the tehri dam in the lesser himalaya and a significant number 15 both in the tethys himalaya and within the flexed crust of the underthrusting indian plate beneath the ganga basin furthermore unlike the reported absence of seismicity in the lower crust of the nepal himalaya and its reappearance in the shallow mantle the kumaon garhwal himalayan crust is found to be seismically active throughout with no event reliably recorded in the shallow mantle of this region during the epoch of our seismic recording online material figures show constrains and distribution of data used in joint inversion and table shows p and s wave corrections obtained from travel time inversion

gradual unlocking of plate boundary controlled initiation of the 2014 iquique earthquake a long foreshock series unlocked the south american plate boundary until eventually initiating the m 8 1 iquique chile earthquake two groups publishing in this issue of nature analyse the seismic context of the iquique earthquake that occurred off the coast of northern chile on 1 april 2014 in a seismic zone that had been quiescent since a significant event in 1877 gavin hayes et al identify areas of remaining or elevated earthquake hazard along the megathrust fault in the region and conclude that the 2014 iquique event was not the earthquake that had been anticipated given that significant sections of the northern chile subduction zone have not ruptured in almost 150 years they suggest that it is likely that future megathrust earthquakes will occur south and potentially north of the 2014 iquique sequence bernd schurr et al show that the april 2014 earthquake broke a central fraction of the northern chile seismic gap the last major segment of the south american plate boundary that had yet to rupture in the past century from july 2013 up to the april earthquake they identify three seismic clusters along this part of the plate boundary each lasting a few weeks with earthquakes of increasing peak magnitudes they conclude that these seismic clusters and their slip transients reflect a gradual weakening of the central part of the seismic gap that was instrumental in initiating the final failure on 1 april 2014 northern chile was struck by a magnitude 8 1 earthquake following a protracted series of foreshocks the integrated plate boundary observatory chile monitored the entire sequence of events providing unprecedented resolution of the build up to the main event and its rupture evolution here we show that the iquique earthquake broke a central fraction of the so called northern chile seismic gap the last major segment of the south american plate boundary that had not ruptured in the past century1 2 since july 2013 three seismic clusters each lasting a few weeks hit this part of the plate boundary with earthquakes of increasing peak magnitudes starting with the second cluster geodetic observations show surface displacements that can be associated with slip on the plate interface these seismic clusters and their slip transients occupied a part of the plate interface that was transitional between a fully locked and a creeping portion leading up to this earthquake the b value of the foreshocks gradually decreased during the years before the earthquake reversing its trend a few days before the iquique earthquake the mainshock finally nucleated at the northern end of the foreshock area which skirted a locked patch and ruptured mainly downdip towards higher locking peak slip was attained immediately downdip of the foreshock region and at the margin of the locked patch we conclude that gradual weakening of the central part of the seismic gap accentuated by the foreshock activity in a zone of intermediate seismic coupling was instrumental in causing final failure distinguishing the iquique earthquake from most great earthquakes finally only one third of the gap was broken and the remaining locked segments now pose a significant increased seismic hazard with the potential to host an earthquake with a magnitude of 8 5

structural interpretation of the great earthquakes of the last millennium in the central himalaya abstract a major question about the himalaya remains open does a great earthquake like the mw 8 1 1934 earthquake release all the strain stored by the tibet india convergence during the preceding interseismic period and only that strain or can it also release a background store of energy that remained unreleased through one or more earlier earthquakes and so potentially engender giant events or a relatively random sequence of events to consider this question the history of the great earthquakes of the last millennium is investigated here by combining data provided by the historical archives of kathmandu trenches through surface ruptures isoseismal damage mapping seismites and the instrumental record in the kathmandu basin the location of the epicenter of the 1934 earthquake was determined from the arrival of high energy p waves that caused sedimentary dikes and ground fractures perpendicular to the epicenter azimuth the epicenter of the mw 7 6 1833 earthquake can therefore be determined analogously from dike orientation and its location to the ne of kathmandu indicates an overlap with the mw 8 1 1934 rupture the 1934 earthquake released strain not released by the 1833 earthquake comparison of the historical records of earthquakes in kathmandu with 14c ages from paleo seismic trenches along the himalayan front suggests that 1 the 1344 kathmandu event ruptured the surface as far away as kumaon and was therefore a giant mw 8 6 earthquake and 2 the 1255 event that destroyed kathmandu is attested by surface ruptures in central and western nepal and by seismites in soft sediment as far away as kumaon geometric and rheologic controls for the different types of ruptures during the medium mw 7 great mw 8 and giant mw 8 4 earthquakes are illustrated in structural cross sections it is found that the epicenters of great himalayan earthquakes are located on the basal thrust farther north or close to the locked zone which is defined from geodetic measurements of regional deformation during the interseismic period this suggests that great earthquakes initiate in a wide transition zone between exclusively brittle and exclusively creeping regimes the extent of which depends on the dip of the main himalayan thrust the succession of the great earthquakes during the last millennium has released all the 20 m millennial himalayan convergence even in the central seismic gap which has been locked since 1505 the millennial seismic release rate is close to the convergence rate nonetheless no evidence of a succession of characteristic earthquakes has been found the 1100 1833 and 1934 earthquakes in the eastern himalaya are characterized neither by constant displacement nor by constant recurrence furthermore some great earthquakes do not release all the strain elastically stored by the himalayan and tibetan upper crust after the 1255 event there was still the potential for a slip of several meters for the mw 8 1 1505 event this suggests a rather random release of seismic energy great earthquakes could occur anytime and in any part of the central himalaya furthermore a future giant earthquake of mw 8 6 cannot be excluded

evidence for a great medieval earthquake 1100 a d in the central himalayas nepal the himalayan orogen has produced three thrust earthquakes with moment magnitude m w 7 8 to 8 5 during the past century yet no surface ruptures associated with these great earthquakes have been documented here we present paleoseismic evidence from east central nepal that since 700 a d a single earthquake ruptured the frontal thrust fault at 1100 a d with a surface displacement of 17 5 3 meters and a lateral extent and size that could have exceeded 240 kilometers and m w 8 8 respectively ruptures associated with m w

strong seismic coupling underneath garhwal kumaun region nw himalaya india abstract no great earthquake has occurred in the garhwal kumaun region of nw himalaya in the past 500 years or more we report results of continuous gps measurements from 28 sites from the region to suggest that the convergence rate in this part of the himalaya is about 18 mm yr which is leading to strain accumulation in the region the main himalayan thrust mht in the frontal part of the himalaya under the outer and southern lesser himalaya is strongly coupled for a width of about 85 km the midcrustal ramp where earthquakes of himalayan seismic belt occur exhibits low coupling strong coupling on the mht beneath the outer and lesser himalaya is homogeneous except in the very shallow updip part of the mht subduction of sediments of the indo gangetic plains or the delhi hardwar ridge does not seem to influence coupling a high rate of strain accumulation which has continued for more than 500 years on a strongly coupled mht makes this one of the most earthquake vulnerable segments of the himalayan arc

the 2016 kumamoto mw 7 0 earthquake a significant event in a fault volcano system the 2016 kumamoto earthquake sequence occurred on the futagawa hinagu fault zone near the aso volcano on kyushu island the sequence was initiated with two major m w 6 0 foreshocks and the mainshock m w 7 0 occurred 25 h after the second major foreshock we combine gps strong motion synthetic aperture radar images and surface offset data in a joint inversion to resolve the kinematic rupture process of the mainshock and coseismic displacement of the foreshocks the joint inversion results reveal a unilateral rupture process for the mainshock involving sequential rupture of four major asperities the slip area of the foreshocks and mainshock and the aftershock loci form a detailed complementary pattern the mainshock rupture terminates near the rim of the caldera leaving a 10 km long gap of aftershocks this area is characterized by high temperature and low shear wave velocity density and resistivity which may be related to the partially melted geothermal condition ductile material property near the volcano may act as a material barrier to the dynamic rupture topographic weight of the caldera increases compressional normal stress on the fault plane which may behave as a stress barrier long term seismic hazard and deformation behaviors related to these two types of barriers are discussed in terms of the associated frictional mechanism significant postseismic creeps observed near the volcano area indicates a velocity strengthening frictional behavior near the rupture termination which confirms that the material barrier mechanism is likely the dominant rupture termination mechanism

electrical structure of the himalaya of central nepal high conductivity around the mid crustal ramp along the mht twelve broadband magnetotelluric mt soundings were performed across the himalaya of central nepal in 1996 in order to determine the electrical structure of the crust and its relation to geological structures and active tectonics the mt impedance tensors were obtained for frequencies between 0 001 and 500 hz the 2 d section derived from joint inversion of te and tm mode after rri and groom bailey decomposition shows high conductivity in the foreland basin 30 m that contrasts with the resistive indian basement 300 m and lesser himalaya 1000 m in addition our mt sounding reveals a major conductive feature beneath the front of the higher himalaya also characterized by intense microseismic activity and the position of a mid crustal ramp along the major active thrust fault mht this high conductivity zone probably reflects metamorphic fluids released during underthrusting of the indian basement and pervading well connected microcracks induced by interseismic stress build up or distributed brittle deformation around the ramp

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

fault slip distribution of the 1999 mw 7 1 hector mine california earthquake estimated from satellite radar and gps measurements we use a combination of satellite radar and gps data to estimate the slip distribution of the 1999 m w 7 1 hector mine earthquake a right lateral strikeslip earthquake that occurred on a northwest southeast striking fault in the southern california mojave desert the data include synthetic aperture radar interferograms insar from both ascending and descending orbits radar amplitude image offset fields sario for both ascending and descending azimuth directions and campaign gps observations from 55 stations provided by agnew et al 2002 we model the fault with nine segments derived from the field mapped fault rupture the sario data and aftershock locations we first estimate the dip of each fault segment as well as a single constant strike slip component across each segment resulting in an average dip of 83 to the northeast and slip of up to 5 6 m then we fix the optimal fault segment dip discretize the fault segments into 1 5 km 1 5 km patches and solve for the variable slip distribution using a nonnegative least squares method that includes an appropriate degree of smoothing our preferred solution has both right lateral strike slip and reverse faulting the estimated geodetic moment is 5 93 10 19 n m m w 7 1 similar to seismological estimates indicating that there are insignificant interseismic and postseismic deformation signals in the data we find strike slip displacements of up to 6 0 m and reverse faulting of up to 1 6 m with the maximum slip located just northwest of the epicenter most of the slip is concentrated northwest and south of the epicenter little slip is found on the northeastern branch of the fault the sario data and our modeling indicate that the amount and extent of surface fault rupture were underestimated in the field

paleoseismic evidence of great surface rupture earthquakes along the indian himalaya 1 toward understanding the relationship between strain accumulation and strain release in the context of the mechanics of the earthquake and mountain building process and quantifying the seismic hazard associated with the globes largest continental thrust system we describe the late quaternary expression and paleoseismic evidence for great surface rupture earthquakes at six sites along the himalayan frontal thrust hft system of india our observations span a distance of 250 km along strike of the hft uplifted and truncated fluvial terrace deposits resulting from the holocene displacements on the hft are preserved along canyons of the ghaggar markanda shajahanpur and kosi rivers dividing the elevation of the bedrock straths at each site by their ages yields estimates of the vertical uplift rate of 4 6 mm yr which when assumed to be the result of slip on an underlying thrust dipping at 20 45 are equivalent to fault slip rates of 6 18 mm yr or shortening rates of 4 16 mm yr trench exposures reveal the hft to fold and break late holocene surface sediments near the cities and villages of chandigarh kala amb rampur ganda lal dhang and ramnagar radiocarbon ages of samples obtained from the displaced sediments indicate surface rupture at each site took place after a d 1200 and before a d 1700 uncertainties attendant to the radiocarbon dating currently do not allow an unambiguous definition of the capping bound on the age of the displacement at each site and hence whether or not the displacements at all sites were contemporaneous trench exposures and vertical separations measured across scarps at rampur ganda lal dhang and ramnagar are interpreted to indicate single event displacements of 11 38 m dividing the observed single event vertical components of displacement by the estimated longer term uplift rates indicates 1330 3250 or more years should be required to accumulate the slip sufficient to produce similar sized displacements surface rupture appears to not have occurred during the historical 1905 kangra mw 7 7 1934 bihar nepal mw 8 1 and 1950 assam mw 8 4 earthquakes which also occurred along the himalayan front yet we observe clear evidence of fault scarps and displacements in young alluvium and progressive and continued offset of fluvial terrace deposits along the hft we suggest on this basis and the size and possible synchroneity of displacements recorded in the trenches that there exists the potential for earthquakes larger than recorded in the historical record and with the potential to rupture lengths of the hft greater than the 250 km we have studied

the 2018 mw6 4 hualien earthquake dynamic slip partitioning reveals the spatial transition from mountain building to subduction abstract the hualien earthquake of 6 february 2018 occurred on the east coast of taiwan in the collision zone between the eurasia and philippine sea plates even though an event with a moderate magnitude of mw6 4 it caused widespread surface rupture in and around the city of hualien leading to significant loss of life and property damages in this study we investigate the rupture process of this earthquake using teleseismic strong motion gps and sar observations we first employ the geodetic data to search for the geometry of multiple fault planes followed by a multi point source inversion using teleseismic and strong motion waveforms to determine the possible change in mechanism during the rupture process we parameterize the fault model with three curved segments to represent the complicated fault system involved in the rupture process an offshore fault was ruptured by the hualien earthquake showing a curved geometry with the dip angle changing from shallow 45 in the north to sub vertical 80 in the south evidence from tomographic velocity model indicates that the offshore fault is the boundary between the subducted coastal range and meilun tableland which is also a plate interface finally a joint inversion is carried out to constrain the complete rupture process the inversion result reveals that the rupture initiated on the offshore fault interface and propagated unilaterally towards south southwest the rupture then jumped on to the sub vertical meilun fault and reached the ground surface significant dynamic slip partitioning occurred during the rupture process near the meilun fault with the majority of thrust and strike slip motions occurring on the off shore fault and the meilun fault respectively the afterslip model is inverted with sar timeseries spanning 6 months after the mainshock the afterslip model shows a strong complementary pattern with co seismic slips the comparison between the 1951 and 2018 hualien earthquake shows that the rupture behavior of the meilun fault does not follow a repeating earthquake pattern thrust motion may occur on the two side walls of the meilun tableland with a teeterboard uplift pattern in different earthquake cycles the change in geometry and slip on the inter plate fault appears to be related to the transition of the coastal range from mountain building to subduction

implications for elastic energy storage in the himalaya from the gorkha 2015 earthquake and other incomplete ruptures of the main himalayan thrust abstract rupture in the 2015 m7 8 gorkha earthquake nucleated at the downdip edge of the main himalayan thrust mht near the transition from interseismic locking to aseismic creep beneath the tibetan plateau and propagated incompletely towards the main frontal thrusts mft despite the imposition of a substantial static strain in the mid decollement afterslip on the mht within a year of the earthquake had decayed to negligible levels earthquakes that incompletely rupture the mht 7

more major earthquakes at the nepal himalaya study on coulomb stress perspective abstract on april 2015 a major earthquake of 7 9 mw occurred in the nepal himalaya followed by 553 earthquakes of local magnitude greater than 4 0 within the first 43 days including another major event of 7 3 mw we resolve the static coulomb failure stress cfs change onto the finite fault models of 7 9 mw after elliott et al 2016 and galezka et al 2015 and its effect on associated receiver faults correlation of aftershocks with the enhanced cfs condition shows that the elliott et al 2016 model explains 60 4 and the galezka et al 2015 model explains about 47 7 of the aftershocks in high stress regions aftershocks were poorly spatially correlated with the enhanced cfs condition after the 7 9 mw main shock and can be explained by correlation with release of seismic energy from the associated secondarily stressed prominent thrust planes and transverse faults stress resolved on the associated receiver faults show increased stress on both transverse and thrust fault systems with the potential of triggering significant aftershocks or subsequent main shocks

the northern cascadia subduction zone at vancouver island seismic structure and tectonic history the structure and tertiary tectonic history of the northern cascadia subduction zone have been delineated by a series of new multichannel seismic lines acquired across the continental shelf to the

bridging the gap between seismically and geodetically detected slow earthquakes 1 recently observed unusual seismic events deep low frequency earthquakes and tremor very low frequency earthquakes and slow slip events form a family of slow earthquakes that are governed by a scaling law different from ordinary earthquakes guided by this scaling law we have observed previously unknown events with source duration of 20 200 s and moment magnitude of 3 4 under the kii peninsula in western japan these events radiate seismic energy in the 2 8 hz band in direct proportion to their seismic moment rate though the constant of proportionality is 4 5 orders of magnitude smaller than for ordinary earthquakes slow earthquakes in this region may be comprised of numerous shear slip events which in aggregate are manifest seismically as tremor or isolated longer events up to 200 s and geodetically as slow slip events

high resolution backprojection at regional distance application to the haiti m7 0 earthquake and comparisons with finite source studies a catastrophic m w7 earthquake ruptured on 12 january 2010 on a complex fault system near port au prince haiti offshore rupture is suggested by aftershock locations and marine geophysics studies but its extent remains difficult to define using geodetic and teleseismic observations here we perform the multitaper multiple signal classification music analysis a high resolution array technique at regional distance with recordings from the venezuela national seismic network to resolve high frequency about 0 4 hz aspects of the earthquake process our results indicate westward rupture with two subevents roughly 35 km apart in comparison a lower frequency finite source inversion with fault geometry based on new geologic and aftershock data shows two slip patches with centroids 21 km apart apparent source time functions from usarray further constrain the intersubevent time delay implying a rupture speed of 3 3 km s the tips of the slip zones coincide with subevents imaged by backprojections the different subevent locations found by backprojection and source inversion suggest spatial complementarity between high and low frequency source radiation consistent with high frequency radiation originating from rupture arrest phases at the edges of main slip areas the centroid moment tensor cmt solution and a geodetic only inversion have similar moment indicating most of the moment released is captured by geodetic observations and no additional rupture is required beyond where it is imaged in our preferred model our results demonstrate the contribution of backprojections of regional seismic array data for earthquakes down to m 7 especially when incomplete coverage of seismic and geodetic data implies large uncertainties in source inversions

the 2015 gorkha earthquake a large event illuminating the main himalayan thrust fault the 2015 gorkha earthquake sequence provides an outstanding opportunity to better characterize the geometry of the main himalayan thrust mht to overcome limitations due to unaccounted lateral heterogeneities we perform centroid moment tensor inversions in a 3 d earth model for the main shock and largest aftershocks in parallel we recompute s top and p to s receiver functions from the hi climb data set inverted centroid locations fall within a low velocity zone at 10 15 km depth and corresponding to the subhorizontal portion of the mht that ruptured during the gorkha earthquake north of the main shock hypocenter receiver functions indicate a north dipping feature that likely corresponds to the midcrustal ramp connecting the flat portion to the deep part of the mht our analysis of the main shock indicates that long period energy emanated updip of high frequency radiation sources previously inferred this frequency dependent rupture process might be explained by different factors such as fault geometry and the presence of fluids

underplating in the himalaya tibet collision zone revealed by the hi climb experiment we studied the formation of the himalayan mountain range and the tibetan plateau by investigating their lithospheric structure using an 800 kilometer long densely spaced seismic array we have constructed an image of the crust and upper mantle beneath the himalayas and the southern tibetan plateau the image reveals in a continuous fashion the main himalayan thrust fault as it extends from a shallow depth under nepal to the mid crust under southern tibet indian crust can be traced to 31 n the crust mantle interface beneath tibet is anisotropic indicating shearing during its formation the dipping mantle fabric suggests that the indian mantle is subducting in a diffuse fashion along several evolving subparallel structures

estimating the return times of great himalayan earthquakes in eastern nepal evidence from the patu and bardibas strands of the main frontal thrust the return times of large himalayan earthquakes are poorly constrained despite historical devastation of cities along the mountain range definitive links between events and specific segments of the main frontal thrust mft are not established and paleoseismological records have not documented the occurrence of several similar events at the same location in east central nepal however recently discovered primary surface ruptures of that megathrust in the a d 1255 and 1934 earthquakes are associated with flights of tectonically uplifted terraces we present here a refined longer slip history of the mft s two overlapping strands patu and bardibas thrusts in that region based on updated geomorphic neotectonic mapping of active faulting two 1 3 km long shallow seismic profiles and logging of two river cut cliffs three paleoseismological trenches and several pits with constraints from 74 detrital charcoals and 14 cosmogenic nuclide ages the amount of hanging wall uplift on the patu thrust since 3650 450 years requires three more events than the two aforementioned the uplift rate 8 5 1 5mm yr thrust dip 25 5 n and apparent characteristic behavior imply 12 17 5m of slip per event on the bardibas thrust discrete pulses of colluvial deposition resulting from the coseismic growth of a flexural fold scarp suggest the occurrence of six or seven paleo earthquakes in the last 4500 50 years the coeval rupture of both strands during great himalayan earthquakes implies that in eastern nepal the late holocene return times of such earthquakes probably ranged between 750 140 and 870 350 years

gps measurements of present day convergence across the nepal himalaya the high elevations of the himalaya and tibet result from the continuing collision between india and asia which started more than 60 million years ago1 4 from geological and seismic studies of the slip rate of faults in asia5 it is believed that approximately one third of the present day convergence rate between india and asia 58 4mmyr 1 is responsible for the shortening uplift and moderate seismicity of the himalaya great earthquakes also occur infrequently in this region releasing in minutes the elastic strain accumulated near the boundary zone over several centuries and accounting for most of the advance of the himalaya over the plains of india the recurrence time for these great earthquakes is determined by the rate of slip of india beneath tibet which has hitherto been estimated indirectly from global plate motions6 from the slip rates of faults in asia7 8 from seismic productivity9 and from the advance of sediments on the northern ganges plain10 here we report geodetic measurements using the global positioning system gps of the rate of contraction across the himalaya which we find to be 17 52 2 mm yr 1 from the form of the deformation field we estimate the rate of slip of india beneath tibet to be 20 5 2 mmyr 1 strain sufficient to drive one or more great himalayan earthquakes with slip similar to that accompanying the magnitude 8 1 bihar nepal 1934 earthquake may currently be available in western nepal

lower edge of locked main himalayan thrust unzipped by the 2015 gorkha earthquake faults are unlocked by earthquakes analysis of seismic data from the 2015 nepal earthquake shows that only part of the main himalayan thrust fault was unzipped by the quake leaving much of the fault locked and ready to slip in a future event

a window into the complexity of the dynamic rupture of the 2011 mw 9 tohoku oki earthquake the 2011 mw 9 tohoku oki earthquake recorded by over 1000 near field stations and multiple large aperture arrays is by far the best recorded earthquake in the history of seismology and provides unique opportunities to address fundamental issues in earthquake source dynamics here we conduct a high resolution array analysis based on recordings from the usarray and the european network the mutually consistent results from both arrays reveal rupture complexity with unprecedented resolution involving phases of diverse rupture speed and intermittent high frequency bursts within slow speed phases which suggests spatially heterogeneous material properties the earthquake initially propagates down dip with a slow initiation phase followed by sustained propagation at speeds of 3 km s the rupture then slows down to 1 5 km s for 60 seconds a rich sequence of bursts is generated along the down dip rim of this slow and roughly circular rupture front before the end of the slow phase an extremely fast rupture front detaches at about 5 km s towards the north finally a rupture front propagates towards the south running at about 2 5 km s for over 100 km key features of the rupture process are confirmed by the strong motion data recorded by k net and kik net the energetic high frequency radiation episodes within a slow rupture phase suggests a patchy image of the brittle ductile transition zone composed of discrete brittle asperities within a ductile matrix the high frequency is generated mainly at the down dip edge of the principal slip regions constrained by geodesy suggesting a variation along dip of the mechanical properties of the mega thrust fault or their spatial heterogeneity that affects rise time

elastodynamic analysis for slow tectonic loading with spontaneous rupture episodes on faults with rate and state dependent friction we present an efficient and rigorous numerical procedure for calculating the elastodynamic response of a fault subjected to slow tectonic loading processes of long duration within which there are episodes of rapid earthquake failure this is done for a general class of rate and state dependent friction laws with positive direct velocity effect the algorithm allows us to treat accurately within a single computational procedure loading intervals of thousands of years and to calculate for each earthquake episode initially aseismic accelerating slip prior to dynamic rupture the rupture propagation itself rapid post seismic deformation which follows and also ongoing creep slippage throughout the loading period in velocity strengthening fault regions the methodology is presented using the two dimensional 2 d antiplane spectral formulation and can be readily extended to the 2 d in plane and 3 d spectral formulations and with certain modifications to the space time boundary integral formulations as well as to their discretized development using finite difference or finite element methods the methodology can be used to address a number of important issues such as fault operation under low overall stress interaction of dynamic rupture propagation with pore pressure development patterns of rupture propagation in events nucleated naturally as a part of a sequence the earthquake nucleation process earthquake sequences on faults with heterogeneous frictional properties and or normal stress and others the procedure is illustrated for a 2 d crustal strike slip fault model with depth variable properties for lower values of the state evolution distance of the friction law small events appear the nucleation phases of the small and large events are very similar suggesting that the size of an event is determined by the conditions on the fault segments the event is propagating into rather than by the nucleation process itself we demonstrate the importance of incorporating slow tectonic loading with elastodynamics by evaluating two simplified approaches one with the slow tectonic loading but no wave effects and the other with all dynamic effects included but much higher loading rate

medieval pulse of great earthquakes in the central himalaya viewing past activities on the frontal thrust the himalaya has experienced three great earthquakes during the last century1934 nepal bihar 1950 upper assam and arguably the 1905 kangra focus here is on the central himalayan segment between the 1905 and the 1934 ruptures where previous studies have identified a great earthquake between thirteenth and sixteenth centuries historical data suggest damaging earthquakes in a d 1255 1344 1505 1803 and 1833 although their sources and magnitudes remain debated we present new evidence for a great earthquake from a trench across the base of a 13m high scarp near ramnagar at the himalayan frontal thrust the section exposed four south verging fault strands and a backthrust offsetting a broad spectrum of lithounits including colluvial deposits age data suggest that the last great earthquake in the central himalaya most likely occurred between a d 1259 and 1433 while evidence for this rupture is unmistakable the stratigraphic clues imply an earlier event which can most tentatively be placed between a d 1050 and 1250 the postulated existence of this earlier event however requires further validation if the two earthquake scenario is realistic then the successive ruptures may have occurred in close intervals and were sourced on adjacent segments that overlapped at the trench site rupture s identified in the trench closely correlate with two damaging earthquakes of 1255 and 1344 reported from nepal the present study suggests that the frontal thrust in central himalaya may have remained seismically inactive during the last similar to 700years considering this long elapsed time a great earthquake may be due in the region

surface ruptures of large himalayan earthquakes in western nepal evidence along a reactivated strand of the main boundary thrust the chronology of the seismic ruptures along the active faults of western nepal remains almost unconstrained despite their high seismogenic potential we present here a slip history of one of these structures a 120 km long reactivated segment of the main boundary thrust named the surkhet gorahi fault this slip history is based on geomorphologic and neotectonic mapping of active faults deduced from the analysis of a high resolution total station digital elevation model and 15 detrital charcoals radiocarbon ages constraining the age of deposition or abandonment of 4 alluvial terraces of the bheri river in botechaur our results show that the last two earthquakes occurred on this fault after 1860 and 640 bp respectively and accommodated slip greater than 8 m each a value corresponding to the incremental vertical offset of the terraces such events released a significant part of the slip deficit accumulated on the main himalayan thrust fault however given the geometry of this fault system as well as the date of occurrence of the last events the ruptures could be associated with major earthquakes also rupturing the main frontal thrust such as the great 1505 earthquake

magnitude calibration of north indian earthquakes this article is concerned primarily with the evaluation of the size and location of northern indian and southern tibetan earthquakes during the last 200 yr it draws attention to the problems of assessing intensity of early and more recent earthquakes in a built environment which is different from that for which the intensity scale has been constructed and to the way in which isoseismals are drawn through a re evaluation of intensities and a reassessment of isoseismals a formula for the estimation of surface wave magnitude using isoseismal radii is derived this formula is used to estimate the surface wave magnitudes of 16 earthquakes that occurred in the region between 1803 and 1900 this study shows that it is possible to calculate accurate surface wave magnitudes for earthquakes that occurred before the advent of the scale and that there is no need to resort to empirical formulae for the assessment of the size and seismic moment release of pre 20th century earthquakes also derived are formulae for the conversion of ms to m0 in total locations surface wave magnitudes and m0 estimates are presented for 43 important events that occurred in the region between 1803 and 1974 eight of which were in the lower crust or were subcrustal we find that the m0 ms scaling for india yields smaller ms than the global relation and that the methodology used can help to evaluate more realistic slip rates as well as to address other issues related to earthquake hazard in northern india

a bayesian method for incorporating self similarity into earthquake slip inversions distributions of coseismic slip help illuminate many properties of earthquakes including fault geometry stress changes frictional properties and potential future hazard slip inversions take observations and calculate slip at depth but there are a number of commonly adopted assumptions such as minimizing the second spatial derivative of slip the laplacian that have little physical basis and potentially bias the result in light of growing evidence that fault slip shows fractal properties we suggest that this information should be incorporated into slip inversions as a form of regularization instead of laplacian smoothing we have developed a bayesian approach to efficiently solve for slip incorporating von karman regularization in synthetic tests our approach retrieves fractal slip better than laplacian regularization as expected but even performs comparably or better when the input slip is not fractal we apply this to the 2014 mw 6 0 napa valley earthquake on a two segment fault using interferometric synthetic aperture radar insar and global positioning system gps data we find the von karman and laplacian inversions give similar slip magnitude but in different locations and the von karman solution has much tighter confidence bounds on slip than the laplacian solution differences in earthquake slip due to the regularization technique could have important implications for the interpretation and modeling of stress changes on the causative and neighboring faults we therefore recommend that choice of regularization method should be routinely made explicit and justified and that von karman regularization is a better default than laplacian