mentawai fault zone off sumatra a new key to the geodynamics of western indonesia the geodynamic evolution of the western part of the sunda arc is controlled by the change from frontal subduction of the indo australian plate along java to oblique subduction along sumatra this obliquity gives rise to the sumatra fault zone that links the accretionary zone of the andaman sea to the sunda strait previous studies have shown a decrease of displacement rate of the movement along the fault zone from north to south consequently it has been proposed that the area between the subduction zone and the fault zone i e the sumatra sliver platelet was deformed an oceanographic cruise on the indonesian ship r v baruna jaya iii was designed to study this area seismic reflection data reveal the existence of a major strike slip fault parallel to the sumatra fault zone that we called the mentawai fault zone located in the fore arc area just east of the mentawai islands it is at least 600 km long thus the sumatra sliver plate appears to be composed of several strips that move toward the northwest to accommodate the oblique subduction

decorrelation in interferometric radar echoes a radar interferometric technique for topographic mapping of surfaces implemented utilizing a single synthetic aperture radar sar system in a nearly repeating orbit is discussed the authors characterize the various sources contributing to the echo correlation statistics and isolate the term which most closely describes surficial change they then examine the application of this approach to topographic mapping of vegetated surfaces which may be expected to possess varying backscatter over time it is found that there is decorrelation increasing with time but that digital terrain model generation remains feasible the authors present such a map of a forested area in oregon which also includes some nearly unvegetated lava flows such a technique could provide a global digital terrain map

implications of the 26 december 2004 sumatra andaman earthquake on tsunami forecast and assessment models for great subduction zone earthquakes results from different tsunami forecasting and hazard assessment mod els are compared with observed tsunami wave heights from the 26 december 2004 indian ocean tsunami forecast models are based on initial earthquake information and are used to estimate tsunami wave heights during propagation an empirical forecast relationship based only on seismic moment provides a close estimate to the observed mean regional and maximum local tsunami runup heights for the 2004 indian ocean tsunami but underestimates mean regional tsunami heights at azimuths in line with the tsunami beaming pattern e g sri lanka thailand standard forecast models developed from subfault discretization of earthquake rupture in which deep ocean sea level observations are used to constrain slip are also tested forecast models of this type use tsunami time series measurements at points in the deep ocean as a proxy for the 2004 indian ocean tsunami a transect of deep ocean tsunami amplitudes recorded by satellite altimetry is used to constrain slip along four subfaults of the m 9 sumatra andaman earthquake this proxy model performs well in comparison to observed tsunami wave heights travel times and inundation patterns at banda aceh hypothetical tsunami hazard assessments models based on end member estimates for average slip and rupture length m w 9 0 9 3 are compared with tsunami observations using average slip low end member and rupture length high end member m w 9 14 consistent with many seismic geodetic and tsunami inversions adequately estimates tsunami runup in most regions except the extreme runup in the western aceh province the high slip that occurred in the southern part of the rupture zone linked to runup in this location is a larger fluctuation than expected from standard stochastic slip models in addition excess moment release 9 deduced from geodetic studies in comparison to seismic moment estimates may gen erate additional tsunami energy if the exponential time constant of slip is less than approximately 1 hr overall there is significant variation in assessed runup heights caused by quantifiable uncertainty in both first order source parameters e g rupture length slip length scaling and spatiotemporal complexity of earthquake rupture

earthquake supercycles inferred from sea level changes recorded in the corals of west sumatra records of relative sea level change extracted from corals of the mentawai islands sumatra imply that this 700 kilometer long section of the sunda megathrust has generated broadly similar sequences of great earthquakes about every two centuries for at least the past 700 years the moment magnitude 8 4 earthquake of september 2007 represents the first in a series of large partial failures of the mentawai section that will probably be completed within the next several decades

effects of kinematic constraints on teleseismic finite source rupture inversions great peruvian earthquakes of 23 june 2001 and 15 august 2007 two great underthrusting earthquakes that occurred along the coast of peru in 2001 and 2007 involve spatiotemporal slip distributions that differ from the predominantly unilateral or bilateral rupture expansion of many great events commonly used finite source rupture model parameterizations with specified rupture velocity and or short duration of slip at each grid point applied to the seismic data for these two events lead to incorrect slip distributions or inaccurate estimation of rupture velocities as a result of intrinsic kinematic constraints imposed on the model slip distributions guided by large aperture array back projections of teleseismic broadband p wave signals that image slip locations without imposing a priori kinematic constraints on the rupture process we exploit the availability of large global broadband body and surface wave data sets to consider the effects of varying the kinematic constraints in teleseismic finite source waveform inversions by allowing longer than usual rupture durations at each point on the fault using a flexible subfault source time function parameterization we find that the anomalous attributes of the 2001 and 2007 peru earthquake ruptures are readily recognized and accounted for by compound rupture models the great 23 june 2001 m w 8 4 8 4 earthquake involved an initial modest size event that appears to have triggered a much larger secondary event about 120 km away that developed an overall slip distribution with significant slip located back along the megathrust in the vicinity of the initial rupture the great 15 august 2007 m w 8 0 8 0 earthquake was also a composite event with a modest size initial rupture followed by a 60 sec delayed larger rupture that initiated 50 60 km away and spread up dip and bilaterally when back projections indicate greater rupture complexity than captured in a simple slip pulse type rupture model one should allow for possible long subfault slip duration or composite triggered sequences and not overly constrain the earthquake slip distribution

