Geochemical and mineralogical studies of chlorine-rich amphibole and biotite from the 2.5 Ga mid-crustal basement beneath the 1993 Killari earthquake region,Maharashtra: Evidence for mantle metasomatism beneath the Deccan Trap

Fluid driven metasomatism and mass transfer from the earth’s mantle have played an important role in the evolution of the lower continental crust in many geodynamically active areas. The epicentral region of the disastrous 1993 Killari earthquake (M 6.2), concealed below a thick suite of Deccan volcanics in central India, appear to be one such region. In connection with the study of seismotectonics of the earthquake prone Deccan volcanic region, we have carried out systematic and detailed geochemical and mineralogical investigation on core samples from the basement, obtained from the 617m deep KLR-1 borehole, drilled in the epicentral region of Killari. Our investigations indicate that the basement, concealed below 338m thick Deccan volcanics, is made up of CO₂, Cl, FeO and CaO-rich, high density (2.82 g/cm³) — high velocity (avg. Vp: 6.2 km/s) moderately retrogressed upper amphibolite to granulite facies mid crustal rocks, which were subjected to pervasive Ca-metasomatism due to infiltration of mantle fluids. Graniticgneissic layer, typical of the upper crust, seems to be totally absent from this earthquake region. Chondrite normalized trace and rare earth elemental patterns display negative Eu anomalies together with LILE enrichment. Similarly, spider diagrams for incompatible elements show depletion in Zr, Hf, Y, Ta and Nb relative to the primitive mantle, indicating possible alterations of such relatively immobile elements at relatively high temperatures. Selective enrichment is also observed in transitional elements like Cu and Zn, indicating the possible role of chlorine in metal transport. The present study suggests that regional metasomatism beneath the Deccan Traps, which apparently alters the basic fabric of the rock during recrystallisation and makes it weak, may have a link with the nucleation of large earthquakes.     

Uplifted crust in parts of western India

During northward movement of the Indian sub-continent, after its breakup from the Gondwanaland in the Late Cretaceous, the western part of India traversed over the Reunion plume. The Saurashtra peninsula and the Cambay Basin are two important geological regions in this part. Two and half dimensional density models, based on the crustal seismic structure, were generated to establish a relationship between these two regions. These models indicate that the crust is 32–33 km thick in the eastern Saurashtra and the northern part of the Cambay Basin. The shallower crust is in a triangular region formed by the extension of the western limb of the Proterozoic Aravalli trend in Saurashtra, its eastern limb and the Narmada fault in the south. Compared to 36–37 km thick crust to the west and 38–40 km to the east of this region the crust in the above triangular region is uplifted by 4 to 6 km. This uplift took place either after the deposition of Mesozoic sediments or was concomitant with the rise of Reunion plume prior to the extrusion of the Deccan volcanics as the region was close to the axis of the plume.     

Structure,tectonics and thermal state of the Lithosphere beneath intraplate seismic region of Latur,central India: An appraisal

Recent surge in intraplate seismicity has led to detailed geological and geophysical investigations, covering different continental segments of India including seismogenic region of Latur. A synthesis of such data sets to understand the prevailing tectonic and thermal state of the Lithosphere beneath Latur region, that witnessed a large scale human loss due to 1993 seismic activity, has revealed shallow surfacing of denser deeper crustal segments which may have resulted due to ongoing active subsurface tectonic activity like uplift and erosion since geological past. Below this region, Moho temperature exceeds 500°C, heat flow input from the mantle is quite high (29–35 mW/m²) and the asthenosphere is shallow (∼100±10 km). It is suggested that stress generated by ongoing upliftment and related subcrustal thermal anomaly is concentrating in this denser and stronger mafic crust within which earthquakes tend to nucleate. In all likelihood, the seismic activity witnessed in the region may stem from the deep crustal/lithospheric dynamics rather than the role of fluids at the hypocentral depth.     

