Modeling a tsunami from the Nicoya,Costa Rica,seismic gap and its potential impact in Puntarenas

Although subduction zones around the world are known to be the source of earthquakes and/or tsunamis, not all segments of these plate boundaries generate destructive earthquakes and catastrophic tsunamis. Costa Rica, in Central America, has subduction zones on both the Pacific and the Caribbean coasts and, even though large earthquakes (Mw = 7.4–7.8) occur in these convergent margins, they do not produce destructive tsunamis. The reason for this is that the seismogenic zones of the segments of the subduction zones that produce large earthquakes in Costa Rica are located beneath land (Nicoya peninsula, Osa peninsula and south of Limón) and not off shore as in most subduction zones around the world. To illustrate this particularity of Costa Rican subduction zones, we show in this work the case for the largest rupture area in Costa Rica (under the Nicoya peninsula), capable of producing Mw ～ 7.8 earthquakes, but the tsunamis it triggers are small and present little potential for damage even to the largest port city in Costa Rica.The Nicoya seismic gap, in NW Costa Rica, has passed its ～50-year interseismic period and therefore a large earthquake will have to occur there in the near future. The last large earthquake, in 1950 generated a tsunami which slightly affected the southwest coast of the Nicoya Peninsula. We present here a simulation to study the possible consequences that a tsunami generated by the next Nicoya earthquake could have for the city of Puntarenas. Puntarenas has a population of approximately eleven thousand people and is located on a 7.5 km long sand bar with a maximum height of 2 m above the mean sea level. This condition makes Puntarenas vulnerable to tsunamis.     

The Aguacaliente Fault,source of the Cartago 1910 destructive earthquake (Costa Rica)

On 4 May 1910, the most destructive earthquake in the history of Costa Rica (Ms 6.4) destroyed the city of Cartago, a major city located in the Valle Central of Costa Rica. Using both palaeo‐seismological and morphotectonic analyses, we have found evidence that points to the Aguacaliente Fault (AF) as the source of this earthquake. This structure is a N100° E trending, strike‐slip fault situated to the south of Cartago and within a wide band of deformation. We excavated two trenches near Bermejo, south of Cartago. We found evidence of three surface ruptures within the last 1000 years on this fault. The age of the most recent rupture is consistent with the Cartago 1910 earthquake. The AF is a seismogenic source capable of producing large earthquakes (Mw 6.5–6.9) with an estimated recurrence interval of about 500 years.     

Modeling of source parameters of the 15 December 2015 Deogarh earthquake of M<Subscript>w</Subscript> 4.0

We herein present source parameters and focal mechanism of a rare cratonic upper crustal earthquake of M_w4.0, which occurred at 8 km depth (centroid depth) below a region near Deogarh, Jharkhand. For our study, we used broadband waveform data from a seismic network of 15 three-component seismographs in the eastern Indian craton. The average seismic moment, moment magnitude and source radius are estimated to be 1.1 × 10¹⁵ N-m, 4.0 and 180.6 m, respectively. The high average stress drop of 14.27 MPa could be attributed to its lower-crustal origin. The mean corner frequency is calculated to be 4.1 Hz. To study the source mechanism, we perform a deviatoric constrained full waveform moment tensor inversion of multiple point sources on the band-passed (0.06 – 0.14 Hz) broadband displacement data of the Deogarh event, using ISOLA software. The best fit is obtained for the source at 8 km centroid depth, with a moment magnitude 3.7, and a right-lateral strike-slip mechanism with strike 162°, dip 72° and rake 169°. The P-axis orients N24°E, which is parallel to the direction of the absolute plate motion direction of Indian plate, while T-axis orients E-W, which is parallel to the strike of the pre-existing Damodar Graben (DG) of Gondwana age. The occurrence of this earthquake is attributed to the neotectonic reactivation of a fault associated with the E-W trending DG shear zone.     

Source parameters of 1<Superscript>st</Superscript> April 2015 Chamoli earthquake (Mw 4.8) vis-à-vis seismotectonics of the region

