WEB-IS (integrated system): an overall view

WEB-IS, Web-based Integrated System, allows remote, interactive visualization of large-scale 3-D data over the Internet, along with data analysis and data mining. In this paper, we discuss the overall structure of WEB-IS. Up until now we have developed three sub-modules geared towards geophysical problems. WEB-IS1 allows geoscientists to navigate through their 3-D geophysical data, such as seismic structures or numerical simulations, and interactively analyze the statistics or apply data-mining techniques, such as cluster analysis. WEB-IS2 lets a user control Amira (a powerful 3-D visualization package) remotely and analyze, render and view large datasets across the Internet. WEB-IS3 is an imaging service that enables the user to control the scale of features to view through interactive zooming. In the near future, we propose to integrate the three components together through a middleware framework called NaradaBrokering (iNtegrated Asynchronous Real-time Adaptive Distributed Architecture, a distributed messaging infrastructure that can be used to intelligently route data between the originators and registered consumers) without regard for time or location. As a result, WEB-IS will improve its scalability and acquire properties of fault-tolerance. WEB-IS uses a combination of Java, C++, and through the use of NaradaBrokering will seamlessly integrate the server-side processing and user interaction utilities on the client. The server takes care of the processor intensive tasks, such as visualization and data mining, and sends either the resulting bitmap image or statistical results to the middleware across the Internet for viewing. WEB-IS is an easy-to-use service, which will eventually help geoscientists collaborate from different sites in a natural manner. It will be very useful in the next 10 years because of the increasing number of space missions and geophysical campaigns.     

Remote data analysis, visualization and problem solving environment (PSE) based on wavelets in the geosciences

We discuss here the issues faced by earth scientists in analyzing and visualizing large datasets over a GRID-like setup from a client-server perspective. We approach this problem by using a remote, web-based visualization and data analysis framework, called WEB-IS, and by employing second-generation wavelets as a means for reducing the amount of data transferred and for extracting coherent features in complex geophysical flows and surface faulting patterns. As an example, we describe how onboard processors on satellites can function as a server for beaming down extracted information to the client computer on the ground, thus exemplifying WEB-IS as a viable middleware on a GRID network for geosciences.Reviewed by: Prof. I.D. Clemens, Prof. L. Hanyk     

Multi-resolution clustering analysis and 3-D visualization of multitudinous synthetic earthquakes

We use modern and novel techniques to study the problems associated with detection and analysis of multitudinous seismic events, which form the background for isolated great earthquakes. This new approach involves multivariate analysis of low and large magnitude events recorded in space over a couple of centuries in time. We propose here the deployment of the clustering scheme both for extracting small local structures and large-scale trends in synthetic data obtained from four numerically simulated models with: uniform properties (U), a Parkfield-type asperity (A), fractal brittle properties (F), and multi-size-heterogeneity fault zone (M). The mutual nearest neighbor (mnn) clustering scheme allows for extraction of multi-resolutional seismic anomalies in both the spatio-temporal and multi-dimensional feature space. We demonstrate that the large earthquakes are correlated with a certain pathway of smaller events. Visualization of the anomalies by using a recently introduced visualization package Amira reveals clearly the spatio-temporal relationships between clusters of small, medium and large earthquakes, indicating significant stress relaxation across the entire fault region. We demonstrate that this mnn scheme can extract distinct clusters of the smallest events, which precede and follow a singularly large magnitude earthquake. These clusters form larger spatio-temporal structures comprising a series of large earthquakes. The link between the large and medium magnitude events is not so clearly understood. Short-ranged correlations are dominated by strong spatio-temporal anomalies, thus reflecting the global seismic properties of the entire fault zone.Electronic Supplementary Material Supplementary material is available for this article if you access the article at . A link in the frame on the left on that page takes you directly to the supplementary material.     

Visualization of multi-scale dynamics of hydrous cold plumes at subduction zones

