The landslides and tsunamis of the 30th of December 2002 in Stromboli analysed through numerical simulations

On the 30th of December 2002 two tsunamis were generated only 7 min apart in Stromboli, southern Tyrrhenian Sea, Italy. They represented the peak of a volcanic crisis that started 2 days before with a large emission of lava flows from a lateral vent that opened some hundreds of meters below the summit craters. Both tsunamis were produced by landslides that detached from the Sciara del Fuoco. This is a morphological scar and is the result of the last collapse of the northwestern flank of the volcanic edifice, that occurred less than 5 ka b.p. The first tsunami was due to a submarine mass movement that started very close to the coastline and that involved about 20×10⁶ m³ of material. The second tsunami was engendered by a subaerial landslide that detached at about 500 m above sea level and that involved a volume estimated at 4–9×10⁶ m³. The latter landslide can be seen as the retrogressive continuation of the first failure. The tsunamis were not perceived as distinct events by most people. They attacked all the coasts of Stromboli within a few minutes and arrived at the neighbouring island of Panarea, 20 km SSW of Stromboli, in less than 5 min. The tsunamis caused severe damage at Stromboli.In this work, the two tsunamis are studied by means of numerical simulations that use two distinct models, one for the landslides and one for the water waves. The motion of the sliding bodies is computed by means of a Lagrangian approach that partitions the mass into a set of blocks: we use both one-dimensional and two-dimensional schemes. The landslide model calculates the instantaneous rate of the vertical displacement of the sea surface caused by the motion of the underwater slide. This is included in the governing equations of the tsunami, which are solved by means of a finite-element (FE) technique. The tsunami is computed on two different grids formed by triangular elements, one covering the near-field around Stromboli and the other also including the island of Panarea.The simulations show that the main tsunamigenic potential of the slides is restricted to the first tens of seconds of their motion when they interact with the shallow-water coastal area, and that it diminishes drastically in deep water. The simulations explain how the tsunamis that are generated in the Sciara del Fuoco area, are able to attack the entire coastline of Stromboli with larger effects on the northern coast than on the southern. Strong refraction and bending of the tsunami fronts is due to the large near-shore bathymetric gradient, which is also responsible for the trapping of the waves and for the persistence of the oscillations. Further, the first tsunami produces large waves and runup heights comparable with the observations. The simulated second tsunami is only slightly smaller, though it was induced by a mass that is approximately one third of the first. The arrival of the first tsunami is negative, in accordance with most eyewitness reports. Conversely, the leading wave of the second tsunami is positive.     

Sector collapse at Kick 'em Jenny submarine volcano (Lesser Antilles): numerical simulation and landslide behaviour

Kick 'em Jenny volcano is the only known active submarine volcano in the Lesser Antilles. It lies within a horseshoe-shaped structure open to the west northwest, toward the deep Grenada Basin. A detailed bathymetric survey of the basin slope at Kick 'em Jenny and resulting high-resolution digital elevation model allowed the identification of a major submarine landslide deposit. This deposit is thought to result from a single sector collapse event at Kick 'em Jenny and to be linked to the formation of the horseshoe-shaped structure. We estimated the volume and the leading-edge runout of the landslide to be ca. 4.4 km³ and 14 km, respectively. We modelled a sector collapse event of a proto Kick 'em Jenny volcano using VolcFlow, a finite difference code based on depth-integrated mass and momentum equations. Our models show that the landslide can be simulated by either a Coulomb-type rheology with low basal friction angles (5.5°–6.5°) and a significant internal friction angle (above 17.5°) or, with better results, by a Bingham rheology with low Bingham kinematic viscosity (0 < ν _B < 30 m²/s) and high shear strength (130 < γ ≤ 180 m²/s²). The models and the short runout distance suggest that the landslide travelled as a stiff cohesive flow affected by minimal granular disaggregation and slumping on a non-lubricated surface. The main submarine landslide deposit can therefore be considered as a submarine mass slide deposit that behaved like a slump.     

FDM Simulation of Seismic Waves, Ocean Acoustic Waves, and Tsunamis Based on Tsunami-Coupled Equations of Motion

We have developed a new, unified modeling technique for the total simulation of seismic waves, ocean acoustic waves, and tsunamis resulting from earthquakes, based on a finite difference method simulation of the 3D equations of motion. Using the equilibrium between the pressure gradient and gravity in these equations, tsunami propagation is naturally incorporated in the simulation based on the equations of motion. The performance of the parallel computation for the newly developed tsunami-coupled equations using a domain partitioning procedure shows a high efficiency coefficient with a large number of CPU cores. The simulation results show how the near-field term associated with seismic waves produced by shallow earthquakes leads to a permanent coseismic deformation of the ground surface, which gives rise to the initial tsunami on the sea surface. Propagation of the tsunami along the sea surface as a gravity wave, and ocean acoustic waves in seawater with high-frequency multiple P-wave reflections between the free surface and sea bottom, are also clearly demonstrated by the present simulations. We find a good agreement in the tsunami waveform between our results and those obtained by other simulations based on an analytical model and the Navier–Stokes equations, demonstrating the effectiveness of the tsunami-coupling simulation model. Based on this simulation, we show that the ratio of the amplitude of ocean acoustic waves to the height of the tsunami, both of which are produced by the earthquake, strongly depends on the rise time of the earthquake rupture. This ratio can be used to obtain a more detailed understanding of the source rupture processes of subduction zone earthquakes, and for implementing an improved tsunami alert system for slow tsunami earthquakes.     

Multiple edifice failures, debris avalanches and associated eruptions in the Holocene history of Shiveluch volcano, Kamchatka, Russia

 Investigation of well-exposed volcaniclastic deposits of Shiveluch volcano indicates that large-scale failures have occurred at least eight times in its history: approximately 10,000, 5700, 3700, 2600, 1600, 1000, 600 ¹⁴C BP and 1964 AD. The volcano was stable during the Late Pleistocene, when a large cone was formed (Old Shiveluch), and became unstable in the Holocene when repetitive collapses of a portion of the edifice (Young Shiveluch) generated debris avalanches. The transition in stability was connected with a change in composition of the erupting magma (increased SiO₂ from ca. 55–56% to 60–62%) that resulted in an abrupt increase of viscosity and the production of lava domes. Each failure was triggered by a disturbance of the volcanic edifice related to the ascent of a new batch of viscous magma. The failures occurred before magma intruded into the upper part of the edifice, suggesting that the trigger mechanism was indirectly associated with magma and involved shaking by a moderate to large volcanic earthquake and/or enhancement of edifice pore pressure due to pressurised juvenile gas. The failures typically included: (a) a retrogressive landslide involving backward rotation of slide blocks; (b) fragmentation of the leading blocks and their transformation into a debris avalanche, while the trailing slide blocks decelerate and soon come to rest; and (c) long-distance runout of the avalanche as a transient wave of debris with yield strength that glides on a thin weak layer of mixed facies developed at the avalanche base. All the failures of Young Shiveluch were immediately followed by explosive eruptions that developed along a similar pattern. The slope failure was the first event, followed by a plinian eruption accompanied by partial fountain collapse and the emplacement of pumice flows. In several cases the slope failure depressurised the hydrothermal system to cause phreatic explosions that preceded the magmatic eruption. The collapse-induced plinian eruptions were moderate-sized and ordinary events in the history of the volcano. No evidence for directed blasts was found associated with any of the slope failures. Received: 28 June 1998 / Accepted: 28 March 1999     

Application of Analytical Modeling to the Far-field Investigation of Tsunami and Oceanic Rayleigh Waves from the 1998 Papua New Guinea Earthquake

We analyze far-field Rayleigh and tsunami waves generated by the 1998 Papua New Guinea (PNG) earthquake. Using the normal mode theory and Thomson-Haskell matrix formalism we calculate synthetic mareograms of oceanic surface waves excited by finite-dimensional line source and propagated in a flat, multilayered oceanic structure. Assuming that the source of destructive sea waves was the main shock of the PNG event and based on the expression for seismic wave displacement in the far-field zone, we compute the energy of the seismic and tsunami waves and the E_z /E_ts ratio. The results of our modeling are generally consistent with those obtained empirically for events with magnitude 7. Also, treating the results of a submarine slide as a single solitary wave and using the theoretical arguments of Striem and Miloh (1976) we estimate the energy of the tsunami induced by a landslide. The difference between the energy of the seismic tsunami and of the aseismic one is about one order of magnitude.The results of our theoretical modeling show that surface sea waves in the far-field zone account well for seismic origin, although additional tsunami energy from a landslide source could be required to explain the local massive tsunami in the Sissano Lagoon.     

