天水及其南北地区温泉分布的地质—地球物理特征

根据天水及其南北地区温泉分布众多的现象,从分析地质、地球物理场着手,探讨温泉形成的区域与深部地质构造背景以及热缘机制。区内近南北向隐伏断裂发育,中酸性侵入岩、碱性玄武喷发岩构成近南北向岩浆构造带;在深部软流圈埋深变浅,岩石圈减薄、上地幔热物质上涌引起莫霍面上隆;中地壳低速、高导层同步出现;地表热流值增高,温泉呈网格状分布。该区温泉是甘肃省隆起断裂对流型地热资源有利开发地段。     

Flow rates in the axial hot springs of the East Pacific Rise (21°N): Implications for the heat budget and the formation of massive sulfide deposits

The first in-situ measurements of flow rates in submarine hot springs (temperatures between 275° and 350°C) were made in the hot springs of the East Pacific Rise at 21°N during November, 1981. The flow rates ranged between 0.7 and 2.4 m/s. We estimate that the total rates at which heat and mass are transported from the three hot spring areas—OBS, National Geographic, and Southwest Vents—are ca. 2.2×10⁸ W and 150 kg of fluid/s respectively. The very high rate of heat loss precludes the existence of hot springs as steady state features of the East Pacific Rise at 21°N. We estimate that hot springs are active for a maximum of ca. 40,000 years. Chemical and/or mechanical clogging of the vents can drastically reduce the lifetime of the hot springs.Using a simple model of plume development [1], we calculated that less than 3% of the sulfide particles that are entrained in the hydrothermal plumes at 21°N settle from the plumes before dispersal by a lateral submarine current at a height of ca. 250 m above the seafloor [2]. The time that is required to form the sulfide mounds that surround the hot springs by the accumulation of sulfide particles that settle from the plumes is ca. 70–85 years. The current rate at which sulfide particles settle from the hydrothermal plumes and accumulate near the hot springs appears to be too small to lead to the formation of large massive sulfide deposits such as the Mavrouni deposit on Cyprus [3], within the time limits that are suggested by our thermal energy balance calculations.     

Plume Mode of Thermal Convection in the Earth’s Mantle

In our previous works, based on numerical models, it was shown that under certain conditions a hot material can rise in portions in the tails of thermal mantle plumes. The spectrum of these pulsations can correspond to the observed spectra of catastrophic hotspot eruptions. Since most of the existing numerical models of thermal convection for the mantle of the present Earth do not reveal these pulsations, in this work, we analyze the physical cause and initiation conditions of pulsations of thermal plumes. The results of a numerical solution of the thermal convection equations for a material with varying parameters in the extended Boussinesq approximation are presented. It is shown how the structure of the convection is transformed with the increase of convection intensity. At the Rayleigh numbers Ra > 10⁶, convection becomes unsteady, and the configuration of the ascending and descending flows changes. The new flow emerging at the mantle bottom acquires a mushroom shape with a head and a tail. After the rise of the plume’s head to the surface, the tail remains in the mantle in the form of a quasi-stationary hot steam. It turns out that at Ra ~ 5 × 10⁷, the thermal mantle plume becomes pulsating and its tail is in fact a heated channel through which the hot material rises in successive portions. At the Rayleigh numbers Ra > 5 × 10⁸, the tail of the thermal plume breaks and the plume becomes a regular conveyor of separate ascending portions of the hot material, which are referred to as thermals. Thus, thermal convection with pulsating plumes takes place at the transitional stage from the regime of quasi-stationary plumes to the regime of thermals.     

On the possibility of crater formation associated with an ascending plume

A hypothetical possibility of a qualitative explanation of large crater formation on the surfaces of the Moon and Mercury is discussed in terms of the concept of thermal mantle plumes. Prerequisites to this hypothesis are revealed under the assumption that the model equation of state of SiO₂ exhibiting an anomaly (a negative coefficient of thermal expansion) in the range of states approximately corresponding to average conditions typical of mantles of minor planets is applicable, in a first approximation, to mantle material. The anomaly reduces the buoyancy of hot plume material in such a way that, under conditions of moderate overheating, only relatively high columns comparable in size to the mantle are capable of ascending from the mantle bottom to the crust; allows cold peripheral material surrounding the hot column to be pushed away; causes compaction of the vertical zone of the contact of the column with the surrounding medium at the first stages after the plume ascent; and leads to compaction of the deep mantle due to the long-term heat supply. Such phenomena can lead to vertical craterlike deformations of the crust in areas of ascending large plumes whose presence can be supposed at early stages of the existence of minor planets. Significant implications of such an anomaly for geophysical processes can also be postulated.     