Recent increases in the computational power of high-performance computing systems have led to a large gap between the high-resolution runs of numerical simulations—typically approaching 50–100 million tracers and 1–5 million grid points in two dimensions—and the modest resolution of 1–2 million pixels for conventional display devices. This technical problem is further compounded by the variety of fields produced by numerical simulations and the limited bandwidth available through the Internet in the course of collaborative ventures. We have developed a visualization system using the paradigm of web-based inquiry to address these mounting problems. We have employed, as a case study, a problem involving two-dimensional multi-scale dynamics of hydrous cold plumes at subduction zones. A Lagrangian marker method, in which the number of markers varies dynamically, is used to delineate the many different fields, such as temperature, viscosity, strain, and chemical composition. We found commercially available software to be insufficient for our visualization needs and so we were driven to develop a new set of tools tailored to high-resolution, multi-aspect, multi-scale simulations, and adaptable to many other applications in which large datasets involving tens of millions of tracers with many different fields are prevalent. In order to address this gap in visualization techniques, we have developed solutions for remote visualization and for local visualization. Our remote visualization solution is a web-based, zoomable image service (WEB-IS) that requires minimal bandwidth while allowing the user to explore our data through time, across many thermo–physical properties, and through different spatial scales. For local visualization, we found it optimal to use bandwidth-intensive, high-resolution display walls for performing parallel visualization in order to best comprehend the causal and temporal relationships between the multiple physical and chemical properties in a simulation.Electronic Supplementary Material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s10069-004-0017-2     

Historical seismicity on the southern margin of the Siberian craton: new data

We discuss historical evidence for seismicity on the southern margin of the Siberian craton collected from old local newspapers. The reported earthquakes vary in magnitude from M = 2.5 to 4.5, and their macroseismic locations agree well with the regional tectonic framework. The new data prove seismic activity in the area and can be used in seismic risk assessment.     

Precursory activity in seismic belts containing the focal region of forthcoming large intraplate earthquake

Seismic activity in the region surrounding the foci is investigated for three severe earthquakes (two with a magnitude of 6.1 and one with a magnitude of 5.3) which have occurred in Japan in recent years. The most conspicuous feature commonly noticed is precursory activation of seismic belts which include the focal regions of main shocks. The repetition of the same pattern in the space-time distribution of earthquake occurrence along the seismic belt is also observed for each case. The precursory activity of seismic belts terminates in rather a short period and, after that, the area around the focus of the forthcoming large earthquake becomes quiescent, which demonstrates the appearance of the seismic gap of the second kind (Mogi, 1979). The periods of seismic quiescence for the cases investigated in this paper are longer than those which are given by the regression relationship between earthquake magnitude and precursor time proposed for example, by Sekiya (1977). However, our definition of anomalous seismic activity is clear, and it is possible to give a physical meaning to it as an increase in the local stress field in the seismic belt. We propose that a kind of coupling between intraplate tectonic blocks, analogous to interplate coupling in the subduction region, is responsible for the formation of the stress field relevant to these earthquakes. Although this is at present only one of the possible viewpoints on the formation of the focal region of large intraplate earthquakes, it may be worthwhile to study various precursory phenomena in-connection with this hypothesis.     

Mapping seismic risk: the current crisis

The seismic risks to which populations are exposed should be estimated reliably for mitigation and preparation of response to disastrous earthquakes. Three parameters need to be known: Population numbers, properties of the built environment, and the seismic hazard. If we focus on large cities, we can say that at least one of these is known satisfactorily, namely the population, but not the other two. In the developing world, the numbers of buildings in a city are known only approximately, their distribution into building types (resistance to shaking) has to be assumed, and the distribution of types throughout the city is unknown. Recent verification of the world seismic hazard map has shown that it is grossly misleading: Instrumental measurements of accelerations due to six earthquakes were about three times larger, on average, than the maximum likely accelerations shown on the map; the macroseismic intensities reported for the last 60 earthquakes with M ≥ 7.5 were all significantly larger than expected, based on the hazard map (by 2.3 intensity units for the 12 deadliest earthquakes); and calculations of losses of life based on the hazard map underestimate the losses sustained in the 12 recent earthquakes with more than 1,000 fatalities by two to three orders of magnitude. This means that the seismic risk in most of the approximately 1,000 large cities at risk in the developing world is unknown. To remedy this intolerable situation, models for the built environment in cities need to be constructed, using cost-effective analyses of satellite images, and worst case scenario estimates of the losses in case of the nearest maximum credible earthquake.     

