Open access peer-reviewed chapter

Quantitative Approximation of Geothermal Potential of Bakreswar Geothermal Area in Eastern India

Written By

Chiranjit Maji, Hirok Chaudhuri and Saroj Khutia

Submitted: January 25th, 2021 Reviewed: February 3rd, 2021 Published: March 9th, 2021

DOI: 10.5772/intechopen.96367

Chapter metrics overview

165 Chapter Downloads

View Full Metrics


Proper utilization of geothermal energy for power generation is still overlooked in India even after having enough potential as much as the equivalent to its other nonconventional energy resources. The source of geothermal energy is the decay of the radio-nuclei present inside the Earth’s crust apart from the primordial heat source. The noble gas 4He is also produced during the radioactive disintegration process. Therefore, measuring the amount of 4He gas along with some other geochemical parameters in an Indian geothermal area, the potential of the reservoir can be evaluated. Mathematical calculations relating to the radioactive disintegration to estimate the geothermal potential of Bakreswar geothermal reservoir utilizing the concept of the 4He exploration technique has been described here. The study showed that the heat (radiogenic) energy generated by the radioactive decay of 232Th, 238U, and 235U inside the reservoir was evaluated as 38 MW. This value raises to 76 MW when primordial heat is included. The detail calculations suggest that a Kalina cycle based binary power plant using ammonia–water mixture as working fluid is supposed to be installed at the identified locations with a drilling depth of about 1,100 m and the plant would be capable of delivering the power of 9.88 MW to 40.26 MW.


  • hot springs
  • radioactive disintegration
  • helium generation
  • geothermal power
  • geothermal power plant

1. Introduction

The origin of the geothermal energy is connected with the internal structure of the planet and the physiochemical processes occurring therein. According to the current knowledge, geothermal energy is unevenly distributed throughout the globe near the surface to the deep interior of the Earth [1, 2]. Depending upon the accessibility as well as the opportunities for the utilization of modern technology, many nations in the world are exploiting this natural energy resources for the commercial production of electric power [3, 4]. Geothermal energy hence, geothermal areas are generally defined through the parameter, geothermal gradient, which is the rate of the increment of the temperature profile of underneath bedrock of the Earth. The average (global) value of the geothermal gradient is typically 30 °C/km in the continental crust and 100 °C/km in the oceanic crust [1, 5]. However, in geothermal areas, its values are well above (>40 °C/km) the global average value [6]. It is so because of the magmatic intrusion. This intrusion is nothing but the molten magma, trapped within the Earth’s crust at a depth of 5–10 km beneath the surface. This may still in a fluid state or the process of solidification and releasing heat constantly [2, 7, 8]. According to the origin of geothermal energy, it is categorized into two. One was from a relic of the Earth’s accretion process, in which huge energy was trapped within the Earth’s interior (∼4.5 billion years ago) [7]. This one is named as the primordial heat source. Another one is the radiogenic heat source, which is produced by the natural decay process of long-lived radioisotopes such as 238U, 235U, 232Th, and 40K. These nuclei, of which the half-life (T12)are comparable to the age of our planet, are found with significant abundance within the crust of the geothermal areas [9, 10]. A considerable amount of heat is contributed from the natural radioactive decay process. Eq. (1) to Eq. (3) represent the physicochemical processes and the heat energy released from the naturally occurring radioactive disintegration in each of the complete decay chain [11, 12, 13, 14]. Moreover, it shows the produced crustal He (4He) atoms and neutrinos during each decay process.


Moreover, within the deep Earth, the production rate of He from 232Th and 238U [and 235U] radio-nuclei are encountered to be 2.43×10100atoms/m3/sand 1.03×108atoms/m3/s,respectively [10]. The fact of characteristics heat–helium coherence at any geothermal system (under the deep reservoir) is interpreted by such physicochemical processes [15]. This radiogenic heat, which is one of the main sources of the Earth’s internal heat, powers all geodynamic processes underneath [16]. Generally, geothermal heat is transferred from the aquifer (reservoirs) to the Earth’s surface by the conduction and convection process. Here, the geothermal fluid (meteoric water) acts as the carrier [1] and the radiogenic He, being highly diffusive gas, generated from the host mineral and mixes by diffusion with the fluid that circulates into the deep Earth [17]. The reservoir, which is nothing but a volume consisting of hot permeable rocks, is usually sandwiched by capping of impermeable rocks. And it is favourably connected to a recharge (surficial) area [18], from which geothermal fluids percolated to recharge the aquifer cyclically [18, 19]. The circulating fluids, to which the heat is transferred from the reservoir, escape through fracture and features from the deep reservoir and manifest through geysers, fumaroles, hot springs, etc. [8, 19]. Moreover, through diffusion and advection process, the radioactive inert gases (like 222Rn & 220Rn) including the stable and inert gases (such as He, Ar) are spontaneously migrating upward from the deep Earth to the superimposing atmosphere [20, 21, 22]. This process, known as ‘Earth degassing’, is non-uniform over space & time [23, 24]. The prominent signature of this degassing is generally noticed along active faults, fractures, oceanic ridges, geothermal fields, and even deep wells [25, 26, 27].

