Techno-Economic Assessment of Geothermal Resources in the Variscan Basement of the Northern Upper Rhine Graben

Deep geothermal energy represents an essential component of the future energy supply because the resources greatly exceed the demand, and the base load capability can compensate for temporal fluctuations in wind and solar power. By far, the largest amount of heat is contained in the crystalline basement, accessible almost everywhere through deep drilling. An interdisciplinary approach for a techno-economic resource assessment was applied to provide stakeholders with a more reliable basis for decision-making in the Northern Upper Rhine Graben. This approach incorporated data from various sources such as boreholes, outcrops, geophysical surveys, geomechanical models, and operating geothermal power plants. Emphasis was placed on resources tied to large-scale fault zones, exhibiting preferential hydraulic properties. Uncertainties in the calculations were quantified using Monte Carlo simulations. The resource base in the basement of the Northern Upper Rhine Graben is about 830 PWhth, of which about 8.2–16.1 PWhth could potentially be extracted with current technologies in enhanced geothermal systems. A comprehensive sensitivity analysis was carried out, examining the economic influence of nine parameters. Considering the reference power plant and the energy prices before the Ukraine war began in February 2022, the geothermal reserves are about 65% of the resources. However, given the massive recent increase in heat and electricity prices, a higher percentage is also possible. Finally, a socio-economic-environmental assessment shows that in the Northern Upper Rhine Graben, geothermal resources largely coincide with favorable conditions at the surface, making the region a preferred target for geothermal utilization.


INTRODUCTION
A rapid reduction in greenhouse gas emissions is essential to mitigate the global effects of manmade climate change. To achieve this aim, major parts of the energy supply need to be shifted from fossil fuels to renewable sources in the coming years (IPCC, 2022). Up to now, German governments and investors have focused primarily on wind and solar power. However, the significant discrepancy between energy supply and demand throughout the year makes the construction of large energy storages necessary (Bä r et al., 2015;Bussar et al., 2016;Welsch et al., 2018). In comparison, geothermal energy enables a base-load capable energy supply such that peak loads can be compensated effectively (Tester et al., 2006;Huenges et al., 2013;Stober & Bucher, 2021). In particular for heat supply, where the dependency on fossil energy sources is still very high both in the residential and commercial sectors, there is a considerable need to expand renewables (BMWi, 2021), much of which can be covered by deep geothermal power plants (Moeck, 2022). At present, geothermal energy accounts for only a small share of the energy mix in Germany, but the expansion of this technology is essential for a successfull energy transition (Bracke & Huenges, 2022).
The Upper Rhine Graben (URG) is an ideal region for geothermal utilization due to the significantly increased temperature gradient and several potential reservoir horizons (Pribnow & Schellschmidt, 2000;Sass et al., 2011;Stober & Bucher, 2015;Freymark et al., 2017). Here, enhanced/engineered geothermal systems (EGS) are generally developed where the permeability of the fractured reservoir is improved by hydro-mechanical and/or hydrochemical stimulation. The reservoir volume is constrained by the hydraulically active fracture network during operation and can be mapped using, e.g., induced seismicity. In total, 17 deep geothermal projects have been implemented in the URG since the late 1970s, seven of which are currently in operation (Bruchsal, Insheim, Landau, Riehen, Rittershoffen, Soultz-sous-Forê ts, Weinheim), providing 50 MW th of heat and 10 MW el of electricity (Frey et al., 2022a). Thus, only a small portion of the total geothermal resources are currently being utilized (e.g., Kock & Kaltschmitt, 2012).
While aquifers are the preferred targets of geothermal drilling, such hydrothermal reservoirs hold only about 1% of the total resources (Paschen et al., 2003). In contrast, by far the largest amount of heat is contained in the Variscan basement, which is characterized by very high temperatures and an increased density of fractures and faults in the URG (Genter & Traineau, 1996;Sausse & Genter, 2005;Dezayes et al., 2010;Glaas et al., 2021;Bossennec et al., 2022;Frey et al., 2022b). The fault and fracture network may have reasonable natural permeability or has the potential to be enhanced by stimulation. However, exploration has been subject to significant uncertainties due to the small number of wells that penetrate the entire sedimentary cover and reach the basement. Moreover, the depth resolution or overall quality of most 2D seismic lines is insufficient to image basement structures accurately.
In the past 2 decades, several studies have been conducted to quantify geothermal resources in the URG basement, but their strong geologic simplification of the basement usually allows only for regional considerations (Paschen et al., 2003;Sass et al., 2011;Kock & Kaltschmitt, 2012;GeORG Projektteam, 2013;Jain et al., 2015). In the framework of the Hessen 3D 2.0 (Bä r et al., 2021b) and DGE rollout projects (www.nweurope.eu/dge-roll out), gravity and magnetic data were additionally integrated to develop the most detailed 3D model of the crystalline crust in the Northern URG to date . These results are combined with comprehensive thermal, petrophysical, geomechanical, and financial information to obtain a refined techno-economic resource assessment. This model provides stakeholders with a valuable basis for decision-making during site and reservoir selection for geothermal drilling. The deep geothermal resources of the sedimentary reservoir horizons in the Northern URG will be covered in a subsequent publication by van der Vaart et al. (in prep.).