Inelastic displacement demand of RC buildings subjected to earthquakes generated by intermediate-depth Vrancea seismic source

Evaluating inelastic displacement demand of structures exposed to seismic hazard is required for the design of new buildings as well as for seismic risk assessment of existing structures. Most of the buildings are designed to withstand strong earthquakes by responding in the nonlinear range. Having special parts of the structure designed to develop a stable hysteretic behaviour allows the structure to deform in order to accommodate the displacement demand imposed by strong ground motions. This paper is centred on finding a correspondence between the maximum elastic and inelastic displacement responses of the single degree of freedom (SDOF) systems subjected to earthquakes generated by Vrancea seismic source. Vrancea intermediate-depth earthquakes are responsible for the seismic hazard throughout Romanian territory. They have distinctive features, such as large displacement demand and large predominant periods, which makes Romania a special seismic environment. Using a database of Romanian and Japanese strong ground motions generated by intermediate-depth earthquakes and performing nonlinear dynamic analysis on the SDOF oscillators following the Takeda model, this study estimates the inelastic to elastic displacement ratio of reinforced concrete systems. Soil conditions, epicentral distance and magnitude influence on inelastic response is analysed using constant ductility response spectra. The main findings of the study are: the local increase of the inelastic to elastic displacement ratio for type C soil (Eurocode 8 classification) for large magnitude earthquakes and the significant effect of soil conditions on the inelastic response of the SDOF systems. The inelastic amplification was evaluated using a functional form depending on system ductility, soil conditions and earthquake magnitude.

     

Strong earthquakes associated with high amplitude daily geomagnetic variations

In this work, we attempt to quantify forces that result from the interaction between the induced Sq-variation currents in the Earth’s lithosphere and the regional Earth’s magnetic field, in order to assess its influence on the tectonic stress field and on seismic activity. The study area is the Sinai Peninsula, a seismically active region where both seismic and magnetic data are available. We show that both short-term and long-term magnetic changes correlate with the seismic activity extending to this area in other previous studies. We also analyze a set of large earthquakes and magnetic data from observatories around the world to deduce a relationship between earthquake magnitude and maximum distance up to which precursory variations of the magnetic field are observed.     

Toward petascale earthquake simulations

Earthquakes are among the most complex terrestrial phenomena, and modeling of earthquake dynamics is one of the most challenging computational problems in science. Computational capabilities have advanced to a state where we can perform wavefield simulations for realistic three-dimensional earth models, and gain more insights into the earthquakes that threaten California and many areas of the world. The Southern California Earthquake Center initiated a major earthquake research program called TeraShake to perform physics-based numerical simulations of earthquake processes for large geographical regions, at high resolution, and for high frequencies. For a large scale simulation such as TeraShake, optimization problems tend to emerge that are not significant in smaller scale simulations. This involves both large parallel computation and also massive data management and visualization coordination. In this paper, we describe how we performed single-processor optimization of the TeraShake AWM application, optimization of the I/O handling, and optimization of initialization. We also look at the challenges presented by run-time data archive management and visualization. The improvements made to the TeraShake AWM code enabled execution on the 40k IBM Blue Gene processors and have created a community code that can be used by seismologists to perform petascale earthquake simulations.     

Web-based service of a visualization package “Amira” for the geosciences

Amira is a powerful three-dimensional visualization package that has been employed recently by the science and engineering communities to gain insight into their data. We discuss a new paradigm for the use of Amira in the Earth sciences that relies on the client-server paradigm. We have developed a module called WEB-IS2, which provides web-based access to Amira. This tool allows Earth scientists to manipulate Amira controls remotely and to analyze, render and view large datasets through the Internet without regard for time or location. This could have important ramifications for GRID computing.     

Probabilistic seismic hazard analysis for Sumatra, Indonesia and across the Southern Malaysian Peninsula

