Field survey of the 1993 Hokkaido Nansei-Oki earthquake tsunami

Runup data in Hokkaido and in three prefectures in the Tohoku District are described with a few witnessed arrival times and with comments of tide records. The highest runup of 31.7 m was found at the bottom of a narrow valley on the west coast of Okushiri Island. In order to explain high runups of 20 m at Hamatsumae in the sheltered area, roles of edge waves, refraction of the Okushiri Spur and tsunami generation by causes other than the major fault motion should be understood. An early arrival of the tsunami on the west coast of Hokkaido suggests another tsunami generation mechanism in addition to the major fault motion.     

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.     

香港海啸监测及警报系统的发展

地震监测、海啸数值模拟和海平面监测是监测和预报海啸的主要工具。为了有效监测南海北部可能发生的地震海啸,香港天文台（HKO）正在香港筹建一个宽频地震站,同时通过太平洋海啸警报及减灾系统（PTWS）的框架取得美国加州综合地震网（CISN）显示系统的实时地震信息,并通过世界气象组织（WM0）的全球通信系统（GTS）接收南海和西北太平洋的验潮站和海啸浮标数据以监测海面的波动情况。香港天文台通过联合国教科文组织（UNESCO）政府间海洋学委员会（IOC）取得海啸漫滩模式交换计划（TIME）下的海啸数值模式,把香港本地的高分辨率水深和地形数据融合在模式之内,并利用这个模式计算南海多处地区在不同地震情景下的海啸传播,为海啸预报提供重要的参考数据。     

南海地震与海啸

地震海啸的形成要具备3个条件:一是有深海盆地,可以容纳巨量海水;二是海底地形隆起与拗陷反差强烈;三是存在倾滑型活断层,可发生6级以上倾滑型的地震。查南海及其周边地形地貌,北西南三面都有宽阔的大陆架,中部又是平坦的深海平原,都不具备发生地震海啸的条件,惟独东侧马尼拉海沟才具备产生地震海啸的条件。南海地壳属于大洋型地壳与大陆型地壳之间的过度类型。其断裂构造非常发育,不同地段具有明显差异。北部为拉张型,南部为挤压型,西部为剪切型,东部为俯冲型,中部是扩张型。按断裂展布方向可分为NE向、NW向、EW向、SN向4组;按断裂切割深度,可分为岩石圈断裂、地壳断裂、基底断裂和盖层断裂。这些断裂多数为活动断裂,而东缘俯冲型断裂又是发震断裂。从地震分布、震源机制解分析,南海北、西、南以及中部都不具备引发地震海啸的条件,只有台南—菲律宾地震带东西两侧的贝尼奥夫带发生的倾滑型或具倾滑分量的走滑型6级以上地震,才有可能引发海啸,并可能对南海及我国东南沿海诸省以及港澳地区产生影响。     

利用PETSc和FEPG编制海啸数值模拟程序的研究

海啸的数值模拟是海啸研究的一个重要领域, 它对于帮助理解海啸的基本物理特性和预防减灾具有重要意义。 海啸数值模拟程序的编制是一项繁杂的工作, 该文介绍了利用两种功能强大的通用软件(PETSc和FEPG)来进行海啸数值模拟程序编写的研究。 PETSc和FEPG采用有限差分、 有限元、 有限体积等多种离散方法, 可以对数值问题给出稳定的求解。 该文采用了直角坐标系下的非线性浅水波方程作为海啸波传播的控制方程进行离散求解, 并将其结果与TUNAMI N1模型进行了比较, 表明这两种方法方便而且有效。     

Demarcation and Scaling of Long-term Seismogenesis