The dynamical and thermal structure of deep mantle plumes

If the interpretation of the D″ layer at the base of the mantle as a thermal boundary layer, with a temperature increment in the order of 800 K, is correct, then the formation of deep-mantle plumes to vent material from it appears inevitable. We demonstrate quantitatively that the strong temperature dependence of viscosity guides the upward flow into long-lived chimneys that are ～ 20 km in diameter near the base of the mantle and decrease in width with progressive upward softening and partial melting of plume material. The speed of flow up the axis of the plume is correspondingly fast; 1.6 m y^?1 at the base and 4.8 m y^?1 at 670 km depth. Thermal diffusive spreading of a heated plume is compensated by a slow horizontal convergence of mantle material toward the chimney in response to the lower pressure there. This convergence, which contributes only a small increment to the flux of material up the plume, also serves to throttle the flow in the chimney. The global plume mass flux necessary to transport 1.6 × 10¹² W of core heat upward through the mantle is 1.8 × 10⁶ kg s^?1. At its base, plume material is probably still significantly below its solidus or eutectic temperature, but substantial partial melting is very likely as it rises. We speculate that a small fraction of this fluid component eventually emerges at the surface in “hot spots”, with the fate of the remainder being unknown. The behaviour and properties of D″ and of plumes are closely coupled. Not only are plumes a necessary consequence of a thermal boundary layer, but their existence is impossible without that layer.     

Plume capture by divergent plate motions: implications for the distribution of hotspots,geochemistry of mid-ocean ridge basalts,and estimates of the heat flux at the core–mantle boundary

The coexistence of stationary mantle plumes with plate-scale flow is problematic in geodynamics. We present results from laboratory experiments aimed at understanding the effects of an imposed large-scale circulation on thermal convection at high Rayleigh number (10⁶≤Ra≤10⁹) in a fluid with a temperature-dependent viscosity. In a large tank, a layer of corn syrup is heated from below while being stirred by large-scale flow due to the opposing motions of a pair of conveyor belts immersed in the syrup at the top of the tank. Three regimes are observed, depending on the ratio V of the imposed horizontal flow velocity to the rise velocity of plumes ascending from the hot boundary, and on the ratio λ of the viscosity of the interior fluid to the viscosity of the hottest fluid in contact with the bottom boundary. When V≪1 and λ≥1, large-scale circulation has a negligible effect on convection and the heat flux is due to the formation and rise of randomly spaced plumes. When V>10 and λ>100, plume formation is suppressed entirely, and the heat flux is carried by a sheet-like upwelling located in the center of the tank. At intermediate V, and depending on λ, established plume conduits are advected along the bottom boundary and ascending plumes are focused towards the central upwelling. Heat transfer across the layer occurs through a combination of ascending plumes and large-scale flow. Scaling analyses show that the bottom boundary layer thickness and, in turn, the basal heat flux q depend on the Peclet number, Pe, and λ. When λ>10, q∝Pe^1/2 and when λ→1, q∝(Peλ)^1/3, consistent with classical scalings. When applied to the Earth, our results suggest that plate-driven mantle flow focuses ascending plumes towards upwellings in the central Pacific and Africa as well as into mid-ocean ridges. Furthermore, plumes may be captured by strong upwelling flow beneath fast-spreading ridges. This behavior may explain why hotspots are more abundant near slow-spreading ridges than fast-spreading ridges and may also explain some observed variations of mid-ocean ridge basalt (MORB) geochemistry with spreading rate. Moreover, our results suggest that a potentially significant fraction of the core heat flux is due to plumes that are drawn into upwelling flows beneath ridges and not observed as hotspots.     

川西高温水热活动区深部热结构的地球物理分析

