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 901 records · Page 50Linked to original sources
USGS2001

Time-variation of hydrothermal discharge at selected sites in the Western United States: Implications for monitoring

We compiled time series of hydrothermal discharge consisting of 3593 chloride- or heat-flux measurements from 24 sites in the Yellowstone region, the northern Oregon Cascades, Lassen Volcanic National Park and vicinity, and Long Valley, California. At all of these sites the hydrothermal phenomena are believed to be as yet unaffected by human activity, though much of the data collection was driven by mandates to collect environmental-baseline data in acticipation of geothermal development. The time series average 19 years in length and some of the Yellowstone sites have been monitored intermittently for over 30 years. Many sites show strong seasonality but few show clear long-term trends, and at most sites statistically significant decadal-scale trends are absent. Thus, the data provide robust estimates of advective heat flow ranging from ~130 MW in the north-central Oregon Cascades to ~6100 MW in the Yellowstone region, and also document Yellowstone hydrothermal chloride and arsenic fluxes of 1740 and 15-20 g/s, respectively. The discharge time series show little sensitivity to regional tectonic events such as earthquakes or inflation/deflation cycles. Most long-term monitoring to date has focused on high-chloride springs and low-temperature fumaroles. The relative stability of these features suggests that discharge measurements done as part of volcano-monitoring programs should focus instead on high-temperature fumaroles, which may be more immediately linked to the magmatic heat source.

Journal of Volcanology and Geothermal Research↗
USGS2003

Summary of recent research in Long Valley Caldera, California

Since 1978, volcanic unrest in the form of earthquakes and ground deformation has persisted in the Long Valley caldera and adjacent parts of the Sierra Nevada. The papers in this special volume focus on periods of accelerated seismicity and deformation in 1980, 1983, 1989-1990, and 1997-1998 to delineate relations between geologic, tectonic, and hydrologic processes. The results distinguish between earthquake sequences that result from relaxation of existing stress accumulation through brittle failure and those in which brittle failure is driven by active intrusion. They also indicate that in addition to a relatively shallow (7-10-km) source beneath the resurgent dome, there exists a deeper (???15-km) source beneath the south moat. Analysis of microgravimety and deformation data indicates that the composition of the shallower source may involve a combination of silicic magma and hydrothermal fluid. Pressure and temperature fluctuations in wells have accompanied periods of crustal unrest, and additional pressure and temperature changes accompanying ongoing geothermal power production have resulted in land subsidence. The completion in 1998 of a 3000-m-deep drill hole on the resurgent dome has provided useful information on present and past periods of circulation of water at temperatures of 100-200??C within the crystalline basement rocks that underlie the post-caldera volcanics. The well is now being converted to a permanent geophysical monitoring station. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2003

Fluid flow in the resurgent dome of Long Valley Caldera: Implications from thermal data and deep electrical sounding

Temperatures of 100??C are measured at 3 km depth in a well located on the resurgent dome in the center of Long Valley Caldera, California, despite an assumed >800??C magma chamber at 6-8 km depth. Local downflow of cold meteoric water as a process for cooling the resurgent dome is ruled out by a Pecle??t-number analysis of temperature logs. These analyses reveal zones with fluid circulation at the upper and lower boundaries of the Bishop Tuff, and an upflow zone in the metasedimentary rocks. Vertical Darcy velocities range from 10 to 70 cm a-1. A 21-km-long geoelectrical profile across the caldera provides resistivity values to the order of 100 to >103 ??m down to a depth of 6 km, as well as variations of self-potential. Interpretation of the electrical data with respect to hydrothermal fluid movement confirms that there is no downflow beneath the resurgent dome. To explain the unexpectedly low temperatures in the resurgent dome, we challenge the common view that the caldera as a whole is a regime of high temperatures and the resurgent dome is a local cold anomaly. Instead, we suggest that the caldera was cooled to normal thermal conditions by vigorous hydrothermal activity in the past, and that a present-day hot water flow system is responsible for local hot anomalies, such as Hot Creek and the area of the Casa Diablo geothermal power plant. The source of hot water has been associated with recent shallow intrusions into the West Moat. The focus of planning for future power plants should be to locate this present-day flow system instead of relying on heat from the old magma chamber. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS1978

