Application of a modified pattern informatics method to forecasting the locations of future large earthquakes in the central Japan

We propose a modification of the Pattern Informatics (PI) method that has been developed for forecasting the locations of future large earthquakes. This forecast is based on analyzing the space–time patterns of past earthquakes to find possible locations where future large earthquakes are expected to occur. A characteristic of our modification is that the effect of errors in the locations of past earthquakes on the output forecast is reduced. We apply the modified and original methods to seismicity in the central part of Japan and compared the forecast performances. We also invoke the Relative Intensity (RI) of seismic activity and randomized catalogs to constitute null hypotheses. We do statistical tests using the Molchan and Relative Operating Characteristic (ROC) diagrams and the log-likelihoods and show that the forecast for using the modified PI method is generally better than the competing original-PI forecast and the forecasts from the null hypotheses. Using the bootstrap technique with Monte-Carlo simulations, we further confirm that earthquake sequences simulated based on the modified-PI forecast can be statistically the same as the real earthquake sequence so that the forecast is acceptable. The main and innovative science in this paper is the modification of the PI method and the demonstration of its applicability, showing a considerable promise as an intermediate-term earthquake forecasting tool.     

Stress transfer in earthquakes, hazard estimation and ensemble forecasting: Inferences from numerical simulations

Observations indicate that earthquake faults occur in topologically complex, multi-scale networks driven by plate tectonic forces. We present realistic numerical simulations, involving data-mining, pattern recognition, theoretical analyses and ensemble forecasting techniques, to understand how the observable space–time earthquake patterns are related to the fundamentally inaccessible and unobservable dynamics. Numerical simulations can also help us to understand how the different scales involved in earthquake physics interact and influence the resulting dynamics. Our simulations indicate that elastic interactions (stress transfer) combined with the nonlinearity in the frictional failure threshold law lead to the self-organization of the statistical dynamics, producing 1) statistical distributions for magnitudes and frequencies of earthquakes that have characteristics similar to those possessed by the Gutenberg–Richter magnitude–frequency distributions observed in nature; and 2) clear examples of stress transfer among fault activity described by stress shadows, in which an earthquake on one group of faults reduces the Coulomb failure stress on other faults, thereby delaying activity on those faults. In this paper, we describe the current state of modeling and simulation efforts for Virtual California, a model for all the major active strike slip faults in California. Noting that the Working Group on California Earthquake Probabilities (WGCEP) uses statistical distributions to produce earthquake forecast probabilities, we demonstrate that Virtual California provides a powerful tool for testing the applicability and reliability of the WGCEP statistical methods. Furthermore, we show how the simulations can be used to develop statistical earthquake forecasting techniques that are complementary to the methods used by the WGCEP, but improve upon those methods in a number of important ways. In doing so, we distinguish between the “official” forecasts of the WGCEP, and the “research-quality” forecasts that we discuss here. Finally, we provide a brief discussion of future problems and issues related to the development of ensemble earthquake hazard estimation and forecasting techniques.     

Reply to the comments of S. Karato on “Petrofabrics and seismic properties of garnet peridotites from the UHP Sulu terrane (China)” by Xu et al. [Tectonophysics 421 (2006) 111–127]

We thank Karato for his interest in our recent paper in Tectonophysics [Xu, Z.Q., Wang, Q., Ji, S.C.. Chen, J., Zeng, L.S., Yang, J.S., Chen, F.Y., Liang, F.H., Wenk, H.R., 2006. Petrofabrics and seismic properties of garnet peridotite from the UHP Sulu terrane (China): Implications for olivine deformation mechanism in a cold and dry subducting continental slab. Tectonophysics 421, 111–127]. However, Karato is incorrect to take water or water-related species (herein referred to generically as ‘water’) as an almost unique cause for any olivine fabric transition in nature no matter what geological conditions are. In this Reply, we will first give a succinct introduction to the context of the ongoing debate, and then provide some additional arguments in favour of our conclusion that “the olivine C-type fabric of the Zhimafang garnet-bearing peridotites was formed under the conditions of high pressure, low temperature and low water content during the subduction of cold and dry supracrustal rocks of the Yangtze plate”. The map of olivine slip systems as a function of differential stress and water content, drawn by Karato and coworkers, cannot be reliably extrapolated to Earth’s interior until relevant scaling laws have been established. The role of high water content on the fabric transitions is more likely to modify recrystallization kinetics than to switch the dominant slip systems.     