The ground motion hazard for Sumatra and the Malaysian peninsula is calculated in a probabilistic framework, using procedures developed for the US National Seismic Hazard Maps. We constructed regional earthquake source models and used standard published and modified attenuation equations to calculate peak ground acceleration at 2% and 10% probability of exceedance in 50 years for rock site conditions. We developed or modified earthquake catalogs and declustered these catalogs to include only independent earthquakes. The resulting catalogs were used to define four source zones that characterize earthquakes in four tectonic environments: subduction zone interface earthquakes, subduction zone deep intraslab earthquakes, strike-slip transform earthquakes, and intraplate earthquakes. The recurrence rates and sizes of historical earthquakes on known faults and across zones were also determined from this modified catalog. In addition to the source zones, our seismic source model considers two major faults that are known historically to generate large earthquakes: the Sumatran subduction zone and the Sumatran transform fault. Several published studies were used to describe earthquakes along these faults during historical and pre-historical time, as well as to identify segmentation models of faults. Peak horizontal ground accelerations were calculated using ground motion prediction relations that were developed from seismic data obtained from the crustal interplate environment, crustal intraplate environment, along the subduction zone interface, and from deep intraslab earthquakes. Most of these relations, however, have not been developed for large distances that are needed for calculating the hazard across the Malaysian peninsula, and none were developed for earthquake ground motions generated in an interplate tectonic environment that are propagated into an intraplate tectonic environment. For the interplate and intraplate crustal earthquakes, we have applied ground-motion prediction relations that are consistent with California (interplate) and India (intraplate) strong motion data that we collected for distances beyond 200 km. For the subduction zone equations, we recognized that the published relationships at large distances were not consistent with global earthquake data that we collected and modified the relations to be compatible with the global subduction zone ground motions. In this analysis, we have used alternative source and attenuation models and weighted them to account for our uncertainty in which model is most appropriate for Sumatra or for the Malaysian peninsula. The resulting peak horizontal ground accelerations for 2% probability of exceedance in 50 years range from over 100% g to about 10% g across Sumatra and generally less than 20% g across most of the Malaysian peninsula. The ground motions at 10% probability of exceedance in 50 years are typically about 60% of the ground motions derived for a hazard level at 2% probability of exceedance in 50 years. The largest contributors to hazard are from the Sumatran faults.     

Implementation of Sugeno: ANFIS for forecasting the seismic moment of large earthquakes over Indo-Himalayan region

The earthquake is known to be an unpredictable geophysical phenomenon. Only few seismic indicators and assumptions of earthquakes can be predicted with probable certainty. This study attempts to analyze the earthquakes over the Indo-Himalayan Border region including Bhutan, Bangladesh, Nepal, China and India during the period from 1995 to 2015. Bangladesh, Bhutan and China borders experience fewer earthquakes than Nepal and India border regions. However, Indo-China rim has inconsistency and vast range in its magnitude. Bangladesh though is a small country with respect to others, but it experiences earthquakes comparable to Bhutan. Nepal experiences highest number of earthquakes. In the last 20 years around 800 records have been observed with moment magnitude > 4.0 Richter scale, while very few records (around 10–12) have been observed for large earthquakes having moment magnitude > 6.0 Richter scale over the region. In this study adaptive neuro-fuzzy inference system has been implemented to assess the predictability of seismic moment associated with large earthquakes having the moment magnitude between 6.0 and 8.0 Richter scales using different combination of epochs, technique and membership functions. The Gaussian membership function with hybrid technique and 40 epochs is observed to be the reasonable model on the basis of the selected spatial and temporal scale. The forecast error in terms of root-mean-square error with the stopping criterion 0.001 has been observed to be 0.006 in case of large earthquakes (> 6.5 Richter scale), that is, forecast accuracy of 99.4%. The model bias of 0.6% may be due to inadequate number of large earthquakes having moment magnitude > 6.5 Richter scale over the region.     

The Earthquake Threat in Southwestern British Columbia: A Geologic Perspective