川西高温水热活动区是地热学研究的典型地区之一,该区温泉、热泉点分布广泛,这些温泉大多沿北西-南东向的金沙江断裂、德格-乡城断裂、甘孜-理塘断裂、鲜水河三条主断裂呈条带状分布.对此高温水热活动区开展地热学分析,是研究青藏高原"东构造结"深部地球动力学过程、开发利用川西高原地热资源的重要基础,具有重要的科学研究意义和实际应用价值.本文针对实测氦同位素所发现的壳幔热结构比,依据重力、航空磁法、地震资料,采用反演莫霍面、居里面深度等方法,对川西高温水热活动区的深部热结构进行了分析.同时对中高温热泉密集出露的巴塘、理塘以及康定水热活动区的典型壳幔热流构与深部背景场进行了对比研究.结果表明,地表氦同位素热流估算的深部热流结构与地球物理异常结果保持较好的对应关系.川西地区热流结构中地壳、地幔的热贡献比较接近：巴塘地区平均Q_c/Q=51.38%、理塘地区平均Q_c/Q=54.39%、康定地区平均Q_c/Q=42.42%.川西地区表现出了"温壳温幔"型的深部热结构,但温地壳背景下形成较高地表热流的原因有待进一步研究.     

地幔对流的实验研究:非立柱状幔柱和地幔涡旋

地幔对流的物理模拟实验结果表明 ,在地幔介质和温度非均匀分布的复杂条件下 ,热卷流 (地幔柱 )往往由立柱状转变为非立柱状 (含斜柱状、涡旋状等 )。在忽略科里奥利力的情况下 ,板块的下插和滞积下沉、岩石圈根的存在以及地幔介质粘度的非均匀分布等都可能构成不同形状的障碍 -导流体 ,导致地幔的涡旋运动。软流圈中的水平涡旋环带属于对数螺线型 ,环带旋转半径及线速度逐渐减小 ,最终在旋转中心处下沉 ,而旋转角速度大致保持恒定     

Making and Breaking of a Continent: Following the Scent of Geodynamic Imprints on the African Continent Using Electromagnetics

The African continent inherits a long history of continental accretion and breakup. The stage of “making” a continent goes back to the Archean, when the first continental masses formed cratons which mostly remained stable ever since. Subsequent collision of weaker continental masses was followed by several extension and compression episodes that resulted in the formation of super-continents. After the assemblage of Gondwana, a period of predominantly “breaking” , i.e., the breakup of super-continents, took over. The modern-day African continent exhibits different types of margins; continental rifting occurs side by side with recent collision. Since the late 1960s, magnetotelluric (MT) experiments have played an important role in studies of the electrical conductivity structure of Africa. The early results significantly shaped the MT community’s understanding of continental-scale conductivity belts and basic characteristics of cratons and mobile belts on both crustal and lithospheric mantle scales for some decades. Modern MT studies in Africa have generally supported earlier results with high resistivities observed on cratons and low resistivities observed across mobile belts. Advances in instrumentation, data processing and interpretation resulted in higher-resolution images of the lithosphere, which in consequence induce an improved understanding of tectonic processes and geological prerequisites for the occurrence of natural resources. The high electrical conductivity of mobile belts and their relation to reactivated fault and detachment zones were often interpreted to characterize mobile belts as tectonic weak zones, which can accommodate stress and constitute zones along which continents can break. Recent breaking of the African continent can be studied on land across the East African rift; however, the lack of amphibian MT experiments across today’s margins does not allow for good resolution of remnants of continental breakup processes. Naturally, the regions and the focus of the MT studies in Africa are diverse, but they all contribute to the story of making and breaking a continent.     

Mineralogic variability of the uppermost mantle along mid-ocean ridges

Modal analyses of 273 different peridotites representing 43 dredge stations in the Atlantic, Caribbean, and Indian Oceans define three separate melting trends. Peridotites dredged in the vicinity of “mantle plumes” or hot spots have the most depleted compositions in terms of basaltic components, while peridotites dredged at locations removed from such regions are systematically less depleted. The modal data correlate well with mineral compositions, with the peridotites most depleted in pyroxene also having the most refractory mineral compositions. This demonstrates that they are the probable residues of variable degrees of mantle melting. Further, there is a good correlation between the modal compositions of the peridotites and the major element composition of spatially associated dredged basalts. This demonstrates for the first time that the two must be directly related, as is frequently postulated. The high degree of depletion of the peridotites in basaltic major element components in the vicinity of some documented mantle plumes provides direct evidence for a thermal anomaly in such regions—justifying their frequent designation as “hot spots”. The high incompatible element concentrations in these “plume” basalts, however, are contrary to what is expected for such high degrees of melting, and thus require either selective contributions from locally more abundant enriched veins and/or contamination by a volatile-rich metasomatic front from depth.     

On the rise of strongly tilted mantle plume tails

The rise of an initially horizontal, buoyant cylinder of fluid through a denser fluid at low Reynolds number is used to look at the ascent of strongly tilted mantle plumes through the mantle. Such ascents are characterized by (1) the growth of instabilities and (2) the development of a thermal wake downstream. Three-dimensional numerical experiments were carried out to examine these features. An hybrid particle-in-cell finite element method was used to look at the rise of non-diffusing cylinders and, a standard finite element method was used to look at the diffusing case. First the experiments show that the timescale of the fastest growing instability vary with the Rayleigh number and the viscosity ratio. In particular the growth rate decreases as the Rayleigh number decreases, in agreement with our analysis of the laboratory experiments of Kerr et al. (2008). Second the experiments show that the length of the thermal wake increases with the Rayleigh number but the change in viscosity has almost no influence on the wake length. Applied to strongly tilted mantle plumes we conclude that such plumes cannot be unstable given the plume timescales. We also discuss the application of this conclusion to weakly tilted plumes. Besides, this study allows to predict that mantle plumes are unlikely to have developed a significant thermal wake by the time they reach the surface. Finally, the resolution that is required to allow for the growth of mantle plume tails by combined diffusion and thermal entrainment is shown to represent a challenge for the large scale mantle convection simulations.     

Geothermal data analysis at the high-temperature hydrothermal area in Western Sichuan

The western Sichuan hydrothermal area is located at the northeastern margin of the eastern syntaxis of the Qinghai-Tibet Plateau, which is also the eastern end of the Mediterranean-Himalayan geothermal activity zone. There are 248 warm or hot springs in this area, and 11 have temperatures beyond the local boiling temperature. Most of these hot springs are distributed along the Jinshajiang, Dege-Xiangcheng, Ganzi-Litang, and Xianshuihe faults, forming a NW-SE hydrothermal belt. A geothermal analysis of this high-temperature hydrothermal area is an important basis for understanding the deep geodynamic process of the eastern syntaxis of the Qinghai-Tibet Plateau. In addition, this study offers an a priori view to utilize geothermal resources, which is important in both scientific research and application. We use gravity, magnetic, seismic, and helium isotope data to analyze the crust-mantle heat flow ratio and deep geothermal structure. The results show that the background terrestrial heat flow descends from southwest to northeast. The crustal heat ratio is not more than 60%. The high temperature hydrothermal active is related to crustal dynamics processes. Along the Batang-Litang-Kangding line, the Moho depth increases eastward, which is consistent with the changing Qc/Qm(crustal/mantle heat flow) ratio trend. The geoid in the hydrothermal zone is 4–6 km higher than the surroundings, forming a local "platform". The NW-SE striking local tensile stress zone and uplift structure in the upper and middle crust corresponds with the surface hydrothermal active zone. There is an average Curie Point Depth(CPD) of 19.5–22.5 km in Batang, Litang, and Kangding. The local shear-wave(S-wave) velocity is relatively low in the middle and lower crust. The S-wave shows a low velocity trap(Vs3.2 km s.1) at 15–30 km, which is considered a high-temperature partial melting magma, the crustal source of the hydrothermal active zone. We conclude that the hydrothermal system in this area can be divided into Batang-type and Kangding-type, both of which rely on a crustal heating cycle of atmospheric precipitation and surface water along the fracture zone. The heat is derived from the middle and lower crust: groundwater penetrates the deep faults bringing geothermal energy back to the surface and forming high-temperature springs.     

Plume-lithosphere interactions in the ocean basins: constraints from the source mineralogy

Trace element relationships of near-primary alkalic lavas from La Grille volcano, Grande Comore, in the Indian Ocean, as well as those of the Honolulu volcanic series, Oahu, Hawaii, show that their sources contain amphibole and/or phlogopite. Small amounts of each mineral (2% amphibole in the source of La Grille and 0.5% phlogopite plus some amphibole in the source of the Honolulu volcanics) and a range of absolute degrees of partial melting from 1 to 5% for both series are consistent with the observed trace element variation. Amphibole and phlogopite are not stable at the temperatures of convecting upper mantle or upwelling thermal plumes from the deep mantle; however, they are stable at pressure-temperature conditions of the oceanic lithospheric mantle. Therefore, the presence of amphibole and/or phlogopite in the magma source region of volcanics is strong evidence for lithospheric melting, and we conclude that the La Grille and the Honolulu series formed by melting of the oceanic lithospheric mantle.

