AG_57.04.14_SURESH_finalonline:Layout 6 ANNALS OF GEOPHYSICS, 57, 4, 2014, S0431; doi:10.4401/ag-6540 S0431 Crustal velocity structure of the Deccan Volcanic Province, Indian Peninsula, from observed surface wave dispersion Gaddale Suresh1,*, Satbir S. Teotia2, Sankar N. Bhattacharya1 1 India Meteorological Department, Seismology Division, New Delhi, India 2 Kurukshetra University, Department of Geophysics, Kurukshetra, India ABSTRACT Through inversion of fundamental mode group velocities of Love and Rayleigh waves, we study the crustal and subcrustal structure across the central Deccan Volcanic Province (DVP), which is one of the world’s largest terrestrial flood basalts. Our analysis is based on broadband seis- mograms recorded at seismological station Bhopal (BHPL) in the central India from earthquakes located near west coast of India, with an average epicentral distance about 768 km. The recording station and epicentral zone are situated respectively on the northern and southern edges of DVP with wave paths across central DVP. The period of group velocity data ranges from 5 to 60 s for Rayleigh waves and 5 to 45 s for Love waves. Using the genetic algorithm, the observed data have been inverted to ob- tain the crust and subcrustal velocity structure along the wavepaths. Using this procedure, a similar velocity structure was also obtained ear- lier for the northwestern DVP, which is in the west of the present study region. Comparison of results show that the crustal thickness decreases westward from central DVP (39.6 km) to northwestern DVP (37.8 km) along with the decrease of thickness of upper crust; while the thickness of lower crust remains nearly same. From east to west S-wave velocity in the upper crust decreases by 2 to 3 per cent, while P-wave velocity in the whole crust and subcrust decreases by 3 to 6 per cent. The P- and S-wave velocities are positively correlated with crustal thickness and negatively correlated with earth’s heat flow. It appears that the elevated crustal and subcrustal temperature in the western side is the main factor for low ve- locities on this side. 1. Introduction The Deccan Volcanic Province (DVP) is the most remarkable continental flood basalts provinces of the Indian Peninsula and formed around 65 My, when this region passed over the Reunion hotspot during its northward movement [Duncan and Pyle 1988]. The DVP has emerged as most interesting subjects of re- search particularly for its enormous size. The west coast and Narmada-Son lineament are the probable source regions from which the Deccan lavas were poured out through tensional fractures [Sen 2001]. The lava flowed great distances towards east and south to form the largest terrestrial flood cover extending nearly half a million square km between latitudes 16-24°N and longitudes 70-77°E (Figure 1). The flood basalts are known as Deccan traps. There are several studies in the past to evaluate the lithospheric velocity structure of the Indian Penin- sula using body waves and surface waves from re- gional and teleseismic earthquakes. Dube et al. [1973] used body wave records of observatories situated in the Peninsular India from 40 earthquakes occurred in the Koyna region (Figure 1). Bhattacharya [1981] in- verted the surface wave group velocities along paths in the central India and obtained the model IP11. Bhat- tacharya [1991] found that IP11 also satisfies the sur- face wave group velocity data across the northwestern part of the Indian Peninsula. Furthermore, reviewing the surface wave studies in the Indian Peninsula, Bhat- tacharya [1992] noted that the short period (0.6 to 2.5 s) surface wave data across the Koyna region (in south- ern DVP) are satisfied by IP12 (Figure 2); IP12 is same as IP11 except the top 1 km of latter is replaced by a Deccan trap layer. Singh et al. [1999] obtained velocity structure using surface waves recorded at regional dis- tances from 1997 Jabalpur earthquake in the central India (Figure 2). From the same earthquake, Rao et al. [2002] used body waves at regional distances and re- fined the velocity model in a grid search method through ray tracing. Figure 2 summarizes the S-wave velocity structure of central India from these studies and shows that main difference is in the depth of the interface between layers of upper and lower crust. The crustal S-wave velocity model of Rao et al. [2002] agrees with that of Singh et al. [1999]; however, there are differences in Moho depth and P-wave velocities Article history Received April 10, 2014; accepted June 10, 2014. Subject classification: Deccan Volcanic Province, Surface waves, Crustal structure. between these two models. The Moho depth of Singh et al. [1999] is nearly same as that of IP11 or IP12 [Bhattacharya 1992]. The crustal structure of northwestern DVP was obtained by Prajapati et al. [2011] using observed sur- face wave dispersion of Rayleigh and Love waves along a few wavepaths across this region. The mean wavepath of this study is shown in Figure 1; we shall refer this re- gion crossed by the wavepaths in northwestern DVP and studied by Prajapati et al. [2011] as northwestern DVP. The authors also evaluated the crustal structure of the adjoining continental shelf. In the present study, we shall evaluate the velocity structure of crust and subcrust (a layer below the crust) of central DVP which is east of northwestern DVP (Figure 1). We shall note the changes in the structure of DVP as we shift west- ward from central to coastal DVP region. Study of these changes will improve our understanding of the geology, tectonics and geodynamics processes of DVP. 2. Data and measurement of group velocity To select surface wavepaths across central-DVP, we consider the broadband records of the station Bhopal (BHPL), which is located in the northern part of DVP and analyze earthquakes occurred in the Koyna region in the west coast of India (Figure 1). The station BHPL is located on hard exposed bed-rock and has low ambi- SURESH ET AL. 2 Figure 1. Major tectonic features of the Deccan Volcanic Province (DVP). Shaded area indicates the DVP. Wavepaths from the epicentres of Koyna earthquakes to Bhopal (BHPL) have been shown in solid lines (present study) and are through central DVP. The dotted line is the mean wave path considered for evaluating the structure of northwest DVP [Prajapati et al. 2011]. Figure 2. S-wave velocity structure of the Indian Peninsula in vari- ous studies: Singh-etal [Singh et al. 1999], IP12 [Bhattacharya 1992], Rao-etal [Rao et al. 2002], northwestern DVP [Prajapati et al. 2011]. Theoretical group velocity curves have been drawn in Figure 3 for these structures. IP12 is modified from IP11 [Bhattacharya 1981, 1992], whose top 1 km is replaced by a Deccan trap layer. 3 ent ground noise. We consider well recorded surface waves from earthquakes of magnitudes between 4.2 and 5.4 (Table 1). A three component velocity broad- band seismogram of BHPL for the April 17, 2006, Koyna earthquake (mb 4.4) is shown in Figure 3 indi- cating well recorded surface waves which are promi- nent compared to the body waves. The hypocentral parameters of the earthquakes are determined using SEISAN analysis software [Havskov and Ottemoeller 1999] utilizing the phase data of regional seismological observatories. All the events were located using IP12 velocity model [Bhattacharya 1992]. Table 1 lists the hypocentral parameters of earthquakes considered; the wavepaths from all these events to the recording station at Bhopal are through the central portion of DVP (Fig- ure 1). Bhopal (23.241°N, 77.424°E, 520 m) is a part of National Seismological Network of India Meteorolog- ical Department (IMD), equipped with three compo- nent broadband sensor, Nanometrics Trillium-240 (with flat 3 dB velocity response from 240 s to 35 Hz) and accelerometer TSA-400S (with flat 3 dB accelera- tion response from DC to 225 Hz) connected to the six channel 24 bit digitizer (Taurus with Trident) and in- ternal time offsets are adjusted by GPS receiver auto- matically. The station is connected to Central Receiving Station in Delhi and Tsunami Warning Center at IN- COIS, Hyderabad through VSAT communication fa- cilities for near real time transfer of continuous recorded data and to perform scheduled quality checks [Dattatrayam et al. 2009]. The station is also equipped with broadband seismometer STS-2 connected to Quanterra digitizer (Q680LVG) [Bhattacharya and Dat- tatrayam 2000]. The average epicentral distance from Koyna earthquakes to BHPL is about 768 km. The group velocities are obtained through fre- quency-time-analysis (FTAN) as described in Dzie- wonskei et al. [1969] and Bhattacharya [1981, 1983]. For each event the 3-component 100 SPS (samples per sec- ond) data of BHPL is (1) decimated to 5 SPS, (2) changed to ground motion using the seismograph response and (3) converted to vertical, radial and transverse compo- nents. The vertical component is used for Rayleigh wave and transverse component for Love wave. In FTAN, the seismogram is passed through narrow Gaussian filter centered at each period. The time of the peak of the envelope to the amplitudes of the filtered seismogram gives the group arrival time (tg) and the group velocity = Δ/tg, where Δ is the epicentral distance. Such filtering makes it easier to separate different modes. The error in hypocenter determination is minimized by locating the earthquakes using regional structure and using phase data from regional observatories [Bhattacharya and Dattatrayam 2000]. Consequently, the error in epicen- ter is less than 5 km and the error in origin time is within 1 s; considering average epicentral distance of 768 km, these two errors cause an error less than 0.03 km/s in group velocity measurement from a single event. However, these errors are reduced by consider- ing the mean of the group velocities from many events; CRUSTAL STRUCTURE OF DVP, INDIAN PENINSULA Event no. Date (yyyy/mm/dd) Origin time hh:mm:ss Longitude (°E) Latitude (°N) Depth (km) Magnitude Used for* 1 1998/02/11 01:08:49.1 73.965 17.321 33.0 4.2 (mb) LQ, LR 2 2000/09/05 00:32:43.1 73.790 17.332 10.0 5.4 (mb) LQ, LR 3 2005/03/14 09:43:49.8 73.873 17.288 26.8 4.7 (mb) LQ, LR 4 2006/04/17 16:39:59.4 73.805 17.083 11.8 4.4 (mb) LQ, LR 5 2007/11/24 10:57:47.0 73.744 17.141 9.0 4.8 (Mw) LQ, LR 6 2008/09/16 21:47:13.3 73.734 17.287 15.0 5.0 (Mw) LQ, LR 7 2009/11/14 13:03:32.8 73.792 17.102 14.8 4.7 (Mw) LQ, LR 8 2009/12/12 11:51:24.8 73.613 17.112 9.7 5.1 (Mw) LQ, LR Table 1. List of earthquakes used for evaluating group velocities. *LQ= Love wave, LR=Rayleigh wave. Figure 3. Broadband waveform of BHPL from the April 17, 2006, Koyna earthquake (Δ = 7.01°, mb = 4.4, Table 1, event no. 4). Wave- form starts at 16:44:55.78. Although the amplitudes of body waves are weak but Rayleigh and Love waves are strong. an estimate for the uncertainty may be seen from the standard deviations of the group velocities. The group velocities are for periods ranging from 5 to 60 s for fun- damental mode Rayleigh wave and from 5 to 45 s for fundamental mode Love waves (Table 2, Figures 4 and 5). Because of this limited distance and moderate mag- nitude of usable earthquakes, no reliable surface wave dispersion measurements could be possible for longer periods. The standard deviations are generally higher for Rayleigh wave than for Love wave. Perhaps this can be explained by the source mechanisms of the earth- quakes, which caused amplitudes of Love wave rela- tively stronger than those of Rayleigh wave at the azimuth of the station from the sources. 3. Comparison of observed group velocity data In Figure 4 the observed mean group velocity data is compared with the theoretical dispersion curves for the previously evaluated models (Figure 2) in the Indian Peninsula discussed in the Introduction. The theoretical curves have been drawn following Bhattacharya [1987]. These theoretical curves are close to the observed data between 25 and 50 s. For Love wave, the dispersion curve for IP12 agrees with the observed data even at lower periods while that for the model of Singh et al. [1999] is close to Rayleigh wave data in this period range. At lower periods (5-25 s), the observed group ve- locities across central DVP are significantly higher than those of the dispersion curves of northwestern DVP; this implies a relatively higher wave speed structure in the upper crust for central DVP. At periods above 25 s, the Rayleigh and Love wave group velocities are only slightly more to that of northwestern DVP, which shows no significant change in lower crust and upper mantle shear velocity structure between them. How- ever, this observation indicates that the crust and upper mantle velocities in central DVP are on average faster than those of northwestern DVP. 4. Inversion and lithospheric structure The observed dispersion is a nonlinear function of the medium parameters viz., shear-, compressional- wave velocities, density and thickness of each layer. In order to obtain these model parameters through joint inversion of Rayleigh and Love wave group velocities, we apply a nonlinear inversion through Genetic Algo- rithms (GA) [Lomax and Snieder 1995, Yamanaka and Ishida 1996, Suresh et al. 2008]. The GA is configured to SURESH ET AL. 4 Period (s) Rayleigh waves Love waves Group velocity (km/s) S.D. (km/s) Number of measurements Group velocity (km/s) S.D. (km/s) Number of measurements 5.000 3.064 0.034 8 3.356 0.024 8 5.741 3.074 0.051 8 3.365 0.017 8 6.591 3.105 0.042 8 3.407 0.026 8 7.568 3.156 0.037 8 3.428 0.022 8 8.689 3.123 0.035 8 3.461 0.025 8 9.976 3.143 0.035 8 3.469 0.024 8 11.454 3.127 0.069 8 3.458 0.021 8 13.151 3.102 0.046 8 3.468 0.033 8 15.100 3.079 0.051 8 3.458 0.041 8 17.337 3.102 0.051 8 3.484 0.019 8 19.905 3.143 0.105 8 3.529 0.042 8 22.854 3.168 0.143 8 3.569 0.077 8 26.240 3.363 0.078 8 3.629 0.050 8 30.128 3.484 0.106 8 3.673 0.051 8 34.592 3.636 0.068 7 3.769 0.102 7 39.716 3.797 0.136 4 3.998 0.113 5 45.600 3.900 0.150 3 4.079 0.132 5 52.356 3.989 0.173 3 60.113 4.097 0.092 3 Table 2. Observed group velocities across central DVP. 5 find acceptable solutions within a user defined search space. During this process, the solution at each gener- ation is improved with misfit between observed and theoretical group velocities becoming less or equal to the preceding generation. Here, the misfit { is consid- ered similar to Prajapati et al. [2011] and is given by where, di R is the absolute difference between the ob- served and theoretical group velocity of a given struc- ture, vi R is the standard deviation of the observed group velocity at period Ti (i = 1, 2, …, NR) for Rayleigh waves. dj L and vj L are the corresponding values at pe- riod Tj (j = 1, 2, …, NL) for Love waves. The misfit { gives a measure of closeness of the theoretical curve to the observed data with respect to standard deviation. A model is acceptable if { < 1 and in such a case theoret- d , ,max i R i R j L j L { d = vv = G CRUSTAL STRUCTURE OF DVP, INDIAN PENINSULA Figure 4. The theoretical group velocity curves for different structures (Figure 3) of the Indian Peninsula are compared with the present mean observed data across central DVP. Figure 5. Observed mean group velocity data (open triangles) with ± one standard deviation. The theoretical curve has been drawn for the structure of central DVP (Table 4). (1) ical curves of this model are within the limits of one standard deviation around each mean observation. We keep the search space extended down to the depth 100 km; this depth limit is due to limit of the longer periods of the observed data. While performing the inversion, we have evaluated the model parameters of the top 5 layers consisting of 4 layers of the crust and a layer of subcrust, whose base is fixed at 100 km depth as in Prajapati et al. [2011]. Below 100 km the structure is same as in IP11 [Bhattacharya 1981]. Such 4-layered crust has also been obtained by forward modeling of gravity data derived from surface measurements on the Indian Peninsula [Arora et al. 2012]. The 14 variables in our GA search are: (1-4) thickness of top 4 layers of the crust, (5-9) VS of the top 5 layers and (10-14) VP/VS of top 5 layers. Here, VP is P-wave velocity and VS is S-wave velocity. The search space of these 14 variables is wide around the model parameters of those of northwest- ern DVP (Table 3). The density in each layer has been kept constant because it has the least effect on the dis- persion curves and the density is considered in each layer as in northwestern DVP (Table 3). The GA begins with a random initial population of K random models within solution limits. We con- sider a population of size K = 60 at each generation, where we create a new population with two elite mem- bers, 46 members through crossover and 12 members by mutation; the elite members replace the worst mod- els in the current generation with the best individuals of previous generation, so that the best individuals are not lost [Yamanaka and Ishida 1996]. Each GA opera- tion gives best solution within 150 generations and the best solution is accepted. Such operation is repeated making a list of 40 accepted solutions and the misfit value of all these solutions are between 0.74 and 0.83. The mean and standard deviations (SD) for each of varying parameters of all the 40 accepted solutions are obtained (Table 4). The SD generally shows the error in the results. VP is obtained from VS and VP/VS for each solution; the mean and SD of VP are also obtained. This model with mean values is considered as the model for central DVP. The misfit for this model is 0.77; this mis- fit value shows that the theoretical curves are within the SD of the observed data. 5. Discussion Central DVP mostly lies in the central India. Com- paring present result with other studies for different areas in the central India [Bhattacharya 1981, Singh et al. 1999, Rao et al. 2002] (Figure 2), we find VS varies between 3.55 and 3.65 km/s in the upper crust and be- tween 3.85 and 3.92 km/s in the lower crust; in the sub- crust VS varies between 4.60 and 4.68 km/s. Thus in the central India, the variation of VS around mean value is 2 per cent in the upper crust and these variations are 1 per cent in lower crust and subcrust. There is large vari- SURESH ET AL. 6 Layer Depth of top layer (km) Thickness of layer (km) P-wave velocity (km/s) S-wave velocity (km/s) Density (gm/cm3) 1. Deccan trap 0.0 1.6 4.79 2.74 2.40 2. Shallow crust 1.6 2.5 5.78 3.29 2.50 3. Upper crust 4.1 12.7 6.12 3.54 2.60 4. Lower crust 16.8 21.0 6.58 3.89 3.01 5. Subcrust 37.8 62.2 8.09 4.57 3.30 Layer Thickness (km) VP (km/s) VS (km/s) Density (g/cm3) 1. Deccan Trap 0.9±0.1 4.97±0.20 2.79± 0.04 2.40 2. Shallow crust 3.3±0.3 6.07±0.07 3.38± 0.05 2.50 3. Upper crust 14.8±1.3 6.44±0.10 3.65± 0.02 2.60 4. Lower crust 20.6±1.3 6.83±0.20 3.87± 0.02 3.01 5. Subcrust 60.4±0.2 8.37±0.08 4.68± 0.01 3.30 Table 3. The crust and subcrustal structure for northwestern DVP down to 100 km [Prajapati et al. 2011]. Below 100 km the structure is same as that of IP11 [Bhattacharya 1981]. Table 4. Crustal and Subcrustal Structure for central DVP down to 100 km. Below 100 km the structure is same as that of IP11 [Bhattacharya 1981]. Thickness of crust = 39.6 ± 0.2 km. The value after ± shows corresponding standard deviation. 7 ation in the depth of lower boundary of the upper crust in different study areas of central India and this depth varies between 13.8 and 20.4 km; however, the Moho depth varies only between 38.7 and 40.2 km and the crustal thickness of 39.6 km obtained in this study lies in this range. Figure 6 shows the comparison of velocity struc- ture of central DVP (Table 4) and northwestern DVP (Table 3). The crustal thickness is 39.6 km at central DVP and it decreases westward becoming 37.8 km at northwestern DVP. The crustal thickness decreases along with decrease of thicknesses of shallow and upper crust. The thickness of lower crust is nearly con- stant from central DVP (20.6 km) to northwestern DVP (21.0 km) and further westward to the continental shelf (20.9 km). In the continental shelf, Prajapati et al. [2011] found radial anisotropy in the crust caused by exten- sional tectonic deformation; it was found that the mean of SH- and SV-wave velocities are close to S-wave ve- locity of the isotropic crust of the adjoining north- western DVP. Such radial anisotropy is not noted either in northwestern DVP or in central DVP. The observa- tion of eastward decrease of trap thickness from 1.6 to 0.9 km is in agreement with other studies [Hari- narayana et al. 2007]. From central DVP in the east to northwestern DVP in the west, VS decreases by 2 to 3 per cent in the shallow and upper crust. Westward decrease of VP is also seen by 3 to 6 per cent in all layers of the crust and subcrust. These differences are above the standard de- viations of the evaluated velocities (Table 4) and thus the changes noted are reliable. Heat flow of central DVP is low compared to northwestern DVP [Sen 2001, Verma 2003]. In the Peninsular India, Sharma et al. [1991] also noted a significant decrease of subcrustal VP in the regions of high heat flow; from this correlation he inferred that low subcrustal VP is associated with el- evated temperature of the subcrustal region. Pei et al. [2004] has shown by quantitative analysis of the ob- served results that subcrustal VP is positively correlated with crustal thickness and negatively correlated with the earth’s heat flow. Here we note such correlation with VP not only in the subcrustal regions but also in all layers of the crust (Tables 3 and 4). The past volcanic activity was more on the west than on central part of DVP [Sen 2001]; thus low VP in northwestern DVP also agrees with the past volcanic activity as found by Pei et al. [2004]. In continuation of the study of Sharma et al. [1991], we may infer that the elevated temperature of the crust and subcrust in northwestern DVP is the main factor for causing low velocity in this part of DVP. 6. Conclusions We make the following inferences from the pres- ent study across central DVP and previous study across northwestern DVP [Prajapati et al. 2011]. (1) The crustal thickness of central DVP is 39.6 km, which remains in the range between 38.7 and 40.2 km obtained in other areas of the central India. How- ever, the thickness of upper crust differs significantly at other areas of the central India. (2) The crustal thickness decreases westward from central DVP (39.6 km) to northwestern DVP (37.8 km) along with the decrease of thickness of the upper crust. However, the lower crust has nearly the same thickness. (3) In a previous study, crustal radial anisotropy was found in the continental shelf adjoining north- western DVP. Such anisotropy is not noted in the con- tinental part of DVP. (4) In the shallow and upper crust of DVP, VS de- creases by 2 to 3 per cent from east to west, while VP in the whole crust and subcrust decreases by 3 to 6 per cent. The variations of VP and VS are positively corre- lated with the crustal thickness and negatively corre- lated with the earth’s heat flow. It appears that the elevated crustal subcrustal temperature in the western side is the main factor for low velocities on this side. Acknowledgements. The first author is thankful to the Di- rector General of Meteorology and Head (Seismology), India Me- teorological Department for providing necessary facilities to carry out the study and permission to publish the paper. Comments from anonymous reviewer were useful in improving the manuscript. References Arora, K., V.M. Tiwari, B. Singh, D.C. Mishra and I. Grevemeyer (2012). Three dimensional lithospheric CRUSTAL STRUCTURE OF DVP, INDIAN PENINSULA Figure 6. Comparison of velocity structures of central DVP (solid thick line) and northwestern DVP (dashed line). For central DVP (present study), ± one standard deviation of the velocities and thick- ness are also shown by shaded areas. structure of the western continental margin of India constrained from gravity modelling: implica- tion for tectonic evolution, Geophys. J. Int., 190, 131- 150. Bhattacharya, S.N. (1981). Observation and inversion of surface wave group velocities across Central India, B. Seismol. Soc. Am., 71, 1489-1501. Bhattacharya, S.N. (1983). Higher order accuracy in multiple filter technique, B. Seismol. Soc. Am., 73, 1395-1406. Bhattacharya, S.N. (1987). Reduction of the deep layers in surface wave computation, Geophys. J. R. Astro. Soc., 88, 97-109. Bhattacharya, S.N. (1991). Surface wave and lithos- pheric structure across the northwestern part of the Indian Peninsula, Pure Appl. Geophys., 135, 53-59. Bhattacharya, S.N. (1992). Crustal and upper mantle structure of India from surface wave dispersion, Current Science, 62, 94-100. Bhattacharya, S.N., and R.S. Dattatrayam (2000). Re- cent advances in seismic instrumentation and data interpretation in India, Current Science, 79, 1347- 1358. Dattatrayam, R.S., A.K. Bhatnagar, G. Suresh, P.R. Baidya and J.L. Gautam (2009). Real time earth- quake monitoring for early warning of Tsunamis in Indian Ocean region, Mausam, Diamond Jubilee Volume, 253-264. Dube, R.K., J.C. Bhayana and H.M. Chaudhary (1973). Crustal structure of the Peninsular India, Pure Appl. Geophys., 109, 1718-1727. Duncan, R.A., and D.G. Pyle (1988). Rapid eruption of the Deccan Traps at the Cretaceous/Tertiary bound- ary, Nature, 333, 841-843. Dziewonskei, A., S. Bloch and M. Landisman (1969). A technique for the analysis of transient seismic sig- nals, B. Seismol. Soc. Am., 59, 427-444. Harinarayana, T., B.P.K. Patro, K. Veeraswamy, C. Manoj, K. Naganjaneyulu, D.N. Murthy and G. Virupakshi (2007). Regional geoelectric structure beneath Deccan Volcanic Province of the Indian sub- continent using magnetotellurics, Tectonophysics, 445, 66-80. Havskov, J., and L. Ottemoller (1999). SeisAn Earth- quake analysis software, Seismol. Res. Lett., 70, 532- 534. Lomax, A., and R. Snieder (1995). Finding sets of ac- ceptable solutions with a genetic algorithm with ap- plication to surface wave group dispersion in Europe, Gephys. Res. Lett., 21, 2617-2620. Pei, S.P., Z.H. Xu and S.Y. Wang (2004). Discussion on origin of Pn velocity variation in China and adjacent region, Acta Seismologica Sinica, 17, 1-11. Prajapati, S., G. Suresh and S.N. Bhattacharya (2011). Crustal Structure of the Northwest Deccan Vol- canic Province, India, and the Adjoining Continen- tal Shelf from Observed Surface-Wave Dispersion, B. Seismol. Soc. Am., 101, 1488-1495. Rao, P.N., T. Tsukuda, M. Kosuga, S.C. Bhatia and G. Suresh (2002). Deep lower crustal earthquakes in Central India: inferences from analysis of regional broadband data of the year 1997 May 21, Jabalpur earthquake, Geophys. J. Int., 148, 132-138. Sen, G. (2001). Generation of Deccan Trap magmas, Proc. Indian Acad. Sci. (Earth Planet. Sci.), 110, 409- 439. Sharma, S.R., A. Sundar, V.K. Rao and D.V. Ramana (1991). Surface heat flow and Pn velocity distribu- tion in Peninsular India, J. Geodynamics, 13, 67-76. Singh, S.K., R.S. Dattatrayam, N.M. Shapiro, P. Mandal, J.F. Pacheco and R.K. Midha (1999). Crustal and Upper Mantle Structure of Peninsular India and Source Parameters of the 21 May 1997, Jabalpur Earthquake (Mw = 5.8): Results from a New Re- gional Broadband Network, B. Seismol. Soc. Am., 89, 1631-1641. Suresh, G., S. Jain and S.N. Bhattacharya (2008). Litho- sphere of Indus block in the northwest Indian sub- continent through genetic algorithm inversion of surface-wave dispersion, B. Seismol. Soc. Am., 98, 1750-1755. Verma, S.K. (2003). Geology, geophysics and geody- namics of extensive KCR volcanism during Pro- terozoic time in India, Long abstracts of 8th Int. Kimberlite Conference, Victoria, Canada, FLA0217. Yamanaka, H., and H. Ishida (1996). Application of Ge- netic Algorithms to an Inversion of Surface-Wave Dispersion Data, B. Seismol. Soc. Am., 86, 436-444. *Corresponding author: Gaddale Suresh, India Meteorological Department, Seismology Division, New Delhi, India; email: gsureshimd@yahoo.com. © 2014 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. SURESH ET AL. 8 << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice