A sustainable clean energy source for mitigating CO 2 emissions: numerical simulation of Hamit granitoid, Central Anatolian Massif

on fluid injection rates and wellbore length emphasizes the importance of fluid resident time for effective heat extraction. Closed-loop systems pose fewer geomechanical risks than fractured systems and can be developed through site selection, system design, and monitoring. Geothermal wellbore casing material must withstand high temperatures, corrosive environments, and should have low thermal conductivity. The HG exhibits the highest heat generation capacity among Anatolian granitoid intrusions and offers potential for sustainable energy development through EGS, thereby reducing CO 2 emissions.


Introduction
As the population and energy needs grow, so do the environmental problems of today.Türkiye's energy demand would increase by 93% and CO 2 emissions would increase by 25% if current consumption trends continue through 2040 (International Energy Agency (IEA) 2014, 2021).Countries supporting the Paris Agreement have already planned for nearly 30 years to reach their net-zero CO 2 emissions target.The netzero path assumes 90% decarbonization of the power sector in Türkiye by 2040 (World Bank Group 2022).The only way to achieve this target is to increase the share of renewables in the primary energy mix.In 2020, only ¼ of the energy generated in Türkiye came from renewable energy sources (International Energy Agency (IEA) 2019, 2020).Renewable energy is an important component of sustainable development, as it is a cost-effective, environmentally friendly, and decentralized alternative to conventional energy production Within the framework of renewable energy, geothermal energy is the only energy model that does not require climate adaptation and that we can use continuously.
Heat comes from a variety of sources, including gravitational energy released during formation, the decay of naturally occurring radioactive isotopes, and the original heat from the formation of the world.The temperature rises on an average of 3 °C every 100 m from the surface to the center of the planet.The earth has a virtually unlimited heat potential, as 99% of the earth is hotter than 1000 °C (Bächler 2003).Hydrothermal reservoirs are those that are naturally hot enough and permeable, while enhanced geothermal systems are those that are naturally hot enough but developed by hydraulic excitation.Hydrothermal reservoirs are geothermal resources that are characterized by their innate ability to store and transmit heat due to naturally occurring conditions within the Earth's crust.These reservoirs are composed of porous and permeable rock formations that contain hot water or steam, which can be readily accessed through drilling.The permeability of the rock is a crucial characteristic, as it allows for the movement of geothermal fluids that can be brought to the surface to generate electricity or for direct use in heating.The term "hydrothermal" refers to the presence of water that can carry thermal energy from depth to the surface.These reservoirs typically occur in regions with volcanic or tectonic activity, where geological processes have created fractures and pathways that facilitate the flow of geothermal fluids.That is, the first approach to construct fracture networks to create conductive fluid pathways in rocks with low permeability in deep geothermal reservoirs is enhanced geothermal systems (EGS), in which a high-pressure fluid is injected (Ellsworth 2013;Olasolo et al. 2016).Considering the sustainability of improved fracture permeability, homogeneous distribution of flow between multiple parallel fractures, and induced seismicity, closed-loop geothermal systems offer numerous advantages over their open-loop counterparts.Closed-loop geothermal systems circulate heat transfer fluids through an underground network of pipes to exchange heat with the ground.They can be designed as either vertical or horizontal loops, depending on site-specific constraints and requirements.They have higher efficiency because the heat transfer fluid is enclosed in a closed loop, minimizing heat loss.The heat transfer fluid absorbs the thermal energy of the earth and high enthalpy geothermal systems can support power generation.Thermo-hydraulic-mechanical (THM) modeling is essential for closed-loop geothermal systems because it provides a more holistic approach to system design and ensures long-term reliability and efficiency.
Türkiye has many hydrothermal provinces due to the high thermal regime, with the most important high enthalpy provinces located in western and central Anatolia (Baba et al. 2014).In addition, there are many alkaline magmatic intrusions in the Anatolian plateau that have the potential to host enhanced geothermal systems (EGS) resulting from plate melting of a series of subduction zones dipping north of the Hercynian-Hellenic arc (Ayzit et al. 2022;Chandrasekharam et al. 2022a).Granitoids in regions with a similar tectonic regime and deeper Curie point depths around the world (Soultz, Cornwall) host enhanced geothermal systems (Koelbel and Genter 2017; Letcher 2020).The Alpine Mountain Belt includes plutons of different ages and compositions that are abundant in Anatolia.Plutonic rocks with high radioactive element content, directly related to alkaline magmatism, are present along the major suture zones.Türkiye has great potential to enrich its geothermal resources with radiogenic granitoids.The EGS potential should be evaluated to exploit, in particular, alkaline magmatic intrusions associated with Page 3 of 18 35 Vol.: (0123456789) the Izmir-Ankara-Erzincan Suture Zone (IAESZ) as a future renewable energy source.In the present paper, it is assessed the power generation capacity of Hamit granitoid (HG), with an outcropping area of about 120 km 2 located within the west of the Central Anatolian Massif (CAM), defines the technology involved in harnessing the heat from these granitoids and a detailed thermo-hydraulic-mechanical (THM) model is developed for a closed-loop geothermal well with geological settings.