the tsunami of 2007 september 12 bengkulu province sumatra indonesia post tsunami field survey and numerical modelling summary the m w 8 4 earthquake on 2007 september 12 offshore of the bengkulu province of sumatra indonesia generated a moderate tsunami with run up heights of up to 4 m as measured by indonesian and international researchers in the days following the earthquake the tsunami was observed along 250 km of coastline and caused damage at several locations the largest wave heights and most severe inundation were observed about 50 km to the northwest of bengkulu elsewhere the effects were less severe with the exception of substantial inundation at a site 150 km to the south in addition to presenting the field data we conduct a modelling study and compare the run up heights and water level predictions from four different seismic deformation models used to initialize a tsunami propagation and inundation model our comparative results suggest that for this event the estimates of fault parameters available immediately after determination of the earthquake size and location predicted the near field run up heights and distribution and far field wave equally well as estimates obtained using more detailed descriptions of the seafloor deformation as typically available hours or days after an event we conclude that while detailed slip patterns can be important to the near field run up distribution simple fault models can be used to rapidly assess the likely near and far field tsunami effects of a particular earthquake

inversion of tsunami waveforms for the estimation of heterogeneous fault motion of large submarine earthquakes the 1968 tokachi oki and 1983 japan sea earthquakes waveform inversion of tsunami data provides independent and reliable information on earthquake source processes the propagation effects of tsunamis are accounted for by computing green s functions using a finite difference method and using the known bathymetry which is much better constrained than the seismic velocity structure the observed tsunami waveforms are corrected for the response of the tide gage system by using the results of in situ measurements using green s functions computed from bathymetric data on a 2 5 km grid and the observed waveforms digitized at 1 min intervals the spatial resolution of the faulting heterogeneity by the present method is as small as 30 km this method is applied to the 1968 tokachi oki and 1983 japan sea earthquakes the slip is concentrated at the northwestern part of the fault plane for the 1968 event for the 1983 event the slip on the fault just north of the epicenter is twice that of the other parts these asperity distributions are similar to those inferred from seismic wave analysis however the recovered total seismic moments are smaller the present method is shown to be a useful technique for estimating heterogeneous fault motion especially for historical events with no or poor seismic records and tsunami earthquakes in which the excitations of seismic and tsunami waves are different

source process of the september 12 2007 mw 8 4 southern sumatra earthquake from tsunami tide gauge record inversion 1 we infer the slip distribution and average rupture velocity of the magnitude mw 8 4 september 12 2007 southern sumatra earthquake from available tide gauge records of the ensuing tsunami we select 9 waveforms recorded along the west coast of sumatra and in the indian ocean slip distribution and rupture velocity are determined simultaneously by means of a non linear inversion method we find high slip values 10 m into a patch 100 km long and 50 km large between 20 and 30 km of depth about 100 km north west from the epicenter we conclude this earthquake did not rupture the whole area of the 1833 event indicating some slip has still to occur our estimate of rupture velocity is of 2 1 0 4 km sec the relatively large depth of the main slip patch is the likely explanation for the low damaging observed tsunami

distribution of slip at the northern sumatran fault system we model spatial variations in horizontal displacements of 117 geodetic sites measured during annual surveys in 1989 1996 with the global positioning system gps as elastic strain across a locked strike slip fault to infer the contemporary slip rate locking depth and location of the sumatran fault sf in northern sumatra 1 s 3 n gps derived slip rate estimates increase slightly northward from 23 plus or minus 3 mm yr at 0 8 deg s to 26 plus or minus 2mm yr at 2 7 n they agree with geologic estimates north of the equator but at 0 5 s they are about 10 mm yr higher strain appears to be distributed asymmetrically about the fault south of 2 n about 5 mm yr of shear is required within the offshore forearc west of the fault to achieve a closer agreement of fault locations inferred from gps velocities with geologically identified traces of the sf locking depth estimates are on the order of 10 20 km the western branch of the major fault bifurcation near 1 n slips at a rate five times higher than the eastern branch the two main strands of the fault at the northwestern tip of sumatra 5 5 n appear to be nearly free of horizontal strain significant slip must occur away from the two strands probably further east at two other geologically active branches the banda aceh embayment is extruded to the northwest at a rate of 5 plus or minus 2 mm yr within the estimated velocity uncertainties of several mm yr fault normal deformation along the sf is insignificant almost strain free the northern part of the back arc basin is part of a rigid sunda shelf while the northern forearc is subjected to 8 plus or minus 5 x 10 exp 8 yr of extension nearly parallel to the arc

