Analysis of Izmit aftershocks 25 days before the November 12th 1999 Düzce earthquake, Turkey

We investigate spatial clustering of 2414 aftershocks along the Izmit M_w = 7.4 August 17, 1999 earthquake rupture zone. 25 days prior to the Düzce earthquake M_w = 7.2 (November 12, 1999), we analyze two spatial clusters, namely Sakarya (SC) and Karadere–Düzce (KDC). We determine the earthquake frequency–magnitude distribution (b-value) for both clusters. We find two high b-value zones in SC and one high b-value zone in KDC which are in agreement with large coseismic surface displacements along the Izmit rupture. The b-values are significantly lower at the eastern end of the Izmit rupture where the Düzce mainshock occurred. These low b-values at depth are correlated with low postseismic slip rate and positive Coloumb stress change along KDC. Since low b-values are hypothesized with high stress levels, we propose that at the depth of the Düzce hypocenter (12.5 km), earthquakes are triggered at higher stresses compared to shallower crustal earthquake. The decrease in b-value from the Karadere segment towards the Düzce Basin supports this low b-value high stress hypothesis at the eastern end of the Izmit rupture. Consequently, we detect three asperity regions which are correlated with high b-value zones along the Izmit rupture. According to aftershock distribution the half of the Düzce fault segment was active before the 12 November 1999 Düzce mainshock. This part is correlated with low b-values which mean high stress concentration in the Düzce Basin. This high density aftershock activity presumably helped to trigger the Düzce event (M_w = 7.2) after the Izmit M_w 7.4 mainshock.     

The rupture process during the 1999 Düzce, Turkey, earthquake from joint inversion of teleseismic and strong-motion data

Teleseismic and strong-motion data are inverted to determine the rupture process during the November 1999 Düzce earthquake in NW Turkey. The fault geometry, rise time and rupture velocity are determined from the aftershock distribution and preliminary inversions of the teleseismic data. Joint inversion of the teleseismic and strong-motion data is then carried out for the slip distribution. We obtain the strike 264°, dip 64°, rake −172°, seismic moment 5.0×10¹⁹ N m (M_w 7.1), and average stress drop 7 MPa. This earthquake was characterized by bilateral fault rupture and asymmetric slip distribution. Two asperities (areas of large slip) are identified, the eastern one being 1.5 times larger than the western one. The derived slip distribution is consistent with the aftershock distribution, surface rupture and damage. The point of rupture initiation in this Düzce earthquake coincided with the eastern tip of the aftershock distribution of the August 1999 Izmit earthquake.     

Poroelastic relaxation and aftershocks of the 2001 Bhuj earthquake, India

We analyse aftershocks of the 26 January 2001 Bhuj earthquake, India, that were recorded for 10 weeks following the mainshock. We calculate undrained or instantaneous pore pressure and change in Coulomb stress due to the earthquake and their poroelastic relaxation in the following 10 weeks period. Almost all aftershocks occurred in the region of coseismic dilatation. In the subsequent period, pore pressure increased through relaxation in the dilatation region which further modified coseismic Coulomb stress. Maximum increase in pore pressure is estimated to be about 0.7 MPa in 60 days time following the mainshock. Correlation between the zones of increased pore pressure and postseismic Coulomb stress with that of aftershocks, suggests a definite role of fluid diffusion in their delayed triggering.     

昆仑山大地震震后形变反映的地壳岩石流变特性

2001年昆仑山地震是我国近50年来最大的地震。用3种模型对震后中国地震局跨断层布设的4个GPS站点记录到显著的震后形变进行了模拟研究。结果表明：单纯的上地壳10 km为滞弹性的模型不能解释震后形变的幅度;30 km的弹性上地壳覆盖在 40 km 的柔性下地壳上的松弛模型可以解释变形速率指数衰减的主要特征;而两个模型的结合不但能解释整 体指数衰减的特征,而且还能更好地拟合震后最初几周的较高形变速率。结果表明,在昆仑山断层两侧存在着流变性质的 差异。     

