Assessment of seismic risk scenarios for Bucharest,Romania

In this paper, seismic risk scenarios for Bucharest, the capital city of Romania, are proposed and assessed. Bucharest has one of the highest seismic risk levels in Europe, and this is due to a combination of relatively high seismic hazard and a building stock built mainly before the devastating Vrancea 1977 earthquake. In this study, the seismic risk of Bucharest is assessed using the most recent information regarding the characteristics of the residential building stock. The ground motion amplitudes are evaluated starting from random fields obtained by coupling a ground motion model derived for the Vrancea intermediate-depth seismic source with a spatial correlation model. The seismic risk evaluation method applied in this study is based on the well-known macroseismic method. For several structural typologies, the vulnerability parameters are evaluated based on a damage survey performed on 18,000 buildings in Bucharest after the March 1977 earthquake. Subsequently, the risk metrics are compared with those from other studies in the literature that apply a different risk assessment methodology in order to gain a better view of the uncertainties associated with a seismic risk study at city level. Finally, the impact of several Vrancea intermediate-depth earthquake scenarios is evaluated and the results show that the earthquake which has the closest epicenter to Bucharest appears to be the most damaging.

     

Assessment of seismic risk scenarios for Bucharest,Romania

In this paper, seismic risk scenarios for Bucharest, the capital city of Romania, are proposed and assessed. Bucharest has one of the highest seismic risk levels in Europe, and this is due to a combination of relatively high seismic hazard and a building stock built mainly before the devastating Vrancea 1977 earthquake. In this study, the seismic risk of Bucharest is assessed using the most recent information regarding the characteristics of the residential building stock. The ground motion amplitudes are evaluated starting from random fields obtained by coupling a ground motion model derived for the Vrancea intermediate-depth seismic source with a spatial correlation model. The seismic risk evaluation method applied in this study is based on the well-known macroseismic method. For several structural typologies, the vulnerability parameters are evaluated based on a damage survey performed on 18,000 buildings in Bucharest after the March 1977 earthquake. Subsequently, the risk metrics are compared with those from other studies in the literature that apply a different risk assessment methodology in order to gain a better view of the uncertainties associated with a seismic risk study at city level. Finally, the impact of several Vrancea intermediate-depth earthquake scenarios is evaluated and the results show that the earthquake which has the closest epicenter to Bucharest appears to be the most damaging.     

Tectonic stress,seismicity, and seismic hazard in the southeastern Carpathians

Intermediate-depth earthquakes in the Vrancea region occur in response to stress generation due to descending lithosphere beneath the southeastern Carpathians. In this article, tectonic stress and seismicity are analyzed in the region on the basis of a vast body of observations. We show a correlation between the location of intermediate-depth earthquakes and the predicted localization of maximum shear stress in the lithosphere. A probabilistic seismic hazard assessment (PSHA) for the region is presented in terms of various ground motion parameters on the utilization of Fourier amplitude spectra used in engineering practice and risk assessment (peak ground acceleration, response spectra amplitude, and seismic intensity). We review the PSHA carried out in the region, and present new PSHA results for the eastern and southern parts of Romania. Our seismic hazard assessment is based on the information about the features of earthquake ground motion excitation, seismic wave propagation (attenuation), and site effect in the region. Spectral models and characteristics of site-response on earthquake ground motions are obtained from the regional ground motion data including several hundred records of small and large earthquakes. Results of the probabilistic seismic hazard assessment are consistent with the features of observed earthquake effects in the southeastern Carpathians and show that geological factors play an important part in the distribution of the earthquake ground motion parameters.     

3-D seismic structure of the Kachchh,Gujarat, and its implications for the earthquake hazard mitigation

