USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geothermics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11Linked to original sources
USGS1990

Geothermal systems within the Mammoth Corridor in Yellowstone National Park and the adjacent Corwin Springs KGRA

A study of potential impacts of geothermal development in the Corwin Springs KGRA north of Yellowstone Park on thermal springs within the Park is being conducted by the U.S. Geological Survey. Thermal waters in the KGRA and at Mammoth Hot Springs, located 13 km inside the Park boundary, are high in bicarbonate and sulfate and are actively depositing travertine. These similarities and the existence of numerous regional-scale structural and stratigraphic features that could provide conduits for fluid flow at depth indicate a possible cause for concern. The objectives of this study include delineations of any hydrologic connections between these thermal waters, the level of impact of geothermal development in the event of such connections, and mitigation measures to minimize or eliminate adverse impacts. The study involves a number of geochemical, geophysical, geologic, and hydrologic techniques, but does not include any test drilling. Preliminary results suggest that thermal waters at Bear Creek Springs may contain a component of water derived from Mammoth but that thermal waters at La Duke Hot Spring do not. The total rate of thermal water that discharges in the area proposed for geothermal development (near La Duke) has been determined; restricting the net production of thermal water to rates less than this total could provide a satisfactory margin of safety for development.

Conference Paper↗
USGS1991

Origins of acid fluids in geothermal reservoirs

Acid fluids in geothermal reservoirs are rare. Their occurrence in geothermal systems associated with recent volcanism (Tatun, Sumikawa, Miravalles) probably indicates that the geothermal reservoir fluid was derived from volcanic fluid incompletely neutralized by reaction with feldspars and micas. Superheated steam containing HCl (Larderello, The Geysers) forms acid where it condenses or mixes with liquid at moderate temperatures (<300??C). The origin of steam with HCl is reaction of NaCl solid with rock minerals at high temperatures (>325??C). Cryptoacidity occurs at Los Humeros where HCl acidity is formed and neutralized without reaching the surface.

Conference Paper↗
USGS1998

Tectonic controls on fault-zone permeability in a geothermal reservoir at Dixie Valley, Nevada

To determine factors controlling permeability variations within and adjacent to a fault-hosted geothermal reservoir at Dixie Valley, Nevada, we conducted borehole televiewer observations of wellbore failure (breakouts and cooling cracks) together with hydraulic fracturing stress measurements in six wells drilled into the Stillwater fault zone at depths of 2 to 3 km. Measurements in highly permeable wells penetrating the main geothermal reservoir indicate that the local orientation of the least horizontal principal stress, Shmin, is nearly optimal for normal faulting on the Stillwater fault. Hydraulic fracturing tests from these wells further show that the magnitude of Shmin is low enough to lead to frictional failure on the Stillwater and nearby subparallel faults, suggesting that fault slip is responsible for the high reservoir productivity. Similar measurements were conducted in two wells penetrating a relatively impermeable segment of the Stillwater fault zone, located approx. 8 and 20 km southwest of the geothermal reservoir (wells 66-21 and 45-14, respectively). The orientation of Shmin in well 66-21 is near optimal for normal faulting on the nearby Stillwater fault, but the magnitude of Shmin is too high to result in incipient frictional failure. In contrast, although the magnitude of Shmin in well 45-14 is low enough to lead to normal faulting on optimally oriented faults, the orientation of the Stillwater fault near this well is rotated by approx. 40?? from the optimal orientation for normal faulting. This misorientation, coupled with an apparent increase in the magnitude of the greatest horizontal principal stress in going from the producing to nonproducing wells, acts to inhibit frictional failure on the Stillwater fault zone in proximity to well 45-14. Taken together, data from the nonproducing and producing wells thus suggest that a necessary condition for high reservoir permeability is that the Stillwater fault zone be critically stressed for frictional failure in the current stress field.

Conference Paper↗
USGS2002

Fracture Development within the Karaha-Telaga Bodas Geothermal Field, Indonesia