GEOLOGICAL FRAMEWORK
The URG is a passive continental rift system that developed from the Eocene onwards due to the changing lithospheric stress field in the alpine foreland (e.g., Ziegler et al., 1995). Both location and orientation of the URG are governed by the reactivation of large fault and shear zones in the basement that formed since the Paleozoic (Schumacher, 2002;Edel et al., 2007;Grimmer et al., 2017). The complex multiphase Cenozoic rift evolution caused a distinct asymmetry both parallel and perpendicular to the strike direction, evident in the significant variations in sediment thickness throughout the URG (e.g., Doebl & Olbrecht, 1974). In the first main rifting phase, crustal extension was dominant, dividing the crust into numerous isolated fault blocks, bounded by large-scale normal faults. During the Miocene, a counterclockwise rotation of r Hmax led to the reactivation of these faults in a strike-slip or oblique sense (Buchner, 1981;Behrmann et al., 2003;Dè zes et al., 2004). The URG represents a convection-dominated geothermal play system, where local thermal anomalies mainly result from the upwelling of deep groundwaters along fault zones with increased fracture permeability in the fault damage zones (Moeck, 2014).
The basement architecture in the URG region was mainly established during the Variscan orogeny (Kossmat, 1927;Behr et al., 1984), involving the juxtaposition of Laurussia, Gondwana and several micro-terranes (Armorican Terrane Assemblage) with the main collisional phase during the Carboniferous (e.g., Kroner et al., 2008). In the Northern URG, the Mid-German Crystalline High (MGCH) makes up most of the basement (Fig. 1), which is interpreted as a former active continental margin in the north of Armorican Terrane Assem-blage (Franke, 2000;Zeh & Gerdes, 2010). The crystalline Odenwald at the northeastern margin of the URG is the largest outcrop of the MGCH and reveals substantial lateral heterogeneity in composition and structure (Krohe & Willner, 1995;Stein, 2001). In the western (Bergsträ ßer) Odenwald, subduction-related igneous rocks predominate that intruded into meta-volcano and meta-sedimentary host rocks (Krohe, 1991;Altherr et al., 1999). From north (Frankenstein Complex) to south (Southern Plutons) a transition from predominantly mafic to felsic plutonic rocks can be observed. Similarly, the age of the intrusives decreases from about 360 Ma to about 325-330 Ma toward the south (Kreuzer & Harre, 1975;Kirsch et al., 1988; Anthes & Reis-  Frey et al. (2021). Dotted signature indicates areas of heterogeneous composition and uncertain interpretation. A detailed image of the fault network in the Northern URG is given in Figure 2. BO Bö llsteiner Odenwald, FC Frankenstein complex, FG Flasergranitoid zone, SG Flasergranitoid zone. chmann, 2001). Only a few basement outcrops exist in isolated quarries along the western rift margin, and less than 10 deep wells in the URG reached the basement so far (Fig. 1). However, additional information is provided by gravity and magnetic data, which were used to trace distinct lithologic units hidden beneath the sedimentary cover (for details see Frey et al., 2021).
The MGCH is separated from the Northern Phyllite Zone (NPZ) to the north by a major continental thrust zone (Blundell et al., 1992;Brun et al., 1992), whose location is not precisely known yet. The NPZ, exposed in the southern Taunus and Hunsrü ck, comprised a tectonic mé lange of greenschist-facies meta-sediments and meta-volcanics (Klü gel, 1997). To the south, the MGCH merges into the Saxothuringian Zone, which consists of a metamorphic early Paleozoic shelf sequence overlying Neoproterozoic gneiss (Falk et al., 1995).