kuril islands tsunami of november 2006 1 impact at crescent city by distant scattering 1 a numerical model for the global tsunami computation constructed by kowalik et al 2005 2007a is applied to the tsunami of november 15 2006 in the northern pacific with spatial resolution of one minute numerical results are compared to sea level data collected by pacific dart buoys the tide gauge at crescent city cc recorded an initial tsunami wave of about 20 cm amplitude and a second larger energy packet arriving 2 hours later the first energy input into the cc harbor was the primary direct wave traveling over the deep waters of the north pacific interactions with submarine ridges and numerous seamounts located in the tsunami path were a larger source of tsunami energy than the direct wave travel time for these amplified energy fluxes is longer than for the direct wave prime sources for the larger fluxes at cc are interactions with koko guyot and hess rise tsunami waves travel next over the mendocino escarpment where the tsunami energy flux is concentrated owing to refraction and directed toward cc local tsunami amplification over the shelf break and shelf are important as well in many locations along the north pacific coast the first arriving signal or forerunner has lower amplitude than the main signal which often is delayed understanding this temporal distribution is important for an application to tsunami warning and prediction as a tsunami hazard mitigation tool we propose that along with the sea level records which are often quite noisy an energy flux for prediction of the delayed tsunami signals be used

coseismic slip and afterslip of the great mw 9 15 sumatra andaman earthquake of 2004 we determine coseismic and the first month postseismic deformation associated with the sumatra andaman earthquake of 26 december 2004 from near field global positioning system gps surveys in northwestern sumatra and along the nicobar andaman islands continuous and campaign gps measurements from thailand and malaysia and in situ and remotely sensed observations of the vertical motion of coral reefs the coseismic model shows that the sunda subduction mega thrust ruptured over a distance of about 1500 km and a width of less than 150 km releasing a total moment of 6 7 7 0 10 22 n m equivalent to a magnitude m w 9 15 the latitudinal distribution of released moment in our model has three distinct peaks at about 4 n 7 n and 9 n which compares well to the latitudinal variations seen in the seismic inversion and of the analysis of radiated t waves our coseismic model is also consistent with interpretation of normal modes and with the amplitude of very long period surface waves the tsunami predicted from this model fits rela tively well the altimetric measurements made by the jason and topex satellites neither slow nor delayed slip is needed to explain the normal modes and the tsunami wave the near field geodetic data that encompass both coseismic deformation and up to 40 days of postseismic deformation require that slip must have continued on the plate interface after the 500 sec long seismic rupture the postseismic geodetic moment of about 2 4 10 22 nm mw 8 8 is equal to about 30 5 of the coseismic moment release evolution of postseismic deformation is consistent with rate strengthening frictional afterslip online material summary of geodetic data used in this study

radar interferometry and its application to changes in the earth s surface geophysical applications of radar interferometry to measure changes in the earth s surface have exploded in the early 1990s this new geodetic technique calculates the interference pattern caused by the difference in phase between two images acquired by a spaceborne synthetic aperture radar at two distinct times the resulting interferogram is a contour map of the change in distance between the ground and the radar instrument these maps provide an unsurpassed spatial sampling density 100 pixels km 2 a competitive precision 1 cm and a useful observation cadence 1 pass month 1 they record movements in the crust perturbations in the atmosphere dielectric modifications in the soil and relief in the topography they are also sensitive to technical effects such as relative variations in the radar s trajectory or variations in its frequency standard we describe how all these phenomena contribute to an interferogram then a practical summary explains the techniques for calculating and manipulating interferograms from various radar instruments including the four satellites currently in orbit ers 1 ers 2 jers 1 and radarsat the next chapter suggests some guidelines for interpreting an interferogram as a geophysical measurement respecting the limits of the technique assessing its uncertainty recognizing artifacts and discriminating different types of signal we then review the geophysical applications published to date most of which study deformation related to earthquakes volcanoes and glaciers using ers 1 data we also show examples of monitoring natural hazards and environmental alterations related to landslides subsidence and agriculture in addition we consider subtler geophysical signals such as postseismic relaxation tidal loading of coastal areas and interseismic strain accumulation we conclude with our perspectives on the future of radar interferometry the objective of the review is for the reader to develop the physical understanding necessary to calculate an interferogram and the geophysical intuition necessary to interpret it