It is notable that geothermal energy sources are still overlooked in India for power generation even after the existence of a lot of potential resources, which are seen in twelve geothermal zones of the country [19]. However, several of them could be well utilized for the generation of power by means of developing geothermal power plants. For the sake of investigation, the hot spring site at Bakreswar in West Bengal, India, was selected as shown in Figure 1. Now, knowing the amount of by-product, He gas which is ultimately reaching the surface through the fracture, fissure and hot springs vents, etc., the associated heat energy (radiogenic) produced inside the reservoir can be estimated. The energy released per unit time from underneath bedrock at the study area was calculated by means of measuring the average amount of He emanated from Agni Kunda hot spring at Bakreswar. Here mainly the decay series of 238U, 235U, and 232Th were considered, and the amount of heat energy contributed due to each series was evaluated. Here the question may arise that each decay series [Eq. (1) to (3)] takes a long period (in geological time scale) to complete its disintegration process and release a certain amount of heat and He discretely. But, heat and He generated due to each series were utilized to calculate the amount of heat production at the said reservoir at a certain instant of time. However, He emanation at the study area shows stable activities for a long-time-interval (5 years), as established by [19]. Therefore, He generation is also stabilized for a long period, i.e., He generation due to every radioactive decay series and emanation of the said gas is in an equilibrium condition. Therefore, no He is being stored at the reservoir at the instant, and, therefore, the He emanation could be considered to be equal to the generation of the same due to the radioactive disintegration process.

Figure 1.

Location of the study area Bakreswar in the map of India (modified after [19]).


2. The study area: Bakreswar geothermal province

A cluster of seven1 hot springs is scattered over Bakreswar geothermal area within a confined zone of the surface area of about 3350 sq. m [19, 28]. The area, which is a geologically complex, heterogeneous, and extensively faulted region, is situated at the eastern end of the SONATA (Son–Narmada–Tapi) geothermal province (Figure 2, window a) [19, 28]. The area lies in the West Bengal Basin (WBB), the extension of the Chotanagpur Gneissic Complex [29]. Furthermore, it is linked with a 1.2-km-long shear zone, which is characteristic by 50 m wide breccia/cherty quartzite aligned through the almost north–south trend-line [26] (Figure 2, window b). The springs here are connected with the extinct Rajmahal volcanic activity (115 Ma), and hence, are associated with the Precambrian granitic rocks [30] (Figure 2, window b). The highly permeable and porous subsurface of the site is facilitated due to the presence of the brecciated, highly sheared, and mylonitized rock here [26]. The association of the study area with the eastern edges of two major fault systems (the ONGC fault and the SONATA fault) made the region to be in a stressed state [31, 32]. This region is characterized by a very high geothermal gradient (∼90 °C/km) and a high heat flow rate (∼ 230 mW/m2) [19, 33]. The presence of a high heat-conducting zone in this area is confirmed by electrical resistivity studies. This conducting zone, which starts at a depth of around 2.8 km and goes down up to a depth of 4 km [34], is supposed to act as the heat feeder to the fault system linked with the Bakreswar hot springs. It is to be noted that the crustal thickness at the study area is only 24 km, whereas the average of the same throughout the country is 38 km [29, 35]. Besides, the average density of the crustal substance here is relatively low. Therefore, inert volatiles like He and 222Rn gases can easily transmit to permeate through crustal constraints due to the presence of the thinner lithospheric overburden here. As a result, the spring and the soil gases here are dominated by the presence of high 222Rn and He flux [29]. High 222Rn and He gases are continuously conveyed and dispersed into the atmosphere via molecular diffusion and the formation of micro-bubbles at the hot spring vents. Here, temperature and He emanation profile of some sites of Bakreswar geothermal provinces are also tabulated in Table 1 for a reference to attain the brief geophysical properties of the study area. Moreover, the reservoir temperature of the geothermal system underneath Bakreswar was predicted to be 100 ± 5 °C (at ∼1 km depth) by [38]. The same was estimated to be in the range of 130 °C to 175 °C (by Na/K ratio) and 110 °C to 124 °C (by TSiO2) by [39]. Furthermore, the range of the reservoir temperature was also evaluated as 212 °C to 124 °C, 118 °C to 120 °C, and 126 °C to 130 °C by means of silica geothermometry by [32, 38, 40] respectively. The audio-magnetotelluric (AMT) studies of the sub-surface beneath the Bakreswar geothermal area were conducted by [41] (Figure 3). The rapid relaxation inversion (RRI) for both transverse-electric (TE) and transverse-magnetic (TM) modes was carried out to figure out the resistivity profile of the subsurface of the site. Here, the suitable locations for drilling for the installation of a future geothermal power plant at the study area were identified by the authors using the result of that AMT survey [personal communication] and the same is discussed later.