The paper presents a detailed analysis of 1^st April 2015 earthquake, whose epicenter (30.16° N, 79.28° E) was located near Simtoli village of Chamoli district, Uttarakhand. The focal depth is refined to 7 km by the grid search technique using moment tensor inversion. The source parameters of the earthquake as estimated by spectral analysis method suggested the source radius of ~1.0 km, seismic moment as 1.99E+23 dyne-cm with moment magnitude (Mw) of 4.8 and stress drop of 69 bar. The fault plane solution inferred using full waveform inversion indicated two nodal planes, the northeast dipping plane having strike 334° and dip 5° and the southwest dipping plane with dip 86° and strike 118°. The parallelism of the nodal plane striking 334° with dip 5° as indicated in depth cross sections of the tectonic elements suggested the north dipping Main Boundary Thrust (MBT) to be the causative fault for this earthquake. Spatio-temporal distribution of earthquakes during the period 1960-2015 showed seismic quiescence during 2006-2010 and migration of seismicity towards south.     

The seismogenic zone in the Central Costa Rican Pacific margin: high-quality hypocentres from an amphibious network

Transition from subduction of normal to thickened oceanic crust occurs in the central portion of the Costa Rican margin, where large interplate earthquakes (M ~ 7) and abundant interseismic seismicity have been associated with subduction of bathymetric highs. We relocated ~1,300 earthquakes recorded for 6 months by a combined on- and offshore seismological network using probabilistic earthquake relocation in a 3D P-wave velocity model. Most of the seismicity originated at the seismogenic zone of the plate boundary, appearing as an 18° dipping, planar cluster from 15 to 25–30 km depth, beneath the continental shelf. Several reverse focal mechanisms were resolved within the cluster. The upper limit of this interseismic interplate seismicity seems to be controlled primarily by the overlying-plate thickness and coherency, which in turn is governed by the erosional processes and fluid release and escape at temperatures lower than ~100 to 120 °C along the plate boundary. The downdip limit of the stick–slip behaviour collocates with relative low temperatures of ~150 to 200 °C, suggesting that it is controlled by serpentinization of the mantle wedge. The distribution of the interseismic interplate seismicity is locally modified by the presence of subducted seamounts at different depths. Unlike in northern Costa Rica, rupture of large earthquakes in the last two decades seems to coincide with the area defined by the interseismic interplate seismicity.     

Relocation of aftershocks, focal mechanisms and stress inversion: Implications toward the seismo-tectonics of the causative fault zone of Mw7.6 2001 Bhuj earthquake (India)

The HYPODD relocation of 1172 aftershocks, recorded on 8–17 three-component digital seismographs, delineate a distinct south dipping E–W trending aftershock zone extending up to 35 km depth, which involves a crustal volume of 40 km × 60 km × 35 km. The relocated focal depths delineate the presence of three fault segments and variation in the brittle–ductile transition depths amongst the individual faults as the earthquake foci in the both western and eastern ends are confined up to 28 km depth whilst in the central aftershock zone they are limited up to 35 km depth. The FPFIT focal mechanism solutions of 444 aftershocks (using 8–12 first motions) suggest that the focal mechanisms ranged between pure reverse and pure strike slip except some pure dip slip solutions. Stress inversion performed using the P and T axes of the selected focal mechanisms reveals an N181°E oriented maximum principal stress with a very shallow dip (= 14°). The stress inversions of different depth bins of the P and T axes of selected aftershocks suggest a heterogeneous stress regime at 0–30 km depth range with a dominant consistent N–S orientation of the P-axes over the aftershock zone, which could be attributed to the existence of varied nature and orientation of fractures and faults as revealed by the relocated aftershocks.     

Widespread fluid expulsion along the seafloor of the Costa Rica convergent margin

Active fluid and gas transport were measured and observed along more than 200 km of the convergent margin of Costa Rica during cruise SO144-2 aboard RV Sonne . Ten profiles were run with the TV-sled OFOS, eight of which detected the dense occurrence of cold vent sites. This discovery shows that seafloor fluid expulsion is widely spread along the Pacific margin of Costa Rica. Surficial evidence of fluid expulsion is indicated by the appearance of chemosynthetic vent organisms such as bacterial mats, vesicomyid, solemyid and mytilid bivalves and tubeworms. Numerous active vents were indicated by elevated methane concentrations (≤ 200 nmol L^–) in the bottom water. Although fluid-venting activity was known previously from a small area south of Nicoya Peninsula, the present study documents active seepage at landslides, headwall scarps related to seamount subduction, morphological intersections of faults and mid-slope mud volcanoes.     

The 20 May 2016 Petermann Ranges earthquake: centroid location,magnitude and focal mechanism from full waveform modelling