Combined Effects of Tectonic and Landslide-Generated Tsunami Runup at Seward, Alaska During the M W 9.2 1964 Earthquake

We apply a recently developed and validated numerical model of tsunami propagation and runup to study the inundation of Resurrection Bay and the town of Seward by the 1964 Alaska tsunami. Seward was hit by both tectonic and landslide-generated tsunami waves during the $M_{\rm W}$ 9.2 1964 megathrust earthquake. The earthquake triggered a series of submarine mass failures around the fjord, which resulted in landsliding of part of the coastline into the water, along with the loss of the port facilities. These submarine mass failures generated local waves in the bay within 5?min of the beginning of strong ground motion. Recent studies estimate the total volume of underwater slide material that moved in Resurrection Bay to be about 211?million m³ (Haeussler et?al. in Submarine mass movements and their consequences, pp 269?C278, 2007). The first tectonic tsunami wave arrived in Resurrection Bay about 30?min after the main shock and was about the same height as the local landslide-generated waves. Our previous numerical study, which focused only on the local landslide-generated waves in Resurrection Bay, demonstrated that they were produced by a number of different slope failures, and estimated relative contributions of different submarine slide complexes into tsunami amplitudes (Suleimani et?al. in Pure Appl Geophys 166:131?C152, 2009). This work extends the previous study by calculating tsunami inundation in Resurrection Bay caused by the combined impact of landslide-generated waves and the tectonic tsunami, and comparing the composite inundation area with observations. To simulate landslide tsunami runup in Seward, we use a viscous slide model of Jiang and LeBlond (J Phys Oceanogr 24(3):559?C572, 1994) coupled with nonlinear shallow water equations. The input data set includes a high resolution multibeam bathymetry and LIDAR topography grid of Resurrection Bay, and an initial thickness of slide material based on pre- and post-earthquake bathymetry difference maps. For simulation of tectonic tsunami runup, we derive the 1964 coseismic deformations from detailed slip distribution in the rupture area, and use them as an initial condition for propagation of the tectonic tsunami. The numerical model employs nonlinear shallow water equations formulated for depth-averaged water fluxes, and calculates a temporal position of the shoreline using a free-surface moving boundary algorithm. We find that the calculated tsunami runup in Seward caused first by local submarine landslide-generated waves, and later by a tectonic tsunami, is in good agreement with observations of the inundation zone. The analysis of inundation caused by two different tsunami sources improves our understanding of their relative contributions, and supports tsunami risk mitigation in south-central Alaska. The record of the 1964 earthquake, tsunami, and submarine landslides, combined with the high-resolution topography and bathymetry of Resurrection Bay make it an ideal location for studying tectonic tsunamis in coastal regions susceptible to underwater landslides.     

Numerical Study of Tsunami Generated by Multiple Submarine Slope Failures in Resurrection Bay, Alaska, during the M W 9.2 1964 Earthquake

We use a viscous slide model of Jiang and LeBlond (1994) coupled with nonlinear shallow water equations to study tsunami waves in Resurrection Bay, in south-central Alaska. The town of Seward, located at the head of Resurrection Bay, was hit hard by both tectonic and local landslide-generated tsunami waves during the M _W 9.2 1964 earthquake with an epicenter located about 150 km northeast of Seward. Recent studies have estimated the total volume of underwater slide material that moved in Resurrection Bay during the earthquake to be about 211 million m³. Resurrection Bay is a glacial fjord with large tidal ranges and sediments accumulating on steep underwater slopes at a high rate. Also, it is located in a seismically active region above the Aleutian megathrust. All these factors make the town vulnerable to locally generated waves produced by underwater slope failures. Therefore it is crucial to assess the tsunami hazard related to local landslide-generated tsunamis in Resurrection Bay in order to conduct comprehensive tsunami inundation mapping at Seward. We use numerical modeling to recreate the landslides and tsunami waves of the 1964 earthquake to test the hypothesis that the local tsunami in Resurrection Bay has been produced by a number of different slope failures. We find that numerical results are in good agreement with the observational data, and the model could be employed to evaluate landslide tsunami hazard in Alaska fjords for the purposes of tsunami hazard mitigation.     

Re-examination of the Source Mechanism of the 1998 Papua New Guinea Earthquake and Tsunami