A space-time envelope of minor seismicity related to major shallow earthquakes is identified from observations of the long-term Precursory Scale Increase () phenomenon, which quantifies the three-stage faulting model of seismogenesis. The envelope, which includes the source area of the major earthquake, is here demarcated for 47 earthquakes from four regions, with tectonic regimes ranging from subduction to continental collision and continental transform. The earthquakes range in magnitude from 5.8 to 8.2, and include the 24 most recent mainshocks of magnitude 6.4 and larger in the San Andreas system of California, the Hellenic Arc region of Greece, and the New Zealand region, together with the six most recent mainshocks of magnitude 7.4 and larger in the Pacific Arc region of Japan. Also included are the destructive earthquakes that occurred at Kobe, Japan (1995, magnitude 7.2), Izmit, Turkey (1999, magnitude 7.4), and W.Tottori, Japan (2000, magnitude 7.3). The space (A _P ) in the space-time envelope is optimised with respect to the scale increase, while the time (T _P ) is the interval between the onset of the scale increase and the occurrence of the earthquake. A strong correlation is found between the envelope A _P T _P and the magnitude of the earthquake; hence the conclusion that the set of precursory earthquakes contained in the envelope is intrinsic to the seismogenic process. Yet A _P and T _P are correlated only weakly with each other, suggesting that A _P is affected by differences in statical conditions, such as geological structure and lithology, and T _P by differences in dynamical conditions, such as plate velocity. Among other scaling relations, predictive regressions are found between, on the one hand, the magnitude level of the precursory seismicity, and on the other hand, both T _P and the major earthquake magnitude. Hence the method, as here applied to retrospective analysis, is potentially adaptable to long-range forecasting of the place, time and magnitude of major earthquakes.     

Identification of 1771 Meiwa Tsunami deposits using a combination of radiocarbon dating and oxygen isotope microprofiling of emerged massive Porites boulders

The Meiwa Tsunami, one of the largest tsunamis recorded in historical documents in Japan struck Ishigaki Island and neighboring islands of the Ryukyu Arc in April 1771 AD, killing more than 12 000 people. An enormous number of massive Porites coral boulders are scattered on the shore and in the reef moat of eastern Ishigaki Island. Although these boulders likely were cast ashore by the Meiwa Tsunami, a detailed examination has not yet been conducted. When the marine reservoir effect is taken into account, one of mode values of calibrated radiocarbon dates possibly corresponds to the time of the 1771 event. However, the range of calibrated radiocarbon ages indicates that the transport of the boulders cannot be ascribed solely to the tsunami. Oxygen isotope microprofiling, which indicates sea-surface temperature variation, was therefore conducted to further investigate the mechanism of transport. The results suggest that the skeletal growth of most coral colonies was interrupted in summer or autumn; hence, tropical storms and typhoons are also very likely to be agents of transport. Thus, by combining radiocarbon dating with oxygen isotope microprofiling to investigate Porites coral boulders, it is possible to separate paleotsunami boulders from those transported by storm events as far as tsunamis occurred during the non-storm season.     

Review of the 2018 Palu Tsunami:Disaster Behavior and Reflections on Disaster Management

The 2018 Palu M_W7.5 earthquake and tsunami attracted geophysicists'' attention for its strike slip focal mechanism and magnitude. We inspected the details of this disaster and discussed its particularity and possible causations. The submarine landslide and special terrain conditions could have contributed to the unusual size of the tsunami. The early warning system and the post disaster response is also reviewed. Efficient social warnings and broadcast systems along with good maintenance is essential. We also found that enhancing publics scientific literacy is the most important way to reduce disaster damage and casualties. Moreover, social conditions and rebuilding difficulties post tsunami are related as reference resources for future disaster management strategies.     

Numerical Simulations of Tsunami Waves and Currents for Southern Vancouver Island from a Cascadia Megathrust Earthquake