rupture process of the 2004 sumatra andaman earthquake the 26 december 2004 sumatra andaman earthquake initiated slowly with small slip and a slow rupture speed for the first 40 to 60 seconds then the rupture expanded at a speed of about 2 5 kilometers per second toward the north northwest extending 1200 to 1300 kilometers along the andaman trough peak displacements reached 15 meters along a 600 kilometer segment of the plate boundary offshore of northwestern sumatra and the southern nicobar islands slip was less in the northern 400 to 500 kilometers of the aftershock zone and at least some slip in that region may have occurred on a time scale beyond the seismic band

the great sumatra andaman earthquake of 26 december 2004 the two largest earthquakes of the past 40 years ruptured a 1600 kilometer long portion of the fault boundary between the indo australian and southeastern eurasian plates on 26 december 2004 seismic moment magnitude mw 9 1 to 9 3 and 28 march 2005 mw 8 6 the first event generated a tsunami that caused more than 283 000 deaths fault slip of up to 15 meters occurred near banda aceh sumatra but to the north along the nicobar and andaman islands rapid slip was much smaller tsunami and geodetic observations indicate that additional slow slip occurred in the north over a time scale of 50 minutes or longer

heterogeneous coupling of the sumatran megathrust constrained by geodetic and paleogeodetic measurements geodetic and paleogeodetic measurements of interseismic strain above the sumatran portion of the sunda subduction zone reveal a heterogeneous pattern of coupling annual banding in corals provides vertical rates of deformation spanning the last half of the 20th century and repeated gps surveys between 1991 and 2001 and continuous measurements at gps stations operated since 2002 provide horizontal velocities near the equator the megathrust is locked over a narrow width of only a few tens of kilometers in contrast the locked fault zone is up to about 175 km wide in areas where great interplate earthquakes have occurred in the past formal inversion of the data reveals that these strongly coupled patches are roughly coincident with asperities that ruptured during these events the correlation is most spectacular for rupture of the m w 8 7 nias simeulue earthquake of 2005 which released half of the moment deficit that had accumulated since its previous rupture in 1861 suggesting that this earthquake was overdue beneath the mentawai islands strong coupling is observed within the overlapping rupture areas of the great earthquakes of 1797 and 1833 the accumulated slip deficit since these events is slowly reaching the amount of slip that occurred during the 1833 earthquake but already exceeds the slip that occurred during the 1797 earthquake thus rerupture of part of the mentawai patch in september 2007 was not a surprise in contrast coupling is low below the batu islands near the equator and around enggano island at about 5 s where only moderate earthquakes m w 8 0 have occurred in the past two centuries the correlation of large seismic asperities with patches that are locked during the interseismic period suggests that they are persistent features this interpretation is reinforced by the fact that the large locked patches and great ruptures occur beneath persistent geomorphologic features the largest outer arc islands depth and convergence rate dependent temperature might influence the pattern of coupling through its effect on the rheology of the plate interface but other influences are required to account for the observed along strike heterogeneity of coupling in particular subduction of the investigator fracture zone could be the cause for the low coupling near the equator

tsunami sources of the november 2006 and january 2007 great kuril earthquakes abstract along the northern kuril trench a great interplate earthquake occurred in november 2006 followed by an outer rise normal fault event in january 2007 surface wave magnitudes m s 7 8 and 8 2 indicate that the 2007 event was larger while the global cmt solutions m w 8 3 and 8 1 indicate that the 2006 event was larger tsunamis from both events were recorded at tide gauge stations in japan russia and the united states as well as deep ocean assessment and reporting of tsunamis dart systems and cabled tsunami sensors installed on deep ocean bottom inversion of 52 tsunami waveforms indicates that the 2006 tsunami source was about 200 km long extending from the epicenter to northeast the largest slips of 4 7 m are estimated on the northeastern subfaults for the 2007 event inversion of 32 waveforms shows that the tsunami source was about 120 km long with the largest slip of 3 5 m on the northeastern subfault the agreements between the observed and synthetic tsunami waveforms are generally good for both events not only within the inversion time windows but also for the later phases which were not used in inversions the slip distributions yield the seismic moment of m w 8 1 for the 2006 event and m w 7 9 for the 2007 event the seismic moment of the 2006 event was larger than the 2007 event from comparison of tsunami data regardless of fault model

synthetic aperture radar interferometry to measure earth s surface topography and its deformation synthetic aperture radar interferometry insar from earth orbiting spacecraft provides a new tool to map global topography and deformation of the earth s surface radar images taken from slightly different viewing directions allow the construction of digital elevation models of meter scale accuracy these data sets aid in the analysis and interpretation of tectonic and volcanic landscapes if the earth s surface deformed between two radar image acquisitions a map of the surface dis placement with tens of meters resolution and subcentimeter accuracy can be con structed this review gives a basic overview of insar for earth scientists and presents a selection of geologic applications that demonstrate the unique capabilities of insar for mapping the topography and deformation of the earth

neotectonics of the sumatran fault indonesia the 1900 km long trench parallel sumatran fault accommodates a significant amount of the right lateral component of oblique convergence between the eurasian and indian australian plates from 10 n to 7 s our detailed map of the fault compiled from topographic maps and stereographic aerial photographs shows that unlike many other great strike slip faults the sumatran fault is highly segmented cross strike width of step overs between the 19 major subaerial segments is commonly many kilometers the influence of these step overs on historical seismic source dimensions suggests that the dimensions of future events will also be influenced by fault geometry geomorphic offsets along the fault range as high as 20 km and may represent the total offset across the fault if this is so other structures must have accommodated much of the dextral component of oblique convergence during the past few million years our analysis of stretching of the forearc region near the southern tip of sumatra constrains the combined dextral slip on the sumatran and mentawai faults to be no more than 100 km in the past few million years the shape and location of the sumatran fault and the active volcanic arc are highly correlated with the shape and character of the underlying subducting oceanic lithosphere nonetheless active volcanic centers of the sumatran volcanic arc have not influenced noticeably the geometry of the active sumatran fault on the basis of its geologic history and pattern of deformation we divide the sumatran plate margin into northern central and southern domains we support previous proposals that the geometry and character of the subducting investigator fracture zone are affecting the shape and evolution of the sumatran fault system within the central domain the southern domain is the most regular the sumatran fault there comprises six right stepping segments this pattern indicates that the overall trend of the fault deviates 4 clockwise from the slip vector between the two blocks it separates the regularity of this section and its association with the portion of the subduction zone that generated the giant m w 9 earthquake of 1833 suggest that a geometrically simple subducting slab results in both simple strike slip faulting and unusually large subduction earthquakes

tsunami source of the 2004 sumatra andaman earthquake inferred from tide gauge and satellite data tsunami source of the 2004 sumatra andaman earthquake was esti mated from a joint inversion of tsunami waveforms recorded on tide gauges and sea surface heights captured by satellite altimetry measurements the earthquake the largest in the past 40 years caused devastating tsunami damage to countries around the indian ocean but the tsunami source in particular its northern end was not well resolved although aftershocks and crustal deformation extended from off north western sumatra island through the nicobar islands to the andaman islands some seismic wave analyses indicated a shorter source length several hundred kilometers we used tsunami waveforms recorded at 12 tide gauge stations around the source and the sea surface heights measured by three satellites jason1 topex poseidon and envisat we numerically computed tsunami propagation using realistic bathym etry more than 66 000 depth points were digitized from nautical charts and combined with the etopo2 data inversion of satellite data indicates that the tsunami source extended to the andaman islands with a total length of 1 400 km but such a model produces much larger tsunami waveforms than observed at indian tide gauge stations inversion of tide gauge records and the joint inversion indicate that the tsunami source was about 900 km long the largest slip about 13 to 25 m was located off sumatra island and the second largest slip up to 7 m near the nicobar islands the slip distribution is similar for different rupture velocities and rise times with a slow velocity of 1 km sec and a rise time of 3 min yielding the largest variance reduction

source parameters of the great sumatran megathrust earthquakes of 1797 and 1833 inferred from coral microatolls large uplifts and tilts occurred on the sumatran outer arc islands between 0 5 and 3 3 s during great historical earthquakes in 1797 and 1833 as judged from relative sea level changes recorded by annually banded coral heads coral data for these two earthquakes are most complete along a 160 km length of the mentawai islands between 3 2 and 2 s uplift there was as great as 0 8 m in 1797 and 2 8 m in 1833 uplift in 1797 extended 370 km between 3 2 and 0 5 s the pattern and magnitude of uplift imply megathrust ruptures corresponding to moment magnitudes m w in the range 8 5 to 8 7 the region of uplift in 1833 ranges from 2 to at least 3 2 s and judging from historical reports of shaking and tsunamis perhaps as far as 5 s the patterns and magnitude of uplift and tilt in 1833 are similar to those experienced farther north between 0 5 and 3 n during the giant nias simeulue megathrust earthquake of 2005 the outer arc islands rose as much as 3 m and tilted toward the mainland elastic dislocation forward modeling of the coral data yields megathrust ruptures with moment magnitudes ranging from 8 6 to 8 9 sparse accounts at padang along the mainland west coast at latitude 1 s imply tsunami runups of at least 5 m in 1797 and 3 4 m in 1833 tsunamis simulated from the pattern of coral uplift are roughly consistent with these reports the tsunami modeling further indicates that the indian ocean tsunamis of both 1797 and 1833 unlike that of 2004 were directed mainly south of the indian subcontinent between about 0 7 and 2 1 s the lack of vintage 1797 and 1833 coral heads in the intertidal zone demonstrates that interseismic submergence has now nearly equals coseismic emergence that accompanied those earthquakes the interseismic strains accumulated along this reach of the megathrust have thus approached or exceeded the levels relieved in 1797 and 1833

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

slip distribution of the 2007 bengkulu earthquake inferred from tsunami waveforms and insar data 1 joint inversion study using tsunami waveforms and insar data provides a new way to understand the magnitude and spatial extent of subduction zone earthquakes a great earthquake mw 8 5 occurred on 12 september 2007 off the west coast of bengkulu indonesia the tsunami generated by the event was recorded by tide gauge stations around the indian ocean and by two bottom pressure sensors in the deep sea the ground surface displacements produced by the earthquake on pagai islands and on sumatra island were observed by interferometric synthetic aperture radar insar we estimated the slip distribution of the earthquake by joint inversion incorporating a spatial smoothness constraint using tsunami waveforms and insar data the total released seismic moment calculated from the slip distribution is 6 7 1021 n m mw 8 5 consistent with the seismic moment of the global cmt solution 6 71 1021 n m the maximum observed tsunami heights along the coast of bengkulu agree with those computed from the slip distribution the slip amount of the 2007 earthquake is smaller than the amount of potential slip that has been accumulated since after the 1797 and 1833 events the premise is that averaged over long periods of time the entire fault must slip equal amounts therefore the 2007 great earthquake could be followed by several great earthquakes that will rupture the plate interface until the potential slip that has been accumulated is completely released

partial rupture of a locked patch of the sumatra megathrust during the 2007 earthquake sequence the great sumatra andaman earthquake and tsunami of 2004 was a dramatic reminder of the importance of understanding the seismic and tsunami hazards of subduction zones 1 2 3 4 in march 2005 the sunda megathrust ruptured again producing an event 5 of moment magnitude mw 8 6 south of the 2004 rupture area which was the site of a similar event in 1861 ref 6 concern was then focused on the mentawai area where large earthquakes had occurred in 1797 mw 8 8 and 1833 mw 9 0 6 7 two earthquakes one of mw 8 4 and twelve hours later one of mw 7 9 indeed occurred there on 12 september 2007 here we show that these earthquakes ruptured only a fraction of the area ruptured in 1833 and consist of distinct asperities within a patch of the megathrust that had remained locked in the interseismic period this indicates that the same portion of a megathrust can rupture in different patterns depending on whether asperities break as isolated seismic events or cooperate to produce a larger rupture this variability probably arises from the influence of non permanent barriers zones with locally lower pre stress due to the past earthquakes the stress state of the portion of the sunda megathrust that had ruptured in 1833 and 1797 was probably not adequate for the development of a single large rupture in 2007 the moment released in 2007 amounts to only a fraction both of that released in 1833 and of the deficit of moment that had accumulated as a result of interseismic strain since 1833 the potential for a large megathrust event in the mentawai area thus remains large

deformation and slip along the sunda megathrust in the great 2005 nias simeulue earthquake seismic rupture produced spectacular tectonic deformation above a 400 kilometer strip of the sunda megathrust offshore northern sumatra in march 2005 measurements from coral microatolls and global positioning system stations reveal trench parallel belts of uplift up to 3 meters high on the outer arc islands above the rupture and a 1 meter deep subsidence trough farther from the trench surface deformation reflects more than 11 meters of fault slip under the islands and a pronounced lessening of slip trenchward a saddle in megathrust slip separates the northwestern edge of the 2005 rupture from the great 2004 sumatra andaman rupture the southeastern edge abuts a predominantly aseismic section of the megathrust near the equator

the global reach of the 26 december 2004 sumatra tsunami numerical model simulations combined with tide gauge and satellite altimetry data reveal that wave amplitudes directionality and global propagation patterns of the 26 december 2004 sumatra tsunami were primarily determined by the orientation and intensity of the offshore seismic line source and subsequently by the trapping effect of mid ocean ridge topographic waveguides

coseismic gravity changes of the 2011 tohoku oki earthquake from satellite gravimetry 1 the massive tohoku oki earthquake of a moment magnitude mw of 9 0 occurred on 11 march 2011 off the pacific coast of the northeastern japan the mass redistribution in and around the focal region associated with this earthquake was studied using the gravity changes detected by gravity recovery and climate experiment grace satellite after the 2004 sumatra andaman and the 2010 central chile maule earthquakes the present study presents the third case of clear detection of coseismic gravity changes by grace the observed gravity changes were dominated by decrease over the back arc region of 7 gal or less this reflects to a large extent coseismic crustal dilatation of the landward plate they agree well with the changes calculated with the green s function for the realistic earth using fault parameters inferred from coseismic crustal movements the spatial patterns of the gravity changes of these earthquakes are very similar because they are all shallow angle reverse faulting at convergent plate boundaries we found linear relationship between gravity decreases and seismic moments