Applications of the VLF induction method for studying some volcanic processes of Kilauea volcano, Hawaii

The very low-frequency (VLF) induction method has found exceptional utility in studying various volcanic processes of Kilauea volcano, Hawaii because: (1) significant anomalies result exclusively from ionically conductive magma or still-hot intrusions (> 800??C) and the attendant electrolytically conductive hot groundwater; (2) basalt flows forming the bulk of Kilauea have very high resistivities at shallow depths that result in low geologic noise levels and relatively deep depths of investigation (???100 m); and (3) the azimuths to two of the usable transmitters (NLK and NPM) are aligned favorably with most of the principal geologic features. Measurements of the tilt angle and ellipticity of the polarization ellipse of the magnetic field, using a simple, hand-held receiver, have been used to: (1) delineate the lateral extent of shallow, partially solidified lava lakes, active lava tubes, and recent intrusive dikes; (2) obtain an indication of the attitude of some recent dikes; (3) show that many eruptive fissures cool faster than their intrusive counterparts; (4) show that some fumarolic areas are underlain by shallow, highly altered, and conductive zones; and (5) provide control information for interpreting data obtained with other electrical techniques. Complementary measurements of scalar apparent resistivity and surface impedance phase, using a new attachment for the VLF receiver, have substantially increased the utility of VLF studies in Kilauea. They provide better lateral resolution of conductors and reduce the ambiguity in interpretation. Notwithstanding recent advances in theoretical modeling techniques, the excellent quality of some of the data warrants extension of interpretive techniques, particularly for quantitatively characterizing the configuration and conductivity of small-dimension bodies. These VLF induction methods should have wide application to studies of active volcanic regions in other parts of the world and could provide some insights into the workings of larger-scaled geothermal systems. ?? 1978.

Journal of Volcanology and Geothermal Research↗
USGS1987

Concentrations of arsenic, antimony, and boron in steam and steam condensate at The Geysers, California

Studies at The Geysers Geothermal Field, California indicate that under some circumstances elements that are transported in the vapor phase can become enriched in the liquid phase. Waters from two condensate traps (steam traps) on steam lines at The Geysers are enriched with arsenic, antimony, and boron compared to the concentrations of these elements in coexisting steam. Concentrations of boron in condensate-trap waters were as high as 160 mg/L, arsenic as high as 35 mg/L, and antimony as high as 200 ??g/L. Enrichment of arsenic, antimony, and boron is at least partially controlled by the partitioning of these elements into the liquid phase, according to their vapor-liquid distribution coefficients, after they are transported in steam. Several of the elements that are most soluble in steam, including arsenic and antimony, are part of the trace-element suite that characterizes precious-metal epithermal ore deposits. ?? 1987.

Journal of Volcanology and Geothermal Research↗
USGS1987

The hydrothermal system of the Calabozos caldera, central Chilean Andes

Active thermal springs associated with the late Pleistocene Calabozos caldera complex occur in two groups: the Colorado group which issues along structures related to caldera collapse and resurgence, and the Puesto Calabozos group, a nearby cluster that is chemically distinct and probably unrelated to the Colorado springs. Most of the Colorado group can be related to a hypothetical parent water containing ???400 ppm Cl at ???250??C by dilution with ???50% of cold meteoric water. The thermal springs in the most deeply eroded part of the caldera were derived from the same parent water by boiling. The hydrothermal system has probably been active for at least as long as 300,000 years, based on geologic evidence and calculations of paleo-heat flow. There is no evidence for economic mineralization at shallow depth. The Calabozos hydrothermal system would be an attractive geothermal prospect were its location not so remote. ?? 1987.

Journal of Volcanology and Geothermal Research↗
USGS1989

The evolution of thermal structure and water chemistry in Lake Nyos

We collected a time series of physical and chemical data to gain a better understanding of the dynamics of Lake Nyos. Measurements of water and gas chemistry, and temperature made during January, March, and May 1987 are compared to data taken in September 1986 just after the initial CO2 gas release. There is no pattern of change in overall heat content of the lake, although heat input to bottom waters (185-208 m) has occurred at a rate of 1600 mW m-2. This increase in heat content translates to a change from 23.38 to 24.12??C at 200 m and can be explained by geothermal heat flow and addition of thermal spring water. Concentrations of Ca2+, Mg2+, Na+, K+, Fe2+ and alkalinity have increased only in bottom waters. In situ lake processes such as sulfate and iron reduction are unable to account for the changes in alkalinity. Observed chemical changes are consistent with a scenario where slightly thermal soda water is being input to the bottom of the lake. Measurements of pCO2 at depth ranged from 18 to 28% of saturation and exhibited horizontal variability. Overall recharge of CO2 in bottom waters is negligible. Mainly because of increasing ion concentrations in bottom water, total stability of the water column increased 33% from 48,800 J m-2 in September 1986 to 64,700 J m-2 in May 1987. As long as CO2 concentrations remain the same, this level of stability is higher than could be disrupted by common limnologic or meteorologic processes. There is thermal and chemical evidence that a buildup of dissolved iron and CO2 in bottom waters must have preceded the August 1986 gas release. In addition, a survey of all crater lakes in Cameroon indicates that only Lakes Nyos and Monoun contain high concentrations of dissolved iron and CO2. Thus there is a low probability of any other Cameroonian lake releasing a substantial volume of CO2. ?? 1989.

Journal of Volcanology and Geothermal Research↗
USGS1990

Deformation monitoring at Nevado del Ruiz, Colombia - October 1985 - March 1988

Deformation studies began at Nevado del Ruiz 23 days before the devastating 13 November 1985 eruption, at least 12 months after precursory seismicity and fumarolic activity began. The late start in geodetic monitoring, limited number of stations in the pre-eruption network, and inconsistent patterns in the observed deformation limit conclusions about intrusive activity in the months and weeks prior to the eruption. However, the data require that the magma source of the devastating eruption was either deeper than 7 km or, if shallow, recovered the same volume and position within one week of the eruption. Geodetic monitoring resumed 1 week after the eruption and, by April 1986, included 11 tilt-leveling stations, 38 EDM lines, and 7 short leveling lines - a network capable of detecting emplacement or movement of magma volumes as small as 3 MCM (3 ?? 106 m3) to a depth of 2-3 km (using a point-source model), 10 MCM to 7 km, 50 MCM to 10 km, and 200 MCM to 15 km beneath Ruiz. In addition, 4 telemetered tiltmeters provided the capability of detecting, in real time, the fairly rapid ascent of much smaller magma bodies. Stations established to detect instability of the summit ice cap after the eruption were discontinued in early 1986. The data collected from the geodetic networks have higher than normal variance but demonstrate that little or no cumulative deformation of Ruiz occurred from October 1985 through March 1988. Thus, little, if any, magma intruded above 5 km beneath the summit during or after the 13 November 1985 eruption. This lack of significant intrusive activity agrees with the surprisingly low seismic energy release under Ruiz and makes direct degassing of a large batholith an improbable explanation of the large sulfur flux to date at Ruiz. Part of the variance in the geodetic data results from real but noncumulative deformation that may in part be pressure-buffered by a fairly large geothermal water-gas mixture for which abundant physical evidence exists. Part of the noncumulative deformation, some of the fairly dispersed and low-level seismicity under Ruiz, and some phreatic events appear to correlate with seasonal precipitation patterns. Hence rain/snow-loading and groundwater interaction may cause deformation events and possibly help trigger some phreatic explosions and seismic events at Ruiz and, as search of the literature reveals, at other volcanoes in metastable states. ?? 1990.

Journal of Volcanology and Geothermal Research↗
USGS1992

Thermal waters along the Konocti Bay fault zone, Lake County, California: a re-evaluation

The Konocti Bay fault zone (KBFZ), initially regarded by some as a promising target for liquid-dominated geothermal systems, has been a disappointment. At least five exploratory wells were drilled in the vicinity of the KBFZ, but none were successful. Although the Na-K-Ca and Na-Li geothermometers indicate that the thermal waters discharging in the vicinity of Howard and Seigler Springs may have equilibrated at temperatures greater than 200??C, the spring temperatures and fluid discharges are low. Most thermal waters along the KBFZ contain >100 mg/l Mg. High concentrations of dissolved magnesium are usually indicative of relatively cool hydrothermal systems. Dissolution of serpentine at shallow depths may contribute dissolved silica and magnesium to rising thermal waters. Most thermal waters are saturated with respect to amorphous silica at the measured spring temperature. Silica geothermometers and mixing models are useless because the dissolved silica concentration is not controlled by the solubility of either quartz or chalcedony. Cation geothermometry indicates the possibility of a high-temperature fluid (> 200??C) only in the vicinity of Howard and Seigler Springs. However, even if the fluid temperature is as high as that indicated by the geothermometers, the permeability may be low. Deuterium and oxygen-18 values of the thermal waters indicate that they recharged locally and became enriched in oxygen-18 by exchange with rock. Diluting meteoric water and the thermal water appear to have the same deuterium value. Lack of tritium in the diluted spring waters suggest that the diluting water is old. ?? 1992.

Journal of Volcanology and Geothermal Research↗
USGS1993

Relations of ammonium minerals at several hydrothermal systems in the western U.S.

Ammonium bound to silicate and sulfate minerals has recently been located at several major hydrothermal systems in the western U.S. utilizing newly-discovered near-infrared spectral properties. Knowledge of the origin and mineralogic relations of ammonium minerals at known hydrothermal systems is critical for the proper interpretation of remote sensing data and for testing of possible links to mineralization. Submicroscopic analysis of ammonium minerals from two mercury- and gold-bearing hot-springs deposits at Ivanhoe, Nevada and McLaughlin, California shows that the ammonium feldspar, buddingtonite, occurs as fine-grained euhedral crystals coating larger sulfide and quartz crystals. Ammonium feldspar seems to precipitate relatively late in the crystallization sequence and shows evidence for replacement of NH4+ by K+ or other monovalent cations. Some buddingtonite is observed in close association with mercury, but not with gold. Ammonioalunite is found in a variety of isolated crystal forms at both deposits. Nitrogen isotopic values for ammonium-bearing minerals show a 14??? range in composition, precluding assignment of a specific provenance to the nitrogen. The correlations of nitrogen isotopic values with depth and ammonium content suggest some loss of nitrogen in the oxidizing supergene environment, possibly as a metastable mineral. The high ammonium content in these hydrothermal systems, the close association to mercury, and the small crystal size of the ammonium-bearing minerals all suggest that ammonium may be transported in a late-stage vapor phase or as an organic volatile. Such a process could lead to the formation of a non-carbonaceous organic aureole above a buried geothermal source. The discovery of a 10-km outcrop of ammonium minerals confirms that significant substitution of ammonium in minerals is possible over an extensive area and that remote sensing is a feasible means to detect such aureoles. ?? 1993.

Journal of Volcanology and Geothermal Research↗
USGS2018

New data yield new geologic insights at the Fallon FORGE site, Carson Sink Region, Nevada

The geologic structure beneath the Fallon Frontier Observatory for Research in Geothermal Energy (FORGE) site represents a record of the Mesozoic through Cenozoic tectonism, volcanism, and sedimentation that has affected the Carson Sink local to Fallon, NV. A robust dataset confirms that the lithologic sequence consists of Quaternary through Miocene sedimentary and volcanic rocks resting non-conformably on Mesozoic crystalline basement. The basement consists of four lithologic units; 1) Triassic-Jurassic low-to-medium grade meta-rhyolites, 2) Jurassic low-to-medium grade quartzites and other metamorphosed marine sedimentary rocks, and 3) Jurassic low-to-medium grade meta-basalt and -basaltic-andesite lavas, all intruded by 4) Jurassic-Cretaceous quartz monzonite. The geologic section dips ~20-25° west, tilting that was accommodated by a predominant system of north-to-north-northeast striking, east and west moderately-to-steeply dipping normal faults. The above, relatively broad-scale characteristics of the geologic framework of the Fallon site have been developed throughout the ~3 year duration of the project. With the collection of new data and re-analysis of existing data in Phase 2B of the Fallon FORGE project, our detailed understanding of the relatively finer-scale aspects of the geologic framework, including aspects of the stratigraphic sequence and the locations and attitudes of individual faults have evolved and an updated 3D geologic map has been developed. Here, we compare the Phase 1 3D geologic map to the Phase 2B 3D geologic map and demonstrate the evolution of our understanding of the geologic framework of the Fallon site and the value of the new data that was collected in Phase 2B in developing this updated framework.

Geothermal Resources Council Transactions↗
USGS2023

Predicting large hydrothermal systems

We train five models using two machine learning (ML) regression algorithms (i.e., linear regression and XGBoost) to predict hydrothermal upflow in the Great Basin. Feature data are extracted from datasets supporting the INnovative Geothermal Exploration through Novel Investigations Of Undiscovered Systems project (INGENIOUS). The label data (the reported convective signals) are extracted from measured thermal gradients in wells by comparing the total estimated heat flow at the wells to the modeled background conductive heat flow. That is, the reported convective signal is the difference between the background conductive heat flow and the well heat flow. The reported convective signals contain outliers that may affect upflow prediction, so the influence of outliers is tested by constructing models for two cases: 1) using all the data (i.e., -91 to 11,105 mW/m2), and 2) truncating the range of labels to include only reported convective signals between -25 and 200 mW/m2. Because hydrothermal systems are sparse, models that predict high convective signal in smaller areas better match the natural frequency of hydrothermal systems. Early results demonstrate that XGBoost outperforms linear regression. For XGBoost using the truncated range of labels, half of the high reported signals are within < 3 % of the highest predictions. For XGBoost using the entire range of labels, half of the high reported signals are in < 13 % of the highest predictions. While this implies that the truncated regression is superior, the all-data model better predicts the locations of power-producing systems (i.e., the operating power plants are in a smaller fraction of the study area given by the highest predictions). Even though the models generally predict greater hydrothermal upflow for higher reported convective signals than for lower reported convective signals, both XGBoost models consistently underpredict the magnitude of higher signals. This behavior is attributed to low resolution/granularity of input features compared with the scale of a hydrothermal upflow zone (a few km or less across). Trouble estimating exact values while still reliably predicting high versus low convective signals suggests that a future strategy such as ranked ordinal regression (e.g., classifying into ordered bins for low, medium, high, and very high convective signal) might fit better models, since doing so reduces problems introduced by outliers while preserving the property of larger versus smaller signals.

Geothermal Resources Council Transactions↗
USGS2014

Maximum magnitude earthquakes induced by fluid injection

