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

the importance of earthquake induced landslides to long term slope erosion and slope failure hazards in seismically active regions abstract this paper describes a general method for determining the amount of earthquake induced landsliding that occurs in a seismically active region over time this determination can be used as a quantitative measure of the long term hazard from seismically triggered landslides as well as a measure of the importance of this process to regional slope erosion rates and landscape evolution the method uses data from historical earthquakes to relate total volume of landslide material dislodged by an earthquake to the magnitude m and seismic moment m0 of the earthquake from worldwide data a linear regression relation between landslide volume v and m0 is determined as v m0 1018 9 0 13 where v is measured in m3 and m0 is in dyn cm to determine the amount of earthquake generated landsliding over time this relation is combined with data on seismic moment release for a particular region which may be derived from either earthquake history or fault slip data the form of the m0 v relation allows the rate of production of earthquake induced landslides over time to be determined from total rate of seismic moment release without regard to the distribution of individual events thus simplifying and generalizing the determination application of the method to twelve seismically active regions with areas ranging from 13 275 to 2 308 000 km2 shows that erosion rates from earthquake induced landslides vary significantly from region to region of the regions studied the highest rates were determined for the island of hawaii new zealand western new guinea and the san francisco bay region of california significantly lower rates were determined for iran tibet the sierra nevada great basin region of california and central japan for the time period from 715 ad to the present intermediate rates were determined for peru southern california onshore california turkey and central japan for the time period from 1586 ad to the present to determine the relative long term importance of seismically triggered landslides these erosion rates are compared to erosion rates calculated for other slope processes and to rates calculated from fluvial sediment discharge comparisons with other slope processes indicate that earthquake induced landslides are the predominant agents of slope erosion on the island of hawaii in the san francisco bay region and in western new guinea for hawaii the san francisco bay region and sierra nevada great basin region of california the erosion rates calculated for earthquake induced landslides also exceed the regional erosion rates calculated from fluvial sediment discharge

orogeny and orography the effects of erosion on the structure of mountain belts a numerical model of the coupled processes of tectonic deformation and surface erosion in convergent orogens is developed to investigate the nature of the interaction between these processes crustal deformation is calculated by a two dimensional finite element model of deformation in response to subduction and accretion of continental crust erosion operates on the uplifted surface of this model through fluvial incision which is taken to be proportional to stream power the relative importance of the tectonic and erosion processes is given by a dimensionless erosion number relating convergence velocity rock erodibility and precipitation rate this number determines the time required for a system to reach steady state and the final topographic shape and size of a mountain belt fundamental characteristics of the model orogens include asymmetric topography with shallower slopes facing the subducting plate and an asymmetric pattern of exhumation with the deepest levels of exhumation opposite to subduction these characteristics are modified when the regional climate exhibits a dominant wind direction and orographically enhanced precipitation on one side of the mountain belt the two possible cases are dominant wind in the direction of motion of the subducting plate and dominant wind direction in the opposite direction of the subducting plate velocity models of the former case predict a broad zone of exhumation with maximum exhumation in the orogen interior models of the latter case predict a focused zone of exhumation at the margin of the orogen and at high erosion number a reversal in the topographic asymmetry natural examples of these two cases are presented the southern alps of new zealand exhibits the climate and exhumation asymmetry characteristic of wind in the direction opposite to motion of the subducting plate the asymmetry of topography suggests that erosion is not efficient enough to have reversed the topographic asymmetry the contrasting example of dominant wind in the direction of subduction motion is provided by the olympic mountains of washington state in this case exhumation of deep levels of the cascadia accretionary wedge shows a broad domal pattern consistent with the observed orographic precipitation

earthquakes as multiscale dynamic ruptures with heterogeneous fracture surface energy 1 we propose a model of the wide scale growth of dynamic rupture during an earthquake based on our multiscale simulation of a planar crack in a three dimensional homogeneous elastic space a simple slip weakening law governs the fracture friction processes and its characteristic parameters slip weakening distance and fracture surface energy have multiscale heterogeneous distributions we consider a set of randomly distributed circular patches whose diameter is proportional to the fracture surface energy each patch represents an asperity between irregular fault surfaces and the size number relation of the patches obeys power law statistics we assess rupture propagation from a small instability using a boundary integral equation method with a renormalization technique although most events stop shortly after their initiation some grow triggering neighboring patches of similar size small and large events show statistically self similar properties of rupture growth and stop spontaneously without requiring a special stopping mechanism the rupture velocity locally exceeds the shear wave speed but globally remains subshear speed due to the increase of the average fracture energy as the rupture grows the relation between size and frequency of events is a power law which is explained by the triggering probability between patches as a consequence of statistically self similar random triggering growth we observe a distinct main phase in seismic waves similar to those of natural earthquakes but we cannot estimate the final size of the event from the initial part of the seismic waves if this is true for the real earthquakes predicting the size of a future earthquake would be quite difficult