Ten moderate to large (magnitude 6–7) earthquakes have occurred in southwestern British Columbia and northwestern Washington in the last 130 years. A future large earthquake close to Vancouver, Victoria, or Seattle would cause tens of billions of dollars damage and would seriously impact the economies of Canada and the United States. An improved understanding of seismic hazards and risk in the region has been gained in recent years by using geologic data to extend the short period of instrumented seismicity. Geologic studies have demonstrated that historically unprecedented, magnitude 8 to 9 earthquakes have struck the coastal Pacific Northwest on average once every 500 years over the last several thousand years; another earthquake of this size can be expected in the future. Geologic data also provide insights into the likely damaging effects of future large earthquakes in the region. Much of the earthquake damage will result directly from ground shaking, but damage can also be expected from secondary phenomena, including liquefaction, landslides, and tsunamis. Vancouver is at great risk from earthquakes because important infrastructure, including energy and transportation lifelines, probably would be damaged or destroyed by landslides and liquefaction-induced ground failure.     

Three-dimensional earthquake locations and upper crustal structure of the Sannio-Matese region (southern Italy)

The Sannio-Matese region is one of the most seismically active regions of Italy and has been struck by large historical earthquakes. At present, the area is characterized by low magnitude background seismicity and small seismic sequences following M4 main events. In this paper, we show Vp and Vp/Vs models and 3D locations for a complete set of earthquakes occurring in the period 1991–2001. We observe a significant crustal heterogeneity, with large scale east-verging high Vp fault-related-folds, stacked by the Pliocene compression. The relocated earthquakes cluster along a 70° east-dipping, NW-striking plane located at the border of the high Vp thrust units. Normal fault earthquakes related to the young and active extension occur within these high Vp zones, interpreted as high strength material. We expect large future earthquakes to occur within these high Vp zones actually characterized by low magnitude seismicity at their borders.     

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.     

Earthquake imprints on a lacustrine deltaic system: The Kürk Delta along the East Anatolian Fault (Turkey)

Deltas contain sedimentary records that are not only indicative of water‐level changes, but also particularly sensitive to earthquake shaking typically resulting in soft‐sediment‐deformation structures. The Kürk lacustrine delta lies at the south‐western extremity of Lake Hazar in eastern Turkey and is adjacent to the seismogenic East Anatolian Fault, which has generated earthquakes of magnitude 7. This study re‐evaluates water‐level changes and earthquake shaking that have affected the Kürk Delta, combining geophysical data (seismic‐reflection profiles and side‐scan sonar), remote sensing images, historical data, onland outcrops and offshore coring. The history of water‐level changes provides a temporal framework for the depositional record. In addition to the common soft‐sediment deformation documented previously, onland outcrops reveal a record of deformation (fracturing, tilt and clastic dykes) linked to large earthquake‐induced liquefactions and lateral spreading. The recurrent liquefaction structures can be used to obtain a palaeoseismological record. Five event horizons were identified that could be linked to historical earthquakes occurring in the last 1000 years along the East Anatolian Fault. Sedimentary cores sampling the most recent subaqueous sedimentation revealed the occurrence of another type of earthquake indicator. Based on radionuclide dating (¹³⁷Cs and ²¹⁰Pb), two major sedimentary events were attributed to the ad 1874 to 1875 East Anatolian Fault earthquake sequence. Their sedimentological characteristics were determined by X‐ray imagery, X‐ray diffraction, loss‐on‐ignition, grain‐size distribution and geophysical measurements. The events are interpreted to be hyperpycnal deposits linked to post‐seismic sediment reworking of earthquake‐triggered landslides.     

Generation of isoseismal maps for Mexico using velocity and acceleration data

Acceleration and velocity data of large earthquakes recorded by the Mexican National Broadband Seismological Network were used to generate isoseismal maps for Mexico. The seismic data consist of 99 events recorded at 27 seismic stations located in the southern and northern regions of Mexico from 2004 to 2009. The magnitude (Mw) of these events ranged from 4.1 to 7.1. Peak ground velocity values (PGV) and peak ground accelerations (PGA) were estimated, and velocity-derived peak ground accelerations (PGAv) were calculated. No important differences between PGAv and PGA values were found; thus, both parameters were used in the isoseismal determination. The generated synthetic isoseismal maps were compared with those from an existing catalog for large earthquakes in Mexico. Using empiric relations between PGA and MMI (modified Mercalli intensity) and PGV and MMI, the obtained simulated intensity maps showed similar behavior to those reported in the catalog. The results indicate that the PGAv values can be used to determine intensities when acceleration records (PGA) are unavailable.