2-D Crustal thermal structure along Thuadara-Sindad DSS profile across Narmada-Son lineament,central India

Central India is traversed by a WSW-ENE trending Narmada-Son lineament (NSL) which is characterized by the presence of numerous hot springs, feeder dykes for Deccan Traps and seismicity all along its length. It is divided in two parts by the Barwani-Sukta Fault (BSF). To the west of this fault a graben exists, whereas to the east the basement is uplifted between Narmada North Fault (NNF) and Narmada South Fault (NSF). The present work deals with the 2-D thermal modeling to delineate the crustal thermal structure of the western part of NSL region along the Thuadara-Sindad Deep Seismic Sounding (DSS) profile which runs almost in the N-S direction across the NSL. Numerical results of the model reveal that the conductive surface heat flow value in the region under consideration varies between 45 and 47mW/m². Out of which 23mW/m² is the contribution from the mantle heat flow and the remaining from within the crust. The Curie depth is found to vary between 46 and 47 km and is in close agreement with the earlier reported Curie depth estimated from the analysis of MAGSAT data. The Moho temperature varies between 470 and 500°C. This study suggests that this western part of central Indian region is characterized by low mantle heat flow which in turn makes the lower crust brittle and amenable to the occurrence of deep focused earthquakes such as Satpura (1938) earthquake.     

Gravity anomalies and associated tectonic features over the Indian Peninsular Shield and adjoining ocean basins

A combined gravity map over the Indian Peninsular Shield (IPS) and adjoining oceans brings out well the inter-relationships between the older tectonic features of the continent and the adjoining younger oceanic features. The NW–SE, NE–SW and N–S Precambrian trends of the IPS are reflected in the structural trends of the Arabian Sea and the Bay of Bengal suggesting their probable reactivation. The Simple Bouguer anomaly map shows consistent increase in gravity value from the continent to the deep ocean basins, which is attributed to isostatic compensation due to variations in the crustal thickness. A crustal density model computed along a profile across this region suggests a thick crust of 35–40 km under the continent, which reduces to 22/20–24 km under the Bay of Bengal with thick sediments of 8–10 km underlain by crustal layers of density 2720 and 2900/2840 kg/m³. Large crustal thickness and trends of the gravity anomalies may suggest a transitional crust in the Bay of Bengal up to 150–200 km from the east coast. The crustal thickness under the Laxmi ridge and east of it in the Arabian Sea is 20 and 14 km, respectively, with 5–6 km thick Tertiary and Mesozoic sediments separated by a thin layer of Deccan Trap. Crustal layers of densities 2750 and 2950 kg/m³ underlie sediments. The crustal density model in this part of the Arabian Sea (east of Laxmi ridge) and the structural trends similar to the Indian Peninsular Shield suggest a continent–ocean transitional crust (COTC). The COTC may represent down dropped and submerged parts of the Indian crust evolved at the time of break-up along the west coast of India and passage of Reunion hotspot over India during late Cretaceous. The crustal model under this part also shows an underplated lower crust and a low density upper mantle, extending over the continent across the west coast of India, which appears to be related to the Deccan volcanism. The crustal thickness under the western Arabian Sea (west of the Laxmi ridge) reduces to 8–9 km with crustal layers of densities 2650 and 2870 kg/m³ representing an oceanic crust.     

Highly heterogeneous Precambrian basement under the central Deccan Traps, India: Direct evidence from xenoliths in dykes