Several pieces of studies on the January 26, 2001, Bhuj earthquake (Mw 7.6) revealed that the mainshock was triggered on the hidden unmapped fault in the western part of Indian stable continental region that caused a huge loss in the entire Kachchh rift basin of Gujarat, India. Occurrences of infrequent earthquakes of Mw 7.6 due to existence of hidden and unmapped faults on the surface have become one of the key issues for geoscientific research, which need to be addressed for evolving plausible earthquake hazard mitigation model. In this study, we have carried out a detailed autopsy of the 2001 Bhuj earthquake source zone by applying three-dimensional (3-D) local earthquake tomography (LET) method to a completely new data set consisting of 576 local earthquakes recorded between November 2006 and April 2009 by a seismic network consisting of 22 numbers of three-component broadband digital seismograph stations. In the present study, a total of 7560 arrival times of P-wave (3820) and S-wave (3740) recorded at least 4 seismograph stations were inverted to assimilate 3-D P-wave velocity (Vp), S-wave velocity (Vs), and Poisson’s ratio (σ) structures beneath the 2001 Bhuj earthquake source zone for reliable interpretation of the imaged anomalies and its bearing on earthquake hazard of the region. The source zone is located near the triple junction formed by juxtapositions of three Indian, Arabian, and Iranian tectonic plates that might have facilitated the process of brittle failure at a depth of 25 km beneath the KRB, Gujarat, which caused a gigantic loss to both property and persons of the region. There may be several hidden seismogenic faults around the epicentral zone of the 2001 Bhuj earthquake in the area, which are detectable using 3-D tomography to minimize earthquake hazard for a region. We infer that the use of detailed 3-D seismic tomography may offer potential information on hidden and unmapped faults beneath the plate interior to unravel the genesis of such big damaging earthquakes. This study may help in evolving a comprehensive earthquake risk mitigation model for regions of analogous geotectonic settings, elsewhere in the world.     

Influence of large dam on seismic hazard in low seismic region of Ulu Padas Area,Northern Borneo

The seismic hazard assessment of a site that lies in the low seismic region affected by the future existence of a large dam has been given less attention in many studies. Moreover, this condition is not addressed directly in the current seismic codes. This paper explains the importance of such information in mitigating the seismic hazard properly. Ulu Padas Area in Northern Borneo is used as an example for a case study of a site classified as a low seismic region. It is located close to the border of Malaysia, Brunei Darussalam, and Indonesia and may have a large dam in the future as the region lies in hilly geography with river flow. This study conducts probabilistic and deterministic seismic hazard analyses, and reservoir-triggered seismicity of a site affected by the future existence of a large dam. The result shows that the spectrum acceleration of the maximum design earthquake for the investigated site in the Ulu Padas Area in Northern Borneo is taken from the reservoir-triggered seismicity earthquake at short periods and from the current condition at longer periods.     

Ground zoning against seismic hazard in Athens, Greece

The 1999, Ms=5.9, Athens earthquake caused serious structural damage to buildings in the western part of Athens, Greece. This paper presents the ground zoning against seismic hazard proposed shortly after the earthquake in order to aid reconstruction of the area. Existing engineering geological and geotechnical data were combined with local observations to provide a unified set of classification criteria, consistent with provisions of the Greek Seismic Code EAK. The accuracy and the possible limitations of this zoning procedure are addressed through comparison with observed damage distribution as well as results from seismic ground response analyses performed at sites with well established soil profiles. There is clear evidence that the proposed zones correspond to geological formations exhibiting grossly different seismic response with regard to the design of common engineering structures. However, the mostly qualitative nature of the guidelines for ground categorisation provided by EAK and the general lack of systematic, site-specific geotechnical data for the whole area induce uncertainties in the definition of the seismic design actions for the different zones. These objective uncertainties certainly demand increased conservatism but do not limit application of the proposed methodology for first aid, preliminary planning in the event of destructive earthquakes.     

Seismicity, seismotectonics and seismic hazard of Italy

A first generation of probabilistic seismic hazard maps of the Italian country are presented. They are based on seismogenic zoning deriving from a kinematic model of the structural tectonic units and on an earthquake catalogue with the foreshock and aftershock events filtered out. The following ground motion parameters have been investigated and mapped using attenuation equations based on strong-motion recordings of Italian earthquakes: peak ground acceleration and velocity; Arias intensity; strong motion duration; and the pseudovelocity and pseudoacceleration spectral values at 14 fixed frequencies both for the vertical and the largest horizontal component. A Poissonian model of earthquake occurrence is assumed as a default and the hazard maps are presented in terms of ground motion values expected to be exceeded at a 10% probability level in 50 years (return period 475 years) according to the requirement of Eurocode 8 for the seismic classification of national territories, as well as in terms of exceedance probabilities of selected ground motion values. Finally, as a tentative study, the use of hybrid methods (implementing both seismogenic zones and structures), renewal processes (including earthquake forecasting) and the influence of site effects (as the basis for the planning of earthquake scenarios) were explored.     

Estimation of seismic hazard in Gujarat region, India