— Simulation of tsunami propagation and runup of the 1998 Papua New Guinea (PNG) earthquake tsunami using the detailed bathymetry measured by JAMSTEC and adding bathymetric data at depths less than 60 m is carried out, reproducing the tsunami energy focus into Warapu and Arop along the Sissano Lagoon. However, the computed runup heights in the lagoon are still lower than those measured. Even if the error in estimating the fault parameters is taken into consideration, computational results are similar. Analysis by the wave ray method using several scenarios of the source size of the tsunami and location by the wave ray method suggests that a source characterized by small size in water 1,000-m deep approximately 25 km offshore the lagoon, best fits the arrival determined from the interviews with eyewitnesses. A two-layer numerical model simulating the interaction of the tsunami with a landslide is employed to study the behavior of a landslide-generated tsunami with different size sand depths of the initial slide just outside the lagoon. A landslide model with a volume of 4–8 × 10⁹ m³ is selected as the best in order to reproduce the distribution of the measured tsunami runup in the lagoon. The simulation of a tsunami generated in two stages, fault and landslide, could show good agreement with the runup heights and distribution of the arrival time, but a time gap of around 10 minutes remains, suggesting that a tsunami generated by the mainshock at 6:49 PM local time is too small for people to notice, and the following tsunami triggered by landslide or mass movement near the lagoon about ten minutes after the mainshock attacked the coast and caused the huge damage.     

Multiple sector collapses at stromboli volcano, Italy: how they work

This paper demonstrates that four large sector collapses have affected the NW flank of the Stromboli volcano in the past 13 ka, alternating with growth phases, in order to contribute to the evaluation of the critical conditions which trigger lateral collapses, a reconstruction of the geometry of each collapse of the volcano edifice in the four stages that preceded the relative collapse events is also presented, and a computation of the landslide volume. This reconstruction is based on new field data plotted in three dimensions. Prior to the initial 13-ka collapse, the volcano was 1125±100 m high above sea level. The collapse had a volume of 2.23±0.87 km³, whereas the pre-collapse volcano volume was 218.8±7.7 km³. The next edifice that failed was 900±70 m high a.s.l. The collapse volume was 1±0.54 km³, with a precollapse volcano volume of 201.4±5.4 km³. The edifice then grew to 1000±60 m a.s.l. The third collapse had a volume of 1.08±0.39 km³ and occurred within a volcano with a volume of 209.1±4.6 km³. This was followed by a new growth phase followed by the last collapse with a volume of 0.73±0.22 km³. The volcano volume was about the same as the present one. The present active crater zone is at 780 m a.s.l. in the first three collapses, sliding surfaces cut the main magma conduit. In the last collapse, the upper scarp coincided with the conduit location. Dyking along a main NE-trending weakness zone across the volcano summit exerted a lateral force for collapse inception. The decrease of the landslide volumes with the age, and the concentric scarps of the four collapses, suggest that the younger sliding planes tended to become more superficial and to decrease the areal extent. This is interpreted as due to: (a) successively weaker eruptive products from dominantly lavas to dominantly pyroclastics; and (b) the feedback effects between collapses and dykes that injected along the lateral segments of the first collapse slide plane.     

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.     

The 1998 debris avalanche at Casita volcano, Nicaragua — investigation of structural deformation as the cause of slope instability using remote sensing

During Hurricane Mitch in 1998, a debris avalanche occurred at Casita volcano, Nicaragua, resulting in a lahar that killed approximately 2500 people. The failure that initiated the avalanche developed at a pre-existing cliff, part of the headwall of a gravitational slide of approximately 1.8 km² in plan view that cuts the southern flank of the volcano. Structural analysis, primarily based on a high-resolution DEM, has shown that this slide is caused by edifice deformation. Casita's eastern side is spreading radially outwards, forming a convex–concave profile and steepening original slopes. This deformation is possibly facilitated by millennia of persistent hydrothermal alteration of the volcano's core. The gravity slide has some typical features of smaller slumps, such as steep headwalls, an inner flatter area and a pronounced basal bulge fronted by thrusts. The headwall is the source of the 1998 avalanche, as well as several previous mass movements. Edifice deformation has led to extensive fracturing of the hydrothermally altered andesitic source rock, increasing instability further. Field evidence indicates that the gravity slide is still actively deforming, and with steep headscarps remaining, the hazard of future avalanches is increasing. The analysis presented here shows how small but highly damaging landslides can occur during the deformation of a volcanic edifice. We show that identification of instability is possible with remote sensing data and minimal reconnaissance work, implying the possibility of similar efficient and cost-effective analysis at other volcanoes known to host extensive hydrothermal systems. We demonstrate this with a simple structural analysis of two similar stratovolcanoes, Orosí (Costa Rica) and Maderas (Nicaragua).     