Current status of multimodel superensemble and operational NWP forecast of the Indian summer monsoon

In the last thirty years great strides have been made by large-scale operational numerical weather prediction models towards improving skills for the medium range time-scale of 7 days. This paper illustrates the use of these current forecasts towards the construction of a consensus multimodel forecast product called the superensemble. This procedure utilizes 120 of the recent-past forecasts from these models to arrive at the training phase statistics. These statistics are described by roughly 10⁷ weights. Use of these weights provides the possibility for real-time medium range forecasts with the superensemble. We show the recent status of this procedure towards real-time forecasts for the Asian summer monsoon. The member models of our suite include ECMWF, NCEP/EMC, JMA, NOGAPS (US Navy), BMRC, RPN (Canada) and an FSU global spectral forecast model. We show in this paper the skill scores for day 1 through day 6 of forecasts from standard variables such as winds, temperature, 500 hPa geopotential height, sea level pressure and precipitation. In all cases we noted that the superensemble carries a higher skill compared to each of the member models and their ensemble mean. The skill matrices we use include the RMS errors, the anomaly correlations and equitable threat scores. For many of these forecasts the improvements of skill for the superensemble over the best model was found to be quite substantial. This real-time product is being provided to many interested research groups. The FSU multimodel superensemble, in real-time, stands out for providing the least errors among all of the operational large scale models.     

Using remotely sensed data to modify wind forcing in operational storm surge forecasting

Storm surges are abnormal coastal sea level events caused by meteorological conditions such as tropical cyclones. They have the potential to cause widespread loss of life and financial damage and have done so on many occasions in the past. Accurate and timely forecasts are necessary to help mitigate the risks posed by these events. Operational forecasting models use discretisations of the governing equations for fluid flow to model the sea surface, which is then forced by surface stresses derived from a model wind and pressure fields. The wind fields are typically idealised and generated parametrically. In this study, wind field datasets derived from remotely sensed data are used to modify the model parametric wind forcing and investigate potential improvement to operational forecasting. We examine two methods for using analysis wind fields derived from remotely sensed observations of three hurricanes. Our first method simply replaces the parametric wind fields with its corresponding analysis wind field for a period of time. Our second method does this also but takes it further by attempting to use some of the information present in the analysis wind field to estimate future wind fields. We find that our methods do yield some forecast improvement, most notably for our second method where we get improvements of up to 0.29 m on average. Importantly, the spatial structure of the surge is changed in some places such that locations that were previously forecast small surges had their water levels increased. These results were validated by tide gauge data.     

订正与集成多模式的中国季度降水预测

针对两个最新换代的季度集合预测系统对中国季度降水预测中存在的系统缺陷,应用改进的贝叶斯联合概率模型（BJP）加以订正。对订正后的单一模式概率预测应用一种混合模型贝叶斯模型平均（BMA）方法加以集成,以综合各模式的优势来提高中国季度降水预测技巧。结果表明:BJP模型可有效地消除集合模式预测的系统偏差,同时大幅提高了概率预测的可靠性。经过订正的欧洲中尺度天气预报中心的 System4预测在许多季度在中国的很大区域范围内都显示出了一定的预测技巧;而澳洲气象局的POAMA2.4预测只在个别季度局部范围内具有技巧。使用BMA对订正后的单一模式预测进行集成可显著提高对中国季度降水预测的精度,相比单一模式预测,技巧得分为正值的网格百分率分别提高了13.3%和20.0%。     

地震预测:从芦山地震到大陆地震

