Measures for groundwater security during and after the Hanshin-Awaji earthquake (1995) and the Great East Japan earthquake (2011), Japan

Many big earthquakes have occurred in the tectonic regions of the world, especially in Japan. Earthquakes often cause damage to crucial life services such as water, gas and electricity supply systems and even the sewage system in urban and rural areas. The most severe problem for people affected by earthquakes is access to water for their drinking/cooking and toilet flushing. Securing safe water for daily life in an earthquake emergency requires the establishment of countermeasures, especially in a mega city like Tokyo. This paper described some examples of groundwater use in earthquake emergencies, with reference to reports, books and newspapers published in Japan. The consensus is that groundwater, as a source of water, plays a major role in earthquake emergencies, especially where the accessibility of wells coincides with the emergency need. It is also important to introduce a registration system for citizen-owned and company wells that can form the basis of a cooperative during a disaster; such a registration system was implemented by many Japanese local governments after the Hanshin-Awaji Earthquake in 1995 and the Great East Japan Earthquake in 2011, and is one of the most effective countermeasures for groundwater use in an earthquake emergency. Emphasis is also placed the importance of establishing of a continuous monitoring system of groundwater conditions for both quantity and quality during non-emergency periods.     

Surface displacements in Japan before the 11 March 2011 M9.0 Tohoku-Oki earthquake

Daily resolution data retrieved from the 1243 ground-based Global Positioning System (GPS) stations in Japan are utilized to expose surface displacements before the destructive M9 Tohoku-Oki earthquake (March 11, 2011). Variations in the residual GPS data, in which effects of the long-term plate movements, short-term noise and frequency-dependent variations have been removed through a band-pass filter via the Hilbert–Huang transform, are compared with parameters of the focal mechanism associated with the Tohoku-Oki earthquake for validation. Analytical results show that the southward movements, which were deduced from the residual displacements and agree with the strike of the rupture fault, became evident on the 65th day before the Tohoku-Oki earthquake. This observation suggests that the shear stress played an important role in the seismic incubation period. The westward movements, which are consistent with the angle of the maximum horizontal compressive stress, covered entire Japan and formed an impeded area (142°E, 42°N) about 75 km away from the epicenter on the 47th day prior to the earthquake. The horizontal displacements integrated with the vertical movements from the residual GPS data are very useful to construct comprehensive images in diagnosing the surface deformation from destructive earthquakes along the subduction zone.     

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.     

3-D crustal-scale gravity model of the San Rafael Block and Payenia volcanic province in Mendoza,Argentina

The San Rafael Block(SRB)is part of one of the main retroarc volcanic provinces in southern Central Andes in Mendoza,Argentina.This block is located in the Andean foothills between the orogenic front and foreland basement uplifts of late Miocene age.In order to analyze the geochronological evolution of the Quaternary volcanism in the region,several geologic and geophysical studies have been conducted.Nevertheless,the crust,where the SRB is located,has not been well characterized yet.Based on gravimetric and magnetic data,together with isostatic and elastic thickness analyses,we modeled the crustal structure of the area.Information obtained has allowed us to understand the crust where the SRB and the Payenia volcanic province are located.Bouguer anomalies indicate that the SRB presents higher densities to the North of Cerro Nevado and Moho calculations suggest depths for this block between 40 and 50 km.Determinations of elastic thickness would indicate that the crust supporting the San Rafael Block presents values of approximately 10 km,being enough to support the block loading.However,in the Payenia region,elastic thickness values are close to zero due to the regional temperature increase.     

Long-traveling landslides in deep snow conditions induced by the 2011 Nagano Prefecture earthquake,Japan

Earthquakes in mountainous areas may produce many landslides that involve abundant snow, but few observations have been made of these hazardous phenomena. The 12 March 2011 north Nagano Prefecture earthquake (M_JMA 6.7) occurred in a mountainous part of Japan that typically has an annual snow cover of more than 2 m, and it induced many snowy landslides. Some of these traveled relatively long distances. We examined the snowy Tatsunokuchi landslide to reconstruct the landsliding processes over deep snow. We infer that the Tatsunokuchi landslide occurred by collapse of a rock debris mass of 5?×?10⁴ m³ that plunged into the abundant snow, forming a mixture of snow and rock debris, which then traveled on top of the snow. Later, the displaced mass included a large amount of snow which was pushed forward at the front and to the sides. The velocity of the landslide was estimated to be approximately 14 m/s. It appears that the displaced mass, having only a small proportion of rock debris, had a low enough density to travel easily on top of the snow. Our observations suggest that there was much liquid water at the base of the displaced mass shortly after the event. Our results suggest that landslides may damage wider areas than expected if they travel over deep snow.     

电阻率三维复杂结构的快速反演

运用共轭梯度迭代算法解三维反问题中的线性方程组,并结合求Jacobi矩阵 G 的Rodi算法,则每次反演迭代仅需一次正演计算,大大加快了计算速度,实现了直流电阻率三维快速反演。另外,由于避免了存储 G 和 G ^T G 所需的庞大存储量,以及在三维反演中加入光滑约束,有利于精细网格下的复杂模型反演。     

Explanation of liquefaction in after shock of the 2011 great east Japan earthquake using numerical analysis

During the 2011 Great East Japan Earthquake, severe liquefaction occurred in reclaimed ground in Urayasu city, Chiba prefecture. This liquefaction provided important lessons for us to re-recognize the liquefaction mechanism. A distinct feature of the liquefaction in this earthquake is that severe liquefaction happened not only in the main shock but also in an aftershock with a maximum acceleration of 25 gal. In some areas, liquefaction happened in the aftershock is even more serious than that happened in the main shock. In this paper, focus is placed on the characteristic features in the occurrence of liquefaction and consequent ground settlement. Based on the observed data, a series of dynamic–static analyses, considering not only the earthquake loading but also static loading during the consolidation after the earthquake shocks, are conducted in a sequential way just the same as the scenario in the earthquake. The calculation is conducted with 3D soil–water coupling finite element–finite difference analyses based on a cyclic elasto-plastic constitutive model. From the results of analyses, it is recognized that small sequential earthquakes, which cannot cause liquefaction of a ground in an independent earthquake vibration, cannot be neglected when the ground has already experienced liquefaction after a major vibration. In addition, the aftershock has great influence on the long-term settlement of low permeability soil layer. The observed and predicted liquefaction and settlements are compared and discussed carefully. It is confirmed that the numerical method used in this study can describe the ground behavior correctly under repeated earthquake shocks.     

First-order rupture features of the 2011 MW 9.0 Tohoku (Japan) earthquake from surface waves

This study analyzed the rupture directivity of the 2011 Tohoku earthquake by using 100-s Rayleigh-wave travel-times, influenced by the finite source, to derive the fault parameters of the earthquake. The results demonstrated that the earthquake exhibited a slow rupture propagation with a rupture velocity of approximately 1.5–2.0 km/s and asymmetric bilateral faulting. The two rupture directions were N60°E and N127°E, with rupture lengths of approximately 276 km and 231 km, respectively. The rupture toward N60°E had a source duration of approximately 183 s, longer than that toward N127°E (approximately 156 s). Overall, the entire source duration of the earthquake faulting lasted approximately 183 s. Regarding historical seismicity in eastern Japan, the 2011 Tohoku earthquake not only ruptured a locked area in which large earthquakes have rarely occurred, but also ruptured the source regions of several historical earthquakes. With the exception of its slow rupture velocity and generation of a tsunami, the rupture features of the 2011 Tohoku earthquake were inconsistent with those of typical tsunami earthquakes.     

Geostatistical integration using 2-D electrical resistivity and 3-D gravity methods for detecting cavities in a karst area

