Temporal behavior of the background seismicity rate in central Japan, 1998 to mid-2003

We studied the temporal behavior of the background shallow seismicity rate in 700 circular areas across inland Japan. To search for and test the significance of the possible rate changes in background seismicity, we developed an efficient computational method that applies the space–time ETAS model proposed by Ogata in 1998 to the areas. Also, we conducted Monte Carlo tests using a simulated catalog to validate the model we applied. Our first finding was that the activation anomalies were found so frequently that the constant background seismicity hypothesis may not be appropriate and/or the triggered event model with constraints on the parameters may not adequately describe the observed seismicity. However, quiescence occasionally occurs merely by chance. Another outcome of our study was that we could automatically find several anomalous background seismicity rate changes associated with the occurrence of large earthquakes. Very significant seismic activation was found before the M6.1 Mt. Iwate earthquake of 1998. Also, possible seismic quiescence was found in an area 150 km southwest of the focal region of the M7.3 Western Tottori earthquake of 2000. The seismicity rate in the area recovered after the mainshock.     

Crustal structure of the Eastern Alps along the TRANSALP profile from wide-angle seismic tomography

The objective of the TRANSALP project is an investigation of the Eastern Alps with regard to their deep structure and dynamic evolution. The core of the project is a 340-km-long seismic profile at 12°E between Munich and Venice. This paper deals with the P-wave velocity distribution as derived from active source travel time tomography. Our database consists of Vibroseis and explosion seismic travel times recorded at up to 100 seismological stations distributed in a 30-km-wide corridor along the profile. In order to derive a velocity and reflector model, we simultaneously inverted refractions and reflections using a derivative of a damped least squares approach for local earthquake tomography. 8000 travel time picks from dense Vibroseis recordings provide the basis for high resolution in the upper crust. Explosion seismic wide-angle reflection travel times constrain both deeper crustal velocities and structure of the crust–mantle boundary with low resolution. In the resulting model, the Adriatic crust shows significantly higher P-wave velocities than the European crust. The European Moho is dipping south at an angle of 7°. The Adriatic Moho dips north with a gentle inclination at shallower depths. This geometry suggests S-directed subduction. Azimuthal variations of the first-break velocities as well as observations of shear wave splitting reveal strong anisotropy in the Tauern Window. We explain this finding by foliations and laminations generated by lateral extrusion. Based on the P-wave model we also localized almost 100 local earthquakes recorded during the 2-month acquisition campaign in 1999. Seismicity patterns in the North seem related to the Inn valley shear zone, and to thrusting of Austroalpine units over European basement. The alignment of deep seismicity in the Trento-Vicenza region with the top of the Adriatic lower crust corroborates the suggestion of a deep thrust fault in the Southern Alps.     

Upper mantle structure of the Northern Eurasia from peaceful nuclear explosion data

Several long-range seismic profiles were carried out in Russia with Peaceful Nuclear Explosions (PNE). The data from 25 PNEs recorded along these profiles were used to compile a 3-D upper mantle velocity model for the central part of the Northern Eurasia. 2-D crust and upper mantle models were also constructed for all profiles using a common methodology for wavefield interpretation. Five basic boundaries were traced over the study area: N1 boundary (velocity level, V = 8.35 km/s; depth interval, D = 60–130 km), N2 (V = 8.4 km/s; D = 100–140 km), L (V = 8.5 km/s; D = 180–240 km) and H (V = 8.6 km/s; D = 300–330 km) and structural maps were compiled for each boundary. Together these boundaries describe a 3-D upper mantle model for northern Eurasia. A map characterised the velocity distribution in the uppermost mantle down to a depth of 60 km is also presented. Mostly horizontal inhomogeneity is observed in the uppermost mantle, and the velocities range from the average 8.0–8.1 km/s to 8.3–8.4 km/s in some blocks of the Siberian Craton. At a depth of 100–200 km, the local high velocity blocks disappear and only three large anomalies are observed: lower velocities in West Siberia and higher velocities in the East-European platform and in the central part of the Siberian Craton. In contrast, the depths to the H boundary are greater beneath the craton and lower beneath in the West Siberian Platform. A correlation between tectonics, geophysical fields and crustal structure is observed. In general, the old and cold cratons have higher velocities in the mantle than the young platforms with higher heat flows.Structural peculiarities of the upper mantle are difficult to describe in form of classical lithosphere–asthenosphere system. The asthenosphere cannot be traced from the seismic data; in contrary the lithosphere is suggested to be rheologically stratified. All the lithospheric boundaries are not simple discontinuities, they are heterogeneous (thin layering) zones which generate multiphase reflections. Many of them may be a result of fluids concentrated at some critical P–T conditions which produce rheologically weak zones. The most visible rheological variations are observed at depths of around 100 and 250 km.     

Seismic lamination and anisotropy of the Lower Continental Crust

Seismic lamination in the lower crust associated with marked anisotropy has been observed at various locations. Three of these locations were investigated by specially designed experiments in the near vertical and in the wide-angle range, that is the Urach and the Black Forrest area, both belonging to the Moldanubian, a collapsed Variscan terrane in southern Germany, and in the Donbas Basin, a rift inside the East European (Ukrainian) craton. In these three cases, a firm relationship between lower crust seismic lamination and anisotropy is found. There are more cases of lower-crustal lamination and anisotropy, e.g. from the Basin and Range province (western US) and from central Tibet, not revealed by seismic wide-angle measurements, but by teleseismic receiver function studies with a P–S conversion at the Moho. Other cases of lamination and anisotropy are from exhumed lower crustal rocks in Calabria (southern Italy), and Val Sesia and Val Strona (Ivrea area, Northern Italy). We demonstrate that rocks in the lower continental crust, apart from differing in composition, differ from the upper mantle both in terms of seismic lamination (observed in the near-vertical range) and in the type of anisotropy. Compared to upper mantle rocks exhibiting mainly orthorhombic symmetry, the symmetry of the rocks constituting the lower crust is either axial or orthorhombic and basically a result of preferred crystallographic orientation of major minerals (biotite, muscovite, hornblende). We argue that the generation of seismic lamination and anisotropy in the lower crust is a consequence of the same tectonic process, that is, ductile deformation in a warm and low-viscosity lower crust. This process takes place preferably in areas of extension. Heterogeneous rock units are formed that are generally felsic in composition, but that contain intercalations of mafic intrusions. The latter have acted as heat sources and provide the necessary seismic impedance contrasts. The observed seismic anisotropy is attributed to lattice preferred orientation (LPO) of major minerals, in particular of mica and hornblende, but also of olivine. A transversely isotropic symmetry system, such as expected for sub-horizontal layering, is found in only half of the field studies. Azimuthal anisotropy is encountered in the rest of the cases. This indicates differences in the horizontal components of tectonic strain, which finally give rise to differences in the evolution of the rock fabric.     

Seismic potential of Southern Italy

To improve estimates of the long-term average seismic potential of the slowly straining South Central Mediterranean plate boundary zone, we integrate constraints on tectonic style and deformation rates from geodetic and geologic data with the traditional constraints from seismicity catalogs. We express seismic potential (long-term average earthquake recurrence rates as a function of magnitude) in the form of truncated Gutenberg–Richter distributions for seven seismotectonic source zones. Seismic coupling seems to be large or even complete in most zones. An exception is the southern Tyrrhenian thrust zone, where most of the African–European convergence is accommodated. Here aseismic deformation is estimated to range from at least 25% along the western part to almost 100% aseismic slip around the Aeolian Islands. Even so, seismic potential of this zone has previously been significantly underestimated, due to the low levels of recorded past seismicity. By contrast, the series of 19 M6–7 earthquakes that hit Calabria in the 18th and 19th century released tectonic strain rates accumulated over time spans up to several times the catalog duration, and seismic potential is revised downward. The southern Tyrrhenian thrust zone and the extensional Calabrian faults, as well as the northeastern Sicilian transtensional zone between them (which includes the Messina Straits, where a destructive M7 event occurred in 1908), all have a similar seismic potential with minimum recurrence times of M ≥ 6.5 of 150–220 years. This potential is lower than that of the Southern Apennines (M ≥ 6.5 recurring every 60 to 140 years), but higher than that of southeastern Sicily (minimum M ≥ 6.5 recurrence times of 400 years). The high seismicity levels recorded in southeastern Sicily indicate some clustering and are most compatible with a tectonic scenario where the Ionian deforms internally, and motions at the Calabrian Trench are small. The estimated seismic potential for the Calabrian Trench and Central and Western Sicily are the lowest (minimum M ≥ 6.5 recurrence times of 550–800 years). Most zones are probably capable of generating earthquakes up to magnitudes 7–7.5, with the exception of Central and Western Sicily where maximum events sizes most likely do not exceed 7.     

Geophysical methods applied to fault characterization and earthquake potential assessment in the Lower Tagus Valley, Portugal