Elastic and visco-elastic stress triggering in the South Iceland Seismic Zone due to large earthquakes since 1706

Damaging earthquakes in the South Iceland Seismic Zone (SISZ) occur fairly regularly and often as a series of events with a few days only between individual events. Tolerably reliable information on epicentre locations and mechanisms are available for 13 M ≥ 6 events between 1706 and 2000. For these events, we computed the co- and post-seismic stress fields, hereby approximating the SISZ by a mixed elastic/visco-elastic layered half-space. The horizontal shear stress and the Coulomb stress changes were analysed to detect possible trigger mechanisms, which may aid future earthquake mitigation efforts. We tested several criteria but must conclude that the start of an earthquake series in the SISZ cannot be explained by triggering through previous events. Inside an individual series, however, one may infer triggering. Our results are in contradiction with the findings in other regions of the world. The reason might be related to the fact that the SISZ is not a mature fault zone, in which old faults are re-activated if a certain stress level threshold is passed. In addition, uncertainties in the model parameters as well as the neglect of horizontal variations in the model and of possible stress transfer due to volcanic activity further complicate the evaluation of our results and need to be taken into account in future studies.     

Seismic behavior in central Taiwan: Response to stress evolution following the 1999 Mw 7.6 Chi-Chi earthquake

Following the 1999 M_w 7.6 Chi-Chi earthquake, a large amount of seismicity occurred in the Nantou region of central Taiwan. Among the seismic activities, eight M_w ⩾ 5.8 earthquakes took place following the Chi-Chi earthquake, whereas only four earthquakes with comparable magnitudes took place from 1900 to 1998. Since the seismicity rate during the Chi-Chi postseismic period has never returned to the background level, such seismicity activation cannot simply be attributed to modified Omori’s Law decay. In this work, we attempted to associate seismic activities with stress evolution. Based on our work, it appears that the spatial distribution of the consequent seismicity can be associated with increasing coseismic stress. On the contrary, the stress changes imparted by the afterslip; lower crust–upper mantle viscoelastic relaxation; and sequent events resulted in a stress drop in most of the study region. Understanding seismogenic mechanisms in terms of stress evolution would be beneficial to seismic hazard mitigation.     

Mantle flow-induced crustal thinning in the area between the easternmost part of the Anatolian plate and the Arabian Foreland (E Turkey) deduced from the geological and geophysical data

Eastern Anatolia consisting of an amalgamation of fragments of oceanic and continental lithosphere is a current active intercontinental contractional zone that is still being squeezed and shortened between the Arabian and Eurasian plates. This collisional and contractional zone is being accompanied by the tectonic escape of most of the Anatolian plate to the west by major strike-slip faulting on the right-lateral North Anatolian Transform Fault Zone (NATFZ) and left-lateral East Anatolian Transform Fault Zone (EATFZ) which meet at Karlıova forming an east-pointing cusp. The present-day crust in the area between the easternmost part of the Anatolian plate and the Arabian Foreland gets thinner from north (ca 44 km) to south (ca 36 km) relative to its eastern (EAHP) and western sides (central Anatolian region). This thinner crustal area is characterized by shallow CPD (12–16 km), very low Pn velocities (< 7.8 km/s) and high Sn attenuation which indicate partially molten to eroded mantle lid or occurrence of asthenospheric mantle beneath the crust. Northernmost margin of the Arabian Foreland in the south of the Bitlis–Pötürge metamorphic gap area is represented by moderate CPD (16–18 km) relative to its eastern and western sides, and low Pn velocities (8 km/s). We infer from the geophysical data that the lithospheric mantle gets thinner towards the Bitlis–Pötürge metamorphic gap area in the northern margin of the Arabian Foreland which has been most probably caused by mechanical removal of the lithospheric mantle during mantle invasion to the north following the slab breakoff beneath the Bitlis–Pötürge Suture Zone. Mantle flow-driven rapid extrusion and counterclockwise rotation of the Anatolian plate gave rise to stretching and hence crustal thinning in the area between the easternmost part of the Anatolian plate and the Arabian Foreland which is currently dominated by wrench tectonics.     