incipient strain partitioning in a slate belt evidence from the early variscan monts d arree slate belt brittany france abstract partitioning of strain is a fundamental process during mountain building it commonly causes a compartmentalisation of a bulk regional strain into deformational domains with contrasting strain characteristics and largely oriented parallel to the orogenic grain the monts d arree slate belt brittany france offers an opportunity to study strain partitioning in a slate belt deformed in an overall transpressional regime the slate belt consists of highly deformed low grade siliciclastic metasediments of upper silurian to lower devonian age the deformation occurred during an early variscan nappe stacking event bretonian phase an extensive structural analysis has demonstrated that the slate belt reflects the initial stages of strain partitioning the slate belt primarily reflects coaxial contraction dominated deformation it resulted in nw verging folding and a pervasive cleavage development giving rise to a pronounced mechanical anisotropy during the later stages of deformation incipient strain partitioning lead to the development of punctuated strain heterogeneities consistently reflecting dextral belt parallel strike slip strain these structures are not organised in networks or domains incipient strain partitioning in the monts d arree slate belt did not reach the stage of compartmentalisation by the formation of an interlinked discontinuity network or wrench dominated deformational domains

three dimensional numerical simulations of crustal systems undergoing orogeny and subjected to surface processes as several modeling studies indicate the structural expression and dynamic behavior of orogenic mountain belts are dictated not only by their rheological properties or by far field tectonic motion but also by the efficiency of erosion and sedimentation acting on its surface until recently numerical investigations have been mainly limited to 2 d studies because of the high computational cost required by 3 d models here we have efficiently coupled the landscape evolution model cascade with the 3 d thermomechanically coupled tectonics code fantom details of the coupling algorithms between both codes are given we present results of numerical experiments designed to study the response of viscous plastic crustal materials subjected to convergence and to surface processes including both erosion and sedimentation in particular we focus on the equilibration of both the tectonic structures and on the surface morphology of the orogen we show that increasing the efficiency of fluvial erosion increases the frontal thrust angle which in turn decreases the width of the orogen in addition the maximum summit elevation of the orogen during transient evolution is significantly higher in those models showcasing surface processes than those that do not this illustrates the strong coupling between tectonics and surface processes we also demonstrate that an along strike gradient of erosion efficiency can have a major impact upon the landscape morphology and the tectonic structure and deformation of the orogen in both the across strike and along strike directions overall our results suggest that surface processes by enhancing localization of deformation can act as a positive forcing to topographic building

the role of sediment in controlling steady state bedrock channel slope implications of the saltation abrasion incision model abstract the saltation abrasion model is a mechanistic model for river incision into bedrock by saltating bedload which we have previously derived and used experimental data to constrain all parameter values here we develop a method for applying the saltation abrasion model at a landscape scale and use the model as a reference for evaluating the behavior of a wide range of alternative incision models in order to consider the implications of the saltation abrasion model as well as other models for predicting topographic steady state channel slope to determine the single valued discharge that best represents the effects of the full discharge distribution in transporting sediment and wearing bedrock we assume all runoff can be partitioned between a low flow and a high flow discharge in which all bedload sediment transport occurs during high flow we then use the gauged discharge record and measurements of channel characteristics at a reference field site and find that the optimum discharge has a moderate magnitude and frequency due to the constraints of the threshold of grain motion and bed alluviation by high relative sediment supply incision models can be classified according to which of the effects of sediment on bedrock incision are accounted for using the predictions of the saltation abrasion model as a reference we find that the threshold of motion is the most important effect that should be represented explicitly followed in order of decreasing importance by the cover effect the tools effect and the threshold of suspension effect models that lack the threshold of motion over predict incision rate for low shear stresses and under predict the steady state channel slope for low to moderate rock uplift rates and rock strengths models that lack the cover effect over predict incision rate for high sediment supply rates and fail to represent the degree of freedom in slope adjustment provided by partial bed coverage models that lack the tools effect over predict incision rate for low sediment supply rates and do not allow for the possibility that incision rate can decline for increases in shear stress above a peak value overall the saltation abrasion model predicts that steady state channel slope is most sensitive to changes in grain size such that the effect of variations in rock uplift rate and rock strength may affect slope indirectly through their possible but as yet poorly understood influence on the size distribution of sediments delivered to channel networks by hillslopes