Volcanic tsunamis and prehistoric cultural transitions in Cook Inlet,Alaska

The 1883 eruption of Augustine Volcano produced a tsunami when a debris avalanche traveled into the waters of Cook Inlet. Older debris avalanches and coeval paleotsunami deposits from sites around Cook Inlet record several older volcanic tsunamis. A debris avalanche into the sea on the west side of Augustine Island ca. 450 years ago produced a wave that affected areas 17 m above high tide on Augustine Island. A large volcanic tsunami was generated by a debris avalanche on the east side of Augustine Island ca. 1600 yr BP, and affected areas more than 7 m above high tide at distances of 80 km from the volcano on the Kenai Peninsula. A tsunami deposit dated to ca. 3600 yr BP is tentatively correlated with a southward directed collapse of the summit of Redoubt Volcano, although little is known about the magnitude of the tsunami. The 1600 yr BP tsunami from Augustine Volcano occurred about the same time as the collapse of the well-developed Kachemak culture in the southern Cook Inlet area, suggesting a link between volcanic tsunamis and prehistoric cultural changes in this region of Alaska.     

Lituya Bay Landslide Impact Generated Mega-Tsunami 50th Anniversary

On July 10, 1958, an earthquake M_w 8.3 along the Fairweather fault triggered a major subaerial landslide into Gilbert Inlet at the head of Lituya Bay on the southern coast of Alaska. The landslide impacted the water at high speed generating a giant tsunami and the highest wave runup in recorded history. The mega-tsunami runup to an elevation of 524 m caused total forest destruction and erosion down to bedrock on a spur ridge in direct prolongation of the slide axis. A cross section of Gilbert Inlet was rebuilt at 1:675 scale in a two-dimensional physical laboratory model based on the generalized Froude similarity. A pneumatic landslide tsunami generator was used to generate a high-speed granular slide with controlled impact characteristics. State-of-the-art laser measurement techniques such as particle image velocimetry (PIV) and laser distance sensors (LDS) were applied to the decisive initial phase with landslide impact and wave generation as well as the runup on the headland. PIV provided instantaneous velocity vector fields in a large area of interest and gave insight into kinematics of wave generation and runup. The entire process of a high-speed granular landslide impact may be subdivided into two main stages: (a) Landslide impact and penetration with flow separation, cavity formation and wave generation, and (b) air cavity collapse with landslide run-out and debris detrainment causing massive phase mixing. Formation of a large air cavity — similar to an asteroid impact — in the back of the landslide is highlighted. A three-dimenional pneumatic landslide tsunami generator was designed, constructed and successfully deployed in the tsunami wave basin at OSU. The Lituya Bay landslide was reproduced in a three-dimensional physical model at 1:400 scale. The landslide surface velocities distribution was measured with PIV. The measured tsunami amplitude and runup heights serve as benchmark for analytical and numerical models.     

Seismicity and magma supply rate of the 1998 failed eruption at Iwate volcano,Japan

Iwate volcano, Japan, showed significant volcanic activity including earthquake swarms and volcano inflation from the beginning of 1998. A large earthquake of magnitude 6.1 hit the south-west of the volcano on September 3. Although a 1 km² fumarole field formed, blighting plants on the ridge in the western part of the volcano in the spring of 1999, no magmatic eruptions occurred. We reconcile the spatio-temporal distributions of volcanic pressure sources determined by previously reported studies in which GPS, strain and tilt data from dense geodetic station networks are analyzed (Miura et al. Earth Planet Space 52:1003–1008, 2000; Sato and Hamaguchi J Volcanol Geotherm Res 155:244–262, 2006). We calculate the magma supply rates from their results and compare them with the occurrence rates of volcanic earthquakes. The results show that the magma supply rates are almost constant or even decrease with time while the earthquake occurrence rate increases with time. This contrast in their temporal changes is interpreted to result from stress accumulation in the volcanic edifice caused by constant magma supply without effusion of magma to the surface. We further show that data showing slight acceleration in strain can be best explained by magma ascent at a constant velocity, and that there is no evidence for increased magma buoyancy resulting from gas bubble growth. This consideration supports the interpretation that the magma stayed at 2 km depth and horizontally migrated. These findings relating magma supply rate and seismicity to magma ascent process are clues to understanding why no magmatic eruption occurred at Iwate volcano in 1998.     