Geological assessment
The Anatolian plateau, located on the western Alpine-Himalayan orogenic belt, is a complex collision of the Gondwanan and Laurasian microcontinents from the Late Cretaceous to the Paleogene (Fig. 1) and is composed of several tectonic zones, which are mainly the Pontides, Anatolide Taurides, and the Central Anatolian Massif (Şengör and Yilmaz 1981; Okay and Tüysüz 1999) (Fig. 2a).
The CAM is roughly triangular-shaped and restricted by the IAESZ to the north, the left-lateral Ecemiş fault to the east, and the right-lateral Tuz Lake fault to the west (Fig. 2b).The Central Anatolian Massif has been affected by uplifting, blockfaulting, and young intracratonic over-thrusts during the closure of the Neotethyan Ocean along two northdipping subduction zones beneath the Pontide arc at the beginning of the Late Cretaceous (Kadıoglu et al. 2006;Lefebvre et al. 2013).

Petrology of the north-western margin of CAM
The Central Anatolian Massif consists of the metamorphic crystalline massifs/complexes of Kırşehir, Akdağ and Niğde (CACC).These are Paleozoic -Mesozoic, medium-to high-grade metamorphic rocks overlain by Late Cretaceous ophiolitic units and intruded by numerous plutons ranging in age from Late Cretaceous to Miocene.The magmatic activity of the Central Anatolia region is characterized by calcareous subalkaline transitional rocks and alkaline plutonic rocks that have been influenced by crustal assimilation associated with fractional crystallization processes of similar ages in a collision-related tectonic setting since the Maastrichtian (Seymen 1983;Ilbeyli 1998).
The magmatic activity is occurred from 110 to 65 Ma, according to the age of emplacement of the plutons described in the literature (Boztuğ and Harlavan 2008;Göncüoğlu 1986;Güleç 1994;Isik et al. 2008;Köksal et al. 2012;Whitney et al. 2003).The oldest granite supersuit is forming as an outer magmatic belt in the western and northern CAM.As move towards the center, monzogranite supersuits and syenite supersuits are occurred, respectively (Kadıoglu et al. 2006).
The north-western boundary of the massif is approximately a 120-km-long strip with massive outcrops Bekrekdağ, Cefalıkdağ, Çelebi, Baranadağ and Hamit that stretches from northeast of Kırıkkale in the north to Kırşehir in the south, bending from NE-SW (Table 1).
The most common mineral phases of the Central Anatolian plutons are feldspars (plagioclase and alkali feldspar).Alkali feldspar can account for about 90% of the Hamit pluton's alkali feldspar syenite.Alkali feldspar syenite, quartz syenite, and foid monzosyenite are the most common rock types in the  Hamit granitoid according to Streckesisen (Streckeisen 1976) classification and the foid-bearing rocks (nepheline monzosyenite and pseudoleucite monzosyenite) are the dominating lithology.Hamit pluton has its source the lithospheric mantle that is metasomatized by subduction fluids.The foid-bearing micro-syenitic dykes cut the nepheline monzosyenite and pseudoleucite monzosyenite in E-W and NE-SW directions.The foid-bearing microsyenitic dykes are mostly phonolite including plagioclase, clinopyroxene, alkali feldspar, and nepheline phenocrysts, and present a well-developed porphyritic texture.Quartz is a late phase that fills the gaps between the other phases and is intergrown with plagioclase as well as alkali feldspar (Ilbeyli 1998).