Karaha-Telaga Bodas is a partially vapor-dominated geothermal system located in an active volcano in western Java. More than 2 dozen geothermal wells have been drilled to depths of 3 km. Detailed paragenetic and fluid-inclusion studies have defined liquid-dominated, transitional and vapor-dominated stages in the evolution of this system. The liquid-dominated stage was initiated by shallow magma intrusion into the base of the volcanic cone. Lava and pyroclastic flows capped a geothermal system. The uppermost andesite flows were only weakly fractured due to the insulating effect of the intervening altered pyroclastics, which absorbed the deformation. Shear and tensile fractures were filled with carbonates at shallow depths and by quartz, epidote and actinolite at depths and temperatures over 1km and 300??C. The system underwent numerous local cycles of overpressuring, which are marked by subhorizontal tensile fractures, anastomosing tensile fractures and implosion breccias. The development of the liquid system was interrupted by a catastrophic drop in fluid pressures. As the fluids boiled in response to this pressure drop, chalcedony and quartz were deposited in fractures having the largest apertures and steep dips. The orientations of these fractures indicate that the escaping overpressured fluids used the shortest possible paths to the surface. Vapor-dominated conditions were initiated within a vertical chimney over the still hot intrusion. As pressures declined these conditions spread outward. Downward migration of the chimney occurred as the intrusion cooled and the brittle-ductile transition migrated to greater depths. Condensate that formed at the top of the vapor-dominated zone percolated downward and lowsalinity meteoric water entered the marginal parts of the system. Calcite, anhydrite, and fluorite precipitated in fractures upon heating. A progressive sealing of the fractures occurred, resulting in the downward migration of the cap rock. In response to decreasing pore pressures in the expanding vapor zone, the fracture system within the vapor-dominated reservoir progressively collapsed, leaving only residual permeability, with apertures supported by asperities or propping breccia. In places, the fractures have completely collapsed where normal stresses acting on the fracture walls exceeded the compressive strength of the wall rock.

Conference Paper↗
USGS2002

Mantle helium and carbon isotopes in Separation Creek geothermal springs, Three Sisters area, central Oregon: Evidence for renewed volcanic activity or a long term steady state system?

Here we present the helium and carbon isotope results from the initial study of a fluid chemistry-monitoring program started in the summer of 2001 near the South Sister volcano in central Oregon. The Separation Creek area which is several miles due west of the volcano is the locus of strong crustal uplift currently occurring at a rate of 4-5 cm/yr (Wicks, et. al. , 2001). Helium [R C /R A = 7.44 and 8.61 R A (R C /R A = ( 3 He/ 4 He) sample-air corrected /( 3 He/ 4 He) air ))] and carbon (δ 13 C = -11.59 to -9.03 ‰ vs PDB) isotope data and CO 2 / 3 He (5 and 9 x 10 9 ) show that bubbling cold springs in the Separation Creek area near South Sister volcano carry a strong mantle signal, indicating the presence of fresh basaltic magma in the volcanic plumbing system. There is no evidence though, to directly relate this signal to the crustal uplift that is currently taking place in the area, which started in 1998. The geothermal system in the area is apparently much longer lived and shows no significant changes in chemistry compared to data from the early 1990s. Hot springs in the area, which are relatively far removed from the volcanic edifice, do not carry a strong mantle signal in helium isotope ratios (2.79 to 5.08 R A ), unlike the cold springs, and also do not show any significant changes in helium isotope ratios compared to literature data for the same springs of over two decades ago. The cold springs of the Separation Creek area form a very diffuse but significant low temperature geothermal system, that should, due to its close vicinity to the center of up uplift, be more sensitive to changes in the deeper volcanic plumbing system than the far removed hot springs and therefore require much more study and consideration when dealing with volcano monitoring in the Cascade range or possibly with geothermal exploration in general.

Oregon↗
USGS2003

Reservoir depletion at The Geysers geothermal area, California, shown by four-dimensional seismic tomography

Intensive geothermal exploitation at The Geysers geothermal area, California, induces myriads of small-magnitude earthquakes that are monitored by a dense, permanent, local seismometer network. Using this network, tomographic inversions were performed for the three-dimensional Vp and Vp/Vs structure of the reservoir for April 1991, February 1993, December 1994, October 1996, and August 1998. The extensive low-Vp/Vs anomaly that occupies the reservoir grew in strength from a maximum of 9% to a maximum of 13.4% during the 7-year study period. This is attributed to depletion of pore liquid water in the reservoir and replacement with steam. This decreases Vp by increasing compressibility, and increases Vs because of reduction in pore pressure and the drying of argillaceous minerals, e.g., illite, which increase the shear modulus. These effects serendipitously combine to lower Vp/Vs, resulting in a strong overall effect that provides a convenient tool for monitoring reservoir depletion. Variations in the Vp and Vs fields indicate that water depletion is the dominant process in the central part of the exploited reservoir, and pressure reduction and mineral drying in the northwest and southeast parts of the reservoir. The rate at which the Vp/Vs anomaly grew in strength in the period 1991-1998 suggests most of the original anomaly was caused by exploitation. Continuous monitoring of Vp, Vs, and Vp/Vs is an effective geothermal reservoir depletion monitoring tool and can potentially provide information about depletion in parts of the reservoir that have not been drilled.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Geothermometer calculations for geothermal assessment

