Correlation Fractal and Multifractal Characteristics of Seismic Activity in the Taiwan Area, China

Based on the analysis of newly collected data of plate tectonics, distribution of active faults and crustal deformation, the Taiwan area is divided into two seismic regions and six seismic belts. Then, correlation fractal dimensions of all the regions and belts are calculated, and the fractal characteristics of hypocenteral distribution can be quantitatively analyzed. Finally, multifractal dimensions Dq and f(α) are calculated by using the earthquake catalog of the past 11 years in the Taiwan area. This study indicates that (1) there exists a favorable corresponding relationship between spatial images of seismic activity described with correlation fractal dimension analysis and tectonic settings; (2) the temporal structure of earthquakes is not single but multifractal fractal, and the pattern of Dq variation with time is a good indicator for predicting strong earthquake events.     

Preliminary analysis of the self-similarity of the aftershocks of the Japanese earthquake on March 11, 2011

The quantitative parameters of the self-similarity of the aftershocks of the Japanese earthquake on March 11, 2011 were obtained. The parameters p in the Omori law (1.06), b in the Gutenberg-Richter law (0.61), and the fractal dimension (D) of the earthquake epicenters (1.52) were determined. Self-similarity is manifested in a range of two orders of temporal and spatial scales and four units of magnitude. The stability in time of parameter p and spatial variations in b and p parameters were revealed.     

四川地区断层空间分布的多重分形特征

断层系统多重分形避免了容量维Ｄ０的不足，考虑了每一条断层对整体分维的贡献，反映了断层的空间分布特征。根据四川地区断层分布图，运用多重分形理论检验了断层空间分布的自相似性，分区计算了四川地区各带的多重分维值。结果表明，断层破裂系统为分形结构，具有很好的自相似性，相关系数均达０．９９。不同带分维值的较大差异，反映不同构造区域的断层分布特点。最后，对分维的含义进行了深入探讨，对油气勘探的有利区域进行了分析。     

An assessment of seismicity parameters in northwest Himalaya and adjoining regions

In this study, an assessment of seismicity parameters in the northwest Himalaya and adjoining regions using an earthquake catalog from India Meteorological Department covering a period from June 1, 1998 to June 30, 2011 has been carried out. The spatial distributions of seismicity parameters, namely magnitude of completeness, M _C, a value, b value, and correlation fractal dimension, D _C, are estimated for the studied region. The M _C, a, and b values are found to be 2.5, 4.601, and 0.83, respectively. Despite significant gaps, the spatial distributions of a and b values are seen to follow similar trend and are found scattering in between Main Boundary Thrust (MBT) and South Tibet Detachment, adjoining areas of Mahendragarh-Dehradun Fault (MDF), Delhi-Haridwar Ridge (DHR) and Moradabad Fault (MF), and the southern flank of Karakoram Fault and Indus-Tsangpo Suture Zone. The estimated spatial distribution of b and a values is within 90 % of confidence level, thereby indicating non-uniform stress accumulation or higher rock fracturing density in the studied region caused by strong tectonization following several earthquakes. Negative correlation between low b value and high D _C is observed predominantly in the region between the MBT and Munsiari Thrust or Main Central Thrust-I of Garhwal and Kumaon Himalaya, adjoining zones of MDF, DHR, and MF of Indo-Gangetic plain, and the eastern flank of the studied region, suggesting the presence of asperities in the zone. At the same time, active creeping process can be inferred in between the MBT and Main Central Thrust of Garhwal Himalaya and the surrounding areas of Shimla region of the Himalayan arc to the northwestern part of the studied region from the positive correlation between b value and D _C. The results indicate that the structural heterogeneity caused by different stress accumulation and rock fracturing densities exists due to continuous tectonic adjustments between different geomorphic features of the studied region. An attempt has also been made to classify the studied region into smaller seismic zones by observing the spatial patterns of b value and D _C that are fractal properties of the observed seismicity, along with the prevalent fault networks.     

Effects of errors and biases on the scaling of earthquake spatial pattern: application to the 2004 Sumatra–Andaman sequence