late holocene paleoseismicity of the southern part of the chelungpu fault in central taiwan evidence from the chushan excavation site the geomorphic expression of the frontal western foothills in central taiwan is usually defined by a late holocene scarp that ranges from tens to hundreds of meters in height this scarp is the product of displacement on a near surface 20 30 east dipping thrust fault the chelungpu fault which ruptured during the chichi earthquake the large scarp height may correspond directly to the accumulation of successive surface ruptures the chushan site is located on the southern part of this fault where the chichi earthquake produced a scarp with a 1 7 m vertical offset for the total vertical separation based on core boring estimates the vertical displace ment on both sides is 7 m along a 24 east dipping thrust fault the results from our paleoseismic analysis indicate that five large paleoearthquake events have caused the large offsets during the past 2 ka the radiocarbon age constraints of the paleoearth quakes suggest a clustering of 540 790 cal yr b p e2 710 950 cal yr b p e3 1380 1700 cal yr b p e4 1710 1930 cal yr b p e5 and the 1999 chichi earth quake events e3 and e4 have not been reported in previous studies and we did not observe event e1 300 430 cal yr b p at the site based on displacement and fault segmentation from the geologic features we argue that the two new events may have occurred along the northern part of the chelungpu fault the vertical slip rate is estimated to be at least 3 9 0 2 mm yr over the past 2 ka which is similar to the long term estimation through a calculation of late pleistocene holocene terrace el evations on the hanging wall

insar detects increase in surface subsidence caused by an arctic tundra fire wildfire is a major disturbance in the arctic tundra and boreal forests having a significant impact on soil hydrology carbon cycling and permafrost dynamics this study explores the use of the microwave interferometric synthetic aperture radar insar technique to map and quantify ground surface subsidence caused by the anaktuvuk river fire on the north slope of alaska we detected an increase of up to 8 cm of thaw season ground subsidence after the fire which is due to a combination of thickened active layer and permafrost thaw subsidence our results illustrate the effectiveness and potential of using insar to quantify fire impacts on the arctic tundra especially in regions underlain by ice rich permafrost our study also suggests that surface subsidence is a more comprehensive indicator of fire impacts on ice rich permafrost terrain than changes in active layer thickness alone

the growth of geological structures by repeated earthquakes 1 conceptual framework in many places earthquakes with similar characteristics have been shown to recur if this is common then relatively small deformations associated with individual earthquake cycles should accumulate over time to create geological structures following this paradigm we show that existing models developed to describe leveling line changes associated with the seismic cycle can be adapted to explain geological features associated with a fault in these models an elastic layer containing the fault overlies a viscous half space with a different density fault motion associated with an earthquake results in immediate deformation followed by a long period of readjustment as stresses relax in the viscous layer and isostatic equilibrium is restored deformation is also caused as a result of the loading and unloading due to sediment deposition and erosion in this paper the parameters that control the growth of dip slip structures are identified we find that the flexural rigidity of the crust or the apparent elastic thickness provides the main control of the width of a structure the loading due to erosion and deposition of sediment determines the ratio of uplift to subsidence between the two sides of the fault the flexure due to sediment load is much more important in this respect than whether the fault is normal or reverse in character we find that in general real structures are associated with apparent elastic thicknesses of 4 km or less and thus with very low flexural rigidities

