Why does near ridge extensional seismicity occur primarily in the Indian Ocean?

The possibility that thermoelastic stresses due to plate cooling contribute significantly to the stress field and seismicity in young oceanic lithosphere has been a subject of considerable recent interest. This effect is suggested by three key observations: a decrease in seismicity with lithospheric age, the fact that focal mechanisms show extension perpendicular to the spreading direction, and a depth stratification of mechanism types. A difficulty with this idea is that although thermoelastic stresses should be comparable in different regions, the intraplate seismicity seems to occur in local concentrations. In particular, the ridge-parallel extensional seismicity occurs preferentially in the Central Indian Ocean region.We explore the possibility that much of the data favoring thermoelastic stresses can be interpreted in terms of stresses resulting from individual plate geometry and local boundary effects. In particular, the dramatic concentration of extensional seismicity in the Central Indian Ocean region is consistent with finite element results for the intraplate stress incorporating the effects of the Himalayan collision and the various subduction zones. The ridge parallel extensional stresses show a decrease with age similar to that of the seismicity. As earthquakes in this area provide a major portion of the data for both ridge-parallel extension and depth stratification, these effects may be due more to the regional stress. We thus propose that thermoelastic stresses provide a low level “background” in all plates, but that the dominant effect is that of individual plate geometry and local boundary effects.     

冲绳板块应力场数值模拟及其动力学特征div>

冲绳板块位于菲律宾海板块向欧亚板块俯冲形成的西太平洋边缘活动带上,构造应力场图像及其动力学机制表现得相当复杂．采用伪三维有限元方法,以WSM2008 观测应力场数据的应力取向和应力型两方面指标作为主要约束,对冲绳板块构造应力场进行了数值模拟．通过对计算结果的分析,对模型涉及的各种作用力作出了估计．在此基础上,对冲绳板块岩石层的状态,以及该地区的板块动力学特征进行了探讨,并得到了以下一些初步认识:① 软流层静压推力控制着该地区构造应力场的基本形态;② 冲绳海槽的演化过程,例如该地区的岩石层减薄与其下地幔流的上升等,也在很大程度上影响了该地区的板内应力场空间分布特征;③ 琉球俯冲带边界力的作用是分段的,不同区段作用力对板内应力场的影响有所不同．      

Apparent stress and stress drop for intraplate earthquakes and tectonic stress in the plates

The magnitude of shear stress in the lithosphere is bounded from below by the apparent stress and stress drop during intraplate earthquakes. Apparent stresses and stress drops for a number of mid-plate earthquakes are calculated from the earthquake magnitude, SH wave amplitude spectra, and estimates of the length of the fault zone. Apparent stresses vary between 0.1 and 2 bars, ifm _b is used as a measure of seismic energy, and stress drops lie between 2 and 70 bars. There is no systematic difference in either apparent stress or stress drop between these intraplate events and typical plate boundary earthquakes. These bounds on intraplate shear stresses are consistent with the inference from current models of plate tectonic driving forces that regional stress differences in the plates are typically on the order of 100 bars. The highest stress drops measured for midplate earthquakes under this model represent nearly total release of local tectonic stress.     

Experimental investigation on local shear stress and turbulence intensities over a rough non-uniform bed with and without sediment using 2D particle image velocimetry

The current work focuses on locally resolving velocities,turbulence,and shear stresses over a rough bed with locally non-uniform character.A nonporous subsurface layer and fixed interfacial sublayer of gravel and sand were water-worked to a nature-like bed form and additionally sealed in a hydraulic flume.Two-dimensional Particle Image Velocimetry(2 D-PIV) was applied in the vertical plane of the experimental flume axis.Runs with clear water and weak sediment transport were done under slightly supercritical flow to ensure sediment transport conditions without formation of considerable sediment deposits or dunes.The study design included analyzing the double-averaged flow parameters of the entire measurement domain and investigating the flow development at 14 consecutive vertical subsections.Local geometrical variabilities as well the presence of sediment were mainly reflected in the vertical velocity component.Whereas the vertical velocity decreased over the entire depth in presence of sediment transport,the streamwise velocity profile was reduced only within the interfacial sublayer.In the region with decelerating flow conditions,however,the streamwise velocity profile systematically increased along the entire depth extent.The increase in the main velocity(reduction of flow resistance)correlated with a decrease of the turbulent shear and main normal stresses.Therefore,effects of rough bed smoothening and drag force reduction were experimentally documented within the interfacial sublayer due to mobile sediment.Moreover,the current study leads to the conclusion that in nonuniform flows the maximum Reynolds stress values are a better predictor for the bed shear stress than the linearly extrapolated Reynolds stress profile.This is an important finding because,in natural flows,uniform conditions are rare.     

The mechanical state of the lithosphere in the Altiplano-Puna segment of the Andes

Information about topography, the shape of the geoid, seismicity, Neogene deformation and volcanism in the region of Altiplano-Puna of western South America is used to analyse the state of stress across the convergent plate margin in terms of the effects of topography and simple models of its compensation. An average elevation near 4 km is consistent with compensation by a yet unresolved combination of crustal root and hot uppermost mantle producing a geoid high of 22–27 meters, average horizontal compressive stress (in excess of a reference sea level lithostatic value) of 390 bars in a 150 km thick lithosphere, and an average shear stress of 170 bars along a 30° dipping interplate boundary. The basis for these estimates is evidence for a neutral to extensional stress regime within the high plateau contrasted with a compressional regime on the eastern slopes and along the interplate boundary itself. Comparison with other plateaus in a convergent plate tectonic setting suggests an evolutionary sequence from compressional to extensional tectonics as elevation of the plateau increases.     

Topographic stress and rock fracture: a two‐dimensional numerical model for arbitrary topography and preliminary comparison with borehole observations

Theoretical calculations indicate that elastic stresses induced by surface topography may be large enough in some landscapes to fracture rocks, which in turn could influence slope stability, erosion rates, and bedrock hydrologic properties. These calculations typically have involved idealized topographic profiles, with few direct comparisons of predicted topographic stresses and observed fractures at specific field sites. We use a numerical model to calculate the stresses induced by measured topographic profiles and compare the calculated stress field with fractures observed in shallow boreholes. The model uses a boundary element method to calculate the stress distribution beneath an arbitrary topographic profile in the presence of ambient tectonic stress. When applied to a topographic profile across the Susquehanna Shale Hills Critical Zone Observatory in central Pennsylvania, the model predicts where shear fractures would occur based on a Mohr–Coulomb criterion, with considerable differences in profiles of stresses with depth beneath ridgetops and valley floors. We calculate the minimum cohesion required to prevent shear failure, C_min, as a proxy for the potential for fracturing or reactivation of existing fractures. We compare depth profiles of C_min with structural analyses of image logs from four boreholes located on the valley floor, and find that fracture abundance declines sharply with depth in the uppermost 15 m of the bedrock, consistent with the modeled profile of C_min. In contrast, C_min increases with depth at comparable depths below ridgetops, suggesting that ridgetop fracture abundance patterns may differ if topographic stresses are indeed important. Thus, the present results are consistent with the hypothesis that topography can influence subsurface rock fracture patterns and provide a basis for further observational tests. Copyright © 2014 John Wiley & Sons, Ltd.     

Break-up spots: Could the Pacific open as a consequence of plate kinematics?

The South Central Pacific is the location of an abnormal concentration of intraplate volcanism. Noting that this volcanism is present from the Kermadec Tonga trench to the Easter microplate and forms a wide east–west channel, we propose to explain its occurrence in relation to the Pacific plate geometry and kinematics. We construct 2D numerical models of stress and strain within the Pacific plate using its velocity field and boundary conditions. The models indicate a shear band, associated to a change from compressional stresses to the south to tensional stresses to the north, which develop after 10 Myr between the Australian plate corner and the Easter microplate. We propose that the Central Pacific intraplate volcanism is related to this process, and may represent the first step of a future plate re-organization which will eventually break the Pacific plate in a southern and a northern plate due to intraplate stresses. Present-day intraplate volcanism would define break up spots of the future border.     

On the origin of compressional intraplate stresses in South America

Intraplate stresses in middle South America are not negligible. We report thrust-faulting mechanisms for five intraplate earthquakes, which indicate a dominant horizontal deviatoric compressional stress oriented in a NW-SE direction. We conclude that this state of stress is due to forces connected with spreading on the Mid Atlantic Ridge and resistive forces exerted by the Caribbean plate to the north and the Nazca plate to the west. The existence and nature of the resistive forces is inferred from earthquake mechanisms and geological evidence presented in other studies. All the available intraplate stress data for Nazca and South America indicate that both plates are under deviatoric compression generated at spreading centers. The absence of tensional earthquake focal mechanisms, particularly in the Nazca plate near the trench, suggests that the forces associated with the gravitational sinking of subducted lithosphere are locally compensated. We present a simple numerical calculation of a non-subducting plate to show how the compressional deviatoric stresses in middle South America can be used to estimate an upper bound of about 10²¹ P for the viscosity of the mantle.     