The seismic hazard in the Gujarat region has been evaluated. The scenario hazard maps showing the spatial distribution of various parameters like peak ground acceleration, characteristics site frequency and spectral acceleration for different periods have been presented. These parameters have been extracted from the simulated earthquake strong ground motions. The expected damage to buildings from future large earthquakes in Gujarat region has been estimated. It has been observed that the seismic hazard of Kachchh region is more in comparison with Saurashtra and mainland. All the cities of Kachchh can expect peak acceleration in excess of 500?cm/s² at surface in case of future large earthquakes from major faults in Kachchh region. The cities of Saurashtra can expect accelerations of less than 200?cm/s² at surface. The mainland Gujarat is having the lowest seismic hazard as compared with other two regions of Gujarat. The expected accelerations are less than 50?cm/s² at most of the places. The single- and double-story buildings in Kachchh region are at highest risk as they can expect large accelerations corresponding to natural periods of such small structures. Such structures are relatively safe in mainland region. The buildings of 3?C4 stories and tall structures that exist mostly in cities of Saurashtra and mainland can expect accelerations in excess of 100?cm/s² during a large earthquake in Kachchh region. It has been found that a total of 0.11 million buildings in Rajkot taluka of Saurashtra are vulnerable to total damage. In Kachchh region, 0.37 million buildings are vulnerable. Most vulnerable talukas are Bhuj, Anjar, Rapar, Bhachau, and Mandvi in Kachchh district and Rajkot, Junagadh, Jamnagar, Surendernagar and Porbandar in Saurashtra. In mainland region, buildings in Bharuch taluka are more vulnerable due to proximity to active Narmada-Son geo-fracture. The scenario hazard maps presented in this study for moderate as well as large earthquakes in the region may be used to augment the information available in the probabilistic seismic hazard maps of the region.     

Integrated seismic hazard evaluation and disaster management approach for Turkey

A significant proportion of the urban areas in Turkey is subject to high seismic risk. An important step for seismic risk mitigation is to define the hazard and damage after an earthquake. This paper proposes an integrated seismic hazard assessment and disaster management processes for Turkey. The proposed methodology utilizes information technologies in its seismic assessment component that provides fast results for assessment. First, image process methodology by using satellite images was implemented in the seismic assessment process for fast evaluation right after an earthquake. Second, the seismic assessment process was integrated with disaster management process. As a result, through integrated seismic hazard evaluation and disaster management procedure, an effective, fast and dependable estimation of loss for Turkey was planned.     

Global strain rates in western to central Himalayas and their implications in seismic hazard assessment

The Himalayas has experienced varying rates of earthquake occurrence in the past in its seismo-tectonically distinguished segments which may be attributed to different physical processes of accumulation of stress and its release, and due diligence is required for its inclusion for working out the seismic hazard. The present paper intends to revisit the various earthquake occurrence models applied to Himalayas and examines it in the light of recent damaging earthquakes in Himalayan belt. Due to discordant seismicity of Himalayas, three types of regions have been considered to estimate larger return period events. The regions selected are (1) the North-West Himalayan Fold and Thrust Belt which is seismically very active, (2) the Garhwal Himalaya which has never experienced large earthquake although sufficient stress exists and (3) the Nepal region which is very seismically active region due to unlocked rupture and frequently experienced large earthquake events. The seismicity parameters have been revisited using two earthquake recurrence models namely constant seismicity and constant moment release. For constant moment release model, the strain rates have been derived from global strain rate model and are converted into seismic moment of earthquake events considering the geometry of the finite source and the rates being consumed fully by the contemporary seismicity. Probability of earthquake occurrence with time has been estimated for each region using both models and compared assuming Poissonian distribution. The results show that seismicity for North-West region is observed to be relatively less when estimated using constant seismicity model which implies that either the occupied accumulated stress is not being unconfined in the form of earthquakes or the compiled earthquake catalogue is insufficient. Similar trend has been observed for seismic gap area but with lesser difference reported from both methods. However, for the Nepal region, the estimated seismicity by the two methods has been found to be relatively less when estimated using constant moment release model which implies that in the Nepal region, accumulated strain is releasing in the form of large earthquake occurrence event. The partial release in second event of May 2015 of similar size shows that the physical process is trying to release the energy with large earthquake event. If it would have been in other regions like that of seismic gap region, the fault may not have released the energy and may be inviting even bigger event in future. It is, therefore, necessary to look into the seismicity from strain rates also for its due interpretation in terms of predicting the seismic hazard in various segments of Himalayas.     