Two-dimensional (2-D) electrical resistivity and three-dimensional (3-D) gravity explorations were undertaken to estimate the 3-D distribution of karst cavities at the area of Yongweol-ri in Korea where ground subsidence has occurred. Although the gravity method is a low-cost way of analyzing a 3-D continuous structure, its vertical resolving power is poor. In contrast, the electrical resistivity method can provide a 2-D and/or 3-D subsurface structure with a much higher lateral and vertical resolution than the gravity method. Accordingly, geostatistical methods and density information were used to enhance the 2-D resistivity structure revealed by the electrical resistivity method into a 3-D structure. The assumptions are, first, that each data set senses the same underlying geological structure in terms of different material properties and, secondly, that two different material properties are correlated locally or globally throughout the entire target area. As a result, the distribution of limestone cavities can be estimated under the assumption that they are mostly filled with groundwater and clayey soils and have abnormally low levels of resistivity and density. The estimated distribution corresponds with the grouting, borehole imaging and monitoring data. In this example, it can be seen that the integration analysis of 2-D electrical resistivity and 3-D gravity methods is a very powerful tool for 3-D subsurface imaging and that the method can provide enhanced imaging capabilities for 3-D cavities.     

Impact of the 11 March 2011, Great East Japan earthquake and tsunami on the chemical industry

The Great East Japan earthquake and tsunami damaged or destroyed many industrial facilities housing or processing hazardous substances, such as refineries, petrochemical facilities and other types of chemical industry. This showed that also generally well prepared countries are at risk of suffering natural hazard triggered technological (Natech) accidents. An analysis of data collected from open sources and through interviews with authorities was performed to understand the main reasons for the industrial damage and downtime as well as the extent of hazardous-materials releases and the associated impact on society. The analysis of the data set confirmed the findings from other studies with respect to main damage and failure modes, as well as hazardous-materials release paths. In addition, gaps in Natech risk management were identified. Based on the data analysis and interviews lessons learned in support of a more far-reaching Natech risk management are presented.     

基于线性回归的重力数据三维RMS成像

这里提出了一种基于线性回归的RMS成像方法,通过该方法可以分析判断重力异常源所处位置及赋存范围的方法。经过线性回归,得相关系数、斜率和截距三个参数,利用这三个参数计算RMS并对场源进行成像分析,以得到场源的位置及赋存状况。为了验证算法的可行性,通过简单模型和实测重力异常数据进行试验分析,验证了该方法可以有效地显示出异常地质体的位置及空间赋存状态。     

前陆盆地挠曲沉降和沉积过程3D模型研究

前陆盆地是在造山带负荷作用下岩石圈发生挠曲沉降而形成的,并且被主要从造山带搬运的沉积物所充填。为了更好地理解和认识前陆盆地的形成演化机制,特别是受控于周缘多个造山带活动所形成的前陆盆地的演化机制,本文通过建立前陆盆地挠曲沉降与沉积过程的3-D模型,模拟展示了造山带逆冲推覆作用、岩石圈挠曲沉降响应及在山盆体系中由于动力地形变化而导致的河流体系的发育变化及其产生的剥蚀和沉积过程。模型的建立和实验完整体现了逆冲推覆、弹性挠曲沉降和沉积物搬运这三者之间的耦合机制,为全面深入研究前陆盆地动力学提供了理论依据和方法。     

Regional differences in cognitive dissonance in evacuation behavior at the time of the 2011 Japan earthquake and tsunami

This paper constructs an evacuation decision-making model that takes cognitive dissonance into consideration. The purpose of this construction is to clarify the psychological mechanism for the evacuation behavior of residents during an emergency, based on Akerlof and Dickens (Akerlof and Dickens Am Econ Rev 72:307–319, 1982) "The economic consequences of cognitive dissonance". Specifically, we empirically explore people’s psychological mechanism (e.g., cognitive dissonance) for evacuation behavior when a tsunami disaster occurs. As a result, we show that the level of anxiety depends on the area where residents live and that the average anxiety of residents is mostly correlated with the level of damage of past disasters, and that it is affected also by the ages of residents. Since the level of anxiety largely affects an individual’s evacuation behavior, this result can indicate for what kinds of people intervention and assistance are required based on the level of anxiety. A high level of anxiety basically promotes evacuation. Since our results show that anxiety is increased by the experience of tsunamis, education having people virtually experience tsunamis may increase evacuation rates efficiently.

     