The 1700 great Cascadia earthquake (M = 9) generated widespread tsunami waves that affected the entire Pacific Ocean and caused damage as distant as Japan. Similar catastrophic waves may be generated by a future Cascadia megathrust earthquake. We use three rupture scenarios for this earthquake in numerical experiments to study propagation of tsunami waves off the west coast of North America and to predict tsunami heights and currents in several bays and harbours on southern Vancouver Island, British Columbia, including Ucluelet, located on the west coast of the island, and Victoria and Esquimalt harbours inside Juan de Fuca Strait. The earthquake scenarios are: an 1100-km long rupture over the entire length of the subduction zone and separate ruptures of its northern or southern segments. As expected, the southern earthquake scenario has a limited effect over most of the Vancouver Island coast, with waves in the harbours not exceeding 1 m. The other two scenarios produce large tsunami waves, higher than 16 m at one location near Ucluelet and over 4 m inside Esquimalt and Victoria harbours, and very strong currents that reach 17 m/s in narrow channels and near headlands. Because the assumed rupture scenarios are based on a previous earthquake, direct use of the model results to estimate the effect of a future earthquake requires appropriate qualification.     

Power-law Distributions of Offspring and Generation Numbers in Branching Models of Earthquake Triggering

We consider a general stochastic branching process,which is relevant to earthquakes as well as to many other systems, and we study the distributions of the total number of offsprings (direct and indirect aftershocks in seismicity) and of the total number of generations before extinction. We apply our results to a branching model of triggered seismicity, the ETAS (epidemic-type aftershock sequence) model. The ETAS model assumes that each earthquake can trigger other earthquakes (aftershocks). An aftershock sequence results in this model from the cascade of aftershocks of each past earthquake. Due to the large fluctuations of the number of aftershocks triggered directly by any earthquake (fertility), there is a large variability of the total number of aftershocks from one sequence to another, for the same mainshock magnitude. We study the regime in which the distribution of fertilities is characterized by a power law ~1/¹⁺. For earthquakes we expect such a power-distribution of fertilities with =b/ based on the Gutenberg-Richter magnitude distribution ~ 10^–bm and on the increase ~ 10^–m of the number of aftershocks with the mainshock magnitude m. We derive the asymptotic distributions p_r(r) and p_g(g) of the total number r of offsprings and of the total number g of generations until extinction following a mainshock. In the regime < 2 for which the distribution of fertilities has an infinite variance, we find This should be compared with the distributions obtained for standard branching processes with finite variance. These predictions are checked by numerical simulations. Our results apply directly to the ETAS model whose preferred values =0.8–1 and b=1 puts it in the regime where the distribution of fertilities has an infinite variance. More generally, our results apply to any stochastic branching process with a power-law distribution of offsprings per mother     

Deduction of groundwater flow regime in a basaltic aquifer using geochemical and isotopic data: The Golan Heights, Israel case study

Groundwater flow-paths through shallow-perch and deep-regional basaltic aquifers at the Golan Heights, Israel, are reconstructed by using groundwater chemical and isotopic compositions. Groundwater chemical composition, which changes gradually along flow-paths due to mineral dissolution and water–rock interaction, is used to distinguish between shallow-perched and deep-regional aquifers. Groundwater replenishment areas of several springs are identified based on the regional depletion in rainwater δ¹⁸O values as a function of elevation (−0.25‰ per 100 m). Tritium concentrations assist in distinguishing between pre-bomb and post-bomb recharged rainwater.

It was found that waters emerging through the larger springs are lower in δ¹⁸O than surrounding meteoric water and poor in tritium; thus, they are inferred to originate in high-elevation regions up to 20 km away from their discharge points and at least several decades ago. These results verify the numerically simulated groundwater flow field proposed in a previous study, which considered the geological configuration, water mass balance and hydraulic head spatial distribution.     

Scaling of response spectrum and duration for aftershocks

