seismic evidence for an earthquake nucleation phase near source observations show that earthquakes initiate with a distinctive seismic nucleation phase that is characterized by a low rate of moment release relative to the rest of the event this phase was observed for the 30 earthquakes having moment magnitudes 2 6 to 8 1 and the size and duration of this phase scale with the eventual size of the earthquake during the nucleation phase moment release was irregular and appears to have been confined to a limited region of the fault it was characteristically followed by quadratic growth in the moment rate as rupture began to propagate away from the nucleation zone these observations suggest that the nucleation process exerts a strong influence on the size of the eventual earthquake

scaling of the critical slip distance for seismic faulting with shear strain in fault zones theoretical and experimentally based laws for seismic faulting contain a critical slip distance1 5 dc which is the slip over which strength breaks down during earthquake nucleation on an earthquake generating fault this distance plays a key role in determining the rupture nucleation dimension6 the amount of premonitory and post seismic slip7 10 and the maximum seismic ground acceleration1 11 in laboratory friction experiments d c has been related to the size of surface contact junctions2 5 12 thus the discrepancy between laboratory measurements of dc 10 5m and values obtained from modelling earthquakes 10 2m has been attributed to differences in roughness between laboratory surfaces and natural faults5 this interpretation predicts a dependence of dc on the particle size of fault gouge2 breccia and wear material but not on shear strain here we present experimental results showing that dc scales with shear strain in simulated fault gouge our data suggest a new physical interpretation for the critical slip distance in which dc is controlled by the thickness of the zone of localized shear strain as gouge zones of mature faults are commonly 102 103 m thick12 17 whereas laboratory gouge layers are 1 10 mm thick our data offer an alternative interpretation of the discrepancy between laboratory and field based estimates of dc

creep dike intrusion and magma chamber deflation model for the 2000 miyake eruption and the izu islands earthquakes 1 analysis of global positioning system data shows shrinkage of miyake island and the widening between nijima and kozu islands during the period of the miyake island volcanic activity and the ensuing izu islands earthquakes in 2000 the estimated time evolution of a model consisting of a dike creeping faults and a mogi source suggests that a crack opening on miyake island occurred immediately after the start of the seismic activities on 26 june and ended within several days at the west coast of miyake island after the seabed eruption on 27 june magma migrated from miyake island to kozu and nijima islands within several days the estimated volume of intruded magma totals around 1 2 109 m3 associated with the magma intrusion between miyake and kozu islands left lateral and right lateral creep motions occurred in regions off the west coast of miyake island and near kozu island the accumulated moment energy is equivalent to an earthquake of mw 6 6 and 6 6 for right lateral and left lateral creeping faults respectively the estimated magma chamber continued deflation beneath the southwestern part of miyake island from 26 june totaling around 0 12 109 m3 in volume change in addition to the collapse volume of 0 6 109 m3 at the summit of mount oyama on miyake island the volume change on miyake island can be compensated by the migrated magma toward kozu island from the deflation source beneath miyake island the deflation speed of the magma chamber beneath miyake island decreased and increased before and after the eruption of 14 15 july and 18 august suggesting a change in balance of mass influx and draining out rate of the magma chamber

propagation of slow slip leading up to the 2011 mw 9 0 tohoku oki earthquake many large earthquakes are preceded by one or more foreshocks but it is unclear how these foreshocks relate to the nucleation process of the mainshock on the basis of an earthquake catalog created using a waveform correlation technique we identified two distinct sequences of foreshocks migrating at rates of 2 to 10 kilometers per day along the trench axis toward the epicenter of the 2011 moment magnitude mw 9 0 tohoku oki earthquake in japan the time history of quasi static slip along the plate interface based on small repeating earthquakes that were part of the migrating seismicity suggests that two sequences involved slow slip transients propagating toward the initial rupture point the second sequence which involved large slip rates may have caused substantial stress loading prompting the unstable dynamic rupture of the mainshock

characterization of nucleation during laboratory earthquakes 1 we observe the nucleation phase of in plane ruptures in the laboratory we show that the nucleation is composed of two distinct phases a quasi static and an acceleration stage followed by dynamic propagation we propose an empirical model which describes the rupture length evolution the quasi static phase is described by an exponential growth while the acceleration phase is described by an inverse power law of time the transition from quasi static to accelerating rupture is related to the critical nucleation length which scales inversely with normal stress in accordance with theoretical predictions and to a critical surfacic power which may be an intrinsic property of the interface finally we discuss these results in the frame of previous studies and propose a scaling up to natural earthquake dimensions