Stress regime over the region
Extraction of heat from such highly radiogenic granitoids requires the creation of a network of fractures by the hydro fracking process.To create a horizontal fracture network, the stress regime acting on the rock mass needs to be evaluated.The structural geology of the Central Anatolia Massif is explained by the synchronized movement of the subduction zone beneath the Pontides and an intra oceanic subduction zone in the Neotethys.
The Central Anatolian Massif ~ E-W tension forces are still active since the subduction segment in the ~ N-S direction has retreated to the west since the Cretaceous (Hinsbergen et al. 2016).Central Anatolian Massif's plutons are mostly controlled by NW, WNW, and NE direction tectonic structures however the Central Anatolian magmatic belt is thought to have originated along with an NNE trending, E-dipping subduction zone (Lefebvre et al. 2013).
Hamit pluton is located between the dextral Delice-Kozaklı Fault Zone and sinistral Savcılı Thrust Zone which are NW and WNW trending, respectively (Fig. 3).As Işık et al. (2014) have stated the Savcılı Thrust Zone is delineated as a WNW-striking and SW-dipping fault carrying mostly N10 directional slicken lines.The Delice-Kozaklı Fault Zone is the right-lateral transpressional is presenting that 200 of clockwise rotation.The northwest Central Anatolian Massif is presenting a slight counter-clockwise rotation move which is 6.10 ± 3.60 (Lefebvre et al. 2013).Plutonic rocks are overlaid by circa 1 km of Cenozoic sediments at low elevations of the basin.

Curie point depth (CPD) of the Hamit granitoid
Temperature affects the properties of materials that are paramagnetic and ferromagnetic, but not those that are diamagnetic.Since there is little thermal motion at very low temperatures, the alignment of these dipoles is greatly simplified.The effective field of the material disappears abruptly at a certain temperature, and its ferromagnetic property is completely replaced by the paramagnetic one (Nagata 1961).The Curie temperature is described by the change from Vol:. ( 1234567890) ferromagnetic to paramagnetic equilibrium, which was discovered by Pierre Curie.The Curie temperature does not correspond to the temperature at which magnetic susceptibility disappears, but rather to a second-order phase transition with an inflection point at which susceptibility is potentially unbounded (Petrovský and Geodaetica 2005).Therefore, the depth associated with this temperature is a measure of the magnitude of the magnetic sources.
The Curie temperatures of different ferromagnetic minerals vary, and the magnetic titanomagnetite (Fe 2 .4 Ti 06 O 4 ) mineral, which is particularly prevalent in the crust, is used as the foundation for calculating the Curie point depth.Generally, aeromagnetic data is used to calculate Curie point depth and at a temperature of 580 °C reaches the curie point, where it begins to acquire paramagnetic properties (Haggerty 1978;Schlinger 1985;Frost and Shive 1986).
Depending on the local geology and rock composition, the Curie-point temperature changes from one area to another.Therefore, it is typical to expect shallow Curie-point depths in areas with geothermal potential and thin crust.The Central Anatolia lithosphere's thermal structure, particularly its viscoelastic strength, can be characterized using knowledge of the Curie isotherm's depth (Quintero et al. 2019).
As a result of the Curie point depth studies carried out in Central Anatolia, it was calculated as an average of 15.4 km between 7.9 and 22.6 km.The Hamit granitoid was calculated in the range of approximately 12-14 km using aeromagnetic anomalies in the region (Ates et al. 2005).The minimum Curie point depth is observed in the basin west of the Hamit granitoid and the shallow Curie point depth direction is approximately northeast (Fig. 4).The northern portion of Central Anatolia appears to have a Curie point depth that is 2 km shallower than the southern portion (Ayzit et al. 2022;Chandrasekharam et al. 2022aChandrasekharam et al. , 2022b;;Uğur and Yahya 2011).

