Integrated approach using petrophysical, gravity, and magnetic data to evaluate the geothermal resources at the Rahat Volcanic Field, Saudi Arabia

It is necessary to develop and explore geothermal resources to achieve sustainable development and clean renewable energy around the Globe. Geothermal energy is crucial to the future energy supply to meet the environmentally friendly energy demand of the World. The Rahat Volcanic Field (Kingdom of Saudi Arabia) is the oldest and lengthiest Cenozoic Volcanic Field in the World. It is a dominantly a mature mafic Volcanic Field that holds three major geological events; the historic eruption (1256 C.E.), the five fingers (∼4500–1500 BP), and the seismic swarm (1999 C.E.). These incidents were studied by utilizing geological information and geophysical data sets. Geophysical and geostatistical research includes gravity and magnetic survey data, including different log curves and major elements, obtained from water samples as well as of volcanic rocks obtained by X-ray fluorescence (XRF). To gain an understanding of the subsurface thermal structure, these datasets were analyzed. The primary goal of this study is to identify the prominent potential geothermal resources with the help of an available data set. Findings suggest that beneath the historic eruption site along with the fissure eruption, on the western side, there is a geothermal anomaly with a surface footprint of about 35 km2. Analyzing gravity and magnetic data as well as density and magnetic susceptibility variations in rock samples led to the mapping of this anomaly. It has been inferred through integrated study that statistical analysis of major elements will be helpful to validate the results of the outcome.