Geothermal exploration programs have relied on the calculation of geothermometers from hot spring chemistry as an early estimation of geothermal reservoir temperatures. Calibration of the geothermometers has evolved from experimental determinations of mineral solubility as a function of temperature to calibration from analyses of water chemistry from known depths and temperatures in thermal wells. Most of the geothermometers were calibrated from analyses of sodium-chloride type waters, and the application of some geothermometers should be restricted to waters of the chemical types that were used in their calibration. Chemical analyses must be determined to be reliable before they are used to calculate geothermometers. The USGS Geothermal Resource Assessment will rely on the silica geothermometer developed by Giggenbach that approximates the transition between chalcedony at 20??C and quartz at 200??C. Above 200??C, the assessment will rely on the quartz geothermometer. In addition, the assessment will also rely on the potassium-magnesium geothermometer.

Conference Paper↗
USGS1976

Audio-magnetotelluric methods in reconnaissance geothermal exploration

An audio-magnetotelluric (AMT) system has been developed by the U.S. Geological Survey for low-cost reconnaissance exploration of geothermal regions. This is an electromagnetic sounding technique in which the scalar or Cagniard resistivity is computed at 12 frequencies logarithmically spaced from 7.5 to 18 600 Hz. Our system uses natural source fields except at the two upper frequencies of 10 200 and 18 600 Hz where artificial VLF sources are available. As a reconnaissance technique we use AMT surveys in conjunction with regional gravity, magnetic, and telluric surveys. The exploration depth is a function of the resistivities of the lithologic section, but typically ranges from the surface to 0.2 km in low-resistivity areas and to greater than 2 km in high-resistivity regions. Results of the initial reconnaissance AMT surveys provide a rational basis for deciding on the extent of costlier follow-up surveys. As part of the U.S. Geological Survey geothermal program, surveys were conducted in Long Valley and Surprise Valley, California; the Vale, Ore-Weiser, Idaho region; and Bruneau-Grand View, Raft River, and Island Park regions of Idaho. AMT surveys in five additional known geothermal resource areas (KGRA's) have been scheduled for completion by May 1975. In the Raft River and Bruneau-Grand View regions and Long Valley, follow-up electrical surveys substantiated the effectiveness of the AMT technique for reconnaissance surveying.

Conference Paper↗
USGS1982

Geothermal systems of the Cascade Range

In the central and southern Cascade Range, plate convergence is oblique, and Quaternary volcanism produces mostly basalt and mafic andesite; large andesite-dacite composite volcanoes and silicic dome fields occur in restricted areas of long-lived igneous activity. To the north, plate convergence is normal, producing widely spaced centers in which mafic lavas are minor. Most Cascade volcanoes are short-lived and unlikely to be underlain at shallow levels by large magma bodies that could support high-temperature geothermal systems. Such systems are known, however, near Meager Mountain, at Newberry Volcano, and near Lassen Peak. Persistent fumaroles occur on several major composite volcanoes, but drilling to date has been insufficient to determine whether exploitable geothermal reservoirs occur at depth. Thermal springs away from the major volcanic centers are few and generally inconspicuous. However, significant geothermal systems along and west of the Cascade Range may well be masked by abundant cold ground water.

California, Oregon, Washington↗
USGS2015

Natural or Induced: Identifying Natural and Induced Swarms from Pre-production and Co-production Microseismic Catalogs at the Coso Geothermal Field

Increased levels of seismicity coinciding with injection of reservoir fluids have prompted interest in methods to distinguish induced from natural seismicity. Discrimination between induced and natural seismicity is especially difficult in areas that have high levels of natural seismicity, such as the geothermal fields at the Salton Sea and Coso, both in California. Both areas show swarm-like sequences that could be related to natural, deep fluid migration as part of the natural hydrothermal system. Therefore, swarms often have spatio-temporal patterns that resemble fluid-induced seismicity, and might possibly share other characteristics. The Coso Geothermal Field and its surroundings is one of the most seismically active areas in California with a large proportion of its activity occurring as seismic swarms. Here we analyze clustered seismicity in and surrounding the currently produced reservoir comparatively for pre-production and co-production periods. We perform a cluster analysis, based on the inter-event distance in a space-time-energy domain to identify notable earthquake sequences. For each event j, the closest previous event i is identified and their relationship categorized. If this nearest neighbor&rsquo;s distance is below a threshold based on the local minimum of the bimodal distribution of nearest neighbor distances, then the event j is included in the cluster as a child to this parent event i. If it is above the threshold, event j begins a new cluster. This process identifies subsets of events whose nearest neighbor distances and relative timing qualify as a cluster as well as a characterizing the parent-child relationships among events in the cluster. We apply this method to three different catalogs: (1) a two-year microseismic survey of the Coso geothermal area that was acquired before exploration drilling in the area began; (2) the HYS_catalog_2013 that contains 52,000 double-difference relocated events and covers the years 1981 to 2013; and (3) a catalog of 57,000 events with absolute locations from the local network recorded between 2002 and 2007. Using this method we identify 10 clusters of more than 20 events each in the pre-production survey and more than 200 distinct seismicity clusters that each contain at least 20 and up to more than 1000 earthquakes in the more extensive catalogs. The cluster identification method used yields a hierarchy of links between multiple generations of parent and offspring events. We analyze different topological parameters of this hierarchy to better characterize and thus differentiate natural swarms from induced clustered seismicity and also to identify aftershock sequences of notable mainshocks. We find that the branching characteristic given by the average number of child events per parent event is significantly different for clusters below than for clusters around the produced field.

Conference Paper↗
USGS1978

Helium in soil gases of the Roosevelt Hot Springs Known Geothermal Resource Ares, Beaver County, Utah

Soil samples were collected in two parallel traverses across the Dome fault zone of the Roosevelt Hot Springs Known Geothermal Resource Area. The samples were sealed in air-tight aluminum cans, and the soil gas was allowed to equilibrate with the atmospheric air in the cans. Gas from the cans was analyzed by mass spectrometry. Samples collected over faults contained anomalously high concentrations of helium. Samples collected close to a geothermal well 884 m deep contained more helium than samples collected near another geothermal well 1370 m deep.

Utah↗
USGS1976

Geothermal energy in the United States; Part II, Assessment of resources

Geothermal energy-from heat deep inside the Earth- is a vast potential source of power. This article is the second part of a series on geothermal energy, the first part of which was in volume 8, number 1, of the Earthquake Information Bulletin (January-February 1976). Part 1 of this series described the categories of the geothermal resource base.

Earthquake Information Bulletin (USGS)↗
USGS2016

Analysis and interpretation of stress indicators in deviated wells of the Coso Geothermal Field

Characterizing the tectonic stress field is an integral part of the development of hydrothermal systems and especially for enhanced geothermal systems (EGS). With a well characterized stress field the propensity of fault slip on faults with known location and orientation can be identified. Faults that are critically oriented for faulting with respect to the stress field are known to provide natural fluid pathways. A high slip tendency makes a fault a likely candidate for reactivation during the creation of an EGS. Similarly, the stress state provides insight for the potential of larger, damaging earthquakes should extensive portions of well-oriented, larger faults be reactivated. The analysis of stress indicators such as drilling-induced fractures and borehole breakouts is the main tool to infer information on the stress state of a geothermal reservoir. The standard procedure is applicable to sub-vertical wellbore sections and highly deviated sections have to be discarded. However, in order to save costs and reduce the environmental impact most recent wells are directionally drilled with deviations that require appropriate consideration of the deviated trajectory. Here we present an analysis scheme applicable to arbitrary well trajectories or a combination of wells to infer the stress state. Through the sampling of the stress tensor along several directions additional information on the stress regime and even relative stress magnitudes can be obtained. We apply this method on image logs from the pair of wells 58-10 and 58A-10 that were drilled from the same well pad. Both wells have image logs of about 2km of their trajectories that are separated by less than 300m. For both wells we obtain a mean orientation of S Hmax of N23° with large standard deviations of locations of stress indicators of 24° and 26°, respectively. While the local stress direction is highly variable along both wells with dominant wavelengths from around 50 to 500m, the mean directions are very consistent and also agree with previous stress estimates in the eastern part of the Coso Geothermal Field. In order to obtain a reliable estimation of the stress orientation in this setting, it is necessary to sample the stress field on an interval long to capture several of the dominant wavelengths.

California↗
USGS2014

Micro-seismicity and seismic moment release within the Coso Geothermal Field, California

We relocate 16 years of seismicity in the Coso Geothermal Field (CGF) using differential travel times and simultaneously invert for seismic velocities to improve our knowledge of the subsurface geologic and hydrologic structure. We expand on our previous results by doubling the number of relocated events from April 1996 through May 2012 using a new field-wide 3-D velocity model. Relocated micro-seismicity sharpens in many portions of the active geothermal reservoir, likely defining large-scale fault zones and fluid pressure compartment boundaries. However, a significant fraction of seismicity remains diffuse and does not cluster into sharply defined structures, suggesting that permeability is maintained within the reservoir through distributed brittle failure. The seismic velocity structure reveals heterogeneous distributions of compressional (Vp) and shear (Vs) wave speed, with Vs generally higher in the Main Field and East Flank and Vp remaining relatively uniform across the CGF, but with significant local variations. The Vp/Vs ratio appears to outline the two main producing compartments of the reservoir at depths below mean ground level of approximately 1 to 2.5 km, with a ridge of relatively high Vp/Vs separating the Main Field from the East Flank. Detailed analyses of spatial and temporal variations in earthquake relocations and cumulative seismic moment release in the East Flank reveal three regions with persistently high rates of seismic activity. Two of these regions exhibit sharp, stationary boundaries at the margins of the East Flank that likely represent barriers to fluid flow and advective heat transport. However, seismicity and moment release in a third region at the northern end of the East Flank spread over time to form an elongated NE to SW structure, roughly parallel both to an elongated cluster of seismicity at the southern end of the East Flank and to regional fault traces mapped at the surface. Our results indicate that high-precision relocations of micro-seismicity and simultaneous velocity inversions in conjunction with mapping of seismic moment release can provide useful insights into subsurface structural features and hydrologic compartmentalization within the Coso Geothermal Field.

California↗
USGS2019

Satellite observations of surface deformation at the Coso Geothermal Field, California

Surface deformation time series and rates are identified at the Coso Geothermal Field (CGF) and surrounding areas by applying interferometric synthetic aperture radar (InSAR) to satellite scenes from Envisat (June 2004 ̶ October 2010) and Sentinel (November 2014 – April 2018). The measurements are done in the line of sight (LOS) to each satellite, within an area of size ~450 km2, at the locations of hundreds of thousands permanent and distributed scatterers. Thirty descending (satellite moves north to south) and 45 ascending (south to north) images were used from Envisat, and 63 descending and 65 ascending from Sentinel. A decomposition into average vertical and east horizontal components is also performed in more than 35,000 100-m pixels where both types of LOS measurements are available. The main observations at CGF include: (1) a subsidence area of size ~70 km2, with a maximum subsidence of –27.6 mm/year for the Envisat period and lower maximum subsidence of –19.1 mm/year for the Sentinel period; (2) eastward movements in the western part of the subsidence area, with Envisat maximum of +23.9 mm/year and a lower Sentinel maximum of +15.9 mm/year; (3) westward displacements in the eastern part of the subsidence area, with Envisat maximum of ̶ 14.2 mm/year and Sentinel maximum of –11.9 mm/year; (4) very good agreement of the InSAR observations with leveling survey data; (5) earthquake clusters in the subsidence area and hypocentral cross-sections showing clusters at various depths and migration in time; and (6) good predictions of the overall geothermal resource, based on poroelastic modeling using both leveling and InSAR data. The ultimate goal of the project is to provide geothermal operators with tools that can be used in reservoir management.

California↗
USGS2020

Geothermal play fairway analysis of the Sou Hills, northern Nevada: A major quaternary accommodation zone in the Great Basin region

To facilitate discovery and development of blind geothermal systems in the Great Basin region, as well as assessment of known systems with surface hot springs, a play fairway approach was developed to evaluate and integrate multiple geologic and geophysical parameters for permeability and heat. Phase 1 of the project produced a geothermal potential map of 96,000 km2 of Nevada. This analysis yielded 24 highly favorable locations with particularly high fairway scores, from which five promising sites were selected for detailed studies. The Sou Hills was chosen as a detailed study area due to a complex structural setting, plethora of Quaternary normal faults, and relatively high rates of both regional strain and slip on Quaternary faults. Our analysis of the Sou Hills utilized: 1) detailed geologic mapping of ~60 km2, 2) reconnaissance mapping of >200 km2, 3) analysis of Quaternary faults, 4) detailed logging of cuttings from four, previously drilled wells (~2,000 m) and integration of ~5,500 m of existing logs from nine other wells, 5) a new gravity survey (355 stations), 6) LiDAR acquisition for 290 km2, 7) a shallow temperature survey (82 stations), 8) interpretation of seven seismic reflection profiles, 9) slip and dilation tendency analyses, and 10) geochemical analyses of 23 water samples. Integration of these datasets shows that the Sou Hills occupies a major accommodation zone between oppositely dipping Quaternary normal fault systems. As such, the area is characterized by multiple, closely-spaced, west- and east-dipping Quaternary normal faults, many of which are favorably oriented for slip and dilation in the current stress regime. Seismic reflection and gravity data show an extensional anticline in northernmost Dixie Valley directly south of the Sou Hills. The anticline marks a zone of multiple intersecting, oppositely dipping normal faults and represents a particularly favorable site for potentially hosting a geothermal system.