Radiogenic characteristics and EGS potential of the Hamit granitoid
Granitoids are the primary hosts of the key radiogenic elements U, Th, and K in the crust.Granitoids have a direct impact on the upper crust's surface heat flow and geothermal regime (Jaupart et al. 2016;Mareschal and Jaupart 2013;McLaren et al. 2006).
The quantity and distribution of radiogenic components in granitoids are used to estimate heat generation capacity and geothermal source potential (Ayzit et al. 2022;Pleitavino et al. 2021).Rybach's Eq. (Pleitavino et al. 2021) was used to calculate the heat generation rate (A).
Each equation should be shown on its own line and carry a number (in consecutive order) on the right margin, in squared brackets, as shown below.where ρ (g/cm 3 ) is the density of the rock (average density for granitoid is 2.7 g/cm 3 ), and C U , C Th and C K are the concentrations of uranium (ppm), thorium (ppm), and potassium (%), respectively.C K is expressed as a percentage of elemental potassium (or K 2 O multiplied by 0.83).The surface heat flow values were calculated using the equation (Lachenbruch 1968).
where Q is the heat flow at the surface, Q 0 is an initial value of heat flow unrelated to the specific decay of radioactive element at a certain time, D is the thickness of rock over which the distribution of radioactive element is just about homogeneous, and A is the radioactive heat production.Since the thin crustal thickness (~ 30 km) is observed in the north-western margin of the Central Anatolian Massif (Tezel et al. 2013), the background heat flow value of 40 mW/m 2 is considered around Hamit granitoid.
The Kırşehir Massif region in the west of Central Anatolia has several plutons that are extremely radiogenic due to the contents of high concentrations of uranium, thorium, and potassium minerals (Table 2).
The heat generation capacity of the Hamit granitoid can reach up to 21.48 (μW/m 3 ), which is much higher than the average production value of the continental crust (5 μW/m 3 ).Also, the heat flow value for this granitoid is approximately 254.8 mW/m 2 (D = 10 km), which is the highest heat flow value calculated in Türkiye so far. (1) Due to their high heat-generating capability and enabling the production potential of base-load power for societal usage, radiogenic high heat-generating granitoids are being considered as a future energy source.Considering the heat flow values, heat generation values, the probable volume of the  granitoids, and the curie depth temperature of this region, the Hamit pluton is a great site for initiating the EGS project.Firstly, the closed-loop geothermal system concept that can be applied to the Hamit granitoid region should be discussed with its mechanisms (Fig. 5).

The reservoir numerical simulation of closed-loop geothermal system
Several modeling techniques have been developed and applied to the analysis and design of closedloop geothermal systems (Schulz 2008;Fallah et al. 2021;Wang et al. 2021a, b;Dornadula et al. 2023).These techniques include analytical models, numerical models, and hybrid approaches.Analytical models provide an expedited understanding of system behavior through simplified mathematical equations, while numerical models are more detailed and flexible, accounting for the complexity of heat transfer and fluid flow processes in the subsurface.Hybrid approaches combine the strengths of analytical and numerical models to optimize system performance and minimize computational effort.Hydrothermal modeling plays a crucial role in the design and optimization of closed loop geothermal systems (Budiono et al. 2022).This approach focuses on simulating the heat transfer processes between the heat transfer fluid and the surrounding ground, considering factors such as thermal conductivity, heat capacity, and fluid flow rates.Recent advances in hydrothermal modeling have led to provide incorporating the effects of groundwater flow, accounting for complex geological structures, and implementing advanced numerical methods to solve governing equations.These advancements have allowed for more accurate system designs, leading to increased efficiency and reliability of closed loop geothermal systems.
Thermo-hydro-mechanical (THM) modeling is essential for closed loop geothermal systems because it provides a comprehensive understanding of the complex interactions occurring in the subsurface environment (Dornadula et al. 2023;Rinaldi et al. 2015;Rutqvist 2012;Rutqvist et al. 2008Rutqvist et al. , 2002)).In these systems, heat transfer, fluid flow, and mechanical stresses are interdependent, resulting in coupled processes that can influence the overall performance and stability of the geothermal installation.