自从1990年以来,通过对青藏高原的调查和研究,认识到下地壳流动同步形成盆地和造山带,并受控于相关洋盆地幔软流圈向大陆的顺层流动和底辟作用。下地壳不均匀流动通过韧脆性中地壳热能-应变能转换孕育地震,部分发震能量通过上地壳脆性断层释放。在地震孕育过程中通常会伴生跨年度干旱和异常降雨,构成热灾害链。近5年内青藏高原东部连续发生汶川、玉树、芦山大地震,形成于从亚东流经羊八井、安多、玉树并分支流向汶川和芦山—康定的下地壳"热河"的仰冲式和侧冲式撞击作用。从2008年9月以来连续发表5篇论文,根据地壳热构造和热灾害链的时空结构对芦山地震的三要素进行了长期和中期预测。2008年9月预测从2013年开始可能发生大地震,2012年9月将鲜水河—安宁河—小江异常热流构造带5年内将发生多个7级地震的首个大震锁定在芦山或西昌。芦山地震只释放了亚东—羊八井—安多—玉树—鲜水河—安宁河—小江"热河"剩余热能中的一小部分,在西昌—会理—昭通地区、道孚—康定地区、通海—石屏地区近5年内很可能发生4个7级左右的地震。此外,华北典型的热灾害链结构表明震情严峻,环渤海地区近3年内很可能发生大地震。从地震热流体撞击机理与地震异常之间的关联性出发,提出了动态立体监测及短临预测地震的思路和方法。     

论地震数值预报——关于我国地震预报研究发展战略的思考

随着中国工业化和社会现代化进程的加快 ,中国城市化进程必将进一步加速。如何减轻地震灾害的问题正变得日益严峻。尽管地震预测是一个举世公认的国际性科学难题 ,但在强化各种减轻地震灾害措施的同时 ,仍须大力推进地震预测研究。为此 ,需要打破长期徘徊在以地震前兆异常监测为基础的经验性预测局面 ,把注意力尽快转向研究以动力学为基础的数值预报。以GPS为代表的空间对地观测技术 ,岩石圈巨型高分辨率宽频带流动地震台阵观测技术以及数值模拟技术已经为地震数值预报研究提供了前所未有的技术基础。以地震数值预报为目标的GPS阵列地壳形变连续观测 ,高分辨率地壳上地幔结构探测 ,地壳动力学 ,地震孕育和破裂过程的理论、模拟试验和实际观测 ,数据同化和计算软件的开发应成为今后研究的重点。现在的问题是 ,需要积极借助数值天气预报的经验 ,强化多学科 ,多部门的组织协调 ,尽早在有条件的地区开展地震动力学数值预报的科学试验。地震数值预报研究必将极大地促进中国地震科学基础研究和地震预报的进一步发展。     

Assessing performances of pattern informatics method: a retrospective analysis for Iran and Italy

Pattern informatics (PI) algorithm, which was introduced at the beginning of past decade, uses instrumental earthquake catalogs to investigate the time-dependent rate of seismicity in the study area and, based on the information from past events, calculates the probabilities for the occurrence of future large earthquakes. The main measure in this method is the number of events above a specified magnitude threshold M _c that is counted over a gridded area. PI has been applied in several regions of the world and different variants of the method have been developed over the past decade. Hence, the problem of formally evaluating and comparing the performances of the different PI variants needs to be addressed from an operational perspective, in order to identify the preferred application scheme and as well as to improve the performances of the method. In this study, PI is applied for the first time to the retrospective analysis of the earthquake catalogs of Iran and Italy, so as to check whether this method could forecast the past large events in these two regions with different level of data completeness and complex seismotectonic setting. The original PI algorithm and one of its modified variants, as well as the relative intensity (RI) model, are used to check the stability and statistical significance of the obtained results. In order to assess and compare the obtained results, the performances of the different PI variants are analyzed considering different evaluation strategies, which turn out to provide significantly different scores even for the same algorithm variant. We show that a critical point in the assessment of the obtained results is related with the definition and quantification of the space uncertainty of the issued forecasts, that is, with the extent of the territory where large earthquakes are to be expected. Accordingly, we emphasize the need for an appropriate definition of the evaluation strategies, clearly and unambiguously indicating the area where a large earthquake has to be expected. The study shows that, with respect to application in Iran and Italy, the performances of PI algorithm (both original and modified variants) are highly dependent on the selected evaluation strategy and do not provide better information than the simple RI model, which does not account for temporal properties of seismicity evolution. The overall performances can be improved by introducing specific thresholds that discard the less active cells; however, being based on some posterior optimization, a rigorous prospective testing is required to assess the forecasting capability of the method. In this paper, we aim to set up the rules for such testing, including advance definition of the evaluation strategy.     