Co and postseismic characteristics of Indian sub-continent in response to the 2004 Sumatra earthquake

The M_w 9.3 Sumatra earthquake of December 26, 2004 caused extensive coseismic displacements globally, measurements of which were made essentially using modern geodetic techniques. This earthquake induced considerable perturbation in stress distribution as far as ∼8000 km away from the epicenteral region, which is tending to relax to its normal rates as seen from postseismic transient deformation. The monitoring of crustal displacements from strategically located sites using GPS provides coseismic as well as postseismic deformation that facilitates the understanding of the fault geometry, elastic thickness, postseismic relaxation mechanisms, rheology and earthquake recurrence time interval.We investigated coseismic and postseismic GPS derived displacements in Indian region together with the GPS data collected from Andaman and Sumatra region. It is found that while EW displacements are significantly large in peninsular India, those in the region to the north of Central India Tectonic Zone (CITZ) are relatively small. We could delineate the postseismic transients from position time series and interpreted them in terms of viscoelastic relaxation. It is inferred that the postseismic deformation is characterized by a power-law viscoelastic flow in the mantle. In Indian peninsula region, the timescale parameter of the exponential decay (τ = 250 days) would require an extremely low viscosity for the upper mantle. Relying on the prevailing coseismic and postseismic displacement fields, the present study also reflects upon the contemporary litho-tectonics of the Indian sub-continent.     

The North Maladeta Fault (Spanish Central Pyrenees) as the Vielha 1923 earthquake seismic source: Recent activity revealed by geomorphological and geophysical research

The Spanish Central Pyrenees have been the scenario of at least two damaging earthquakes in the last 800 years. Analysis of macroseismic data of the most recent one, the Vielha earthquake (19 November 1923), has led to the identification of the North Maladeta Fault (NMF) as the seismic source of the event. This E–W trending fault defines the northern boundary of the Maladeta Batholith and corresponds to a segment of the Alpine Gavarnie thrust fault. Our study shows that the NMF offsets a reference Neogene peneplain. The maximum observed vertical displacement is  730 m, with the northern downthrown sector slightly tilting towards the South. This offset provides evidence of normal faulting and together with the presence of tectonic faceted spurs allowed us to geomorphically identify a fault trace of 17.5 km. This length suggests that a maximum earthquake of M_w = 6.5 ± 0.66 could occur in the area. The geomorphological study was improved with a resistivity model obtained at Prüedo, where a unique detritic Late Miocene sequence crops out adjacent to the NMF. The section is made up of 13 audiomagnetotelluric soundings along a 1.5 km transect perpendicular to the fault trace at Prüedo and reveals the structure in depth, allowing us to interpret the Late Miocene deposits as tectonically trapped basin deposits associated with normal faulting of the NMF. The indirect age of these deposits has been constrained between 11.1 and 8.7 Ma, which represents a minimum age for the elevated Pyrenean peneplain in this part of the Pyrenees. Therefore, we propose the maximum vertical dip-slip rate for the NMF to be between 0.06 and 0.08 mm/a. Normal faulting in this area is attributed to the vertical lithospheric stress associated with the thickened Pyrenean crust.     

汶川地震的发生对周围断层稳定性影响的数值模拟