characteristic features of local breakdown near a crack tip in the transition zone from nucleation to unstable rupture during stick slip shear failure abstract characteristic features of local breakdown near a propagating crack tip during slip failure nucleation and its transition process to unstable rupture under mode ii conditions are experimentally studied using a rock sample with a simulated fault the size of which is large compared with the size of the breakdown zone in addition distinctive features in the nucleation zone are compared with those in the zone of steady dynamic rupture propagation it is found that the local breakdown stress drop the local dynamic stress drop the local stress increase immediately before the slip weakening process the dimensionless parameter s defined by das and aki and the shear fracture energy are functions of crack length in the transition zone from nucleation to unstable rupture i e they all increase with crack growth in the nucleation zone this is contrasted with the observations in the zone of steady dynamic rupture propagation where these physical parameters do not depend on crack length an increase in the magnitudes of these parameters with crack growth in the nucleation zone results in increasing resistance to crack extension with crack growth in the nucleation zone the reason for the increase in crack growth resistance in this zone is that a slip failure nucleus is formed at a point or zone where the local strength and the shear fracture energy are at a minimum this basic physical information can provide the key to earthquake prediction since the nucleation itself is a precursory phenomenon the ratio of the size of the breakdown zone to the crack length decreases abruptly beyond the critical crack length rupture velocity v accelerates with crack growth and reaches the shear wave velocity or super shear and v is found to be expressed empirically as a function of crack length in the nucleation zone a positive correlation is found between s and v in the nucleation zone v increases with increasing s however no such correlation was observed in the zone of steady dynamic rupture propagation where v seems to have reached the shear wave velocity of super shear for s 0 5 this is consistent with das and aki s numerical result the interrelationships between shear stress slip displacement slip velocity and acceleration near the crack tip during the breakdown in the nucleation zone are revealed and they are compared with those in the zone of dynamic rupture propagation the present experiments suggest that the relationship between shear stress and slip displacement is more fundamental during the breakdown process than the relationship between shear stress and slip velocity a model is presented to describe the breakdown process in the nucleation zone

smoothing and re roughening processes the geometric evolution of a single fault zone abstract the geometry of a fault zone exerts a major control on earthquake rupture processes and source parameters observations previously compiled from multiple faults suggest that fault surface shape evolves with displacement but the specific processes driving the evolution of fault geometry within a single fault zone are not well understood here we characterize the deformation history and geometry of an extraordinarily well exposed fault using maps of cross sectional exposures constructed with the structure from motion photogrammetric method the la quinta fault located in southern california experienced at least three phases of deformation multiple layers of ultracataclasite formed during the most recent phase crosscutting relations between the layers define the evolution of the structures and demonstrate that new layers formed successively during the deformation history wear processes such as grain plucking from one layer into a younger layer and truncation of asperities at layer edges indicate that the layers were slip zones and the contacts between them slip surfaces slip surfaces that were not reactivated or modified after they were abandoned exhibit self affine geometry preserving the fault roughness from different stages of faulting roughness varies little between surfaces except the last slip zone to form in the fault which is the smoothest this layer contains a distinct mineral assemblage indicating that the composition of the fault rock exerts a control on roughness in contrast the similar roughness of the older slip zones which have comparable mineralogy but clearly crosscut one another suggests that as the fault matured the roughness of the active slip surface stayed approximately constant wear processes affected these layers so for roughness to stay constant the roughening and smoothing effects of fault slip must have been approximately balanced these observations suggest fault surface evolution occurs by nucleation of new surfaces and wear by competing smoothing and re roughening processes

