Statistical analysis of time-dependent earthquake occurrence and its impact on hazard in the low seismicity region Lower Rhine Embayment

The time-dependence of earthquake occurrence is mostly ignored in standard seismic hazard assessment even though earthquake clustering is well known. In this work, we attempt to quantify the impact of more realistic dynamics on the seismic hazard estimations. We include the time and space dependences between earthquakes into the hazard analysis via Monte Carlo simulations. Our target region is the Lower Rhine Embayment, a low seismicity area in Germany. Including aftershock sequences by using the epidemic type aftershock-sequence (ETAS) model, we find that on average the hypothesis of uncorrelated random earthquake activity underestimates the hazard by 5–10 per cent. Furthermore, we show that aftershock activity of past large earthquakes can locally increase the hazard even centuries later. We also analyse the impact of the so-called long-term behaviour, assuming a quasi-periodic occurrence of main events on a major fault in that region. We found that a significant impact on hazard is only expected for the special case of a very regular recurrence of the main shocks.     

Resolution and significance assessment of precursory changes in mean earthquake magnitudes

The Loma Prieta earthquake (magnitude 7.0), which occurred in October 1989 in central California, was preceded by a period during which the mean magnitude of background seismicity in a small region near the eventual epicentre was abnormally low. This period may have begun as early as 1979, and it continued until mid-1988, after which the mean magnitude increased to a higher than normal value until the main earthquake. These changes were observed in the seismicity of an area 40  km in radius, centred on the Loma Prieta epicentre, and are consistent with the predictions of fracture mechanics studies. The 1988 change correlates with a reported change in long-term strain.
  A procedure has been developed for resolving such temporal changes in seismicity using CUSUM statistics. It demonstrates that the anomaly was highly significant, on the basis of analyses of two independent catalogues. There was also a significant anomaly before the 1994 Northridge earthquake.
  The hypothesis that large earthquakes are preceded by periods in which the mean magnitude of background activity is abnormally low, in the immediate vicinity of the eventual epicentre, is a tantalizing one. The analysis tool examined here may be useful for resolving such changes. Care needs to be taken, however, in routine surveillance of earthquake populations that contain large aftershock sequences.     

Locally triggered seismicity in the central Swiss Alps following the large rainfall event of August 2005

In 2005 August, an unusual series of 47 earthquakes occurred over a 12-hr period in central Switzerland. The earthquakes occurred at the end of 3-d period of intensive rainfall, with over 300 mm of precipitation. The highest seismicity occurred as two distinct clusters in the region of Muotatal and Riemenstalden, Switzerland, a well-known Karst area that received a particularly large amount of rainfall. The large increase in seismicity, compared to the background, and the short time delay between the onset of the intense rainfall and the seismicity strongly suggest that earthquakes were triggered by rainfall. In our preferred model, an increase in fluid pressure at the surface due to a large amount of rain leads to a local increase in pore fluid pressure at depth. The increase in pore fluid pressure will reduce the shear strength of a porous medium by counteracting normal stress and, at the end, provoke failure. The series of triggered earthquakes in central Switzerland occurred in regions that have been seismically active in the past, showing similar hypocentre locations and magnitudes. This suggests that these earthquakes occurred on existing faults that were critically stressed. We modelled the intense rainfall as a step increase in fluid pressure at the surface that migrates to greater depths following the solution of the one-dimensional diffusion equation in a homogeneous half space. This allowed us to estimate the hydraulic diffusivity by plotting triggered seismicity in a time–depth plot. We found values of hydraulic diffusivity in the range from 0.01 to 0.5 m² s⁻¹ for our study area. These values are in good agreement with previous studies on earthquakes that were triggered by fluids, supporting the idea that the observed earthquake series was triggered by the large amount of rainfall.     

A test of the precursory accelerating moment release model on some recent New Zealand earthquakes