Analysis of numerous case histories of earthquake sequences induced by fluid injection at depth reveals that the maximum magnitude appears to be limited according to the total volume of fluid injected. Similarly, the maximum seismic moment seems to have an upper bound proportional to the total volume of injected fluid. Activities involving fluid injection include (1) hydraulic fracturing of shale formations or coal seams to extract gas and oil, (2) disposal of wastewater from these gas and oil activities by injection into deep aquifers, and (3) the development of enhanced geothermal systems by injecting water into hot, low-permeability rock. Of these three operations, wastewater disposal is observed to be associated with the largest earthquakes, with maximum magnitudes sometimes exceeding 5. To estimate the maximum earthquake that could be induced by a given fluid injection project, the rock mass is assumed to be fully saturated, brittle, to respond to injection with a sequence of earthquakes localized to the region weakened by the pore pressure increase of the injection operation and to have a Gutenberg-Richter magnitude distribution with a b value of 1. If these assumptions correctly describe the circumstances of the largest earthquake, then the maximum seismic moment is limited to the volume of injected liquid times the modulus of rigidity. Observations from the available case histories of earthquakes induced by fluid injection are consistent with this bound on seismic moment. In view of the uncertainties in this analysis, however, this should not be regarded as an absolute physical limit.

Journal of Geophysical Research B: Solid Earth↗
USGS2014

Soil CO2 emissions as a proxy for heat and mass flow assessment, Taupō Volcanic Zone, New Zealand

The quantification of heat and mass flow between deep reservoirs and the surface is important for understanding magmatic and hydrothermal systems. Here, we use high-resolution measurement of carbon dioxide flux (φCO 2 ) and heat flow at the surface to characterize the mass (CO 2 and steam) and heat released to the atmosphere from two magma-hydrothermal systems. Our soil gas and heat flow surveys at Rotokawa and White Island in the Taupō Volcanic Zone, New Zealand, include over 3000 direct measurements of φCO 2 and soil temperature and 60 carbon isotopic values on soil gases. Carbon dioxide flux was separated into background and magmatic/hydrothermal populations based on the measured values and isotopic characterization. Total CO 2 emission rates (ΣCO 2 ) of 441 ± 84 t d −1 and 124 ± 18 t d −1 were calculated for Rotokawa (2.9 km 2 ) and for the crater floor at White Island (0.3 km 2 ), respectively. The total CO 2 emissions differ from previously published values by +386 t d −1 at Rotokawa and +25 t d −1 at White Island, demonstrating that earlier research underestimated emissions by 700% (Rotokawa) and 25% (White Island). These differences suggest that soil CO 2 emissions facilitate more robust estimates of the thermal energy and mass flux in geothermal systems than traditional approaches. Combining the magmatic/hydrothermal-sourced CO 2 emission (constrained using stable isotopes) with reservoir H 2 O:CO 2 mass ratios and the enthalpy of evaporation, the surface expression of thermal energy release for the Rotokawa hydrothermal system (226 MW t ) is 10 times greater than the White Island crater floor (22.5 MW t ).

Taupō Volcanic Zone↗
USGS2016

Seismic imaging of the metamorphism of young sediment into new crystalline crust in the actively rifting Imperial Valley, California

Plate-boundary rifting between transform faults is opening the Imperial Valley of southern California and the rift is rapidly filling with sediment from the Colorado River. Three 65–90 km long seismic refraction profiles across and along the valley, acquired as part of the 2011 Salton Seismic Imaging Project, were analyzed to constrain upper crustal structure and the transition from sediment to underlying crystalline rock. Both first arrival travel-time tomography and frequency-domain full-waveform inversion were applied to provide P-wave velocity models down to ∼7 km depth. The valley margins are fault-bounded, beyond which thinner sediment has been deposited on preexisting crystalline rocks. Within the central basin, seismic velocity increases continuously from ∼1.8 km/s sediment at the surface to >6 km/s crystalline rock with no sharp discontinuity. Borehole data show young sediment is progressively metamorphosed into crystalline rock. The seismic velocity gradient with depth decreases approximately at the 4 km/s contour, which coincides with changes in the porosity and density gradient in borehole core samples. This change occurs at ∼3 km depth in most of the valley, but at only ∼1.5 km depth in the Salton Sea geothermal field. We interpret progressive metamorphism caused by high heat flow to be creating new crystalline crust throughout the valley at a rate comparable to the ≥2 km/Myr sedimentation rate. The newly formed crystalline crust extends to at least 7–8 km depth, and it is shallower and faster where heat flow is higher. Most of the active seismicity occurs within this new crust.