Evaluation of real-time flash flood forecasts for Haiti during the passage of Hurricane Tomas, November 4–6, 2010

The January 2010 earthquake that devastated Haiti left its population ever more vulnerable to rainfall-induced flash floods. A flash flood guidance system has been implemented to provide real-time information on the potential of small (~70 km²) basins for flash flooding throughout Haiti. This system has components for satellite rainfall ingest and adjustment on the basis of rain gauge information, dynamic soil water deficit estimation, ingest of operational mesoscale model quantitative precipitation forecasts, and estimation of the times of channel flow at bankfull. The result of the system integration is the estimation of the flash flood guidance (FFG) for a given basin and for a given duration. FFG is the amount of rain of a given duration over a small basin that causes minor flooding in the outlet of the basin. Amounts predicted or nowcasted that are higher than the FFG indicate basins with potential for flash flooding. In preparation for Hurricane Tomas’ landfall in early November 2010, the FFG system was used to generate 36-h forecasts of flash flood occurrence based on rainfall forecasts of the nested high-resolution North American Model of the National Centers for Environmental Prediction. Assessment of the forecast flood maps and forecast precipitation indicates the utility and value of the forecasts in understanding the spatial distribution of the expected flooding for mitigation and disaster management. It also highlights the need for explicit uncertainty characterization of forecast risk products due to large uncertainties in quantitative precipitation forecasts on hydrologic basin scales.     

Reply to comments by H. Mashima on ‘Evolution of the eastern margin of Korea: constraints on the opening of the East Sea (Japan Sea)’ by Kim et al. [Tectonophysics 436 (2007) 37–55]

We thank H. Mashima for his interest in our recent article in Tectonophysics [Kim, H.J., Lee, G.H., Jou, H.T., Cho, H.M., Yoo, H.S., Park, G.T., Kim, J.S., 2007, Evolution of the eastern margin of Korea: Constraints on the opening of the East Sea (Japan Sea). Tectonophysics 436, 37–55.] and welcome the opportunity to respond to his comments. In our article we suggested that the southern part of the East Sea (Japan Sea) opened principally in the southeast direction in response to the northwestward subduction of the Pacific Plate beneath the Japan Arc. In contrast, Mashima claims that the opening of the East Sea was achieved in the south–southeast direction. However, there are many crucial things in his comments that we find scientifically unconvincing and misleading. In this reply, we give a detailed response to his comments.     

Operational tropical cyclone intensity prediction—an empirical technique

One very specific operational requirement of the Tropical Cyclone (TC) Programme of the Regional Specialized Meteorological Centre, New Delhi is to provide 12-hourly forecasts valid up to 48 h (preferably 72 h) on the intensity of cyclones over the southern Indian Seas. In this paper, a simple empirical model for predicting the intensity of TCs occurring in the Bay of Bengal is proposed. The model parameter has been determined from a database assembled on 30 recent cyclones, and the model itself is based on the assumption that a TC intensifies exponentially. A method for correcting the forecast during subsequent observation hours (6- or 12-h intervals) is also presented. The results show that the forecast skill for forecasts of up to 48 h is reasonably good. The absolute mean errors are less than 12 knots for 48-h forecasts, with the forecast skill decreasing with time. With the incorporation of a correction procedure based on the latest observations, some improvement in the forecast skill can be obtained. The model is expected to be useful to operational forecasters.     

Short term earthquake forecasts at Koyna,India