Crustal or mantle xenoliths are not common in evolved, tholeiitic flood basalts that cover huge areas of the Precambrian shields. Yet, the occasional occurrences provide the most direct and unequivocal evidence on basement composition. Few xenolith occurrences are known from the Deccan Traps, India, and inferences about the Deccan basement have necessarily depended on geophysical studies and geochemistry of Deccan lavas and intrusions. Here, we report two basalt dykes (Rajmane and Talwade dykes) from the central Deccan Traps that are extremely rich in crustal xenoliths of great lithological variety (gneisses, quartzites, granite mylonite, felsic granulite, carbonate rock, tuff). Because the dykes are parallel and only 4 km apart, and only a few kilometres long, the xenoliths provide clear evidence for high small-scale lithological heterogeneity and strong tectonic deformation in the Precambrian Indian crust beneath. Measured ⁸⁷Sr/⁸⁶Sr ratios in the xenoliths range from 0.70935 (carbonate) to 0.78479 (granite mylonite). The Rajmane dyke sampled away from any of the xenoliths shows a present-day ⁸⁷Sr/⁸⁶Sr ratio of 0.70465 and initial (at 66 Ma) ratio of 0.70445. The dyke is subalkalic and fairly evolved (Mg No. = 44.1) and broadly similar in its Sr-isotopic and elemental composition to some of the lavas of the Mahabaleshwar Formation. The xenoliths are comparable lithologically and geochemically to basement rocks from the Archaean Dharwar craton forming much of southern India. As several lines of evidence suggest, the Dharwar craton may extend at least 350–400 km north under the Deccan lava cover. This is significant for Precambrian crustal evolution of India besides continental reconstructions.     

Crustal structure of the Gujarat region,India: New constraints from the analysis of teleseismic receiver functions

In this study, receiver function analysis is carried out at 32 broadband stations spread all over the Gujarat region, located in the western part of India to image the sedimentary structure and investigate the crustal composition for the entire region. The powerful Genetic Algorithm technique is applied to the receiver functions to derive S-velocity structure beneath each site. A detail image in terms of basement depths and Moho thickness for the entire Gujarat region is obtained for the first time. Gujarat comprises of three distinct regions: Kachchh, Saurashtra and Mainland. In Kachchh region, depth of the basement varies from around 1.5 km in the eastern part to 6 km in the western part and around 2–3 km in the northern part to 4–5 km in the southern part. In the Saurashtra region, there is not much variation in the depth of the basement and is between 3 km and 4 km. In Gujarat mainland part, the basement depth is 5–8 km in the Cambay basin and western edge of Narmada basin. In other parts of the mainland, it is 3–4 km. The depth of Moho beneath each site is obtained using stacking algorithm approach. The Moho is at shallower depth (26–30 km) in the western part of Kachchh region. In the eastern part and epicentral zone of the 2001 Bhuj earthquake, large variation in the Moho depths is noticed (36–46 km). In the Saurashtra region, the crust is more thick in the northern part. It varies from 36–38 km in the southern part to 42–44 km in the northern part. In the mainland region, the crust is more thick (40–44 km) in the northern and southern part and is shallow in Cambay and Narmada basins (32–36 km). The large variations of Poisson’s ratio across Gujarat region may be interpreted as heterogeneity in crustal composition. High values of σ (∼0.30) at many sites in Kachchh and few sites in Saurashtra and Mainland regions may be related to the existence of high-velocity lower crust with a mafic/ultramafic composition and, locally, to the presence of partial melt. The existing tectono-sedimentary models proposed by various researchers were also examined.     

Seismically constrained two-dimensional crustal thermal structure of the cambay basin

The temperature field within the crust is closely related to tectonic history as well as many other geological processes inside the earth. Therefore, knowledge of the crustal thermal structure of a region is of great importance for its tectonophysical studies. This work deals with the two-dimensional thermal modelling to delineate the crustal thermal structure along a 230 km long Deep Seismic Sounding (DSS) profile in the north Cambay basin. In this work P-wave velocities obtained from the DSS studies have been converted into heat generation values for the computation of temperature distribution. The model result reveals the Curie isotherm at a depth of ≈22 km and Moho temperature at around 900‡C.     

Anomalous crustal and lithospheric mantle structure of southern part of the Vindhyan Basin and its geodynamic implications