contiguous rupture areas of two nankai trough earthquakes revealed by high resolution tsunami waveform inversion 1 we have developed a new method for inverting tsunami waveforms that reveals considerable detail in megathrust slip during subduction zone earthquakes previous methods have ensured compliance with the shallow water theory used to compute tsunami waveforms by using large subfaults that generate only long wavelength seafloor deformation we show that a better approach is to use small subfaults coupled with a smoothing criterion that ensures compliance with the shallow water approximation in an application of the method to historical earthquakes in the nankai trough southwestern japan we find that the areas with slip 1 m for the earthquakes of 1944 and 1946 which ruptured adjacent segments of the subduction zone are separated by a sharp non overlapping boundary this establishes that interseismic accumulation of strain energy extends very close to the boundary between rupture zones and strongly suggests that this boundary is associated with a physical barrier to rupture

asperity distributions and large earthquake occurrence in subduction zones abstract ruff l j 1992 asperity distributions and large earthquake occurrence in subduction zones in t mikumo k aki m ohnaka l j ruff and p k p spudich editors earthquake source physics and earthquake precursors tectonophysics 211 61 83 plate tectonics and the seismic gap hypothesis provide the framework for long term earthquake forecasting of plate boundary earthquakes unfortunately detailed examination reveals that earthquake recurrence times and rupture length vary between successive earthquake cycles in the same subduction zone furthermore larger coseismic slip is commonly associated with larger rupture length hence large earthquake occurrence in subduction zones is characterized by variability in 1 recurrence times 2 rupture length and 3 coseismic slip these facts plus many other observations indicate that there are significant spatial variations in the strength of the plate interface one simple description of these variations and their role in the earthquake cycle is the asperity model where the large strong regions of the plate interface are called asperities and the large earthquakes occur when the large asperities break the asperity model of earthquake occurrence is able to qualitatively explain several features of large plate boundary earthquakes to go beyond general qualitative notions i pose the following scientific test are the observed asperity distributions and a simple model of their interaction self consistent with the above three observed features of large earthquake occurrence the distribution of the major asperities along plate boundary segments has now been determined for several subduction zones rupture process studies of adjacent large and great earthquakes have provided reliable estimates of the along strike asperity lengths and separations for several adjacent asperities in the kurile islands colombia and peru subduction zones the simplest mechanical model for asperity interaction is to idealize two adjacent asperities as frictional sliders that are connected by main springs to the upper plate by a coupling spring to each other and maintain frictional contact with a conveyer belt the lower plate that moves with a constant velocity an earthquake occurs when the net force on the asperity frictional slider reaches some specified level the failure force and spring constants are determined by the observed asperity distribution and simple models of elastic interaction two different macroscopic failure criteria are used this simple mechanical model displays a remarkable range of behavior from simple to complex when the two asperities are identical in all their properties sequences of identical earthquakes are produced for the more realistic case of non identical asperities earthquake sequences show great variety using system variables from the observed asperity distributions the earthquake sequences typically display 1 variable recurrence times 2 variable rupture length i e a combination of single asperity and double asperity failures and for one of the failure criteria 3 larger coseismic slip for double asperity failures statistical summaries of thousands of simulated earthquake sequences for asperity pairs in the kuriles colombia and peru subduction zones are broadly consistent with the observed features of large earthquake occurrence in these subduction zones the main conclusion is that the asperity model provides a self consistent explanation for fault zone heterogeneity the rupture process and recurrence times and rupture mode of large earthquake sequences via a simple model for adjacent asperity interaction in addition a conclusion independent of any particular model for fault zone heterogeneity is that simple deterministic models of fault zone interaction can explain complex patterns of large earthquake occurrence in subduction zones

seismic history and seismotectonics of the sunda arc historic records of the last 300 years reveal two great interplate earthquakes 1833 mw 8 1861 mw 8 8 which ruptured major segments of the sumatra fore arc in western indonesia a significant percentage of interplate slip along this portion of the plate boundary appears to occur seismically the ends of these rupture zones are coincident with clusters of large ms 7 and moderate 6 ms 7 shocks and with heterogeneities in the plate interface as inferred from geologic and geophysical data the northern extent of the 1833 rupture zone is coincident with a group of moderate and large earthquakes a structural arch in the fore arc and a fracture zone on the subducted plate the southern portion is coincident with a very prominent cluster and the along strike extension of a sediment filled trough on the subducted plate prior to this study sumatra was characterized as relatively aseismic as inferred from the lack of great earthquakes in the instrumental record of this century java and the lesser sunda islands had major ms 6 earthquakes in the historic record but none are of the same magnitude as the great events near sumatra there is no clear evidence that shocks in the java fore arc were interplate events yet geophysical data indicate recent uplift of the trench slope break and a change in strike of the trench axis where the northernmost flank of the roo rise a prominent bathymetrie high appears to be interacting with the fore arc this feature was formed at a ridge crest and is more buoyant and thus more difficult to subduct than the surrounding seafloor there is no increase in seismicity associated with this tectonism however that may indicate a lack of contact between the crystalline portions of the plates at this time we infer that a majority of slip on the plate interface near java occurs aseismically the variation in occurrence of great interplate earthquakes along the sunda arc may be interpreted in terms of a model that attributes variation in interplate coupling to the age of the subducted lithosphere great interplate earthquakes occur near sumatra where the age of the youngest crust is 46 ma the characteristics of the sunda arc and analogies with pacific arcs imply that the entire length of sumatra has the potential to produce great thrust earthquakes its seismic potential should be considered high as there has been no recurrence of the great events of 1833 and 1861 in contrast no such events have been reported off java and the lesser sunda islands where the oldest crust is 152 ma the plate interface near java and the lesser sunda islands should be considered to have low seismic potential

