Using numerical modelling to evaluate tsunami hazard near the Kuril Island

A sequence of computer experiments is used to study questions concerning the tsunami problem as a quantitative estimate of tsunami danger, detailed geographical tsunami classification, determination of the parameters of critical tsunami waves, and the conditions of their development. We call a wave critical, if its impact on the coast is most hazardous.Using the Middle Kuril Island as an example, we present the results of a computer experiment which includes determining the wavefields on the shelf and estimating the effects connected with the deep-water Bussol and Diana Straits.Numerical simulation of tsunami waves of different sources permits the assessment of the extent of tsunami danger in different areas of the coastal zone of Simushir Island, depending on the location of the focus zone and their geometry.The major singularities of the wavefield arise in the zones of the deep-water straits. The distribution of the amplification factors is determined by both the global parameters of the wavefields and the local properties of individual harbours. The results obtained for a particular harbour in the northern part of Simushir Island, formed the basis for the quantitative estimate of tsunami danger for this area.     

Na depletion in modern adakites via melt / rock reaction within the sub-arc mantle

Interaction between slab-derived melt and mantle peridotite and the role of slab melt as a metasomatizing agent in the sub-arc mantle is being increasingly recognized. Adakite, the slab melt erupted on the surface, usually exhibits anomalously high MgO, CaO, Cr and Ni contents that indicate interaction with mantle peridotitite. Here we note that Cenozoic adakites have Na₂O contents below 5.8 wt.% with ∼95% samples lower than 5.0 wt.%, and are generally depleted in this component relative to experimental basalt partial melts (mostly beyond 5.0 wt.% and up to 9.0 wt.% Na₂O) produced under 1.5-3.0 GPa conditions that are most relevant to adakite production. We interpret the adakite Na depletion to be also a consequence of the melt / rock reaction that takes place within the hot mantle wedge. During ascent and reaction with mantle peridotite, primary adakite melts gain mantle components MgO, CaO, Cr and Ni but lose Na₂O, SiO₂ and perhaps K₂O to the mantle, leading to Na-rich mantle metasomatism. Selective assimilation of predominately mantle clinopyroxene, some spinel and minor olivine at high T/P has been considered to be an important process in producing high-Mg adakites from primary low-Mg slab melts [Killian, R., Stern, C. R., 2002. Constraints on the interaction between slab melts and the mantle wedge from adakitic glass in peridotite xenoliths. Eur. J. Mineral. 14, 25-36]. In such a process, Na depletion in the assimilated melt is the result of dilution due to the increase in melt mass. Phase relationships in the reaction system siliceous melt + peridotite and quantitative calculation suggest that assimilation of mantle clinopyroxene, olivine and spinel and fractional crystallization of sodic amphibole and orthopyroxene, under conditions of moderate T/P and increasing melt mass, is also an important process that modifies the composition of adakites and causes the Na depletion.     

Aftershocks series monitoring of the September 18, 2004 M = 4.6 earthquake at the western Pyrenees: A case of reservoir-triggered seismicity?

On September 18, 2004, a 4.6 mbLg earthquake was widely felt in the region around Pamplona, at the western Pyrenees. Preliminary locations reported an epicenter less than 20 km ESE of Pamplona and close to the Itoiz reservoir, which started impounding in January 2004. The area apparently lacks of significant seismic activity in recent times. After the main shock, which was preceded by series of foreshocks reaching magnitudes of 3.3 mbLg, a dense temporal network of 13 seismic stations was deployed there to monitor the aftershocks series and to constrain the hypocentral pattern. Aftershock determinations obtained with a double-difference algorithm define a narrow epicentral zone of less than 10 km², ESE–WNW oriented. The events are mainly concentrated between 3 and 9 km depth. Focal solutions were computed for the main event and 12 aftershocks including the highest secondary one of 3.8 mbLg. They show mainly normal faulting with some strike-slip component and one of the nodal planes oriented NW–SE and dipping to the NE. Cross-correlation techniques applied to detect and associate events with similar waveforms, provided up to 33 families relating the 67% of the 326 relocated aftershocks. Families show event clusters grouped by periods and migrating from NW to SE. Interestingly, the narrow epicentral zone inferred here is located less than 4 km away from the 111-m high Itoiz dam. These hypocentral results, and the correlation observed between fluctuations of the reservoir water level and the seismic activity, favour the explanation of this foreshock–aftershock series as a rapid response case of reservoir-triggered seismicity, burst by the first impoundment of the Itoiz reservoir. The region is folded and affected by shallow dipping thrusts, and the Itoiz reservoir is located on the hangingwall of a low angle southward verging thrust, which might be a case sensible to water level fluctuations. However, continued seismic monitoring in the coming years is mandatory in this area to infer more reliable seismotectonic and hazard assessments.     

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.     

Evidence for a seismic activity mainly constituted of hybrid events at Cayambe volcano,Ecuador. Interpretation in a iced-domes volcano context