3D Model and Parametrization
The resource assessment approach of this study is based on a plausible 3D model of the crystalline basement in the Northern URG. The software Petrel was used to merge existing models, mainly from the GeORG, Hesse 3D and DGE Rollout projects (Sass et al., 2011;GeORG Projektteam, 2013;Bä r et al., 2021b;Frey et al., 2021;. Furthermore, additional borehole data and seismic profiles, e.g., the reprocessed DEKORP 9 N line (Homuth et al., 2021), were included. The final model spans over 110 km from north to south, 60 km from east to west, and extends to a depth of 6 km. In previous geothermal potential studies (Paschen et al. 2003;Tester et al, 2006), 10 km was often assumed as the maximum depth of technically feasible utilization, which leads to far too optimistic estimates. In fact, almost no geological data are available below 6 km, preventing any reliable statements about the hydraulic, thermal and geomechanical conditions in this depth range. The 3D model was converted into an irregular grid with a cell size of typically 500 9 500 9 500 m 3 . At the top of the basement, the cell height varies to handle the anticipated subsurface topography (Fig. 2a) as accurately as possible.
Supplementary to the structural model of the basement, information on the subsurface temperature distribution is essential. Therefore, the Germany-wide temperature model of Agemar et al. (2012) was used (Fig. 2b), which covers both the sedimentary horizons and the basement with a spatial resolution of 2 9 2 km. The model used data from approximately 10,500 wells of varying quality. A 3D universal kriging approach was applied to calculate the temperature model, providing indications of the measurement and kriging uncertainties. In the basement, where almost no temperature information is available in the model (mostly below 4-5 km depth), a constant temperature increase of 30°C/km was assumed to extrapolate temperature values. This value agrees with observed geothermal gradients below the convection-dominated depth interval, e.g., in Soultz-sous-Forê ts . It should be noted that local convection cells, lithological changes, or fault zones are poorly resolved in this purely interpolative model. However, the measured regional-scale thermal anomalies are well represented. Temperature models with a higher resolution exist for specific parts of the URG (Bä r et al., 2021b), but because they do not cover the entire model domain, they were not considered further here. The reference temperature was set to the mean surface temperature issued by the German Weather Service for 1991-2020 (DWD, 2021) (about 11°C in the Northern URG). Additionally, information about the inverted rock density was taken from Frey et al. (2021) (Fig. 2c). Petrophysical data, mainly the specific heat capacity, were adapted from the databases of Bä r et al. (2020) and Weinert et al. (2020). Information on the 3D stress tensor in the Northern URG, provided by Ahlers et al. (2021), was used to determine the reactivation potential of fault zones (see details below).

Technical Assessment
The volumetric (Ôheat in placeÕ) approach, proposed by Muffler and Cataldi (1978), was adopted in this study and further refined to quantify the geothermal resources of the basement. This rather simplistic method involves a limited number of parameters, allowing the uncertainties to be reasonably quantified. In comparison, assessment schemes based on numerical simulation require several assumptions about the reservoir (properties, pore pressure, temperature, etc.), the well comple-tion design, the pumping rate, and the operational parameters of the power plant (e.g., Jain et al., 2015). However, especially the hydraulic properties can vary over several orders of magnitude in the basement (Manning & Ingebritsen, 1999;Stober & Bucher, 2007;Achtziger-Zupančič et al., 2017). This makes any simulation result hardly reliable if either the fracture network properties or the local hydraulic parameters are not known from borehole geophysics or hydraulic well testing.
The resource base Q total (J or Wh), also referred to as heat in place or theoretical potential, describes the total amount of thermal energy stored in the Earth's crust (in this case, down to 6 km depth) relative to the ambient surface temperature T 0 (°C). Q total of the individual model cells is calculated as: where V total is the volume (m 3 ) of each model cell, q is the rock density (kg/m 3 ), c is the specific heat capacity (J/kg/K) and T r is the reservoir temperature (°C). Due to the fundamentally low porosity of the crystalline basement, the thermal properties of the pore fluid are not regarded. The geothermal resources Q rec (J or Wh), also known as recoverable heat or technical potential (Muffler and Cataldi, 1978), are defined as the share of the total thermal energy stored in the reservoir that can be extracted considering the current technologies, thus: Consequently, the thermal recovery factor R (unitless) is one of the key factors that influence the potential productivity of geothermal projects. It depends on various parameters, such as the geometry of the reservoir, permeability, reinjection and abandonment temperature (temperature drop in the reservoir at which operation becomes uneconomical), flow rate, and well design (Garg & Combs, 2015;Agemar et al., 2018). The most reliable infor- mation on the recovery factor is derived from production data of existing geothermal projects. Accordingly, the recovery factor can be as high as 25% for hydrothermal systems (Tester et al., 2006). In fractured reservoirs, R is generally smaller due to the lower permeability. In active geothermal fields in the USA, such as The Geysers, Coso, and Dixie Valley, the recovery factor ranges from about 5-20% (Williams, 2007). In contrast, Grant (2016) reported recovery factors of only 0.2-2% for the research EGS sites Fenton Hill (USA), Rosemanowes (United Kingdom), Hijori (Japan), and Cooper Basin (Australia). Recovery factors from EGS plants in the URG are currently not publicly available. In the present study, a range of R from 0 to 2% is therefore assumed for the base case, i.e., without the influence of large-scale fault zones. The possibility of zero recovery is included in the assessment to account for the risk of too low productivity.
Fault zones in the crystalline basement are generally considered to be effective fluid conduits (Caine et al., 1996;Evans et al., 1997;Bense et al., 2013). A two to three order of magnitude increase in permeability compared to the host rock is primarily observed in the damage zone, which can be several hundred meters wide depending on the displacement. Consequently, an improved heat recovery is expected when targeting major fault zones (e.g., Jolie et al., 2015). Here, a maximum recovery factor of up to 20% was assumed, which is considered the upper limit for fracture-dominated systems (Tester et al., 2006;Williams et al., 2008). Furthermore, the relative location of the fault planes within the regional stress field was taken as a proxy for maximum heat recovery (Agemar et al., 2017). For this purpose, the slip and dilation tendency were calculated for each modeled fault (Morris et al., 1996;Ferrill & Morris, 2003), which were then used as linear weightings for the recovery factor. The slip tendency T s (Fig. 3a) describes the ratio of the shear stress s to the normal stress r n on a surface (stress components are given in N/m 2 ), thus: The dilations tendency T d (Fig. 3b) is calculated from the normal stress as well as the largest and smallest principal stresses r 1 and r 3 , thus: In addition, e.g., Stober and Bucher (2007) showed that the permeability of the basement is strongly dependent on the lithology. They reported results from hydraulic tests in the Black Forest showing that permeability in gneiss can be up to two orders of magnitude lower than in granite under comparable stress conditions. The geomechanical behavior of the foliated metamorphic rocks during reservoir stimulation and production is more complex to predict than for igneous rocks. There is an increased tendency of re-closing fractures, which are formed primarily along the foliation. To account for this, the recovery factor of the metamorphic basement was downweighed by a factor of 0.1 based on the inversion results and interpretations of Frey et al. (2021).
The potential amount of electricity E (J or Wh) that could be generated in the model area is then obtained by multiplying Q rec with the unitless electrical system efficiency g (Beardsmore et al., 2010). A distinction needs to be made between gross and net system efficiency to consider the parasitic electricity consumption of the pumps, thus: According to Eyerer et al. (2020), g gross and g net can be described as a function of the fluid temperature T fl (°C), thus: g gross ¼ 13:59 Á ln T fl ð Þ À 62:38 ð6Þ The mean electrical power output P el (W) of a power plant is calculated by dividing the total gross electricity produced by the project lifetime T life (years) and the annual full load hours h full , thus: Monte Carlo Simulation The resource assessment described above involves several parameters subject to varying degrees of uncertainty. Probability density functions (PDF) can be derived directly from the data basis for some parameters, such as temperature, density, or specific heat capacity. For the recovery factor, ranges of possible values have to be assumed since only a limited database is available (see assumed ranges above). The associated data publication of this study (Frey et al., 2022d) contains the parameterized model used for the calculations, including the standard deviations or ranges of the parameters.
To quantify the risks of geothermal exploration, it is important to combine the individual uncertainties of the input data into a cumulative uncertainty of the target variable. For this purpose, a Monte Carlo simulation was implemented, which is widely used in geosciences (Sambridge & Mosegaard, 2002;Garg & Combs, 2010;Caers, 2011). The basic principle is that random samples of the input parameters are drawn from the previously defined PDFs. The generated values are then used to perform the above-described calculations. This process is re-peated 100,000 times to obtain a significant distribution, which can then be analyzed statistically.

Economical Assessment
The share of geothermal resources that can be exploited economically at present is referred to as geothermal reserves. These can be determined by comparing the levelized costs of energy (LCOE) with the average revenues earned from the sale of electricity and heat. LCOE describes the ratio of the power plantÕs net present value (NPV), including the total investment c invest , the operation and maintenance costs c O&M , and the potential revenues r, to the NPV of the net electricity generated over the project lifetime (Short et al., 1995), thus: where t is the time step (years) and d the annual effective discount rate (%). Due to the comparatively high risk of geothermal projects, 8.25% was assumed (e.g., Eyerer et al., 2020).
Compilations of the investment costs of geothermal power plants are provided by various sources (Stefá nsson, 2002;Sanyal, 2004;Beckers et al., 2014;Limberger et al., 2014;Schlagermann, 2014;Adams et al., 2021). The cost components described therein are in some cases site-and power plant-specific and therefore not directly applicable for each region. To compare data from different years, the Producer Price Index was used to convert prices to be valid for January 2022 (Destatis, 2022b). After consulting with German industry experts (e.g., Menzel 2022), a generalized cost overview was compiled (Table 1). Accordingly, the total invest-ments consist of the costs for the exploration, the drilling including stimulation, hydraulic tests and geophysical logs, the submersible pump, the power plant at the surface, the management, engineering, planning, and public relations work, the seismological monitoring network, and the insurance of exploration risks. In the following, the major cost elements are discussed in more detail.
Drilling and well completion are the largest cost figures and can reach up to 70% of the total investment for EGS projects. Empirical well cost equations, e.g., published by Legarth (2003), Bloomfield and Laney (2005), Tester et al. (2006), Lukawski et al. (2014) and Schlagermann (2014), generally show an exponential or quadratic increase in cost with depth. For this study, the empirical equation of Schlagermann (2014) was used (Table 1), which focuses on the German market. A simple doublet system with one production well and one injection well was assumed, as this is until now still the standard for EGS projects in the URG. Geothermal wells are usually drilled from one plat- form and thus need to be deviated to develop a sufficiently large reservoir volume. The resulting measured depth z md (m) of a well was approximated as: where z tvd is the true vertical depth and s the horizontal distance of the wells (both given in m).
Apart from borehole construction, significant costs are related to surface installations of the power plant. Because the reservoir temperatures are generally below 200°C in the URG, only binary plants based on the Organic Rankine Cycle (ORC) or Kalina Cycle are applicable. Literature data on the specific costs of binary power plants display a wide range, from about 1000 e/kW to over 5000 e/kW (Rettig et al., 2011;Quoilin et al., 2013;Franco & Vaccaro, 2014;Heberle & Brü ggemann, 2015;Lemmens, 2015;Tartiè re & Astolfi, 2017). However, high values are mainly associated with small ORCs (< < 1000 kW), and costs are expected to decrease with increasing electrical power. A step function was defined accordingly (Table 1). In analogy to the power plant in Landau, an electricity-driven combined heat and power (CHP) process with a power to heat ratio of 0.5 was assumed, in which the residual heat from the ORC can be fed into a district heating network. Such a system additionally requires a heat exchanger and piping. These investments can vary greatly depending on the plant specifications. Thus, additional costs of 25% were estimated for simplification.
While capital expenditures are considered overnight expenses, operation and maintenance costs are recurring annually. They include the costs of personnel, well and surface facility maintenance, seismic monitoring, and insurance. These amount to about 3% of the total investment costs for EGS projects. Potential revenues arise from the sale of heat from CHP operations.

Sensitivity Analysis
The LCOE and thus the geothermal reserves are influenced by various technical and economic factors. To quantify the respective effect, a comprehensive sensitivity study was conducted. For this purpose, nine parameters were selected and were assumed to have significant impact on the LCOE (Table 2). Based on the experience of existing geothermal projects, values for a reference power plant were defined. These parameters were then varied systematically in equal steps (À 40%, À 20%, + 20%, + 40%). Excluded from this are the well depth and the lateral reservoir extent, as the variation of these parameters is predetermined by the fixed model discretization. Figure 4a illustrates the total geothermal resource base in the crystalline basement of the Northern URG to a depth of 6 km. The combined heat in place is approximately 830 PWh, with significant spatial variations. While the resource base exceeds 100 TWh per 500 9 500 m 2 raster cell at the western and northeastern rift margin, it only ranges from 20 to 40 TWh south of Heidelberg. These variations can be explained primarily by differences in the depth to the basement. Accordingly, the thickness of the basement wedge ranges from less than 1 km to more than 5 km within the asymmetric URG (cf. Fig. 2a). The calculation uncertainties are shown as relative standard deviations in Figure 4b. With more than 14%, they are largest at the southeastern margin of the URG. Here, the greater depth of the basement leads to higher standard deviations of the input parameters and thus also of the calculation result.

Geothermal Resource Base and Resources
To avoid the dominating effect of the varying basement thickness down to 6 km depth, three distinct depth slices (3-4 km, 4-5 km, and 5-6 km) of the geothermal resource base were extracted from the model (Fig. 5). The heat in place shown therein is primarily dependent on the distribution of temperature and petrophysical properties. In general, the heat content increases with depth as the temperature also increases. Between 3 and 4 km depth, the resource base is very heterogeneously distributed, mainly because the basement is not yet reached or only partially reached at this depth interval in the eastern URG (see Fig. 2). Between 4-5 km and 5-6 km depth, the resource base varies by about ± 10% from the average. Consistently high heat in place values are present in the area southwest of Mannheim due to the local positive temperature anomaly.  Figure 6 shows three alternative scenarios of potentially recoverable geothermal resources considering the currently available EGS technologies. Scenario A (Fig. 6a) describes the base case, where a maximum R of 2% is assumed for the entire basement, as observed in EGS projects worldwide. The combined resources, in this case, are about 8.2 PWh th , allowing a net electricity generation of 0.7 PWh el . The highest values of up to 2 TWh th per 500 9 500 m 2 are again located at the western and northeastern edges of the Northern URG. Scenario B (Fig. 6b) integrates the influence of large-scale fault zones and the regional stress field into the calculation. As a result, the resources amount to about 16.1 PWh th , with a potential electricity production of 1.3 PWh el . The highest values of up to 5 TWh per 500 9 500 m 2 are reached along favorably oriented faults and at fault intersections. The area along the southwestern rift margin exhibits particularly high resources due to the increased fault density and the positive local temperature anomaly.
In addition, there are significant geothermal resources around Mannheim and west of Darmstadt. Scenario C (Fig. 6c) also considered the basement lithology. It was assumed that less favorable hydrogeological conditions characterize the predominantly metamorphic and mafic basement compared to the granitoid basement (Stober & Bucher, 2007). Accordingly, the northwestern and southwestern rift margins are less suitable for deep geothermal exploitation. In this case, the total combined resources are about 9.8 PWh th , corresponding to a net electricity production of 0.8 PWh el . The uncertainties are shown as relative standard deviations in Figure 6d-f. The general pattern is similar to Figure 4b, but the standard deviation is substantially higher, ranging from 20 to 60% of the resources. For scenario B, the standard deviation is additionally increased along the modeled fault zones and for scenario C along the boundaries of the major basement units.

LCOE and Geothermal Reserves
An economic evaluation was performed based on the calculated recoverable geothermal resources as displayed in Figure 6c and the properties of the reference power plant in Table 2. From this, the cost-supply diagram in Figure 7a is obtained, where the LCOE is plotted against the respective cumulative net electricity production. The curve approaches the geothermal resources for large LCOE. The economically exploitable fraction is obtained by comparing the LCOE with the mean revenues for electricity. Assuming a project lifetime of 25 years (Eyerer et al., 2020), these amount to about 21.6 ct/ kWh el , as the feed-in tariff of 25.2 ct/kWh el is paid for 20 years, and afterward a market revenue of about 7.1 ct/kWh el can be expected (AEE, 2013). Accordingly, the geothermal reserves for the reference case amount to about 300 TWh el or 65% of the resources. Note that these assumed revenues are valid for the situation prior to the energy crisis resulting from the Ukraine war. In the first months of 2022, a significant increase in energy prices oc-curred, but due to the extreme volatility of the market, no reliable forecast of prices can currently be made.
The technical and the economic-financial framework may differ considerably from the reference configuration for specific projects. For this reason, the influence of nine parameters on the LCOE and geothermal reserves was systematically investigated in a sensitivity analysis (Fig. 7b). Detailed cost-supply diagrams for each parameter are given in Figure 8. In general, there is a high degree of agreement with the results of Eyerer et al. (2020), who conducted a similar analysis for hydrothermal reservoirs in Germany.
The maximum drilling depth has the highest positive effect on the total available geothermal resources. It is considered that with deeper wells, also the open hole section becomes longer and therefore more heat can be produced from the reservoir. Furthermore, the fluid temperature increases with depth and thus also the heat content. If the basement is developed to a depth of 6 km instead of 5 km in the entire Northern URG, the reserves are about 110% higher. However, while drilling costs, risk, and technical challenges increase exponentially with depth, the permeability of the basement decreases, making wells deeper than 5 km still rather unlikely in the near future. Recovery factor and the lateral extent of the reservoir, which essentially determine the total amount of heat that a doublet can potentially produce, also positively influence the reserves. They strongly depend, among other factors, on the hydro-mechanical properties of the fracture network, well design, and stimulation effectiveness. A higher recovery factor primarily increases the available recoverable resources and decreases the LCOE, leading to an improved ratio of reserves to resources. The extent of the reservoir does not affect the overall amount of recoverable heat, but with larger reservoir volume, the ratio of reserves to resources increases. Accordingly, geothermal reserves can be significantly increased in the future with improved reservoir enhancement technologies or customized well design patterns with more wells than just the standard doublet and multilaterals in the reservoir section.
Average electricity revenues during the project period are particularly sensitive to political changes in the energy market. Feed-in tariffs are regularly reassessed, hence reductions are possible in the coming years. Conversely, energy prices are generally very volatile, so electricity revenues can fall or  Table 2 for details of these parameters.. rise steeply on the free energy market. An increase in average revenues of about 20% would result in about 15% higher reserves, while a 20% decrease would result in 50% lower reserves. In addition, the impact of investment costs is significant. A 20% increase in costs means halving the reserves compared to the reference case. Potentials for cost reductions exist for the ORC and Kalina technologies on the one hand and drilling costs on the other hand.
Interestingly, project lifetime seems to have a similar impact on reserves as investment costs. The reason for this is that with increasing time, the average revenues decrease since feed-in tariffs are only paid for a fixed period. Therefore, from an economic point of view, limiting the project lifetime to 20 years is recommended. The discount rate reflects the risk assessment of geothermal projects. As more power plants are successfully built, this rate may decrease, resulting in lower power production costs and higher reserves. The thermal full load hours and thermal revenues for the assumed powerto-heat ratio of 0.5 only slightly affect the reserves.

Limitations of the Volumetric Approach
For the first time, a techno-economic assessment of geothermal resources was carried out for the Variscan basement in the Northern URG, which produced a high-resolution map (500 9 500 m 2 cell size) of the resources including uncertainties (Fig. 6). That said, it is important to note that the approach still significantly simplifies the geothermal system. The purely static model neglects all temporal changes in the reservoir during operation, as they are difficult to quantify on the regional scale. This includes the temperature decline during production, the replenishment of the resource by conduction, convection and radiogenic heat production, as well as changes in permeability due to mineral scaling or thermal contraction of the host rock. As a result, the recovery factor is not constant over the project lifetime (Grant, 2018), but reliable data on the temporal behavior in EGS are still rare.
Furthermore, the fracture network properties are excluded from the regional model, which greatly impact the flow rate in geothermal wells and thus affect the profitability of the power plant. A simplified well design was assumed, but depending on the local structural framework, the actual well path and completion may vary significantly. Finally, the cost components and other properties for deep geothermal power plants compiled in Tables 1 and 2 are also subject to large variations. A sensitivity analysis therefore compared the individual influence on the LCOE and geothermal reserves.
Consequently, the results presented are mainly indicative of areas with high potential for geothermal use and cannot replace a detailed site-specific investigation prior to drilling. The resource maps are an important basis for a first site selection and the planning of tailored exploration measures. In this context, the integration of 3D seismic datasets is crucial to locate all potentially permeable faults by tracing them from sedimentary horizons down into the basement. Additional information for the reservoir characterization can be provided, for example, by electromagnetic and gravity surveys as well as thermal and hydrochemical monitoring (Bä r et al., 2021a). Based on a detailed local 3D model, the temporal thermo-hydro-mechanical (THM) processes relevant for the reservoir stimulation and sustainable operation of the geothermal doublet should then be simulated by high-resolution coupled numerical models. Thereby, the petrophysical rock properties can be derived from comprehensive laboratory investigations and corrected for reservoir conditions (e.g., Bä r et al., 2020). Hydrogeological properties of the reservoir may be approximated by combining the few hydraulic test data in the crystalline basement with discrete fracture network models (e.g., Koike et al., 2015;Mahmoodpour et al., 2021). Information on the fracture network properties is provided by sparse boreholes (Afshari et al., 2019;Glaas et al., 2021) but more importantly by analog studies (Bossennec et al., 2021;Bossennec et al., 2022;Frey et al., 2022b). Of particular interest is the influence of the fault damage zone thickness and permeability on the productivity of a geothermal doublet, which can be systematically investigated in high-resolution local models. With a better understanding of the hydraulic reservoir characteristics, it is then possible to estimate the time period over which geothermal heat production can potentially be sustained.

Implications for Geothermal Exploration and Utilization
The presented resource assessment demonstrated that the crystalline basement in the Northern URG is an attractive reservoir for deep geothermal exploitation due to the vast resource base of 830 PWh th . However, only a fraction (8.2-16.1 PWh th ) of it is presumably recoverable with currently available technologies. As expected, the largest resources are associated with major fault zones with a high reactivation potential. Note that in some cases, even fault zones with intermediate slip and dilation tendency may exhibit comparatively high permeabilities, as for example, in Bruchsal. Additional indicators for the potential permeability of fault zones are listed by Agemar et al. (2017) but were not further integrated in the regional-scale study. However, these indicators might be of greater use for local exploration when defining potential drilling targets.
Depending on the hydraulic far-field connectivity, reservoir depletion usually progresses faster than heat resupply through pure conduction and radiogenic heat production. Therefore, after the project lifetime, a rest period is required to recharge the geothermal resources for sustainable utilization of the basement. Conservative estimates suggest a period of about 1000 years (Paschen et al., 2003;Eyerer et al., 2020), which would correspond to a maximum annual heat production of about 8-16 TWh th /year. For comparison, the total annual heat demand in the Northern URG is about 30.6 TWh th /year (Strozyk et al., 2021). Pritchett (1998), on the other hand, reports that 90% of the heat will be renewed after three times the production period. Thus, with a 25-year lifetime, an annual heat production of about 80-160 TWh th /year would be possible. Tester et al. (2006) argue that production can be considered renewable if less than 10% of the heat is extracted annually. This corresponds to a maximum heat production of 0.8-1.6 PWh th /year. More precise insights about the sustainability of geothermal energy will only be provided by longterm experience from commercial EGS projects, but the experiences so far show that the reservoirs in the URG are very well connected to the far-field and can thus be quickly recharged by natural brine circulation on a large scale.
The LCOE analysis has shown that about 65% of the geothermal resources in the Northern URG basement could be exploited economically in the reference case (considering the energy price increase in 2022, the percentage is likely to be even higher). This assessment depends on both the reservoir characteristics and specific investments of the power plants. Drilling is the largest cost factor in EGS projects. The more wells are drilled, the more pro-ven the technologies are and the better the subsurface is explored, resulting in lower risk and shorter drilling times. Therefore, e.g., Limberger et al. (2014) assume that these costs will decrease significantly in the next decades, which is, however, hardly quantifiable. Another way to reduce drilling costs is to drill multiple wells from one drill pad to take advantage of economies of scale. This is particularly useful if a single doublet proves the resources.
Moreover, the effective stimulation of the reservoir plays an essential role. It essentially controls the reservoir volume, the connectivity of the fracture network, and the recovery factor (e.g., Schulte et al., 2010). In this context, the development of novel stimulation methods is required to minimize the induced seismicity (e.g., Amann et al., 2018). In situ experiments addressing these objectives are to be conducted at the planned GeoLaB underground geoscience laboratory on the eastern margin of the URG (Schä tzler et al., 2020).
This study investigates an exclusive geothermal utilization of the basement. Additionally, several sedimentary formations in the Northern URG are suitable for deep geothermal exploitation, namely the Permo-Carboniferous, Buntsandstein, Muschelkalk, and Tertiary sandstones. In particular in the eastern rift area, where the resource base in the basement is low (see Fig. 4), these horizons represent the preferred exploration targets. However, their geothermal potential is not further quantified herein but will be the subject of a subsequent publication (van der Vaart et al., in prep.). Some geothermal projects in the URG pursue a multireservoir approach, for example, in Insheim, Landau and Rittershoffen, where both sedimentary horizons and the basement are exploited together. This reduces the exploration risk, and a higher flow rate as well as a higher heat production can be achieved at the same drilling depth. Co-production of heat and electricity currently provides the most economic benefit, as the feed-in tariff for electricity is relatively high, and the sale of heat provides an additional revenue stream. For direct use of heat from the basement, the LCOH usually exceed 0.03 e/ kWh th , which was higher than average prices in the Northern URG before the Ukraine war. However, as mentioned above, energy prices have increased by up to 90% in the first months of 2022 (Destatis, 2022a), making the direct use of heat from the basement an economically viable solution as of now.
In addition to conventional EGS projects, interest in closed-loop geothermal systems (CLGS) has increased in recent years, where heat transfer from the reservoir to the working fluid is based solely on conduction (Beckers et al., 2022). Various designs exist in this regard, including coaxial pipes or large subsurface heat exchangers consisting of multiple lateral horizontal wells (Winsloe et al., 2021). These approaches include several advantages in that they are independent of subsurface hydraulic properties, no stimulation is required, less corrosion and scaling occur, and plants can be scaled as required. CLGS can also be retrofitted into existing geothermal, oil, and gas wells. Disadvantages include high drilling costs and comparatively rapid temperature decline during operation. The companies Eavor and GreenFire Energy have already realized initial demonstration projects in North America (Higgins et al., 2019;Toews et al., 2020). Nevertheless, the technical and economic feasibility on a full industrial scale has yet to be proven. In the URG, there is no direct competition to the classical EGS projects since the CLGS target mostly undisturbed rock units, which are not the primary drill sites for open systems.

Socio-economic-environmental Potential
Apart from the geological, technical and financial parameters, the conditions at the surface must always be considered when developing geothermal power plants, as these can have both positive and negative effects on the feasibility of such projects. Frey et al. (2022c) mapped the socio-economic-environmental potential for deep geothermal energy in the URG. The three considered dimensions consist in total of nine individual indicators. The social dimension includes the population-related aspects of population density, heat demand, social level, and acceptance of renewable energy forms. The economic dimension comprises the gross domestic product, the existing district heating network infrastructure and the public debt rate. Finally, the environmental dimension examines the aspects of land access and greenhouse gas emissions. The different quantifiable components were first normalized and then joined into a composite index via an unweighted generalized mean (for a detailed description of the approach see Frey et al., 2022c).
In the Northern URG, the socio-economic-environmental index exhibits a wide distribution (Fig. 9), with a clear trend between urban and rural areas. The highest potential is found in the large cities and their peripheral districts, which can be attributed to the concentration of population and industry as well as the extensive heat supply infrastructure (e.g., district heating networks of Frankfurt, Darmstadt, Mannheim and Heidelberg). Here, natural and landscape conservation areas can be limiting factors that make access to geothermal resources difficult or in rare cases completely impossible, since drilling below these areas is still possible by deviated wells at sufficient depth. These are mainly concentrated along the Rhine. The more sparsely populated and agriculturally dominated areas west of the Rhine have the lowest overall socio-economic-environmental potential.
The comparison of the recoverable resources and mean LCOE (Fig. 9a, b) with the socio-economic-environmental index (Fig. 9c) reveals a high degree of correlation. In particular in the Rhine-Neckar area, favorable conditions in the subsurface and on the surface are combined. Likewise, the areas along the southwestern margin of the Northern URG and around Darmstadt are positively rated. Although the socio-economic conditions are good in the Frankfurt metropolitan region, the crystalline basement is less suitable for geothermal utilization.
Again, this regional evaluation is mostly indicative for locations with an overall high potential. In particular, the social acceptance of geothermal plants can vary considerably at the local level and depends primarily on the quality of communication and the integration of the population by the project developers. The specific decision-making for a drill site represents a multidimensional problem that has to be solved via analytical hierarchical processes (Raos et al., 2022). This involves carefully weighing the conditions in the subsurface against the conditions at the surface.

CONCLUSIONS
A detailed techno-economic resource assessment was conducted for the crystalline basement in the Northern URG. The following conclusions are drawn from these investigations: The Northern URG basement exhibits a vast geothermal resource base of approximately 830 PWh th , of which between 8.2 and 16.1 PWh th are potentially recoverable with current technologies. Geothermal reserves amount to at least 65% of the resources. The substantial energy price increases at the beginning of 2022 result in an even higher share of resources that could be used economically. However, due to the energy market volatility, it is currently impossible to make a reliable forecast of price developments. The geothermal resources are mainly tied to fault zones with high slip and dilation tendencies located in predominantly granitoid basement units. The largest uncertainties are found along the eastern rift margin, where the top of basement is deepest. The impact of various CHP plant properties on the LCOE was investigated using a sensitivity analysis. Significant influences were found for the maximum drilling depth, the recovery factor, the reservoir development, the investment costs, the project lifetime, the discount rate and the average revenue for electricity. In contrast, the thermal full load hours and the average revenue from the sale of heat have a lower significance. The socio-economic-environmental potential for deep geothermal energy reveals a clear urbanrural trend. A comparison with the resources shows a high degree of agreement in the Rhine-Main area, west of Darmstadt, the southeastern areas around Mannheim, and the southwestern rift margin close to Landau The techno-economic and socio-economic-environmental analyses are primarily indicative of areas with high potential. They are not a substitute for local targeted exploration activities and multidimensional decision-making processes. In-depth studies of the dynamic behavior of reservoir properties during production need to be carried out using local coupled THM models

ACKNOWLEDGMENTS
We thank the three anonymous reviewers for their very helpful comments on this manuscript. The Figure 9. Comparison of (a) the recoverable geothermal resources (per 500 9 500 m 2 raster cell, Scenario C, see Fig. 6), (b) the LCEO for the reference power plant (see. Table 2) and (c) the socio-economical-environmental potential in the Northern URG. The unitless index ranges from 0 (worst conditions) to 100 (ideal conditions). research was funded by the Interreg NWE Program through the Roll-out of Deep Geothermal Energy in North-West Europe (DGE-ROLLOUT Project, w ww.nweurope.eu/DGE-Rollout). The Interreg NWE Program is part of the European Cohesion Policy and is financed by the European Regional Development Fund (ERDF). We acknowledge support by the Deutsche Forschungsgemeinschaft (DFG-German Research Foundation) and the Open Access Publishing Fund of Technical University of Darmstadt, through the DEAL project. Petrel was kindly provided by Schlumberger for this research.

FUNDING
Open Access funding enabled and organized by Projekt DEAL.

DATA AVAILABILITY
The discretized and parameterized model used for the geothermal resource assessment can be found in the TUDatalib repository, https://doi.org/ 10.48328/tudatalib-898.

Conflict of Interest
The authors have no relevant financial or non-financial interests to disclose.

OPEN ACCESS
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecom mons.org/licenses/by/4.0/.