Some Aspects of Energy Balance and Tsunami Generation by Earthquakes and Landslides

— Tsunamis are generated by displacement or motion of large volumes of water. While there are several documented cases of tsunami generation by volcanic eruptions and landslides, most observed tsunamis are attributed to earthquakes. Kinematic models of tsunami generation by earthquakes — where specified fault size and slip determine seafloor and sea-surface vertical motion — quantitatively explain far-field tsunami wave records. On the other hand, submarine landslides in subduction zones and other tectonic settings can generate large tsunamis that are hazardous along near-source coasts. Furthermore, the ongoing exploration of the oceans has found evidence for large paleo-landslides in many places, not just subduction zones. Thus, we want to know the relative contribution of faulting and landslides to tsunami generation. For earthquakes, only a small fraction of the minimum earthquake energy (less than 1% for typical parameter choices for shallow underthrusting earthquakes) can be converted into tsunami wave energy; yet, this is enough energy to generate terrible tsunamis. For submarine landslides, tsunami wave generation and landslide motion interact in a dynamic coupling. The dynamic problem of a 2-D translational slider block on a constant-angle slope can be solved using a Green's function approach for the wave transients. The key result is that the largest waves are generated when the ratio of initial water depth above the block to downslope vertical drop of the block H ₀ /W sin δ is less than 1. The conversion factor of gravitational energy into tsunami wave energy varies from 0% for a slow-velocity slide in deep water, to about 50% for a fast-velocity slide in shallow water and a motion abruptly truncated. To compare maximum tsunami wave amplitudes in the source region, great earthquakes produce amplitudes of a few meters at a wavelength fixed by the fault width of 100 km or so. For submarine landslides, tsunami wave heights — as measured by b, block height — are small for most of the parameter regime. However, for low initial dynamic friction and values of H ₀ /W sin δ less than 1, tsunami wave heights in the downslope and upslope directions reach b and b/4, respectively.Wavelengths of these large waves scale with block width. For significant submarine slides, the value of b can range from meters up to the kilometer scale. Thus, the extreme case of efficient tsunami generation by landslides produces dramatic hazards scenarios.     

Tsunami Excitation by Submarine Slides in Shallow-water Approximation

-- Landslide-induced tsunamis are receiving increased attention since there is evidence that recent large devastating events have been caused by underwater mass failures. Normally, numerical models are used to simulate tsunami excitation, most of which are based on shallow water, known also as long wave, approximation to the full equations of hydrodynamics. Analytical studies may handle only simplified problems, but help understand the basic features of physical processes. This paper is an analytical investigation of long-water waves excited by rigid bodies sliding on the sea bottom, based on the shallow-water approximation, which is here derived by properly scaling Euler equations for an inviscid, incompressible and irrotational ocean. In one-dimensional (1-D) cases (where motion depends only on one horizontal coordinate), under the further assumptions of small-height slide, which permits the recourse to linear theory, and of flat ocean floor, a solution for arbitrary body shape and velocity is deduced by applying the Duhamel theorem. It is also shown that this theorem can be advantageously used to obtain a general solution in case of a non-flat ocean floor, when the sea bottom follows a special power law, that can be adapted to study reasonable bottom profiles. The characteristics of the excited tsunamis are then evaluated by computing solutions in numerous examples, with special focus on wave pattern and wave evolution. The energy of the wave system is shown to depend on time: it grows expectedly in the initial phase of tsunami generation, when the moving body transfers energy to the water, but it may also diminish later, implying that a certain amount of energy may pass back from water waves to the slide.     

喀拉喀托火山喷发过程与巽他海啸成因

介绍了2018年12月22日发生的巽他海峡喀拉喀托火山喷发的过程及其火山监测情况,并提取了欧空局哨兵1A遥感卫星在火山喷发前后的遥感影像,通过遥感影像的对比分析获得了火山锥体的坍塌范围;使用3DAnalyst软件模块对坍塌部分的DEM影像进行分析,计算出海平面以上火山锥体的坍塌体积约为54000000m3,海面以下崩塌锥体体积更加巨大,崩塌导致火山周围水体发生激荡形成海浪,海浪相干传播至周边海岸附近引发了巽他海峡海啸;海啸灾害主要发生在印尼万丹省西冷县和板底兰县西部沿海,与根据火山喷发引发海啸的传播路径推测的受灾地区基本一致。     