Modeling of strong motion generation areas of the 2011 Tohoku,Japan earthquake using modified semi-empirical technique

Modification in the semi-empirical technique for the simulation of strong ground motion has been introduced to incorporate the strong motion generation areas (SMGA) in the modeled rupture plane. Strong motion generation areas identified within the rupture plane of the Tohoku earthquake of March 11, 2011 (M _w = 9.0), have been modeled using this modified technique. Two different source models having four and five SMGAs, respectively, are considered for modeling purpose. Strong motion records using modified semi-empirical technique have been simulated at two near-field stations located at epicentral distance of 137 and 140 km, respectively, using two different source models. Comparison of the observed and simulated acceleration waveforms is made in terms of root mean square error (RMSE) at both stations. Minimum root mean square error of the waveform comparison has been obtained at both the stations for source model having five SMGAs. Simulations from same rupture model have been made at other four stations lying at epicentral distance between 154 and 249 km. Comparison of observed and simulated records has been made in terms of RMSE in acceleration records, velocity records and response spectra at each six station. Simulations have been made at six other stations to obtain distribution of peak ground acceleration and peak ground velocity with hypocentral distance. Peak ground acceleration and velocity from simulated and observed records are compared at twelve stations surrounding the source of Tohoku earthquake. Comparison of waveforms and parameters extracted from observed and simulated strong motion records confirms the efficacy of the developed modified technique to model earthquake characterized by SMGAs.     

非结构网格上的三维浅水流动数值模型

针对当前复杂环境水流模拟的需求,建立了新型的基于特征型高分辨率数值算法的三维非结构网格浅水动力模型。模型采用有限体积法离散sigma坐标下的三维浅水方程,运用Roe黎曼近似解评估水平界面通量。模型网格拟合边界能力强,可根据需要局部加密;格式数值性能优良,具有守恒性、单调迎风性、高数值分辨率等特性。同时,应用干湿判别法处理动边界,以适应浅滩地形漫/露过程模拟的需要。封闭水池内部风生环流、干河床上溃坝过程和长江口实际潮流场的模拟从不同侧面展示了模型的特点,结果表明它能够准确地预测水流的三维流动结构,而且计算简单高效,具有良好的数值稳定性。     

3-D seismic structure of the Kachchh,Gujarat, and its implications for the earthquake hazard mitigation

Several pieces of studies on the January 26, 2001, Bhuj earthquake (Mw 7.6) revealed that the mainshock was triggered on the hidden unmapped fault in the western part of Indian stable continental region that caused a huge loss in the entire Kachchh rift basin of Gujarat, India. Occurrences of infrequent earthquakes of Mw 7.6 due to existence of hidden and unmapped faults on the surface have become one of the key issues for geoscientific research, which need to be addressed for evolving plausible earthquake hazard mitigation model. In this study, we have carried out a detailed autopsy of the 2001 Bhuj earthquake source zone by applying three-dimensional (3-D) local earthquake tomography (LET) method to a completely new data set consisting of 576 local earthquakes recorded between November 2006 and April 2009 by a seismic network consisting of 22 numbers of three-component broadband digital seismograph stations. In the present study, a total of 7560 arrival times of P-wave (3820) and S-wave (3740) recorded at least 4 seismograph stations were inverted to assimilate 3-D P-wave velocity (Vp), S-wave velocity (Vs), and Poisson’s ratio (σ) structures beneath the 2001 Bhuj earthquake source zone for reliable interpretation of the imaged anomalies and its bearing on earthquake hazard of the region. The source zone is located near the triple junction formed by juxtapositions of three Indian, Arabian, and Iranian tectonic plates that might have facilitated the process of brittle failure at a depth of 25 km beneath the KRB, Gujarat, which caused a gigantic loss to both property and persons of the region. There may be several hidden seismogenic faults around the epicentral zone of the 2001 Bhuj earthquake in the area, which are detectable using 3-D tomography to minimize earthquake hazard for a region. We infer that the use of detailed 3-D seismic tomography may offer potential information on hidden and unmapped faults beneath the plate interior to unravel the genesis of such big damaging earthquakes. This study may help in evolving a comprehensive earthquake risk mitigation model for regions of analogous geotectonic settings, elsewhere in the world.     

