Integrated geophysical interpretation for the area located at the eastern part of Ismailia Canal, Greater Cairo, Egypt

The integrated geophysical interpretation for the different geophysical tools such as resistivity and gravity is usually used to define the structural elements, stratigraphic units, groundwater potentiality, and depth to the basement rocks. In the present work, gravity and resistivity data were utilized for detecting the groundwater aquifer and structural elements, as well as the upper and lower surfaces of the subsurface basaltic sheet in an area located at the eastern side of Ismailia Canal, northeastern Greater Cairo, Egypt. Two hundred and ten gravity stations were measured using an Autograv instrument through a grid pattern of 50?×?50 m. The different required corrections were carried out, such as drift, elevation, tide, and latitude corrections. The final corrected data represented by the Bouguer anomaly map were filtered using high- and low-pass filters into regional and residual gravity anomaly maps. The resulting residual gravity anomaly map was used for gravity modeling to calculate the depths to the upper and lower surfaces of the basaltic sheet. The resulting gravity models indicated that the depths to the upper surface of the basaltic sheet are ranged between 26 and 314 m, where the shallower depths were found around the southern and eastern parts. The depths to the lower surface of the basaltic sheet are varied from 86 to 338 m, and the thickness of the basaltic sheet is ranged from 24 to 127 m, where the biggest thicknesses were found around the southern and northern parts of the study area. Forty-two vertical electrical soundings (VES) were carried out using Schlumberger configuration with AB/2 spacings ranged from 1.5 to 500 m. 1D quantitative interpretation was carried out through manual and analytical interpretations. The VES data were also inverted assuming a 3D resistivity distribution. The results from the 3D resistivity inversion indicated that the subsurface section consists of sand, sandstone, and sandy–clays of Miocene deposits overlying the basalts. Such basaltic features (of Oligocene age) are underlain by Gabal Ahmar Formation of Oligocene deposits, which are composed of sand and sandstone. Therefore, two aquifers were deduced in the area. The first is the Miocene aquifer (shallower) and the other is the Oligocene aquifer (deeper).     

Subsurface structural trends of the offshore Nile Delta area, Egypt: evidences from gravity and magnetic data

Subsurface structural trends and tectonics affecting the offshore Nile Delta area, Egypt, have been studied through the interpretations of gravity and magnetic data. Reduced to the pole, regional–residual separation, Tilt derivative and Euler deconvolution techniques are applied for the processing and interpretations of the magnetic and gravity data. The average depth of the sedimentary cover, estimated from the two-dimensional power spectrum technique ranges between 8 km and 13 km. The interpretation of the gravity and magnetic data indicates that the study area is affected by many subsurface structural trends. The NW–SE is the major trend related to El-Temsah and Misfaq-Bardwil trend. The NE–SW direction is the second dominant trend, related to the Rosetta trend. Other trends defined through the interpretation of gravity and magnetic data include: the N–S direction, related to the Baltim fault trend, the E–W direction, related to the Neogene hinge line and the NNE–SSW related to the Gulf of Aqaba. Accessory trends include the ENE–WSW, WNW–ESE and finally the NNW–SSW.     

Exploration for iron ore in Agbado-Okudu,Kogi State,Nigeria