Modeling methodology
The method consists of a numerical model for a closed-loop geothermal system in the Hamit granitoid region, focusing on the fully coupled thermohydro-mechanical (THM) processes and employing the finite element technique for discretization using COMSOL Multiphysics (Multiphysics 2020).The development of a numerical model provides a comprehensive understanding of the complex interplay between fluid flow, heat transfer, and mechanical deformation in the geothermal system.
The numerical model incorporates the Hamit granitoid region's geological characteristics, such as the subsurface rock formations, thermal gradients, and in-situ stress conditions, to ensure a realistic representation of the geothermal system's behavior.Additionally, the model accounts for the region-specific thermal and mechanical properties of the working fluid and the surrounding rock formations, which significantly influence the overall system efficiency and lifespan.This comprehensive approach enables the identification and optimization of key design parameters such as well length, and fluid injection and production rates, ultimately contributing to improved system performance and sustainability.
To understand the effect of THM processes during a closed loop geothermal well installation model at Hamit granitoid region, a three-dimensional geometry is considered which is 35 km long, 35 km wide and 5 km deep as shown in Fig. 6a.The size of the geometry is considered to include the stress redistribution impact by the geothermal wells on the nearby naturally faulted zone.Geological settings, petrophysical, petrographical and thermophysical parameters considered for this modeling are shown in Table 3.To shed light on the design of geothermal system at region, closed loop geothermal well technology is adopted.The wells are of different length and in order to enhance the heat exchange between the rock mass the flowing fluid.Figure 6b shows three components of the closed loop geothermal well.Section A in mainly used for injecting cold water and four different lengths of section B is considered to enhance the fluid resident time in the wellbore.Section C is used Page 9 of 18 35 Vol.: (0123456789) for producing the hot water.The length of section A is 2500 m whereas that of section C is 3000 m.Four different lengths of section B are considered in this study: 500 m, 1 km, 2 km and 5 km.Finite element formulation is employed to address mass conservation, energy conservation, and stress equations concurrently.Comprehensive information on governing equations is available in Mahmoodpour et al. (2022) and Singh et al. (2023).To solve these equations, COMSOL Multiphysics v5.6 (Multiphysics 2020) is utilized.Roller boundary conditions are designated to the side walls, displacement restriction is applied to the bottom boundary, and the top boundary is given free displacement.Elastic geomechanical stresses are assumed.The side boundaries allow fluid flow, while the top and bottom boundaries are  designated as no-flow for mass.A 0.254 W/m 2 constant heat flux is applied to the bottom surface, and the top surface is defined as a no-heat-flow boundary.The model's thermal gradient is 60 °C/km, with a surface temperature of 20 °C.An initial temperature field is developed based on the geothermal gradient and site's earth surface temperature.Hydrostatic pressure distribution is used to approximate the initial pore pressure.It is assumed that the geothermal system is saturated with water before injection or production commences.Initial stress values for the region are sourced from Çiftçi (2013).
Water's thermophysical properties, the working fluid for this geothermal site, are described in Singh et al. (2023).The numerical model is validated in Mahmoodpour et al. (2022), and a satisfactory match with operational data is obtained (Mahmoodpour et al. 2021).Unstructured meshing is adopted to discretize the geothermal reservoir.To counteract numerical diffusion due to mesh resolution, the wellbore's maximum mesh size is set to 1 m, significantly smaller than the fault zone aperture (10 m).Free tetrahedral meshes are used to discretize the entire domain, free triangular meshes for the fault zone, and edge elements for the wellbore geometry.The number of mesh elements varies depending on the type of geothermal system.All models are solved for a 30-year operational period.

Modeling results
Based on a global sensitivity analysis authors identified three key parameters influencing the energy extraction from a discretely fractured geothermal system: fracture aperture, rock matrix permeability and the wellbore radius (Mahmoodpour et al. 2022).However, for a closed-loop geothermal system, fracture aperture and rock matrix are not critically impacting the heat exchange.Therefore, wellbore radius or the injection flow rate is the main controlling parameter.Furthermore, wellbore length may also impact the heat extraction rate as longer wells will allow longer fluid resident time for the cold fluid.Therefore, a sensitivity analysis is performed based on the fluid injection rates and the wellbore length in this study.Temperature isosurface evolution in the vicinity of the closed loop is shown in Fig. 7a-c at three different times.Heat exchange is allowed only in sections A and B to extract heat from the rock whereas section C is fully insulated.A corresponding mean of combined thermoelastic and poroelastic stresses is shown in Fig. 8a-c.
The mean total stress in the vicinity of the wellbore is shown in Fig. 8a-c indicating that a small region is influenced by the heat transfer.A sensitivity study is performed for injection rate and wellbore length for the studied model and results are shown in Fig. 9.It is evident that the longer the well trajectory, the higher resident time will be available for water which ultimately results in higher production temperature.Another point to notice is that for the 500 m section B well trajectory, an increase in mass flux decreases the final production temperature at the end of 30 years.However, for longer well trajectories, higher production temperature is observed.This behavior could be attributed to the same reason that higher fluid resident time provides higher heat extraction capacity for the water.Considering a closed-loop geothermal system at Hamit granitoid region will limit the geomechanical impacts such as fracture propagation and unwanted seismic activities for the entire operational time.Permeability changes in fractured geothermal systems can be both beneficial and detrimental to the overall performance of the system.Increased permeability can enhance convective heat transfer and improve system efficiency, while decreased permeability can reduce heat extraction rates and limit the system's capacity.Thermal, hydraulic, and mechanical factors can all contribute to changes in permeability, necessitating a comprehensive understanding of the reservoir's properties and behavior.However, in case of a closed-loop geothermal model as shown in Fig. 6 does not show any noticeable permeability change in the system.

Environmental effects of closed-loop geothermal systems and factors affecting its efficiency
Closed-loop geothermal systems offer several advantages over traditional open-loop systems, including reduced environmental impact, increased operational stability, and the ability to function independently of naturally occurring geothermal reservoirs.This design reduces the risk of contamination, as no groundwater is directly utilized or discharged.Moreover, closed-loop systems are often more suitable for regions with limited water resources, as they require minimal water input and have a lower environmental impact.Despite their benefits, closed-loop geothermal systems also exhibit certain limitations.The upfront costs for installation can be higher than traditional HVAC (Heating, ventilation, and air conditioning) systems, as drilling and excavation are required for loop installation.Additionally, the performance of closed loop systems can be influenced by ground conditions and soil properties, requiring careful site selection and design.Lastly, closed loop systems may require periodic maintenance to ensure optimal performance and efficiency over their operational lifespan (Schulz 2008).
Fractured geothermal systems involve the extraction of heat from underground reservoirs with high permeability and natural fractures, allowing for more efficient heat transfer through convection.In these systems, water or other working fluids flow through the fractures, where they come into direct contact with hot rock formations.The fluid absorbs the heat and transports it to the surface, where it can be used for power generation or other applications (Schulz 2008;Rutqvist 2012;Rutqvist et al. 2008Rutqvist et al. , 2002;;Mahmoodpour et al. 2022;Singh et al. 2023).The convective heat transfer process in fractured systems is significantly faster than conductive heat transfers in closed-loop systems, leading to enhanced heat generation capacity.However, for a closed-loop geothermal system, water is circulated through the pipes, and heat is transferred from the surrounding ground to the fluid via conduction.This process is relatively slow, as it depends on the thermal conductivity of the ground and the temperature gradient between the ground and the hat transfer fluid (HTF).As a result, the overall heat generation capacity of closed-loop systems is limited by the rate of conductive heat transfer.
The geomechanical impacts of closed-loop geothermal systems are generally less severe than those of fractured systems.The primary concern in closedloop systems is the potential for ground deformation and subsidence due to the removal of heat from the subsurface.This can result in changes in stress distribution and soil compaction, which may affect the structural integrity of buildings and infrastructure near the geothermal installation.However, these geomechanical impacts are typically limited in scope and can be mitigated through careful site selection, system design, and monitoring.It should be noted that close to the wellbore trajectory, the total stress reaches a higher value due to higher injection rates.To minimize this stress, a smaller injection rate may be suitable given that the total economic output is justified.In contrast, fractured geothermal systems are subject to a range of geomechanical impacts due to their reliance on fluid flow through fractures and the associated changes in stress and temperature as investigated by Dornadula et al. (2023).These impacts can include induced seismicity, fracture propagation, and permeability changes, which can affect the performance and sustainability of the geothermal reservoir.Fracture propagation is another geomechanical impact associated with fractured geothermal systems.As fluid flows through the fractures, it can cause them to propagate and extend, potentially resulting in changes to the reservoir's permeability, heat transfer efficiency, and overall performance.Induced seismicity is a potential concern in fractured geothermal systems, as fluid injection and extraction can alter the stress state of the subsurface, potentially leading to the reactivation of existing faults or the creation of new ones.
The choice of casing material for a geothermal wellbore is critical for several reasons.Firstly, the casing material must be able to withstand the high temperatures and corrosive environments encountered in geothermal reservoirs.This is essential to ensure the structural integrity of the wellbore and to prevent the risk of casing failure, which could lead to catastrophic consequences, such as loss of well control or contamination of the surrounding environment.In addition to withstanding the harsh conditions, the casing material should also exhibit low thermal conductivity to minimize heat loss from the geothermal fluids as they are brought to the surface.This is particularly important in high-temperature geothermal systems, where maximizing energy extraction efficiency is a priority.Ignoring the heat transfer through the wellbore casing could result in several negative consequences.One significant issue would be a reduction in the overall efficiency of the geothermal energy extraction process.If the casing material conducts heat readily, a considerable amount of thermal energy will be lost to the surrounding rock formations, reducing the temperature of the geothermal fluids as they travel to the surface.This, in turn, lowers the efficiency of energy conversion systems, such as turbines and heat exchangers, leading to a decreased power output.Furthermore, excessive heat transfer through the casing could also cause thermal stresses, potentially compromising the structural integrity of the casing material.This could lead to premature failure of the casing and increased maintenance costs, making the geothermal energy extraction process less economically viable.These behaviors are necessary to keep in check while designing a closed-loop geothermal wellbore as demonstrated in this study.