大地震发生后, 研究地震的发生对周围断层的影响尤为重要.利用川西-藏东地区三维粘弹性有限元模型,考虑地表高程和粘弹性松弛等因素的影响,研究主要断裂带库仑应力累积速率和汶川地震的发生对周围断层的影响.结果表明:(1)龙门山断裂带年累积速率为0.28×10-3~0.35×10-3 MPa/a,这种较小的累积速率与龙门山断裂带强震较长复发间隔一致;(2)汶川地震的发生除造成震源区应力减小外, 还造成断裂带北东段不同程度的应力增加, 这与震后余震的分布基本吻合;(3)鲜水河断裂北西段、东昆仑、龙日坝、岷江以及虎牙断裂库仑应力水平增加显著,且汶川地震对于玉树地震的发生有微弱的加载效应;(4)汶川地震的发生造成鲜水河断裂带强震复发间隔缩短约52~104 a,是值得关注的强震危险区.      

两次于田M_S7.3地震间应力触发作用及2014年于田地震的发生对周缘断层的影响

基于青藏高原及邻区的三维粘弹性有限元模型,讨论2008年于田MS7.3级地震与2014年于田MS7.3级地震之间的关系,并研究2014年于田MS7.3级地震的发生造成周围断层的库仑破裂应力变化。初步结果表明:1)2008年于田MS7.3级地震在2014年于田MS7.3级地震震中滑动方向上产生的库仑破裂应力变化高于地震触发的阈值0.01 MPa,存在明显的触发作用。在视摩擦系数分别取0.4和0.6时,震源区同震库仑破裂应力变化为0.0167 MPa和0.0170 MPa;而考虑粘弹性松弛作用时产生的库仑应力增加量分别为0.0187 MPa和0.0194 MPa。结合断裂带构造应力年累计速率的结果,2008年于田地震的发生造成2014年于田地震提前21.4~24.9 a;2)在较短的时间尺度内,对于距离相近的两次地震之间,同震产生的应力变化远大于粘弹性松弛效应产生的变化;3)2014年于田MS7.3级地震的发生造成阿尔金断裂中北段、玛尼—玉树断裂中段、东昆仑断裂西段、柴达木北缘断裂东段、西秦岭北缘断裂西段等不同程度的加载效应,地震危险性有所增强。其中阿尔金断裂中段库仑应力增加最为明显,最大达2.8×10–3 MPa;玛尼—玉树断裂中段次之,应力增加量最大达5.6×10–4 MPa;东昆仑断裂西段应力增加量最大达4.75×10–4 MPa。而玛尼—玉树断裂西段库仑破裂应力最大卸载量达3.6×10–3 MPa。     

Influence of the precollisional stage on subduction dynamics and the buried crust thermal state: Insights from numerical simulations

At continental subduction initiation, the continental crust buoyancy may induce, first, a convergence slowdown, and second, a compressive stress increase that could lead to the forearc lithosphere rupture. Both processes could influence the slab surface P–T conditions, favoring on one side crust partial melting or on the opposite the formation of ultra-high pressure/low temperature (UHP-LT) mineral. We quantify these two effects by performing numerical simulations of subduction. Water transfers are computed as a function of slab dehydration/overlying mantle hydration reactions, and a strength decrease is imposed for hydrated mantle rocks. The model starts with an old oceanic plate ( 100 Ma) subducting for 145.5 Myr with a 5 cm/yr convergence rate. The arc lithosphere is thermally thinned between 100 km and 310 km away from the trench, due to small-scale convection occuring in the water-saturated mantle wedge. We test the influence of convergence slowdown by carrying on subduction with a decreased convergence rate (≤ 2 cm/yr). Surprisingly, the subduction slowdown yields not only a strong slab warming at great depth (> 80 km), but also a significant cooling of the forearc lithosphere at shallower depth. The convergence slowdown increases the subducted crust temperature at 90 km depth to 705 ± 62 °C, depending on the convergence rate reduction, and might thus favor the oceanic crust partial melting in presence of water. For subduction velocities ≤ 1 cm/yr, slab breakoff is triggered 20–32 Myr after slowdown onset, due to a drastic slab thermal weakening in the vicinity of the interplate plane base. At last, the rupture of the weakened forearc is simulated by imposing in the thinnest part of the overlying lithosphere a dipping weakness plane. For convergence with rates ≥ 1 cm/yr, the thinned forearc first shortens, then starts subducting along the slab surface. The forearc lithosphere subduction stops the slab surface warming by hot asthenosphere corner flow, and decreases in a first stage the slab surface temperature to 630 ± 20 °C at 80 km depth, in agreement with P–T range inferred from natural records of UHP-LT metamorphism. The subducted crust temperature is further reduced to 405 ± 10 °C for the crust directly buried below the subducting forearc. Such a cold thermal state at great depth has never been sampled in collision zones, suggesting that forearc subduction might not be always required to explain UHP-LT metamorphsim.     