Ground velocity records of the 20 May 2016 Petermann Ranges earthquake are used to calculate its centroid-moment-tensor in the 3?D heterogeneous Earth model AuSREM. The global-centroid-moment-tensor reported a depth of 12?km, which is the shallowest allowed depth in the algorithm. Solutions from other global and local agencies indicate that the event occurred within the top 12?km of the crust, but the locations vary laterally by up to 100?km. We perform a centroid-moment-tensor inversion through a spatiotemporal grid search in 3?D allowing for time shifts around the origin time. Our 3?D grid encompasses the locations of all proposed global solutions. The inversion produces an ensemble of solutions that constrain the depth, lateral location of the centroid, and strike, dip and rake of the fault. The centroid location stands out with a clear peak in the correlation between real and synthetic data for a depth of 1?km at longitude 129.8° and latitude –25.6°. A collection of acceptable solutions at this centroid location, produced by different time shifts, constrain the fault strike to be 304?±?4° or 138?±?1°. The two nodal planes have dip angles of 64?±?5° and 26?±?4° and rake angles of 96?±?2° and 77?±?5°, respectively. The southwest-dipping nodal plane with the dip angle of 64° could be seen as part of a near vertical splay fault system at the end of the Woodroffe Thrust. The other nodal plane could be interpreted as a conjugate fault rupturing perpendicular to the splay structure. We speculate that the latter is more likely, since the hypocentres reported by several agencies, including the Geoscience Australia, as well as the majority of aftershocks are all located to the northeast of our preferred centroid location. Our best estimate for the moment magnitude of this event is 5.9. The optimum centroid is located on the 20?km surface rupture caused by the earthquake. Given the estimated magnitude, the long surface rupture requires only ～4?km of rupture down dip, which is in agreement with the shallow centroid depth we obtained.     

Mélange formation by subduction erosion: the case of the Osa mélange in southern Costa Rica

The convergent plate margin off the Osa peninsula in southern Costa Rica is characterized by the indentation of the Cocos ridge at 4–5 Ma. The indentation causes the uplift of the Osa mélange which we interpret to represent an exhumed major channel for the transport of tectonically eroded material down into the subduction zone. We present evidence that, similar to the Nicoya segment of the Costa Rica convergent margin, subduction erosion rather than accretion has been the dominant process along the plate boundary. The composition of the Osa mélange is dominated by tectonized material of the upper-plate Nicoya ophiolite complex (basalt, radiolarite, limestone). Strong deformation is concentrated in numerous discrete shear zones and produced the layered fabric of large rock volumes, which partly experienced temperatures > 200°C. We thus interpret the Osa mélange to be a product of subduction erosion at the base of the outer arc wedge structure.     

The ultramafic flora of the Santa Elena peninsula,Costa Rica: A biogeochemical reconnaissance

The Santa Elena peninsula in the northwest of Costa Rica protrudes about 30 km westwards into the Pacific Ocean, and measures about 8–16 km in a north–south direction. Several geological studies have been carried out since 1953, showing that much of the peninsula is made up of peridotite, cut by mafic dykes. Only one previous brief examination appears to have been made of the vegetation in relation to the composition of the soils. We present here the results of a survey of some soils and plants of the eastern part of the peridotite massif, in which 73 plant specimens representing 51 identified species were collected and analyzed. The soils sampled all showed extreme ultramafic characteristics: Fe 10–16%, Mg 4–16%, Ca 0.5–1.4%; Ni 3000–7500 mg/kg, Cr 1400–3650 mg/kg, Co 150–325 mg/kg. The plants collected include several from genera such as Arrabidaea, Chamaesyce, Helicteres, Hyptis, Lippia, Oxalis, Polygala, Turnera and Waltheria that are also represented on ultramafics elsewhere in the Americas (e.g. Cuba, Puerto Rico, Brazil). Few of the species appear to be endemic to Costa Rica or to the ultramafics of Santa Elena. None of the specimens collected exhibited hyperaccumulation of nickel, the highest Ni concentration being 275 mg/kg in Buchnera pusilla.     

The geologic history of the Caribbean-Cocos plate boundary with special reference to the Nicoya ophiolite complex (Costa Rica) and D.S.D.P. results (Legs 67 and 84 off Guatemala): A synthesis

The Pre-Upper Senonian basement of Costa Rica crops out in the Santa Elena and Nicoya peninsulas. From south to north and from base to top the basement includes: the Esperanza, Matapalo and Santa Elena units. The Esperanza unit is Albian-Santonian in age and consists mainly of pillow basalt and massive basalt flows. The Matapalo unit includes Callovian to Cenomanian radiolarite and includes massive basalt flows, basalt, and dolerite basement. The Santa Elena unit contains ultramafic and mafic rocks in which harzburgite is the major component. The most important tectonic features of the Nicoya Complex are the large Santa Elena and Matapalo nappes. Nappe emplacement was from north to south during upper Santonian time. The sedimentary cover of the Nicoya Complex comprises:

1. (1) the Campanian El Viejo Formation that consists of shallow-water sediments in the north (Santa Elena Peninsula) and the Campanian-Maastrichtian Sabana Grande Formation of deep-water origin in the South (Nicoya Peninsula);

2. (2) Paleocene strata indicating deposition in a deep-water environment comprises the Rivas, Las Palmas and Samara Formations;

3. (3) a post-upper Eocene (?) sequence that consists of the shallow-water Barra Honda and Montezuma Formations.

Two unconformities are significant geological features of the upper-Senonian to Tertiary history of Costa Rica. The lower one is at the base of the Sabana Grande Formation and marks a major change in the geologic conditions (basalt is scarce in the Campanian-Tertiary series); the upper unconformity at the base of the Barra Honda and Montezuma Formations is not as major as the lower one. During post-Campanian time, normal faulting occurred in two stages separated by a strong erosional phase.The geology of the landward slope of the adjacent Middle America Trench is outlined by interpreting multifold seismic reflection records off the west coast of Costa Rica and the DSDP Legs 67 and 84 transects off Guatemala. The western Caribbean plate boundary may have been under extensional stress for the last 75 m.y. The strong landward-dipping reflectors of the Middle America Trench landward slope off Guatemala could be equivalent to the on-land pre-Campanian overthrusts of Costa Rica. The available data are consistent with the Convergent Extensional margin concept.     

A Source Study of the 9 September 1998 (Mw 5.6) Castelluccio Earthquake in Southern Italy using Teleseismic and Strong Motion Data

Teleseismic and strong motion data are used to derive the source parameters of the September 9, 1998 Castelluccio earthquake (M5.6). Teleseismic body-wave modeling indicates normal faulting (Plane 1: strike = 328°, dip = 50°, rake = −75°; Plane 2: strike = 126°, dip = 42°, rake = −107°) along NW–SE striking planes. Both nodal planes of the computed focal mechanism are tested for their capability of reproducing the mainshock acceleration time histories at three strong motion stations. Synthetic accelerograms are estimated using the stochastic method for finite sources in combination with the H/V ratios technique for the incorporation of the site effect. Our preferred model, which provided the best fit between synthetic and observed waveforms and corresponding Fourier amplitude spectra, consists of a NE-dipping normal fault with dimensions 8 × 7 km. The rupture nucleation point, which is assumed to coincide with the hypocenter location, was confined to the southeasternmost, deepest part of the fault. Our results are in good agreement with the so far released information regarding the aftershock sequence of the examined event, as well as the general seismotectonic knowledge on the broader epicentral area.     

Mineralogy and geochemistry of El Dorado epithermal gold deposit, El Sauce district, central-northern Chile

Summary The El Dorado Au-Cu deposit is located in an extensive intra-caldera zone of hydrothermal alteration affecting Upper Cretaceous andesites of the Los Elquinos Formation at La Serena (≈ 29°47′S Lat., 70°43′W Long., Chile). Quartz-sulfide veins of economic potential are hosted by N25W and N20E fault structures associated with quartz-illite alteration (+supergene kaolinite). The main ore minerals in the deposit are pyrite, chalcopyrite ± fahlore (As/(As + Sb): 0.06−0.98), with electrum, sphalerite, galena, bournonite-seligmanite (As/(As + Sb): 0.21−0.31), marcasite, pyrrhotite being accessory phases. Electrum, with an Ag content between 32 and 37 at.%, occurs interstitial to pyrite aggregates or along pyrite fractures. Pyrite commonly exhibits chemical zonation with some zones up to 1.96 at.% As. Electron probe microanalyses of pyrite indicate that As-rich zones do not exhibit detectable Au values. Fluid inclusion microthermometry shows homogenization temperatures between 130 and 352 °C and salinities between 1.6 and 6.9 wt.% NaCl eq. Isotope data for quartz, ankerite and phyllosilicates and estimated temperatures show that δ¹⁸O and δD for the hydrothermal fluids were between 3 and 10‰ and between −95 and −75‰, respectively. These results suggest the mineralizing fluids were a mixture of meteoric and magmatic waters. An epithermal intermediate-sulfidation model is proposed for the formation of the El Dorado deposit. Author’s present address: J. Carrillo-Rosúa, Dpto. de Didáctica de las Ciencias Experimentales, Universidad de Granada, Campus de Cartuja, 18071, Granada, Spain     