Earthquake activity is monitored in real time at the Koyna reservoir in western India, beginning from August 2005 and successful short term forecasts have been made of M ∼ 4 earthquakes. The basis of these forecasts is the observation of nucleation that precedes such earthquakes. Here we report that a total of 29 earthquakes in the magnitude range of 3.5 to 5.1 occurred in the region during the period of August 2005 through May 2010. These earthquakes could broadly be put in three zones. Zone-A has been most active accounting for 18 earthquakes, while 5 earthquakes in Zone-B and 6 in Zone-C have occurred. Earthquakes in Zone-A are preceded by well defined nucleation, while it is not the case with zones B and C. This indicates the complexity of the earthquakes processes and the fact that even in a small seismically active area of only 20 km × 30 km earthquake forecast is difficult.     

Modification of the pattern informatics method for forecasting large earthquake events using complex eigenfactors

Recent studies have shown that real-valued principal component analysis can be applied to earthquake fault systems for forecasting and prediction. In addition, theoretical analysis indicates that earthquake stresses may obey a wave-like equation, having solutions with inverse frequencies for a given fault similar to those that characterize the time intervals between the largest events on the fault. It is therefore desirable to apply complex principal component analysis to develop earthquake forecast algorithms. In this paper we modify the Pattern Informatics method of earthquake forecasting to take advantage of the wave-like properties of seismic stresses and utilize the Hilbert transform to create complex eigenvectors out of measured time series. We show that Pattern Informatics analyses using complex eigenvectors create short-term forecast hot-spot maps that differ from hot-spot maps created using only real-valued data and suggest methods of analyzing the differences and calculating the information gain.     

A risky business: Mining, rent and the neoliberalization of “risk”

Natural resource investment in the mining sector is often mediated through conflicts over rent distribution between corporate capital and landowner states. Recent rounds of neoliberal policy promoted by the World Bank have highlighted the need for landowner states to offer incentives in order to attract “high risk” capital investment. In Sub-Saharan Africa, in particular, countries have been pushed to offer attractive fiscal terms to capital, thereby lowering the proportion traditionally called rent. This paper examines how the concept of “risk” has been mobilized to legitimate such skewed distributional arrangements. While certain conceptions of social and ecological “risk” have been prevalent in political and social theoretic discourse on mining, such focus elides the overwhelming contemporary power of our notion of “neoliberal risk” – or the financial/market risks – in actually setting the distributional terms of mineral investment. We illustrate our argument by examining the nexus of World Bank mining policy promotion and Tanzanian policy in the late 1990s meant to attract foreign direct investment in gold production. In closing, we suggest that just as “risk” is used to legitimate attractive fiscal terms for investment, recent events highlight how skewed distribution of benefits may set into motion risks that corporate capital had not bargained for.     

Reply to Comment on “Analyses of Seismic Deformation at the Kibyra Roman Stadium,Southwest Turkey”

In their comment, Elitez and Yaltırak (2013) criticize our paper (Karabacak et al. 2013 ) by arguing that it contains misinterpretations and unreliable data for a fault cutting the Roman Stadium in the ancient city of Kibyra. However their comments are not based on strong arguments to disprove our geological and archeological field observations or previous data. Here we present additional data supported by new figures which validate the faulting both in the bedrock and in the stadium floor. We reject their comments and uphold our original viewpoint on the faulting at the Stadium.     

Experimental superensemble forecasts of tropical cyclones over the Bay of Bengal

This study entails the implementation of an experimental real time forecast capability for tropical cyclones over the Bay of Bengal basin of North Indian Ocean. This work is being built on the experience gained from a number of recent studies using the concept of superensemble developed at the Florida State University (FSU). Real time hurricane forecasts are one of the major components of superensemble modeling at FSU. The superensemble approach of training followed by real time forecasts produces the best forecasts for tracks and intensity (up to 5 days) of Atlantic hurricanes and Pacific typhoons. Improvements in track forecasts of about 25–35% compared to current operational forecast models has been noted over the Atlantic Ocean basin. The intensity forecasts for hurricanes are only marginally better than the best models. In this paper, we address tropical cyclone forecasts over the Bay of Bengal for the years 1996–2000. The main result from this study is that the position and intensity errors for tropical cyclone forecasts over the Bay of Bengal from the multimodel superensemble are generally less than those of all of the participating models during 1- to 3-day forecasts. Some of the major tropical cyclones, such as the November 1996 Andhra Pradesh cyclone and October 1999 Orissa super cyclone were well handled by this superensemble approach. A conclusion from this study is that the proposed approach may be a viable way to construct improved forecasts of Bay of Bengal tropical cyclone positions and intensity.     