This study discusses the scaling properties of the spatial distribution of the December 26, 2004, Sumatra aftershocks. We estimate the spatial correlation dimension D ₂ of the epicentral distribution of aftershocks recorded by a local network operated by Geological Survey of India. We estimate the value of D ₂ for five blocks in the source area by using generalized correlation integral approach. We assess its bias due to finite data points, scaling range, effects of location errors, and boundary effects theoretically and apply it to real data sets. The correlation dimension was computed both for real as well as synthetic data sets that include randomly generated point sets obtained using uniform distributions and mimicking the number of events and outlines of the effective areas filled with epicenters. On comparing the results from the real data and random point sets from simulations, we found the lower limit of bias in D ₂ estimates from limited data sets to be 0.26. Thus, the spatial variation in correlation dimensions among different blocks using local data sets cannot be directly compared unless the influence of bias in the real aftershock data set is taken into account. They cannot also be used to infer the geometry of the faults. We also discuss the results in order to add constraints on the use of synthetic data and of different approaches for uncertainty analysis on spatial variation of D ₂. A difference in D ₂ values, rather than their absolute values, among small blocks is of interest to local data sets, which are correlated with their seismic b values. Taking into account the possible errors and biases, the average D ₂ values vary from 1.05 to 1.57 in the Andaman–Nicobar region. The relative change in D ₂ values can be interpreted in terms of clustering and diffuse seismic activity associated with the low and high D ₂ values, respectively. Overall, a relatively high D ₂ and low b value is consistent with high-magnitude, diffuse activity in space in the source region of the 2004 Sumatra earthquake.

     

Seismic source characteristics in Kachchh and Saurashtra regions of Western India: b-value and fractal dimension mapping of aftershock sequences

Seismic source characteristics in the Kachchh rift basin and Saurashtra horst tectonic blocks in the stable continental region (SCR) of western peninsular India are studied using the earthquake catalog data for the period 2006–2011 recorded by a 52-station broadband seismic network known as Gujarat State Network (GSNet) running by Institute of Seismological Research (ISR), Gujarat. These data are mainly the aftershock sequences of three mainshocks, the 2001 Bhuj earthquake (M _w 7.7) in the Kachchh rift basin, and the 2007 and 2011 Talala earthquakes (M _w ≥ 5.0) in the Saurashtra horst. Two important seismological parameters, the frequency–magnitude relation (b-value) and the fractal correlation dimension (D _c) of the hypocenters, are estimated. The b-value and the D _c maps indicate a difference in seismic characteristics of these two tectonic regions. The average b-value in Kachchh region is 1.2 ± 0.05 and that in the Saurashtra region 0.7 ± 0.04. The average D _c in Kachchh is 2.64 ± 0.01 and in Saurashtra 2.46 ± 0.01. The hypocenters in Kachchh rift basin cluster at a depth range 20–35 km and that in Saurashtra at 5–10 km. The b-value and D _c cross sections image the seismogenic structures that shed new light on seismotectonics of these two tectonic regions. The mainshock sources at depth are identified as lower b-value or stressed zones at the fault end. Crustal heterogeneities are well reflected in the maps as well as in the cross sections. We also find a positive correlation between b- and D _c-values in both the tectonic regions.

     

The perimeter-area fractal model and its application to geology

Perimeters and areas of similarly shaped fractal geometries in two-dimensional space are related to one another by power-law relationships. The exponents obtained from these power laws are associated with, but do not necessarily provide, unbiased estimates of the fractal dimensions of the perimeters and areas. The exponent (D_AL) obtained from perimeter-area analysis can be used only as a reliable estimate of the dimension of the perimeter (D_L) if the dimension of the measured area is D_A=2. If D_A<2, then the exponent D_AL=2D_L/D_A>D_L. Similar relations hold true for area and volumes of three-dimensional fractal geometries. The newly derived results are used for characterizing Au associated alteration zones in porphyry systems in the Mitchell-Sulphurets mineral district, northwestern British Columbia.     

Evaluation of earthquake potential along the Northern Anatolian Fault Zone in the Marmara Sea using comparisons of GPS strain and seismotectonic parameters

Seismotectonic parameters including the Gutenberg-Richter b-value and multifractal dimensions D₂ and D₁₅ of seismicity patterns (both spatial and temporal) were compared to GPS-derived maximum shear and dilatation strains measured in the Marmara Sea region of western Turkey along the Northern Anatolian Fault Zone (NAFZ). Comparisons of seismotectonic parameters and GPS-derived maximum shear and dilatation strain along the NAFZ in the vicinity of the 1999 M7.4 Izmit earthquake reveal a positive correlation (r = 0.5, p = 0.05) between average dilatation and the Gutenberg-Richter b-value. Significant negative correlation (r = − 0.56, p = 0.03 and r = − 0.56, p = 0.02) was also observed between the spatial fractal dimension D₂ and GPS-derived maximum geodetic and shear strain. This relationship suggests that, as maximum geodetic and shear strains increase, seismicity becomes increasingly clustered.Anomalous interrelationships are observed in the Marmara Sea region prior to the Izmit event along a bend in the NAFZ near the eastern end of the Marmara Sea known as the Northern Boundary Fault (NBF). An asperity is located near the northwest end of the NBF. Along the 50-km length of the NBF, GPS strains become slightly compressive. The correlation between b-value and GPS-derived dilatation suggests that regions in compression have increased probability of larger magnitude rupture. The NBF appears to serve as an impediment to the transfer of strain from east to west along the NAFZ. Recurrence times for large earthquakes along the NBF are larger than in surrounding areas. Temporal clustering of seismicity in the vicinity of the NBF may represent foreshocks of an impending rupture.     