Figure 2.

(a) Regional and (b) Local geological maps of the study area (modified after [19,26]).

Sl. No.Test Site
(distance from Agni Kunda)
Sample typeTemperature (°C)He Conc. (vol %)References
1Bakreswar Agni Kunda (0 m)HSG69.01.72[19]
2Bakreswar Khar Kunda (16 m)HSG68.01.36[36]
3Bakreswar Bhairab Kunda (7 m)HSG62.01.12[28]
4Bakreswar Brahma Kunda (20 m)HSG46.01.26[28]
5Bakreswar Surya Kunda (18 m)HSG63.00.31[28]
6Bakreswar Reserve Tank (5 m)HSG52.00.91[28]
7PWD Bungalow at Bakreswar (987 m)SG (1 m depth)32.0 (Ambient)0.35[37]
8PWD Bungalow at Bakreswar (988 m)SG (3 m depth)31.0 (Ambient)0.02[37]
9PWD Bungalow at Bakreswar (990)AA (1 m height)33.0 (Ambient)0.05[29]
(10 km)
AA (1 m height)28.0 (Ambient)0.07[29]
11Mallarpur (43 km)BG (100 m depth)58.01.20RWA

Table 1.

Temperature and He emanation profile of some sites at Bakreswar geothermal province.

Note: HSG = Hot spring gas; SG = Soil gas; AA = Ambient air; BG = Borehole gas; Conc. = Concentration; RWA = Recent work by the authors.

Figure 3.

(A) AEW traverses (AMT sites) on the map of the study area; (B) 2D RRI (rapid relaxation inversion) along traverse AEW1; (C) 2D RRI along traverse AEW2; (D) 2D RRI along traverse AEW4 (modified after [41] [personal communication] and [42]).


3. Methodology

3.1 Experimental techniques

In view of continuous monitoring of gases emanated from the hot spring Agni Kunda at the spring site of Bakreswar, a field laboratory was established. In this regard, a giant inverted SS funnel was placed under hot water at Agni Kunda at a position where gas out flux was significantly high, to trap hot spring gases which were comprised of He, Ar, O2, N2, CH4, CO2, 222Rn, etc. A portable and programmable μ-GC (micro-gas chromatograph) CP 490 (make Agilent, Netherland) comprised of a μ-thermal conductivity detector, was utilized to detect the relative concentration of different gases present in the spring gas. Here, ultra-pure (>99.998 vol%) H2gas was utilized as the carrier gas for running the equipment. The entire measurement was carried out in-round-the clock (24 × 7) measurement fashion for a continuous five-year (1st August 2005 to 31st July 2010). The back-up power supply was maintained to keep a stable and continuous power supply in case of a power failure. The schematic diagram of the experimental set-up is illustrated in Figure 4. Further details of the above-mentioned experimental procedure are already described by [19]. Here, the average value of the He concentration (vol%) of 5 years of continuous measurement was adopted in our study. Moreover, the flow rate of the emanated gases from the spring vent was measured by means of collecting the spring gases in a gas container (5 litres) from the main channel of the incoming gas line. The gas collection procedure was kept running up to a certain time until its pressure makes equilibrium with the pressure (1.58 atm) at the spring vent underwater. This type of measurement was done once every month for a continuous five years, and the average value of those was considered as the final value of the gas flow rate under consideration. The ambient temperature of the study area was monitored for the same interval at the time of measurement of gas flow rates.