EGS technology as a renewable energy source for mitigating CO 2 emissions
EGS technology is gaining importance due to the advancement in drilling technology and the availability of radiogenic granitoids in several countries.

Several countries have implemented this technology
Page 13 of 18 35 Vol.: (0123456789) as shown in Table 4 and generating both electricity and heat.Based on the field, a theoretical and modeling investigation conducted during the exploration stages of the Cooper Basin EGS project, it was estimated that 1 km 3 of such high radiogenic granitoids can generate about 79 × 10 6 kWh of electricity (Chandrasekharam et al. 2015;Somerville and Wyborn 1994;Cooper et al. 2010).The Hamit granitoid, being a batholith, with a surface area of 120 km 2 , extends to several kilometers below the surface.Thus 1 km thick granitoid at 3 km depth with a volume of 120 km 3 will be able to generate about 9 billion kWh of electricity.At least 2% of this energy can be extractable (Somerville and Wyborn 1994) which is equal to 189.6 million kWh which reduces CO 2 emissions by about 189.6 million kg.Since Hamit pluton is located over the region with high heat flow and with a curie point depth of ~ 12 km (Uğur et al. 2014;Aydın et al. 2005), the leveled cost of electricity generated from such granitoid will be comparable (8 euro cents/kWh) with the cost of electricity generated from fossil fuels based energy sources (7 euro cents/kWh) (Baba and Chandrasekharam 2022).These costs are based on drilling to a depth of about 5 km.But since the depth of drilling is sallow in the case of the western Anatolian region (Baba and Chandrasekharam 2022) and Hamit granitoid, in particular, the cost of electricity will be much lower than 7 euro cents/kWh.If a closed-loop heat extraction method is adopted instead of creating a network of fractures in the granitoids at 3 km depth, as is followed traditionally in the earlier EGS projects like that followed in Soultz, France, the unit cost of power can further be reduced (Chandrasekharam et al. 2022a).

Discussion
The scenarios produced as a result of the research show that the world's population, energy demand and therefore CO 2 emissions will increase rapidly day by day.Fossil fuel-based power plants are the main source of CO 2 emissions in Türkiye, which can be reduced by developing renewable energy technologies.Hydrothermal source can save about 32,211 million kg (~ 32 million tonnes) of CO 2 by replacing an equivalent amount of fossil fuels (about 400 million tonnes) as the main energy source (Ayzit et al. 2022;Chandrasekharam et al. 2022a).
On the other hand, since 1965, Europe has followed the strategy of consuming natural gas and renewables instead of coal, even though natural gas resources are limited.Although Europe has reduced its energy supplies from Russia since 2020, Russia currently supplies 25% of its oil needs, 46% of its coal needs, and 40% of its gas needs.Especially since February 2022, the serious tensions with Russia, with which it was not possible to reach an agreement on political issues, had a negative impact on energy supplies, and this situation had a negative effect on European energy prices.This situation showed that countries need to invest more in renewable energy to become more independent and pave the way for technological research in renewable energy.Considering climate adaptation, 24/7 applicability, small land footprint (Table 5), ability to provide baseload power, and highest efficiency in power generation, it is obvious that geothermal energy technologies are the fastest and most sustainable solution.
EGS has the opportunity to provide renewable energy to the Earth because it is an adaptable, endless, and local source.The EGS operating process can provide CO 2 sequestration, uses little water, emits less CO 2 , and provides reliable baseload power as well as a "brine mining" environment for extracting critical raw materials (e.g., Li, Sr, B, REE, etc.) from the waste stream.Deep-well drilling technology and material improvements could reduce the cost of EGS technology, and geothermal energy could eventually contribute significantly to a net-zero energy system (Fig. 10).EGS energy costs (Table 5), currently around 6-9 cents/kWh, could very likely be much lower given the thermal conditions in western and central Anatolia.It was announced that the highest geothermal base temperature in Türkiye was reached Unfortunately, it is inevitable that the problem of water security will be the main problem of the world in the coming years.Seawater can be desalinated and the problem of water security is solved by using renewable energy sources with low carbon emissions.However, after the enhanced geothermal systems have been improved with technology, wastewater and desalination may be integrated into the same system and a new generation of renewable energy technology can be developed more economically and cleaner for sustainable development.Enhanced geothermal systems may satisfy the need for water-energy-foodenvironment security as long as the earth's core continues to cool.

