^{1}

^{*}

^{1}

Geostatistics belongs to the wide class of statistical methods. It is used and applied to analyze and predict the values associated with spatial or spatio-temporal phenomena such as seismological events. Thus, in addition to its adaptability to perform spatial data analysis based on the principles of variography, many other methods more representative of the spatial distributions of the epicentres have been experimented. The unidirectional or isotropic fractal variogram associates the criterion of variance with the fractal dimension which is a better descriptive parameter of the spatial organization of events. The same analysis procedure carried out through a directional or azimuthal variogram introduces the context of a preferential direction of earthquake evolution. Moreover, b-value is a privileged seismological parameter that can be studied alone or in combination with other more advanced geostatistical analysis factors such as fractals, fractal dimensions and the anisotropic variogram. All these concepts were used as methodology for a protocol of analysis of the catalog of the Algeria seismicity including 1919 events for the period extending from 1673 to 2010. This catalog is provided by the compilation of partial catalogs synthesized by the CRAAG (Research Centre in Astronomy, Astrophysics and Geophysics).

Many natural phenomena are identifiable by measurements or recordings. Thus, the seismological data can be processed using stochastic approach or on the basis of existing uncertainties. More generally, statistical methods like Geostatistics are proving to be highly adapted tools to make it possible to provide an efficient assessment of the analysis parameters. The seismological data recorded in a catalog are identifiable through three types of dimensions: the space dimension (i.e. longitude and latitude) the magnitude and the time of occurrence (Wierner, 1996; Ogata & Katsura, 1993) . In the strictly geostatistical context, these data are assimilated to time series obeying the laws of distribution, clustering and scale invariance. The seismological data recorded in a catalog are identifiable through three types of dimensions, a dimension of space (longitude and latitude), magnitude and time of occurrence (Wiemer & Wiss, 2002) . On the purely geostatistical level, these data are assimilated to time series whose distribution laws, clustering mode and scale invariance status will be determined in the fractal case. These characteristics are elucidated only in relation to the parameters of analysis imposed by the objectives of the study. Various statistical techniques for studying local or regional seismicity have been applied to catalog data across a restricted geographic area or country. The frequency-magnitude distribution provided by the associated Gutenberg-Richter relationship (Gutenberg & Richter, 1942) and the fractal dimension concept (Hirata, 1989; Legrand, 2002) is considered as a convenient and objective approach in the seismicity analysis.

Conventional statistical analysis provides only a limited set of tools because it could be unaware certain no obvious aspects of the data to be analyzed (Cressie, 1993) . So the French mathematician Georges Matheron introduced geostatistics as a set of probabilistic methods designed to understand the structure of variables regionalized and distributed in a real space like seismicity data (Mathéron, 1970; Schaefer et al., 2014) .

Contrary to conventional statistics, methods, geostatistics takes care through adapted mechanisms of study and specific concepts describing the complexity and irregularity of the phenomena (Goovaerts, 1997; Mathéron, 1970) . Thus, the more recent statistical approaches are variogram (Olea, 2009) and directional or azimuthal variogram. Fractal structures, fractal variogram and the principle of spatial heterogeneity are associated with the classical concepts of variance, correlation, and random distribution. All these arguments require the implementation of specific methods combining geostatistical processing, variography and fractal analysis. The materials and methods included in this study aim to demystify the seismic catalog by highlighting the intrinsic parameters of earthquakes. The number of seismic events extracted from the catalog before declustering is around 1919 events reduced to 1879 after eliminating the non-significant smallest events and after the declustering step. Because of their direct connection with sismotectonic activity, the b-value extracted from the Gutenberg-Richter relationship, the variogram representation, the fractal dimensions and the interconnections between these parameters will constitute the framework of this study.