The identification of amphibole ± phlogopite in the source region of both series implies that the metasomatism by fluids or volatile-rich melts occurred prior to melting. The presence of hydrous phases results in a lower solidus temperature of the lithospheric mantle, which can be reached by conductive heating by the thermal plumes. Isotope ratios of the La Grille and the Honolulu series display a restricted range in composition and represent compositional end-members for each island. Larger isotopic variations in shield lavas, represented by the contemporaneous Karthala volcano on Grande Comore and the older Koolau series on Oahu, reflect interaction of the upwelling thermal plumes with the lithospheric mantle rather than the heterogeneity of deep-seated mantle plume sources or entrainment of mantle material in the rising plume. Literature OsSr isotope ratio covariations constrain the process of plume-lithosphere interaction as occurring through mixing of plume melts and low-degree melts from the metasomatized oceanic lithospheric mantle.

The characterization of the lithospheric mantle signature allows the isotopic composition of the deep mantle plume components to be identified, and mixing relationships show that the Karthala and Koolau plume end-members have nearly uniform isotopic compositions. Based on independent arguments, isotopic variations on Heard and Easter islands have been shown to be a result of mixing between deep plume sources having distinct isotopic compositions with lithosphere or shallow asthenospheric mantle. To the extent that these case studies are representative of oceanic island volcanism, they indicate that interaction with oceanic lithospheric mantle plays an important role in the compositions of lavas erupted during the shield-building stage of plume magmatism, and that isotopic compositions of deep mantle plume sources are nearly uniform on the scale that they are sampled by melting.     

动态地幔柱尾管结构

动态地幔柱模式被广泛用于讨论地球科学中的一些重要课题,如巨大火成区的成因,冈瓦纳古陆解体的原因,板块内中小尺度动力过程的驱动因素等．但是这个基于实验研究而建立的模式中,忽略了地幔柱尾管特征及其作用．地幔柱尾管内温度和速度分布是研究地幔柱上升过程的必要条件．本文从控制尾管结构的基本方程出发,给出了一个定常轴对称地幔柱温度和速度分布的近似分析解．从而得到尾管结构的基本特征：影响尾管内温度分布的主导因素是地幔柱的热流通量,而尾管内上涌速度的大小则不仅取决于热流通量,主要是取决于地幔粘度随深度的变化方式．结果表明,对弱地幔柱,尾管的热损失可能是不可忽略的,而对强地幔柱,径向质量传递可能是不可忽略的.     

A longitudinal seismic reflection profile of the Reykjanes Ridge: Part II — Implications for the mantle hot spot hypothesis

A longitudinal seismic reflection profile of the Reykjanes Ridge, together with earthquake seismicity patterns, is interpreted in terms of the mantle plume hypothesis. Between 52°N and 57°N Reykjanes Ridge is cut by about 12 fractures whose trend, inferred from other data, is approximately east-west. North of 57° there is little or no indication of east-west fracturing.The 57°N transition from fractured to unfractured basement occurs about 900 km southwest of the postulated Iceland mantle plume. The fractured province exhibits higher seismicity and rougher basement, on transverse profiles, than does the unfractured province. A similar transition to rougher, more seismic ridge crest also occurs 900 km northeast of Iceland. We propose that flowage of hot, basalt-rich asthenosphere away from the Iceland hot spot keeps the axial lithosphere hot, thin, sparsely fractured, and relatively aseismic out to 900 km from the plume. Similar effects are evident in the vicinity of some other plumes located near spreading axes. Some plumes also exhibit a greater number of earthquakes at some distance from the spreading axis — possibly a reflection of non-axial igneous activity or fracturing due to local, plume-generated stresses.The regional basement slope along the longitudinal profile is about 8 × 10^?4. If this slope represents a balance between viscous and gravity forces in the flow, a viscosity of the order 10¹⁹ poises can be estimated from the Poiseuille equation.A peculiarly flat, opaque reflector was discovered near the Reykjanes axis, about 300 km southwest of Iceland. Several hypotheses are advanced to account for such reflectors by the exceptional volcanic activity associated with high plume discharge.     