Introduction
Over the past few years, the interest in finding productive geothermal reservoirs has increased. Geothermal energy is a good alternative to replace the hydrocarbon (HC)-based economy (Aboud et al., 2021). Although, these geothermal reservoir fields can be used to generate energy, the technologies used to make them suitable are still being developed, and significant improvements are needed in the research area of Harrat Rahat Field (Harrat is the Arabic name of Lava Fields) (Asfahani, 2018). Drilling through the high-temperature zones OPEN ACCESS EDITED BY Nima Mohamadian, Islamic Azad University, Iran is an essential step before establishing any geothermal plant (Ashraf et al., 2021). These zones gather geothermal energy, which is transported through the subsurface level to the ground to generate electricity (Chandreasekharam et al., 2014).
As a result, drilling in these reservoir zones encounters serious geotechnical problems, such as the loss of circulation (Chandreasekharam et al., 2014). It has been observed based on previous studies that loss of circulation is a frequent issue in FIGURE 1 (A) Location and geologic map of the study area (Rahat Volcanic Field) (after Aboud et al., 2022;Downs et al., 2019). (B) Location map of the drilling wells within the study area based on Google Earth Satellite imagery. geothermal drilling operations. Many geothermal wells go to be abandoned due to lost circulation issues and increased overall cost of the project. This issue may arise in the present of fractured or formations with high permeability values. Loss of circulation is also occurred due to fractured rocks or layers. Drilling at great depths or in low-pressure, depleted reservoirs increases the risk of lost circulation (Behnoud far and Hosseini, 2017; Kiran and Salehi, 2020;Saleh et al., 2020). Research indicates that greater than 10% of the typical drilling expenses are accumulated from the lost circulation difficulties (Kiran and Salehi, 2020). Snyder et al. (2019) recommended that around 90% of the circulation is lost in well drilling described for the non-beneficial period in the extended segment of geothermal drilling.
Geothermal resources are useful and contain broad reserves with advantages such as their widespread, vast reserves, eco-friendliness, reduced emissions of carbon, and extreme consumption factors (Zhou et al., 2015). Geothermal energy is categorized as a renewable energy resource because it constantly generates thermal heat within the Earth. A variety of uses can be obtained from this source, including heating systems, industrial production, and electricity generation . Although hydrocarbon resources are the main source of electricity production in Saudi Arabia with a greater percentage of electric power produced being used for cooling purposes, the potential use of envisioned geothermal energy is for producing electric power (Al-Amri et al., 2019).
It has been observed based on the literature survey that very limited research has been performed to assess Saudi Arabia's geothermal energy resources and the ongoing work mainly focused on the western part of the kingdom where hot springs and Volcanic fields are located (Abdelwahed et al., 2016). In this study, the suspected geothermal field is in the Rahat Volcanic Field (RVF), south of the city of Madinah is shown in Figure 1A. Moreover, the geothermal resources database established for Saudi Arabia demonstrates that the most likely geothermal region is found near the western parts of the Red Sea (Aboud et al., 2016). Based on the previous research, it was concluded that geothermal resources like Volcanic fields and hot springs were not encountered in any phenomena related to geothermal features (Rehman, 2005;Rehman and Shash, 2005). Previous studies of geothermal features (e.g., hot springs and volcanic fields) within Saudi Arabia concluded that volcanoes with high heat flow content eruption size present a consistent potential for geothermal power production (Hussein et al., 2013;Chandreasekharam et al., 2014;Lashin et al., 2014). Known geological formations in Saudi Arabia might accommodate geothermal systems of up to 150°C and 300°C (Downs et al., 2018).
In the classification of geological processes, the information about the temperature is considered one of the most important parameters (Norden and Forster, 2006). Earth temperature is an important factor in understanding geothermal systems. Heat production (HP), thermal conductivity (TC), thermal gradient (TG), and heat flow (HF) are all key parameters used to classify and evaluate geothermal resources. In cases where limited temperature data is available from wells, these parameters can be estimated using Earth temperature models and geophysical techniques (Förster, 2001;Ullah et al., 2022). Well-log data analysis can be used to determine the thermal conductivity (TC) of geological formations, which can vary laterally and vertically in sedimentary basins. Understanding the TC and its spatial variations can be helpful in studying geothermal resources and determining heat flow to quantify the subsurface temperature regime. TC can be determined using well logs data as a technique in the absence of measured TC values from core data (Fuchs and Förster, 2014).
Exploration geophysical techniques are mostly applied to identify the subsurface structure, as well as to detect and estimate the geothermal reservoirs at a depth of up to 200 m. Geophysical techniques play a critical role in geothermal energy exploration, helping to identify and characterize the geothermal reservoirs that are necessary for the development of geothermal energy resources (Ashraf et al., 2021). For example, magnetic and gravity techniques are used to map the variation of density and observed gravity at shallow and deep depths. Bouguer anomaly and magnetic susceptibility help to understand the geothermal reservoir and to identify the geothermal zones through the areas of interest marked on the map (Asfahani, 2018). This may reduce the number of wells planned to be drilled for classification and generate consistent interpretations (Ehsan et al., 2019;Ehsan and Gu, 2020).
The main objective of the current research is to estimate the potential of the geothermal resource in the Rahat Volcanic Field (RVF) by an integrated study that includes magnetic, gravity, major element analysis, and well logs data. The study aims to assess the geothermal potential of the Rahat Volcanic Field using gravity, well data, and magnetic techniques. These methods are important for geothermal exploration and can provide insights into the subsurface structure and potential energy resources. The finding of this study will be helpful for further exploration and development of renewable energy resources in the Kingdom of Saudi Arabia.