This report details the result of geophysical exploration for iron ore; which involved vertical magnetic intensity (?Z) and gravity measurements, to delineate the geometry and depth extent of the deposit and acquiring quantitative and qualitative information for pre-drilling purposes in Agbado-Okudu. It is located about 3 km from Jakura along Lokoja-Jakura marble quarry and within low latitude precambrian basement complex district of Kogi State, Nigeria. A total of 517 magnetic measurement points along 16 traverses and 330 gravity reading along 11 profiles on the deposit in northeast–southwest azimuth were undertaken. The magnetic and gravity data enhancement involved linear regression curve fitting and fast Fourier transform, which were used to construct residual magnetic (RM) and gravity (RG) anomalies, analytic signal amplitude, Euler deconvolution at varying spectral indices (SI), power spectrum, and source parameter image (SPI), using the submenu of Geosoft Oasis Montaj software. Interpretation of the RM and RG anomalies revealed a primary causative body which perfectly correlates the positive anomalies and iron ore deposit, in form of a horizontal or gently dipping dyke with strike length of 600 m and average width of 110–130 m, within the gneiss complex in the north and trending south of the area. A secondary causative body associated with the negative anomalies and inferred as a vertical/near vertical thin sheet striking northeast–southwest coincided with the granitic and quartzitic intrusion. The NW–SE and E–W lineament trend conformed Kibarian and Liberian orogeny cycles of generally known structural trends in Nigeria, which shows that the iron ore deposit is structurally controlled. Depths to sources were estimated within range ≤ 2–24 m and 37.5–60 m, regarded as shallow and relatively deep depths, respectively. Ten vertical boreholes ranging in depth between 50 and 100 m are recommended, five of which require a priority attention to ascertain the thickness of the primary causative body.     

Tectonic deformation and subsurface active fault trends of El-Faiyum province,Egypt as deduced by seismicity and potential field data

We utilized the seismicity and the potential field data to study the tectonic deformation and to delineate the seismically active subsurface tectonic trends of El-Faiyum area. To accomplish these goals, we analyzed and interpreted the seismicity data, the reduced to pole total magnetic intensity, and the Bouguer anomaly maps. We also used the spatial distribution of the recent seismic events and the focal mechanism to outline the local seismic zones that control the seismicity of the study area and to determine the sense of the motion along the subsurface active faults. The focal mechanism of the recent seismic events and the interpreted subsurface tectonic faults from the potential field data reflect strike-slip movements with normal components along the subsurface active faults. The gravity and magnetic maps show a NE-SW regional trend with low gravity and magnetic values. The NE-SW regional trend extends across the whole area and could be related to the Pelusium Megashear fault. A NE-SW trend with high gravity and magnetic exists at the northern part and coincides with the Kattania Uplift and the basaltic flows in Gabal Qatrani area. The gravity and magnetic maps also reveal several local anomalies with different polarities, amplitudes, and extensions, which reflect anticlinal and synclinal structures on the basement surface. The seismotectonic map, generated by linking the basement structure map and the spatial distribution of the recent earthquake foci, reveals the dominant tectonic trends and the subsurface active faults.     

Surface magnetic survey of the Souss Basin, southwestern Morocco: Evaluation of the tectonic role postulated for the Agadir and Tarfaya fault zones and the South Atlas flexure

The results of a surface magnetic survey over the Souss Basin, the High Atlas and Anti-Atlas margins, combined with the gravity data of Van den Bosch (1971) and a review of the regional geology permits us to trace structures from the High Atlas into the Souss Basin. The South Atlas flexure is an east—west trending structure imposed during the Cenozoic upon an area with a northeast—southwest, Anti-Atlas structural grain.The major structure is the Tarfaya fault zone, whose location, projected by Choubert and Faure-Muret (1970), lies within a corridor of low residual magnetic anomalies. The continuation of this structure in the High Atlas is the Tizi n'Test fault zone. The Hercynian activity of the Tizi n' Test fault zone (Mattauer et al., 1972) was just one period of reactivation of an older, Precambrian, structure for we believe the Tarfaya fault zone was also active in post-late Triassic time.The Agadir earthquake occurred along a minor fault parallel to the Tarfaya-Tizi n' Test zone and is not related to the South Atlas flexure. Neither does it represent the entry into the African continent of a major ocean linked structure. The Tarfaya-Tizi n' Test fault zone is a much more likely candidate, for its offshore trend appears to match that of the Canary Fracture Zone.     

浑河断裂地球物理场特征及地震活动性