Integrated interpretation of the magnetotelluric and magnetic data from Mahallat geothermal field,Iran

Magnetotelluric (MT) and ground magnetic surveys were conducted on the Mahallat geothermal field situated in Markazi province, central Iran, as a primary part of the explorations and developments of a geothermal energy investigation program in the region. Mahallat region has the greatest geothermal fields in Iran. MT survey was performed in November 2011 on an 8 km profile crossing the hot springs with a total of 17 stations. The 2D inversion of the determinant MT data was performed using a 2D inversion routine based on the Occam approach. The 2D resistivity model obtained from the determinant data shows a low resistivity zone at 800-2000 m depth and a higher resistivity zone above the low resistivity zone, interpreted as geothermal reservoir and cap rock, respectively. It also revealed two major concealed faults which are acting as preferential paths for the circulation of hydrothermal fluids. To obtain more geophysical evidence, a ground magnetic survey with 5000 stations was also performed over an area of 200 km² around the MT profile. Magnetic measurements show a main positive anomaly of about +1000 nT over the study area, which could be interpreted as an intrusive body with the high magnetic susceptibility (i.e. mafic and ultramafic rocks) into the sedimentary host rocks. We interpret the body as the heat source of the geothermal system. Structural index and depth estimation of the anomaly indicate that the intrusive body is similar to a cylinder extending from about one kilometer depth down to greater depths. The results of MT and magnetic investigations indicate a geothermal reservoir which proves the preliminary geological observations to a great extent.     

地球内部压力变化及其构造意义

地球内部的压力远比人们预想的要高得多，地心处的压强为无穷大，在重力分异过程中，地壳地幔减压膨胀（吸热），地核增压收缩（放热），核幔边界是热交换界面，热幔 核幔边界生成。     

Mantle plumes control magnetic reversal frequency

Magnetic reversal frequency correlates inversely with mantle plume activity for the past 150 Ma, as measured by the volume production rate of oceanic plateaus, seamount chains, and continental flood basalts. This inverse correlation is especially striking during the long Cretaceous magnetic normal “superchron”, when mantle plume activity was at a maximum. We suggest that mantle plumes control magnetic reversal frequency by the following sequence of events. Mantle plumes rise from theD″ seismic layer just above the core/mantle boundary, thinningD″ to fuel the plumes. This increases core cooling by allowing heat to be conducted more rapidly across the core/mantle boundary. Outer core convective activity then increases to restore the abnormal heat loss, causing a decrease in magnetic reversal frequency in accord with model predictions for bothα² andαω dynamos. When core convective activity increases above a critical level, a magnetic superchron results. The pulse of plume activity that caused the Cretaceous superchron resulted in a minimum increase in core heat loss of about 1200 GW over the present-day level, which corresponds to an increase in Joule heat production of about 120 GW within the core.     

地幔柱研究进展

地幔柱假设自提出以来,在学术界争议很大,对于其是否存在还没有统一答案.本文对地幔柱的起源做了详细阐述,将质疑者的观点和地幔柱假设自身存在的问题进行罗列,对地幔柱存在的证据进行汇总.从地球化学及岩石学方面对大陆溢流玄武岩和大洋岛玄武岩的分布进行研究,了解地幔柱在地表的形态;从大地构造学角度分析地幔柱如何使超大陆不断合并、...     

Geology, water geochemistry and geothermal potential of the jemez springs area, Canon de San Diego, new Mexico

Studies of the geology, geochemistry of thermal waters, and of one exploratory geothermal well show that two related hot spring systems discharge in Canõn de San Diego at Soda Dam (48°C) and Jemez Springs (72°C). The hot springs discharge from separate strands of the Jemez fault zone which trends northeastward towards the center of Valles Caldera. Exploration drilling to Precambrian basement beneath Jemez Springs encountered a hot aquifer (68°C) at the top of Paleozoic limestone of appropriate temperature and composition to be the local source of the fluids in the surface hot springs at Jemez Springs. Comparisons of the soluble elements Na, Li, Cl, and B, arguments based on isotopic evidence, and chemical geothermometry indicate that the hot spring fluids are derivatives of the deep geothermal fluid within Valles Caldera. No hot aquifer was discovered in or on top of Precambrian basement. It appears that low- to moderate-temperature geothermal reservoirs (< 100°C) of small volume are localized along the Jemez fault zone between Jemez Springs and the margin of Valles Caldera.