地震丛集的分形新方法—物理分形

地震时间丛集的传统分形分析方法把地震事件看作是时间轴上无质量的数学点。这些方法得到的结果主要反映了大量小地震的活动特征。本文发展了一个新的物理分形分析的筛子方法，采用一系列震级筛子来研究不同震级地震的分布特征。对华北和南加洲的研究结果表明，同一地区地震目录中不同震级范围的子集合具有近似相等的分维数，华北Ｄ＝０．３，南加洲Ｄ＝０．４。     

台湾地区地震的空间关联维特征与构造环境研究

采用关联维方法对台湾地区地震活动的空间特征进行了研究。先利用 10 0a来台湾的地震目录计算各个地震区、带的关联维数 ,将地震空间分布的分形特征定量表达出来 ,然后综合分析地震空间分布的关联维数和孕震构造环境之间的关系 ,得出了以下结论 :1)台湾东、西部地震区由于地震属于不同的大地构造单元 ,因此关联维数有较大的差异 ;2 )在各地震区内部的各个地震带由于板块构造、地壳结构、活断层分布上的差异 ,而具有与其构造特征相对应的关联维数 ;3)各地震带内部的各个不同的部位又由于不同的构造应力场 ,而导致地震分布上出现不同的丛集性 ,表现为不同的关联维数。这些结论充分说明通过关联维分析所得到的地震活动的空间图像与地震活动所代表的不同地质构造背景有着良好的对应关系     

Temporal Changes of Shallow Seismic Velocity Around the Karadere-Düzce Branch of the North Anatolian Fault and Strong Ground Motion

We analyze temporal variations of seismic velocity along the Karadere-Düzce branch of the north Anatolian fault using seismograms generated by repeating earthquake clusters in the aftershock zones of the 1999 M_w7.4 İzmit and M_w7.1 Düzce earthquakes. The analysis employs 36 sets of highly repeating earthquakes, each containing 4–18 events. The events in each cluster are relocated by detailed multi-step analysis and are likely to rupture approximately the same fault patch at different times. The decay rates of the repeating events in individual clusters are compatible with the Omori's law for the decay rate of regional aftershocks. A sliding window waveform cross-correlation technique is used to measure travel time differences and evolving decorrelation in waveforms generated by each set of the repeating events. We find clear step-like delays in the direct S and early S-coda waves (sharp seismic velocity reduction) immediately after the Düzce main shock, followed by gradual logarithmic-type recoveries. A gradual increase of seismic velocities is also observed before the Düzce main shock, probably reflecting post-seismic recovery from the earlier İzmit main shock. The temporal behavior is similar at each station for clusters at various source locations, indicating that the temporal changes of material properties occur in the top most portion of the crust. The effects are most prominent at stations situated in the immediate vicinity of the recently ruptured fault zones, and generally decrease with normal distance from the fault. A strong correlation between the co-seismic delays and intensities of the strong ground motion generated by the Düzce main shock implies that the radiated seismic waves produced the velocity reductions in the shallow material.     

Comparing seismicity declustering techniques by means of the joint use of Allan Factor and Morisita index

In this paper, we propose to compare different declustering methods on the basis of the time-correlation and the space-clustering of the residual earthquake catalog after the declustering techniques have been applied. To this aim, we applied two point process clustering measures, the Allan Factor and the Morisita Index, for the identification and quantification of temporal correlation and spatial clustering in point processes, respectively. We used our joint space–time approach to study the earthquake space–time point processes of southern California and Switzerland with surrounding area, declustered by using the method of Gardner and Knopoff (with Grünthal and Uhmhammer window) and that of Reasenberg (with different setting parameters). Our results show that the residual declustered catalog is still characterized by time-correlated structures at long timescales; however, the cutoff timescale that is the lowest timescale above which the time-correlation is visible is higher with the Reasenberg method while is smaller with the Gardner and Knopoff method with Grünthal window. The space-clustering analysis performed by means of the Morisita Index suggests that the declustering technique effectively reduces the spatial clustering of the seismicity of Switzerland, but does not change the spatial properties of the residual seismic catalogue of the southern California.     