preslip and cascade processes initiating laboratory stick slip recent modeling studies have explored whether earthquakes begin with a large aseismic nucleation process or initiate dynamically from the rapid growth of a smaller instability in a cascade up process to explore such a case in the laboratory we study the initiation of dynamic rupture stick slip of a smooth saw cut fault in a 76 mm diameter cylindrical granite laboratory sample at 40 120 mpa confining pressure we use a high dynamic range recording system to directly compare the seismic waves radiated during the stick slip event to those radiated from tiny m 6 discrete seismic events commonly known as acoustic emissions aes that occur in the seconds prior to each large stick slip the seismic moments focal mechanisms locations and timing of the aes all contribute to our understanding of their mechanics and provide us with information about the stick slip nucleation process in a sequence of 10 stick slips the first few microseconds of the signals recorded from stick slip instabilities are nearly indistinguishable from those of premonitory aes in this sense it appears that each stick slip begins as an ae event that rapidly 20 s grows about 2 orders of magnitude in linear dimension and ruptures the entire 150 mm length of the simulated fault we also measure accelerating fault slip in the final seconds before stick slip we estimate that this slip is at least 98 aseismic and that it both weakens the fault and produces aes that will eventually cascade up to initiate the larger dynamic rupture

the long precursory phase of most large interplate earthquakes many earthquakes are preceded by foreshocks however the mechanisms that generate foreshocks and the reason why they occur before some earthquakes and not others are unknown here we use seismic catalogues from the best instrumented areas of the north pacific to analyse the foreshock sequences preceding all earthquakes there between 1999 and 2011 of magnitude larger than 6 5 and at depths shallower than 50 km the data set comprises 31 earthquakes at plate boundaries and 31 in plate interiors we find that there is a remarkable contrast between the foreshock sequences of interplate compared with intraplate earthquakes most large earthquakes at plate interfaces in the north pacific were preceded by accelerating seismic activity in the months to days leading up to the mainshock in contrast foreshocks are much less frequent in intraplate settings we suggest that at plate boundaries the interface between the two plates begins to slowly slip before the interface ruptures in a large earthquake this relatively long precursory phase could help mitigate earthquake risk at plate boundaries

dynamic rupture scenarios of anticipated nankai tonankai earthquakes southwest japan 1 we investigated dynamic rupture scenarios of anticipated megathrust earthquakes on the nankai tonankai subduction zone southwest japan to improve the scenario reliability the model parameters should be constrained by available data or derived from their analysis we employed the three dimensional plate interface geometry and the slip deficit rate on the interface accumulated slip deficit was used to obtain the stress drop distribution of anticipated earthquakes the estimated stress drop distribution is consistent with the seismogenic asperity locations known from the analysis of past earthquakes fault friction constitutive parameters however had to be assumed from indirect observations because they cannot be constrained directly by the data based on various geophysical observations we defined three regions where larger fracture energy is required these are the eastern edge of the tonankai area the western edge of the nankai area and the region between the tonankai and nankai areas beneath the kii peninsula such lateral heterogeneity promoted the segmented rupture along the nankai trough with predefined stress drop and constitutive parameters various rupture scenarios for tonankai and nankai asperities were obtained for different initiation locations in some cases a single segment is ruptured while in other cases all the segments are broken due to dynamic linkage at the segment boundary causing a giant earthquake the initiation location is a critical parameter that controls the rupture propagation across the segment boundary these scenarios will be extremely useful to evaluate deterministically the strong ground motions and tsunami hazards caused by the next major earthquakes in southwest japan

earthquake nucleation on rate and state faults aging and slip laws we compare 2 d quasi static earthquake nucleation on rate and state faults under both aging and slip versions of the state evolution law for both versions mature nucleation zones exhibit 2 primary regimes of growth well above and slightly above steady state corresponding respectively to larger and smaller fault weakening rates well above steady state aging law nucleation takes the form of accelerating slip on a patch of fixed length this length is proportional to b 1 and independent of a where a and b are the constitutive parameters relating changes in slip speed and state to frictional strength under the slip law the nucleation zone is smaller and continually shrinks as slip accelerates the nucleation zone is guaranteed to remain well above steady state only for values of a b that are low by laboratory standards near steady state for both laws the nucleation zone expands the propagating front remains well above steady state giving rise to a simple expression for its effective fracture energy g c this fracture energy controls the propagation style for the aging law g c increases approximately as the square of the logarithm of the velocity jump this causes the nucleation zone to undergo quasi static crack like expansion to a size asymptotically proportional to b b a 2 for the slip law g c increases only as the logarithm of the velocity jump and crack like expansion is not an option instead the nucleation zone grows as an accelerating unidirectional slip pulse under both laws the nucleation front propagates at a velocity larger than the slip speed by roughly b divided by the logarithm of the velocity jump where is the effective elastic shear modulus for this prediction to be consistent with observed propagation speeds of slow slip events in subduction zones appears to require effective normal stresses as low as 1 mpa