A study of microearthquake seismicity and focal mechanisms within the Sea of Marmara (NW Turkey) using ocean bottom seismometers (OBSs)

We have carried out seismological observations within the Sea of Marmara (NW Turkey) in order to investigate the seismicity induced after Gölcük–İzmit (Kocaeli) earthquake (M_w 7.4) of August 17, 1999, using ocean bottom seismometers (OBSs). High-resolution hypocenters and focal mechanisms of microearthquakes have been investigated during this Marmara Sea OBS project involving deployment of 10 OBSs within the Çınarcık (eastern Marmara Sea) and Central-Tekirdağ (western Marmara Sea) basins during April–July 2000. Little was known about microearthquake activity and their source mechanisms in the Marmara Sea. We have detected numerous microearthquakes within the main basins of the Sea of Marmara along the imaged strands of the North Anatolian Fault (NAF). We obtained more than 350 well-constrained hypocenters and nine composite focal mechanisms during 70 days of observation. Microseismicity mainly occurred along the Main Marmara Fault (MMF) in the Marmara Sea. There are a few events along the Southern Shelf. Seismic activity along the Main Marmara Fault is quite high, and focal depth distribution was shallower than 20 km along the western part of this fault, and shallower than 15 km along its eastern part. From high-resolution relative relocation studies of some of the microearthquake clusters, we suggest that the western Main Marmara Fault is subvertical and the eastern Main Marmara Fault dips to south at 45°. Composite focal mechanisms show a strike-slip regime on the western Main Marmara Fault and complex faulting (strike-slip and normal faulting) on the eastern Main Marmara Fault.     

Patterns in the Turkish reinsurance claim data, 1998–2002: An application of the extreme value theory

The devastating 1999 Marmara and Düzce earthquakes led to a significant increase in the earthquake studies in Turkey in geological, engineering and financial aspects. Extreme Value Theory (EVT) has a range of applications from stock market changes to natural disasters like floods and hurricanes. Here EVT is fitted to the ordinary and earthquake reinsurance claims of Turkey.     

从同震和震后形变分析1668年M8.5级郯城地震对周边地震活动性的影响

以山东郯城1668年大地震为例,以前人地表地质调查结果为约束,利用弹性位错理论初步获取了该地震的同震破裂模型;在此基础上,基于粘弹性分层模型分析了该地震的同震和震后形变,同时以主震断层为接收断层计算了库仑应力分布,进一步讨论了地幔不同粘滞性系数对地表形变和库仑应力变化的影响。计算结果显示,该地震是一个右旋走滑为主兼有一定逆冲性质的地震,其同震位移巨大,能量释放较彻底;同震破裂造成震中郯城县西北、东北和南部部分断层库仑应力增加,而震后形变使得这些断层库仑应力进一步增加,在单县、宿迁和日照等地,地震后350 a库仑应力变化量达到+1bar-+1MPa量级;地幔粘滞性系数不同,形变量和库仑应力变化达到稳定的时间不同,但最终趋于稳定的数值基本一致。     

A speculation on the structure of the D″ layer: The growth of anti-crust at the core–mantle boundary through the subduction history of the Earth