The study region is located in the Lower Tagus Valley, central Portugal, and includes a large portion of the densely populated area of Lisbon. It is characterized by a moderate seismicity with a diffuse pattern, with historical earthquakes causing many casualties, serious damage and economic losses. Occurrence of earthquakes in the area indicates the presence of seismogenic structures at depth that are deficiently known due to a thick Cenozoic sedimentary cover. The hidden character of many of the faults in the Lower Tagus Valley requires the use of indirect methodologies for their study. This paper focuses on the application of high-resolution seismic reflection method for the detection of near-surface faulting on two major tectonic structures that are hidden under the recent alluvial cover of the Tagus Valley, and that have been recognized on deep oil-industry seismic reflection profiles and/or inferred from the surface geology. These are a WNW–ESE-trending fault zone located within the Lower Tagus Cenozoic basin, across the Tagus River estuary (Porto Alto fault), and a NNE–SSW-trending reverse fault zone that borders the Cenozoic Basin at the W (Vila Franca de Xira–Lisbon fault). Vertical electrical soundings were also acquired over the seismic profiles and the refraction interpretation of the reflection data was carried out. According to the interpretation of the collected data, a complex fault pattern disrupts the near surface (first 400 m) at Porto Alto, affecting the Upper Neogene and (at least for one fault) the Quaternary, with a normal offset component. The consistency with the previous oil-industry profiles interpretation supports the location and geometry of this fault zone. Concerning the second structure, two major faults were detected north of Vila Franca de Xira, supporting the extension of the Vila Franca de Xira–Lisbon fault zone northwards. One of these faults presents a reverse geometry apparently displacing Holocene alluvium. Vertical offsets of the Holocene sediments detected in the studied geophysical data of Porto Alto and Vila Franca de Xira–Lisbon faults imply minimum slip rates of 0.15–0.30 mm/year, three times larger than previously inferred for active faults in the Lower Tagus Valley and maximum estimates of average return periods of 2000–5000 years for M 6.5–7 co-seismic ruptures.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

Optimization of nonconventional wells under uncertainty using statistical proxies

The determination of the optimal type and placement of a nonconventional well in a heterogeneous reservoir represents a challenging optimization problem. This determination is significantly more complicated if uncertainty in the reservoir geology is included in the optimization. In this study, a genetic algorithm is applied to optimize the deployment of nonconventional wells. Geological uncertainty is accounted for by optimizing over multiple reservoir models (realizations) subject to a prescribed risk attitude. To reduce the excessive computational requirements of the base method, a new statistical proxy (which provides fast estimates of the objective function) based on cluster analysis is introduced into the optimization process. This proxy provides an estimate of the cumulative distribution function (CDF) of the scenario performance, which enables the quantification of proxy uncertainty. Knowledge of the proxy-based performance estimate in conjunction with the proxy CDF enables the systematic selection of the most appropriate scenarios for full simulation. Application of the overall method for the optimization of monobore and dual-lateral well placement demonstrates the performance of the hybrid optimization procedure. Specifically, it is shown that by simulating only 10% or 20% of the scenarios (as determined by application of the proxy), optimization results very close to those achieved by simulating all cases are obtained.     

城市化过程中土地征用风险表征及对策

从土地征用基本概念和基本特征着手，描述了土地征用的主要风险，其表现为土地征用权力滥用、侵害失地农民的权益、诱发集体土地黑市、诱发社会矛盾等方面。分析了产生风险的“公共利益需要”缺乏明确界定、土地征用的补偿标准偏低或滞后、土地征用及相关制度不健全或存在缺陷等主要原因。提出了控制与规避风险的对策与措施：严格限定公益性用地范围，提高征地补偿标准，扩大补偿范围，完善相关制度，严格土地征用程序。     

Quantification of the risk associated with high levels of mercury pollution at Temirtau, Kazakhstan

This research aims to establish the risks of hurnan exposure to modeling approach. The study site was located on the River Nura 1990s, an acetaldehyde production plant located in the industrial high environmental mercury levels, based on a quantitative risk in central Kazakhstan. From the 1950s until its closure in the mid city of Temirtau released mercury-containing wastewater into the river. River sediments are highly contaminated with mercury up to a distance of 25 km and beyond. In addition, a local power station released an estimated 6 million tonnes of fly-ash into the water and the mercury has become quite tightly associated with the ash deposits. River water, fish and agricultural land in the floodplain are also contaminated with mercury, yet the risks posed to the local population have not been evaluated to date. In June and July 2005, we took samples of soil, interior and exterior dust, drinking water, and food from individual households, communal areas and markets. Additionally, water and sediment samples and fish were taken from the river. Interviews were conducted with householders to establish their age and body weight, general living conditions and sources of irrigation and drinking water. A food frequency questionnaire （FFQ） was designed to investigate the frequencies of consumption of several common regional food items, including fish from the fiver and/or local market. Human hair samples were also collected to estimate the mercury bioburden and to enable the validity of the modeling approach to be established. The paper expands on the main pathways of contamination and looks at linkages between exposure pathways and mercury concentrations found in human hair. Uncertainties inherent in risk analysis as well as their influence on the relative importance of different exposure routes are also discussed.