Local seismic hazard assessment in areas of weak to moderate seismicity—case study from Eastern Germany

Generally the seismic hazard of an area of interest is considered independent of time. However, its seismic risk or vulnerability, respectively, increases with the population and developing state of economy of the area. Therefore, many areas of moderate seismic hazard gain increasing importance with respect to seismic hazard and risk analysis. However, these areas mostly have a weak earthquake database, i.e., they are characterised by relative low seismicity and uncertain information concerning historical earthquakes. In a case study for Eastern Thuringia (Germany), acting as example for similar places in the world, seismic hazard is estimated using the probabilistic approach. Because of the lack of earthquakes occurring in the recent past, mainly historical earthquakes have to be used. But for these the actual earthquake sources or active faults, needed for the analysis, are imprecisely known. Therefore, the earthquake locations are represented by areal sources, a common practice. The definition of these sources is performed carefully, because their geometrical shape and size (apart from the earthquake occurrence model) influence the results significantly. Using analysis tools such as density maps of earthquake epicentres, seismic strain and energy release support this. Oversizing of areal sources leads to underestimation of seismic hazard and should therefore be avoided. Large location errors of historical earthquakes on the other hand are represented by several alternative areal sources with final superimposition of the different results. In a very similar way information known from macroseismic observations interpreted as source rather than as site effects are taken into account in order to achieve a seismic hazard assessment as realistic as possible. In very local cases the meaning of source effects exceeds those of site effects very likely. The influence of attenuation parameter variations on the result of estimated local seismic hazard is relatively low. Generally, the results obtained by the seismic hazard assessment coincide well with macroseismic observations from the thoroughly investigated largest earthquake in the region.     

Non-Poisson probabilistic seismic hazard assessment

The development of the new seismic hazard map of metropolitan Tehran is based on probabilistic seismic hazard computation using the non-Poisson recurrence time model. For this model, two maps have been prepared to indicate the earthquake hazard of the region in the form of iso-acceleration contour lines. They display the non-Poisson probabilistic estimates of peak ground accelerations over bedrock for 10 and 63 % probability of exceedance in 50 years. To carry out the non-Poisson seismic hazard analysis, appropriate distributions of interoccurrence times of earthquakes were used for the seismotectonic provinces which the study region is located and then the renewal process was applied. In order to calculate the seismic hazard for different return periods in the probabilistic procedure, the study area encompassed by the 49.5–54.5°E longitudes and 34–37°N latitudes was divided into 0.1° intervals generating 1,350 grid points. PGA values for this region are estimated to be 0.30–0.32 and 0.16–0.17 g for 10 and 63 % probability of exceedance, respectively, in 50 years for bedrock condition.     

Spatial-temporal variability of seismic hazard in peninsular India

This paper examines the variability of seismic activity observed in the case of different geological zones of peninsular India (10°N–26°N; 68°E–90°E) based on earthquake catalog between the period 1842 and 2002 and estimates earthquake hazard for the region. With compilation of earthquake catalog in terms of moment magnitude and establishing broad completeness criteria, we derive the seismicity parameters for each geologic zone of peninsular India using maximum likelihood procedure. The estimated parameters provide the basis for understanding the historical seismicity associated with different geological zones of peninsular India and also provide important inputs for future seismic hazard estimation studies in the region. Based on present investigation, it is clear that earthquake recurrence activity in various geologic zones of peninsular India is distinct and varies considerably between its cratonic and rifting zones. The study identifies the likely hazards due to the possibility of moderate to large earthquakes in peninsular India and also presents the influence of spatial rate variation in the seismic activity of this region. This paper presents the influence of source zone characterization and recurrence rate variation pattern on the maximum earthquake magnitude estimation. The results presented in the paper provide a useful basis for probabilistic seismic hazard studies and microzonation studies in peninsular India.     

Dynamic estimate of seismic danger based on multifractal properties of low-frequency seismic noise