Role of interannual equatorial forcing on the subsurface temperature dipole in the Bay of Bengal during IOD and ENSO events

Role of equatorial forcing on the thermocline variability in the Bay of Bengal (BoB) during positive and negative phases of the Indian Ocean Dipole (IOD) and El Niño Southern Oscillation (ENSO) was investigated using the Regional Ocean Modeling System (ROMS) simulations during 1988 to 2015. Two numerical experiments were carried out for (i) the Indian Ocean Model (IOM) with interannual open boundary conditions and (ii) the BoB Model (BoBM) with climatological boundary conditions. The first mode of Sea Surface Height Anomalies (SSHA) variability showed a west-east dipole nature in both IOM and altimetry observations around 11°N, which was absent in the BoBM. The vertical section of temperature along the same latitude showed a sharp subsurface temperature dipole with a core at ~ 100 m depth. The positive (negative) subsurface temperature anomalies were observed over the whole northeastern BoB during NIOD (PIOD) and LN (EN) composites due to stronger (weaker) second downwelling Kelvin Waves. During the negative phases of IOD and ENSO, the cyclonic eddy on the southwestern BoB strengthened due to intensified southward coastal current along the western BoB and local wind stress. The subsurface temperature dipole was at its peak during October–December (OND) with 1-month lag from IOD and was evident from the Argo observations and other reanalysis datasets as well. A new BoB dipole index (BDI) was defined as the normalized difference of 100-m temperature anomaly and found to be closely related to the frequency of cyclones and the surface chlorophyll-a concentration in the BoB.

     

Stress pattern in two seismogenic sources in Nepal-Himalaya and its vicinity

The composite stereographic projection of orientations of the compression and tension axes using thirty-nine fault-plane solutions of earthquakes from two active seismogenic sources of Nepal and adjoining areas were examined and the nature of stress pattern and their influence on tectonics in the region have been studied. The seismogenic source located in Eastern Nepal region, which has been the site of 1934 Bihar-Nepal great earthquake of M 8.4, is presently experiencing N-S to NE-SW directed compressive stresses. The inferred pattern of compression axes in Western Nepal region suggests a shallow compressive stress, dipping N-S to NE-SW. Approximately similar nature of the stress regime is observed in Western and Eastern regions of Nepal, separated by nearly 700 km; it shows N-S to NNE-SSW direction of compression and underthrusting of the Indian Plate beneath the Himalaya at a shallow angle. Present study indicates that the stress is being released along the strikes of some of the transverse faults present in the region since the compressive stress exerted by the northward movement of the Indian Plate is approximately perpendicular to the Himalayan collision belt. Unilateral stress pattern generated by the northward movement of the Indian Plate in the central part of the Himalaya reveals that the present day collision occurs roughly perpendicular to the local strike of the Himalaya.     

Fractal analysis of earthquake swarms of Vogtland/NW-Bohemia intraplate seismicity

The special type of intraplate microseismicity with swarm-like occurrence of earthquakes within the Vogtland/NW-Bohemian Region is analysed to reveal the nature and the origin of the seismogenic regime. The long-term data set of continuous seismic monitoring since 1962, including more than 26000 events within a range of about 5 units of local magnitude, provides an unique database for statistical investigations. Most earthquakes occur in narrow hypocentral volumes (clusters) within the lower part of the upper crust, but also single event occurrence outside of spatial clusters is observed. Temporal distribution of events is concentrated in clusters (swarms), which last some days until few month in dependence of intensity. Since 1962 three strong swarms occurred (1962, 1985/86, 2000), including two seismic cycles. Spatial clusters are distributed along a fault system of regional extension (Leipzig-Regensburger Störung), which is supposed to act as the joint tectonic fracture zone for the whole seismogenic region. Seismicity is analysed by fractal analysis, suggesting a unifractal behaviour of seismicity and uniform character of seismotectonic regime for the whole region. A tendency of decreasing fractal dimension values is observed for temporal distribution of earthquakes, indicating an increasing degree of temporal clustering from swarm to swarm. Following the idea of earthquake triggering by magma intrusions and related fluid and gas release into the tectonically pre-stressed parts of the crust, a steady increased intensity of intrusion and/or fluid and gas release might account for that observation. Additionally, seismic parameters for Vogtland/NW-Bohemia intraplate seismicity are compared with an adequate data set of mining-induced seismicity in a nearby mine of Lubin/Poland and with synthetic data sets to evaluate parameter estimation. Due to different seismogenic regime of tectonic and induced seismicity, significant differences between b-values and temporal dimension values are observed. Most significant for intraplate seismicity are relatively low fractal dimension values for temporal distribution. That observation reflects the strong degree of temporal earthquake clustering, which might explain the episodic character of earthquake swarms and support the idea of push-like triggering of earthquake avalanches by intruding magma.     