Figure 4.

Schematic diagram of the experimental set-up installed at Bakreswar (modified after [19]).

3.2 Analytical techniques

To move towards the desired direction for calculation, the following steps were adopted.

The volume of He gas (VHe) emanating from the spring per second was estimated as


Where F= average flow rate of He gas emanation (recorded); CHe= relative concentration of He in the gas mixture, which was expelled through the spring vent (recorded). The no. of moles of He gas emanated per unit second from the hot spring was calculated using the real gas equation as stated below:


Where V(i.e.,VHe) = volume of the He gas at temperature T and pressure P; R = universal gas constant = 0.0821 litre atm/ mole K; n= number of mole (to be calculated); ‘a’ and ‘b’ are real gas constants and for He, a = 0.03457 atm litre2/ mole2; b = 0.02370 litre/mole [43]. Solving the Eq. (6) and considering the real root for ‘n’, the corresponding total number of He atom (NHe) was calculated by


Here, NA= Avogadro’s number = 6.022140857×1023[44]. The relative contribution of the individual isotope in the generation of He atoms was calculated according to their relative abundance in the natural resources because the total number of He atom is produced via the radioactive decay series of 238U, 235U, and 232Th. Here, the same was not applicable to the 40 K series as it disintegrates through only β emission. Therefore, for production of He atoms,

The relative contribution of UraniumCU=conc.of UraniumUtotal conc.of UraniumU & Thorium ThE8
The relative contribution of ThoriumCTh=conc. of Thorium Thtotal conc. of Uranium U & Thorium ThE9

The basement of the study area is predominantly composed of granite gneiss belonging to the Precambrian Chotanagpur Gneissic Complex [26, 30]. Here the relative contribution of U and Th were evaluated according to their (average) content in granite type rock material, i.e., 238U [or 235U] content as 4.80 ppm and 232Th content as 21.50 ppm were considered [10, 45]. Moreover, natural Uranium is an admixture of 238U (99.28%) and 235U (0.71%) [10]. Therefore, for production of He atoms by radioactive decay,

The relative contribution of238UCU238=99.28100×CUE10
The relative contribution of235UCU235=0.71100×CUE11

The no. of the He atoms generated (in a unit second) due to the decay of radio nuclei 232Th, 238U and 235U are respectively CTh232×NHe, CU238×NHeand CU235×NHe. According to Eq. (1) to Eq. (3), it is reflected that 6 He atoms and 42.6 MeV/atom heat energy, 8 He atoms and 51.7 MeV/atom heat energy and 7 He atoms and 46.4 MeV/atom heat energy from the decay of 232Th, 238U and 235U are releasing respectively. Therefore, the energy release (per unit second) from 232Th, 238U, and 235U decay can be evaluated respectively by


And the total energy generated due to the decay of all these three radioelements were


An important issue to discuss is that the loss of generated heat energy may be considered to be negligible here as capping of the impermeable and insulating bedrock over the geothermal system prevents the heat transfer by means of conduction and convection [30, 41, 42]. Therefore, the heat energy would be stored inside the geothermal system, which may be subjected to break its dynamical stability after the accumulation of enough energy within it. However, that does not happen as excess heat is drained to the surface, along with the transfer of geothermal fluid through the spring vent [30, 42]. Moreover, here only radiogenic heat is accounted for, and the contribution of energy belonging to primordial heat sources is not included. However, [46] documented that heat from radioactive decay was contributed about half of Earth’s total heat flux, and the rest was accounted for from the primordial heat source of the Earth. Considering the similar concept, we can also assume that the primordial heat source also would contribute as much as heat energy generated by radioactive decay of radio-nuclei at the reservoir of the study area.

Therefore, the heat generated by the primordial source,


Therefore, total energy contributed from the radiogenic and primordial source is


Moreover, If the geothermal gradient dxis considered to be constant at least up to the depth of the reservoir (x) then the depth (x) of the reservoir could be calculated from the below stated linear relationship.


Where, θr= reservoir temperature of the geothermal system; θa= average ambient temperature at the study area.


4. Result and discussion