Aftershocks have the potential to cause collapse of a structure that has been already damaged by the preceding main shock. Seismic safety of a structure should therefore be ascertained through a damage analysis using the anticipated main shock and few larger-aftershock motions. Simulation of aftershock motions needs characterization of the seismic hazard due to aftershocks, and therefore it will be useful to develop a conditional scaling model that can predict the response spectrum of an anticipated aftershock motion consistent with the design spectrum of the main shock motion anticipated at the same station. In this study an attempt is made to develop a conditional scaling model for the pseudo spectral velocity spectrum via linear regression analysis on the aftershock and main shock recordings for the 1999 Chi–Chi earthquake. It is shown that it may be possible to obtain a simpler and approximate version of the conditional model from an unconditional model. Damage-causing potential of a ground motion also depends on its strong motion duration (SMD) and therefore a conditional scaling model is developed for SMD of the aftershock motion in several narrow frequency-bands. The model is developed for the larger-aftershock motions and it is shown that a reasonable replacement of such a model may be obtainable directly from an unconditional model. Finally, a simple weighted averaging scheme is proposed to obtain the composite SMD from the SMDs for different frequency bands by using the pseudo spectral acceleration spectrum of the motion.     

Features and formation processes of multiple deposition layers from the 2004 Indian Ocean Tsunami at Ban Nam Kem,southern Thailand

Multiple‐layered tsunami deposits have been frequently reported from coastal stratigraphic sequences, but the formation processes of these layers remain uncertain. A terrestrial sandy deposit formed by the 2004 Indian Ocean Tsunami was investigated at Ban Nam Kem, southern Thailand. Four internal layers induced by two tsunami waves were identified in the tsunami deposit. Sedimentary structures indicated that two units were formed by run‐up currents caused by the tsunami and the other two units were deposited by the backwash flows. Graded bedding was common in the layers, but inverse grading was observed at limited intervals on the surveyed transects. The characteristics of the multiple‐layered tsunami deposit vary remarkably over a very short distance (<1 m) in response to the local topography. Remarkable asymmetries in thickness and grain‐size distribution are recognized between the run‐up and backwash flow deposits. On the basis of the interpretation of sedimentary structures, the formation process of the multiple‐layered tsunami deposit observed in this study can be explained in a schematic model as the modification of the ideal tsunami sequence by local erosion and the asymmetric hydraulic properties of tsunami waves, such as the maximum shear velocity and the heterogeneity of the flow velocity field.     

Scaling and regionalization of flood flows in British Columbia,Canada

A regionalization of flood data in British Columbia reveals a common scaling with drainage area over the range 0·5×10²<A_d<10⁴ km². This scaling is not a function of flood return period, which implies that simple scaling—consistent with a snowmelt‐dominated flow regime—applies to the province. The observed scale relation takes the form , similar to values reported in previous studies. The scaling relation identified was used to define the regional pattern of hydroclimatic variability for flood flows in British Columbia after discounting the effect of drainage area. The pattern was determined by kriging a scale‐independent runoff factor k for the mean annual flood, 5 year flood and 20 year flood. The analysis permits quantification of uncertainty of the estimates, which can be used in conjunction with the mapped k‐fields to calculate a mean and range for floods with the identified return period for ungauged basins. Owing to the sparsity of data, the precision is relatively poor. The standard error is generally less than 75% of the estimate in the southern half of the province, whereas in the northern half it is often between 75 and 100%. Examination of the relative increase in flood magnitude with increasing return period reveals spatially consistent but statistically insignificant differences. Flood magnitude tends to increase more rapidly in the western regions, where rain events may contribute to flood generation. The relative increase in flood magnitude with return period is consistently lower in the eastern mountain ranges, where snowmelt dominates the flood flow regime. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

Scaling of Fracture and Faulting of Ice on Earth

The scaling properties of fracture and faulting of ice on Earth are reviewed.Numerous evidences for the scaling of fracture and faulting of ice are given,including self-affine fracture surfaces, fractal fracture networks at small(laboratory) and large (geophysical) scales, power law distributions of fracturelengths or of fragment sizes within fault gouges. These scaling laws are discussedin terms of the underlying mechanics. Scaling of the observables associated withfracture and faulting argues for the scale invariance of the fracture and faultingprocesses and indicates that small scales cannot be arbitrarily disconnected fromlarge scales. Consequently, quantitative links between scales cannot be performedthrough classical homogenization procedures. Scaling can also induce scale effectson different mechanical parameters such as fracture energy, strength or stiffness.Although scaling is ubiquitous for the fracture of ice on Earth, important exceptionsexist such as the nucleation of microcracks or the crevassing of glaciers. Theseexceptions are stressed and discussed.     