The growth curve of the continental crust shows that large amounts of continental crust formed in the early part of the Earth history are missing. In order to test a hypothesis that the former crust was subducted to the deep mantle, we performed phase assemblage analysis in the systems of mid-oceanic ridge basalt (MORB), anorthosite, and tonalite–trondhjemite–granite (TTG) down to the core–mantle boundary (CMB) conditions. Results show that all these materials can be subducted to the CMB leading to the development of a compositional layering in the D″ layer. We speculate that there could be five layers of FeO-enriched melt from partial melting of MORB, MORB crust, anorthosite, TTG, and slab or mantle peridotite in ascending order. Although the polymorphic transformation of perovskite to post-perovskite in (Mg,Fe)SiO₃ may explain the seismic discontinuity at the top of the D″ layer (D″ discontinuity), the effects of solid solution on the sharpness of the transformation suggest that the compositional layering is more plausible for the origin of the D″ discontinuity. The D″ layer can be an “anti-crust” made up mostly of TTG + anorthosite derived from the former continental crust. Tectonic style of the anti-crust at the CMB is similar to that at the surface. At both places, chemically distinct layers are density stratified and are also characterized by the processes of accretion, magmatism, and metasomatism.     

基于分裂节点法的地震同震和震后形变数值模拟及其在汶川大地震中的应用

大地震导致的同震及震后效应,对于分析不同地震之间的相互影响及区域地震危险性等有着重要的作用.文中开发了模拟地震同震及震后效应的三维黏弹性有限元程序,通过计算走滑断层震例(概念性模型)引起的同震及震后效应,并与解析/半解析解进行对比,验证了程序的可靠性.同时基于概念性模型,分析了不同介质参数对同震及震后的地表变形的影响....     

Large-strain deformation and strain partitioning in polyphase rocks: Dislocation creep of olivine–magnesiowüstite aggregates

Aggregates composed of olivine and magnesiowüstite have been deformed to large strains at high pressure and temperature to investigate stress and strain partitioning, phase segregation and possible localization of deformation in a polyphase material. Samples with 20 vol.% of natural olivine and 80 vol.% of (Mg_0.7Fe_0.3)O were synthesized and deformed in a gas-medium torsion apparatus at temperatures of 1127 °C and 1250 °C, a confining pressure of 300 MPa and constant angular displacement rates equivalent to constant shear strain rates of 1–3.3 × 10^− 4 s^− 1. The samples deformed homogeneously to total shear strains of up to γ  15. During constant strain rate measurements the flow stress remained approximately stable at 1250 °C while it progressively decreased after the initial yield stress at the lower temperature. Mechanical data, microstructures and textures indicate that both phases were deforming in the dislocation creep regime. The weaker component, magnesiowüstite, controlled the rheological behavior of the bulk material and accommodated most of the strain. Deformation and dynamic recrystallization lead to grain refinement and to textures that were not previously observed in pure magnesiowüstite and may have developed due to the presence of the second phase. At 1127 °C, olivine grains behaved as semi-rigid inclusions rotating in a viscous matrix. At 1250 °C, some olivine grains remained largely undeformed while deformation and recrystallization of other grains oriented for a-slip on (010) resulted in a weak foliation and a texture typical for pure dry olivine aggregates. Both a-slip and c-slip on (010) were activated in olivine even though the nominal stresses were up to 2 orders of magnitude lower than those needed to activate these slip systems in pure olivine at the same conditions.     

A composite geologic and seismic profile beneath the southern Rio Grande rift, New Mexico, based on xenolith mineralogy, temperature, and pressure