浑河断裂是辽宁省内一条重要的北东向地震构造带,断裂两侧的地球物理重磁场存在差异。通过浑河断裂的区域重力场与磁场研究,认为浑河断裂既控制了深、浅部岩石地层的发育,同时浑河断裂的重力场、磁场分布随着深、浅部地层岩性发育的差异而有所变化。浑河断裂的重力场和磁场梯度带主要位于抚顺以西至沈阳一带,浑河断裂历史以来有记录的最大破坏性地震位于重力梯度变化较大的沈阳段。综合浑河断裂地球物理场特征及区域地质资料,浑河断裂虽然活动性一般,但抚顺地区的采煤活动对断裂的活动性有一定影响。近些年抚顺地区频发的沉降、滑坡和地裂缝等地质灾害应当引起对浑河断裂活动性研究的重视。     

芦山地震同震和震后地表形变及重力变化的理论模拟

为了研究芦山地震的孕震过程和震源区的长期构造过程以及解释实测的震后形变和重力资料, 采用分层介质模型, 利用数值模拟的方法, 考虑区域流变系数, 计算了地震引起的地表同震、震后的形变和重力变化以及区域内部分GPS与重力连续观测台站的震后形变和重力变化的时间序列.结果表明: 芦山地震的地表同震形变显示出发震断层明显的逆冲特性; 粘弹性松弛效应引起的震后地表形变和重力变化比同震形变和重力变化的范围明显扩大, 但随着粘滞系数的增加, 变化量明显减小; 观测台站的震后变化时变曲线显示震后形变和重力变化在震后50 a间变化显著, 100 a后基本平缓, 趋于稳定; 模拟计算的GPS台站中除了MEIG台和MYAN台以外, 其余台站的震后观测必须考虑粘弹性松弛的影响.      

Delineation the active subsurface structures at Fayoum-Cairo district, Egypt, using magnetic and magnetotelluric tools

In this work, we present a reconnaissance study to elucidate and delineate subsurface fault structures for an active tectonics area that lies between Cairo and El Fayoum provinces and consider major sources of seismicity in Egypt. Well logging, aeromagnetic, land magnetic, and magnetotelluric data have been used. The well-logging data were used for several drilled wells along W–E direction. The magnetic data were analyzed using trend analysis, 3D magnetic modeling, and Werner deconvolution techniques. The magnetotelluric data were interpreted using 2D (TM–TE) modeling techniques. The results show that there are eight major fault structures having E–W, N–S, and NW–SE directions. These faults extend downward for about 20 km at the Dahshour and Qatrani areas. The epicenter sources of the earthquakes are clustering around the intersections of these structures. The Kattaniya horst structure has been interpreted as a regional structure that exceeds the limits previously determined by geologists. The depth to this horst reaches about 1.8 km at the NW and more than 4.3 km at the southern parts. The interpreted values of magnetic susceptibility at the horst zone indicate that they are ultrabasic/basic intrusion bodies.     

Source mechanism and parameters of the 19 October 2012 earthquake,northern Egyptian continental margin

The 19 October 2012 earthquake (M _L = 5.1) occurred in the northern continental margin of Egypt within the Nile Cone at latitude 32.35° N and longitude 31.27° E. The quake was felt over a wide area in north Egypt and East Mediterranean countries, but no casualties have been reported. This area had experienced the large earthquake (Ms = 6.7) of 12 September 1955. The fault plane solution of the 19 October 2012 earthquake is here presented based on the digital seismograms recorded by the Egyptian National Seismological Network (ENSN) and other regional seismic stations. The analysis is carried out using the well-known techniques of first motion polarities of P-wave and the amplitude ratios of P-, SH-, and SV-waves with lower hemisphere projection. The fault plane solution based on the first P-wave onset demonstrates a left lateral strike-slip faulting mechanism, while the solution based on both P-wave polarities and amplitude ratios of P-, SH-, and SV-waves reveals a reverse fault with strike-slip component trending NW–SE to NE–SW, in conformity with the N–S compression along the Hellenic Arc convergence zone. Following the Brune’s model, the source dynamic parameters for the 19 October 2012 earthquake are estimated as corner frequency = 1.47 Hz, fault radius = 0.7 km, stress drop = 22.1 MPa, seismic moment = 2.80E + 16 Nm, and moment magnitude M _w = 4.9. These parameters may provide important quantitative information for the seismic hazard assessment studies.     