Landslide and tsunami hazard at Yate volcano,Chile as an example of edifice destruction on strike-slip fault zones

The edifice of Yate volcano, a dissected stratocone in the Andean Southern Volcanic Zone, has experienced multiple summit collapses throughout postglacial time restricted to sectors NE and SW of the summit. The largest such historic event occurred on 19th February 1965 when ～6.1–10?×?10⁶ m³ of rock and ice detached from 2,000-m elevation to the SW of the summit and transformed into a debris flow. In the upper part of the flow path, velocities are estimated to have reached 40 m s^?1. After travelling 7,500 m and descending 1,490 m, the flow entered an intermontane lake, Lago Cabrera. A wavemaker of estimated volume 9?±?3?×?10⁶ m³ generated a tsunami with an estimated amplitude of 25 m and a run-up of ～60 m at the west end of the lake where a settlement disappeared with the loss of 27 lives. The landslide followed 15 days of unusually heavy summer rain, which may have caused failure by increasing pore water pressure in rock mechanically weathered through glacial action. The preferential collapse directions at Yate result from the volcano’s construction on the dextral strike-slip Liquiñe-Ofqui fault zone. Movement on the fault during the lifetime of the volcano is thought to have generated internal instabilities in the observed failure orientations, at ～10° to the fault zone in the Riedel shear direction. This mechanically weakened rock may have led to preferentially orientated glacial valleys, generating a feedback mechanism with collapse followed by rapid glacial erosion, accelerating the rate of incision into the edifice through repeated landslides. Debris flows with magnitudes similar to the 1965 event are likely to recur at Yate, with repeat times of the order of 10² years. With a warming climate, increased glacial meltwater due to snowline retreat and increasing rain, at the expense of snow, may accelerate rates of edifice collapse, with implications for landslide hazard and risk at glaciated volcanoes, in particular those in strike-slip tectonic settings where orientated structural instabilities may exist.     

Results of repeated leveling surveys at Newberry Volcano, Oregon, and near Lassen Peak Volcano, California

 Personnel from the U.S. Geological Survey's Cascades Volcano Observatory conducted first-order, class-II leveling surveys near Lassen Peak, California, in 1991 and at Newberry Volcano, Oregon, in 1985, 1986, and 1994. Near Lassen Peak no significant vertical displacements had occurred along either of two traverses, 33 and 44 km long, since second-order surveys in 1932 and 1934. At Newberry, however, the 1994 survey suggests that the volcano's summit area had risen as much as 97±22 mm with respect to a third-order survey in 1931. The 1931 and 1994 surveys measured a 37-km-long, east–west traverse across the entire volcano. The 1985 and 1986 surveys, on the other hand, measured only a 9-km-long traverse across the summit caldera with only one benchmark in common with the 1931 survey. Comparison of the 1985, 1986, and 1994 surveys revealed no significant differential displacements inside the caldera. A possible mechanism for uplift during 1931–1994 is injection of approximately 0.06 km³ of magma at a depth of approximately 10 km beneath the volcano's summit. The average magma supply rate of approximately 1×10^–3 km³/year would be generally consistent with the volcano's growth rate averaged over its 600,000-year history (0.7–1.7×10^–3 km³/year). Received: 10 September 1998 / Accepted: 4 December 1998     

Stochastic coupled simulation of strong motion and tsunami for the 2011 Tohoku,Japan earthquake

This study conducts coupled simulation of strong motion and tsunami using stochastically generated earthquake source models. It is focused upon the 2011 Tohoku, Japan earthquake. The ground motion time-histories are simulated using the multiple-event stochastic finite-fault method, which takes into account multiple local rupture processes in strong motion generation areas. For tsunami simulation, multiple realizations of wave profiles are generated by evaluating nonlinear shallow water equations with run-up. Key objectives of this research are: (i) to investigate the sensitivity of strong motion and tsunami hazard parameters to asperities and strong motion generation areas, and (ii) to quantify the spatial variability and dependency of strong motion and tsunami predictions due to common earthquake sources. The investigations provide valuable insights in understanding the temporal and spatial impact of cascading earthquake hazards. Importantly, the study also develops an integrated strong motion and tsunami simulator, which is capable of capturing earthquake source uncertainty. Such an advanced numerical tool is necessary for assessing the performance of buildings and infrastructure that are subjected to cascading earthquake–tsunami hazards.