The proposal that the moment release rate increases in a systematic way in a large region around a forthcoming large earthquake is tested using three recent, large New Zealand events. The three events, 1993–1995, magnitudes 6.7–7.0, occurred in varied tectonic settings. For all three events, a circular precursory region can be found such that the moment release rate of the included seismicity is modelled significantly better by the proposed accelerating model than by a linear moment release model, although in one case the result is dubious. The 'best' such regions have radii from 122 to 167 km, roughly in accord with previous observations world-wide, but are offset by 50–60 km from the associated main shock epicentre. A grid-search procedure is used to test whether these three earthquakes could have been forecast using the accelerating moment release model. For two of the earthquakes the result is positive in terms of location, but the main shock times are only loosely constrained.     

Regional disparity in the changes of agricultural land use intensity in China during 1980–2002

Based on the cost-benefit data (1980–2002) of farm products and China Agriculture Yearbooks, this paper studies the regional disparity in the changes of the agricultural land use in China during the period 1980–2002 from three aspects such as the degree of intensity, the sown area and the abandoned farmland. The results show that: (1) The degree of intensity of land use in the western region during 1980–2000 has a strong uptrend, but in the eastern and central regions the degree of intensity descends obviously and has shown a continuous downtrend since 1997. (2) The total sown area shrinks notably in the eastern region, while it enlarges constantly in the western region. (3) The sown area in the eastern, central and western regions has gone through a similar cyclic process: down (1980–1985)–up (1985–1991)–down (1991–1994)–up (1994–1999)–down (1999–2002). However, there are obvious differences in amplitude variation and tendency among them. The sown area has shrunk in the eastern region and expanded in the central and western regions especially before 1999. (4) The most cases of abandoned farmland are reported in the central region, the second in the eastern region and the least in the western region. The abandonment phenomena chiefly occurred during 1992–1995 in the eastern region, and during 1998–2002 in the central region.     

Precursory seismic quiescence in the Mudurnu Valley, North Anatolian fault zone, Turkey

The seismicity rate in the Mudurnu Valley of Turkey was studied using an earthquake catalogue that reports events homogeneously down to magnitude 2.3 for the years 1985–1989, and covers the area between latitudes 40.2° and 41.0°N, and longitudes 30.0° and 31.5°E. During this period the only two main shocks, M = 4.0 and M = 4.3, occurred on 1988 September 6 and 1988 December 9 within about 30km of each other. A highly significant seismic quiescence is evident in the area surrounding these main shocks, while the seismicity rate in the rest of the area covered by the catalogue remains constant. the quiescence becomes more pronounced the smaller the area around the main shocks that is studied. the smallest areas that can be studied contain about 60 earthquakes and have dimensions of approximately 25km on each side. the decreases in seismicity rates are 50–80 per cent depending on the volume and period used for defining the quiescence. the quiescence started in 1988 January and lasted about seven months, with approximately 4.5 months of normal activity separating it from the main shock of December. the precursor time of 12 months for an M = 4.3 main shock is similar to those observed in California. It is concluded that it is possible to resolve precursory quiescence before moderate and large earthquakes in the Mudurnu area with the existing seismograph network.     

Fractal clustering of induced seismicity in The Geysers geothermal area, California

We analyse earthquakes recorded at The Geysers geothermal field in California, an area where industrial activity induces seismicity. The seismicity is characterized by the seismic b -value and D , the fractal dimension of earthquake hypocentres measured from sliding windows containing 200 events. We study a group of events strongly clustered around an injection well. Over most of the time period examined we find a positive correlation between b and D . However, during the initiation of injection into a new well we find instead a negative correlation. The differences in correlation are statistically significant at the 1 σ level but only marginally so at the 2 σ level. These results provide evidence for a transient change in the seismic mechanisms operating, and may be explained by a change from conditions of slow stress loading to rapid loading as a result of the build-up of the rate of water injection into the reservoir.     

Tectonics of the Marmara Sea region of Turkey: new evidence from microearthquake fault plane solutions

Summary. Composite and single-event fault plane solutions for microearthquakes in the Izmit Bay area of the Marmara Sea indicate right-lateral strike-slip motion and tension on this extension of the North Anatolian Fault. This interpretation is consistent with teleseismically determined fault-plane solutions obtained for large earthquakes on the Marmara Sea seismic lineation. Consideration of the microplate geometry of north-western Turkey, inferred from seismicity as well as earthquake mechanisms, suggests that the region comprises two seismotectonic units with differing styles of deformation. The (Anatolid) structures of south- and central-western Anatolia are undergoing major extension, whereas the (Pontid) structures of the Marmara Sea region are being sheared, resulting in a mixed regime of both strike-slip and extensional faulting.     