2016年11月13日新西兰凯库拉(Kaikoura)MW7.8地震同震位移和应力数值分析

2016年11月13日新西兰南岛北端凯库拉（Kaikoura）发生了M_W7.8大地震,造成了强烈的地表变形并引发大面积滑坡和海啸的发生.基于美国地质调查局（USGS）断层滑动模型,建立全球同震横向不均匀并行椭球型地球模型,计算了此次新西兰凯库拉大地震产生的同震形变和应力及库仑应力变化.初步计算结果表明：新西兰凯库拉M_W7.8地震造成断层上盘东北向抬升,下盘西南俯冲;引起发震区域同震位移较大,从凯库拉到坎贝拉（Campbell）以及首都惠灵顿（Wellington）整体上东北向抬升,最大同震水平位移1.2 m,垂直位移1.1 m.此次大地震释放了发震断层上积累的压应力,但增加了发震断层两端的挤压力;同时,同震剪应力变化增加了NE-SW向断层发生右旋滑动的危险性;采用此次地震发震断层参数计算得出的最大库仑应力变化增加区域集中在发震断层两端,可达到MPa量级.当分别采用新西兰北岛Awatere断裂系和南岛Wellington断裂系参数计算库仑应力变化时,发现新西兰北岛和南岛震中以南区域的库仑应力均增加,可触发部分余震的发生.     

A comparison of surface and underground array measurements of ambient noise recorded in Naples (Italy)

In this study, we describe two experiments of seismic noise measurements carried out in Naples, Italy. The site allowed measurements to be obtained both at the surface and in a tunnel that is 120-m-deep. The main goal was to compare the seismic response evaluated at the surface to the in-tunnel response, through spectral, polarization, and resonance directivity analyses. In the 1 to 20 Hz frequency band, the noise level was up to 15 dB higher at the surface than in the tunnel. The polarization properties and horizontal-to-vertical spectral ratios appear not to be influenced by the tunnel geometry or by the topography. Some preferential alignments were observed in the polarization azimuths computed at the surface, which are likely to be due to local sources, rather than morphological features. The absence of directivity effects and the low noise levels in the tunnel make this site suitable for installing seismic stations. We also studied how the subsoil structure affects the seismic motion at the surface. The dispersive properties of the Rayleigh waves were investigated using the spatial autocorrelation method. A joint inversion of the dispersion data and the horizontal-to-vertical spectral ratios provided the subsurface Vs profile. The derived model has a low velocity contrast at depth, such as to generate moderate and broad H/V spectral ratio peak amplitude. The normalized spectral ratio appears more appropriate to identify the soil-resonance frequencies.     

Modelling tide–surge interaction effects using finite volume and finite element models of the Irish Sea

An unstructured mesh model of the west coast of Britain, covering the same domain and using topography and open boundary forcing that are identical to a previous validated uniform grid finite difference model of the region, is used to compare the performance of a finite volume (FV) and a finite element (FE) model of the area in determining tide–surge interaction in the region. Initial calculations show that although qualitatively both models give comparable tidal solutions in the region, comparison with observations shows that the FV model tends to under-estimate tidal amplitudes and hence background tidal friction in the eastern Irish Sea. Storm surge elevations in the eastern Irish Sea due to westerly, northerly and southerly uniform wind stresses computed with the FV model tend to be slightly higher than those computed with the FE model, due to differences in background tidal friction. However, both models showed comparable non-linear tide–surge interaction effects for all wind directions, suggesting that they can reproduce the extensive tide–surge interaction processes that occur in the eastern Irish Sea. Following on from this model comparison study, the physical processes contributing to surge generation and tide–surge interaction in the region are examined. Calculations are performed with uniform wind stresses from a range of directions, and the balance of various terms in the hydrodynamic equations is examined. A detailed comparison of the spatial variability of time series of non-linear bottom friction and non-linear momentum advection terms at six adjacent nodes at two locations in water depths of 20 and 6 m showed some spatial variability from one node to another. This suggests that even in the near coastal region, where water depths are of the order of 6 m and the mesh is fine (of order 0.5 km), there is significant spatial variability in the non-linear terms. In addition, distributions of maximum bed stress due to tides and wind forcing in nearshore regions show appreciable spatial variability. This suggests that intensive measurement campaigns and very high-resolution mesh models are required to validate and reproduce the non-linear processes that occur in these regions and to predict extreme bed stresses that can give rise to sediment movement. High-resolution meshes will also be required in pollution transport problems.     