hydrodynamic models of boulder berm deposition the hydrodynamic environment associated with boulder berm deposition at flow expansions during high magnitude flood is considered with respect to the plan geometry of surveyed berms hydraulic computer simulation of a prototype channel expansion and froude scaled physical hydraulic models are used to demonstrate the processes of berm formation all the geometric features observed in the field and temporal scales of berm formation can be reproduced in the laboratory and compared in principle with hydraulic energy diffusers used in engineering practice the experiments indicate that berms form most readily at low froude numbers f 1 5 owing to the deposition of individual bedload clasts and the shearing of transverse bars in the shear flow zone at the edge of the flow separation zone the resultant accretionary berms minimize energy losses and ensure that high flows are transmitted efficiently through channel transitions berm formation may occur with either steady or unsteady supercritical flows f 2 0 through a variety of channel expansions where initial flow already maintains a minimal backwater effect it would appear that formation of berms is not associated with initial flow into dry expansions where no backwater depth is maintained in all cases stable berm formation occurs extremely rapidly and the depositional process is not sensibly affected by steady or unsteady sediment feed rates

rupture preparation process controlled by surface roughness on meter scale laboratory fault abstract we investigate the effect of fault surface roughness on rupture preparation characteristics using meter scale metagabbro specimens we repeatedly conducted the experiments with the same pair of rock specimens to make the fault surface rough we obtained three experimental results under the same experimental conditions 6 7 mpa of normal stress and 0 01 mm s of loading rate but at different roughness conditions smooth moderately roughened and heavily roughened during each experiment we observed many stick slip events preceded by precursory slow slip we investigated when and where slow slip initiated by using the strain gauge data processed by the kalman filter algorithm the observed rupture preparation processes on the smooth fault i e the first experiment among the three showed high repeatability of the spatiotemporal distributions of slow slip initiation local stress measurements revealed that slow slip initiated around the region where the ratio of shear to normal stress was the highest as expected from finite element method fem modeling however the exact location of slow slip initiation was where became locally minimum probably due to the frictional heterogeneity in the experiment on the moderately roughened fault some irregular events were observed though the basic characteristics of other regular events were similar to those on the smooth fault local stress data revealed that the spatiotemporal characteristics of slow slip initiation and the resulting drop for irregular events were different from those for regular ones even under similar stress conditions on the heavily roughened fault the location of slow slip initiation was not consistent with anymore because of the highly heterogeneous static friction on the fault which also decreased the repeatability of spatiotemporal distributions of slow slip initiation these results suggest that fault surface roughness strongly controls the rupture preparation process and generally increases its complexity with the degree of roughness

foreshocks and pre events associated with the nucleation of large earthquakes we propose a theoretical model which can explain the occurrence of foreshocks and pre events associated with the nucleation of large earthquakes we consider a broad weak zone with a locally strong part asperity on a fault plane in a two dimensional framework and examine the transition process from quasistatic nucleation to high speed rupture propagation through numerical simulations the dynamic rupture of the asperity occurs in three different manners as aseismic slip foreshock or pre event depending on the peak stress of the asperity in all cases the rupture of the asperity accelerates the nucleation of main rupture

extended nucleation of the 1999 mw 7 6 izmit earthquake laboratory and theoretical studies suggest that earthquakes are preceded by a phase of developing slip instability in which the fault slips slowly before accelerating to dynamic rupture we report here that one of the best recorded large earthquakes to date the 1999 moment magnitude mw 7 6 izmit turkey earthquake was preceded by a seismic signal of long duration that originated from the hypocenter the signal consisted of a succession of repetitive seismic bursts accelerating with time and increased low frequency seismic noise these observations show that the earthquake was preceded for 44 minutes by a phase of slow slip occurring at the base of the brittle crust this slip accelerated slowly initially and then rapidly accelerated in the 2 minutes preceding the earthquake