Monte Carlo Studies of Relations between Fractal Dimensions in Monofractal Data Sets

—Within the fractal approach to studying the distribution of seismic event locations, different fractal dimension definitions and estimation algorithms are in use. Although one expects that for the same data set, values of different dimensions will be different, it is usually anticipated that the direction of fractal dimension changes among different data sets will be the same for every fractal dimension.¶Mutual relations between the three most popular fractal dimensions, namely the capacity, cluster and correlation dimensions, have been investigated in the present work. The studies were performed on the Monte Carlo generated data sets. The analysis has shown that dependence of the fractal dimensions on epicenter distribution, and relations among the fractal dimensions, are complex and variable. Neither values nor even inequalities among dimension estimates are preserved when different fractal dimensions are used. The correlation and the capacity dimensions seem to be good tools to trace collinear tendencies of eipicenters while the cluster dimension is more appropriate to studying uniform clustering of points.     

天津地区井水位气压系数变化与地震关系的探讨

地下水位气压系数在震前出现异常变化的研究,在地震预测研究中具有积极的作用和现实意义。利用高采样率的数字化资料,采取一阶差分筛选回归分析的方法,对天津地区4口地下水动态观测井的气压系数进行计算分析,并对其与地震的关系进行探讨。研究认为:不同的井孔,由于自身条件及干扰因素的不同,气压系数在震前的异常变化形态及幅度都有所不同;气压系数是随时间和空间而变化的物理量。通过文安5.1级地震震例检验,提出利用气压系数预报地震的可能性,但将之应用到地震监测预报实践中,尚需更多震例检验,并需开展广泛、系统的研究。     

Capturing the essential spatial variability in distributed hydrological modelling: Hydraulic roughness

The spatial variability of each parameter affecting storm runoff must be accounted for in distributed modelling. The objective of the work reported here is to assess the effects of using distributed versus lumped hydraulic roughness coefficients in the modelling of direct surface runoff. A spatially variable data set composed of Manning roughness coefficients is used to model direct surface runoff. To assess the information content (as measured by entropy) of spatially variable data and its significance in distributed modelling, various degrees of smoothing are applied. The error resulting from smoothing the hydraulic roughness coefficients is determined by modelling overland flow using a finite element solution. The Manning roughness coefficients were taken from field measurements of the Manning roughness coefficient at 0.6 m on a 14 m hillslope. These values were then used in a numerical simulation of outflow hydrographs to investigate the dependence of error on spatial variability. Our study focuses on the characteristics of spatial data used in distributed hydrological modelling. The field sites have fractal dimensions of ≈? 1.4, which is close to a Brownian variation. The sampling interval that captures the essential spatial variability of the Manning roughness coefficient does not seem to matter due to its Brownian variation in the field sites. Hence due to the nearly uniform random distribution, measurements at 0.6 m intervals are not necessary and larger intervals would yield results that are just as acceptable provided the mean value together with a uniformly random distribution is maintained for any size of finite element or sampling resolution. Because detailed measurements of hydraulic roughness are not practically available for deterministic catchment modelling, it is important to know that larger sampling resolutions may be used than 0.6 m.     

Fractal analysis of 2-D fracture networks of geothermal reservoirs in south-western Turkey

Natural fracture patterns of producing geothermal formations in south-western Turkey are mapped at different scales. The fractal dimensions of different fracture network properties, such as spatial distribution, density, connectivity, orientation, and length are measured by different methods. Analysis of the natural fracture patterns from giga to microscales identifies the descending behavior of box-counting fractal dimension with respect to the scale. It is observed that the fracture networks represent scale-invariant properties, but fractal dimensions might notably differ when the mass dimension is measured applying different methods. Anisotropic nature of fracture networks is also included in the fractal analysis.     

1983年菏泽5.9级地震前后小震频数及能量空间分布的分形特征