A model study of spin-down and circulation in a cold water dome

A three-dimensional prognostic hydrodynamic model in cross sectional form is used to examine the influence of bottom friction, mixing and topography upon the spin-down and steady-state circulation in a cold water bottom-dome. Parameters characteristic of the Irish Sea or Yellow Sea cold water domes are used. In all calculations, motion is induced by specifying an initial temperature distribution characteristic of the dome, and an associated along frontal flow. The spin-down of the dome is found to be influenced by the coefficient of bottom friction, with a typical time scale of order 10 days, and in general to be independent of the chosen initial vertical profile of along frontal flow. However, in the case in which the along frontal flow is such that the near bed velocity is zero, then bottom stress is also zero, and there is no appreciable spin-down. Calculations showed that the formulation of viscosity and diffusivity had a greater effect upon the steady-state circulation than topography, suggesting that background mixing of tidal origin is important. The lack of topographic influence was due mainly to the formulation of the initial conditions which were taken to be independent of topography. The steady-state circulation was characterized by a cyclonic flow in the surface region, with an anti-cyclonic current near the bed, where frictional effects produced a bottom Ekman layer and an across frontal flow. This gave rise to vertical circulation cells in the frontal region of the dome with prevailing downwelling motion inside the dome. A detailed analysis of the dynamic balance of the various terms in the hydrodynamic equations yielded insight into the processes controlling the steady-state circulation in cold water domes. Responsible Editor: Phil Dyke     

Refined models of gravitational potential energy compared with stress and strain rate patterns in Iberia

This study examines the role of gravitational potential energy (GPE) in generating second-order (spatial scale ∼10² km) variations in the Iberia stress and strain-rate patterns. We present a new map of present-day strain rate field derived from the secular velocity field computed using all available continuously operating Global Navigation Satellite Systems (GNSS) stations in Iberia. The estimated strain rate field is generally consistent with the tectonic framework of the Iberian region, even though sporadic sharp local variations downgrade its correlation with the regional stress patterns. Many of the sharp spatial variations in the strain rate map are consistent with local changes of deformation style determined by prevailing faults. To obtain a more accurate estimate of GPE we use new data on the structure of the crust and apply a thin sheet approach using a 3-D definition of deviatoric stress. The GPE is derived from two isostatically compensated models (GPEd and GPEe compensated by density and elevation adjustment, respectively) and from the truncated geoid (GPEg). The GPE stresses are then summed with the first-order stress field due to the Eurasia–Nubia (EU–NU) convergence and the results compared with both the stress and strain rate data. In agreement with previous studies, we find that the GPE does not significantly change the NW–SE average direction of the most compressive stress (SHmax) imposed by the EU–NU collision, its main effect being to cause spatially changing stress regimes. From the analysis of the different GPE models we find: (1) in the Pyrenees, the tectonic forces have a secondary role when compared to the GPE. In this region, the model that best correlates with observations is the one emphasizing the role of surface elevation as a source of GPE (GPEe); (2) in the Iberian Chain and the Betics, the GPE imposes NE–SW extension consistent with a strike-slip regime and is equally (GPEe) or more (GPEg) important than the tectonic forces. In these regions, both deep heterogeneities associated with mantle convection and elevation are important sources of GPE; (3) in western Iberia, the GPE differences work against dominant tectonic forces by reducing the SHmax magnitude. The GPEg model is the one that best predicts the average strike-slip regime in Galicia; and finally (4) in the Gulf of Cadiz the gravitational potential stresses have a minor role and the style of deformation is clearly controlled by the tectonic forces.     

Earthquake source mechanisms from body-waveform inversion and intraplate tectonics in the northern Indian Ocean