Seismicity and major geologic structures of Tiran and Sanafir islands and their surroundings in the Red Sea

Tiran and Sanafir islands and their surrounding areas are very important due to their location within the Red Sea, which is in a triple junction among the African and Arabian plates and the Sinai Peninsula microplate. Consequently, this area should be studied from a geological point of view, particularly because there is a plan to construct King Salman’s bridge connecting Egypt and Saudi Arabia. Freely available potential field data, such as magnetic and gravity data, were integrated with seismological data from nearby seismic stations to understand the regional structure and seismic activity in the area. Potential field data were analysed using edge detection techniques (Tilt DeRivative and horizontal gradient) for qualitative interpretation and 2D inversion modelling for quantitative interpretation. Seismological data were analysed geostatistically to identify many epicentre locations and estimate the focal depths. The results of the potential field data analysis indicate that Tiran and Sanafir islands are located along a subsurface geological edge trending NW, parallel to the Red Sea. A similarity in potential field anomalies between both islands and the southern Sinai Peninsula indicates that these islands were separated from the Sinai Peninsula during the formation of the Gulf of Aqaba via Dead Sea Transform fault. The analysis of the seismic data indicates that a specific motion characterizes each focal depth solution. The seismic events are related to main structural trends NW–SE and NE–SW.     

汶川地震后龙门山断裂带活动特征

利用地震后2009 ～ 2011年GPS监测数据,获得了龙门山断裂带所在地区2009～2010年、2010～2011年以及2009～2011年GPS测站运动速度场,分析了区域地壳运动总体趋势及形变特征;通过分析龙门山断裂带北段、中段、南段横切剖面的测站运动速度变化,探讨了汶川地震后龙门山断裂带运动特征.分析表明:汶川地震前后,地壳运动总体趋势未变,作顺时旋转;断裂带西侧GPS测站运动速度变大,东侧运动速度变小;龙门山断裂带的断裂性质地震前后都为右旋走滑挤压,断裂带运动速率受汶川地震影响较大,震后运动速率较震前有显著的增加.龙门山断裂带震后各段次级断裂活动不同,中南段以前山断裂运动为主,其它各段以后山断裂运动为主.地震后龙门山断裂带表现出的运动特征主要与地震活动有关.受汶川地震的影响,区域动力学、运动学平衡被打破,龙门山断裂带东侧震后初期弹性回返,表现为低速反向运动.龙门山断裂带西侧震后松弛为拉张区,运动速度加大.地震对断裂带的影响不同,导致断裂带各段及次级断裂表现出不同的运动特征.     

郯庐断裂带中南段重磁特征与深部构造

郯庐断裂带是中国东部地球物理场的分界线。它在磁场上呈现出一条NNE向线性正磁异常带,它的重力场特征表现为规模较大的布格重力异常梯度带。在断裂带的两侧,其区域磁场、区域重力场无论是场值大小,还是异常规模以及各异常的展布形态等诸方面都具有明显的区别。地球物理特征揭示郯庐断裂带深部构造,为研究重点矿产分布规律及地震灾害预报具有重要的实际意义。     

The effect of uniaxial compression in the behaviour of remanent and induced magnetizations for basaltic rocks from the western Sinai Peninsula, Egypt

Seventy oriented basaltic samples were collected from six sites from the Wadi Budra and Farsh El Azraq areas of the west-central part of Sinai. Rock magnetic properties such as Curie temperatures and hysteresis parameters, as well as microscopic observations, point to magnetite as the main carrier of the remanent magnetization. The stress sensitivity of basaltic rocks is relatively low. High stress produces an increase in the remanent and induced magnetizations perpendicular to the applied stress axis and a decrease parallel to it. The change of magnetization during stress action ranges from 1.1% to 3.5% for a stress of 100?bar. The differential total magnetic intensity field with time (within 2?years) was observed through 80 magnetic observation points set up on both sites of the basaltic sheet at the studied area. The observed temporal variations of magnetization can be interpreted as stress loading parallel to the regional stress field in the order of 50?±?20?bar, according to the stress sensitivity of the precursor basalt.     