Unravelling the complexity of Apenninic extensional fault systems: A review of the 2009 L'Aquila earthquake (Central Apennines,Italy)

The 2009 L'Aquila sequence activated a normal fault system 50 km long in the Central Apennines, composed of two main NW-trending faults 12–16 km long: the main high angle L'Aquila segment and the Campotosto listric fault.The M_W 6.1 L'Aquila mainshock nucleated on the Paganica fault at a depth of ∼8.6 km and cut through the upper crust producing coseismic surface slip of up to 10 cm observed along a strike length of ∼13 km. Analysis of historical seismicity and data collected in paleo-seismological trenches suggest that this event filled a >500-year gap. In contrast, the blind Campotosto listric fault is composed of different fault segments displaying abrupt changes in dip at a depth where major events nucleate suggesting a rheological and geometrical control on stress concentration.A foreshock sequence that started around 4 months before the L'Aquila mainshock activated the deepest portion of the Paganica fault and marked the onset of large variations in elastic properties of the crustal volume. The variations have been modelled in terms of dilatancy and diffusion processes, corroborating the hypothesis that fluids play a key role in the nucleation process of extensional faults in the crust.     

The Todagin Creek landslide of October 3, 2006, Northwest British Columbia, Canada

The Todagin Creek landslide is located at 57.61° N 129.98° W in Northwest British Columbia. A seismic station 90 km north of the landslide recorded the event at 1643 hours coordinated universal time (UTC; 0943 hours Pacific daylight time (PDT)) on October 3, 2006. The signal verifies the discovery and relative time bounds provided by a hunting party in the valley. The landslide initiated as a translational rock slide on sedimentary rock dipping down slope at 34° and striking parallel to the valley. The landslide transformed into a debris avalanche and had a total volume estimated at 4 Mm³. An elevation drop of 771 m along a planar length of 1,885 m resulted in a travel angle (fahrb?schung) of 21.3°. The narrowest part of the landslide through the transport zone is 345 m. The widest part of the divergent toe of the landslide reaches a width of 1,010 m. Landslide debris impounded a lake of approximately 32 ha and destroyed an additional 67 ha of forest. The impoundment took 7 to 10 days to fill, with muddied waters observed downstream on October 13. No clear linkage exists with precipitation and temperature records preceding the landslide, but strong diurnal temperature cycles occurred in the days prior to the event. The Todagin Creek area appears to have an affinity for large landslides with the deposits of three other landslides >5 Mm³ observed in the valley.     

Seismotectonic and seismological aspects of the Mostaganem (Western Algeria) May 22, 2014 (Mw 4.9) seismic event