Double-couple point-source parameters for 11 of the largest intraplate earthquakes in the northern Indian Ocean during the last 20 y were determined from a formal inversion of long-period P and SH waveforms. Nine of the events have centroid depths at least 17 km below the seafloor, well into the upper mantle; two have centroid depths as great as 39 km. Using the source mechanisms of these earthquakes, we distinguish two major intraplate tectonic provinces in the northern Indian Ocean. To the west of the Ninetyeast Ridge, in the southern Bay of Bengal, intraplate earthquakes have thrust-faulting mechanisms with P axes oriented N-S. The centroid depths of these earthquakes range from 27 to 39 km below the seafloor. Lithospheric shortening in this region is thus accomplished by thrust faulting in the strong core of the oceanic upper mantle, while other geophysical evidence suggests that shallow sedimentary and crustal layers apparently deform predominantly by folding. In the immediate vicinity of the Ninetyeast Ridge, earthquakes display strike-slip mechanisms with left-lateral motion on planes parallel to the ridge. This type of faulting occurs from at least 10°S to the northern end of the Ninetyeast Ridge near 10°N, where the ridge meets the Sunda Arc. Seismic activity diminishes to the east of the Ninetyeast Ridge, but is also characterized by strike-slip faulting. Despite these variations in deformational style, the inferred orientation of greatest compressive stress in the northern Indian Ocean displays a consistent long-wavelength pattern over a large portion of the Indian plate, varying smoothly from nearly N-S in the Bay of Bengal to NW-SE in the northeastern Indian Ocean. This plate-wide stress pattern and the high level of intraplate seismicity in the northern Indian Ocean are likely the results of substantial resistance, along the Himalayan continental collision zone, to the continued northward motion of the western portion of the Indian plate. Oceanic intraplate earthquakes in other regions, where the level of deviatoric stress associated with the long-wavelength part of the stress field is likely to be smaller, need not be comparably reliable indicators of the plate-wide stress field.     

The pattern of the paleo- and present-day stresses of Northern Tien Shan

The study of neotectonic stresses of the North Tien Shan is carried out within the northern slopes of the Kyrgyz Ridge bordering the Chyua submontane depression. The purpose of the study was to reconstruct the neotectonic stresses from the geological indicators and to compare the obtained results with the data on the present-day stress state of this region. The work is based on the analysis of the field in situ data acquired in the expeditions in 2009 and 2011. For the first time for the Northern Tien Shan region, the general (averaged) neotectonic stresses, which are distinct between the uplifts and depressions, are reconstructed. The deformation of the positive topographic features at the recent tectonic stage takes place in the conditions of thrust faulting with the meriodional horizontal compression axis and subvertical tension axis, while deformation of the negative landforms is dominated by normal faulting with the vertical compression axis and northnortheast trending subhorizontal tension axis. Based on the field data, separate sites with different types of local stress tensors, which are determined by the Lode-Nadai parameter, are revealed. The neotectonic geodynamical regime within the depressions and uplifts and the variations in the type of the stress tensor agree with the characteristics of the present-day stress state reconstructed from seismological data.     

On the nature of intraplate earthquakes

Continental intraplate regions are characterized by uniform stresses over thousands of kilometers. Local stresses, with wavelengths of tens to hundreds of kilometers can accumulate at inhomogeneities lying within these regional fields. A variety of geological structures, herein called local stress concentrators (LSCs), act as elastic inhomogeneities. The temporal buildup of stress depends on the particular structure and its geometrical relationship with the regional stress field. The interaction of the local and the regional stress fields can result in the rotation of the latter over wavelengths of tens to hundreds of kilometers. This rotation can be detected by direct measurement or from seismicity data. Intraplate earthquakes (IPEs) result when the local stresses become comparable with their regional counterparts, i.e., hundreds of megapascals. Globally, most of the seismic energy release associated with IPEs occurs within old rifts which contain LSCs most favorable for stress buildup by stress inversion. Of the various LSCs, stepover en echelon faults are associated the largest IPEs. In low tectonic strain rate regions, IPEs are associated with larger stress drops. With the availability of a variety of LSCs, there is generally an absence of repeat earthquakes. Instead, successive earthquakes occur on different structures, leading to the observation of “roaming” earthquakes. These observations suggest a need for a reevaluation of seismic hazard estimation techniques. This study addresses some of these facets of the nature of IPEs with global examples, including a unique, detailed seismicity and geodetic data set collected in a dozen years following the 2001 M 7.7 Bhuj earthquake in western India.