Tectonically active Vindhyan intracratonic basin situated in central India, forms one of the largest Proterozoic sedimentary basins of the world. Possibility of hydrocarbon occurrences in thick sediments of the southern part of this basin, has led to surge in geological and geophysical investigations by various agencies. An attempt to synthesize such multiparametric data in an integrated manner, has provided a new understanding to the prevailing crustal configuration, thermal regime and nature of its geodynamic evolution. Apparently, this region has been subjected to sustained uplift, erosion and magmatism followed by crustal extension, rifting and subsidence due to episodic thermal interaction of the crust with the hot underlying mantle. Almost 5–6 km thick sedimentation took place in the deep faulted Jabera Basin, either directly over the Bijawar/Mahakoshal group of mafic rocks or high velocity-high density exhumed middle part of the crust. Detailed gravity observations indicate further extension of the basin probably beyond NSL rift in the south. A high heat flow of about 78 mW/m² has also been estimated for this basin, which is characterized by extremely high Moho temperatures (exceeding 1000 °C) and mantle heat flow (56 mW/m²) besides a very thin lithospheric lid of only about 50 km. Many areas of this terrain are thickly underplated by infused magmas and from some segments, granitic–gneissic upper crust has either been completely eroded or now only a thin veneer of such rocks exists due to sustained exhumation of deep seated rocks. A 5–8 km thick retrogressed metasomatized zone, with significantly reduced velocities, has also been identified around mid to lower crustal transition.     

Magmatic underplating beneath the Rajmahal Traps: Gravity signature and derived 3-D configuration

The early Cretaceous thermal perturbation beneath the eastern continental margin of the Indian shield resulted in the eruption of the Rajmahal Traps. To understand the impact of the magmatic process that originated in the deep mantle on the lower crustal level of the eastern Indian shield and adjoining Bengal basin the conspicuous gravity anomalies observed over the region have been modelled integrating with available geophysical information. The 3-D gravity modelling has delineated 10–15 km thick high-density (ρ = 3.02 g/cm³) accreted igneous layer at the base of the crust beneath the Rajmahal Traps. Thickness of this layer varies from 16 km to the west of the Rajmahal towards north to about 12 km near Kharagpur towards south and about 18 km to the east of the Raniganj in the central part of the region. The greater thickness of the magmatic body beneath the central part of the region presents itself as the locus of the potential feeder channel for the Rajmahal Traps. It is suggested that the crustal accretion is the imprint of the mantle thermal perturbation, over which the eastern margin of the eastern Indian shield opened around 117 Ma ago. The nosing of the crustal accretion in the down south suggests the possible imprint of the subsequent magmatic intrusion along the plume path.     

Crustal accretion beneath the Koyna coastal region (India) and Late Cretaceous geodynamics

In 1967 a major earthquake in the Koyna region attracted attention to the hitherto considered stable Indian shield. The region is covered by a thick pile of Deccan lava flows and characterized by several hidden tectonic features and complex geophysical signatures. Although deep seismic sounding studies have provided vital information regarding the crustal structure of the Koyna region, much remains unknown. The two available DSS profiles in the region have been combined along the trend of Bouguer gravity anomalies. Unified 2-D density modelling of the Koyna crust/mantle suggests a ca. 3 km thick and 40 km wide high velocity/high density anomalous layer at the base of the crust along the coastline. The thickness of this anomalous layer decreases gradually towards the east and ahead of the Koyna gravity low the layer ceases to be visible. Based on the seismic and gravity data interpretation in the geodynamical/rheological boundary conditions the anomalous layer is attributed to igneous crustal accretion at the base of the crust. It is suggested that the underplated layer is the imprint of the magmatism caused by the deep mantle plume when the northward migrating Indian plate passed over the Reunion hotspot.     

Implications for the lithospheric structure of Cambay rift zone,western India:Inferences from a magnetotelluric study