本文对菏泽地震为例,从频数及能量两方面对地震不均匀空间颁的分形特征进行研究,并探讨其预报含义。结果表明,地震频数及地震能量空间分布均表现出统计意义上的自相似特征,并且前者的分布结构较后者列为复杂。菏泽地震前几年的“降维”中视为中期预报的一个显著指标,而动态均值则可作为判断“降维”的一个较好标准;ｔ检验表明能量分布的降维较频数分布更为显著,并且其对中强地震孕育过程各阶段的反映似也列为灵敏,降维之后的     

分形方法在地震序列类型早期判别中的应用

着重讨论了地震序列的空间分布随时间变化的特性，利用分形方法提取大震后余震的空间分维信息，从而达到早期判别地震序列类型的目的。以１９６６年以来定位精度较好、余震记录较全的１２个７级以上大震和大同６级强震序列为资料，计算了震后１个月内余震随时间变化的空间分维数。结果发现，不同类型的地震序列分维数不同。一般震后５－６天，主震型序列的分维数在１．１－１．２之间，多震型序列的分维数在１．３以上。利用前７天的     

Magma intrusion as a driving mechanism for the seismic clustering following the 9 May 1989 earthquake swarms at the Canary Islands

On 9 May 1989 a M _L = 5.2 earthquake struck a region between the islands of Tenerife and Gran Canaria. We investigated the time-spatial evolution of seismic patterns affecting the Canary Islands region during 1989–1995, using a quantitative spatial fractal analysis method. This method allows quantitative investigation of subtle trends in seismicity distribution through time. The fractal analysis indicates that epicenters clustered around a large zone during the May 1989 sequence affected narrow zones during 1991–1993, but then larger zones during 1993–1995 with an overall trend to shallower focal depths. The spatial localisation of seismic data and its time evolution appear to be related to magmatic rather than tectonic activity. Spatial clustering properties of seismicity are consistent with a major intrusive episode in 1989, followed by a period of quiescence and renewed deep intrusive activity from 1993 onwards. This interpretation suggests an increasing probability of future volcanic hazard in the region investigated.     

Probabilistic seismic hazard assessment of Eastern Marmara Region

The objective of this study is to evaluate the seismic hazard in Eastern Marmara Region using an improved probabilistic seismic hazard assessment methodology. Two significant improvements over the previous seismic hazard assessment practices are accomplished in this study: advanced seismic source characterization models in terms of source geometry and recurrence relationships are developed, and improved global ground motion models (NGA-W1 models) are employed to represent the ground motion variability. Planar fault segments are defined and a composite magnitude distribution model is used for all seismic sources in the region to properly represent the characteristic behavior of the North Anatolian Fault without the need for an additional background zone. Multi-segment ruptures are considered using the rupture model proposed by the Working Group on California Earthquake Probabilities (2003). Events in the earthquake catalogue are attributed to the fault zones and scenario weights are determined by releasing the accumulated seismic energy. The uniform hazard spectra at 10 % probability of exceedance in 50 years hazard level for different soil conditions (soil and rock) are revealed for specific locations in the region (Adapazar?, Düzce, Gemlik, Izmit, Iznik and Sapanca). Hazard maps of the region for rock site conditions at the selected hazard levels are provided to allow the readers perform site-specific hazard assessment and develop site-specific design spectrum for local site conditions.     

关于地震丛集特征、成因及临界状态的讨论

给出了研究地震临界丛集现象的最新方法,即时间变异诊断方法。讨论了不同类型的地震丛集活动特征和成因,结合中国大陆及邻区M≥ 8. 0大震序列和华北地区自1966年以来M≥6. 0地震序列,研究了丛集状态的自相似性及临界时间分支现象。对发震系统的内部时间给出了定义,并研究了地震丛集的非线性时间结构特征。研究认为: ① 事件的发生意味着系统现存的能量一次跳跃式的释放和状态的一次转换。只有当系统以某种足够随机方式动作时,系统状态才会转换,事件才会发生;② 涨落是系统状态的探测器,在临界状态时系统出现很大的涨落,时间变异系数W= 1。临界丛集是由系统内在随机性决定的,事件发生时间是随机地出现的,未来可能出现不止一次的W= 1的临界点;③ 相继发生的事件(M≥ M₀) 可以反映系统的演化过程和特征,丛集是系统演化过程中的有序性和无序性的综合反映;④ 事件丛集和丛集的事件具有自相似性和分形结构,也具有耗散系统的自组织和映射特点;⑤ 根据时间变异系数,系统内部时间(年龄),时间间隔与所经历时间非线性关系式以及映射和迭代方法可以用来预测未来地震趋势和未来M≥ M₀ 地震发生的可能时间点。     

Fractal analysis of earthquake swarms of Vogtland/NW-Bohemia intraplate seismicity

The special type of intraplate microseismicity with swarm-like occurrence of earthquakes within the Vogtland/NW-Bohemian Region is analysed to reveal the nature and the origin of the seismogenic regime. The long-term data set of continuous seismic monitoring since 1962, including more than 26000 events within a range of about 5 units of local magnitude, provides an unique database for statistical investigations. Most earthquakes occur in narrow hypocentral volumes (clusters) within the lower part of the upper crust, but also single event occurrence outside of spatial clusters is observed. Temporal distribution of events is concentrated in clusters (swarms), which last some days until few month in dependence of intensity. Since 1962 three strong swarms occurred (1962, 1985/86, 2000), including two seismic cycles. Spatial clusters are distributed along a fault system of regional extension (Leipzig-Regensburger Störung), which is supposed to act as the joint tectonic fracture zone for the whole seismogenic region. Seismicity is analysed by fractal analysis, suggesting a unifractal behaviour of seismicity and uniform character of seismotectonic regime for the whole region. A tendency of decreasing fractal dimension values is observed for temporal distribution of earthquakes, indicating an increasing degree of temporal clustering from swarm to swarm. Following the idea of earthquake triggering by magma intrusions and related fluid and gas release into the tectonically pre-stressed parts of the crust, a steady increased intensity of intrusion and/or fluid and gas release might account for that observation. Additionally, seismic parameters for Vogtland/NW-Bohemia intraplate seismicity are compared with an adequate data set of mining-induced seismicity in a nearby mine of Lubin/Poland and with synthetic data sets to evaluate parameter estimation. Due to different seismogenic regime of tectonic and induced seismicity, significant differences between b-values and temporal dimension values are observed. Most significant for intraplate seismicity are relatively low fractal dimension values for temporal distribution. That observation reflects the strong degree of temporal earthquake clustering, which might explain the episodic character of earthquake swarms and support the idea of push-like triggering of earthquake avalanches by intruding magma.     

The application of the method of reconstruction a phase space with higher dimensions to the study of seismic processes and earthquake prediction

The method of reconstructing a phase space with higher dimensions can be applied to analysing the observaiton data for studying the dynamic behavior of seismic processes and earthquake prediction. It is found that seismic processes are chaotic. So earthquakes could be determinately predicted only within the predictable period. Here, some related problems (e.g. several primary parameters and the relations among them, the relation between the quality of used data and analytical results, the relation between a reductive process of the fractal dimensions and earthquake prediction and so on) are discussed. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 470–476, 1993. This study was sponsored by the Chinese Joint Seismological Science Foundation under Contract NO. 9000007.     

Multifractality in Seismicity Spatial Distributions: Significance and Possible Precursory Applications as Found for Two Cases in Different Tectonic Environments

We explore fractal properties of two observed seismicity distributions prior to the 2003 M _w 7.4 Colima, Mexico and 1992 M _w 7.3 Landers, USA earthquakes, together with several mathematical fractal distributions and two non-fractal ones, in order to estimate minimum reliable sample sizes, determine whether fractality for observed seismicity is essentially different from random uniform distributions, and explore the possibility of extracting premonitory information from fractal characteristics of seismicity before large earthquakes. Sample sizes above 800 events for whole catalogs appear to be sufficient to maintain ordered multifractality and to yield dimension estimates that vary smoothly and reliably. Fractal estimates appear to be best for whole catalogs that include aftershocks. The fractal characteristics of spatial distributions of seismicity are essentially different from those of the uniform random distribution, which is the null hypothesis of a non-fractal distribution with minimum information. The fractal dimensions and afractality measures of seismicity distributions change with time and show distinctive behaviors associated with foreshocks and main events, although these behaviors are different for each example. Results suggest the possibility of a priori identification of foreshocks to large earthquakes. A combination of fractal dimension and afractality measures over time may be helpful in large earthquake premonitory studies.     

