Observations of physical effects from tsunamis of December 30, 2002 at Stromboli volcano, southern Italy

On December 30, 2002, following an intense period of activity of Stromboli volcano (south Tyrrhenian Sea, Italy), complex mass failures occurred on the northwest slope of the mountain which also involved the underwater portion of the volcanic edifice for a total volume of about 2–3×10⁷ m³. Two main landslides occurred within a time separation of 7 min, and both set tsunami waves in motion that hit the coasts of Stromboli causing injuries to three people and severe damage to buildings and structures. The tsunamis also caused damage on the island of Panarea, some 20 km to the SSE from the source. They were observed all over the Aeolian archipelago, at the island of Ustica to the west, along the northern Sicily coasts to the south as well as along the Tyrrhenian coasts of Calabria to the east and in Campania to the north. This paper presents field observations that were made in the days and weeks immediately following the events. The results of the quantitative investigations undertaken in the most affected places, namely along the coasts of Stromboli and on the island of Panarea, are reported in order to highlight the dynamics of the attacking waves and their impact on the physical environment, on the coastal structures and on the coastal residential zone. In Stromboli, the tsunami waves were most violent along the northern and northeastern coastal belt between Punta Frontone and the village of Scari, with maximum runup heights of about 11 m measured on the beach of Spiaggia Longa. Measured runups were observed to decay rapidly with distance from the source, typical of tsunami waves generated by limited-area sources such as landslides.     

Comparison of two great Chile tsunamis in 1960 and 2010 using numerical simulation

Great differences in hazard and losses were shown from two tsunamis, both generated in Chile, one in 1960 and the other in 2010. Numerical simulation was applied to the tsunami analysis. The fault dislocation of the seafloor was assumed to equal to the initial tsunami wave field, which can be calculated by the formula of fault dislocation in the elastic isotropic half-space. The linear long wave theory was used as the tsunami hydrodynamic model, and the finite difference method and leap-frog scheme were selected for solving the equations. The accuracy of the simulated results was verified by the observed data in five tide gauges. By means of two scenario tsunamis, the analytical results show that the earthquake magnitude, bathymetry in rupture zone and rapid release of warning information in 2010 tsunami are the main explanations of the aforementioned great difference.     

The morphology and evolution of the Stromboli 2002–2003 lava flow field: an example of a basaltic flow field emplaced on a steep slope

The use of a hand-held thermal camera during the 2002–2003 Stromboli effusive eruption proved essential in tracking the development of flow field structures and in measuring related eruption parameters, such as the number of active vents and flow lengths. The steep underlying slope on which the flow field was emplaced resulted in a characteristic flow field morphology. This comprised a proximal shield, where flow stacking and inflation caused piling up of lava on the relatively flat ground of the vent zone, that fed a medial–distal lava flow field. This zone was characterized by the formation of lava tubes and tumuli forming a complex network of tumuli and flows linked by tubes. Most of the flow field was emplaced on extremely steep slopes and this had two effects. It caused flows to slide, as well as flow, and flow fronts to fail frequently, persistent flow front crumbling resulted in the production of an extensive debris field. Channel-fed flows were also characterized by development of excavated debris levees in this zone (Calvari et al. 2005). Collapse of lava flow fronts and inflation of the upper proximal lava shield made volume calculation very difficult. Comparison of the final field volume with that expecta by integrating the lava effusion rates through time suggests a loss of ~70% erupted lava by flow front crumbling and accumulation as debris flows below sea level. Derived relationships between effusion rate, flow length, and number of active vents showed systematic and correlated variations with time where spreading of volume between numerous flows caused an otherwise good correlation between effusion rate, flow length to break down. Observations collected during this eruption are useful in helping to understand lava flow processes on steep slopes, as well as in interpreting old lava–debris sequences found in other steep-sided volcanoes subject to effusive activity.     

海啸在冲绳海槽和东海浅水大陆架地形上产生和传播的数值模拟

运用数值模拟的方法对在冲绳海槽产生9.0级地震,并引发海啸的过程和海啸波在东海浅水大陆架地形上的传播过程进行研究.模拟的结果表明,数值模拟产生的波浪符合海啸波的特点,东海浅水大陆架适合海啸波的传播.     