effects of physical fault properties on frictional instabilities produced on simulated faults laboratory studies of large scale simulated faults show that physical properties of the fault specifically normal stress and fault roughness significantly influence the unstable shear failure behavior of the fault in addition the experiments provide insights into important length scaling effects that are useful for assessing concepts such as critical crack length or rupture nucleation dimension stick slip shear failures have been generated along a 2 m long simulated fault in a block of sierra white granite experiments were conducted at normal stresses between 0 6 and 4 mpa with two different prepared roughnesses of the simulated fault a smooth fault with a profilometer measured roughness of 0 2 m and a rough fault with a measured roughness of 80 m high speed records of shear strains and slip velocities clearly show both the propagation of the stick slip failures and details of the local deformations and motions from which slip weakening of the fault at the onset of sliding can be resolved both the dynamic stress drops of the stick slip events and the apparent fracture energies of the events determined from the records of the fault slip weakening increase with increasing normal stress although critical slip weakening displacements are insensitive to normal stress they increase with increased fault roughness such that the apparent fracture energies are also larger on the rough fault than on the smooth fault slip velocities in these experiments are influenced by roughness and are proportional to stress drop although the onset of slip is sharply defined the termination of slip is more gradual without a clear stopping phase measured rupture velocities on the rough fault are lower than those on the smooth fault when a length scale based on rupture nucleation dimensions directly proportional to critical displacements is introduced the data suggest that at fault lengths close to the nucleation lengths rupture velocity and slip velocity increase with increasing fault length

properties of the seismic nucleation phase abstract near source observations show that earthquakes begin abruptly at the p wave arrival but that this beginning is weak with a low moment rate relative to the rest of the main shock we term this initial phase of low moment rate the seismic nucleation phase we have observed the seismic nucleation phase for a set of 48 earthquakes ranging in magnitude from 1 1 8 1 the size and duration of the seismic nucleation phase scale with the total seismic moment of the earthquake suggesting that the process responsible for the seismic nucleation phase carries information about the eventual size of the earthquake the seismic nucleation phase is characteristically followed by quadratic growth in the moment rate consistent with self similar rupture at constant stress drop in this paper we quantify the properties of the seismic nucleation phase and offer several possible explanations for it

hydrological response to earthquakes in the haibara well central japan i groundwater level changes revealed using state space decomposition of atmospheric pressure rainfall and tidal responses summary for the groundwater level observed at the haibara well shizuoka prefecture central japan time series analysis using state space modelling is applied to extract hydrological anomalies related to earthquakes this method can decompose observed groundwater level time series into five components atmospheric pressure tidal and precipitation responses observation noise and residual water level the decomposed responses to atmospheric pressure and precipitation are independently determined and are consistent with the expected response to surface loading in the groundwater level at the haibara well 28 coseismic changes can be discerned during the period from 1981 april to 1997 december there is a threshold in the relationship between earthquake magnitude and the well hypocentre distance above which earthquakes cause coseismic changes in the residual water level all of the coseismic water level changes at the haibara well are decreases although 33 per cent of the estimated coseismic volumetric strain steps are contraction which would be expected to cause water level increases the coseismic changes in groundwater level are more closely proportional to the estimated ground motion than to coseismic volumetric strain steps suggesting that ground motion due to earthquakes is the major cause of the coseismic water level drops and that the contribution from static strain is rather small possible pre or inter earthquake water level changes have occurred at the haibara well and may have been caused by local aseismic crustal deformation

magmatic rifting of pangaea linked to onset of south american plate motion abstract the causes of the transition from amagmatic to magmatic rifting during continental break up are not always clear and have been often linked to the break through of a plume however stress fields recorded in the crust may offer new insights into the relationship between changes in the stress field and the onset of magmatism in this paper stress fields recorded in the crust of north america are used to test possible causes of the break up of pangaea and the transition from amagmatic to magmatic rifting finite element models reveal that the most likely scenario for the break up involves the initial northwest motion of north america at 230 ma followed by a south southeast motion of south america at 200 ma with the initiation of magmatism and finally a weakened area between north and south america sometime soon after 200 ma which resulted in north america being dominated by northwest motion once again it was also determined that plate boundary structure and orientation play a large part in the recorded stress fields and must be taken into consideration when modeling continental break ups and rifting

normal stress distribution of rough surfaces in contact we study numerically the stress distribution on the interface between two thick elastic media bounded by interfaces that include spatially correlated asperities the interface roughness is described using the self affine topography that is observed over a very wide range of scales from fractures to faults we analyse the correlation properties of the normal stress distribution when the rough surfaces have been brought into full contact the self affinity of the rough surfaces is described by a hurst exponent h we find that the normal stress field is also self affine but with a hurst exponent h 1 fluctations of the normal stress are shown to be important especially at local scales with anti persistent correlations