Conclusions
The west of Central Anatolia hosts many plutons that are highly radiogenic due to the high concentration of uranium, thorium, and potassium minerals.These high heat-generating granitoids spread over an area of almost 650 km 2 between Kırşehir and Kırıkkale, in the Kırşehir Massif region.Hamit granitoid has  the highest heat generation potential in the region.
The heat generation capacity of the Hamit granitoid can reach up to 21.48 (μW/m 3 ), which is much higher than the average production value of the continental crust (5 μW/m 3 ).Also, the heat flow value for this granitoid is 254.8 mW/m 2 , which is the highest heat flow value calculated in Türkiye so far.Practically it may not be possible to extract the entire heat in place in the Hamit granitoid.Based on the earlier EGS projects, it is estimated that 2% of this heat can be extracted for use.This amounts to 189.6 million kWh.This will help the country to reduce 189.6 million kg of CO 2 emissions and comply with the Paris 2015 agreement on reducing CO 2 emissions.This energy can be utilized for generating electricity for direct applications like space cooling and heating or greenhouse cultivation.To investigate one such idea, a fully coupled thermo-hydro-mechanical model of a closed loop geothermal model is developed in this study for the Hamit granitoid.In this study, the influence of key parameters on energy extraction from closedloop geothermal systems was investigated.Results showed that wellbore radius, fluid injection rate, and wellbore length are the primary controlling factors, with longer well trajectories resulting in higher production temperatures.The geomechanical impacts of closed-loop systems are generally less severe than those of fractured systems, with ground deformation and subsidence being the primary concerns.However, these impacts can be mitigated through careful site selection, system design, and monitoring.In comparison, fractured systems are subject to induced seismicity, fracture propagation, and permeability changes, which can affect the performance and sustainability of the reservoir.Casing material selection is critical for geothermal wellbores, as it must withstand high temperatures, corrosive environments, and exhibit low thermal conductivity to minimize heat loss.Ignoring heat transfer through the wellbore casing can lead to reduced energy extraction efficiency, lower power output, and increased maintenance costs.Careful consideration of these factors is essential for the successful design and operation of closedloop geothermal systems.This study indicates that a closed loop geothermal system will be a suitable idea to harness geothermal energy from Hamit granitoid and similar geological regions with high geothermal gradient.Thus, a future sustainable development strategy can be adopted by utilizing this green renewable energy by the country.

Fig. 3
Fig. 3 Lithostratigraphic column section of CAM (not to scale); (modified from Lefebvre et al. 2011) and the north-east direction cross-section between Tuzgölü Basin and IAESZ which is a boundary between Pontides and CAM

Fig. 5
Fig. 5 Closed-loop system concept that can be applied for the Hamit granitoid region

Fig. 6 a
Fig. 6 a Three-dimensional geometry of the Hamit granitoid region used for numerical modeling.Three major fault zones are naturally present in this region.A closed loop well design is considered in this study.b A zoomed version of the geometry indicating the well alignment in the modeled geothermal

Fig. 7
Fig. 7 Temperature development in the geothermal reservoir after a 1 year, b 5 years and c 30 years.Since the production well is fully insulated, there is no heat losses along the production well

Fig. 8
Fig. 8 Total mean stress (mean of thermoelastic and poroelastic stress) development in the geothermal reservoir after a 1 year, b 5 years and c 30 years

Fig. 10 a
Fig. 10 a Schematic of a conceptual two-well enhanced geothermal system (EGS) in hot dry rock in a low-permeability crystalline basement formation and b Drilling in hard dry rock that lacks natural cracks and connecting between the wells via closeloop technology