Broad-band and long-period magnetotelluric(MT) data were acquired along an east-west trending traverse of nearly 200 km across the Kachchh,Cambay rift basins,and Aravalli-Delhi fold belt(ADFB),western India.The regional strike analysis of MT data indicated an approximate N59°E geoelectric strike direction under the traverse and it is in fair agreement with the predominant geological strike in the study area.The decomposed transverse electric(TE)-and transverse magnetic(TM)-data modes were inverted using a nonlinear conjugate gradient algorithm to image the electrical lithospheric structure across the Cambay rift basin and its surrounding regions.These studies show a thick(～1-5 km) layer of conductive Tertiary-Mesozoic sediments beneath the Kachchh and Cambay rift basins.The resistive blocks indicate presence of basic/ultrabasic volcanic intrusives,depleted mantle lithosphere,and different Precambrian structural units.The crustal conductor delineated within the ADFB indicates the presence of fluids within the fault zones,sulfide mineralization within polyphase metamorphic rocks,and/or Aravalli-Delhi sediments/metasediments.The observed conductive anomalies beneath the Cambay rift basin indicate the presence of basaltic underplating,volatile(CO_2,H_2 O) enriched melts and channelization of melt fractions/fluids into crustal depths that occurred due to plume-lithosphere interactions.The variations in electrical resistivity observed across the profile indicate that the impact of Reunion plume on lithospheric structures of the Cambay rift basin is more dominant at western continental margin of India(WCMI) and thus support the hypothesis proposed by Campbell Griffiths about the plume-lithosphere interactions.     

Pn wave velocity and Moho geometry in north eastern India

Pn velocity has been computed across the NE India and Moho geometry constrained, using regional earthquake travel times recorded by a network of 30 seismological stations operated during February-May 1993. Using an appropriate velocity model and the arrival times at the network stations, preliminary hypocentres of 16 regional earthquakes provided by NEIC were also improved. The average Pn wave velocity in NE India has been found to be 8.5 ±0.2 km/s. It varies from 8.3 to 8.5 km/s beneath the Shillong Plateau, Mikhir hills and Assam valley, which is significantly higher than those in other parts of India. The crustal thickness in NE India is also high, varying from 45–49 km under the Shillong plateau and the adjoining region to 55–65 km in the convergence zone. The presence of a thick crust and high Pn velocity suggests that the lithosphere in NE India is colder, as also indicated by the observed deeper level (45-51 km) seismicity of the region.     

Mantle-derived fluids in the basement of the Deccan Trap: evidence from stable carbon and oxygen isotopes of carbonates from the Killari borehole basement, Maharashtra, India

We report here for the first time geochemical, mineralogical and stable carbon and oxygen isotopic data on the crystalline basement rocks of the 1993 Killari earthquake region of Maharashtra (India), which is covered by a thick suite of Deccan volcanics. Our results revealed the hitherto unknown amphibolite–granulite nature of the 2.5?Ga basement, which contains 2.00–2.28?wt% of CO₂. The stable carbon (δ¹³C) and oxygen (δ¹⁸O) isotopic measurements on carbonates separated from two basement samples KIL-13 (440.5?m depth) and KIL-20 (499.6?m depth) collected from the KLR-1 borehole drilled in the epicentral region showed the respective values of ?6.23 and ?6.22‰ versus PDB for δ¹³C and 7.94 and 8.11‰ versus SMOW for δ¹⁸O. The samples plot in the primary igneous carbonatite field, indicating the mantle origin of the carbonates, derived through the process of mantle metasomatism from the deep mantle carbon reservoir. This would suggest large-scale crust-mantle thermal fluid interaction beneath the Killari seismogenic region, which is characterized by massive upwarping of the high-velocity mafic crust and retrograde metamorphism.     

南华北盆地群岩石圈热-流变结构