scale dependence of rock friction at high work rate in metre sized rock specimens rock friction starts to decrease at a much smaller work rate than in centimetre sized rock specimens thus demonstrating that rock friction is scale dependent frictional sliding at faults controls the rupture of earthquakes and much effort has gone into the study of the frictional properties of rocks in laboratory experiments this study shows that previously used scaling of results obtained in lab studies investigating the earth s crust didn t accurately take into account the effect of sample size futoshi yamashita et al used metre sized rock specimens not the centimetre sized samples commonly used and demonstrate that friction in the metre size range starts to decrease at a work rate that is an order of magnitude smaller than that of the smaller specimens the authors propose that stress concentrated areas exist in which more gouge materials are produced due to frictional slip resulting in further stress concentrations at these areas as such heterogeneity is common in nature the authors conclude that a natural fault could lose its strength faster than that expected from the properties estimated from small rock samples determination of the frictional properties of rocks is crucial for an understanding of earthquake mechanics because most earthquakes are caused by frictional sliding along faults prior studies using rotary shear apparatus1 2 3 4 5 6 7 8 9 10 11 12 13 revealed a marked decrease in frictional strength which can cause a large stress drop and strong shaking with increasing slip rate and increasing work rate the mechanical work rate per unit area equals the product of the shear stress and the slip rate however those important findings were obtained in experiments using rock specimens with dimensions of only several centimetres which are much smaller than the dimensions of a natural fault of the order of 1 000 metres here we use a large scale biaxial friction apparatus with metre sized rock specimens to investigate scale dependent rock friction the experiments show that rock friction in metre sized rock specimens starts to decrease at a work rate that is one order of magnitude smaller than that in centimetre sized rock specimens mechanical visual and material observations suggest that slip evolved stress heterogeneity on the fault accounts for the difference on the basis of these observations we propose that stress concentrated areas exist in which frictional slip produces more wear materials gouge than in areas outside resulting in further stress concentrations at these areas shear stress on the fault is primarily sustained by stress concentrated areas that undergo a high work rate so those areas should weaken rapidly and cause the macroscopic frictional strength to decrease abruptly to verify this idea we conducted numerical simulations assuming that local friction follows the frictional properties observed on centimetre sized rock specimens the simulations reproduced the macroscopic frictional properties observed on the metre sized rock specimens given that localized stress concentrations commonly occur naturally our results suggest that a natural fault may lose its strength faster than would be expected from the properties estimated from centimetre sized rock samples

a jurassic granite from southern georgia u s a silicic extension related magmatism along the southeastern coastal plain abstract ion microprobe shrimp u pb dating of single zircon grains extracted from granite cuttings from the subsurface of southeastern georgia pierce co yields a well defined crystallization age of documentclass aastex usepackage amsbsy usepackage amsfonts usepackage amssymb usepackage bm usepackage mathrsfs usepackage pifont usepackage stmaryrd usepackage textcomp usepackage portland xspace usepackage amsmath amsxtra usepackage ot2 ot1 fontenc newcommand cyr renewcommand rmdefault wncyr renewcommand sfdefault wncyss renewcommand encodingdefault ot2 normalfont selectfont declaretextfontcommand textcyr cyr pagestyle empty declaremathsizes 10 9 7 6 begin document landscape 159 pm 3 end document ma 2 confirming the presence of mesozoic silicic magmatism in the northern suwannee terrane revision of suwannee terrane subcrop age province boundaries is required by these data the granite is 20 40 m yr younger than subsurface diabases in the same r

a mechanism for preseismic steady rupture fronts observed in laboratory experiments it has been shown that the onset of frictional instability is characterized by a transition from stable quasi static rupture growth to unstable inertially controlled high speed rupture in particular slow rupture fronts propagating at a steady speed v slow of the order of 5 of the s wave speed have been observed prior to the onset of dynamic rupture in recent fault friction laboratory experiments however the precise mechanism governing this v slow stage is unknown here we reproduce this phenomenon in numerical simulations of earthquake sequences that incorporate laboratory derived rate and state friction laws our simulations show that the v slow stage originates from a stress concentration inherited from the coalescence of interseismic slow creep fronts its occurrence is limited to a narrow range of the parameter space but is found in simulations with two commonly used state variable evolution laws in the rate and state formulation the sensitivity of the speed v slow to the model parameters suggests that the propagation speed v slow reported in laboratory experiments may also be sensitive to parameters of friction and stress conditions our results imply that time and space dimensions associated with the propagation of v slow on natural faults can be as much as a few seconds and several hundred meters respectively hence the detection of such preseismic signals may be possible with near field high resolution observations

experimental observation of stick slip instability fronts summary stick slip dynamic instability is a key mechanism governing frictional processes from microscale physics to earthquake faults and landslides yet challenging questions are still open about its nucleation and propagation dynamics here we present novel observations on laboratory experimental faults where spontaneously nucleating fractures are produced describing 1 an initial quasistatic stable rupture front propagating at about 5 per cent the shear wave velocity vs accelerating to subshear and then to intersonic velocity 2 the arisal of a higher degree of complexity when the friction to prestress ratio is increased on the sliding surface the complex behaviour includes stop and go sequences irregular propagation and rerupturing episodes within short time intervals implying rapid restrengthening of the surface and the formation of self healing pulses reproducing experimentally for the first time a behaviour observed on seismic faults

earthquake nucleation on faults with rate and state dependent strength abstract dieterich j h 1992 earthquake nucleation on faults with rate and state dependent strength in t mikumo k aki m ohnaka l j ruff and p k p spudich editors earthquake source physics and earthquake precursors tectonophysics 211 115 134 faults with rate and state dependent constitutive properties reproduce a range of observed fault slip phenomena including spontaneous nucleation of slip instabilities at stresses above some critical stress level and recovery of strength following slip instability calculations with a plane strain fault model with spatially varying properties demonstrate that accelerating slip precedes instability and becomes localized to a fault patch the dimensions of the fault patch follow scaling relations for the minimum critical length for unstable fault slip the critical length is a function of normal stress loading conditions and constitutive parameters which include dc the characteristic slip distance if slip starts on a patch that exceeds the critical size the length of the rapidly accelerating zone tends to shrink to the characteristic size as the time of instability approaches solutions have been obtained for a uniform fixed patch model that are in good agreement with results from the plane strain model over a wide range of conditions above the steady state stress the logarithm of the time to instability linearly decreases as the initial stress increases because nucleation patch length and premonitory displacement are proportional to dc the moment of premonitory slip scales by d3c the scaling of dc is currently an open question unless dc for earthquake faults is significantly greater than that observed on laboratory faults premonitory strain arising from the nucleation process for earthquakes may by too small to detect using current observation methods excluding the possibility that dc in the nucleation zone controls the magnitude of the subsequent earthquake then the source dimensions of the smallest earthquakes in a region provide an upper limit for the size of the nucleation patch

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

foreshocks during the nucleation of stick slip instability 1 we report on laboratory experiments which investigate interactions between aseismic slip stress changes and seismicity on a critically stressed fault during the nucleation of stick slip instability we monitor quasi static and dynamic changes in local shear stress and fault slip with arrays of gages deployed along a simulated strike slip fault 2 m long and 0 4 m deep in a saw cut sample of sierra white granite with 14 piezoelectric sensors we simultaneously monitor seismic signals produced during the nucleation phase and subsequent dynamic rupture we observe localized aseismic fault slip in an approximately meter sized zone in the center of the fault while the ends of the fault remain locked clusters of high frequency foreshocks mw 6 5 to 5 0 can occur in this slowly slipping zone 5 50 ms prior to the initiation of dynamic rupture their occurrence appears to be dependent on the rate at which local shear stress is applied to the fault the meter sized nucleation zone is generally consistent with theoretical estimates but source radii of the foreshocks 2 to 70 mm are 1 to 2 orders of magnitude smaller than the theoretical minimum length scale over which earthquake nucleation can occur we propose that frictional stability and the transition between seismic and aseismic slip are modulated by local stressing rate and that fault sections which would typically slip aseismically may radiate seismic waves if they are rapidly stressed fault behavior of this type may provide physical insight into the mechanics of foreshocks tremor repeating earthquake sequences and a minimum earthquake source dimension