Mediterranean Sea belongs to an active and permanent geodynamic context resulting from the collision of the two African and Eurasian plates inducing a frequent but moderate seismicity. Indeed, the basins delimiting the extent of the Mediterranean form a perimeter zone of separation of these two geodynamic entities. This limit is characterized by a strong geological complexity induced by a sustainable and diversified geodynamics. This geodynamics includes alternating sequences of extensions and compressions. The extension phases lead to the appearance of new oceans whereas the compression phases led to the appearance of subduction zones and to the creation of collision chains. This geological complexity continues until recently in terms of deformation and seismotectonic activity. Overall, these phenomena are local or regional but exert an impact on the deformation field related to the convergence of African and Eurasian tectonic plates (

The coexistence of the compressive and distensile deformation regimes is the driving force of the NW direction convergence of the two plates. However, the speed of convergence is not uniform; it is estimated at between 3 and 8 mm/yr in the eastern Mediterranean and at 2 and 5 mm/yr in the western Mediterranean (

Algeria country is a part of the north-west of the African continent. It also belongs to the eastern Ibero-Maghrebian region. In the regional tectonic context, Algeria is located on the edge of the African and Eurasian plates, and it undergoes the effects induced by the convergence dynamics of these two plates (Kenzi, 1972) . The seismicity of Algeria is located on the northern fringe of the country formed by four morpho-structural domains that are the Tellian Atlas, the High Plateaus, the Saharan Atlas and the northern part of the Saharan Atlas platform (

The origin of this proven seismicity comes from the compressive regime to which the Tell is subjected, with a shortening of direction NNW-SSE. Earthquakes are generally generated by folds, faulty folds and especially by NE-SW oriented faults. The frequency and magnitude of the seismicity are important on the Tellian Atlas while this seismicity is ostensibly attenuated towards the southern Algeria. Moderate earthquakes are recorded in the Saharan Atlas and in the part of the Saharan platform. The frequency and magnitude of the seismicity are more important on the Tellian Atlas but this seismicity is disparate towards the southern Algeria region. However, historically, some moderate earthquakes have been recorded the Saharan platform.

Northern Algeria is characterized throughout its extent by seismicity caused by the confrontation of the African and Eurasian Plates. This seismicity involves the occurrence of weak to moderate earthquakes although sometimes violent earthquakes can occur. These earthquakes generate many disasters like the earthquakes of El-Asnam of 10/10/1980 (Benouar, 1994) and more recently the May 21, 2003 Boumerdes earthquake (Aitouche et al., 2012) . On a larger time scale, a large number of earthquakes were recorded, including those used in the catalog on which this study was based for a period ranging from 1673 to 2010 (

It is not easy to produce a homogeneous and reliable seismic catalog over such a long period of time (1673-2010). The catalog used in this study is available from the Centre for Research in Astronomy, Astrophysics and Geophysics (CRAAG―Algiers). It is the result of a synthesis and a compilation of partial catalogs each attached to a precise interval of time (Benouar, 1994) . In addition to the catalog provided by the CRAAG, several authors contributed to the final version of the seismic catalog as used in this study. So, a number of 1919 seismic events before the declustering sequence and covering the period 1673-2010 for the northern region of Algeria (longitude (Longitude 2.85˚W - 10˚E; latitude 31.28˚S - 37.80˚N) defines the core of the earthquake catalog. The main obstacle to the construction of the catalog lies in an objective conversion of intensities into surface magnitudes.

Northern Algeria is affected all along its Tellian part by a recurring seismicity but diversified in terms of periods of occurrence and magnitude (Bellalem et al., 2008) . Local or regional tectonics materialized by a network of active faults is the main source of this recurrent seismicity.

circled point has an anomaly that differs from all other points. It represents an area with high seismic potential. It is also worth noting the isolation of some epicentres which could indicate a sporadic seismicity due, perhaps, to a fault reactivation. Areas with low seismicity are either not visited by the isoseismal curves or by their relatively wide spacing.

The occurrence of any global seismic event is generally composed of two specific partial events including the major seismic event followed by a sequence of induced earthquakes also known as aftershocks. Under the assumption that the seismological processes generating these two phenomena are statistically Gaussian, they are then supposed to be independent of each other with a random spatial distribute. It then becomes necessary to remove these aftershocks classified as noise. The declustering of the catalog is the first step to perform before any processing operation. The declustering is conceived as a process of separating the catalog data into foreshocks, mainshocks and aftershock events. There are several algorithms for removing the aftershocks like Rosenberg algorithm (1975). Currently, the algorithm proposed by Gardner and Knopoff (1974) is widely used due to its availability and simplicity. This algorithm is based on an empirical identification of the aftershock sequence using time and space windows that vary as a function of the mainshock magnitude.

The Gardner-Knopoff algorithm is based on the windowing approach to identify the mainshock and the aftershocks. For each earthquake listed in the catalog with a magnitude M the associated seismic events are identified as aftershocks if they belong to a specified time interval and within a range of corresponding distance. Thereafter, the time windows will be reset according to the magnitude value of the largest magnitude in the sequence. An approximation of the size of the temporal window T and distance D according to the Gardner and Knopoff algorithm is given by the following equations

T ( Days ) = { 10 0.032 M + 2.7389 if M ≥ 6.5 10 0.5409 M − 0.547 else (1)

D ( Kms ) = 10 0.1238 M + 0.983 (2)

The Algerian seismicity catalog was declustered using Zmap software (Wierner, 2001) running under Matalb language (Wierner, 1996) . The Gardner-Knopoff’s algorithm is implemented with the possibility to vary the declustering parameters. This operation provided the map shown in

The north of Algeria is characterized practically throughout the Tellian part by intense seismicity induced by the collision of the African and Eurasian plates. This seismicity generally results in moderate to weak earthquakes, although severe strong earthquakes may occur. An abnormal high magnitude ( M w = 7.5 ) was cited in the catalog.

epicentres of all earthquakes that occurred from the year 1673 to the year 2010. They illustrate a scattering of the epicentres along the entire Tellian band.

The cumulative number vs. time curve as shown in the figure (small squares) depicts a regular increase in the number of earthquakes for the period considered (1673-2010) with two remarks, however. The first relates to the presence of irregularities comparable to fluctuations in the number of earthquakes in the interval 1850-1950 (figure (small squares) probably due to a regional relaxation of the stress field. This state may be correlated in the same time with the frequency histogram (

For example, the destructive earthquake of magnitude M w = 6.8 which hit the balneal town earthquake of Boumerdes (east of the capital Algiers) on May 21, 2003 is ranked among the strongest seismic events that occur in this region. It was a shallow quake (10 km in depth) which caused important material and human damage in the epicentral region. Shallow earthquakes are caused by fault breaks or continental inter-plate movements. This may explain the large spread and the greater damage caused by some quakes that occurred in the Algerian Tellian region.

The relationship between magnitude and frequency of occurrence of past earthquakes has been described by the Gutenberg-Richter power law (Gutenberg & Richter, 1942) as given by n the following equation

log 10 N ( > M ) = a − b M (3)

where N ( > M ) indicates the cumulative number of earthquakes that have a magnitude greater than M while a and b are specific real positive coefficients which vary in a fixed time and space window (Crownover, 1995) . In a seismotectonic context, the constant a is related to the level of seismicity and to the seismic activity while the parameter b also known as b-value describes the number of the large seismic events that occurred or indirectly the level of the accumulated stress. b-value may be also used for interpretation of geological structures and to understand the faulting systems. The changes of b-value are an indicator of the level of stress in the area of study. Globally the b-value is close to 1.0, but experimentally it can vary from 0.5 to 1.5. Deviations of this parameter are then frequently observed and a mapping of the b-value variations can help with a better seismological investigation. Graphically, the b-value is calculated from the slope of the line representing the Equation (3) as showed in

Since b-value is related to the size of the seismic event distribution, it can thus define a measure of the frequency of weak or strong earthquakes that occurred in the region under study (Wackermagel, 2003; Mc Nutt et al., 2004) . Moreover, b-value is also an indicator of the tectonic activity and consequently of the field of stress.

The map shown in

The concept of fractal analysis has become a widely used tool in the descriptive, studies and understanding of some geophysical problems (Turcotte, 1981) . This theory essentially encompasses so-called fractal geometry and the fractal scaling property. Fractal geometry was initially formulated by the Franco-American mathematician B. Mandelbrot in terms of scale invariance and self-similarity observed in certain objects or phenomena (Mandelbrot, 1982; Crownover, 1995; Feder, 1989) . In addition to these two characteristics, we must add the fundamental parameter of fractal dimension which will play a key role in this study.

Fractal geometry provides a tool for investing at different scales of fractal objects and processes. It also allows the observation of the constitutive elements of the irregular shapes which can not be revealed by the Euclidean conventional geometry. The major factor of the fractal theory is undoubtedly the fractal dimension D F as well as its many variants including the capacity dimension or box counting fractal dimension and the correlation dimension D C R (Klinkenberg, 1994) . A self-similar object is identical or approximately similar to any subset of itself which can be continuously subdivided into ever smaller elements whose structure is a scale-copy of the whole. A seismic regime can be assimilated to a set of earthquake swarms distributed in space and time. The grouping of these swarms gives a representation similar to an aggregate whose size and distribution follow the relationship

N ( r > R ) = C R − D F (4)

where N ( r > R ) is the full number of objects whose size r is greater than a reference size R and D F is the fractal dimension. Equation (4) shows that the value of the fractal dimension D F depends on both the shape of the individual object and the spatial scales that govern the events distribution. By extension, variations in the fractal dimension could be an indicator of the spatial heterogeneity of the distribution. A high value of the fractal dimension implies an increase in the number of objects of small size and hence of the heterogeneity. The fractal dimension thus becomes an indicator of heterogeneity of the spatial distribution of events. Several other approaches to the definition of the fractal dimension have been proposed by various authors.

Seismicity through its different analysis parameters highlights fractal structures in space, time and magnitude distribution as described by the fractal dimension of the epicentre distribution. Among these parameters the b-value, deduced from the Gutenberg-Richter relationship, plays a key role in the analysis of seismological data. Indeed, b-value has the advantage of being correlated with several other parameters of seismicity of which the fractal dimension (Hirata, 1989) . One of these relations was that proposed by Aki (Aki & Richards, 1980) and formulated as follows

D F = 3 b c (5)

where c defines a constant closely related to the size of the seismic events:

・ For small events: c = 1 so D F = 3 b .

・ For intermediate events: c = 1.5 so D F = 2 b .

・ For large events: c = 2 so D F = 3 b / 2 .

Generally an, with some exceptions, Algeria’s seismicity catalog shows intermediate-size earthquakes. The second case cited above may be accepted i.e. D = 2 b , the value of b being that provided by the Gutenberg-Richter law as shown in

There are various methods to compute the fractal dimension of a discrete space or time series apparently randomly distributed. The most popular and the first one used among these methods is the so-called “box counting” or capacity fractal dimension. As introduced by B. Mandelbrot (Mandelbrot, 1982) , it consists in counting the number of small squares (in two dimensions) or cubed (in three dimensions) required to cover the whole of the fractal object under consideration. Therefore, this number is compared with the used box size. The limit of the ratio when the size of a square or cube tends towards 0 defines the box counting dimension noted D C B (Equation (6))

D C B = lim r → 0 log [ N ( r ) ] log ( 1 r ) (6)

where N ( r ) is the number of no, empty boxes (squares or cubes of size r).

The first disadvantage encountered when applying the capacity dimension is that the latter is purely geometric and depends on the density of the number of points of the spatial or temporal series. In the case of a spatial distribution of seismic events, this distribution can be assimilated to a dynamic process. The capacity dimension does not take this dynamic into account, in particular the distance that can separate two samples. To correct this drawback of the box-counting dimension, another empirical dimension called fractal correlation dimension has been proposed. The assessment of this dimension is based on the set of sample points ( Z 1 , Z 2 , ... , Z n ) belonging to the considered fractal object. For a chosen radius r “sufficiently small”, we compute the proportion of pairs ( Z i , Z j ) denoted C ( r ) such as their inter-distance is less than r. This procedure of computing is repeated for many values of r, say ( r 1 , r 2 , .... , r n ) . These values must be chosen sufficiently small. ln [ C ( r ) ] is then regressed by − ln ( r ) . The slope of the regression line is the D C R correlation dimension, defined as

D C R = lim r → 0 [ − ln ( C ( r ) ) ln ( r ) ] (7)

Applied to the data of the seismological catalog this procedure of calculation of the correlation fractal dimension leads to the value 1.4 +/− 0.03 i.e. 1.61 < D C R < 1.67 . It was previously calculated the fractal dimension using the b-value. The result was D F = 1.546 . The two fractal dimensions are therefore substantially equal (

The statistical analysis by application of the so-called statistical variography method consists mainly of the definition of an experimental variogram calculated from the data and then of the selection of a theoretical variogram model adapted to the data to be processed (Yarus & Chambers, 1994) . The variogram belongs to the class of descriptive statistics techniques (Amorese et al., 2010) . It graphically illustrates and characterizes the spatial continuity of the data set i.e. its roughness. The experimental variogram is calculated by averaging one-half the squared differences of the variable under consideration on all pairs of observations with the specified value separation distance and direction (Equation (8))

γ ( u , h ) = 1 2 N ( a , h ) ∑ i = 1 N ( a , h ) [ Z ( i ) − Z ( i + h ) ] 2 (8)

where γ ( u , h ) is the semi-variogram estimated for the distance h, Z ( i ) and Z ( i + h ) are the samples of the variable at the points i and i + h , N ( u , h ) defines the number of pairs of data observed in the direction u separated each other by distance lag h. When the direction u is not specified, the following model is then obtained

N ( h ) = 1 2 N ( h ) ∑ i = 1 N ( h ) [ Z ( i ) − Z ( i + h ) ] 2 (9)

According to the data, the theoretical variogram is selected from a set of other several variograms so as to coincide as much as with the experimental variogram (

The seismological data of this study require the choice of an exponential variogram (

γ ( h ) = { C 0 + C [ 1 − exp ( − 3 h a ) ] if h ≠ 0 0 if h = 0 (10)

where C o ≥ 0 , C > 0 and a > 0 .

The exponential variogram model reaches its sill fairly rapidly and asymptotically. For this model, the range is the value of the distance lag where the variogram reaches 95% of the sill. An example of an exponential variogram model is illustrated by

A variogram is characterized by at least three descriptive parameters cited as follows: the nugget C 0 , the sill identified at the ordinate C + C 0 and the range a which represents the lag distance in which the plateau reaches its horizontal shape (Kouadri et al., 2012; Barnes, 1991) .

A non-objective analysis of each of these parameters may the quality of the empirical variogram interpretation. Concretely each parameter defines a specific characteristic of the variogram. The nugget C 0 or intercept of the variogram could be attributed to measurement errors or a scale spatial variation (

zero distance and the samples will be the same. However, some variables may change suddenly. A difference over a very short distance is then observable.

When analyzing a semivariogram, it appears at a certain distance, a stabilization of the model. The distance at which the model begins to flatten is called range noted a. If sample locations are separated by a distance smaller than the range, then the samples are spatially correlated. In the opposite case, the samples will be spatially uncorrelated. The value of the semivariogram corresponding to thee range is called sill expressed by C + C 0 . Mathematically, the sill is defined

as e lim h → + ∞ γ ( h ) of the random field of data. In another words, the sill is simply variance of the variables.

Depending on the statistical objectives and spatial characteristics of the data under consideration a more expressive variogram may be required. For a uniform distribution of the seismic events, a unidirectional or isotropic variogram may be sufficient to describe and to understand that distribution. However in the case of random or irregular distribution, a directional or anisotropy variogram may be used to assess the degree of anisotropy and to extract the major direction of spatial continuity. Anisotropy can be identified by the covariance function or by the nature of the experimental variogram model used for the geostatistical analysis. In practice, a suspected spatial anisotropy can be carried out by the high variability displayed by the experimental variogram which becomes thus disturbed. A directional variogram has the property of varying in structure and fractal dimension by a preliminary choice of a direction of analysis. This behaviour may describe the anisotropic tendency of the data. It also measures the possible dissimilarities expressed as a function of the separation distance on which the variogram is conceived. Indeed, the variable under consideration can change spatially more rapidly in a direction rather than in another.

Starting from the declustered seismic catalog, the first step is to draw the so-called experimental variogram model. The shape of the graph of the experimental variogram is then compared with a graph of the theoretical variograms in order to extract the mathematical model. As shown in

The analysis of

・ Nugget: C 0 = 0.13160 .

・ Sill: C 0 + C = 0.42820 .

・ Range (estimated): a ≈ 0.64 .

The graph of the exponential semi-variogram describing the spatial distribution of the data of the seismicity catalog of Algeria thus has an intercept represented by the non-zero nugget C 0 . This is a consequence of the reliability to be given to the data of the catalog. Indeed, the conversion of seismic data from the intensity scale to the scale of magnitudes can be a source of errors. Let us add to this probably measurement or recording errors.

The calculated sill C + C 0 refers to the variance of the set of magnitudes. The variance scale shown in

The variogram cloud is a graph that highlights the average dissimilarity between two observations as a function of the distance separating them spatially. The variogram cloud is applied in the analysis of spatial data to better identify the dissimilarity of the couples and to compensate the deficiencies analysis details not provided by the classical variogram. The variogram cloud plot individualizes the contribution of each data pair to the global variogram (Klinkenberg, 1994; Wackermagel, 2003) . Compared to the classical variogram, it allows to have a wider view of the spatial variation and to identify the unusual behavior of certain points. The latter may be aberrant points that are in fact unnecessary in geostatistical analysis.

The construction of a cloud variogram is based on the calculation of the quadratic difference of the samples relatively to lag distance. The possible dissimilarities are then plotted against separation distance of the pairs of samples. The resulting graph forms the variogram cloud for a fixed class of distance.

Semivariogram based analysis must not only measure the dissimilarity of the data against distance but also investigate anisotropic data. When using a conventional semivariogram the variable is studied without taking into account the

directions of investigation. In the isotropic variogram model only the modulus of the vector is used (equation). Directional variography introduces another parameter which is the direction of azimuthal investigation (

This form of representation of an anisotropic semivariogram also called 2D-variogram map is simply a visual image of semi-variance for various directions. The two coordinate axes contain the separation distances or lags. Thus the abscissa axis corresponds to the E-W direction lags and the Y axis to the N-S direction. The center of the image corresponds to the origin of the semivariogram for each direction. The color bar inserted next to the anisotropic semivariance surface denotes a scale of values for this semivariance. The result of the application of the 2D-variogram map on the data of the seismic catalog of Algeria is illustrated by

One of the simplest semivariogram models is the power model given by the following equation

γ ( h ) = a h b (11)

Coefficients a and b are real numbers and they represent respectively a quantification of the rate of variation and a scaling factor. The scaling factor a plays a key role in the shape of the variogram depending on whether 0 < a < 1 or 1 < a < 2 (

γ ( h ) = 1 2 a h 2 H (12)

If we carry out a contraction of scale on the variable h by taking h k H = h ′ , the corresponding variogram is given by

γ ( h ′ ) = 1 2 a ( h k ) 2 H (13)

Relationship (13) is precisely of the same form as that of the initial variogram. A self-similarity state is thus shown in the case of the variogram. A junction between variogram and fractals is therefore established.

Two kinds of geometric dimensions can characterize the spatial distribution of data such as those provided by a seismic catalog. These are the natural or topological dimension D T and the fractal dimension D F . Generally, the fractal dimension is comprised between two consecutive topological dimensions, it is then fractional. Fractal objects are characterized by their apparent complexity,

irregularity and random fluctuations. The “fractality” of an object is an intuitive perception of its irregularity and roughness assessed by the fractal dimension. The more an object is irregular, the higher its fractal dimensions. The relationship between the topological dimension D T and the fractal dimension D F is provided by the physical process called the Brownian Fractional Movement (fBm), the study of which is outside the scope of this study. Equation (13) is a model of fBm indexed by the parameter H such that 0 < H < 1 . The three quantities H, D T and D F are unified by the following relationship

D F = D T + 1 − H (14)

For D T = 1 we obtain

H = 2 − D F (15)

Equation (14) is the basis of the geostatistical analysis by fractal variogram of the seismicity catalog data. Equation (13) can be written as

γ ( h ) = 1 2 a h 2 ( 2 − D F ) (16)

Bi-logarithmic coordinates applied on Equation (16) i.e. log [ γ ( h ) ] against log ( h ) provides the equation of a straight line pe m = 2 ( 2 − D F ) . The fractal dimension D F is then easily calculable:

D F = 2 − m 2 (17)

This procedure will allow the analysis of the spatial anisotropy of the data of the seismic catalog jointly by using the directional semivariogram and the fractal dimension.

Four a directions of analysis of the spatial anisotropy were chosen. These directions correspond to those most commonly used. By applying the previous calculation procedure, the fractal dimensions btained are summarized by

The most popular scaling power law for earthquakes i.e. the Gutenberg-Richter’s relation is the basis of this study. It describes for a given seismogenic area the distribution of earthquake occurrence. Two important seismological parameters are provided by this law namely the a-value and the b-value the a-value and especially the b-value which depends on the tectonic regime induced by the stress field. The Gutenberg-Richter’s b-value was the first seismological parameter calculated and interpreted.

α | 0˚ | 45˚ | 90˚ | 135˚ |
---|---|---|---|---|

D_{F} | 1.845 | 1.911 | 2 | 1.942 |

The b-value is of great importance, because the concept of geometrical self-similarity it has been shown that b-value can be directly related to the fractal dimension. The graphically estimated value of the b-value for the seismicity of Algeria is b = 0.773 ± 0.06 . The analysis of this global value as well as the distribution map of this seismic parameter reveals a moderate, frequent and uniformly distributed seismic activity throughout the north of Algeria.

The concept of fractal and related concepts of self-similarity and scale invariance have been introduced to understand the spatial organization of cluster epicentres that obeys an underlying order known as fractal structure.

The value of the fractal dimension is considered as an index of the complexity that can prevail in a time series such as a data of magnitude. The moderate seismicity of northern of Algeria has suggested using the relationship D F = 2 b i.e. D F = 1.546 . This value of the fractal dimension is intermediate between the two topological dimensions 1 and 2, first confirming average data complexity and non-uniform epicentre distribution. However, the fractal dimension D F does not take into consideration the clustering mode in the spatial distribution of epicentres. Another type of fractal dimension is then proposed: the correlation dimension denoted D C R . Graphically the value D C R = 1.64 ± 0.03 was obtained. Fractal dimension D F and correlation dimension D C R are substantially equivalent. They are substantially equivalent. The clustering mode of earthquakes in northern Algeria thus presents a not very complex distribution.

Basically, the mean and variance parameters are of great importance in the geostatistical analysis of the time series because they are respectively indicators of the homogeneous and isotropic character of the data. With respect to homogeneity, the variographic approach was used in this study because it provided an experimental variogram of exponential type with a non-zero nugget C 0 which gives some reliability to the data despite the likely uncertainties of measurement or recording. The calculated sill C + C 0 means that the data contained in the seismic catalog tend rather towards a spatial heterogeneity.

The study of the spatial anisotropy of the data was carried out by means of the variogram which showed an NW direction of anisotropy i.e. the same direction of the compressive geodynamic regime previously discussed.

The junction between semivariogram and fractal dimension is established through the directional anisotropic variogram. For a chosen range of directions, the fractal dimension was calculated; thus for the direction α = 45 ˚ it was obtained the fractal dimension D F = 1.942 close to the topological dimension D T = 2 . This means that the density of the earthquakes is very high in this direction, which therefore has a high seismic activity.

Despite the constraints imposed by the nature and quality of the seismological data, the catalog of seismicity of Algeria was proved to be a fruitful source of information on the seismological activity. The extraction and processing of this information requires the development of an analysis protocol. Statistical analysis in general and the application of geostatistical methods in particular are tools adapted to these objectives. Despite their relevance, current seismological parameters such as b-value can not suffice to carry out an extended survey of seismological events. They lack spatial dimension that is to say a description of the distribution of epicentres. The concepts of conventional or isotropic variogram and its variants such as the variogram cloud, the variogram surface and the directional variogram are additional analytical tools recommended for catalog data analysis. Fractal theory and fractal geostatistics are among the most recent orientations for the analysis of seismological data in the catalog. The first one provided the concept of fractal variogram across the fractal dimension which thus becomes a seismic parameter suitable for the study of the anisotropy of spatial distribution of epicentres for the study area.

The authors declare no conflicts of interest regarding the publication of this paper.

Aitouche, M.-A., & Djeddi, M. (2018). Conventional and Fractal Variogram Based on Time―Space Analysis of Seismicity Distribution―Case Study: Algeria Seismicity (1673-2010). Journal of Geoscience and Environment Protection, 6, 147-172. https://doi.org/10.4236/gep.2018.611012