Dextral transpression and late Eocene magmatism in the trans-Himalayan Ladakh Batholith (North India): implications for tectono-magmatic evolution of the Indo-Eurasian collisional arc

The trans-Himalayan Ladakh batholith is a result of arc magmatism caused by the northward subduction of the Tethyan oceanic lithosphere below the edge of the Eurasian plate. The batholith dominantly consists of calc-alkaline I-type granitoids which are ferromagnetic in nature with the presence of magnetite as the principal carrier of magnetic susceptibility. The mesoscopic and magnetic fabric are concordant and generally vary from WNW–ESE to ENE–WSW for different intrusions of ferromagnetic granites in different parts of the batholith. Strike of magnetic fabric is roughly parallel with the regional trend of the Ladakh batholith in the present study area and is orthogonal to the direction of India-Eurasia collision. In Khardungla and Changla section, the magnetic fabric is distributed in a sigmoidal manner. It is inferred that this sigmoidal pattern is caused by shearing due to transpression induced by oblique convergence between the two plates. U–Pb zircon geochronology of a rhyolite from the southern parts of the batholith gives a crystallization age of 71.7 ± 0.6 Ma, coeval with ~68 Ma magmatism in the northern parts of the batholith. The central part of the batholith is characterized by S-type two-mica granites, which gives much younger age of magmatism at 35.5 ± 0.5 Ma. The magnetic fabric of these two-mica granites is at a high angle to the regional trend of the batholith. It is proposed that these two-mica granites were emplaced well after the cessation of subduction and arc magmatism, along fractures that developed perpendicular to the regional strike of the batholith due to shearing.     

Moho topography under central Greece and its compensation by Pn time-terms for the accurate location of hypocenters: The example of the Gulf of Corinth 1995 Aigion earthquake

In this paper we expand over the whole of central Greece, the Moho map centered on the Gulf of Corinth from tomographic inversion of PmP traveltime profile data recorded by several tens of temporary stations. Our approach is based on Pn, Moho refracted waves, from a large regional earthquake recorded by both temporary stations and the permanent Hellenic network. The Moho map shows the large Moho depth under the Hellenides belt. It also highlights the shallower Moho domain towards the Aegean Sea south and east of the Corinth Gulf. The domain of shallow Moho is limited along a NE–SW prolongation ahead of the North Anatolian Fault, from the North Aegean Trough to the western tip of the Gulf of Corinth towards the Gulf of Patras. The Pn time-terms provide corrections for the permanent stations that can be used together with the 1D velocity–depth model for a first-order compensation of lateral heterogeneity and contribute to the accurate and fast location of earthquake hypocenters. As a test we relocated the 1995 Aigion earthquake in this way, using only the sparse data of the permanent stations. Hypocentral coordinates then shift close to those derived by a dedicated dense array deployed after the earthquake, implying improvement of the routine location.     

Tsunami Hazard for the Auckland Region and Hauraki Gulf,New Zealand

The Hauraki Gulf is a semi-enclosed sea next to the largest population centre in New Zealand, the Auckland metropolitan region. The potential tsunami hazard is of concern to regional and local planners around the Hauraki Gulf. The Hauraki Gulf has recorded 11 tsunamis and one meteorological tsunami (rissaga) since 1840.The historical tsunami data are relatively sparse, particularly for the largest events in 1868 and 1883. Moreover, local sources may produce damaging tsunamis but none has occurred during recorded history. Therefore numerical modelling of potential tsunami events provides a powerful tool to obtain data for planning purposes. Three main scenarios have been identified for numerical modelling:1. A teletsunami event from an earthquake off the West Coast of South America. Historically this region has produced the largest teletsunamis in the Hauraki Gulf.2. A tsunami generated by a local earthquake along the Kerepehi Fault. This fault bisects the Gulf, has been active during the last century at the southern inland end, and is overlain by a considerable thickness of soft sediment that may amplify the seismic waves.3. A tsunami generated by a volcanic eruption within the Auckland Volcanic Field. This field has involved a series of mainly monogenetic basaltic eruptions over the last 140,000 years. Many of these eruptions have involved phreatomagmatic eruptions around the coastal margins, or within the shallow waters close to Auckland.     

