Layout 6 The effects of zoning methods on forecasting of the next earthquake occurrences, according to semi-Markov models Ramin Sadeghian1,* 1 Bu Ali Sina University, Hamedan, Iran ANNALS OF GEOPHYSICS, 53, 5-6, 2010; doi: 10.4401/ag-4558 ABSTRACT Earthquakes are natural phenomena that can be viewed in three dimensions: time, space and magnitude. Earthquakes can be investigated not only physically, but also mathematically. In this study, semi-Markov models are applied, which can be considered as useful methods to analyze and forecast the occurrence of future earthquakes based on previous earthquake data. In the present study, the target region, Iran, is divided into zones, and each zone is examined as one of the semi-Markov model states. Several methods to determine the levels of forecasting error are then introduced and applied to the target area. The results of the application of these semi-Markov models to investigate and forecast the occurrence of future earthquakes are obtained and analyzed mathematically. A new zoning method is developed and compared with that of Karakaisis, through the proposed forecasting method. Moreover, the effects of the type of zoning and the number of zones on the forecasting error of the next earthquake occurrences are investigated using several algorithms. Editor's note. The present paper revisits and completes the paper by Sadeghian and Jalali-Naini [2008b] cited by the author in the References. Introduction Earthquakes can be examined both mathematically and physically [Sadeghian 2007]. Stochastic processes represent one of the branches of mathematics that can be applied to probabilistic investigations of these phenomena, and one of the stochastic process models that has been frequently applied over recent years is the Markov model, and especially the semi-Markov model. Semi-Markov models belong to those that can model events that have a relationship with previous events [Jalali-Naini 1997, Minh 2001, Rice 2007, Sadeghian and Jalali-Naini 2008a]. For analyzing of the three temporal, spatial, magnitudinal dimensions of earthquakes, Markov models, including semi-Markov models, can be applied [Patwardhan et al. 1980]. Markov models have already been applied for earthquake occurrence analysis [Di Luccio et al. 1997, Console 2001, Console et al. 2002]; moreover, in this regard, a few studies have recommended the use of semi-Markov models [Patwardhan et al. 1980, Altinok and Kolcak 1999, Sadeghian and Jalali-Naini 2008b, Jalali-Naini and Sadeghian 2009]. One of the advantages of the model applied in this study is that it can forecast this triad of dimensions simultaneously, which has not been provided in previous studies [Kelleher et al. 1973, Nanjo et al. 2006, Papazachos 1992, Rikitake 1976]. For every semi-Markov model, and the inspection of each event, several states have to be defined. In this study, the Iran area was chosen as the target for study. Seismologically, this area was divided into various parts by Karakaisis [1994] (K94), and also in the present study. Each part is called a zone, and each zone can be considered as a state of a semi-Markov model. A similar classification has to be used for the magnitudes of the earthquake occurrences too. Using the proposed semi-Markov models, later earthquakes can be forecast probabilistically. Then, to demonstrate whether the forecasting results obtained are mathematically valid or not, several methods are introduced to calculate the differences between the probabilities of the occurrences of actual earthquakes and those forecast, which from now on is called the forecasting error. Nava et al. [2005] presented another method for validation of their model. The forecasting error of the algorithms presented is obtained, and the results are investigated and analyzed. Moreover, a deterministic forecasting algorithm, based on the provided probabilistic one, is presented, and results related to deterministic forecasting are also mentioned. Modeling Applying semi-Markov models in the forecasting of the dimensions of earthquakes was completely examined by Altinok and Kolcak [1999] and Jalali-Naini [1997]. One of the most important elements of the semi-Markov models is the interval transition probability matrix. The probability of a transition from state i to state j in the interval [0,n] requires the process to make at least one transition during that Article history Received February 2, 2010; accepted October 4, 2010. Subject classification: Semi-Markov model, Transition probability matrix, Probabilistic forecasting, Deterministic forecasting, Forecasting error, K94 zoning. 83 interval. The interval transition probability matrix can be determined in matrix form as follows: (1) where n is the number of time intervals, WG (n) is a diagonal matrix where the i-th element is equal to , and is the probability that the waiting time in state i is greater than n, and G is a transition matrix, so that Gij is the probability that if the last step is in state i, the next step will be in state j. In other words, Gij is the probability of transition from state i to state j. Also T(m) is the holding time matrix, which is defined as follows: The interval transition probability F(n) is obtained by a recursive procedure. Since T(0) is equal to zero, F(n) is just obtained for the interval 1≤m≤n. In the case n=0, F (n) is equal to the Kronecker Delta or identity matrix, which is defined as follows: Also, C (m) is the core matrix, which is defined as follows: Where N is the total number of states in the system [Jalali-Naini and Sadeghian 2009]. For the earthquake phenomenon, F (n) can be used to study earthquake hazards and to evaluate their risk. The two interval transition probability matrices that consist of FR(k) n (for region-to-region transitions) and FM(k) n (for magnitude-to-magnitude tran- sitions) can be determined through Equation (1). If the last earthquake occurred in region r0 with a magnitude m0, the matrix of probabilistic forecasting after k; n time periods (i.e. n) is obtained by the following formula: (2) where r is the number of zones of the assumed area considered, m is the number of classes for all magnitudes considered, and k specifies the number of time periods. Therefore, by computing , all of the three dimensions of any earthquake can be taken into account and specified. As in any region, any earthquake with any magnitude can occur, the region and magnitude can be considered as two independent random variables. Consequently, based on probability theory, if Xand Yare two independent random variables, f(x,y)=f(x).f(y), where f is a probability density function [Jalali-Naini 1997, Rice 2007]. The same relationship is also demonstrated by Equation (2). Also is the probability that an earthquake will occur in the region ri with the magnitude mi after the time period k. Figure 1 illustrates this more clearly. In this way, forecasting the dimensions of following earthquakes is possible, through determining the probability forecasting matrices (i.e. ). In this study, the probabilistic forecasting matrix and deterministic forecasting matrix are demonstrated by the and FRM symbols respectively. For model validation, assume that is a probability forecasting matrix, where its elements are the estimated probabilities of the next earthquake occurrences according to our proposed model, during the next k time periods, and FMR (k) is a deterministic forecasting matrix defined as follows: There are three common ways for calculating the forecasting error: mean square error (MSE), mean absolute deviation (MAD) [Rice 2007], and mean absolute percent error (MAPE) [Haj Shirmohammadi 2003]. Let Fi be the forecasting value for period i, and Di be the real data for period i, then the forecasting error (Ei) is calculated by Ei = Di − Fi. These three attributes for comparing forecast precision can be used according to the following relationships: USING SEMI-MARKOV MODELS FOR ZONING 84 ! ! ! ! Figure 1. The ways that the matrix is determined [from Sadeghian and Jalali-Naini 2008b]. ( ) ( ) ( ) ( ) ( ) ( ) ( ); 0, 1, 2, ... F n W n G T m F n m W n C m F n m n 0 0 G m n G m n = + = = + = 7 - - = = / / ( )w ni G ( ) Pr ; 1, ..., .T m t m m nij ij= = =" , (0) 1 0 F i j i jij = = ! ' ( ) ( ); , 1, ..., , 1, ..., .C m m i j N m nG Tij ij ij= = = 1, ...,k =6 1, ...,k =6 1, ...,k =6 ( ) 1, ...,FRM k k =6t 1, ..., ; , ..., ; , ..., FRM ( ) ( ) ( ) k r j m k n k FR k FM k 1 1 r mr m m m0i ji j j =# 6= = = t ( ) 1, ..., ; 1, ..., ; 1, ...,FRM k i r j m k nr mi j = = =6 t ( ) 1, ..., ; 1, ..., ; 1, ...,FRM k i r j m k nr mi j = = =6 t ( ) 1, ...,FRM k k n=6t FRMt ( )FRM kt with a magnitude in th time periodm k Otherwise ( ) ( )FRM k FRM k i jij r mi j =/ If an earthquake occurs in region r1 0 * FRMt MSE E 2i i= N / MAD Ei i= ; ; N / 85 where N is the number of historical data for the forecasting. Note that the MSE and MAD are not scale less, while the MAPE is scale less. Therefore MAPE is more suitable than MSE and MAD to compare several forecasting errors together. Here, two fundamental questions are posed: The first question relates to the probable estimation of earthquake occurrences. In other words, can forecast the next earthquake occurrences successfully? The next question is to find out whether can be used for forecasting the next earthquake occurrences not only probabilistically, but also deterministically. The response to the first question can be found in methods 1, 2 and 3, and the response to the second question is explained in method 3, as: Method 1: MSE and MAD [Rice 2007] Method 2: MAPE [Haj Shirmohammadi 2003] Method 3: An innovative plan that can be used for determining probabilistic forecasting, deterministic forecasting and their forecasting errors (in this study, this is called the zero and one (0-1) method technically). In methods 1 and 2, if the total number of the past data is equal to n, the first n1 data will be used for forecasting the next n2 data, which will be used to determine the forecasting error (i.e. n=n1+n2). With all of the n data having already occurred and being available, the forecasting errors can be calculated. In method 3, the first n1 data are used for forecasting the next n2 data, which is applied for the making of a pattern, and the first n1+n2 data are used for forecasting the next n3 data, which is applied to the determination of the forecasting error (i.e. n=n1+n2+n3). The whole ndata have already occurred and are available; therefore, the forecasting errors can be calculated. In these methods, calculation of the forecasting error can be performed in two forms: by using successive data or by using random data within the set of whole data. In the first form, the first n1 data are used for forecasting the next n2 data, (i.e. n=n1+n2), while in the second form, the n2 data within the main n data are eliminated randomly and then they will be forecast by their past data. Method 1 In this method, the relationships between MSEand MAD are used. For this goal, an algorithm is presented for MSE, which can be used with some changes for both MAD and method 2. Algorithm I: This algorithm can be used to determine the forecasting error by MSE as follows: Step 0: Begin Step 1: n2 (0)=0 Step 2: i=0 Step 3: Use the first n1+n2 (i) data to determine , FRM (n1 is the number of the first data that can be used for forecasting the next data, and n2 (i) is the number of data that occurred during i time periods after the first n1 data). Step 4: Determine , FRM Step 5: (3) Step 6: i=i+1 Step 7: If then go to step 3; otherwise go to step 8 (n is the number of total data) Step 8: Step 9: End. For calculating the forecasting error by MAD, use algorithm (I); just be careful to substitute steps 5 and 8 in algorithm (I) by considering the following steps: Step 5: (4) Step 8: For each forecast, if the past data are more complete, the forecast results are more accurate. Therefore applying after gives more accurate information than . Accordingly, only (1) and FRM (1) are used in Equations (3) and (4), and in the remaining equations. One of the disadvantages of this method is that neither MSEnor MADhave any mathematical interpretation, nor is there any benchmark for their degree of goodness. In other words, if an amount is obtained for them, there is no common mathematical interpretation for every one. For example, by MSE=2, regarding the forecasting error, it is not clear whether 2 is a suitable forecasting error or not. Method 2 In this method, a similar equation to algorithm (I) is used, yet steps 5 and 8 in algorithm (I) are substituted by the following intended steps: Step 5: USING SEMI-MARKOV MODELS FOR ZONING 100MAPE D E i i i= # ; ; N / FRMt FRMt FRMt FRMt ( 1) ( (1) (1)) MSE i FRM FRM 1 2 1 ij j m ij i r + = - == r m# t// ( )n n j n1 2 1j i + 1 = / ( ) MSE i MSE k 1k i = = / ( 1) (1) (1) MAD i FRM FRM 11 ij j m ij i r + = ; ;- == r m# t// ( ) i k MAD MAD 1k i = = / ( )FRM k1t( )FRM k2t ( )FRM k k1 2+t FRMt ( 1) (1) (1) (1) 100MAPE i FRM FRM FRM 11 ij ij ij j m i r + = # ; ;- == r.m t // (5) Step 8: where, as indicated above, MAPE is the mean absolute percentage error. MAPE is a useful equation, as it specifies the error percentage, which is both comprehensible and interpretable for every one. Therefore, this method will be concentrated upon in this study. Method 3 Definition: the i-th maximum order The i-th element in a descending sorted list, where none of the elements are equal to each other, is named the i-th maximum order. This method is an innovative plan that can be used for two goals: a) Determining the forecasting error b) Deterministic forecasting of earthquake occurrences. a) Determining the forecasting error In this section, an algorithm to determine the forecasting error is presented. This algorithm has two parts. The first part of the algorithm is devoted to making a pattern, and the second part is assigned to determine the forecasting error. This algorithm consists of the following steps: Algorithm II: Step 0: Begin Step 1: n2 (0)=0 Step 2: i=0 Step 3: Use first n1 +n2 (i) data for determining , FRM Step 4: Determine , FRM Step 5: Step 6: i= i+1 Step 7: If then go to step 3; otherwise k1 = i and go to step 8 (n3 is the number of data used to determine the forecasting error) Step 8: For i=1 to k1 do Step 10: i=1 Step 11: Mij={i-th maximum order in } Step 12: ([a] obtains the greatest integer number smaller than the real number a) Step 13: All of the elements greater than or equal to 1 in , are replaced by 1 Step 14: Step 15: Step 16: i= i+1 Step 17: If i-th maximum order is available in then go to step 11 otherwise go to step 18 Step 18: If MAPE(t+1) is the first element greater than the MAPE obtained in method 2, the i-th maximum order is the most correct one to be considered as a benchmark to calculate the errors. Step 19: n3 (0)=0 Step 20: i=0 Step 21: Use the first n1 +n2 +n3 (i) data to determine , FRM (n1 is the number of the first data which can be used to forecast the next n2 data, which can be used in making a pattern, and n3 (i) is the number of data that occurred during i time periods after the first n1 +n2 data) Step 22: Determine , FRM Step 23: Step 24: i= i+1 Step 25: If then go to step 21; other- wise k2 = i and go to step 26 Step 26: Mj={t-th maximum order in } Step 27: ([a] obtains the greatest integer number smaller than the real number a) Step 28: All elements greater than or equal to 1 in are replaced by 1 Step 29: USING SEMI-MARKOV MODELS FOR ZONING 86 ( ) MAPE i MAPE k 1k i = = / FRMt FRMt 1 ( 1) (1), 1 ( 1) (1)FRM i FRM FRM i FRM+ = + =t t t t ( )n n j n n1 2 1j i 3+ 1 - = / 2 ( ) ( is a zero matrix)FRM i O O=t 2 ( ) 1 ( ) 1 ( ) 2 ( ) 100MAPE i FRM i FRM i FRM i 11 ij ij ij j m i r = # ; ;- == r.m t // (0) 2 ( ) MAPE k MAPE k 1 1k k1 = = / Step 9: 1 ( ), 1, ...,FRM j j k1=6t ( ) 1 ( ) , 1, ...,FRM j M FRM j j k2 1 ij = =6t t; E ( ) 1, ...,FRM j j k2 1=6t 2 ( )MAPE j = 1 ( ) 1 ( ) 2 ( ) 100, 1, ..., FRM j FRM j FRM j j k 11 1 pq pq pq q m p r = =# 6 ; ;- == r.m t // ( ) 2 ( ) MAPE i k MAPE j 1 1j k1 = = / 1 ( ),FRM j j =6t 1, ..., k1 FRMt FRMt 1 ( 1) (1), 1 ( 1) (1)FRM i FRM F M i FRMR+ = + =t t ( )n n n j n1 1j i 2 3+ + 1 = / 1 ( ), 1, ...,FRM j j k2=6t 2 ( ) 1 ( ) , 1, ...,FRM j M FRM j j k j 2= =6t t; E 2 ( ),FRM jt 1, ...,j k2=6 1 ( ) 1 ( ) 2 ( ) 100 , 1, ..., FRM FRM FRM j k j j j 11 2 pq mr pq pq qp = =# 6 ; ;- == r.m t // 2 ( )MAPE j = 87 Step 30: Step 31: End. Note that in algorithm (II), MAPE can be substituted by either MSE or MAD, but because of the aforementioned reasons, MAPE is preferred. b) Deterministic forecasting of earthquake occurrences Section (a) of method 3, along with the probabilistic forecasting matrixes, such as , can be used for determin- istic forecasting of earthquake occurrences. Considering algorithm (II), and using the following steps, this makes it possible to forecast the earthquakes occurring in the next several time periods, deterministically (i.e. one and zero; so that one specifies the occurrence and zero specifies the lack of occurrence of any earthquake). Algorithm III: Step 0: Begin Step 1: Use the past ndata to determine (n is the number of total data) Step 2: Determine Step 3: (k is the number of predictable time periods in the future) Step 4: Mj={t-th maximum order in } Step 5: ([a] obtains the greatest integer number smaller than the real number a) Step 6: Replace any element(s) greater than or equal to 1 in by 1 and name the resulting matrices Step 7: If FRMDrm(j)=1 then an earthquake in region rwith a magnitude m in j-th time period will occur, otherwise any earthquake in region rwith a magnitude m in j-th time period will not occur Step 8: End. Application In this section, the earthquake occurrences are forecast probabilistically and deterministically using the actual data collected from the Iran area. Then, two zoning methods of the Iran area are introduced and considered, consisting of zoning by K94 and a proposed zoning in the present study, and the forecasting error for these zoning methods is determined. The Iran area was selected as the area of investigation. This area is bounded by longitudes 44.23˚E, 63.33˚E and latitudes 25.05˚N, 39.78˚N. The data were collected from the United States Geological Survey website (http://neic.usgs.gov/neis/epic/epic.html). After filtering and removing the unsuitable data according to Table 1 [Gardner and Knopoff 1974], 3,179 data related to earthquakes that occurred during 1973-2007 are used. The maximum time interval between the times of the earthquake occurrences is 45 days, so by considering each 10 days as one time unit, forecasting during the next time units in each time will be possible. Also, each zone in each zoning method is considered as a state in the region-to- region transitions of a semi-Markov model, and the magnitude of total occurrences is classified into 5 classes for each of the two zoning methods, and each class of magnitudes is considered as a state in the magnitude-to- magnitude transitions of a semi-Markov model. Zoning by K94 was carried out according to occurrences or epicenters, while the zoning in the present study was carried out by fault lines. In other words, K94 considered the historical data of earthquake occurrences according to their longitudes and latitudes, and classified them into 22 classes; while the present study considered the historical data of earthquake occurrences according to the longitudes and latitudes of their fault lines, and obtained their linear regressions, and then classify these linear regressions into 10 classes. Since the epicenters are positioned on the fault lines, zoning by fault lines can be more accurate than by epicenters. a) Zoning by K94 (based on epicenters) Karakaisis [1994] divided the Iran area into 21 zones. In this section, his zoning method is applied. K94 does not considered the center of Iran as a zone; hence the zoning is used to divide Iran into 22 zones, i.e. 21 zones by K94, plus one central part, as shown in Figure 2. Also, the magnitudes of the past occurrences were classified into 5 classes, as follows: (6) USING SEMI-MARKOV MODELS FOR ZONING 2 ( ) MAPE k MAPE j 2 1j k2 = = / FRMt FRMt FRMt ( ) ( ) 1, ...,FRM FRM i i ki1 = =6t t 1 ( ), 1, ...,FRM j j k=6t 2 ( ) 1 ( ) , 1, ...,FRM j M FRM j j k j = =6t t; E ( ), 1, ...,FRM j j k2 =6t ( ), 1, ...,FRM j j kD =6t 10 45 5=` j Table 1. A reference for distinguishing foreshocks and aftershocks from main-shocks [from Sadeghian and Jalali-Naini 2008b]. : 3.6M mb1 # : 3.6 .M mb 4 82 1 1 : 4.8 .M mb 5 43 1 1 : 5.4 .M mb 6 34 1 1 : 6.3M mb5 1 These classes were obtained through the agglomerative nesting (AGNES) method, which is a technique for data clustering [Ghazanfari and Rezaei 2006, Jalali-Naini and Sadeghian 2009]. The minimum of the considered magnitudes is 3.1 MB (body wave magnitude) and their maximum is 7.1 MB. With respect to the K94 zoning in Figure 2 and Equations (1), (2) and (6), the transition probability matrix for both the magnitude-to-magnitude and the region-to- region transitions are obtained as in Tables 7, 8a and 8b, which are included in Appendix 1. By applying the data in Tables 7, 8a and 8b and Equation (1), interval transition probability matrices in both magnitude-to-magnitude and region-to-region transitions are determined. By using these and Equation (2), the probabilistic forecasting matrix for the next 5 time periods (i.e. the next 50 days) after normalizing are obtained (see also Tables 9a and 9b in Appendix 1, and more comments in the modeling section mentioned above). Considering n = 3179, n1 = 3000, n2 = 179, hence, in zoning by K94, the values of the forecasting error of the algorithms mentioned are calculated as follows: In method 3, according to algorithm (II) in the validation section, all of the data are divided into three classes. The data ranging from 1 to 3,000 are used for forecasting the next 104 data, which are used for making a pattern. The data ranging from 3,001 to 3,104, which is equal to 28 time periods (each time unit is equal to 10 days), are used for deterministic forecasting of the data ranging from 3,105 to 3,179, which are the data later used for determining the forecasting error in the case of deterministic forecasting. In this way, according to algorithm (II) in the validation section, t is equal to 5. This value means that the element(s) greater than 5th maximum order in forecasting matrixes are replaced by 1 and the other elements are replaced by 0, then the deterministic forecasting is the nearest forecasting to the real occurrences and its forecasting error is the least. However by considering t= 5 in this zoning method, its MAPE is 2.398%. b) The proposed zoning (based on fault lines) In this zoning method, the estimated linear equations of fault lines (according to their longitudes and latitudes) were obtained. Next, the fault lines and their linear equations were classified into 10 zones according to specific attributes, including the Provinces of Iran, their importance, the length of the fault lines, etc. Therefore, Figure 3 was drawn approximately. The above steps were completely run on the MATLAB software. In this way, the fault lines in the Iran area (shown in Figure 4) [Sadeghian 2007] are divided into 10 zones, as shown in Figure 3. The classes of the magnitudes are considered similar to the K94 zoning method. The matrices GR (i.e. probability matrixof region-to-region transitions) and are determined through a method similar to the K94 zoning method, as with Table 10, in Appendix 2. By using Tables 7 and 10, the interval transition probability matrices in both region-to-region and magnitude- USING SEMI-MARKOV MODELS FOR ZONING 88 Figure 2. The siesmogenic source zones of Iran, as proposed by K94 [from Karakaisis 1994; re-edited by Sadeghian and Jalali-Naini 2008b]. Figure 3. The siesmogenic source zones of Iran proposed by Sadeghian (present study). 0.022MSE = 0.0MAD 5= .0 %MAPE 5 0= ( ) 1, ..., 5FRM i i =6 89 to-magnitude transitions are determined; in addition, using these and Equation (2), the probabilistic forecasting matrices for the next 5 time periods (i.e. the next 50 days) after normalizing are determined. (See Tables 11a and 11b in Appendix 2). In this zoning method, the forecasting errors are determined through a method similar to the K94 zoning method, as follows: In method 3, similar to the zoning by K94, all of data were divided into three sections. The value of tobtained was 4 and its MAPE in deterministic forecasting error is 3.552%. Discussion With respect to the two zoning methods investigated, the percentages of earthquake occurrences in each of the zones are shown in Figures 5 and 6. Moreover, the percentages of earthquake occurrences in each class of magnitudes, which are similar in both methods, are presented in Figure 7. The transition probability matrices in the magnitude-to- magnitude transitions are the same with both of the zoning methods. Table 7 shows that the maximum and minimum probabilities in the magnitude-to-magnitude transitions are related to the transitions from M5 to M2 (1.00) and from M2 to USING SEMI-MARKOV MODELS FOR ZONING Figure 4. The fault lines in the Iran area [from the website of the International Institute of Earthquake Engineering and Seismology (Tehran, Iran), http://www.iiees.ac.ir/iiees/English/Seismology/eng_seis_faults.html]. 0.0 2MSE 11= 0.0 5MAD 44= . %MAPE 4 45= M5 (0.0026), respectively (except for none-zero elements). These results are illustrated in Figure 7. With respect to the zoning by K94, presented in Tables 8a and 8b, it is clear that the maximum and the minimum probabilities in these tables are related to the R8 →R18 (0.333) and R22 → R6, R22 → R8 (0.0028) transitions, respectively. Moreover, in the proposed zoning (Table 10), the maximum and minimum probabilities are related to R2 →R7 (0.355) and R2 →R10 (0.0059) transitions, respectively. A summary of the minimum and maximum transition probabilities is shown in Table 2. The results forecast for the probabilities of earthquake occurrence in the future 5 time periods with the zoning by K94 and the proposed zoning are shown in Tables 9a and 9b, and Tables 11a and 11b, respectively. These tables show that in zoning by K94, the maximum probability of earthquake occurrences pertains to R22M2 in the first time period, and R18M2 in the second, third, fourth and fifth time periods. Similarly, in the proposed zoning, the maximum probability of earthquake occurrences pertains to R5M2 in the first time period, and R7M2 in the second, third, fourth and fifth time periods. Ultimately, the forecasting errors in both of the zoning methods show that both the number of zones and the type of zoning affect the forecasting probabilities and the forecasting errors. MAPE in the proposed zoning is 4.45%, while in zoning by K94, it is 5.0%. Therefore, it can be concluded that the proposed zoning will be more successful and precise, although both of these methods have rather low forecasting errors. A summary of the maximum probabilities in the two types of zoning is shown in Table 3. In the deterministic forecast method, for different t's, the minimum errors of the forecasting are calculated. So, by considering the forecasting errors obtained through method 3, it can be easily determined that the minimum error of forecasting is related to t=5 in the zoning by K94, and to t= 4 in the proposed zoning. The MAPEs in the deterministic forecasting error of zoning by K94 and by the proposed zoning can be obtained through the aforementioned method (i.e. algorithm III). These values are 2.398% in the zoning by K94, and 3.552% in the proposed zoning. These values show that in comparison to the proposed zoning, the K94 zoning method has fewer errors. Of course, this priority has a rational reason, and that is because in method 3, although the sum of the absolute percentage error (SAPE) is the same for both methods, to find the MAPE in the K94 zoning method, the SAPE has to be divided by 22×5, while, in the proposed zoning, it has to be divided by 10×5 (since the number of the classes of magnitudes is five, the number of the K94 zones is 22, while in the proposed zoning this is just 10); therefore, the MAPE obtained through the K94 zoning method is less than that obtained through the proposed method in the present study. Of course, this priority exists in USING SEMI-MARKOV MODELS FOR ZONING 90 Figure 5. The percentage of occurrences in each region, with the zoning by K94 [from Sadeghian and Jalali-Naini 2008b]. Figure 6. The percentage of occurrences in each region, with the zoning by Sadeghian (present study). Figure 7. The percentage of occurrences in each class of magnitudes [from Sadeghian and Jalali-Naini 2008b]. 91 USING SEMI-MARKOV MODELS FOR ZONING Table 6. Summary of the results of the forecasting. Table 2. Summary of minimum and maximum transition probabilities. Table 3. Summary of maximum probabilities in the two types of zoning. Table 4. Deterministic forecasting matrix during the first to fifth time periods for zoning by K94 and by the proposed zoning. Table 5. Actual earthquake occurrences during the following 5 time periods after the last earthquake occurrence. Transition Minimum probability Maximum probability Magnitude-to-magnitude transition probability M2 ! M5 (0.0026) M5 ! M2 (1.00) Region-to-region transition probability for the K94 zoning R22 ! R 6 , R 22 ! R 8 (0.0028) R 8 ! R 18 (0.333) Region-to-region transition probability for the proposed zoning R2 ! R 10 (0.0059) R2 ! R 7 (0.355) Zoning type Period1 Period2 Period3 Period4 Period5 Maximum probabilities in K94 zoning R22 M2 R18 M2 R18 M2 R18 M2 R18 M2 Maximum probabilities in the proposed zoning R5 M2 R7 M2 R7 M2 R7 M2 R7 M2 Zoning type Future time periods 1 to 5 K94 zoning method R16 M2, R18 M2, R19 M2, R20 M2, R22 M2 Proposed zoning method R5 M2, R6 M2, R7 M2, R8 M2 Period K94 zoning method Proposed zoning method 1 R16 M2, R20 M2, R21 M2 R6 M2, R8 M2, 2 R2 M2, R13 M2, R14 M2 R4 M2, R5 M2, R9 M2, 3 R7 M2, R17 M2, R18 M2, R19 M2, R22 M2 R2 M2, R4 M2, R7 M2, R8 M2, 4 R12 M2, R18 M2, R18 M2, R19 M2, R20 M3 R5 M2, R6 M3, R7 M1, R7 M2, 5 R20 M2, R21 M2, R22 M3 R6 M2, R6 M3, R7 M2, Zoning type Forecast (%) Completely correct Region correct, Magnitude wrong Adjacent region correct, Magnitude correct Not forecast K94 zoning 42 16 16 26 Proposed zoning 56 19 19 6 deterministic forecasting just due to the greater number of zones in the K94 zoning method, so it does not exist for probabilistic forecasting. Conclusion In this study, first, the method of applying the semi- Markov model in forecasting these triad dimensions of earthquakes was mentioned. Then, the method of calculating the probabilistic and deterministic forecasting was elaborated on. Moreover, several algorithms for validating the proposed model and obtaining the forecasting errors are introduced. Finally, while introducing two zoning methods, as those of K94 and of the proposed zoning, a case study on the Iran area is carried out and the forecasting errors of the proposed models are obtained for both of these zoning methods. In conclusion, it can noted that both the number of zones and the type of zoning affect the forecast probabilities and the forecasting errors. The forecasting errors show that the MAPE for the K94 zoning is 5.0%, while it is just 4.45% in the proposed zoning. These levels show that the forecasting error can be reduced through manipulating the numbers of zones or the types of zones. With respect to method 3, in calculating the deterministic forecasting errors, this method can also forecast the next earthquake occurrences deterministically. With respect to the collected data in the case study, the last earthquake considered in this study occurred on March 26, 2007, with a magnitude of 4.9 MB, M3, and in the regions of R22 and R5 (in the K94 and the proposed zonings, respectively). Accordingly, the deterministic forecasting during the future 5 time periods (the future 50 days) are determined in the K94 and the proposed zoning by algorithm III, as shown in Table 4. During the 50 days after carrying out the proposed model in the Iran area, some earthquakes occurred for both of the zoning methods. These are presented in Table 5. These results obtained show that in the zoning by K94, 42% of the earthquakes were correctly forecast. The areas of 16% of the earthquakes were correctly forecast, but not their magnitudes. Then 16% of the earthquakes occurred in areas adjacent to the forecast earthquakes, with the correct magnitudes, and 26% of the earthquakes were not forecast at all. In the proposed zoning, 56% of the earthquakes were correctly forecast in all of the three dimensions. The areas of 19% of the earthquakes were correctly forecast, but not their magnitudes. Then 19% of the earthquakes occurred in areas adjacent to the forecast earthquakes, with the correct magnitudes, and 6% of the earthquakes were not forecast at all. The summary of these results is given in Table 6. In addition to the errors obtained and explained above, the statistics mentioned here show that in comparison to the K94 zoning, the proposed zoning provides more accurate forecasting. One of the reasons for the priority of the proposed zoning over that of K94 is that the K94 zoning was not performed correctly (i.e. it did not take into account the center of Iran, which is a very extended area, while this was part of the zones considered in the present study). References Altinok, Y. and D. Kolcak (1999). An application of the semi- Markov model for earthquake occurrences in north Ana- tolia, Turkey, J. Balkan Geophys. Soc., 2 (4), 90-99. Console, R. (2001). Testing earthquake forecast hypotheses, Tectonophysics, 338, 261-268. Console, R., D. Pantosti and G. D'Addezio (2002). Probabilis- tic approach to earthquake prediction, Annals of Geo- physics, 45 (6), 723-731. Di Luccio, F., R. Console, M. Imoto and M. Murru (1997). Analysis of short time-space range seismicity pattern in Italy, Annals of Geophysics, 40 (4), 783-798. Gardner, J.K. and L. Knopoff (1974). Is the sequence of earth- quakes in southern California, with aftershocks removed, Poissonian?, B. Seismol. Soc. Am., 64, 1363-1367. Ghazanfari, M. and M. Rezaei (2006). Introduction on Fuzzy Theory and its Applications, IUST Press, Tehran. Haj Shirmohammadi, A. (2003). Principle of production and inventory planning and control, Arkan Press, Isfahan. Jalali-Naini, G.R. (1997). Stochastic processes and its applica- tions in biology, agriculture, management and industries, Statistics and Information Office Press, Tehran. Jalali-Naini, G.R. and R. Sadeghian (2009). Forecasting time and place of earthquakes using a semi-Markov model, Mathemathical Modelling and Applied Computing, 2, 1. Karakaisis, G.F. (1994). Long-term earthquake prediction in Iran based on the time- and magnitude-predictable model, Phys. Earth Planet. In., 83, 129-145. Kelleher, J., L. Sykes and O. Jack (1973). Possible Criteria for Predicting Earthquake Locations and Their Application to Major Plate Boundaries of the Pacific and the Carib- bean, J.Geophys. Res., 78 (14), 2547-2585. Minh, D.L. (2001). Applied Probability Models, Duxbury, Thomson Learning Press, California. Nanjo, Z.K., R.J. Holliday, C. Chen, B.J. Rundle and L.D. Turcotte (2006). Application of a modified pattern infor- matics method to forecasting the locations of future Large earthquakes in central Japan, Tectonophysics, 424, 351-366. Nava, F.A., C. Herrera, J. Frez and E. Glowacka (2005). Seismic hazard evaluation using Markov chains: application to the Japan area, Pure Appl. Geophys., 162, 1347-1366. Papazachos, B.C. (1992). A time- and magnitude-predictable model for generation of shallow earthquakes in the Aegean area, Pure Appl. Geophys., 138 (2), 287-308. Patwardhan, S.A., B.R. Kulkarni and D. Tocher (1980). A semi- Markov model for characterizing recurrence of great earth- quakes, B. Seismol. Soc. Am., 70 (1), 323-347. USING SEMI-MARKOV MODELS FOR ZONING 92 93 Rice, J.A. (2007). Mathematical Statistics and Data Analysis, Thomson Press, California. Rikitake, T. (1976). Recurrence of great earthquakes at sub- duction zones, Tectonophysics, 35, 335-362. Sadeghian, R. and G.R. Jalali-Naini (2008a). A new probability density function in earthquake occurrences, J. Ind. Engi- neering Int., 4 (6), 1-6. Sadeghian, R. and G.R. Jalali-Naini (2008b). Applying semi- Markov models for forecasting the triple dimensions of next earthquake occurrences, Dynamics of Continuous, Discrete and Impulsive Systems, Series B: Applications & Algorithms, 15, 395-416. Sadeghian, R. (2007). Whatever You Don't Know about Earth- quakes, Shooka Press, Tehran. *Corresponding author: Ramin Sadeghian, Bu Ali Sina University, Hamedan, Iran; e-mail: rsadeghian@basu.ac.ir. © 2010 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. USING SEMI-MARKOV MODELS FOR ZONING USING SEMI-MARKOV MODELS FOR ZONING 94 0. 00 0 0. 02 8 0. 00 0 0. 08 3 0. 00 0 0. 00 0 0. 00 0 0. 02 8 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 04 5 0. 09 1 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 06 9 0. 00 0 0. 03 4 0. 03 4 0. 00 0 0. 00 0 0. 03 4 0. 00 0 0. 00 0 0. 00 0 0. 13 8 0. 01 2 0. 00 6 0. 01 8 0. 07 9 0. 06 7 0. 01 2 0. 02 4 0. 00 6 0. 01 8 0. 00 0 0. 01 2 0. 00 0 0. 00 0 0. 00 0 0. 12 6 0. 26 4 0. 00 0 0. 01 1 0. 00 0 0. 00 0 0. 01 1 0. 01 1 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 04 3 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 04 3 0. 00 9 0. 00 9 0. 00 0 0. 01 8 0. 05 3 0. 00 0 0. 03 5 0. 00 0 0. 03 5 0. 01 8 0. 02 6 0. 00 0 0. 00 0 0. 05 6 0. 05 6 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 05 6 0. 00 0 0. 00 0 0. 01 5 0. 00 0 0. 00 0 0. 09 0 0. 00 0 0. 01 5 0. 03 0 0. 01 5 0. 04 5 0. 00 0 0. 03 0 0. 00 0 0. 00 0 0. 00 0 0. 06 5 0. 03 2 0. 00 0 0. 06 5 0. 00 0 0. 00 0 0. 03 2 0. 00 0 0. 00 0 0. 00 0 0. 01 0 0. 04 0 0. 01 0 0. 01 0 0. 04 0 0. 00 0 0. 03 0 0. 01 0 0. 28 0 0. 00 0 0. 01 8 0. 01 8 0. 01 8 0. 01 8 0. 00 0 0. 05 4 0. 01 8 0. 00 0 0. 03 6 0. 07 1 0. 05 1 0. 01 0 0. 01 0 0. 05 1 0. 03 1 0. 01 0 0. 05 1 0. 00 0 0. 01 0 0. 01 0 0. 02 0 0. 00 0 0. 00 0 0. 00 0 0. 07 3 0. 00 0 0. 00 0 0. 02 4 0. 00 0 0. 00 0 0. 00 0 0. 02 4 0. 00 0 0. 00 0 0. 00 0 0. 04 8 0. 04 8 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 04 8 0. 00 0 0. 00 7 0. 02 1 0. 01 0 0. 03 1 0. 03 1 0. 01 0 0. 02 8 0. 01 0 0. 03 1 0. 00 7 0. 02 8 0. 01 6 0. 00 0 0. 01 2 0. 05 3 0. 01 6 0. 00 0 0. 04 1 0. 01 6 0. 01 6 0. 01 6 0. 01 2 0. 01 5 0. 00 7 0. 01 3 0. 04 6 0. 01 8 0. 01 1 0. 03 7 0. 00 4 0. 02 0 0. 00 0 0. 02 4 0. 01 3 0. 00 3 0. 00 6 0. 05 4 0. 00 6 0. 00 6 0. 03 5 0. 00 3 0. 02 8 0. 01 3 0. 01 6 0. 00 9 0. 00 3 0. 00 3 0. 03 5 0. 00 9 0. 01 5 0. 04 7 0. 00 9 0. 02 6 0. 01 7 0. 03 2 0. 00 4 0. 00 7 0. 00 4 0. 06 7 0. 02 6 0. 00 7 0. 04 9 0. 00 0 0. 01 1 0. 00 7 0. 02 2 0. 01 1 0. 01 1 0. 00 8 0. 05 4 0. 01 7 0. 00 3 0. 03 4 0. 00 3 0. 02 5 0. 01 1 0. 02 3 0. 00 0 0. 11 1 0. 00 0 0. 00 0 0. 08 3 0. 11 1 0. 22 2 0. 08 3 0. 08 3 0. 05 6 0. 11 1 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 04 5 0. 09 1 0. 09 1 0. 09 1 0. 22 7 0. 13 6 0. 18 2 0. 00 0 0. 10 3 0. 03 4 0. 00 0 0. 03 4 0. 06 9 0. 10 3 0. 06 9 0. 10 3 0. 03 4 0. 13 8 0. 01 8 0. 02 4 0. 01 2 0. 01 2 0. 09 1 0. 06 1 0. 14 0 0. 11 0 0. 07 3 0. 06 1 0. 14 0 0. 00 0 0. 02 3 0. 02 3 0. 00 0 0. 08 0 0. 08 0 0. 06 9 0. 09 2 0. 05 7 0. 08 0 0. 06 9 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 04 3 0. 04 3 0. 17 4 0. 00 0 0. 21 7 0. 21 7 0. 21 7 0. 03 5 0. 06 1 0. 00 0 0. 00 9 0. 04 4 0. 07 0 0. 16 7 0. 08 8 0. 11 4 0. 08 8 0. 12 3 0. 00 0 0. 00 0 0. 00 0 0. 00 0 0. 16 7 0. 05 6 0. 33 3 0. 16 7 0. 11 1 0. 00 0 0. 00 0 0. 01 5 0. 06 0 0. 00 0 0. 00 0 0. 07 5 0. 04 5 0. 13 4 0. 13 4 0. 11 9 0. 10 4 0. 07 5 0. 03 2 0. 00 0 0. 03 2 0. 00 0 0. 06 5 0. 09 7 0. 06 5 0. 16 1 0. 19 4 0. 00 0 0. 16 1 0. 02 0 0. 01 0 0. 01 0 0. 01 0 0. 11 0 0. 03 0 0. 08 0 0. 05 0 0. 08 0 0. 05 0 0. 12 0 0. 17 9 0. 03 6 0. 00 0 0. 00 0 0. 01 8 0. 12 5 0. 12 5 0. 00 0 0. 08 9 0. 01 8 0. 16 1 0. 03 1 0. 16 3 0. 03 1 0. 01 0 0. 10 2 0. 06 1 0. 08 2 0. 02 0 0. 09 2 0. 07 1 0. 08 2 0. 00 0 0. 00 0 0. 09 8 0. 00 0 0. 12 2 0. 04 9 0. 17 1 0. 04 9 0. 07 3 0. 14 6 0. 17 1 0. 00 0 0. 09 5 0. 00 0 0. 00 0 0. 19 0 0. 04 8 0. 04 8 0. 19 0 0. 09 5 0. 09 5 0. 09 5 0. 01 0 0. 01 7 0. 01 4 0. 00 7 0. 14 0 0. 10 1 0. 14 3 0. 11 9 0. 08 4 0. 05 6 0. 09 1 0. 02 4 0. 02 0 0. 01 2 0. 00 8 0. 09 3 0. 09 3 0. 16 7 0. 07 7 0. 11 4 0. 08 5 0. 10 6 0. 01 5 0. 02 0 0. 01 8 0. 00 7 0. 09 4 0. 07 4 0. 18 2 0. 09 8 0. 09 8 0. 09 2 0. 10 7 0. 00 9 0. 04 4 0. 01 3 0. 00 6 0. 07 0 0. 08 9 0. 13 0 0. 12 0 0. 07 3 0. 10 8 0. 15 5 0. 01 5 0. 01 7 0. 00 9 0. 00 3 0. 07 6 0. 08 4 0. 16 0 0. 10 2 0. 17 2 0. 08 7 0. 07 3 0. 01 5 0. 01 5 0. 01 1 0. 00 4 0. 08 6 0. 08 6 0. 10 9 0. 10 1 0. 12 0 0. 11 6 0. 13 1 0. 01 1 0. 02 8 0. 00 6 0. 01 4 0. 09 9 0. 05 6 0. 15 2 0. 12 7 0. 12 4 0. 07 6 0. 10 7 0. 04 82 0. 91 57 0. 02 41 0. 01 20 0. 00 00 0. 02 75 0. 83 68 0. 11 08 0. 02 22 0. 00 26 0. 01 12 0. 81 23 0. 13 17 0. 04 48 0. 00 00 0. 02 60 0. 75 32 0. 19 48 0. 01 30 0. 01 30 0. 00 00 1. 00 00 0. 00 00 0. 00 00 0. 00 00 A p p en d ix 1 :T ra ns it io n pr ob ab ili ty m at ri ce s in th e zo ni ng a cc or di ng to K 94 T ab le 7 (t op ). T ra ns it io n pr ob ab ili ty m at ri x of th e m ag ni tu de -t o- m ag ni tu de tr an si ti on s in th e zo ni ng b y K 94 . T ab le 8 a/ b (b ot to m , l ef t/ ri gh t) . P ro ba bi lit y m at ri x of th e re gi on -t o- re gi on tr an si ti on s in th e zo ni ng b y K 94 . 95 USING SEMI-MARKOV MODELS FOR ZONING 0. 00 07 0. 05 36 0. 01 90 0. 00 34 0. 00 00 0. 00 18 0. 06 68 0. 01 36 0. 00 18 0. 00 00 0. 00 07 0. 05 36 0. 01 90 0. 00 34 0. 00 00 0. 00 18 0. 05 33 0. 01 08 0. 00 18 0. 00 00 0. 00 07 0. 04 00 0. 01 42 0. 00 27 0. 00 00 0. 00 18 0. 05 15 0. 01 08 0. 00 18 0. 00 00 0. 00 27 0. 22 66 0. 08 07 0. 01 42 0. 00 00 0. 01 08 0. 36 31 0. 07 32 0. 01 17 0. 00 09 0. 00 07 0. 08 01 0. 02 85 0. 00 54 0. 00 00 0. 00 45 0. 15 36 0. 03 07 0. 00 45 0. 00 00 0. 00 00 0. 01 36 0. 00 47 0. 00 07 0. 00 00 0. 00 09 0. 04 34 0. 00 90 0. 00 18 0. 00 00 0. 00 14 0. 12 01 0. 04 27 0. 00 75 0. 00 00 0. 00 81 0. 27 73 0. 05 60 0. 00 90 0. 00 09 0. 00 00 0. 01 36 0. 00 47 0. 00 07 0. 00 00 0. 00 09 0. 03 61 0. 00 72 0. 00 09 0. 00 00 0. 00 14 0. 10 65 0. 03 80 0. 00 68 0. 00 00 0. 00 45 0. 16 17 0. 03 25 0. 00 54 0. 00 00 0. 00 07 0. 05 36 0. 01 90 0. 00 34 0. 00 00 0. 00 18 0. 06 50 0. 01 26 0. 00 18 0. 00 00 0. 00 14 0. 10 65 0. 03 80 0. 00 68 0. 00 00 0. 00 54 0. 17 98 0. 03 61 0. 00 54 0. 00 00 0. 00 07 0. 04 00 0. 01 42 0. 00 27 0. 00 00 0. 00 36 0. 11 74 0. 02 35 0. 00 36 0. 00 00 0. 00 14 0. 12 01 0. 04 27 0. 00 75 0. 00 00 0. 00 72 0. 22 40 0. 04 52 0. 00 72 0. 00 09 0. 00 00 0. 02 65 0. 00 95 0. 00 20 0. 00 00 0. 00 18 0. 06 96 0. 01 45 0. 00 18 0. 00 00 0. 00 07 0. 06 65 0. 02 37 0. 00 41 0. 00 00 0. 00 18 0. 05 24 0. 01 08 0. 00 18 0. 00 00 0. 00 47 0. 41 32 0. 14 79 0. 02 65 0. 00 00 0. 01 99 0. 66 31 0. 13 28 0. 02 08 0. 00 18 0. 00 27 0. 22 66 0. 08 07 0. 01 42 0. 00 00 0. 01 54 0. 50 68 0. 10 21 0. 01 63 0. 00 09 0. 00 75 0. 63 98 0. 22 86 0. 04 14 0. 00 00 0. 03 07 1. 00 00 0. 20 14 0. 03 16 0. 00 27 0. 00 61 0. 51 97 0. 18 59 0. 03 32 0. 00 00 0. 02 26 0. 75 25 0. 15 18 0. 02 35 0. 00 18 0 .0 06 1 0. 47 96 0. 17 16 0. 03 05 0. 00 00 0. 02 53 0. 83 47 0. 16 80 0. 02 62 0. 00 18 0. 00 41 0. 33 31 0. 11 94 0. 02 17 0. 00 00 0. 01 81 0. 58 27 0. 11 74 0. 01 81 0. 00 18 0 .0 12 2 1. 00 00 0. 35 75 0. 06 45 0. 00 00 0. 02 71 0. 88 62 0. 17 80 0. 02 80 0. 00 18 0. 00 25 0. 07 12 0. 01 09 0. 00 17 0. 00 00 0. 00 25 0. 06 96 0. 00 98 0. 00 16 0. 00 00 0. 00 17 0. 04 69 0. 00 67 0. 00 17 0. 00 00 0. 00 16 0. 04 75 0. 00 66 0. 00 16 0. 00 00 0. 00 17 0. 05 28 0. 00 84 0. 00 17 0. 00 00 0. 00 16 0. 05 81 0. 00 82 0. 00 16 0. 00 00 0. 01 09 0. 35 54 0. 05 36 0. 01 01 0. 00 08 0. 01 06 0. 35 71 0. 05 24 0. 01 06 0. 00 08 0. 00 50 0. 16 60 0. 02 51 0. 00 50 0. 00 08 0. 00 49 0. 16 79 0. 02 46 0. 00 49 0. 00 08 0. 00 17 0. 04 53 0. 00 67 0. 00 17 0. 00 00 0. 00 16 0. 04 42 0. 00 66 0. 00 16 0. 00 00 0. 00 75 0. 24 90 0. 03 77 0. 00 75 0. 00 08 0. 00 74 0. 24 57 0. 03 60 0. 00 74 0. 00 08 0. 00 08 0. 04 02 0. 00 59 0. 00 08 0. 00 00 0. 00 08 0. 04 01 0. 00 57 0. 00 08 0. 00 00 0. 00 42 0. 14 75 0. 02 26 0. 00 42 0. 00 00 0. 00 41 0. 14 33 0. 02 05 0. 00 41 0. 00 00 0. 00 17 0. 06 45 0. 01 01 0. 00 17 0. 00 00 0. 00 16 0. 06 31 0. 00 90 0. 00 16 0. 00 00 0. 00 59 0. 19 53 0. 02 93 0. 00 59 0. 00 08 0. 00 57 0. 19 82 0. 02 87 0. 00 57 0. 00 08 0. 00 34 0. 11 57 0. 01 76 0. 00 34 0. 00 00 0. 00 33 0. 11 55 0. 01 72 0. 00 33 0. 00 00 0. 00 67 0. 22 05 0. 03 35 0. 00 67 0. 00 08 0. 00 66 0. 21 62 0. 03 11 0. 00 66 0. 00 08 0. 00 25 0. 07 88 0. 01 17 0. 00 25 0. 00 00 0. 00 25 0. 08 03 0. 01 15 0. 00 25 0. 00 00 0. 00 17 0. 04 36 0. 00 67 0. 00 17 0. 00 00 0. 00 16 0. 04 59 0. 00 66 0. 00 16 0. 00 00 0. 01 93 0. 63 20 0. 09 64 0. 01 84 0. 00 17 0. 01 88 0. 62 00 0. 09 01 0. 01 80 0. 00 16 0. 01 68 0. 56 08 0. 08 55 0. 01 68 0. 00 17 0. 01 64 0. 55 45 0. 08 11 0. 01 56 0. 00 16 0. 03 02 1. 00 00 0. 15 17 0. 02 93 0. 00 25 0. 03 03 1. 00 00 0. 14 58 0. 02 87 0. 00 25 0. 02 01 0. 67 31 0. 10 23 0. 01 93 0. 00 17 0. 01 97 0. 64 37 0. 09 42 0. 01 88 0. 00 16 0 .0 22 6 0. 76 53 0. 11 65 0. 02 26 0. 00 17 0. 02 21 0. 73 55 0. 10 73 0. 02 13 0. 00 16 0. 01 76 0. 58 09 0. 08 80 0. 01 68 0. 00 17 0. 01 72 0. 57 58 0. 08 35 0. 01 64 0. 00 16 0 .0 24 3 0. 80 13 0. 12 15 0. 02 35 0. 00 17 0. 02 29 0. 76 33 0. 11 14 0. 02 21 0. 00 16 A p p en d ix 1 :T ra ns it io n pr ob ab ili ty m at ri ce s in th e zo ni ng a cc or di ng to K 94 T ab le 9 a (l ef t to r ig ht ). P ro ba bi lis ti c fo re ca st in g m at ri xe s du ri ng t he fi rs t, s ec on d, t hi rd a nd fo u rt h ti m e pe ri od s in t he z on in g by K 94 . USING SEMI-MARKOV MODELS FOR ZONING 96 0. 21 73 0. 03 19 0. 01 60 0. 03 19 0. 06 39 0. 15 34 0. 23 96 0. 19 81 0. 02 88 0. 01 92 0. 06 51 0. 10 06 0. 01 78 0. 00 00 0. 09 47 0. 20 12 0. 35 50 0. 11 83 0. 04 14 0. 00 59 0. 06 52 0. 07 61 0. 28 26 0. 03 26 0. 05 43 0. 16 30 0. 16 30 0. 09 78 0. 03 26 0. 03 26 0. 05 13 0. 06 41 0. 02 56 0. 05 13 0. 05 13 0. 17 95 0. 32 05 0. 23 08 0. 01 28 0. 01 28 0. 08 91 0. 06 93 0. 03 96 0. 02 97 0. 16 34 0. 14 36 0. 18 81 0. 18 81 0. 05 45 0. 03 47 0. 11 11 0. 04 96 0. 02 39 0. 02 56 0. 03 59 0. 23 76 0. 26 32 0. 20 00 0. 02 22 0. 03 08 0. 07 57 0. 05 16 0. 02 09 0. 02 31 0. 05 27 0. 18 66 0. 33 04 0. 19 98 0. 03 40 0. 02 52 0. 10 00 0. 04 84 0. 01 56 0. 02 34 0. 06 72 0. 15 94 0. 30 00 0. 22 19 0. 03 59 0. 02 81 0. 03 77 0. 03 77 0. 03 77 0. 01 89 0. 08 49 0. 19 81 0. 25 47 0. 24 53 0. 05 66 0. 02 83 0. 04 88 0. 06 10 0. 01 22 0. 02 44 0. 04 88 0. 15 85 0. 29 27 0. 31 71 0. 01 22 0. 02 44 0. 00 25 0. 07 12 0. 00 99 0. 00 25 0. 00 00 0. 00 17 0. 04 80 0. 00 66 0. 00 17 0. 00 00 0. 00 17 0. 06 13 0. 00 83 0. 00 17 0. 00 00 0. 01 08 0. 36 34 0. 05 22 0. 01 08 0. 00 08 0. 00 50 0. 17 22 0. 02 48 0. 00 50 0. 00 08 0. 00 17 0. 04 64 0. 00 66 0. 00 17 0. 00 00 0. 00 75 0. 25 00 0. 03 56 0. 00 75 0. 00 08 0. 00 08 0. 04 06 0. 00 58 0. 00 08 0. 00 00 0. 00 41 0. 14 40 0. 02 07 0. 00 41 0. 00 00 0. 00 17 0. 06 37 0. 00 91 0. 00 17 0. 00 00 0. 00 58 0. 20 28 0. 02 90 0. 00 58 0. 00 08 0. 00 33 0. 11 75 0. 01 66 0. 00 33 0. 00 00 0. 00 66 0. 21 85 0. 03 15 0. 00 66 0. 00 08 0. 00 25 0. 08 36 0. 01 16 0. 00 25 0. 00 00 0. 00 17 0. 04 64 0. 00 66 0. 00 17 0. 00 00 0. 01 90 0. 62 83 0. 08 94 0. 01 82 0. 00 17 0. 01 66 0. 56 29 0. 08 03 0. 01 66 0. 00 17 0. 02 98 1. 00 00 0. 14 24 0. 02 98 0. 00 25 0. 01 99 0. 65 23 0. 09 27 0. 01 90 0. 00 17 0 .0 22 4 0. 74 34 0. 10 60 0. 02 15 0. 00 17 0. 01 74 0. 58 28 0. 08 28 0. 01 74 0. 00 17 0 .0 23 2 0. 77 40 0. 11 01 0. 02 32 0. 00 25 A p p en d ix 1 . T ab le 9 b. P ro ba bi lis ti c fo re ca st in g m at ri xe s du ri ng t he fi ft h ti m e pe ri od in t he z on in g by K 94 . A p p en d ix 2 :T ra ns it io n pr ob ab ili ty m at ri ce s in th e pr op os ed z on in g (p re se nt s tu dy ) T ab le 1 0 (t op ). P ro ba bi lit y m at ri x of t he r eg io n- to -r eg io n tr an si ti on s in t he p ro po se d zo ni ng . T ab le 1 1a (b ot to m ). P ro ba bi lis ti c fo re ca st in g m at ri ce s du ri ng th e fir st a nd s ec on d ti m e pe ri od s in th e pr op os ed z on in g. 0. 00 41 0. 31 98 0. 11 41 0. 02 03 0. 00 00 0. 01 03 0. 33 97 0. 06 81 0. 01 08 0. 00 10 0. 00 29 0. 26 01 0. 09 27 0. 01 68 0. 00 00 0. 00 69 0. 22 60 0. 04 56 0. 00 69 0. 00 05 0. 00 17 0. 15 99 0. 05 74 0. 01 04 0. 00 00 0. 00 39 0. 12 55 0. 02 50 0. 00 39 0. 00 05 0. 00 12 0. 11 99 0. 04 29 0. 00 75 0. 00 00 0. 00 25 0. 08 87 0. 01 76 0. 00 29 0. 00 00 0. 01 22 1. 00 00 0. 35 75 0. 06 43 0. 00 00 0. 01 27 0. 42 79 0. 08 58 0. 01 32 0. 00 10 0. 00 58 0. 50 00 0. 17 90 0. 03 19 0. 00 00 0. 02 01 0. 66 23 0. 13 28 0. 02 06 0. 00 15 0. 00 81 0. 65 99 0. 23 58 0. 04 23 0. 00 00 0. 03 04 1. 00 00 0. 20 10 0. 03 14 0. 00 25 0. 00 81 0. 68 02 0. 24 33 0. 04 35 0. 00 00 0. 02 30 0. 76 42 0. 15 34 0. 02 40 0. 00 20 0. 00 23 0. 19 99 0. 07 13 0. 01 27 0. 00 00 0. 00 44 0. 14 36 0. 02 89 0. 00 44 0. 00 05 0. 00 17 0. 14 02 0. 04 98 0. 00 93 0. 00 00 0. 00 29 0. 09 85 0. 01 96 0. 00 29 0. 00 00 97 USING SEMI-MARKOV MODELS FOR ZONING A p p en d ix 2 :T ra ns it io n pr ob ab ili ty m at ri ce s in th e pr op os ed z on in g (p re se nt s tu dy ) T ab le 1 1b ( le ft t o ri gh t) . P ro ba bi lis ti c fo re ca st in g m at ri ce s du ri ng t he t hi rd , f ou rt h an d fif th t im e pe ri od s in t he p ro po se d zo ni ng . 0. 01 06 0. 35 27 0. 05 36 0. 01 01 0. 00 09 0. 01 03 0. 35 02 0. 05 11 0. 00 99 0. 00 09 0. 00 57 0. 19 44 0. 02 95 0. 00 57 0. 00 04 0. 00 56 0. 18 97 0. 02 75 0. 00 56 0. 00 04 0. 00 31 0. 10 38 0. 01 58 0. 00 31 0. 00 04 0. 00 30 0. 09 96 0. 01 46 0. 00 30 0. 00 04 0. 00 26 0. 08 09 0. 01 23 0. 00 22 0. 00 00 0. 00 26 0. 08 28 0. 01 20 0. 00 26 0. 00 04 0. 00 88 0. 28 63 0. 04 35 0. 00 84 0. 00 09 0. 00 73 0. 24 12 0. 03 52 0. 00 69 0. 00 09 0. 02 02 0. 67 06 0. 10 20 0. 01 98 0. 00 18 0. 01 97 0. 65 79 0. 09 57 0. 01 89 0. 00 17 0. 02 99 1. 00 00 0. 15 22 0. 02 90 0. 00 26 0. 03 00 1. 00 00 0. 14 59 0. 02 88 0. 00 26 0. 02 15 0. 72 16 0. 10 99 0. 02 11 0. 00 18 0. 02 15 0. 72 36 0. 10 56 0. 02 06 0. 00 21 0. 00 35 0. 11 65 0. 01 76 0. 00 35 0. 00 04 0. 00 34 0. 11 55 0. 01 67 0. 00 34 0. 00 04 0. 00 31 0. 09 54 0. 01 45 0. 00 26 0. 00 04 0. 00 30 0. 09 57 0. 01 37 0. 00 26 0. 00 04 0. 01 02 0. 34 51 0. 04 92 0. 01 02 0. 00 08 0. 00 55 0. 18 65 0. 02 67 0. 00 55 0. 00 04 0. 00 30 0. 09 71 0. 01 40 0. 00 30 0. 00 04 0. 00 25 0. 08 35 0. 01 19 0. 00 25 0. 00 04 0. 00 68 0. 22 89 0. 03 26 0. 00 68 0. 00 04 0. 01 95 0. 64 99 0. 09 24 0. 01 91 0. 00 17 0. 03 01 1. 00 00 0. 14 24 0. 02 92 0. 00 25 0. 02 16 0. 72 02 0. 10 26 0. 02 12 0. 00 21 0. 00 34 0. 11 32 0. 01 61 0. 00 34 0. 00 04 0. 00 30 0. 09 50 0. 01 36 0. 00 30 0. 00 04