The measured parameters as well as calculated parameters such as number of He moles emanating per second (n), the total number of He atoms emanating per second (NHe), the contributed relative concentration of 232Th (CTh), 238U (CU238) and 235U (CU235) including the energy contributed due to decay of 232Th, ETh232, 238U, EU238and 235U, EU235etc. are listed in Table 2. The Table reflects that the heat energy generated per unit second by the radioactive decay of 232Th, 238U and 235U inside the reservoir are 31.58 MW, 6.3585 MW and 0.0467 MW respectively and together contributed as approximately 38 MW [radiogenic source, ERTotal]. Considering the concept of [46], as mentioned above, total heat energy [ETotal] related to the radiogenic source [ERTotal] and the primordial source [EPTotal] is expected to be 76 MW. Moreover, the values would be likely increased whenever the He emanations through the others hot springs (where He emanation is comparably less than that of Agni Kunda) and the vast surface area (soil gas) at Bakreswar would be included in this estimation. However, this was a little bit difficult as well as complicated due to the technical coerces and geographical constraints. It is notable that Kalina cycle based binary power plant using ammonia–water mixture as working fluid (thermal efficiency: 13–53%) [47], is already proposed to be installed at the spring site [48]. Accordingly, if such type of power plant is supposed to be installed, for say, the plant would be capable of delivering the power of 4.94 MW (minimum thermal efficiency 13%) to 20.13 MW (maximum thermal efficiency 53%) only considering the radiogenic heat source. These values are changed to 9.88 MW and 40.26 MW, respectively, when the primordial heat source is comprised of the radiogenic heat source.

Sl. No.ParametersParameters’ valueReferences [Eq., if any)
1He concentration, CHe1.72 vol%RA
2Flow rate, F3.5 L/minRA
3He emanation per minute, VHe0.0602 L/minEA [Eq. (4)]
4Temperature inside the spring gas trapping funnel, T342 K (69 °C)[19]
5Pressure inside the spring gas trapping funnel, P1.58 atmRA
6Number of moles emanating per second, n56.3877×106EA [Eq. (6)]
7Total number of He atoms emanating per second, NHe33.9623×1018EA [Eq. (7)]
8The relative concentration of 232Th, CTh81.7490×102EA [Eq. (9)]
9The relative contribution of 238U,CU23818.1195×102EA [Eq. (10)]
10The relative contribution of 235U, CU2350.1296×102EA [Eq. (11)]
11Energy contributed due to decay of 232Th, ETh23231.58 MWEA [Eq. (12)]
12Energy contributed due to decay of 238U, EU2386.3585 MWEA [Eq. (13)]
13Energy contributed due to decay of 235U, EU2350.0467 MWEA [Eq. (14)]
14Energy accounted for radiogenic source, ERTotal37.9834 MWEA [Eq. (15)]
15Energy accounted for primordial source, EPTotal37.9834 MWEA [Eq. (16)]
16Total Energy accounted from radiogenic & primordial source, ETotal75.9668 MWEA [Eq. (17)]
17Geothermal gradient, dx90 °C/km[33]
18Reservoir temperature, θr130 °C[40]
19Average ambient temperature, θa26 °CRA
20Depth of the geothermal reservoir, x1155 mEA [Eq. (18)]

Table 2.

Experimental and calculated parameters.

Note: EA = Estimated by the authors; RA = Recorded by the authors

Recently, the reservoir temperature [θr] of the Bakreswar geothermal system was estimated to be 126–130 °C by [40] by means of silica geothermometry. The average ambient temperature [θa] of the study area was recorded to be 26 °C. Our estimation using the mathematical relation [Eq. (18)] shows that the geothermal reservoir at the area is expected to be located at about [x=] 1,111 to 1,155 m beneath the surface. Moreover, the results of audio magnetotelluric (AMT) studies collected from [41] revealed the existence of a deep heat reservoir in the N–W part of Bakreswar. Observation sites 000, 001, 003, 005, 007, 008, 103, 105, 110, as marked in Figure 3, show low resistivity profile, and these are more favourable sites for deep drilling.

It is notable that no such work has been carried out to figure out the potential of Bakreswar geothermal region in terms of power harnessing capability. However, [42] estimated the geo-heat of the site to be 1158 KW-hr (=416.88×107Joule) by means of considering total water discharge through Agni Kunda hot spring and difference between the spring temperature and mean ambient temperature. Further, another approach may be considered here to get a comparative view in this regard. In this connection, the water discharges from the Agni Kunda and Khar Kunda hot spring were measured to be 790 L/min (∼790 kg/min) and 680 L/min (∼680 kg/min) respectively. The amount of the heat energy carried out (per second) by the hot water through these springs would be 10.29–10.70 MJ (equivalent to 10.29–10.70 MW in terms of power). However, the energy carried out by other springs is not included herewith due to technical difficulties. Furthermore, the energy carried out by the gaseous phase of the hot springs and the heat loss through the soil of the vast region is out of scope to be counted.