The high activity level of Hybrid Events (HE) detected beneath the Cayambe volcano since 1989 has been more thoroughly investigated with data from a temporary array. The unusual HE spectral content allows separating a high-frequency signal riding on a low-frequency one, with a probable single source. HEs are interpreted as high frequency VT events, produced by the interaction between magmatic heat and an underground water system fed by thaw water from the summital glacier, which trigger simultaneous low-frequency fluid resonance in the highly fractured adjacent medium. Pure VTs are interpreted as ‘aborted’ HEs occurring probably in the oldest and coldest part of the volcano complex. To cite this article: B. Guillier, J.-L. Chatelain, C. R. Geoscience 338 (2006).     

Great earthquakes of variable magnitude at the Cascadia subduction zone

Comparison of histories of great earthquakes and accompanying tsunamis at eight coastal sites suggests plate-boundary ruptures of varying length, implying great earthquakes of variable magnitude at the Cascadia subduction zone. Inference of rupture length relies on degree of overlap on radiocarbon age ranges for earthquakes and tsunamis, and relative amounts of coseismic subsidence and heights of tsunamis. Written records of a tsunami in Japan provide the most conclusive evidence for rupture of much of the plate boundary during the earthquake of 26 January 1700. Cascadia stratigraphic evidence dating from about 1600 cal yr B.P., similar to that for the 1700 earthquake, implies a similarly long rupture with substantial subsidence and a high tsunami. Correlations are consistent with other long ruptures about 1350 cal yr B.P., 2500 cal yr B.P., 3400 cal yr B.P., 3800 cal yr B.P., 4400 cal yr B.P., and 4900 cal yr B.P. A rupture about 700-1100 cal yr B.P. was limited to the northern and central parts of the subduction zone, and a northern rupture about 2900 cal yr B.P. may have been similarly limited. Times of probable short ruptures in southern Cascadia include about 1100 cal yr B.P., 1700 cal yr B.P., 3200 cal yr B.P., 4200 cal yr B.P., 4600 cal yr B.P., and 4700 cal yr B.P. Rupture patterns suggest that the plate boundary in northern Cascadia usually breaks in long ruptures during the greatest earthquakes. Ruptures in southernmost Cascadia vary in length and recurrence intervals more than ruptures in northern Cascadia.     

New Approach for Estimation of Static and Seismic Active Earth Pressure

To estimate static and seismic active earth pressure (P_ad) on a rigid retaining wall, numerical analyses using different step sizes have been carried out in this paper, based on the modified Culmann line method by considering Coulomb’s planar rupture surface. Equivalent pseudo-static seismic forces are considered in the analysis. A new concept of modified unit weight by considering ground surcharge is introduced under static and seismic conditions. By numerical analysis, area of soil (A) has been estimated to obtain the ratio of A/A₀ where A₀ is θh², θ is the angle between retaining structure and ground surface and h is the vertical height of the wall. This ratio remains constant for a particular type of soil and has been used to estimate the maximum active earth pressure using force diagram. Results are provided in tabular form for easy calculation of the coefficient of static and seismic active earth pressure. Present results by considering the new technique, compares well with the results obtained by earlier researchers.     

Centrifuge modelling-based verification of response spectra design methods for some vulnerable structures

Residential RC framed structures suffered heavily during the 2001 Bhuj earthquake in Gujarat, India. These types of structures also saw severe damage in other earthquakes such as the 1999 Kocaeli earthquake in Turkey and 921 Ji-Ji earthquake in Taiwan. In this paper the seismic response of residential structures was investigated using physical modelling. Idealised soft storey and top heavy, two degrees of freedom (2DOF) portal frame structures were developed and tested on saturated and dry sand models at 25 g using the Schofield Centre 10-m Beam Centrifuge. It was possible to recreate observed field behaviour using these models. As observed in many of the recent earthquakes, soft storey structures were found to be particularly vulnerable to seismic loads. Elastic response spectra methods are often used in the design of simple portal frame structures. The seismic risk of these structures can be significantly increased due to modifications such as removal of a column or addition of heavy water tanks on the roof. The experimental data from the dynamic centrifuge tests on such soft storey or top-heavy models was used to evaluate the predictions obtained from the response spectra. Response spectra were able to predict seismic response during small to moderate intensity earthquakes, but became inaccurate during strong earthquakes and when soil structure interaction effects became important. Re-evaluation of seismic risk of such modified structures is required and time domain analyses suggested by building codes such as IBC, UBC or NEHRP may be more appropriate.     

用现今小震推断洪洞、临汾两次历史大震的震源断层

山西临汾地区是一个历史强震多发区，1303年和1695年发生了洪洞（M＝8）和临汾（M＝73／4）两次特大地震，这两次地震所在区域至今仍在持续不断的小地震活动，具有明显的大震区地震长期活动特征，我们对临汾无线传输地震台网记录的1987－1999年期间发生的1670次中，小地震重新进行了震源定位，根据对这些地震震源位置三维空间分布特征和震源机解制的分析，认为洪洞地震的震源断层应是长80km，埋深5－26km的NNE走向，高倾角的右旋走滑型断层，而临汾地震的震源断层是长70km，埋深5－22km的NWW走向，高倾角的左旋走滑型断层。这与洪洞，临汾两次大震极震区的等震线及该地区应力场的构造环境是吻合的。