黄骅坳陷中区断裂系统分形研究

利用分形几何方法分析了黄骅坳陷中区断裂分布的分形特征,以探讨该区断裂系统的发育规律及其与油气的关系。利用数盒子方法进行分析表明,黄骅坳陷中区断裂系统在标度区间1.75～7km具有很好的统计自相似性。根据断裂分维值等值线可知,海岸线两侧断裂分维值分布规律有着明显不同,显示了两侧断裂系统在形成机制上的差异性。研究区分维值自西向东逐渐增大,在歧口凹陷南北两侧有着最高的分维值。结合黄骅坳陷中区断裂系统的形成机制及油田的分布,初步探讨了断裂分形特征与黄骅坳陷中区断裂系统的形成机制及断裂分形特征与油田分布之间的关系。研究结果表明,断裂系统分维值的变化特征与断裂的发育规律有着很好的相关性,断裂发育程度高的地区其分维值也更高,也将更有利于油气的运聚。     

A comparison of seismic and structural measurements of scaling exponents during tensile subcritical crack growth

The observed fractal nature of both fault length distributions and earthquake magnitude-frequency distributions suggests that there may be a relationship between the structure of active fault systems and the resulting seismicity. In previous theoretical work, a positive correlation between the exponent D from the fracture length distribution, and the seismic or acoustic emission (AE) b-value has been inferred from a simple dislocation model of the seismic source. Here, we present the first experimental evidence for a correlation between D and b from a series of tensile fracture mechanics tests on crystalline rock, carried out in different environmental conditions, both air-dry and water-saturated, and at ambient temperature and pressure. The microseismic acoustic emissions were monitored during subcritical crack growth under controlled conditions of constant stress intensity, K_I, and quantitative analyses of the resulting fracture patterns were carried out on the same specimens. It is found that AE b-values, ranging from 1.0 to 2.3, correlate negatively with the normalized stress intensity K_I/KIC, where K_IC is the fracture toughness of the specimen. The microcrack length distribution exponent D, ranges from 1.0 to 1.7. Fluid presence has a first-order influence on both the AE and structure produced in these experiments. For experiments at low stress intensity or high fluid content, the activation of the stress corrosion mechanism for K_I < K_IC leads to a greater relative proportion both of small cracks and of low amplitude acoustic emissions, reflected in higher values of D and b. The exponent D is found to correlate positively with the AE b-value.     

Doubling of coastal erosion under rising sea level by mid-century in Hawaii

This study discusses the scaling properties of the spatial distribution of the December 26, 2004, Sumatra aftershocks. We estimate the spatial correlation dimension D ₂ of the epicentral distribution of aftershocks recorded by a local network operated by Geological Survey of India. We estimate the value of D ₂ for five blocks in the source area by using generalized correlation integral approach. We assess its bias due to finite data points, scaling range, effects of location errors, and boundary effects theoretically and apply it to real data sets. The correlation dimension was computed both for real as well as synthetic data sets that include randomly generated point sets obtained using uniform distributions and mimicking the number of events and outlines of the effective areas filled with epicenters. On comparing the results from the real data and random point sets from simulations, we found the lower limit of bias in D ₂ estimates from limited data sets to be 0.26. Thus, the spatial variation in correlation dimensions among different blocks using local data sets cannot be directly compared unless the influence of bias in the real aftershock data set is taken into account. They cannot also be used to infer the geometry of the faults. We also discuss the results in order to add constraints on the use of synthetic data and of different approaches for uncertainty analysis on spatial variation of D ₂. A difference in D ₂ values, rather than their absolute values, among small blocks is of interest to local data sets, which are correlated with their seismic b values. Taking into account the possible errors and biases, the average D ₂ values vary from 1.05 to 1.57 in the Andaman–Nicobar region. The relative change in D ₂ values can be interpreted in terms of clustering and diffuse seismic activity associated with the low and high D ₂ values, respectively. Overall, a relatively high D ₂ and low b value is consistent with high-magnitude, diffuse activity in space in the source region of the 2004 Sumatra earthquake.     