Nevada↗
USGS2025

The influence of pre-existing structures on geothermal springs: Inferences from potential field mapping in Surprise Valley, CA and other sites In the northwestern Great Basin

Surprise Valley, located in the northwestern Great Basin, is an asymmetric extensional basin that marks a major tectonic transition between the relatively un-extended volcanic Modoc Plateau to the west, and the Basin and Range to the east that has undergone 10-15% extension. In addition, it sits just north of the Walker Lane which accommodates up to 20% of dextral slip associated with Pacific-North American plate interactions. Thermal springs issue from eight areas within Surprise Valley. Most of these occur within the basin and are not situated on the main basin forming range-front faults. As a result, efforts to resolve the structural setting of the valley’s hydrothermal system have relied on geophysics to characterize basin structure and geology. Extensive efforts to map the basin with ground and airborne magnetics have revealed a >35 km-long linear, intra-basin magnetic high, interpreted as a buried dike swarm. Geothermal springs on the eastern side of the valley, including Seifert hot springs, Leonards hot springs, and Surprise Valley hot springs (SVHS), are all situated along the magnetic high and occur at local breaks and bends in the anomaly, suggesting that fracture permeability is enhanced along the feature and particularly at these discontinuities. Recent studies, including drilling over the anomaly near SVHS that likely intersected dike material, as well as subsequent mapping and sampling of dikes outcropping along the anomaly on the playa surface south of SVHS, confirm (as previously inferred) that mafic intrusives are the principal source of the anomaly. Similar interpretations made in two other valleys (in southern Oregon and northwestern Nevada), where inferred intra-basin dikes appear to be spatially correlated with hot springs or prospective geothermal resource areas, suggest that the impact of pre-existing basement structure on hydrothermal activity may pertain more generally to other hydrothermal settings throughout the Great Basin. If so, efforts to map basement may enhance understanding structural controls on some geothermal systems. Furthermore, similarities across these disparate sites suggest that magmatism may play a much larger role in accommodating extension and influencing basin evolution across the western Great Basin than previously recognized.

Conference Paper↗
USGS2025

Exploring Martian geothermal and liquid water potential with basin modeling

Assessing the potential for geothermal energy and liquid water presence in the Martian subsurface is crucial for future exploration and habitability studies. In this work, we employed comprehensive finite element model simulations adapted specifically for Martian conditions to estimate subsurface temperatures and the potential for liquid water at depth within Martian crater basins. Rock and fluid property values for basin fill were carefully adjusted to match Martian gravity, radiogenic heat generation, and compositional characteristics derived from rover analyses, Martian meteorite samples, and orbital spectroscopy data. Multiple modeling scenarios were explored to systematically evaluate end-member cases across critical variables such as heat flow, lithological composition, and average surface temperature. Sensitivity testing revealed that heat flow and average annual surface temperatures are the most important variables. Results were used in calculations based on a database of Martian craters to estimate the temperature of crater fill at depth. Our model results indicate significant potential for sustained liquid water in the subsurface within sedimentary deposits across a range of crater sizes and latitudes. They further suggest that viable geothermal reservoirs likely exist and are potentially accessible for future Martian missions seeking energy sources or exploring astrobiological hypotheses. This study provides a methodological framework for geothermal and hydrological assessments for the subsurface of Mars, contributing to ongoing planetary exploration strategies.

Conference Paper↗