lifespan of mountain ranges scaled by feedbacks between landsliding and erosion by rivers an important challenge in geomorphology is the reconciliation of the high fluvial incision rates observed in tectonically active mountain ranges with the long term preservation of significant mountain range relief in ancient tectonically inactive orogenic belts river bedrock erosion and sediment transport are widely recognized to be the principal controls on the lifespan of mountain ranges but the factors controlling the rate of erosion and the reasons why they seem to vary significantly as a function of tectonic activity remain controversial here we use computational simulations to show that the key to understanding variations in the rate of erosion between tectonically active and inactive mountain ranges may relate to a bidirectional coupling between bedrock river incision and landslides whereas fluvial incision steepens surrounding hillslopes and increases landslide frequency landsliding affects fluvial erosion rates in two fundamentally distinct ways on the one hand large landslides overwhelm the river transport capacity and cause upstream build up of sediment that protects the river bed from further erosion on the other hand in delivering abrasive agents to the streams landslides help accelerate fluvial erosion our models illustrate how this coupling has fundamentally different implications for rates of fluvial incision in active and inactive mountain ranges the coupling therefore provides a plausible physical explanation for the preservation of significant mountain range relief in old orogenic belts up to several hundred million years after tectonic activity has effectively ceased

new constraints on sediment flux dependent river incision implications for extracting tectonic signals from river profiles we present new field data from rivers draining across active normal faults that incise across the same lithology at the fault have been subjected to similar climatic regimes and tectonic settings and were perturbed by a well documented increase in fault slip rate ca 1 ma in spite of these similarities the rivers exhibit markedly different long profiles and patterns of catchment incision we use channel slope and hydraulic geometry data for each river to calculate bed shear stresses b and show that there is no simple relationship between peak b and the relative uplift rates across the faults u which differ by a factor of four the long term average sediment supply to each channel q s estimated from time averaged catchment erosion rates can explain the b versus u data if bedload modulates bedrock incision rate e in a strongly nonlinear way together these field data allow us for the first time to evaluate theoretical predictions of the role of sediment on river profile evolution and to quantify the magnitude of the effect in natural systems

prolonged seismically induced erosion and the mass balance of a large earthquake large earthquakes deform the earth s surface and drive topographic growth in the frontal zones of mountain belts they also induce widespread mass wasting reducing relief the sum of these two opposing effects is unknown using a time series of landslide maps and suspended sediment transport data we show that the mw7 6 chi chi earthquake in taiwan was followed by a period of enhanced mass wasting and fluvial sediment evacuation peaking at more than five times the background rate and returning progressively to pre earthquake levels in about six years therefore it is now possible to calculate the mass balance and topographic effect of the earthquake the choshui river has removed sediment representing more than 30 of the added rock mass from the epicentral area this has resulted in a reduction of surface uplift by up to 0 25 m or 35 of local elevation change and a reduction of the area where the chi chi earthquake has built topography for other large earthquakes erosion may evolve in similar predictable ways reducing the efficiency of mountain building in fold and thrust belts and the topographic expression of seismogenic faults prolonging the risk of triggered processes and impeding economic regeneration of epicentral areas

a more precise chronology of earthquakes produced by the san andreas fault in southern california improved methods of radiocarbon analysis have enabled us to date more precisely the earthquake ruptures of the san andreas fault that are recorded in the sediments at pallett creek previous dates of these events had 95 confidence errors of 50 100 calendar years new error limits are less than 23 calendar years for all but two of the dated events this greater precision is due to larger sample size longer counting time lower background noise levels more precise conversion of radiocarbon ages to calendric dates and better stratigraphic constraints and statistical techniques the new date ranges with one exception fall within the broader ranges estimated previously but our estimate of the average interval between the latest 10 episodes of faulting is now about 132 years variability about the mean interval is much greater than was suspected previously five of the nine intervals are shorter than a century three of the remaining four intervals are about two to three centuries long despite the wide range of these intervals a pattern in the occurrence of large earthquakes at pallett creek is apparent in the new data the past 10 earthquakes occur in four clusters each of which consists of two or three events earthquakes within the clusters are separated by periods of several decades but the clusters are separated by dormant periods of two to three centuries this pattern may reflect important mechanical aspects of the fault s behavior if this pattern continues into the future the current period of dormancy will probably be greater than two centuries this would mean that the section of the fault represented by the pallett creek site is currently in the middle of one of its longer periods of repose between clusters and sections of the fault farther to the southeast are much more likely to produce the next great earthquake in california the greater precision of dates now available for large earthquakes recorded at the pallett creek site enables speculative correlation of events between paleoseismic sites along the southern half of the san andreas fault a history of great earthquakes with overlapping rupture zones along the mojave section of the fault remains one of the more attractive possibilities