结合南华北盆地群现代地温场资料和深部地震测深资料及岩石热物性参数,对南华北盆地群的热结构进行了研究。结果表明：南华北盆地群平均热流值为53.7 mW/m2,地幔热流为30～34 mW/m2,莫霍面温度为500～550℃,热岩石圈厚度为110～130 km。在此基础上,进行了岩石圈流变模拟,探讨了研究区的岩石圈流变特征及其地球动力学意义。南华北盆地群岩石圈强度为（7.6~23.3）×1012 N/m,具有显著的 “三明治”结构。上地壳表现为脆性变形,中、下地壳为韧性的流动变形。这一分层变形机制决定了南华北盆地群的成盆演化动力学过程。     

Upwarped high velocity mafic crust, subsurface tectonics and causes of intraplate Latur-Killari (M 6.2) and Koyna (M 6.3) earthquakes, India – A comparative study

A unique attempt is made to understand the genesis of intraplate seismicity in the Latur-Killari and Koyna seismogenic regions of India, through derived crustal structure by synthesizing active and passive seismic, magnetotelluric, gravity and heat flow data. It has indicated presence of relatively high velocity/density intermediate granulite (and amphibolite) facies rocks underneath the Deccan volcanic cover caused mainly due to a continuous geodynamic process of uplift and erosion since Precambrian times. These findings have been independently confirmed by detailed borehole geological, geochemical and mineralogical investigations. The crystalline basement rock is found to contain 2 wt% of carbon-di-oxide fluid components. The presence of geodynamic process, associated with thermal anomalies at subcrustal depths, is supported by a high mantle heat flow (29–36 mW/m²) beneath both regions, although some structural and compositional variations may exist as evidenced by P- and S-wave seismic velocities. We suggest that the stress, caused by ongoing uplift and a high mantle heat flow is continuously accumulating in this denser and rheologically stronger mafic crust within which earthquakes tend to nucleate. These stresses appear to dominate over and above those generated by the India–Eurasia collision. The role of fluids in stress generation, as advocated through earlier studies, appears limited.     

Origin of the Felsic and Basaltic Dikes and Flows in the Rajula-Palitana-Sihor Area of the Deccan Traps,Saurashtra, India: A Geochemical and Geochronological Study

Rhyolite, trachyte, pitchstone, and granophyre dikes are associated with mafic dolerite dikes and basaltic flows of the northwestern part of the Deccan flood basalt province in the Saurashtra Peninsula, India. Felsic dikes, exposed in the Rajula area of Saurashtra, are similar in age to the basaltic flows of neighboring Palitana. The ages of both the felsic and mafic rocks straddle the ~65 Ma Cretaceous-Tertiary boundary and correspond to the main Deccan flood basalt episode. Palitana is centered on an elongated gravity high whose major axis is NE-SW, and Rajula is located on its southwestern flank. Unlike the younger Bombay felsic rocks from the western coast of India, which have been explained as partial melts of gabbros in deep crustal sills or previously erupted basalts, the incompatible-element characteristics of the Rajula rocks indicate that the Rajula rhyolites, trachytes, and dacites may have been generated by an almost complete melting of upper crustal rocks at the southwestern flank of the Rajula-Palitana-Sihor magmatic body. High potential temperatures of the Deccan plume, quick migration of the hot basaltic parent magma through lithospheric weak trends, and collection and residence of magma in upper-crustal magma chambers before eruption may have produced the right conditions to melt the upper crust in the vicinity of the Rajula-Palitana-Sihor magma chamber. On the other hand, the andesite located northeast of the magmatic body possibly evolved by assimilation of upper-crustal wall rocks accompanied by 5-10% crystallization of a Rajula-type basalt near the wall of the magma chamber. The Sihor rhyolites may also have been derived from the Sihor basalts through fractional crystallization accompanied by crustal assimilation. The Rajula granophyres, however, do not show any involvement of the upper crust in their genesis. These may have a history similar to that of the Bombay rocks and may have erupted in response to rifting along the Cambay rift.     

Assessment of predominant frequencies using ambient vibration in the Kachchh region of western India: implications for earthquake hazards