Fractal analysis of the fracture strength of lava dome material based on the grain size distribution of block-and-ash flow deposits at Unzen Volcano, Japan

A fractal theory of rock fragmentation is applied to block-and-ash flow deposits from the Fugendake dome, Unzen Volcano, Kyushu, Japan, in order to analyze the material strength and the energy required for size reduction of the source lava dome. Two fractal dimensions h and D_s, which are mutually interchangeable, represent the relative strength and energy for particles reduced to a given size. They can be theoretically estimated from the power relations of a reference grain size to the cumulative mass and number of fragments smaller than that size. The Unzen–Fugendake block-and-ash flow deposits have been further modified by size sorting and secondary fragmentation that occurred during flowage, so that the h value decreases (or D_s value increases) with increasing distance from the source. Coarse, reversely graded deposits are, however, found to retain the original size population relatively well. The D_s values estimated from deposits of this type are compatible with those previously reported from decompression–fragmentation experiments conducted on the same dome material. The employed fractal approach could thus give insights into the potential mode of dome collapse that generates block-and-ash flows.     

新疆东准噶尔地区断裂的多重分形研究

结合东准噶尔地区强应变构造带的空间展布情况和断裂的空间分布特征 ,用多重分形模型测算了各带的多重分维值。研究结果表明 :强应变构造带在不同尺度上具有自相似性 ;断裂的空间分布为多重分形结构 ,断裂的演化有从复杂几何结构的次级断裂组合向单一连续的大断裂过渡而分维值逐渐减小的趋势 ;不同的大地构造单元有不同的分维值 ,而且有靠近板块缝合带 (或海沟 )分维值减小、远离板块缝合带 (或海沟 )分维值增大的趋势     

Do Skeletal Networks Derived from Water Bodies Follow Horton's Laws?

The aim of this short note is to test whether the morphological skeletal network (MSN) of water bodies that resembles a river network follows Horton's laws. A fractal relationship of MSN of a water body is also shown. This investigation shows that the MSN of the Nizamsagar reservoir follows Horton's laws. Furthermore, this reservoir has a fractal dimension (D _m) of 1.92 which was computed by using two morphometric quantities and the fractal dimension of the main skeletal length (d). This value tallies exactly with the fractal dimension (D _f) of the whole MSN computed through box-counting method.     

The use of the fractal dimension to quantify the morphology of irregular-shaped particles

For several decades, sedimentologists have had difficulty in obtaining an efficient index of particle form that can be used to specify adequately irregular morphology of sedimentary particles. Mandelbrot has suggested the use of the fractal dimension as a single value estimate of form, in order to characterize morphologically closed loops of an irregular nature. The concept of fractal dimension derives from Richardson's unpublished suggestion that a stable linear relationship appears when the logarithm of the perimeter estimate of an irregular outline is plotted against the logarithm of the unit of measurement (step length). Decreases in step length result in an increase in perimeter by a constant weight (b) for particles whose morphological variations are the same at all measurement scales (self-similarity). The fractal dimension (D) equals 1.0-(b), where b is the slope coefficient of the best-fitting linear regression of the plot. The value of D lies between 1.0 and 2.0, with increasing values of D correlating with increasing irregularity of the outline. In practice, particle outline morphology is not always self-similar, such that two or possibly more fractal elements can occur for many outlines. Two fractal elements reflect the morphological difference between micro-scale edge textural effects (D₁) and macro-scale particle structural effects (D₂) generated by the presence of crenellate-edge morphology (re-entrants). Fractal calibration on a range of regular/irregular particle outline morphologies, plus examination of carbonate beach, pyroclastic and weathered quartz particles indicates that this type of analysis is best suited for morphological characterization of irregular and crenellate particles. In this respect, fractal analysis appears as the complementary analytical technique to harmonic form analysis in order to achieve an adequate specification of all types of particles on a continuum of irregular to regular morphology.     

Natural fracture profiles,fractal dimension and joint roughness coefficients

Summary Thirteen natural rock profiles (Barton and Choubey, 1977) are analyzed for their fractal properties. Most of the profiles were found to approximate fractal curves but some also showed features of specific wavelengths and amplitudes superimposed on fractal characteristics. The profiles showed fractal dimensions from 1.1 to 1.5 covering a range of selfsimilar and self-affine curves. The analysis results suggest a negative correlation between fractal dimension,D, and amplitude,A. Joint roughness coefficients (JRC) show a positive correlation with amplitude,A, and a negative correlation with fractal dimension,D. A numerical model of fracture closure is used to investigate the effects of different profile characteristics (D, A and sample size) on the nature of dilation and contact area, using the natural profiles and synthetic fractional Brownian motion profiles. Smooth profiles (low JRC, highD, lowA) display many small contact regions whereas rough fractures (high JRC, lowD, highA) display few large contact areas. The agreement with published experimental data supports the suggested correlations between JRC and the fractal parameters,A andD. It is suggested that observed scale effects in JRC and joint dilation can be explained by small differential strain discontinuities across fractures, which originate at the time of fracture formation.     