A new method of dynamic estimate of seismic danger is presented which is based on estimating multifractal properties of low-frequency seismic noise. The efficiency of the method is illustrated by the analysis of seismic noise from broadband seismic network F-net in Japan. The analysis of multifractal properties of low-frequency seismic noise from Japan seismic network F-net since the beginning of 1997 allowed a hypothesis about approaching Japan Islands to a future seismic catastrophe to be formulated at the middle of 2008. The base for such a hypothesis was statistically significant decreasing of multi-fractal singularity spectrum support width mean value. The peculiarities of correlation coefficient estimate within 1 year time window between median values of singularity spectra support width and generalized Hurst exponent allowed to make a decision that starting from July 2010, Japan come to the state of waiting strong earthquake. This prediction of Tohoku mega-earthquake, initially with estimate of lower magnitude as 8.3 only (at the middle of 2008) and further on with estimate of the time beginning of waiting earthquake (from the middle of 2010), was published in advance in a number of scientific articles and abstracts on international conferences. The analysis of seismic noise data after Tohoku mega-earthquake indicates increasing probability of the 2nd strong earthquake within the region where the north part of Philippine Sea plate is approaching island Honshu (Nankai Trough). This region is characterized by relatively low values of singularity spectrum support width which is an indicator of seismic danger. In one paper (Sobolev in Izv Phys Solid Earth 47:1034–1044, 2011), the low-frequency seismic noise at the same range of periods was investigated retrospectively using data from the stations of broadband network IRIS which are located around the epicenter of Tohoku mega-earthquake with a distance up to 1,200 km. It was shown that the variance of the noise and the number of high-amplitude asymmetric impulses were grown dramatically before the event for stations which are located within the radius up to 500 km from the epicenter.     

Continued seismic activity at the Koyna reservoir site, India

Following the impounding of the Shivaji Sager Lake in 1962, tremors became prevalent in the Koyna region, considered previously to be aseismic. During ensuing years the tremor frequency appears to have been dependent on the rate of increase of water level, maximum water level reached, and the period for which high levels were retained. This culminated in a burst of seismic activity from September 1967 to January 1968 following the record water levels in the reservoir and included the earthquake of September 13, 1967 with magnitude 5.5 and the damaging December 10, 1967 earthquake of magnitude 6.0. During the next five years water levels were kept low and no significant earthquakes occurred subsequent to the October 29, 1968 earthquake of magnitude 5.

The reservoir was filled to maximum capacity during September 1973 and this was followed by a conspicuous increase in seismic activity which included an earthquake of magnitude 5.1 on October 17, 1973. However, seismic activity during 1973 was much less severe than that of 1967. This relative decrease in seismicity may indicate that (a) the “threshhold level” for relatively large magnitude earthquakes had increased; (b) a major portion of the accumulated strains had been released; and/or (c) the importance of the longer period of high loading in 1967. Similar observations have been made at other seismically active reservoir sites.     

Comprehensive seismic hazard assessment of Tripura and Mizoram states

Northeast India is one of the most highly seismically active regions in the world with more than seven earthquakes on an average per year of magnitude 5.0 and above. Reliable seismic hazard assessment could provide the necessary design inputs for earthquake resistant design of structures in this region. In this study, deterministic as well as probabilistic methods have been attempted for seismic hazard assessment of Tripura and Mizoram states at bedrock level condition. An updated earthquake catalogue was collected from various national and international seismological agencies for the period from 1731 to 2011. The homogenization, declustering and data completeness analysis of events have been carried out before hazard evaluation. Seismicity parameters have been estimated using G–R relationship for each source zone. Based on the seismicity, tectonic features and fault rupture mechanism, this region was divided into six major subzones. Region specific correlations were used for magnitude conversion for homogenization of earthquake size. Ground motion equations (Atkinson and Boore 2003; Gupta 2010) were validated with the observed PGA (peak ground acceleration) values before use in the hazard evaluation. In this study, the hazard is estimated using linear sources, identified in and around the study area. Results are presented in the form of PGA using both DSHA (deterministic seismic hazard analysis) and PSHA (probabilistic seismic hazard analysis) with 2 and 10% probability of exceedance in 50 years, and spectral acceleration (T = 0. 2 s, 1.0 s) for both the states (2% probability of exceedance in 50 years). The results are important to provide inputs for planning risk reduction strategies, for developing risk acceptance criteria and financial analysis for possible damages in the study area with a comprehensive analysis and higher resolution hazard mapping.     

Recent evaluation of the assessment seismic hazards for Nuweiba, Gulf of Aqaba