mountain erosion over 10 yr 10 k y and 10 m y time scales we used cosmogenic 10 be to measure erosion rates over 10 k y time scales at 32 idaho mountain catchments ranging from small experimental watersheds 0 2 km 2 t o large river basins 35 000 km 2 these long term sediment yields are on average 17 times higher than stream sediment fluxes measured over 10 84 yr but are consistent with 10 m y erosion rates measured by apatite fission tracks our results imply that conventional sediment yield measurements even those made over decades can greatly underestimate long term average rates of sediment delivery and thus overestimate the life spans of engineered reservoirs our observations also suggest that sediment delivery from mountainous terrain is extremely episodic sporadically subjecting mountain stream ecosystems to extensive disturbance

effects of normal stress perturbations on the frictional properties of simulated faults 1 we report on laboratory experiments to investigate the frictional response of creeping faults to sudden changes in normal stress experiments were conducted on layers of quartz powder bare surfaces of westerly granite and layers of a 50 50 mixture of quartz powder and smectite clay powder the tests were carried out at room temperature and controlled humidity using a servo controlled double direct shear configuration normal stress perturbations corresponding to loading and unloading of tectonic fault zones were applied during steady sliding at constant loading rate from 3 to 1000 m s shear strain rates of 1 5 10 3 to 0 5 s 1 sudden changes in normal stress resulted in a linear elastic response of shear stress followed by a transient evolution of friction over a characteristic displacement the transient inelastic response is quantified as ln 0 where is the transient change in shear stress following a step change from initial normal stress 0 to final normal stress we find that is independent of sliding velocity and varies with ambient relative humidity and shear loading history for unloading we document a transition from stable to unstable behavior as a function of net slip in the range 3 to 30 mm shear strains of 1 5 to 15 increased humidity led to higher values of for pure quartz gouge but smaller for the quartz clay gouge the effects of shear displacement and humidity are discussed in the context of particle characteristics and gouge fabric development the extended rate and state dependent friction laws using one state variable and the ruina evolution law with normal stress variation describe our observations

earthquakes on dipping faults the effects of broken symmetry dynamic simulations of earthquakes on dipping faults show asymmetric near source ground motion caused by the asymmetric geometry of such faults the ground motion from a thrust or reverse fault is larger than that of a normal fault by a factor of 2 or more given identical initial stress magnitudes the motion of the hanging wall is larger than that of the footwall in both thrust reverse and normal earthquakes the asymmetry between normal and thrust reverse faults results from time dependent normal stress caused by the interaction of the earthquake generated stress field with earth s free surface the asymmetry between hanging wall and footwall results from the asymmetric mass and geometry on the two sides of the fault

brittle frictional mountain building 1 deformation and mechanical energy budget an active fold and thrust belt is analogous to a wedge of soil or snow that forms in front of a moving bulldozer such wedges exhibit a critical taper and a regional state of stress that is everywhere on the verge of coulomb failure the width of such a critically tapered fold and thrust belt does not depend on its brittle strength or frictional properties but rather on the accretionary influx rate of fresh material at its toe and on the rate of erosion a steady state fold and thrust belt is one in which the accretionary influx is balanced by the erosive efflux rocks are accreted at the toe and then horizontally shortened as they are transported toward the rear those that enter lower in the accreted section are more deeply buried before being uplifted by erosion mass balance and isotropy constrain the kinematics of this large scale deformation enabling us to infer the trajectories residence times and stress strain histories of rocks incorporated into eroding fold and thrust belts a typical rock resides in the steady state taiwan wedge for 2 3 m y before it is uplifted and eroded during its motion through the wedge it experiences strain rates in the range 10 13 to 10 14 s 1 the mechanical energy budget of brittle frictional mountain building is described by the equation where is the rate of work performed on the base and front of the fold and thrust belt by the subducting plate is the rate at which energy is dissipated against friction on the decollement fault is the rate at which energy is dissipated by internal frictional processes within the deforming brittle wedge and is the rate of work performed against gravitational body forces in a reference frame attached to the overriding plate the total mechanical power being supplied by the subducting eurasian plate to the active fold and thrust belt in taiwan is slightly over 3 gw approximately 60 of this work of steady state mountain building is being dissipated against friction on the decollement fault and about another 25 is being dissipated against internal friction this leaves only 15 or roughly half a gigawatt available to do useful work against gravity in general fold and thrust belts with moderate pore fluid pressures are dominated by work done against friction on the decollement fault however those with nearly lithostatic pore fluid pressures may be dominated by work done against gravity internal frictional dissipation is always less than basal frictional dissipation as it is in taiwan an alternative and equivalent description of the mechanical energy balance of a steady state fold and thrust belt is provided by the equation in this version the quantity on the left is the rate at which work is performed on the back of the wedge by the overriding plate and is the rate of work performed against gravity in a reference frame attached to the subducting plate the latter quantity is always positive for any critically tapered fold and thrust belt whose decollement fault dips toward its rear in contradiction to the central premise of the gravity spreading theory of fold and thrust tectonics

variability of earthquake nucleation in continuum models of rate and state faults and implications for aftershock rates using two continuum models of rate and state faults one with a weaker patch and the other with rheological transition from steady state velocity weakening to velocity strengthening friction we simulate several scenarios of spontaneous earthquake nucleation plausible for natural faults investigate their response to static shear stress steps and infer the corresponding aftershock rates overall nucleation processes at weaker patches behave similarly to theories based on spring slider models with some notable deviations in particular nucleation and aftershock rates are affected by normal stress heterogeneity in the nucleation zone nucleation processes at rheological transitions behave differently producing complex slip velocity histories nonmonotonic responses to static stress changes and aftershock rates with pronounced peaks and seismic quiescence for such processes positive stress steps sometimes delay nucleation of seismic events by inducing aseismic transients that relieve stress and postpone seismic slip superposition of the complex aftershock response for spatially variable stress changes results in omori s law for a period of time followed by seismic quiescence such behavior was observed at the base of the seismogenic zone near the 1984 morgan hill earthquake we show that the computed aftershock rates are linked to unperturbed slip velocity evolution in the nucleation zone and construct simplified analytical scenarios that explain some features of the response the qualitative differences that we find between the two nucleation models indicate that aftershock response of rate and state faults to static stress changes would depend on the conditions under which nucleation occurs on natural faults and may be different from predictions based on spring slider models

how rivers react to large earthquakes evidence from central taiwan earthquakes and bedrock river incision are fundamental processes in the evolution of tectonically active landscapes yet little work has focused on understanding how a bedrock river responds to a single large earthquake data from the 1999 m w 7 6 chi chi earthquake in central taiwan show dramatic differences in river response that depend on the proximity to the rupture zone near the fault scarp vertical ground deformation intensifi es river incision on a time scale of years to decades while distal to the fault landslides induced by the earthquake mantle bedrock on the river bed with sediment impeding incision for decades to centuries as the material is evacuated this surprising spatial and temporal variability in bedrock incision caused by earthquakes has implications for the paleoseismic interpretation of bedrock terraces and for the evolution of tectonically active landscapes

three dimensional boundary integral modeling of spontaneous earthquake sequences and aseismic slip 1 fault processes involve complex patterns of seismic events and aseismic slip this work develops a three dimensional 3 d methodology for simulating long term history of spontaneous seismic and aseismic slip on a vertical planar strike slip fault subjected to slow tectonic loading our approach reproduces all stages of earthquake cycles from accelerating slip before dynamic instability to rapid dynamic propagation of earthquake rupture to postseismic slip and to interseismic creep including aseismic transients we use the developed 3 d methodology to study interaction of fault slip with a small patch of higher normal stress over long term slip history for uniform initial prestress dynamic rupture is significantly affected by the stronger patch in the first simulated event but not in subsequent ones the change in behavior is due to redistribution of shear stress by prior slip which demonstrates that distributions of fault strength and stress are related and illustrates the importance of simulating long slip histories even in studies of dynamic rupture despite no long term effect on dynamic rupture the small patch of higher normal stress influences nucleation processes and hence long term slip patterns in the model comparison of the fully dynamic simulations and a widely used quasi dynamic approach shows that the quasi dynamic approach modifies long term slip patterns in addition to resulting in much smaller slip velocities and rupture speeds during dynamic events we show that the response of quasi dynamic formulations with reduced radiation damping terms can be scaled to match the results of the standard quasi dynamic formulation and hence cannot improve the comparison