On Thursday, 22 of May 2014, at 6 h 22 min 0.3.3 s (GMT?+?1) a moderate-sized earthquake struck the Mostaganem, Western Algeria, region. The main shock, recorded by many international and national seismological stations, was preceded by a foreshock, 3 hours before, on May 22, 2014 (Ml?=?4.1) at 3 h 57 min 41.4 s and followed by four well-felt aftershocks (M?>?3.0) that lasted about 1 year. The main shock did not cause loss of lives but serious panic among the population was reported. The main shock, however, caused cracks in walls and roofs, sometimes destroyed, the old non-engineered and precarious adobe dweller corresponding to I₀?=?VI–VII (Msk scale). We used accelerograph records to (i) determine the epicenter location (longitude?=?0.3537 E, latitude?=?35.8598 N, (ii) perform waveforms inversion to calculate the earthquake parameters. The obtained results are, respectively, the seismic moment (M₀)?=?2.71 E + 16, the Mw?=?4.9 and the focal depth?=?6 km. The obtained focal mechanism solution shows reverse faulting with small right lateral component with the following nodal plans: NP1, strike?=?193.5, dip?=?49.5, slip?=?57.6 and NP2, strike?=?57.8, dip?=?50, slip?=?122.1. On the other hand, the seismotectonic framework of the Dahra area exhibits a serie of NE-SW trending “en echelon” faulted folds that may be active as suggested by this study.     

Co‐seismic Faults and Geological Hazards and Incidence of Active Fault of Wenchuan Ms 8.0 Earthquake,Sichuan, China

Abstract: There are two co-seismic faults which developed when the Wenchuan earthquake happened. One occurred along the active fault zone in the central Longmen Mts. and the other in the front of Longmen Mts. The length of which is more than 270 km and about 80 km respectively. The co-seismic fault shows a reverse flexure belt with strike of N45°–60°E in the ground, which caused uplift at its northwest side and subsidence at the southeast. The fault face dips to the northwest with a dip angle ranging from 50° to 60°. The vertical offset of the co-seismic fault ranges 2.5–3.0 m along the Yingxiu-Beichuan co-seismic fault, and 1.5–1.1 m along the Doujiangyan-Hanwang fault. Movement of the co-seismic fault presents obvious segmented features along the active fault zone in central Longmen Mts. For instance, in the section from Yingxiu to Leigu town, thrust without evident slip occurred; while from Beichuan to Qingchuan, thrust and dextral strike-slip take place. Main movement along the front Longmen Mts. shows thrust without slip and segmented features. The area of earthquake intensity more than IX degree and the distribution of secondary geological hazards occurred along the hanging wall of co-seismic faults, and were consistent with the area of aftershock, and its width is less than 40km from co-seismic faults in the hanging wall. The secondary geological hazards, collapses, landslides, debris flows et al., concentrated in the hanging wall of co-seismic fault within 0–20 km from co-seismic fault.     

青藏高原羌塘中部冈玛日—菊花山地区大型逆冲推覆构造的基本特征及形成机制

青藏高原羌塘中部冈玛日地区发现冈玛日-菊花山大型新生代逆冲推覆构造,延伸长度超过200km,总体走向近东西,构成新生代戈木错盆地的北部主干边界。推覆构造在冈玛日一带出露最好,表现为一系列单冲型逆冲断层构成的叠瓦状构造,逆冲断面倾向北,倾角12~20°,逆冲方向190~200°。冈玛日-菊花山逆冲推覆构造的发现,对探讨高原隆升在羌塘地区的表现形式具有重要的科学意义,为探讨新生代戈木错盆地的演化和盆地油气资源远景评价提供了新的构造依据。     

Thermochronological constraints of the exhumation and uplift of the Sierra de Pie de Palo,NW Argentina

The Sierra de Pie de Palo located between 67°30′–68°30′ W and 31°00′–32°00′ S in the Argentine Western Sierras Pampeanas in Argentina is a distinct basement range, which lacks thermochronological data deciphering its exhumation and uplift history below 200 °C. Integrated cooling histories constrained by apatite fission-track data as well as (U–Th)/He measurements of zircon and apatite reveal that the structural evolution of this mountain range commenced during the Late Paleozoic and was mainly controlled by tectonically triggered erosion. Following further erosional controlled exhumation in a more or less extensional regime during the Mesozoic, the modern topography was generated by denudation in the Paleogene during the early stage of the Andean deformation, whereupon deformation propagated towards the west since the Late Mesozoic to Paleogene. This evolution is characterised by a total of 3.7–4.2 km vertical rock uplift and by 1.7–2.2 km exhumation with a rate of 0.03–0.04 mm/a within the Sierra de Pie de Palo since ca. 60 Ma. Onset of uplift of peak level is also referred to that time resulting in a less Pliocene amount of uplift than previously assumed.     

Tectonics and thermal structure of western Satpura,India

Increased seismicity and occurrences of hot springs having surface temperature of 36–58 °C are observed in the central part of India (74–81° E, 20–25° N), where the NE trending Middle Proterozoic Aravalli Mobile Belt meets the ENE trending Satpura Mobile Belt. Earlier Deep Seismic Sounding (DSS) studies along Thuadara-Sendhwa-Sindad profile in the area has showed Mesozoic Sediments up to around 4 km depth covered by Deccan Trap and the Moho depth with a boundary velocity (P_n) of 8.2 km/s. In the present study, surface heat flow of 48 ± 4 mW m^?2 has been estimated based on P_n velocity, which agrees with the value of heat flow of 52 ± 4 mW m^?2 based on Curie point isotherms estimates. The calculated temperature-depth profile shows temperature of 80–120 °C at the basement, which is equivalent to oil window temperature in Mesozoic sediments and around 570–635 °C at Moho depth of 38–43 km and the thermal lithosphere is about 110 km thick, which is comparatively higher than those of adjoining regions. The present study reveals the brittle–ductile transition zone at 14–41 km depth (temperature around 250–600 °C) where earthquake nucleation takes place.