A model of characteristic earthquakes and its implications for regional seismicity

Regional seismicity (i.e. that averaged over large enough areas over long enough periods of time) has a size–frequency relationship, the Gutenberg–Richter law, which differs from that found for some seismic faults, the Characteristic Earthquake relationship. But all seismicity comes in the end from active faults, so the question arises of how one seismicity pattern could emerge from the other. The recently introduced Minimalist Model of Vázquez‐Prada et al. of characteristic earthquakes provides a simple representation of the seismicity originating from a single fault. Here, we show that a Characteristic Earthquake relationship together with a fractal distribution of fault lengths can accurately describe the total seismicity produced in a region. The resulting earthquake catalogue accounts for the addition of both all the characteristic and all the non‐characteristic events triggered in the faults. The global accumulated size–frequency relationship strongly depends on the fault length fractal exponent and, for fractal exponents close to 2, correctly describes a Gutenberg–Richter distribution with a b exponent compatible with real seismicity.     

吉林省水系构成的分形研究

根据对分形水系的新近认识探讨了吉林省水系结构的自相似规律。Horton-Strahler的水系标度定律隐含着水系的位序-规模法则和河流长度-流域面积的异速生长关系,这暗示着α＝ lnR_b／lnR_l是一种等级结构的维数,不能将之与空间结构维数混同;Hack模型的标度因子b＝lnR_l／lnR_a是一种广义的空间维数之比,不能据之确定主河道的分维。基于上述思想,对吉林省10个主要水系的等级结构进行考察,发现气候相对湿润的山区水系的α值高于气候相对干燥的平原地区水系的α值,而平原-干燥区水系的b值高于山地-湿润区水系的b值。从河流发育的地质、地貌背景和气候-水文关系等角度对上述现象进行了的初步解释,并根据洮儿河的异常α值修正了LaBabera-Rosso的水系分维定义。     

Quantification of Aperture and Relations Between Aperture,Normal Stress and Fluid Flow for Natural Single Rock Fractures

Accurate quantification of rock fracture aperture is important in investigating hydro-mechanical properties of rock fractures. Liquefied wood’s metal was used successfully to determine the spatial distribution of aperture with normal stress for natural single rock fractures. A modified 3D box counting method is developed and applied to quantify the spatial variation of rock fracture aperture with normal stress. New functional relations are developed for the following list: (a) Aperture fractal dimension versus effective normal stress; (b) Aperture fractal dimension versus mean aperture; (c) Fluid flow rate per unit hydraulic gradient per unit width versus mean aperture; (d) Fluid flow rate per unit hydraulic gradient per unit width versus aperture fractal dimension. The aperture fractal dimension was found to be a better parameter than mean aperture to correlate to fluid flow rate of natural single rock fractures. A highly refined variogram technique is used to investigate possible existence of aperture anisotropy. It was observed that the scale dependent fractal parameter, K _v, plays a more prominent role than the fractal dimension, D _a1d, on determining the anisotropy pattern of aperture data. A combined factor that represents both D _a1d and K _v, D _a1d × K _v, is suggested to capture the aperture anisotropy.     

广西贺州水岩坝矿田断裂构造分形特征分析

文章基于盒计维数法对水岩坝矿田断裂构造进行了系统的定量分析,结果显示矿田内不同尺度的断裂构造具有良好的自相似性。其中,断层的分维值为1.3475;矿田中部的烂头山矿床节理和手标本微裂隙的平均分维值分别为1.7841和1.3310;在子区分维高次趋势图上,分维高值区总体延伸趋势为北西向,与矿田断层构造总体展布趋势相吻合。平面上,烂头山矿床正常围岩微裂隙、蚀变围岩微裂隙和矿(化)体微裂隙的平均分维值分别为1.2097、1.3310和1.8622,表明矿化程度的高低与其微裂隙分维值的大小成正比;垂向上,烂头山矿床露头节理平均分维值随高程的升高而增大,显示烂头山矿床具有含钨石英脉型的"五层楼"模式特点。此外,水岩坝矿田内已知矿床(点)均位于分维高值区,说明分维值的高低与已知矿床(点)的分布具有耦合性。