Soil moisture assessment in the alluvial basins associated with seismic activity of northwest Anatolia,Turkey, using Landsat TM/ETM satellite imagery / Estimation de l’humidité du sol associée à l’activité sismique dans les bassins alluviaux du nord-ouest de l’Anatolie (Turquie) par l’imagerie satellitaire Landsat TM/ETM

Abstract

Abstract The aim of this study is to detect seasonal surface temperature changes and to estimate soil moisture conditions based on the evaporative cooling principle of damp ground in the alluvial basins of northwest Anatolia, Turkey, using Landsat TM/ETM data. According to analysis of satellite sensor data acquired on different dates, soil moisture is greatest in the spring season in the basins. Soil moisture decreases toward the summer and autumn. The 17 August 1999 earthquake occurred in the high surface temperature (low soil moisture) period, and the 12 November 1999 earthquake occurred in the low surface temperature (high soil moisture) period. It is possible to conclude that the urban-rural settlements and industrial developments on the loose deposits of the Adapazari, ?zmit and Düzce depressions have been affected by the seasonal changes in the local ground conditions.     

地壳垂直形变的分维特征─以京、津、唐复测网为例

本文应用分数维的概念及方法，对京、津、唐地区１９６８─１９７５年的大面积地壳垂直形变复测成果数据处理后的结果进行了分维特征分析。结果表明，地震前该地区地壳垂直形变空间分布的分数维明显降低，地壳的有序性增强。本文认为，地壳垂直形变空间分布分数维能够明确地揭示地壳的有序程度，是一种较好地利用形变测量资料提取地震前兆异常的方法。     

Spatial distribution of aftershocks and background seismicity in central California

We examine the spatial distribution of earthquake hypocenters in four central California areas: the aftershock zones of the (1) 1984 Morgan Hill (2) 1979 Coyote Lake, and (3) 1983 Coalinga earthquakes, as well as (4) the aseismically creeping area around Hollister. The basic tool we use to analyze these data are frequency distributions of interevent distances between earthquakes. These distributions are evaluated on the basis of their deviation from what would be expected if earthquakes occurred randomly in the study areas. We find that both background seismic activity and aftershocks in the study areas exhibit nonrandom spatial distribution. Two major spatial patterns, clustering at small distances and anomalies at larger distances, are observed depending on tectonic setting. While both patterns are seen in the strike-slip environments along the Calaveras fault (Morgan Hill, Coyote Lake, and Hollister), aftershocks of the Coalinga event (a thrust earthquake) seem to be characterized by clustering only. The spatial distribution of earthquakes in areas gradually decreasing in size does not seem to support the hypothesis of a self-similar distribution over the range of scales studied here, regardless of tectonic setting. Spatial distributions are independent of magnitude for the Coalinga aftershocks, but events in strike-slip environments show increasing clustering with increasing magnitude. Finally, earthquake spatial distributions vary in time showing different patterns before, during, and following the end of aftershock sequences.     

天津市面向震害快速评估的房屋和人口空间化研究*

以天津市的大比例尺房屋空间数据和第六次人口普查数据为基础,利用GIS空间分析功能,模拟了天津市房屋和人口空间分布,产出1km格网的房屋数据和分时段人口公里格网数据。经检验,模拟的房屋和人口数据达到了较高精度,能够较准确地表达行政单元内部的房屋和人口空间分布。将模拟结果应用到设定地震震害评估中,与基于平均密度法的评估结果对比,表明基于公里格网评估的结果更为合理,减少了由于房屋和人口分布不均匀造成的误差,能够为地震应急工作提供可靠的数据支撑。