基于小波变换和支持向量机的中国大陆强震预测

将小波变换和支持向量机用于中国大陆年度最大地震震级预测。 先用小波变换把中国大陆年度最大地震序列分解成几个不同尺度水平(频率)的子序列, 然后使用支持向量机对分解后的子序列分别进行预测, 最后通过重构几个子序列的支持向量机预测结果得到最终预测结果, 预测次年中国大陆最大地震震级。 与支持向量机和神经网络方法对比, 结果表明小波变换和支持向量机相结合方法具有更高的预测精度, 预测效果很好, 说明此方法可用于地震时间序列预测。     

Spatial Scaling of Effective Modulus and Correlation of Deformation Near the Critical Point of Fracturing

Many observations point to the lithosphere being metastable and close to a critical mechanical point. Exercises in modelling deformation, past or present, across subsurface reservoirs need to take account of this criticality in an efficient way. Using a renormalization technique, the spatial scaling of effective elastic modulus is derived for 2-D and 3-D bodies close to the critical point of through-going fracturing. The resulting exponent, d_μ, of spatial scaling of effective modulus with size, , takes the values ~ −2.5 and −4.2 in two- and three-dimensional space, respectively. The exponents are compatible with those for scaling of effective modulus with fracture density near the percolation threshold determined by other workers from numerical experiments; the high absolute values are also approximately consistent with empirical data from a) fluctuations in depth of a seismic surface; b) `1/k' scaling of heterogeneities observed in one-dimensional well-log samples; c) spatial correlation of slip displacements induced by water injection. The effective modulus scaling modifies the spatial correlation of components of displacement or strain for a domain close to the critical point of fracturing. This correlation function has been used to geostatistically interpolate components of the strain tensor across subsurface reservoirs with the prime purpose of predicting fracture densities between drilled wells. Simulations of strain distributions appear realistic and can be conditioned to surface depths and observations at wells of fracture-related information such as densities and orientations, welltest permeabilities, changes in well-test permeabilities, etc.     

Scaling properties and spatial organization of snow depth fields in sub‐alpine forest and alpine tundra

This paper reports on a study analysing the spatial distribution functions, the correlation structures, and the power spectral densities of high‐resolution LIDAR snow depths (～1 m) in two adjacent 500 m × 500 m areas in the Colorado Rocky Mountains, one a sub‐alpine forest the other an alpine tundra. It is shown how and why differences in the controlling physical processes induced by variations in vegetation cover and wind patterns lead to the observed differences in spatial organization between the snow depth fields of these environments. In the sub‐alpine forest area, the mean of snow depth increases with elevation, while its standard deviation remains uniform. In the tundra subarea, the mean of snow depth decreases with elevation, while its standard deviation varies over a wide range. The two‐dimensional correlations of snow depth in the forested area indicate little spatial memory and isotropic conditions, while in the tundra they indicate a marked directional bias that is consistent with the predominant wind directions and the location of topographic ridges and depressions. The power spectral densities exhibit a power law behaviour in two frequency intervals separated by a break located at a scale of around 12 m in the forested subarea, and 65 m in the tundra subarea. The spectral exponents obtained indicate that the snow depth fields are highly variable over scales larger than the scale break, while highly correlated below. Based on the observations and on synthetic snow depth fields generated with one‐ and two‐dimensional spectral techniques, it is shown that the scale at which the break occurs increases with the separation distance between snow depth maxima. In addition, the breaks in the forested area coincide with those of the corresponding vegetation height field, while in the tundra subarea they are displaced towards larger scales than those observed in the corresponding vegetation height field. Copyright © 2009 John Wiley & Sons, Ltd.