Tectonic structures across the East African Rift based on the source parameters of the 20 May 1990 M7.2 Sudan earthquake

Earthquakes in Kenya are common along the Kenya Rift Valley because of the slow divergent movement of the rift and hydrothermal processes in the geothermal fields. This implies slow but continuous radiation of seismic energy, which relieves stress in the subsurface rocks. On the contrary, the NW-SE trending rift/fault zones such as the Aswa-Nyangia fault zone and the Muglad-Anza-Lamu rift zone are the likely sites of major earthquakes in Kenya and the East African region. These rift/fault zones have been the sites of a number of strong earthquakes in the past such as the M _w = 7.2 southern Sudan earthquake of 20 May 1990 and aftershocks of M _w = 6.5 and 7.1 on 24 May 1990, the 1937 M _s = 6.1 earthquake north of Lake Turkana close to the Kenya-Ethiopian border, and the 1913 M _s = 6.0 Turkana earthquake, among others. Source parameters of the 20 May 1990 southern Sudan earthquake show that this earthquake consists of only one event on a fault having strike, dip, and rake of 315°, 84°, and ?3°. The fault plane is characterized by a left-lateral strike slip fault mechanism. The focal depth for this earthquake is 12.1 km, seismic moment M _o = 7.65 × 10¹⁹ Nm, and moment magnitude, M _w = 7.19 (?7.2). The fault rupture started 15 s earlier and lasted for 17 s along a fault plane having dimensions of ?60 km × 40 km. The average fault dislocation is 1.1 m, and the stress drop, ?σ, is 1.63 MPa. The distribution of historical earthquakes (M _w ≥ 5) from southern Sudan through central Kenya generally shows a NW-SE alignment of epicenters. On a local scale in Kenya, the NW–SE alignment of epicenters is characterized by earthquakes of local magnitude M _l ≤ 4.0, except the 1928 Subukia earthquake (M _s = 6.9) in central Kenya. This NW–SE alignment of epicenters is consistent with the trend of the Aswa-Nyangia Fault Zone, from southern Sudan through central Kenya and further southwards into the Indian Ocean. We therefore conclude that the NW–SE trending rift/fault zones are sites of strong earthquakes likely to pose the greatest earthquake hazard in Kenya and the East African region in general.     

Basement-controlled faulting of Paleozoic strata in southern Ontario, Canada: new evidence from geophysical lineament mapping

Basement fault reactivation is now recognized as an important control on sedimentation and fault propagation in intracratonic basins. In southern Ontario, the basement consists of complexly structured mid-Proterozoic (ca. 1.2 Ga) crystalline rocks and metasedimentary rocks that are overlain by up to 1500 m of Paleozoic sedimentary strata. Reactivation of basement structures is suspected to control the location of Paleozoic fault and fracture systems, but evaluation has been hindered by a limited understanding of the regional structural characteristics of the buried basement. New aeromagnetic- and gravimetric-lineament mapping presented in this paper better resolves the location of basement discontinuities and provides further evidence for basement controls on the distribution of Paleozoic fault and fracture systems. Lineament mapping was facilitated by reprocessing and digital image enhancement (micro-levelling, regional residual separation, derivative filtering) of existing regional gravity and aeromagnetic datasets. Reprocessed images identify new details of the structural fabric of the basement below southern Ontario and delineate several previously unrecognized aeromagnetic and gravity lineaments and linear zones. Linear zones parallel the projected trends of mid-Proterozoic terrane boundaries identified by field mapping on the exposed shield to the north of the study area, and are interpreted as zones of shearing and basement faulting. Mapped aeromagnetic and gravity lineaments show similar trends to Paleozoic faults and fracture networks and broad zones of seismicity in southern Ontario. These new data support an ‘inheritance model’ for Paleozoic faulting, involving repeated reactivation and upward propagation of basement faults and fractures into overlying cover strata.     