The Kachchh region is the second most seismically active region in India after the Himalaya. One of the disastrous Indian earthquakes of the millennium was the Bhuj earthquake of January 26, 2001, which caused about 14,000 casualties and huge property damage. The main reason for such devastation is due to lack of earthquake awareness and poor construction practices. Hence, an increase in the knowledge and awareness, based on improved seismic hazard assessment, is required to mitigate damage due to an earthquake. Natural predominant ground frequencies have been investigated in the Kachchh region of western India using ambient vibrations. The horizontal-to-vertical spectral ratio technique has been applied to estimate the predominant frequency at 126 sites. The ambient vibration measurements were conducted for about 1 h at each site in the continuous mode recording at 100 samples/s. We have validated the estimated predominant frequency with earthquake data recorded at six broadband stations in the region. It has been observed that geological time period has a significant effect on predominant frequency of the ground. The estimated predominant frequencies vary from 0.24 to 2.25 Hz for the Quaternary, 0.41–2.34 Hz for the Tertiary, 0.32–4.91 Hz for the Cretaceous, and 0.39–8.0 Hz for the Jurassic/Mesozoic. In the Deccan trap, it varies from 1.30 to 3.80 Hz. We found distinct variation of predominant frequencies of sites associated with hard rock and soft soil. The predominant frequencies were related to the thickness of the sediments, which are deduced by other geophysical and geological methods in the region. Our results suggest that frequencies of the region reveals the site characteristics that can be considered for studying the seismic risks to evolve a plan for disaster risk mitigation for the region.     

Temperature-depth profiles in Czechoslovakia and some adjacent areas derived from heat-flow measurements, deep seismic sounding and other geophysical data

The results of seismic measurements along three deep seismic sounding (DSS) profiles on the territory of Czechoslovakia and in adjacent countries have provided sufficient material about the crustal structure and the depth of the Moho discontinuity. These data, together with gravity and aeromagnetic data and the determinations of heat-flow values, were used to select several locations where the temperature—depth profiles were calculated. The Moho temperature of about 500 C beneath the Bohemian Massif increases to 800–1000 C and even more beneath the inner Neogene depressions of the Carpathian system. The regional differences in mantle heat-flow contribution between both these provinces may reach 1 μcal. cm⁻² sec⁻¹; such a variation in energy inflow may then be the driving force for the geological evolution. The geophysical implications of different thermal structure of the crust are discussed. Because of high subsurface temperatures in the Hungarian basin, partial melting at a depth of about 30 km may not be excluded.     

Crustal structure beneath Kottamiya broadband station,northern Egypt,from analysis of teleseismic receiver functions

We analyzed a total of 206 receiver functions beneath Kottamiya broadband station in northern Egypt to study the crustal structure and any azimuthal variations in the crustal thickness. The computed receiver functions are subdivided according to their azimuth into eight subgroups and analyzed separately using a genetic algorithm. The genetic algorithm is more appropriate than conventional linearized inversion schemes in regions where there is little a priori information about local crustal structures such as northern Egypt because it does not strongly depend on an initial model. The study region is located on the unstable shelf of Egypt in the northeastern corner of Africa. Little information about the deep structure of the crust beneath this region is available. For this reason, we have adopted the genetic algorithm to seismic waveform data recorded by Kottamiya broadband station. The crustal thickness varies slightly from 32 to 34 km with an average of 32.25 km, which is consistent with previous studies in the region. The crustal thickness shows a tendency of decrease toward the east and northeast being consistent with the general tectonic setting of the region including the opening of the Red Sea in the Tertiary times. Nonetheless, more teleseismic receiver functions from earthquakes recorded at denser seismic stations in northern Egypt and the southeastern Mediterranean combined with surface wave dispersion data as well as other geophysical investigations are necessary for more detailed imaging of the crustal structure which will deepen our understanding of the current tectonic and seismic activities of the region.