optimum sea surface displacement and fault slip distribution of the 2017 tehuantepec earthquake m w 8 2 in mexico estimated from tsunami waveforms

geochemistry and hydrology of a calcareous fen within the savage fen wetlands complex minnesota usa abstract savage fen is a wetlands complex at the base of north facing bluffs in the minnesota river valley the complex includes 27 8 hectares of calcareous fen that host rare calciphile plants whose populations are declining in minnesota water and sediment compositions in the calcareous fen were studied to gain a better understanding of the hydrologie system that sustains the rare vegetation groundwater in the fen is a calcium magnesium bicarbonate type with circumneutral ph values the groundwater composition is the resuit of interactions among water dissolved and gaseous carbon species carbonates and ion exchangers shallow groundwater is distinguished from deep groundwater by smaller concentrations of chloride sulfate magnesium and sodium and larger concentrations of calcium bicarbonate hydrogen sulfide and ammonium magnesian calcite is the prevalent carbonate in unconsolidated sedimentary fill beneath the fen and is an important source and sink for dissolved calcium magnesium and inorganic carbon calcite concentrations just below the water table are small because aerobic and anaerobic oxidation of organic matter increase the partial pressure of carbon dioxide pco2 decrease ph and cause calcite to dissolve thick calcite accumulations just above the water table in the root zone of calciphile plants resuit from water table fluctuations and attendant changes in pco2 groundwater beneath savage fen recharges in lakes and ponds south of the fen and upwells to the surface within the fen water at the water table is a mixture of upwelling groundwater and water near the surface that flows downslope from higher elevations in the fen changes in oxygen and hydrogen isotope compositions of shallow groundwater indicate that the proportion of upwelling groundwater in shallow groundwater decreases downgradient in the calcareous fen encroachment of reed grasses into the calcareous fen may reflect human caused disturbances in the recharge area

the great kurile earthquake of october 4 1994 tore the slab the fault geometry depth and slip distribution of the kurile earthquake of oct 4 1994 are estimated using seismic waveforms aftershock distribution geodetic measurements and tsunami waveforms previous large earthquakes in the kurile arc had typical underthrusting focal mechanisms due to the local subduction seismic wave inversions of the 1994 event indicate a thrust type mechanism with a large strike slip component this does not represent an underthrust event at the subduction interface the focal depth is estimated at 50 kin just beneath the previously estimated subduction interface we calculate crustal deformation and tsunami waveforms from the two possible fault planes of the focal mechanism solution and find that both fault planes explain the observations we choose the north south striking plane because of the occurrence of a large m w 7 8 intermediate depth earthquake on dec 6 1978 northwest of the 1994 rupture area the lineation of aftershock distributions of these two events suggests that both earthquakes were due to a tearing of the slab perpendicular to the trench the joint inversion of geodetic and tsunami data shows that the largest slip about 17 m occurred at the upper end of the slab

geodetic data inversion using a bayesian information criterion for spatial distribution of fault slip summary we developed a new inversion method to reconstruct static images of seismic sources from geodetic data using akaike s bayesian information criterion abic coseismic surface displacements are generally related with a slip distribution on a fault surface by linear integral equations parametric expansion of the fault slip distribution by a finite number of known basis functions yields a set of observation equations expressed in a simple vector form incorporating prior constraints on the smoothness of slip distribution with the observation equations we construct a bayesian model with unknown hyperparameters the optimal values of the hyperparameters which control the structure of the bayesian model are objectively determined from observed data by using abic once the values of hyperparameters are determined we can use the maximum likelihood method to find the optimal distribution of fault slip we examined the validity of this method through a numerical experiment using theoretical data with random noise we analysed geodetic data associated with the 1946 nankaido earthquake ms 8 2 by using this method the result shows that the fault slip distribution of this earthquake has two main peaks of 4 and 6 m located off kii peninsula and muroto promontory these two high slip areas are clearly separated by a low slip zone extending along kii strait such a slip distribution corresponds with the fact that the rupture process of this earthquake in the western part is notably different from that in the eastern part