Gravity and magnetic data analysis of block-35, Jiza' basin,Eastern Yemen

The Jiza' basin is located in the eastern part of Yemen, trending generally in the E–W direction. It is filled with Middle Jurassic to recent sediments, which increase in thickness approximately from 3,000 m to more than 9,000 m. In this study, block-35 of this sedimentary basin is selected to detect the major subsurface geological and structural features characterizing this basin and controlling its hydrocarbon potentials. To achieve these goals, the available detailed gravity and magnetic data, scale 1:100,000, were intensively subjected to different kinds of processing and interpretation steps. Also, the available seismic reflection sections and deep wells data were used to confirm the interpretation. The results indicated three average depth levels; 12.5, 2.4, and 0.65 km for the deep, intermediate, and shallow gravity sources and 5.1 and 0.65 km for the deep and shallow magnetic sources. Accordingly, the residual and regional anomaly maps were constructed. These maps revealed a number of high and low structures (horsts and grabens and half grabens), ranging in depth from 0.5 km to less than 4.5 km and trending mainly in the ENE, NW, and NE directions. However, the analytical signal for both gravity and magnetic data also showed locations, dimensions, and approximate depths of the shallow and near surface anomaly sources. The interpretation of the gravity and magnetic anomalies in the area indicated that the NW, NNW, ENE, and NE trends characterize the shallow to deep gravity anomaly sources; however, the NE, NW, and NNE trends characterize the magnetic anomaly sources, mainly the basement. Two-dimensional geologic models were also constructed for three long gravity anomaly profiles that confirmed and tied with the available deep wells data and previously interpreted seismic sections. These models show the basement surface and the overlying sedimentary section as well as the associated faults.     

Tectonics and crustal structures related to Bhuj earthquake of January 26, 2001: based on gravity and magnetic surveys constrained from seismic and seismological studies

Gravity and magnetic data of the Kachchh basin and surrounding regions have delineated major E–W and NW–SE oriented lineaments and faults, which are even extending up to plate boundaries in the north Arabian Sea and western boundary of the Indian plate, respectively. The epicentral zone of Bhuj earthquake and its aftershocks is located over the junction of Rann of Kachchh and median uplifts viz. Kachchh mainland and Wagad uplifts, which are separated by thrust faults. Gravity data with constraints from the results of the seismic studies along a profile suggest that the basement is uplifted towards the north along thrust faults dipping 40–60° south. Similarly gravity and magnetic modeling along a profile across Wagad uplift suggest south dipping (50–60°) basement contacts separating rocks of high susceptibility and density towards the north. One of these contacts coincides with the fault plane of the Bhuj earthquake as inferred from seismological studies and its projection on the surface coincides with the E–W oriented north Wagad thrust fault. A circular gravity high in contact with the fault in northern part of the Wagad uplift along with high amplitude magnetic anomaly suggests plug type mafic intrusive in this region. Several such gravity anomalies are observed over the island belt in the Rann of Kachchh indicating their association with mafic intrusions. The contact of these intrusives with the country rock demarcates shallow crustal inhomogeneities, which provides excellent sites for the accumulation of regional stress. A regional gravity anomaly map based on the concept of isostasy presents two centers of gravity lows of −11 to −13 mGal (10⁻⁵ m/s²) representing mass deficiency in the epicentral region. Their best-fit model constrained from the receiver function analysis and seismic refraction studies suggest crustal root of 7–8 km (deep crustal inhomogeneity) under them for a standard density contrast of −400 kg/m³. It is, therefore, suggested that significant amount of stress get concentrated in this region due to (a) buoyant crustal root, (b) regional stress due to plate tectonic forces, and (c) mafic intrusives as stress concentrators and the same might be responsible for the frequent and large magnitude earthquakes in this region including the Bhuj earthquake of January 26, 2001.