RUSSIAN JOURNAL OF EARTH SCIENCES VOL. 8, ES6003, doi:10.2205/2006ES000212, 2006
Numerical modeling of the generation of long waves by a dynamic seismic source and their propagation in the Black SeaR. Kh. Mazova and E. A. Tresvyatskaya Nizhniy Novgorod Technical University, Nizhniy Novgorod, RussiaContents
Abstract[1] Numerical modeling of generation and propagation of tsunami waves caused by underwater earthquakes in the Black Sea basin is carried out using the model of dynamic seismic source. A number of historical tsunamis are considered as well as possible catastrophic tsunami with large magnitudes from the sources whose location is close to the historical events considered here. The process of tsunami source formation and propagation of tsunami waves in the entire basin of the sea up to 10-meter isobath was considered. Kinematic vertical motion of key-blocks, into which the source of the earthquake is divided, was used to model the seismic source. The obtained results were compared with the data of pressure gauges and results of other authors. The distribution of maximal heights of tsunami waves along the entire Russian coast of the Black Sea was obtained for each calculation. Pressure gauge records were plotted for a number of points at the eastern coast located close to eight Russian cities such as Anapa, where the probability of earthquakes and tsunami is high in the nearest future. Introduction[2] Numerical modeling of tsunami at the coast of the Black Sea including the estimate of tsunami hazard and seismic hazard of the Russian coast (see, for example, [Garagash and Lobkovsky, 2000; Solovieva and Kuzin, 2005]) is a pressing objective of research in the last decades. High industrial potential of the region (large ports, terminals of gas and oil pipelines), and the largest recreation zone that exists here determine the importance of this problem. The urgency of such calculations is related, in particular, to the problem of operation of the marine part of gas pipeline Russia-Turkey (the "Blue Flow" project) connecting the territories of these countries over the bottom of the Black Sea, which should operate in the conditions of increased seismic and landslide hazard of the Russian and Turkish slopes of the Black Sea [Garagash and Lobkovsky, 2000; Solovieva and Kuzin, 2005]. [3] Although seismicity of the Caucasus and Crimea coasts is considered moderate (see, for example, [Solovieva and Kuzin, 2005]), according to the chart of the maximal vibrations over the territory of North Caucasus, the coast of the Black Sea from Anapa to Sochi falls within the zone of seventh grade of intensity, which approximately corresponds to magnitude 6 and horizontal velocities from 0.1 to 0.2 m/s [Lobkovsky and Baranov, 1984; Solovieva and Kuzin, 2005]. Even more, detailed seismological observations using bottom stations of high sensitivity in the Jubga region where the terminal of the gas pipeline Russia-Turkey is located, revealed a very high level of seismic activity at the level of micro- and weak earthquakes. Such seismicity is caused by the regional tectonics since the Black Sea is located in the boundary region of the interaction between the Arabian and Eurasian tectonic plates, which results in the fact that in the northeastern part of the Black Sea Depression, significant stresses are formed in the Earth's crust, which later relax in the form of earthquakes in the zones of shelf and continental slope as well as in the fractures of the costal zone [Solovieva and Kuzin, 2005]. [4] In addition to the seismic displacements of the sea bottom, development of large landslide displacements of large blocks of sedimentary formations are possible over the underwater slopes in their upper parts induced by seismicity of the bottom [Garagash and Lobkovsky, 2000; Lobkovsky and Baranov, 1984]. Such landslides can also be the sources of tsunamis (cf. with [Garder et al., 1993; Gusyakov V.K., 1978; Novikova and Ostrovsky, 1978]), which were mentioned in historical tsunami data in this region (between Sochi and Tuapse in 1909 [Solovieva and Kuzin, 2005; Zaitsev et al., 2002]). [5] As was mentioned in [Solovieva and Kuzin, 2005], the peculiarity of tsunami propagation in the basin of the Black Sea is trapping of wave energy in the zone of the shelf and continental slope, which results in the fact that for the sources in the western part of the Black Sea the wave energy is almost not transferred to the eastern part of the sea, and on the contrary, tsunamis generated in the eastern part of the sea are strongly attenuated when the waves approach its western coast. [6] The Anapa region is one of the seismic active regions at the Russian coast of the Black Sea. Underwater earthquakes repeatedly occurred and tsunamis were recorded here. The recurrence period of tsunami is approximately 40 years [Solovieva and Kuzin, 2005]. The last of the recorded earthquakes and tsunamis in this region occurred in July 1966, thus numerical modeling of the possible tsunami in the nearest time at the coast of the Black Sea based on the new model of seismic source and characteristic of tsunami at the coast in the Anapa region, where the well known recreation center for children is located, seems were actual.
Formulation of the Problem
[9] Here,
u(x, t) is horizontal velocity component;
h (x, y, t) is perturbation of
the free surface;
H is maximal depth of the basin;
B(x, y, t) is variation in the bottom
of the basin;
H [10] Dynamic seismic source in the form of key-blocks, which can rise (or descend) to a given distance with different velocity, was used for tsunami wave generation [Lobkovsky and Baranov, 1984]. At the initial time moment, parameters of dynamic seismic source (coordinates, velocity of the vertical motion of key-block) are specified in the calculation area. The boundaries are located either in the land regions or the water surface. Numerical Modeling of Generation and Propagation of Tsunami WavesEquations of Motion[11] We shall consider generation of tsunami waves in the shallow water approximation (see, e.g. [Gusyakov V.K., 1978; Pelinovsky, 1982]. We consider nonlinear equations (1), where function B(x, y, t) describes the vertical motion of key-blocks in the seismic source. The law of the motion of key-blocks changes depending on the version of the calculation. This dynamic model allows us to model various variations of the sea bottom without significant complication of the algorithm.
where x, y are spatial coordinates, t is time, u(x, t), v(y, t) are horizontal components of velocity, h(x, y, t) is perturbation of the free surface with respect to its level at rest, H is maximal depth of the basin, B is variation in the sea bottom (the account for the characteristics of dynamic seismic source), g is acceleration due to gravity. Initial and Boundary Conditions[12] Boundary conditions in the problems of tsunami wave propagation differ depending on the peculiarity of the problems considered in the study. For example, if propagation of tsunami waves is considered in the basin of the sea, it is possible to specify the condition of complete reflection (specifying zero normal components of velocity un, vn, and derivative hn ) or to specify "free" boundary conditions. In the latter case, the wave passes through the boundary without distortions. However, a more complete description of the interaction between the wave and coast is needed in the solution of the problems of tsunami wave propagation in the shelf zone [Pelinovsky, 1982; Pelinovsky and Mazova, 1992].[13] We consider that before the generation, the fluid that fills the region with initial form of the bottom is at rest, i.e. the velocity and perturbation of the free surface are zero: Boundary conditions (full reflection) at 10-meter isobath are written as:
The conditions of free passing wave are written as
where n+1 denotes the time layer following layer n, and c = (g (H + h))1/2. [14] This model makes possible to correlate the wave with the generating source taking into account the characteristics of displacements and velocities at each point of variable sea bottom. The method allows us to combine the processes of numerical modeling of the generation and propagation of tsunami waves taking into account the effects of real basin bathymetry residual variations in bottom topography and superposition of subsequent perturbations on the wave propagating in the ocean. The seismic characteristics of the earthquake source influence the wave only by means of function B(x, y, t). [15] The numerical scheme described in and adjusted for the solution of the problems with dynamic source [Garagash et al., 2003] was used to calculate the generation and propagation of tsunami waves.
[16] Bathymetry of the Black Sea with a resolution of 1.5 km was used in the modeling. The
number of nodes in the numerical scheme is equal to 505
Results of Numerical ModelingCalculation of Historical Tsunami[17] Solovieva and Kuzin, [2005] considered the main earthquakes and tsunamis at the Russian coast of the Black Sea. Events were selected from these data, for which the dates of earthquakes, coordinates of the earthquake epicenters, and run-up data were known. This makes possible adequate modeling of the past events for analyzing possible dynamic displacements within the seismic source. The results of the calculations were checked with the data of observations on run-up at the coast for each event and with the results of calculations performed by other authors [Solovieva and Kuzin, 2005; Zaitsev et al., 2002]. Thus, we selected three seismic events on 26 June 1927; 12 September 1927; and 12 July 1966 with coordinates of the sources 34.4oE, 44.4oN; 34.5oE, 44.5oN; and 37.3oE, 44.7oN, respectively (see Table 1). [18] We obtained approximate sizes of tsunami source and altitude, at which water could be displaced in the source using Iida relation (see, for example [Pelinovsky, 1982]) for estimating the parameters of tsunami source at given magnitude of the earthquake M.
where R is radius of the source, and h is displacement of water over the seismic source. It is noteworthy that owing to incompressibility of the fluid, the displacement at the surface would repeat exactly the bottom displacement in the seismic source. Thus, formation of the tsunami source would occur during the time of dynamic displacements in the seismic source, in our case, during t seconds of the rising of the block. The size of the tsunami source R allows us to specify the approximate size of the block (or blocks) comprising the seismic source, while the displacement of the wave surface over the seismic source would allow us to calculate the altitude of the elevation of the block (or blocks) in the seismic source. Calculation 1.[19] Let us consider earthquake on 26 June 1927 (M = 5.8) with coordinates of the epicenter: 34.4oE, 44.4oN. Using relations (4) we get that the radius of the source is R = 6.3 km, and the altitude of the displacement of water surface is h = 0.12 m.[20] Thus, we assume that approximate parameters of the seismic source are: D = 12.6 km is the size of the key-block; B = 0.2 m is the altitude of the key-block.
[22] Patterns of the locations of wave fronts for four time moments are shown in Figures 4a, 4b, 4c, and 4d: in 200 s after the beginning of tsunami wave generation by the seismic source, after 1200 s, after 2200 s, and after 8800 s. It is clearly seen that after 20 min (1200 s) the wave front reached Yalta and Feodosiya, and approximately after 40 min, the wave front reached the Turkish coast of the Black Sea. Almost in two and a half hours, the wave reached all points at the Black Sea coast.
[25] It is clearly seen from the pressure gauge series that at each of the three points, tsunami started with sea level rise. For example, in Eupatoria, the first sea level rise was 2 cm, and only the fifth wave caused a rise of 9 cm. Fluctuations of the sea level continued for a long time and ranged from - 8 to 9 cm. In Yalta, the height of the first wave was 8 cm, while the second was already equal to 15 cm. After this, the fluctuations of the sea level decreased and continued for a long time ranging from - 7 to 8 cm. In Feodosiya, the height of the first wave was 4 cm, and for a long time the fluctuations of the level ranged from - 10 to 8 cm. It is clearly seen that the distribution of maximal wave heights along the Russian coast of the Black Sea has a small scatter (Figures 5 and 6). Calculation 2.[26] The earthquake on 12 September 1927 (M = 6) corresponds to the following coordinates of the epicenter: 34.5oE, 44.4oN. The coordinates of the source are close to the case considered above, however, the recorded magnitude of the earthquake was higher (M = 6). The parameters of the tsunami source calculated from relations (4) are the following: radius of the tsunami source is R = 7.9 km, while the elevation of water at the source is h = 0.18 m.[27] Thus, we assume approximate parameters of the seismic source: D = 15.8 km is the size of the key-block; B = 0.3 m is the altitude of the key-block.
[29] Locations of tsunami wave fronts at four time moments: after 220 s, after 2200 s, after 3550 s, and after 7100 s are shown in Figure 8.
[31] Pressure gauge series calculated at three points of the Black Sea coast: Yalta, Novorossiysk, and Tuapse are shown in Figure 10. The choice of the towns is related to the possibility of testing the results of the calculations on the basis of the available data of measurements and calculations. Calculation 3.[32] The earthquake on 12 July 1966 (M = 5.8) corresponds to the following coordinates of the epicenter: 37.3oE, 44.7oN. The parameters of the tsunami source calculated from relations (4) are the following: radius of the tsunami source is R = 6.3 km, while the elevation of water at the source is h = 0.12 m.[33] Thus, we assume approximate parameters of the seismic source: D = 12.6 km is the size of the key-block; B = 0.22 m is the altitude of the key-block.
[35] Locations of tsunami wave fronts and spatial distribution of sea level at four time moments after wave propagation from this source are shown in Figure 12: after 310 s, after 1200 s, after 3500 s, and after 9100 s. It is clearly seen that already in 20 min after the earthquake, the waves would reach Anapa, Novorossiysk, and Gelendzhik. In half an hour, the waves would reach Tuapse, and in 40 min they would reach Sochi. Less than one hour would be necessary for the waves to reach the Turkish coast of the Black Sea (see Figures 12 a-d).
[37] One can see that localization of the earthquake source near Anapa yields in a peak of maximal wave heights at the coast in the region of Anapa, Novorossiysk, and Gelendzhik. The greatest heights reach 42 cm in Gelendzhik (see Figure 15a, and Table 1). [38] It is clearly seen that the results obtained by numerical modeling using the dynamic keyboard model of subduction [Zaitsev et al., 2002] agree well with the data of field observations and calculations of other authors almost at all points. However, to our opinion, a certain difference from the calculations made by other authors for the earthquakes on 26 June 1927 in Eupatoria and on 12 September 1927 in Yalta is related to a better physical basis of the earthquake source model. Calculation of Possible Catastrophic Tsunami at the Black Sea Coast Generated at Hypothetical Seismic Sources[39] Let us consider generation and propagation of a tsunami wave over the basin of the Black Sea for six sources with equal size, the same localization, and assume that in all six cases (scenarios) the magnitude of the possible catastrophic earthquake is M = 8.5. The motion of blocks in the source is considered with the same velocity 0.54 km h -1. The displacement is the same: up or down. The location of the source corresponds to the coordinates of the earthquake on 26 June 1927: 34.4oE-33.6oE; 44.3oN-43.5oN. From relations (4) we get the size of the tsunami source R = 141.2 km and height of water displacement at the tsunami source h = 2.5 m. Rising of block (blocks) in the seismic source occurs during 30 s, then the recalculated elevation of the block at the bottom would be equal to 4.5 m. We assume that the size of the seismic source is approximately D = 282 km.
[41] 1. The source consists of one block moving upwards (Figure 16a); [42] 2. The source consists of one block moving downwards (Figure 16b); [43] 3. The source consists of two blocks simultaneously moving upwards and downwards and the negative motion in the source is oriented to the sea side (Figures 16c, 16d); [44] 4. The source consists of two blocks simultaneously moving upwards and downwards and the negative motion in the source is oriented to the coast side (Figures 16e, 16f). Calculation 1.[45] Let us, first, consider the motion of the block as a whole (Figure 16a). We shall consider propagation of waves in the basin of the Black Sea. Distribution of maximal wave heights at 10-meter isobath for the entire Black Sea coast of Russia for this calculation is shown in Figure 15a.
[47] Analysis of pressure gauge series indicates that at all points the wave would reach the coast in 15-50 minutes, and only in Yalta, the first front would be observed almost instantaneously, in 3-5 m. after the beginning of the seismic process. The greatest values of run-up and run-down in Yalta would be equal almost to three meters, after this the amplitudes would decrease but however, for a long time the maxima of the approaching wave trains would remain within 1.7-1.2 m. In Feodosiya and Eupatoria, the wave heights would reach 1.5-1.8 m, while the maximal run-down values would reach 1.5 m. In Gelendzhik, Novorossiysk, Tuapse, and Sochi, the first maximal waves of run-up and run-down would be in the interval 0.5-1 m. It is worth noting that since we consider only local tsunamis and the source is located close to the coast, the first wave that reaches the coast (10-meter isobath), would be a run-up wave at all points of the coast. Calculation 2.[48] Let us consider the same formulation of the problem, but the block as a whole moves downwards (Figure 16b). Distribution of maximal wave heights at 10-meter isobath for the entire Black Sea coast of Russia for this calculation is shown in Figure 15b.
[50] It is interesting to note that it follows from the comparison of Figures 19a and 19b that the patterns of the distribution of wave heights for the maximal and minimal values calculated along the coast are overturned. It is seen well that maximal values of run-up and run-down are observed in the region left of Yalta. Calculation 3.[51] Let us consider a seismic source (Figure 16c) consisting of two blocks of equal size. The blocks move simultaneously upwards and downwards with equal speed of 0.54 km h -1. The calculated values of vertical displacement of blocks are equal to 4.5 m and - 4.5 m, respectively.[52] The distributions of maximal wave heights at 10-meter isobath along the entire Russian coast of the Black Sea for this calculation are shown in Figure 15c.
[54] It is clearly seen that the wave arrived to Eupatoria with a negative height, i.e. the water, first, flowed off the coast (60 cm in the vertical direction), and then returned with a height of 1.3 m. At all other points, tsunami started with run-up. The maximal run-up reaching 3.5 m is possible in Yalta. Calculation 4.[55] In this calculation we consider an "overturned" seismic source (relative to Calculation 3) oriented as shown in Figure 16d. Negative displacement is oriented in the region of Gelendzhik, Tuapse, and Sochi. The blocks move simultaneously upwards and downwards with equal speed of 0.54 km h -1. The calculated values of vertical displacement of blocks are equal to 4.5 m and - 4.5 m, respectively.
Calculation 5.[57] In this calculation, in the source of the earthquake of the same size as in the previous calculations the negative displacement is oriented into the open sea (Figure 16e). The blocks move simultaneously upwards and downwards with equal velocity of 0.54 km h -1. The displacement of the blocks is the same with the same velocity.[58] Distribution of maximal wave heights at 10-meter isobath is shown in Figure 15e.
Calculation 6.[60] In this calculation, the source is similar to the previous one; however, the negative displacement is oriented to the coast (Figure 16f). The blocks move simultaneously upwards and downwards with equal velocity of 0.54 km h -1. The displacement of the blocks is the same with the same velocity as in the previous calculation. Localization of the source and magnitude of the earthquake M = 8.5 is the same as in the previous calculations (see Figure 15f).
Conclusion[64] In this work, we carried out numerical modeling of the generation of long gravity waves in water by a dynamic seismic source and their propagation for the specific historical events in the Black Sea. It was shown that for a given magnitude of earthquake the source provides a good agreement with the data of pressure gauges recorded for the given event. The performed calculation of possible catastrophic tsunami at the Black Sea coast from hypothetical seismic sources showed that if the source is localized in the region with coordinates 34.4oE-33.6oE; 44.3oN-43.5oN, the regions of Yalta and Novorossiysk would be subjected to greatest hazard, while the regions of Sochi and Gelendzhik would rest in the calm zone. [65] The analysis demonstrated that for local tsunamis at the same earthquake magnitude and conservation of the source parameters and its localization the possible values of tsunami run-up and run-down at the coast depend significantly in the sequence of the motion of key-blocks in the seismic source. [66] It is clear that the performed analysis is only an example of using the dynamic keyboard model for calculating possible consequences after earthquakes with magnitude M = 8.5. Complete analysis of the possible catastrophic consequences of tsunami wave impact caused by strong earthquakes requires numerical modeling of the generation and propagation of tsunami waves at the location of the earthquake source in the most dangerous zones of fractures of tectonic plates, which will be studied in our future work. Acknowledgments[67] This work was supported by Russian Foundation for Basic Research no. 05-05-64685 and EU Project 502247 (COMSHELFRISKS). ReferencesGaragash, I. A., and L. I. Lobkovsky (2000), Geomechanical evaluation of slide processes and their monitoring at slopes of the Black Sea in connection with realization of the "Blue Flow", 6 Int. Sci.-Indust. Conf. "Modern Methods and Equipment for Oceanological Research", Nauka, Moscow. Garagash, I. A., L. I. Lobkovsky, O. R. Kozyrev, and R. Kh. Mazova (2003), Generation and Runup of Tsunami Wave at Moving of Submarine Slide, Oceanology, 43, 185. Garder, O. I., I. S. Dolina, E. N. Pelinovsky, A. A. Poplavsky, and V. E. Fridman (1993), Generation of tsunami waves by gravity lithodynamic processes, Issledovaniya Tsunami (in Russian), (5), 50. Gusyakov, V. K. (1978), Review of Works on the Problem of Initiation of Tsunami Wave, in: Methods of Calculation of Generation and Propagation of Tsunami (in Russian), p. 18, Nauka, Moscow. Lobkovsky, L. I., and B. V. Baranov (1984), Keyboard model of strong earthquakes in the island arches and active continental boundaries, Dokl. Acad. Sci. USSR, 275, (4), 843. Novikova, L. E., and L. A. Ostrovsky (1978), Excitation of tsunamy waves by movements of ocean bottom., in: Methods of Calculation of Generation and Propagation of Tsunami (in Russian), p. 88, Nauka, Moscow. Pelinovsky, E. N. (1982), Nonlinear dynamics of tsunami waves, 227 pp., IPF AN SSSR, Gorky. Pelinovsky, E. N., and R. Kh. Mazova (1992), Exact analytical solution of nonlinear Problem of tsunami wave run-up on slops with different profiles, Sci. Natural Hazards, 6, 227, doi:10.1007/BF00129510. [CrossRef] Solovieva, O. N., and I. P. Kuzin (2005), Seismicity and tsunami in the northeastern part of the Black Sea, Oceanology, 45, 826. Yanushauskas, A. I. (1978), Generation of waves by underwater sources with alternative intensity, in: Methods of Calculation of Generation and Propagation of Tsunami (in Russian), p. 100, Nauka, Moscow. Zaitsev, A. I., A. S. Kozelkov, A. A. Kurkin, E. N. Pelinovsky, T. G. Talipova, and A. S. Yalchiner (2002), Tsunami modeling in the Black Sea, Izv. AIN RF, Series: Applied Mathematics and Mechanics (in Russian), 3, 27. Received 21 November 2006; accepted 15 December 2006; published 24 January 2007. Keywords: tsunami danger evaluation, seismic danger, region tectonics, shelf zone and continental slope zone, dynamic seismic source. Index Terms: 3025 Marine Geology and Geophysics: Marine seismics; 3285 Mathematical Geophysics: Wave propagation; 4255 Oceanography: General: Numerical modeling; 4564 Oceanography: Physical: Tsunamis and storm surges; 7209 Seismology: Earthquake dynamics. ![]() Citation: 2006), Numerical modeling of the generation of long waves by a dynamic seismic source and their propagation in the Black Sea, Russ. J. Earth Sci., 8, ES6003, doi:10.2205/2006ES000212. (Copyright 2006 by the Russian Journal of Earth SciencesPowered by TeXWeb (Win32, v.2.0). |