The Event of 26th of December 2004 – The Biggest Earthquake in the World in the Last 40 years

A great earthquake occurred at 00:58:49 (UTC) on Sunday, December 26, 2004 off the northwest coast of Sumatra, Indonesia. Its revised moment magnitude was M 9.3 making it in the top four largest earthquakes in the world since 1900 and the largest since the Alaskan 1964 event. The earthquake caused tsunami waves which killed more than 300,000 people in Southern Asia and Africa. There were 31 earthquakes with magnitudes between 5.5 and 7.3 in the 48-h period after the main event, and it seemed that seismicity migrated northwards along the 1200 km fault (http: //www.ga.gov.au). Similar size events occurred in that location off Sumatra in the 19th century, but no evidence of written records of their tsunami effects in Australia is found. The devastating megathrust earthquake of 26 December 2004 occurred on the interface of the Indo-Australian and Euro-Asian plates where the first plate subducts beneath the overriding second plate and the Indo-Australian plate begins its descent into the mantle. In the epicentral region, the Indo-Australian plate moves toward the northeast at a rate of about 7 cm/year relative to the Euro-Asian plate resulting in oblique convergence and partitioning into thrust-faulting. From the size of the earthquake, it is likely that the displacement on the fault plane was up to fifteen meters. As with the recent event, megathrust earthquakes often generate large tsunamis that cause damage over a much wider area than is directly affected by ground shaking near the earthquake’s rupture. The subduction zone continues further south of the Indonesian archipelago and that area is also a potential risk of producing a megathrust event that may affect coastal parts of northwest Australia. The tragic events of Boxing Day 2004 highlighted the importance of establishing a tsunami warning system for the Indian Ocean like the one for the Pacific. Issues like more and better instrumentation, and a long-term program to educate people in the region about the dangers of tsunamis, were identified as priorities. Of particular interest is the time for identifying and issuing alerts for such devastating earthquakes with possibility to reduce it in future for warning purposes.     

The contemporaneous emission of low-K and high-K trachybasalts and the role of the NE Rift during the 2002 eruptive event,Mt. Etna,Italy

Mount Etna volcano erupted almost simultaneously on its northeastern and southern flanks between October 27 and November 3, 2002. The eruption on the northeastern flank lasted for 8 days, while on the southern flank it continued for 3 months. The northeastern flank eruption was characterized by the opening of a long eruptive fracture system between 2,900 and 1,900 m.a.s.l. A detailed survey indicates that the fractures’ direction shifted during the opening from N10W (at the NE Crater, 2,900 m) to N45E (at its lowest portion, 1,900 m) and that distinct magma groups were erupted at distinct fracture segments. Based on their petrological features, three distinct groups of rocks have been identified. The first group, high-potassium porphyritic (HKP), is made up of porphyritic lavas with a Porphyritic Index (P.I.) of 20–32 and K₂O content higher than 2 wt%. The second group is represented by lavas and tephra with low modal phenocryst abundance (P.I. < 20) named here oligo-phyric (low-phyric), and K₂O content higher than 2 wt% (HKO, high-potassium oligophyric). The third group, low-potassium oligophyric (LKO), consists of tephra with oligophyric texture (P.I. < 20) but K₂O content < 2 wt%. K-rich magmas (HKP and HKO) are similar to the magma erupted on the southern flank, and geochemical variations within these groups can be accounted for by a variable degree of fractionation from a single parent magma. The K-poor magma (LKO), erupted only in the upper segment of the fracture, cannot be placed on the same liquid line of descent of the HK groups, and it is similar to the magmas that fed the activity of Etna volcano prior to the eruption of 1971. This is the first time since then that a magma of this composition has been documented at Mt. Etna, thus providing a strong indication for the existence of distinct batches of magma whose rise and differentiation are independent from the main conduit system. The evolution of this eruption provides evidence that the NE Rift plays a very active role in the activity of Mt. Etna volcano, and that its extensional tectonics allows the intrusion and residence of magma bodies at various depths, which can therefore differentiate independently from the main open conduit system.     

日本Mw9.0级地震海啸数值模拟与启示

2011年3月11日13时46分日本东北部海域发生Mw9.0级特大地震,地震诱发了海啸.本文依据USGS的震源机制解,进行了地震海啸的数值模拟,并与美国国家海洋与大气管理局布设的海底压力计记录的水深数据对比分析,结果表明数值模拟结果可信.同时,进一步分析了海啸造成巨大损失的原因,并对未来我国海啸防灾减灾工作给出了几点建...     

Characterization of shape and terminal velocity of tephra particles erupted during the 2002 eruption of Etna volcano,Italy

In this paper, we present a complete morphological characterization of the ash particles erupted on 18 December 2002 from Etna volcano, Italy. The work is based on the acquisition and processing of bidimensional digital images carried out by scanning electron microscopy (SEM) to obtain shape parameters by image analysis. We measure aspect ratio (AR), form factor (FF), compactness (CC), and rectangularity (RT) of 2065 ash particles with size between 0.026 and 1.122 mm. We evaluate the variation of these parameters as a function of the grain-size. Ash particles with a diameter of <0.125 mm vary from mostly equant to very equant, ash particles between 0.125 and 0.250 mm have an intermediate shape, and particles with diameters >0.250 mm are subelongate. We find that, on average, particles with a diameter of <0.250 mm are subrounded, particles between 0.250 and 0.50 mm are subangular, and particles >0.50 mm are angular. Using this morphological analysis and an empirical relation between the drag coefficient (C _D) and Reynolds number (R _e) of Wilson and Huang (Earth Planet Sci Lett 44:311–324, 1979), we calculate the terminal settling velocities (V _WH). The comparisons between these velocities and those calculated with the formula of Kunii and Levenspiel (Fluidization engineering. Wiley, New York, pp 97, 1969) (V _KL), which considers ash particles as spheres, show that V _KL are in average 1.28 greater than V _WH. Hence, we quantify the systematic error on the spatial distribution of the mass computed around the volcano carried out by tephra dispersal models when the particles are assumed to be spherical.     

The Sumatra tsunami of 26 December 2004 as observed in the North Pacific and North Atlantic oceans

The M _w=9.3 megathrust earthquake of December 26, 2004 off the coast of Sumatra in the Indian Ocean generated a catastrophic tsunami that caused widespread damage in coastal areas and left more than 226,000 people dead or missing. The Sumatra tsunami was accurately recorded by a large number of tide gauges throughout the world's oceans. This paper examines the amplitudes, frequencies and wave train structure of tsunami waves recorded by tide gauges located more than 20,000 km from the source area along the Pacific and Atlantic coasts of North America.     

A comparative study of buried structure in soil subjected to blast load using 2D and 3D numerical simulations

The response of underground structures subjected to subsurface blast is an important topic in protective engineering. Due to various constraints, pertinent experimental data are extremely scarce. Adequately detailed numerical simulation thus becomes a desirable alternative. However, the physical processes involved in the explosion and blast wave propagation are very complex, hence a realistic and detailed reproduction of the phenomena would require sophisticated numerical models for the loading and material responses. In this paper, a fully coupled numerical model is used to simulate the response of a buried concrete structure under subsurface blast, with emphasis on the comparative performance of 2D and 3D modeling schemes. The explosive charge, soil medium and the RC structure are all incorporated in a single model system. The SPH (smooth particle hydrodynamics) technique is employed to model the explosive charge and the close-in zones where large deformation takes place, while the normal FEM is used to model the remaining soil region and the buried structure. Results show that the 2D model can provide reasonably accurate results concerning the crater size, blast loading on the structure, and the critical response in the front wall. The response in the remaining part of the structure shows noticeable differences between the 2D and 3D models. Based on the simulation results, the characteristics of the in-structure shock environment are also discussed in terms of the shock response spectra.     

Identification and mapping of riverbed sediment facies in the Columbia River through integration of field observations and numerical simulations

In the Hanford Reach of the Columbia River, a thin layer of recent alluvium overlies the sedimentary formations that comprise the unconfined groundwater aquifer. Experimental and modelling studies have demonstrated that this alluvial layer exerts significant control on the exchange of groundwater and surface water (hydrologic exchange flux), and is associated with elevated levels of biogeochemical activity. This layer is also observed to be strongly heterogeneous, and quantifying the spatial distribution of properties over the range of scales of interest is challenging. Facies are elements of a sediment classification scheme that groups complex geologic materials into a set of discrete classes according to distinguishing features. Facies classifications have been used as a framework for assigning heterogeneous material properties to grid cells of numerical models of flow and reactive transport in subsurface media. The usefulness of such an approach hinges on being able to relate facies to quantitative properties needed for flow and reactive transport modelling, and on being able to map facies over the domain of interest using readily available information. Although aquifer facies have been used in various modelling contexts, application of this concept to riverbed sediments is relatively new. Here, we describe an approach for categorizing and mapping recent alluvial (riverbed) sediments based on the integration of diverse observations with numerical simulations of river hydrodynamics. The facies have distinct distributions of sediment texture that correspond to variations in hydraulic properties, and therefore provide a useful framework for assigning heterogeneous properties in numerical simulations of hydrologic exchange flows and biogeochemical processes.     

Directional damage due to near-fault and site effects in the M6.4 Changureh–Avaj earthquake of 22 June 2002

The Changureh–Avaj earthquake (M = 6.4) occurred 220 km northwest of Tehran on June 22, 2002. According to the official report, 226 people were killed, 1300 injured and 33,000 houses were seriously damaged or collapsed. Two villages located in the near-fault region with a population more than 200 people were completely destroyed. The damage survey was conducted within 30 km of the epicenter. During damage investigation, an observable pattern of damage in toppled masonry and concrete block fences was recognized. The data of damaged fences were documented and analyzed to study the effect of direction and distance from the epicenter. From the 109 fences, 85 fences were built within 15° of either north or east, providing an opportunity to study the effect of orientation. The results of analysis revealed that the east–west oriented fences toppled 2.6 times more than those oriented north–south. More than 70% of toppled fences were oriented in east–west direction, with more falling towards the north than to the south. This preferred direction of damage for toppled fences was observed in all four quadrants about the epicenter and at distance up to 15 km. In the northern quadrants almost all toppled fences were oriented east–west. Two possibilities, namely, near-fault effect and 2D or 3D site response, are examined to find an explanation for the observed damage pattern at the sites. Since the preferred direction of damage inferred from toppled fences within 15 km of epicenter was found to be parallel to the direction of rupture propagation, the near-fault forwarded rupture directivity was investigated as the most possible reason for the observed phenomenon at the sites. However, there were no records of strong motions available in near-fault zone. Thus, the near-fault ground motions were theoretically simulated by using developed local source parameters and heterogeneous slip model. Next, the principal axes were calculated by using eigenvector analysis of simulated horizontal motions. The preferred direction of damage in over 90% of the sites is well predicted by the directional dependence of near-fault ground motion. Nevertheless, the preferred direction of damage in Avaj station, where the main shock was recorded, demonstrates totally different result from the one predicted by the source. The shear-wave polarization is investigated as an alternative reason at the sites. The polarization of motion is found to be the greatest at frequencies with peaks in maximized and H/V spectra, suggesting that directional motion is associated with local site and geologic conditions. The dominant direction of motion due to source mechanism and local site conditions was identified and compared with the preferred direction of damage observed at the sites. It is found that the near-fault source effects can strongly control the motion directionality and even dictate its dominant direction to the one caused by local geologic or site conditions. Finally, the outcome provides field evidence of directional damage due to near-fault effects and local site conditions. This may aid in predicting the directions of damaging earthquake motions, and serve as an important factor in the design of critical facilities.     

青海湖热力状况的模拟与未来情景之研究

运用基于湖泊热量收支与湍流扩散的湖泊热力学模型模拟了青海湖近３０年来湖水热力状况,内容包括湖泊水面温度,温度沿深度的垂直分布,冬季结冰与融冰的起迄时间,冰盖厚度;积雪深度与天数等。在此基础上,运用４个ＧＣＭｓ模型输出的ＣＯ２倍增情景下该地区的气候状况,评价了湖泊热力状况在未来的可能变化情况。