A complete understanding of the processes of crustal growth and recycling in the earth remains elusive, in part because data on rock composition at depth is scarce. Seismic velocities can provide additional information about lithospheric composition and structure, however, the relationship between velocity and rock type is not unique. The diverse xenolith suite from the Potrillo volcanic field in the southern Rio Grande rift, together with velocity models derived from reflection and refraction data in the area, offers an opportunity to place constraints on the composition of the crust and upper mantle from the surface to depths of  60 km. In this work, we calculate seismic velocities of crustal and mantle xenoliths using modal mineralogy, mineral compositions, pressure and temperature estimates, and elasticity data. The pressure, temperature, and velocity estimates from xenoliths are then combined with sonic logs and stratigraphy estimated from drill cores and surface geology to produce a geologic and velocity profile through the crust and upper mantle. Lower crustal xenoliths include garnet ± sillimanite granulite, two-pyroxene granulite, charnokite, and anorthosite. Metagabbro and amphibolite account for only a small fraction of the lower crustal xenoliths, suggesting that a basaltic underplate at the crust–mantle boundary is not present beneath the southern Rio Grande rift. Abundant mid-crustal felsic to mafic igneous xenoliths, however, suggest that plutonic rocks are common in the middle crust and were intraplated rather than underplated during the Cenozoic. Calculated velocities for garnet granulite are between  6.9 and 8.0 km/s, depending on garnet content. Granulites are strongly foliated and lineated and should be seismically anisotropic. These results suggest that velocities > 7.0 km/s and a layered structure, which are often attributed to underplated mafic rocks, can also be characteristic of alternating garnet-rich and garnet-poor metasedimentary rocks. Because the lower crust appears to be composed largely of metasedimentary granulite, which requires deep burial of upper crustal materials, we suggest the initial construction of the continental crust beneath the Potrillo volcanic field occurred by thickening of supracrustal material in the absence of large scale magmatic accretion. Mantle xenoliths include spinel lherzolite and harzburgite, dunite, and clinopyroxenite. Calculated P-wave velocities for peridotites range from 7.75 km/s to 7.89 km/s, with an average of 7.82 km/s. This velocity is in good agreement with refraction and reflection studies that report P_n velocities of 7.6–7.8 km/s throughout most of the Rio Grande rift. These calculations suggest that the low P_n velocities compared to average uppermost mantle are the result of relatively high temperatures and low pressures due to thin crust, as well as a fertile, Fe-rich, bulk upper mantle composition. Partial melt or metasomatic hydration of the mantle lithosphere are not needed to produce the observed P_n velocities.     

The large aftershocks triggered by the 2011 Mw 9.0 Tohoku-Oki earthquake,Japan

The Mw 9.0 Tohoku-Oki earthquake that occurred off the Pacific coast of Japan on March 11, 2011, was followed by thousands of aftershocks, both near the plate interface and in the crust of inland eastern Japan. In this paper, we report on two large, shallow crustal earthquakes that occurred near the Ibaraki-Fukushima prefecture border, where the background seismicity was low prior to the 2011 Tohoku-Oki earthquake. Using densely spaced geodetic observations (GPS and InSAR datasets), we found that two large aftershocks in the Iwaki and Kita-Ibarake regions (hereafter referred to as the Iwaki earthquake and the Kita-Ibarake earthquake) produced 2.1 m and 0.44 m of motion in the line-of-sight (LOS), respectively. The azimuth-offset method was used to obtain the preliminary location of the fault traces. The InSAR-based maximum offset and trace of the faults that produced the Iwaki earthquake are consistent with field observations. The fault location and geometry of these two earthquakes are constrained by a rectangular dislocation model in a multilayered elastic half-space, which indicates that the maximum slips for the two earthquakes are 3.28 m and 0.98 m, respectively. The Coulomb stress changes were calculated for the faults following the 2011 Mw 9.0 Tohoku-Oki earthquake based on the modeled slip along the fault planes. The resulting Coulomb stress changes indicate that the stresses on the faults increased by up to 1.1 MPa and 0.7 MPa in the Iwaki and Kita-Ibarake regions, respectively, suggesting that the Tohoku-Oki earthquake triggered the two aftershocks, supporting the results of seismic tomography.