Source process and tectonic implications of the great 1975 North Atlantic earthquake

Summary. The Atlantic segment of the Africa–Europe plate boundary has usually been interpreted as a transform boundary on the basis of the bathymetric expression of the Gloria fault and dextral strike-slip first-motion mechanisms aligned along the Azores–Gibraltar line of seismicity. The 1975 May 26 earthquake ( M _s=7.9) was assumed to fit into this framework because it occurred in the general area of this line and has a similar first-motion focal mechanism (strike=288°, dip=72°, slip angle=184°). However, several anomalies cast doubt on this picture: the event is abnormally large for an oceanic transform event; a sizeable tsunami was excited; the aftershock area is unusually small for such a large event; and most significantly, the epicentre is 200 km south of the presumed plate boundary. The Rayleigh wave radiation pattern indicates a change in focal mechanism to one with a significant dip-slip component. The short duration of the source time history (20 s, as deconvolved from long-period P -waves), the lack of directivity in the Rayleigh waves, and the small one-day aftershock area suggest a fault length less than 80 km. One nodal plane of the earthquake is approximately aligned with the trace of an ancient fracture zone.
We have compared the Pasadena 1-90 record of the 1975 earthquake to that of the 1941 North Atlantic strike-slip earthquake (200 km to the NNW) and confirmed the large size of the 1941 event ( M =8.2). The non-colinear relationship of the 1975 and 1941 events suggests that there is no well-defined plate boundary between the Azores and Gibraltar. This interpretation is supported by the intraplate nature of both the 1975 event and the large 1969 thrust event 650 km to the east. This study also implies that the largest oceanic strike-slip earthquakes occur in old lithosphere in a transitional tectonic regime.     

Observation of Coriolis coupled modes below 1 mHz

We present observations of spectral energy at toroidal mode frequencies in vertical seismic recordings of the 1998 Balleny Islands earthquake. Since toroidal modes on a spherically symmetric, nonrotating Earth have horizontally polarized particle motion these observations call for an explanation. We first rule out local and instrumental effects as being responsible for the verticalcomponent signal of the toroidal modes ₀ T ₃ (0.59 mHz) and ₀ T ₄ (0.77 mHz). The global effects that we consider are general heterogeneous mantle structure, ellipticity of figure and rotation. We find that rotation through Coriolis coupling of loworder spheroidal and toroidal oscillations is the dominant mechanism.     

Seismicity and seismic gap in the Lesser Antilles arc and earthquake hazard in Guadeloupe

Summary. A seismic study of the Lesser Antilles arc has been carried out, first for the period 1950–1978, for which we can use local seismic networks to draw maps of instrumental seismicity, then for the period 1530–1950, for which we have catalogues of felt earthquakes. The striking feature of the spatial distribution of foci is the cluster of epicentres in the northern half of the arc; all large earthquakes ( M > 7.5) are located north of 14° latitude. Seismicity cross-sections through the arc show a variable dipping subduction zone along the arc; the deep seismic zone is steeper in the centre of the arc than on the extremity.
The time-space diagram for historical seismicity, and the evidence of a seismic gap at the east of Guadeloupe lead us to consider the northern half arc as a likely site for a large earthquake in the near future.
The seismic slip rate calculated from all major earthquakes since 1530 is of much greater value than that obtained from recent plate tectonic models, suggesting that the recurrence rate of earthquakes is more than many hundreds of years with a possible aseismic creep.     

New seismogenic source and deep structures revealed by the 1999 Chia-yi earthquake sequence in southwestern Taiwan