California↗
USGS2018

Mainshock-aftershock clustering in volcanic regions

Earthquakes break their general Poissonean behavior through two types of seismic bursts: swarms and mainshock-aftershock sequences. The former is commonly thought to dominate in volcanic and geothermal regions, but aftershock production, including within swarms, is not well studied in volcanic regions. Here we compare mainshock-aftershock clustering in active volcanic regions in Japan to nearby nonvolcanic regions. We find that aftershock production is similar in both areas by two separate metrics: (1) Both volcanic and nonvolcanic regions have similar proportions of areas that cluster into mainshock-aftershock sequences. (2) Volcanic areas with mainshock-aftershock sequences have aftershock productivity at least as high as nonvolcanic regions. We also find that volcano-tectonic events that are precursors to an eruption are more common at volcanoes without mainshock-aftershock clusters than at volcanoes with well-defined mainshock-aftershock clusters. This last finding hints at a strategy to identify volcanic systems where cataloged earthquakes are good predictors of behavior.

Geophysical Research Letters↗
USGS2017

Quantifying the heterogeneity of the tectonic stress field using borehole data

The heterogeneity of the tectonic stress field is a fundamental property which influences earthquake dynamics and subsurface engineering. Self-similar scaling of stress heterogeneities is frequently assumed to explain characteristics of earthquakes such as the magnitude-frequency relation. However, observational evidence for such scaling of the stress field heterogeneity is scarce. We analyze the local stress orientations using image logs of two closely spaced boreholes in the Coso Geothermal Field with sub-vertical and deviated trajectories, respectively, each spanning about 2 km in depth. Both the mean and the standard deviation of stress orientation indicators (borehole breakouts, drilling-induced fractures and petal-centerline fractures) determined from each borehole agree to the limit of the resolution of our method although measurements at specific depths may not. We find that the standard deviation in these boreholes strongly depends on the interval length analyzed, generally increasing up to a wellbore log length of about 600 m and constant for longer intervals. We find the same behavior in global data from the World Stress Map. This suggests that the standard deviation of stress indicators characterizes the heterogeneity of the tectonic stress field rather than the quality of the stress measurement. A large standard deviation of a stress measurement might be an expression of strong crustal heterogeneity rather than of an unreliable stress determination. Robust characterization of stress heterogeneity requires logs that sample stress indicators along a representative sample volume of at least 1 km.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Occurrence of arsenic in ground water of the Middle Rio Grande Basin, central New Mexico

Chemical data from more than 400 ground-water sites in the Middle Rio Grande Basin of central New Mexico indicate that arsenic concentrations exceed the U.S. Environmental Protection Agency drinking-water standard of 10 micrograms per liter across broad areas of the Santa Fe Group aquifer system, which is currently the almost exclusive source of drinking-water supply for residents of the basin. Identification of sources of arsenic to ground water of the basin is complicated by multiple sources of ground-water recharge that differ substantially in chemical composition. Establishment of a clear hydrologic framework for the basin was useful in interpreting the significance of patterns in arsenic concentration. This investigation indicates that there are two main sources of high-arsenic water to the Middle Rio Grande Basin. One primary source is related to silicic volcanism in the Jemez Mountains to the north, where dilute recharge water likely flows through rocks that have been altered by contact with geothermal fluids. The other primary source is mineralized water of deep origin that mixes with shallower ground water in several locations around the basin, particularly along major structural features. Ground water that has not been affected by either of these two high-arsenic sources generally has low arsenic concentrations. In some areas of the basin, values of pH exceeding about 8.5 appear to contribute to elevated arsenic concentrations through desorption of arsenic from metal oxides.

New Mexico↗