The Gulf of Aqaba is considered seismically as one of the most active zones of the Dead Sea Transform region. The main shock of the 1995 Gulf of Aqaba earthquake sequence is considered as the largest shock in the (surface wave magnitude Ms?=?7.2) since the sixteenth century. The present study is a trial to detect the probabilistic seismic hazard analysis (PSHA) for Nuweiba site. Data used for this study was a combination of both historical and recent instrumental data. Results of the hazard assessment, expressed as in the worst case scenario, reveal that Nuweiba is exposed to the occurrence of a maximum credible earthquake of magnitude $ m_{{\max }} ~ = ~7.4 \pm 0.31 $ , at hypocentral distance of 15.6?±?10 km. For structure with the return period of 100 years, with a 90% probability of exceedance, the maximum expected earthquake magnitude (M_L) is 5.9 in this lifetime. The possibility of the maximum peak ground acceleration at the Nuweiba site is 0.41 g. Results of the hazard assessment can be used as an input data to assess the seismic risk for site of interest.     

Probabilistic seismic hazard assessment in the Constantine region,Northeast of Algeria

This paper presents a seismic hazard evaluation and develops an earthquake catalogue for the Constantine region over the period from 1357 to 2014. The study contributes to the improvement of seismic risk management by evaluating the seismic hazards in Northeast Algeria. A regional seismicity analysis was conducted based on reliable earthquake data obtained from various agencies (CRAAG, IGN, USGS and ISC). All magnitudes (M _l, m _b) and intensities (I ₀, I _MM, I _MSK and I _EMS) were converted to M _s magnitudes using the appropriate relationships. Earthquake hazard maps were created for the Constantine region. These maps were estimated in terms of spectral acceleration (SA) at periods of 0.1, 0.2, 0.5, 0.7, 0.9, 1.0, 1.5 and 2.0 s. Five seismogenic zones are proposed. This new method differs from the conventional method because it incorporates earthquake magnitude uncertainty and mixed datasets containing large historical events and recent data. The method can be used to estimate the b value of the Gutenberg-Richter relationship, annual activity rate λ(M) of an event and maximum possible magnitude M _max using incomplete and heterogeneous data files. In addition, an earthquake is considered a Poisson with an annual activity rate λ and with a doubly truncated exponential earthquake magnitude distribution. Map of seismic hazard and an earthquake catalogue, graphs and maps were created using geographic information systems (GIS), the Z-map code version 6 and Crisis software 2012.     

A conception of casualty control based seismic design for buildings

Earthquake disaster is still the number one among all natural disasters, particularly, in terms of destructive power in causing deaths. Can earthquake engineers control seismic casualties through the seismic design of buildings? For this purpose, a conception of casualty control based seismic design is presented and a “two-step decision-making” method is proposed for determining the optimum seismic design intensity (or ground-motion) for controlling both seismic death and economic losses. The key problems in establishing the model are to determine the appropriate socially acceptable level of earthquake mortality and establish the corresponding objective function and /or constraint conditions in determining the optimum seismic design intensity. Ten different grades of socially acceptable mortality are suggested and the final socially acceptable mortality level was proposed for seismic design based on a questionnaire that was distributed nationwide in China. Finally, the method was applied to eight cities with different seismic hazard in China and the effects of various grades of acceptable earthquake mortality on seismic design intensity are analyzed.     

Finite element models to represent seismic activity of the Indian plate

Quantification of seismic activity is one of the most challenging problems faced by earthquake engineers in probabilistic seismic hazard analysis. Currently, this problem has been attempted using empirical approaches which are based on the regional earthquake recurrence relations from the available earthquake catalogue. However, at a specified site of engineering interest, these empirical models are associated with large number of uncertainties due to lack of sufficient data. Due to these uncertainties, engineers need to develop mechanistic models to quantify seismic activity. A wide range of techniques for modeling continental plates provides useful insights on the mechanics of plates and their seismic activity. Among the different continental plates, the Indian plate experiences diffused seismicity. In India, although Himalaya is regarded as a plate boundary and active region, the seismicity database indicates that there are other regions in the Indian shield reporting sporadic seismic activity. It is expected that mechanistic models of Indian plate, based on finite element method, simulate stress fields that quantify the seismic potential of active regions in India. This article explores the development of a finite element model for Indian plate by observing the simulated stress field for various boundary conditions, geological and rheological conditions. The study observes that the magnitude and direction of stresses in the plate is sensitive to these conditions. The numerical analysis of the models shows that the simulated stress field represents the active seismic zones in India.