In a tectonically active setting large earthquakes are always threats; however, they may also be useful in elucidating the subsurface geology. Instrumentally recorded seismicity is, therefore, widely utilized to extend our knowledge into the deeper crust, especially where basement is involved. It is because the earthquakes are triggered by underground stress changes that usually corresponding to the framework of geological structures. Hidden faults, therefore, can be recognized and their extension as well as orientation can be estimated. Both above are of relevance for assessment on seismic hazard of a region, since the active faults are supposed to be re-activated and cause large earthquakes. In this study, we analysed the 1999 October 22 earthquake sequence that occurred in southwestern Taiwan. Two major seismicity clusters were identified with spatial distribution between depths of 10 and 16 km. One cluster is nearly vertical and striking 032°, corresponding to the strike-slip Meishan fault (MSF) that generated the 1906 surface rupture. Another cluster strikes 190° and dips 64° to the west, which is interpreted as west-vergent reverse fault, in contrast to previous expectation of east vergence. Our analysis of the focal solutions of all the larger earthquakes in the 1999 sequence with the 3-D distribution of all the earthquakes over the period 1990–2004 allows us reinterpret the structural framework and suggest previously unreognized seismogenic sources in this area. We accordingly suggest: (1) multiple detachment faults are present in southwestern Taiwan coastal plain and (2) additional seismogenic sources consist of tear faults and backthrust faults in addition to sources associated with west-vergent fold-and-thrust belt.     

基于信息扩散理论的中国西部地区地震风险评估

利用全国地震目录数据库,系统梳理建国以来中国西部地区5.0级以上地震信息;利用信息扩散理论,从年最大震级和各级地震灾害年频次两个方面对西部12省、市、自治区的地震灾害进行风险评估。在此基础上,将风险估计离散值进行曲线拟合,系统评价中国西部地震最为多发5省区的地震风险分布趋势及差异,定量揭示中国西部地区地震风险,为西部地区地震防灾减灾提供技术支持和决策建议。     

Application of the negative binomial to earthquake occurrences in the Alpide-Himalayan belt

Summary . Shallow focus earthquakes ( h ≤ 60 km) of magnitude range M = 4.0–6.0, occurred during 1954–75 in various high seismicity zones of the Alpide-Himalayan belt have been tested by the Poisson and the negative binomial laws. When the clustering of events make the simple Poisson model inapplicable in most of the high seismicity zones of the Alpide-Himalayan belt it has been shown that the negative binomial entries provide an excellent model for describing the earthquake occurrences. The chi-square ( X ²) test is employed for testing the actual observations with theoretical distributions.     

Characteristics of change in the Antarctic sea ice area and its relation to SST in the tropical Pacific

１ＩｎｔｒｏｄｕｃｔｉｏｎＢｅｆｏｒｅｔｈｅ１９８０ｓ,Ａｎｔａｒｃｔｉｃｉｎｖｅｓｔｉｇａｔｉｏｎｓａｎｄｔｈｅｉｒｒｅｓｅａｒｃｈｄａｔａｗｅｒｅｒｅｌａｔｉｖｅｌｙｆｅｗ．Ｓｉｎｃｅｔｈｅ１９８０ｓ,ａｌｏｎｇｗｉｔｈｔｈｅｄｅｖｅｌｏｐｍ...     

南极重力场特征及界面的计算与解释

本文根据ＯＳＵ９１Ａ地球重力场模型和地形资料，计算了南极洲的空间重力异常和布格重力异常，分析了空间重力异常变化剧烈的原因及其与高程的相关性，同时根据布格重力异常用两层界面的反演方法计算了冰盖厚度和地壳厚度，冰盖厚度较大的地区位于极区的东南部，而极区周围和西南部地区厚度较小。地壳较厚的地区位于极区的东南部，最大达５６ｋｍ，西南部地区地壳较薄，最小值为８ｋｍ。     

Postseismic deformation at El Asnam (Algeria) in the seismotectonic context of northwestern Algeria