5. Conclusion

Using a simple technique by means of He exploration study at the field site, the probable energy generated inside the reservoir was estimated here. Considering the combined source of heat generation inside the reservoir system, the energy was expected to be generated from the source of power of 38 to 76 MW using the appropriate technology. The utilization of proper technology for power generation could facilitate to build a Kalina cycle based geothermal power plant (using ammonia–water mixture as working fluid) of power harnessing capability of 9.88 MW to 40.26 MW at the study area. Moreover, the values would be likely increased whenever, the He emanations through the others hot springs (where He emanation is comparably less than that of Agni kunda) and through the vast surface area (soil gas) at Bakreswar would be included in this estimation. However, this was a little bit difficult as well as complicated due to the technical coerces and geographical constraints. Furthermore, the deep drilling (production & injection well) of the proposed power plant to be rooted upto a depth of approximately 1,100 m at a location near to the hot spring area as indicated in the Figure 3. However, a detail geophysical survey may also be required for selecting the appropriate and exact location for drilling as well as the measurement of the horizontal (length & width-wise) and vertical (depth-wise) dimension of the geothermal reservoir at the area. The same would be subjected to accurately calculate the possible capacity of the power plant to be installed at the site.



The authors owe a debt of gratitude to the National Institute of Technology Durgapur (NIT Durgapur) and the Ministry of Human Resource Development (Presently known as Ministry of Education), Govt. of India for providing the financial as well as Institutional support in all respect for carrying out such type of the research activities at the field site, Bakreswar.


  1. 1. Barbier E. Geothermal energy technology and current status: an overview. Renewable and Sustainable Energy Reviews. 2002:6:3–65
  2. 2. Billings MP. Structural geology. New Delhi: Prentice-Hall of India Pvt. Ltd; 1990. 326 p
  3. 3. Breeze P. Power Generation Technologies. 3rd ed. Cambridge: Newnes; 2019.449 p. DOI:10.1016/b978-0-08-102631-1.00012-2
  4. 4. Mburu M. Geothermal Energy Utilisation, [Internet] 2009. Available from:[Accessed: 2021-01-24]
  5. 5. Cann JR. Metamorphism in the ocean crust. In: Talwani M, Harrison CG, Hayes D, editors. Deep drilling results in the Atlantic Ocean: ocean crust. 1st ed. Washington: American Geophysical Union; 1979. 230–238 p. DOI: 10.1029/ME002p0230
  6. 6. Armstead HCH. Geothermal energy–its past, present and future contribution to the energy needs of man. London: E & FN Spon Press; 1983. 4 p
  7. 7. Csányi Ľ, Krištof V, Kušnír S, Katin M, Marci M. Geothermal Energy: Intensive Programme Renewable Energy Sources, Železná Ruda-Špičák, University of West Bohemia, Czech Republic [Internet] 2010. Available from:[Accessed: 2021-01-23]
  8. 8. Dickson MH and Fanelli M. Geothermal background. In: Dickson MH, Fanelli M, editors. Geothermal energy: utilization and technology. 2nd ed. France: UNESCO Publishing; 2003. p. 1–27
  9. 9. Beus AA. Geochemistry of the lithosphere. 2nd ed. Moscow: MIR Publishers; 1976. 366 p
  10. 10. Pereira EB. Some problems concerning the migration and distribution of He–4 and Radon–222 in the upper sediments of the crust–a theoretical model; and the development of a quadrupole ion filter for measuring He at the soil air interface [thesis]. Houston: Rice University; 1980
  11. 11. Cupps VR. Heavy metal clocks, U–Pb and Th–Pb dating models: radioactive dating, part 7; Acts & facts 44(5) [Internet] 2015. Available from: 2021-01-23]
  12. 12. Humphreys DR, Austin SA, Baumgardner JR, Snelling AA. Recently measured Helium diffusion rate for Zircon suggests inconsistency with U-Pb age for Fenton Hill Granodiorite. Transactions of the American Geophysical Union (AGU Fall Meeting Abstracts). 2003;2003:V32C-1047
  13. 13. Humphreys DR, Austin SA, Baumgardner JR, Snelling AA. Helium Diffusion Age of 6,000 Years Supports Accelerated Nuclear Decay. Creation Research Society Quarterly. 2004:41:1–16
  14. 14. Faure G. Principles of Isotope Geology. 2nd ed. New York: John Wiley & Sons; 1986. 589 p
  15. 15. Kennedy BM, Fischer TP, Shuster DL. Heat and Helium in Geothermal Systems. In: Proceeding of Twenty–Fifth Workshop on Geothermal Reservoir Engineering (SGP–TR–165); 24–26 January 2000; California. 2000. p. 167–173
  16. 16. Joshua EO, Alabi OO. Pattern of radiogenic heat production in rock samples of southwestern Nigeria. Journal of Earth Sciences and Geotechnical Engineering. 2012:2(2):25–38
  17. 17. Nicholson K. Geothermal Fluid. 1st ed. Berlin: Springer Verlag; 1993. 263 p. DOI: 10.1007/978-3-642-77844-5
  18. 18. Hochstein MP. Classification and assessment of geothermal resources. In: Dickson MH, Fanelli M, editors. Small geothermal resources: a guide to development and utilization. 1st ed. New York: UNITAR/UNDP Centre on Small Energy Resources; 1990. p. 31–59
  19. 19. Chaudhuri H, Maji C, Seal K, Pal S, Mandal MK. Exploration of geothermal activity using time series analysis of subsurface gases data from Bakreswar hot springs area, eastern India. Arabian Journal of Geosciences. 2018:11(12):324. DOI:
  20. 20. Sathaye KJ, Smye AJ, Jordan JS, Hesse MA. Noble gases preserve history of retentive continental crust in the bravo dome natural CO2 field, New Mexico. Earth and Planetary Science Letters. 2016:443:32–40
  21. 21. Etiope G, Martinelli G. Migration of carrier and trace gases in the geosphere: an overview. Physics of the Earth and Planetary Interiors. 2002:129:185–204
  22. 22. Zhao X, Fritzel TLB, Quinodoz HAM, Bethke CM, Torgersen T. Controls on the distribution and isotopic composition of helium in deep ground-water flows. Geology. 1998:26(4):291–294
  23. 23. King CY. Gas geochemistry applied to earthquake prediction: an overview. Journal of Geophysical Research: Solid Earth.1986:91:12269–12281
  24. 24. Damon PE, Kulp JL. Inert gases and the evolution of the atmosphere. Geochimica et Cosmochimica Acta. 1958:13:280–292
  25. 25. King CY, King BS, Evans WC. Spatial radon anomalies on active faults in California. Applied Geochemistry. 1996:11:497–510. DOI: 10.1016/0883-2927(96)00003-0
  26. 26. Nagar RK, Vishwanathan G, Sagar S, Sankaranarayanan A. Geological, geophysical and geochemical investigations in Bakreswar-Tantloi thermal field, Birbhum and Santhal Parganas Districts, West Bengal and Bihar, India. In: Pitale UL, Padhi RN, editors. Geothermal Energy in India. 1st ed. Kolkata: Geological Society of India (Special Publication 45); 1996. p. 87–98
  27. 27. Wakita H, Nakamura Y, Kita I, Fujii N, Notsu K. Hydrogen release: new Indicator of fault activity. Science. 1980:210:188–190
  28. 28. Chaudhuri H, Sinha B, Chandrasekharam D. Helium from Geothermal Sources. In: Proceedings of World Geothermal Congress 2015; 19–25 April 2015; Melbourne. Bonn: International Geothermal Association; 2015. p. 1–14. Available from:[Accessed 2021-01-25]
  29. 29. Ghose D, Chowdhury DP, Sinha B. Large-scale helium escape from earth surface around Bakreswar-Tantloi geothermal area in Birbhum District, West Bengal, and Dumka District, Jharkhand, India. Current Science. 2002:82(8):993–996
  30. 30. Majumdar RK, Majumdar N, Mukherjee AL. Geoelectric investigations in Bakreswar geothermal area, West Bengal, India. Journal of Applied Geophysics. 2000:45:187–202
  31. 31. Majumdar N, Majumdar RK, Mukherjee AL, Bhattacharya SK and Jani RA. Seasonal variations in the isotopes of oxygen and hydrogen in geothermal waters from Bakreswar and Tantloi, eastern India: implications for groundwater characterization. Journal of Asian Earth Sciences. 2005:25:269–278
  32. 32. Shanker R. Thermal and crustal structure of SONATA–azoneof mid continental rifting in Indian shield. Journal of Geological Society of India. 1991:37:211–220
  33. 33. Shanker R. Heat-flow of India and discussion on its geological and economic significance. Indian Minerals. 1988:42:89–110
  34. 34. Shalivahan S, Sinharay RK, Bhattacharya BB. Electrical conductivity structure over geothermal province of Bakreswar, Eastern India. In: Proceedings of 17th IAGA WG 1.2 Workshop on Electromagnet Induction in the Earth 2004; 18–23 October 2004; Hyderabad. Hyderabad: IAGA WG; 2005. p. 18–23
  35. 35. Mukhopadhyay M, Verma RK, Ashraf MH. Gravity field and structures of the Rajmahal Hills: example of the paleo-Mesozoic continental margin in eastern India. Tectonophysics. 1986:131(3):353–367
  36. 36. Chaudhuri H, Ghose D, Bhandari RK, Sen P, Sinha B. The enigma of helium. Acta Geodaetica et Geophysica Hungarica. 2010:45:452–470. DOI: 10.1556/AGeod.45.2010.4.5
  37. 37. Chaudhuri H, Seal K, Maji C, Pal S, Mandal MK. The unrevealed facts on helium resources of India. Arabian Journal of Geosciences. 2019:12:216. DOI: 10.1007/s12517-019-4369-1
  38. 38. Majumdar N, Mukherjee AL, Majumdar RK. Mixing hydrology and chemical equilibria in Bakreswar geothermal area, Eastern India. Journal of Volcanology and Geothermal Research. 2009:183:201–212
  39. 39. Mukhopadhyay DK, Sarolkar PB. Geochemical appraisal of Bakreshwar-tantloi hot springs, West Bengal and Jharkhand, India. In: Proceedings of 37th Workshop on Geothermal Reservoir Engineering; 30 January-01 February 2012; Stanford. Stanford: Stanford Geothermal Program; 2012. p. 1010–1014
  40. 40. Gupta T, Maji C, Gupta RK, Das R, Chaudhuri H. Harnessing geothermal energy at Bakreswar geothermal area. In: Chaudhury S, Hazra AK, Balachandran S, editors. Green energy and sustainable environment. 1st ed. Kolkata: Akshar Publication; 2016. p. 57–70
  41. 41. Sinharay RK, Srivastava S, Bhattacharya BB. Audiomagnetotelluric studies to trace the hydrological system of thermal fluid flow of Bakreswar Hot Spring, Eastern India: A case history. Geophys. 2010:75(5):B187-B195 [personal communication]. DOI:
  42. 42. Majumdar RK, Majumdar N, Mukherjee AL. 2010 Geological, geochemical and geoelectric studies for hydrological characterization and assessment of Bakreswar thermal springs in hard rock areas of Birbhum District, West Bengal, India. In: Proceedings of 8th Biennial International Conference and Exposition on Petroleum Geophysics; 1–3 February 2010; Hyderabad. 2010. p. 186–193
  43. 43. Fishbane PM, Gasiorowicz S, Thornton ST. Physics for Scientists and Engineers. 3rd ed. New York: Pearson Prentice–Hall; 2005. 240 p
  44. 44. Fox RF, Hill TP. An Exact Value for Avogadro’s Number. American Scientist. 2017:95:104–107
  45. 45. Brown AA. Formation of High Helium Gases: A Guide for Explorationist. AAPG [Internet]. 2010. Available from:[Accessed 2021-01-24]
  46. 46. Gando et al. Partial radiogenic heat model for earth revealed by Geoneutrino measurements. Nature Geoscience. 2011:4:647–651
  47. 47. DiPippo R. Second Law assessment of binary plants generating power from low-temperature geothermal fluids. Geothermics. 2004:33:565–586
  48. 48. Das R, Gupta RK, Gupta T, Maji C, Chaudhuri H. Study on Geothermal Power Generation Techniques Related to Bakreswar-Tantloi Geothermal Area. In: Proceedings ofInternational Conference on Renewable Energy – Extension & Outreach; 20–21 March 2016; Santiniketan. Available from:[Accessed: 2021-01-24]


  • He emanated from other six hot springs and through the soil (soil gas) of the geothermal area are not included in the estimation due to the lack of other necessary information.

Written By

Chiranjit Maji, Hirok Chaudhuri and Saroj Khutia

Submitted: January 25th, 2021 Reviewed: February 3rd, 2021 Published: March 9th, 2021