Low-latitude “dusty events” vs. high-latitude “icy Heinrich events”

It has been proposed that tropical events could have participated in the triggering of the classic, high-latitude, iceberg-discharge Heinrich events (HE). We explore low-latitude Heinrich events equivalents at high resolution, in a piston core recovered from the tropical north-western African margin. They are characterized by an increase of total dust, lacustrine diatoms and fibrous lacustrine clay minerals. Thus, low-latitude events clearly reflect severe aridity events that occurred over Africa at the Saharan latitudes, probably induced by southward shifts of the Inter Tropical Convergence Zone. At a first approximation, it seems that there is more likely synchronicity between the high-latitude Heinrich Events (HEs) and low-latitude events (LLE), rather than asynchronous behaviours.     

Deep structure under Yellowstone National Park U.S.A.: A continental “hot spot”

In order to understand the origin of long-lived loci of volcanism (sometimes called “hot spots”) and their possible role in global tectonic processes, it is essential to know their deep structure. Even though some work has been done on the crustal, upper-mantle, and deep-mantle structure under some of these “hot spots”, the picture is far from clear. In an attempt to study the structure under the Yellowstone National Park U.S.A., which is considered to be such a “hot spot”, we recorded teleseisms using 26 telemetered seismic stations and three groups of portable stations. The network was operated within a 150 km radius centered on the Yellowstone caldera, the major, Quaternary volcanic feature of the Yellowstone region. Teleseismic delays of about 1.5 sec are found inside the caldera, and the delays remain high over a 100 km wide area around the caldera. The spatial distribution and magnitude of the delays indicate the presence of a large body of low-velocity material with horizontal dimensions corresponding approximately to the caldera size (40 km × 80 km) near the surface and extending to a depth of 200–250 km under the caldera. Using ray-tracing and inversion techniques, it is estimated that the compressional velocity inside the anomalous body is lower than in the surrounding rock by about 15% in the upper crust and by 5% in the lower crust and upper mantle. It is postulated that the body is partly composed of molten rock with a high degree of partial melting at shallow depths and is responsible for the observed Yellowstone volcanism. The large size of the partially molten body, taken together with its location at the head of a 350 km zone of volcanic propagation along the axis of the Snake River Plain, indicates that the volcanism associated with Yellowstone has its origin below the lithosphere and is relatively stationary with respect to plate motion. Using our techniques, we are unable to detect any measurable velocity contrast in the mantle beneath the low-velocity body, and, hence, we are unable to determine whether the Yellowstone melting anomaly is associated with a deep heat source or with any deep phenomenon such as a convection plume, chemical plume, or gravitational anchor.     

Shear‐wave splitting and earthquake forecasting

Seismic shear‐wave splitting (SWS) monitors the low‐level deformation of fluid‐saturated microcracked rock. We report evidence of systematic SWS changes, recorded above small earthquakes, monitoring the accumulation of stress before earthquakes that allows the time and magnitude of impending large earthquakes to be stress‐forecast. The effects have been seen with hindsight before some 15 earthquakes ranging in magnitude from an M1.7 seismic swarm event in Iceland to the Ms7.7 Chi‐Chi Earthquake in Taiwan, including a successfully stress‐forecast of a M5.0 earthquake in SW Iceland. Characteristic increases in SWS time‐delays are observed before large earthquakes, which abruptly change to deceases shortly before the earthquake occurs. There is a linear relationship between magnitudes and logarithms of durations of both increases and decreases in SWS time‐delays before large impending earthquakes. However, suitably persistent swarms of small earthquakes are too scarce for routine stress‐forecasting. Reliable forecasting requires controlled‐source cross‐hole seismics between neighbouring boreholes in stress‐monitoring sites (SMS). It would be possible to stress‐forecast damaging earthquakes worldwide by a global network of SMS in real time.