三维梁-颗粒模型在岩石材料破坏模拟中的应用

用三维梁-颗粒模型BPM3D(beam-particlemodelinthreedimensions)对岩石类非均质脆性材料的力学性质和破坏过程进行了数值模拟。梁-颗粒模型是在离散单元法基础上,结合有限单元法中的网格模型提出的用于模拟岩石类材料损伤破坏过程的数值模型。在模型中,材料在细观层次上被离散为颗粒单元集合体,相邻颗粒单元由有限单元法中的弹脆性梁单元联结。梁单元的力学性质均按韦伯(Weibull)分布随机赋值,以模拟岩石类材料力学参数的空间变异性。材料内部裂纹通过断开梁单元来模拟。通过自动生成的非均质材料模型对岩石类材料的破坏机理进行研究。岩石类非均质脆性材料在单轴压缩状态下破坏过程细观数值模拟结果显示,岩石材料宏观破坏是由于其内部细观裂纹产生、扩展、贯通的结果。通过数值模拟结果之间的对比分析,揭示出岩石试样宏观破坏模式随细观层次上韦伯分布参数的变化而不同。与实际矿柱破坏形态的对比分析表明了模型的适用性。根据数值模拟结果对岩石类非均质材料的破坏机理进行了探讨。     

重力异常二维正演中的无网格方法

无网格法是一类新型数值算法,具有精度高、高阶形函数构造与物性加载便利等特点,在计算力学领域应用广泛。将无网格方法（PIM、RPIM及EFGM）用于重力异常场二维正演计算：首先从重力异常二维变分问题出发,利用Galerkin法结合高斯积分公式推导了对应的无网格离散系统矩阵表达式;其次通过数值试验得出了RPIM-MQ、RPIM-exp及EFGM-exp形状参数的建议值,最后比较分析了最优形状参数下不同无网格法的计算效果。结果表明：无网格法适用于介质物性分布变化较大的重力异常二维正演,exp函数形状参数α_c最优取值区间为[1.5,1.7],β建议值为0.6,MQ函数q取值区间为-4.1~1.9;EFGM较PIM及RPIM具有更高的计算精度。     

长江口斜压诊断模式三维流场数值模拟

建立了σ坐标系下长江口斜压诊断模式三维流场数学模型,采用k-kl二方程紊流闭合模型求解垂向涡粘系数,计算域内恒定、非均匀盐度场反映了计算域密度斜压效应。验证结果表明,斜压诊断模式的模拟成果较切合实际。模型成功复演了长江口水域由往复流向旋转流过渡的流场特性;潮波从四条汊道传入后在分汊口附近相遇形成多个会潮点,其位置随着径流的大小、潮汐的强弱、汊道的形态等因素而变化。     

Using a genetic algorithm for 3-D inversion of gravity data in Fuerteventura (Canary Islands)

The use of genetic algorithms in geophysical inverse problems is a relatively recent development and offers many advantages in dealing with the non-linearity inherent in such applications. We have implemented a genetic algorithm to efficiently invert a set of gravity data. Employing several fixed density contrasts, this algorithm determines the geometry of the sources of the anomaly gravity field in a 3-D context. The genetic algorithms, based on Darwins theory of evolution, seek the optimum solution from an initial population of models, working with a set of parameters by means of modifications in successive iterations or generations. This searching method traditionally consists of three operators (selection, crossover and mutation) acting on each generation, but we have added a further one, which smoothes the obtained models. In this way, we have designed an efficient inversion gravity method, confirmed by both a synthetic example and a real data set from the island of Fuerteventura. In the latter case, we identify crustal structures related to the origin and evolution of the island. The results show a clear correlation between the sources of gravity field in the model and the three volcanic complexes recognized in Fuerteventura by other geological studies.