scaling of the shear rupture process from nucleation to dynamic propagation implications of geometric irregularity of the rupturing surfaces a series of systematic high resolution laboratory experiments have been performed on the nucleation of propagating slip failure on preexisting faults having different surface roughnesses to demonstrate how the size scale and duration of shear rupture nucleation are affected by geometric irregularity of the rupturing surfaces on the basis of the experimental results it has been discussed theoretically how consistently scale dependent physical quantities inherent in shear rupture are scaled the experiments led to conclusive results that the nucleation process consists of two phases phase i an initial quasi static phase and phase ii a subsequent accelerating phase and that the nucleation process is greatly affected by geometric irregularities on the rupturing surfaces in phase i the rupture grows at a slow steady speed which is independent of the rupture growth length l in contrast during phase ii the rupture develops at accelerating speeds v which increase with an increase in l obeying a power law v vs l c n where vs is the shear wave velocity c is the characteristic length representing the geometric irregularity of the fault surfaces and and n are constants 8 87 10 29 and n 7 31 scale dependency of scale dependent physical quantities including the nucleation zone size and its duration is commonly ascribed to scale dependency of the slip dependent constitutive law parameter dc which is in turn governed by c it has been discussed that a unified comprehension can be provided for shear rupture of any size scale if the constitutive law for shear rupture is formulated as a slip dependent law

slow and fast ruptures on a laboratory fault controlled by loading characteristics recent geodetic observations indicate that a single fault may slip slowly and silently or fast and seismically in different circumstances we report laboratory experiments that demonstrate how the mode of faulting can alternate between fast stick slip events 70 to 500 mm s sliding rates and slow silent events 0 001 to 30 mm s as a result of loading conditions rather than friction properties or stiffness of the loading machine the 760 mm long granite sample is close to the critical nucleation length scale for unstable slip so small variations in nucleation properties result in measurable differences in slip events slow events occur when instability cannot fully nucleate before reaching the sample ends dynamic events occur after long healing times or abrupt increases in loading rate which suggests that these factors shrink the spatial and temporal extents of the nucleation zone arrays of slip strain and ground motion sensors installed on the sample allow us to quantify seismic coupling and study details of premonitory slip and afterslip we find that seismic coupling decreases when slip rates fall below about 70 mm s the slow slip events we observe are primarily aseismic less than 1 of the seismic coupling of faster events and produce swarms of very small m 6 to m 8 events these mechanical and seismic interactions suggest that faults with transitional behavior where creep small earthquakes and tremor are often observed could become seismically coupled if loaded rapidly either by a slow slip front or dynamic rupture of an earthquake that nucleated elsewhere

the dynamics of the onset of frictional slip the way in which a frictional interface fails is critical to our fundamental understanding of failure processes in fields ranging from engineering to the study of earthquakes frictional motion is initiated by rupture fronts that propagate within the thin interface that separates two sheared bodies by measuring the shear and normal stresses along the interface together with the subsequent rapid real contact area dynamics we find that the ratio of shear stress to normal stress can locally far exceed the static friction coefficient without precipitating slip moreover different modes of rupture selected by the system correspond to distinct regimes of the local stress ratio these results indicate the key role of nonuniformity to frictional stability and dynamics with implications for the prediction selection and arrest of different modes of earthquakes

large nucleation before large earthquakes is sometimes skipped due to cascade up implications from a rate and state simulation of faults with hierarchical asperities 1 does a large earthquake have a large quasi static preparation a hierarchical asperity model in which a large tough patch patch l contains smaller fragile patches patch s enables a large earthquake to start with only small preparation because of cascade up rupture growth we realized such a model in a rate and state framework by heterogeneous distributions of the state evolution distance and its consequences are investigated by earthquake sequence simulations focus is put on elementary processes the interaction between two scales one patch l containing one patch s is simulated with their size ratio as a parameter if patch s is larger than the nucleation size of patch l the system falls into a limit cycle consisting of only one earthquake that starts with patch s nucleation and grows dynamically to rupture entire patch l if patch s is considerably smaller than the nucleation size of patch l small earthquakes never dynamically cascade up and the large earthquakes are initiated by large quasi static nucleation in between large earthquakes start in various ways by large nucleation dynamic cascade up or delayed cascade up in the final stage of quasi static nucleation the preseismic moment release rate increases roughly inversely proportional to the time to failure tf with its amplitude depending on the nucleation size for a patch s rupture to cascade up strength in the adjacent region must have been reduced manifested by a higher creep velocity before the patch s nucleation starts following 1 tf acceleration large nucleation sometimes has a precursory small earthquake characterized by larger afterslip than nonprecursory ones