We use a combination of seismicity. tectonic features, focal mechanisms, seismic strain and postseismic movement to study the western part of North Algeria, the El Asnam region and its surrounding area in particular. A seismotectonic map of this part of Algeria, delimited by the Mediterranean Sea in the north and the Tellian mountains in the south, was built from available geological and seismological data. An examination of this map shows that the most significant earthquakes are concentrated along tectonic features and quaternary basins elongated in an east-west direction, suggesting NNW-SSE compressional movements. During the large El Asnam earthquake of 1980 October 10, M _w= 7.1, vertical movement was measured along a 40 km northeast-southwest thrust fault. These movements were determined geodetically in 1981 with reference to a basic network previously measured in 1976. In order to control postseismic movement and to ensure the monitoring of the seismic area, a dense geodetic network has been regularly measured since 1986, both in planemetry and altimetry. The results of the altimetric remeasurements show significant vertical movements. The elevation changes of the benchmarks have been deduced from precise levelling measurements: a remarkable uplift (5.1 ± 1.9 mm yr⁻¹) of the northwestern block, during the 1986-91 period has been observed, whereas the southeastern block is seen to be relatively stable. The Sar El Marouf anticline, situated along the central segment of the El Asnam surface breaks, appears to be growing with a maximum postseismic slip rate of (9.6 ± 1.4 mm yr⁻¹). The mean uplift rates computed for the northwestern block support the view that the 1954 earthquake did not occur on the same reverse fault as the 1980 event.     

Seismic risk in southern Europe through to India examined using Gumbel's third distribution of extreme values

Summary. A technique is described for the analysis of seismicity using Gumbel's third asymptotic distribution of extreme values. Seismicity of southern Europe through to India, nominally for the period 1900–74, is subdivided in a cellular manner, without recourse to tectonic discrimination between regions, and a covariance analysis on the three parameters of Gumbel's distribution is performed for each cell of seismicity. The results indicate that the upper bound to the magnitude of earthquake occurrence is often uncertain although it is discernible, while curvature of the earthquake occurrence distributions is usually established. Uncertainties in the forecasts of largest earthquakes, with a return period of 75 yr, are distinctly improved by taking into account the large and negative covariance which is measured between the curvature and upper bound to earthquake magnitude for the observed seismicity. These results are then used to map seismic risk for southern Europe through to India.     

近30年南极海冰的变化特征

采用NCEP的1973-2002年南极海冰密集度资料,对近30年南极海冰冰密集度的季节变化、年际变化及其与南极海冰涛动指数的长期变化关系进行了分析研究。结果表明,南极海冰的季节变化特点是海冰融化速度远大于凝结速度,而北极海冰融化速度与凝结速度基本相同。南极海冰存在着明显的年际变化,海冰面积指数呈增加趋势,年平均倾向率为28/10a。而北极海冰年际变化则相反,呈减少趋势,年平均面积指数的倾向率-3.5/10a。南极海冰涛动指数能代表南极地区近1/3的海水变化,是南极海冰变化的重要指数,具有10年、3-5年和2年左右的准振荡周期。     

Spatio-temporal patterns of land use change along the Bohai Rim in China during 1985–2005

Based on TM image data and other survey materials, this paper analyzed the spatiotemporal patterns of land use change in the Bohai Rim during 1985–2005. The findings of this study are summarized as follows: (1) Land use pattern changed dramatically during 1985–2005. Industrial and residential land in urban and rural areas increased by 643,946 hm2, of which urban construction land had the largest and fastest increase of 294,953 hm2 at an annual rate of 3.72%. (2) The outward migration of rural population did not prevent the expansion of residential land in rural areas by 184,869 hm2. This increase reveals that construction of rural residences makes seriously wasteful and inefficient use of land. (3) Arable land, woodland and grassland decreased at a rate of –0.02%, –0.12% and –1.32% annually, while unused land shrank by 157,444 hm2 at an annual rate of –1.69%. (4) The change of land use types showed marked fluctuations over the two stages (1985–1995 and 1995–2005). In particular, arable land, woodland and unused land experienced an inversed trend of change. (5) There was a significant interaction between arable land and woodland. Industrial construction land in urban and rural areas showed a net trend of increase during the earlier period, but only adjustment to its internal structure during the second period. The loss of arable land to the construction of factories, mines and residences took place mainly in the fringe areas of large and medium-sized cities, along the routes of major roads, as well as in the economically developed coastal areas in the east. Such changes are closely related to the spatial differentiation of the level of urbanization and industrialization in the region.