Petrophysical analysis
Petrophysics analysis is the study of the physical and chemical properties of rocks. In RVF, seven different well log data were available for geothermal analysis. The well GR-01, Hirma-1, Hirma-2, RAH-10, RAH-11, RAH-22, and RAH-38 were utilized for the given studies as shown in Figure 1B. The analysis was performed with help of a log curve on all seven wells. The well data of the Rahat Volcanic Field is shown in Table 1 and details of the log curve that is present in the well are shown in Table 2. The minimum variation of temperature is an encounter in all the wells because of their shallow depth. Table 3 demonstrate the statistical analysis of the well based on given log curve data.
The core samples were usually used for the calculation of heat production (HP) in the laboratory with the help of gamma ray spectrometry (He et al., 2008). Researchers used an NGR logging technique and airborne spectrometry GR method for the estimation of HP (Bücker and Rybach, 1996;Salem et al., 2005;Zhou et al., 2015). The calculation of HP at each depth is very helpful for GR (Asfahani and Abdul-Hadi, 2001). So, we used the natural gamma ray for the calculation of HP with the help of the following equation (Asfahani and Abdul-Hadi, 2001 Where HP is heat production and its unit microwatt per cubic meter, and GR is gamma ray log whereas API stands for American Petroleum Institute. Thermal energy flows from the earth's interior to the surface. A well drilled in the subsurface shows a persistent rise in the temperature with depth due to geothermal gradient (GG) (Liu et al., 2016). The temperature rise usually represents in terms of a geothermal gradient and increases per kilometer of depth . The subsurface temperature also increases eventually due to the presence of a radioactive element that increases the temperature without any influence of depth. The formula for calculating the geothermal gradient is as follows (Cooper and Cowan, 2006).
Whereas T˚formation is a formation temperature and T˚surface is the temperature of a surface. Thermal conductivity (TC) is typically estimated in the laboratory from core samples (drill cuttings) using various steady-state techniques (Si et al., 2014;Rezaei et al., 2015;Prol-Ledesma and Morán-Zenteno, 2019). Normally in worldwide, core data may not always be readily available in the targeted zone. In such cases, well logs can be used to estimate TC and obtain the well profile. To address this issue, a study by Radwan (2021) used petrophysical parameters derived from well log data to assess TC. This study used well log-derived petrophysical parameters to assess TC.

Hydrochemical characteristics of the geothermal fluids
The study involved the collection of 33 water samples on 09 March 2022 from wells located in the Rahat Volcanic Field. These samples were gathered systematically and analyzed for their major elements. Table 4 presents the Sr. No., Latitude, Longitude, Sample ID, Temperature, Water Depth, and major elements (Ca +2 , K +1 , Mg +2 , Na +1 , Si) data for each water sample. By analyzing the chemical composition of geothermal fluids and monitoring changes in their properties over time, scientists and engineers can develop strategies for the sustainable use of geothermal resources and mitigate potential environmental impacts. The samples were analyzed at the Center of Excellence in Environmental Studies, King Abdulaziz University, Kingdom of Saudi Arabia, using inductively coupled plasma atomic emission spectroscopy (ICP-AES).

Major elements analysis of rock samples
The major elements data set consists of 691 rock samples and was obtained from the USGS. The major elements in each sample were analyzed at the GeoAnalytical Laboratory in Pullman, Washington, US through X-ray fluorescence (XRF) analysis (Downs, 2019; https://www. sciencebase.gov/catalog/item/5c01bcc7e4b0815414cc723a). Analysis of geothermal sources in water from several chemicals such as various oxides and elements as iron and calcium released from several anthropogenic sources is becoming a worldwide concern (Rapant et al., 2017). Analysis of geothermal resources with the help of a water sample gives a clue about geothermal reservoirs in the subsurface (Dai and Chen, 2008). The chemical contents of water including calcium, sodium, potassium, magnesium, and silicates, are measured based on an available data set.

Gravity and magnetic data
The United States Geological Survey (USGS) acquired the geological and geophysical data in the Northern Rahat Volcanic Field within the framework of their 1 to 75,000 scale geological mapping. The area that is covered in the data set is approximately 3600 km 2 . Gravity data was acquired using 300 gravity stations with the help of a gravimeter and differential GPS (Langenheim, 2018; https://www.sciencebase.gov/catalog/item/ 5b0d8c1ce4b0c39c934b04c1). The gravity survey was carried out with the help of fixed wing light and the data is recorded in the eastto-west direction, which covers both sides of the Precambrian basement outcrop. The distance between the base camp and to stations was approximately 1-1.5 km. Different types of data correction such as latitude, free air, Bouguer, and terrain correction are applied to the provided gravity data (Langenheim et al., 2019).
Moreover, magnetic surveys were acquired  by Arabian Geophysical and Surveying (ARGAS) company supervised by USGS (Langenheim, 2018). The data sets were accessible for the total field in the public domain. These acquired data were utilized to analyze and interpret the subsurface density and susceptibility maps. The reduction and corrections are applied to the gravity data for interpretation purposes (Aboud et al., 2015;. It should be noted that changes in the density or susceptibility of rocks are related to heat or thermal change (Salem et al., 2005). At a temperature of 570°C (Curie temperature), certain magnetic materials undergo a sharp change in their magnetic properties, and above this temperature certain materials lose their permanent magnetic properties (Aboud et al., 2016).

Petrophysical interpretation
Petrophysical and statistical analysis techniques were normally used to evaluate the geothermal reservoirs. For example, the temperature of the Earth is a valuable variable in the examination and characteristics of numerous essential geothermic estimations such as HP, TC, GG, and HF, when the observed well temperature data is not reliable (Gegenhuber, 2011). The statistical analysis of each well was performed with the help of Python language.
Natural Gamma Ray (NGR) log is extremely useful for understanding the lithology and facies, which are compulsory factors that help in the recognition of fundamental parameters of geothermal resources (Ashraf et al., 2021;Aboud et al., 2022). NGR log has been utilized to calculate the HP and was significantly employed in recognizing the thickness of the geological units  (Asfahani and Abdul-Hadi, 2001;Asfahani, 2018). Thermal conductivity and geothermal gradients are important factors that can be applied to contact the flowing of a fluid region in the wells data (Haffen et al., 2013). The evaluation of tectonics can be decided by using the difference in HF, which can give the impression of frictional heating (Zheng et al., 2016). The Not a Number (NaN) is

FIGURE 2
Detailed view of well log profile of RAH-10 well. Track 1 (correlation) shows the caliper log. Track 2 (resistivity) demonstrates the variation of the resistivity log. Track 3 showed the natural gamma ray log and density log the value of the density log is merged with the natural gamma ray log. used where the data has not been recorded. Table 3 demonstrates the detailed statistical analysis of the log curve in which the name of the log curve is shown in column 1, the total number of values is counted in column 2, and the mean, standard deviation, and minimum value of all log curve are created in column 3,4 and 5 respectively. The maximum value is shown in column 9.
Temperature contrast can be produced by several factors such as time reaction, position, depth, and geology, as shown by the temperature logging curves, which exhibit an increasing alignment with the logging rates in depth. Due to differences in speed, temperature inclination, and duration, it is also possible to measure the temperature curve among many logs using the probe (Garcia-Gutierrez et al., 2002;Rezaei et al., 2015). Based on the well data available, the evaluation of the result is carried out on all the wells. The best result shows that the average temperature for well RAH-10 is 32.88°C. There is no temperature variation shown in the well. The estimated temperature does not represent the accurate temperature of the formation because of the substantial turbulences in the drillings. A low Gamma Ray (GR) log value, which typically ranges from 0 to 50 counts per second (CPS), indicates that there is less radioactive material present in the formation being measured. Gamma ray readings can vary depending on the lithology or rock type. Clay typically contains a high amount of clay minerals like illite, kaolinite, and smectite, which show a low gamma ray response. Basalt, on the other hand, is a volcanic rock composed of minerals such as plagioclase, pyroxene, and olivine, which have a high gamma ray response. Therefore, the gamma ray log values in basalt are generally higher than in clay. The gamma ray value is low in clay and higher in basalt (Rider, 1986). HP is calculated with the help of the GR log. HP logs the fluctuation in the overall log profile which demonstrates that there is a chance of a geothermal reservoir present in the subsurface. The geothermal gradient and heat flow show no variation all over the log, a constant value is observed in the log profile as shown in Figure 2. The result of GG and HF show no presence of geothermal resources in the subsurface. In the lithology section of the log curve, three lithologies are identified based on the behavior of the log curve:  (Middlemost, 1975) showing the geochemical classification of the Northern Harrat Rahat region basalts as alkalic and alkalic-sub-alkalic basalts. Source of data: (Downs, 2019). (C) Heat map of the major oxide elements shows the correlation among different major oxides. Source of data: (Downs, 2019).

Frontiers in Earth Science
frontiersin.org 08 Clay, Basalt, and Weathered Basalt. We chose to display the results of RAH-10 ( Figure 2) because this well is characterized by continuous recording and no missing data is observed during the measurement. Several parameters can indicate the potential for a geothermal reservoir, including high porosity values, variations in gamma ray values, high thermal conductivity values, and specific rock types such as granite or basalt. For example, studies have shown that high porosity values were a key factor in the development of the Longwangmiao geothermal reservoir in China, while variations in GR values were useful in identifying different lithologies in the Menengai geothermal field in Kenya. Additionally, high thermal conductivity values may suggest the presence of fluids transporting heat (Xiao et al., 2020;Nyakundi et al., 2021).
Based on the well log data of RAH-10, the depth interval between 12 m and 45 m shows potential for geothermal resources. The gamma ray and porosity logs for this particular interval indicate higher values. If the geothermal resources in this area have high heat productivity, it implies that they have the potential to produce significant amounts of thermal energy over an extended period.
Generally, geothermal gradient depends on a geological formation's thermal conductivity (the productivity with which that formation transfer heat or in the case of the earth, permits loss of heat). Table 5 represents some ranges of thermal conductivity for typical lithology (Garcia-Gutierrez et al., 2002). The wells data interpretation of RVF indicates the lithology is mainly composed of clay, basalt, weathered basalt, and silty sand. As is well known, the thermal conductivity of basalt is high, and this strengthens the possibility that we considered it as one of the main geothermal energy resources in the study area.

Hydrochemical characteristics of the geothermal fluids
The lowest and highest values for the most important ion (Ca 2+ , Na + , Mg 2+ , K + , Cl − , Si) statistical analysis results are given in Table 6. Na + and Ca 2+ cations dominate in geothermal fluids. The Na + concentration of the geothermal fluids ranged from 87 to 889 mg/L and the highest Na + concentration detected in the water is approximately 889 mg/L in water samples ES and 15. The highest Ca 2+ concentration (525 mg/L) is observed from wells in the Harrat Rahat Field. Mg 2+ with also high concentrations compared with other cations varied from 12 to 336 mg/L. The lowest concentration K + value lies between (3.9-24.5 mg/L) which was observed in the water sample. The Si-concentrations of the water sample were found to range from 2.2 to 18.3 mg/L which is relatively low compared with other elements.

Major elements analysis of rock sample
The major element (their oxides) being analyzed with XRF equipment has an average SiO 2 content of 50.7812 wt.%. The data shows a trend of increasing cumulative Na 2 O + K 2 O with increasing SiO 2 on the total alkali vs. silica (TAS) diagram, as seen in Figure 3A. The SiO 2 vs. K 2 O diagram ( Figure 3B) with fields of Na 2 O+ K 2 O, sub-alkalic, and low potassium sub-alkalic rocks (Middlemost, 1975) shows the Harrat Al-Madinah rock covering a broad variety of SiO 2 and K 2 O and are unique in being alkali to sub-alkali basalts (Moufti et al., 2012). The rocks have an average value of MgO (4.94%), CaO (7.29 wt.%), and mean Fe 2 O 3 (10.90 wt.%). The Al 2 O 3 contents mean is (16.54 wt.%) showing restricted variation, while TiO 2 average value is 1.98 wt.%, with the highest values being in the basaltic trachyandesite, Trachyandesite Trachydacite. Trachydacite has the highest contents of Na 2 O, their average value is (4.75 wt.%) and K 2 O (1.82 wt.%). As a correlation among oxide elements using various Harker-diagrams and TAS diagrams, the variation of major elements ( Figure 3C) shows good negative correlation (i.e., Total alkali silica (TAS) diagram), weak negative correlation (i.e., Al 2 O 3 ), or positive correlation (i.e., CaO). Table 7 shows the details of statistical analysis of mineral oxides present in the study area.
Heat maps are a useful tool for visualizing the distribution of a particular variable or the correlation between multiple variables. In the context of geothermal studies, a heat map could be used to display the major oxide elements present in a given area, revealing correlations among these elements and showing which combinations are typically found together or apart. This information is important for understanding the origin and evolution of rocks and minerals in the area, as well as the properties of the geothermal system itself. By highlighting patterns and relationships in the data, heat maps can reveal insights that may not be immediately apparent when looking at the raw data alone. In a heat map, darker colors indicate a higher frequency or intensity of the variable being measured, while lighter colors indicate a lower frequency. In geothermal studies, it is common to find rocks and minerals in the Earth's crust that contain different major oxide elements.
Based on XRF studies of rock in the Rahat Volcanic Field, the trend of increasing Fe 2 O 3 (average value=10.90) and TiO 2 (average value=1.98) with decreasing MgO (average value=4.94) argues against significant fractionation of Fe-Ti oxides in the early stages of Silicic volcanic rock. P 2 O 5 (average value = 0.51) also shows a good negative correlation with MgO in the silicic volcanic, consistent with apatite fractionation, which is further supported by the appearance of apatite in the silicic lavas. Moufti and Németh (2016) studies also indicates the presence of silicic volcanic.

Hydrochemical characteristics of rock samples
In the earth's crust, the 6 th abundant element is sodium and found in several minerals like rock salt (NaCl), sodalite, and feldspars. Most of the salts of Na are highly soluble in water in Harrat Volcanic Field. The concentration of sodium that exists in the water sample is in the range of 87-889 mg/L. According to the World Health Organization (WHO), 500 mg/L is considered good for human health as well as the environment (Bermudez et al., 2011). Based on the allowable limit of WHO the samples F4, F5, F6, F7, ES, and 15 exceeds the limit, which is the indication of a geothermal reservoir on the location at which the samples were taken shown in Figure 4A. Chloride is mostly found as a component of salt (Sodium Chloride) and in a few cases, it is present in combination with calcium (Ca + ) and potassium (K + ). The concentration of Calcium and Potassium is in the range (27-525 mg/L) and (3.9-24.5 mg/L), respectively, whereas the WHO limit is 250 mg/L, the samples F5, F6, F7, ES, and 15 exceed the limit as shown in Figure 4B. Magnesium concentration also exceeds in samples F4, F5, F6, F7, F8, ES, and 15 as shown in Figure 4C.

Gravity and magnetic study
The gravity studies in the Rahat region reveal shallow subsurface density variations that are linked to the volcano's structural and magmatic history. The density variations are more significant near the surface compared to deeper depths. A connection exists between high and low gravity values in the volcanic anomaly, as demonstrated in Figure 5A, which shows the elevation map of the Harrat Volcanic Field between depths of 600-1100 m. The graph in Figure 5B was created using acquired gravity data obtained from USGS and shows observed gravity variations. As seen in the graph, the gravity values vary from station to station, indicating localized variations in density or mass distribution within the volcanic anomaly. The observed variation in gravity, as depicted FIGURE 6 (A) shows an elevation contour map created using gravity data. The map reveals a variation in elevation from north to south, with higher elevations observed on the southern side. The map was created using elevation and coordinate data. The low elevation is observed on the northern side of the map. The average elevation obtained by the gravity data lies between 550 and 1100 m in depth. (B) Observed gravity map of the Rahat Volcanic Field shows the variation in gravity. The variation in gravity is increased from south to north. The zone of interest shown in ellipsoid black color usually demonstrates a higher value. (C) Bouguer anomaly map of Rahat Volcanic Field.

Frontiers in Earth Science
frontiersin.org in the Bouguer anomaly profile of the available station, is an inverted image of the elevation graph shown in Figure 5C. This means that regions with higher elevation and lower mass density can have negative Bouguer anomalies, while regions with lower elevation and higher mass density can have positive Bouguer anomalies. The reason for using these parameters profile is to understand the subsurface structure and variation along the acquired survey profile. Understanding the subsurface structural system and how it relates to geothermal resources was accomplished using gravity data. The overall area is characterized by a basin-like structural system (Aboud et al., 2015). Three different types of maps based on gravity data are prepared to understand its behavior on Rahat volcanic field; elevation, observed gravity and Bouguer anomaly map have been generated. There is a correlation between elevation and gravity, whereby an increase in elevation leads to an increase in the distance between an object (such as geothermal bodies) and the Earth's center. As a result, the force of gravity decreases slightly with higher elevation, causing the acceleration due to gravity to be slightly less than 9.81 m/s 2 (Meiburg et al., 2015). Overall, the gravity values are reduced due to an increase in elevation as shown in Figure 6A. Figure 6B shows the observed gravity map of the Rahat Volcanic Field. It displays the low gravity anomalies beneath the historical sites (1256 AD) and the fissure eruptions (641 AD) (Langenheim et al., 2019). The low gravity anomalies are bounded by large gravity anomalies involving thermal variation.
The gravity anomaly area refers to the area where the density of crustal material changes sharply along the horizontal direction, and this indicates the existence of a graben system (Langenheim et al., 2019). It is caused by an uneven distribution of underground rock mass and mineral density, or by the density difference between geological bodies surrounding rocks. The Bouguer gravity anomaly can identify or determine large fault structures arising from local or global tectonic activity, which accordingly indicates the occurrence a potential geothermal zone. The Bouguer gravity map is shown in Figure 6C.
The variation of magnetic anomalies that usually observed in the magnetic plots does not occur unsystematically; it mainly reflects a certain structural geological system. High magnetic values are often associated with rocks that have a high magnetic mineral content, such as basalts. However, the strength of the magnetic signal can vary widely and may be negligible or very weak. In a geothermal FIGURE 7 (A) Magnetic susceptibility variations are shown in samples as the range of magnetic susceptibility is from 0 to 45 SI. (B) Bulk density variation is the overall field the density range lies between 1.5 and 3.0 g/cm 3 . The highest density variation is shown in the range of 2.5-2.9 g/cm 3 , the sample from 0 to 400 mostly shows the higher values which are a clear indicator that in the subsurface some geothermal anomaly is present.

Frontiers in Earth Science
frontiersin.org 13 environment, there is generally an inverse relationship between magnetic susceptibility and high temperatures (Aboud et al., 2022). Figure 7A is displaying the generalized trend of susceptibility at each sample point. Figure 7B is displaying bulk density variation in the overall field.
Magnetic susceptibility and density of rocks play an influential role in determining physical properties of rocks samples (Bérubé et al., 2018). An increase in density can lead to a decrease in thermal conductivity due to the inverse relationship between the two properties.
While magnetic anomalies are often used to infer the age and orientation of the Earth's magnetic field, they are not a direct cause of geothermal or tectonic stress variations. These stresses are primarily driven by processes such as mantle convection, plate tectonics, and the Earth's internal heat budget. However, the distribution of magnetic anomalies can provide insights into the history of these processes and the evolution of the Earth's crust and mantle. A magnetic anomaly map can be used to create a model without requiring complicated modifications. In specific conditions, small magnetic anomalies are strictly correlated with each other. Active groundwater temperatures can weaken the magnetic properties of rocks through thermal alteration. Tectonic activity can also diminish the magnetic characteristics of rocks along the stress path and in the tectonic fracture zone, leading to weakened aero-magnetic properties. Figure 8A shows the overall variation of magnetic susceptibility, with the highest values observed on the northwest side. Figure 8B is a density map of the Rahat Volcanic Field, showing higher density values ranging from 1.5 to 2.9 g/cm 3 .

Conclusion
It has been concluded that no geothermal body is present in the subsurface surrounding all studied wells based on the analysis of the petrophysical data and the temperature log. The range of temperature in all well is between approximately 30 o to 49°C. The calculated parameter like heat productivity is estimated based on a natural gamma ray log which also supports the conclusion of no geothermal body is present in the surrounding. The geothermal reservoir is also estimated with the help of radioactive rock and minerals. Major element analysis (XRF analysis) proves that the rocks present in the subsurface are mainly basalt. Previous studies based on geological information of the study area indicates that the flow of magma in past historical eruption provide a clue for the presence of geothermal resources. Major elements of water samples that are obtained from different regions of the Rahat Volcanic Field demonstrate that some samples show very high value exceeds the limits of WHO in hydrochemical analysis which is an indication of geothermal resources in the subsurface. Samples F5, F6, F7, ES, and 15 indicate the presence of a geothermal reservoir in the subsurface. The magnetic and gravity data variations in density and susceptibility maps can be examined (highest density and susceptibility zones) which represent the geothermal anomalies. Three different types of anomalies are situated beneath the historical eruption, fissure eruption, and north of the swarm area. Such superior density is suggested for high possible geothermal resources.

Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.