USGS Science⌕ Search

SEARCH · USGS Science

Results for “The R Journal”

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 595 records · Page 33Linked to original sources
USGS2016

Automated mapping of persistent ice and snow cover across the western U.S. with Landsat

We implemented an automated approach for mapping persistent ice and snow cover (PISC) across the conterminous western U.S. using all available Landsat TM and ETM+ scenes acquired during the late summer/early fall period between 2010 and 2014. Two separate validation approaches indicate this dataset provides a more accurate representation of glacial ice and perennial snow cover for the region than either the U.S. glacier database derived from US Geological Survey (USGS) Digital Raster Graphics (DRG) maps (based on aerial photography primarily from the 1960s–1980s) or the National Land Cover Database 2011 perennial ice and snow cover class. Our 2010–2014 Landsat-derived dataset indicates 28% less glacier and perennial snow cover than the USGS DRG dataset. There are larger differences between the datasets in some regions, such as the Rocky Mountains of Northwest Wyoming and Southwest Montana, where the Landsat dataset indicates 54% less PISC area. Analysis of Landsat scenes from 1987–1988 and 2008–2010 for three regions using a more conventional, semi-automated approach indicates substantial decreases in glaciers and perennial snow cover that correlate with differences between PISC mapped by the USGS DRG dataset and the automated Landsat-derived dataset. This suggests that most of the differences in PISC between the USGS DRG and the Landsat-derived dataset can be attributed to decreases in PISC, as opposed to differences between mapping techniques. While the dataset produced by the automated Landsat mapping approach is not designed to serve as a conventional glacier inventory that provides glacier outlines and attribute information, it allows for an updated estimate of PISC for the conterminous U.S. as well as for smaller regions. Additionally, the new dataset highlights areas where decreases in PISC have been most significant over the past 25–50 years.

ISPRS Journal of Photogrammetry and Remote Sensing↗
USGS2013

Next generation of global land cover characterization, mapping, and monitoring

Land cover change is increasingly affecting the biophysics, biogeochemistry, and biogeography of the Earth's surface and the atmosphere, with far-reaching consequences to human well-being. However, our scientific understanding of the distribution and dynamics of land cover and land cover change (LCLCC) is limited. Previous global land cover assessments performed using coarse spatial resolution (300 m–1 km) satellite data did not provide enough thematic detail or change information for global change studies and for resource management. High resolution (∼30 m) land cover characterization and monitoring is needed that permits detection of land change at the scale of most human activity and offers the increased flexibility of environmental model parameterization needed for global change studies. However, there are a number of challenges to overcome before producing such data sets including unavailability of consistent global coverage of satellite data, sheer volume of data, unavailability of timely and accurate training and validation data, difficulties in preparing image mosaics, and high performance computing requirements. Integration of remote sensing and information technology is needed for process automation and high-performance computing needs. Recent developments in these areas have created an opportunity for operational high resolution land cover mapping, and monitoring of the world. Here, we report and discuss these advancements and opportunities in producing the next generations of global land cover characterization, mapping, and monitoring at 30-m spatial resolution primarily in the context of United States, Group on Earth Observations Global 30 m land cover initiative (UGLC).

International Journal of Applied Earth Observation↗
USGS2006

Quantitative estimation of minimum offset for multichannel surface-wave survey with actively exciting source

Multichannel analysis of surface waves is a developing method widely used in shallow subsurface investigations. The field procedures and related parameters are very important for successful applications. Among these parameters, the source-receiver offset range is seldom discussed in theory and normally determined by empirical or semi-quantitative methods in current practice. This paper discusses the problem from a theoretical perspective. A formula for quantitatively evaluating a layered homogenous elastic model was developed. The analytical results based on simple models and experimental data demonstrate that the formula is correct for surface wave surveys for near-surface applications. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Applied Geophysics↗
USGS2008

Biogeochemical and ecological impacts of livestock grazing in semi-arid southeastern Utah, USA

Relatively few studies have examined the ecological and biogeochemical effects of livestock grazing in southeastern Utah. In this study, we evaluated how grazing has affected soil organic carbon and nitrogen to a depth of 50 cm in grasslands located in relict and actively-grazed sites in the Canyonlands physiographic section of the Colorado Plateau. We also evaluated differences in plant ground cover and the spatial distribution of soil resources. Results show that areas used by domestic livestock have 20% less plant cover and 100% less soil organic carbon and nitrogen compared to relict sites browsed by native ungulates. In actively grazed sites, domestic livestock grazing also appears to lead to clustered, rather than random, spatial distribution of soil resources. Magnetic susceptibility, a proxy for soil stability in this region, suggests that grazing increases soil erosion leading to an increase in the area of nutrient-depleted bare ground. Overall, these results, combined with previous studies in the region, suggest that livestock grazing affects both plant cover and soil fertility with potential long-term implications for the sustainability of grazing operations in this semi-arid landscape. ?? 2007 Elsevier Ltd. All rights reserved.

Journal of Arid Environments↗
USGS2009

Using packrat middens to assess grazing effects on vegetation change

Research on grazing effects usually compares the same sites through time or grazed and ungrazed sites over the same time period. Both approaches are complicated in arid environments where grazing can have a long undocumented history and landscapes can be spatially heterogenous. This work employs both approaches simultaneously by comparing grazed and ungrazed samples through both time and space using fossil plant macrofossils and pollen from packrat middens. A series of 27 middens, spanning from 995 yr BP to the present, were collected from Glen Canyon in southeastern Utah, USA. These middens detail vegetation change just prior to, and following, the historical introduction of domesticated grazers and also compares assemblages from nearby ungrazable mesas. Pre-grazing middens, and modern middens from ungrazed areas, record more native grasses, native herbs, and native shrubs such as Rhus trilobata, Amelanchier utahensis, and Shepherdia rotundifolia than modern middens from grazed areas. Ordinations demonstrate that site-to-site variability is more important than any temporal changes, making selection of comparable grazed versus ungrazed study treatments difficult. But within similar sites, the changes through time show that grazing lowered the number of taxa recorded, and lessened the pre-existing site differences, homogenizing the resultant plant associations in this desert grassland.

Journal of Arid Environments↗
USGS2009

Microbial characterization of nitrification in a shallow, nitrogen-contaminated aquifer, Cape Cod, Massachusetts and detection of a novel cluster associated with nitrifying Betaproteobacteria

Groundwater nitrification is a poorly characterized process affecting the speciation and transport of nitrogen. Cores from two sites in a plume of contamination were examined using culture-based and molecular techniques targeting nitrification processes. The first site, located beneath a sewage effluent infiltration bed, received treated effluent containing O 2 (> 300 µM) and NH 4 + (51–800 µM). The second site was 2.5 km down-gradient near the leading edge of the ammonium zone within the contaminant plume and featured vertical gradients of O 2 , NH 4 + , and NO 3 − (0–300, 0–500, and 100–200 µM with depth, respectively). Ammonia- and nitrite-oxidizers enumerated by the culture-based MPN method were low in abundance at both sites (1.8 to 350 g − 1 and 33 to 35,000 g − 1 , respectively). Potential nitrifying activity measured in core material in the laboratory was also very low, requiring several weeks for products to accumulate. Molecular analysis of aquifer DNA (nested PCR followed by cloning and 16S rDNA sequencing) detected primarily sequences associated with the Nitrosospira genus throughout the cores at the down-gradient site and a smaller proportion from the Nitrosomonas genus in the deeper anoxic, NH 4 + zone at the down-gradient site. Only a single Nitrosospira sequence was detected beneath the infiltration bed. Furthermore, the majority of Nitrosospira -associated sequences represent an unrecognized cluster. We conclude that an uncharacterized group associated with Nitrosospira dominate at the geochemically stable, down-gradient site, but found little evidence for Betaproteobacteria nitrifiers beneath the infiltration beds where geochemical conditions were more variable.

Journal of Contaminant Hydrology↗
USGS2010

The Lake Ontario zooplankton community before (1987-1991) and after (2001-2005) invasion-induced ecosystem change

We assessed changes in Lake Ontario zooplankton biomass, production, and community composition before (1987–1991) and after (2001–2005) invasion-induced ecosystem changes. The ecosystem changes were associated with establishment of invasive dreissenid mussels and invasive predatory cladocerans ( Bythotrephes and Cercopagis ). Whole-lake total epilimnetic plus metalimnetic zooplankton production declined by approximately half from 42.45 (g dry wt∙m −2 ∙ year −1 ) during 1987–1991 to 21.91 (g dry wt∙m −2 ∙ year −1 ) in 2003 and averaged 21.01 (g dry wt∙m −2 ∙ year −1 ) during 2001–2005. Analysis of two independent data sets indicates that the mean biomass and biomass proportion of cyclopoid copepods declined while the same measures increased for the invasive predatory cladocerans. Changes in means and proportions of all other zooplankton groups were not consistent between the data sets. Cyclopoid copepod biomass and production declined by factors ranging from 3.6 to 5.7. Invasive predatory cladoceran biomass averaged from 5.0% to 8.0% of the total zooplankton biomass. The zooplankton community was otherwise resilient to the invasion-induced disruption as zooplankton species richness and diversity were unaffected. Zooplankton production was likely reduced by declines in primary productivity but may have declined further due to increased predation by alewives and invasive predatory cladocerans. Shifts in zooplankton community structure were consistent with increased predation pressure on cyclopoid copepods by alewives and invasive predatory cladocerans. Predicted declines in the proportion of small cladocerans were not evident. This study represents the first direct comparison of changes in Lake Ontario zooplankton production before and after the invasion-induced disruption and will be important to food web-scale investigations of invasion effects.

Journal of Great Lakes Research↗
USGS2019

Tracking changes in nutrient delivery to western Lake Erie: Approaches to compensate for variability and trends in streamflow

Tracking changes in stream nutrient inputs to Lake Erie over multidecadal time scales depends on the use of statistical methods that can remove the influence of year-to-year variability of streamflow but also explicitly consider the influence of long-term trends in streamflow. The methods introduced in this paper include an extended version of Weighted Regressions on Time, Discharge, and Season (WRTDS) modeling that explicitly considers nonstationary streamflow by incorporating information on changes in the frequency distribution of daily measured streamflow (discharge) over time. Soluble reactive phosphorus (SRP) trends in annual flow-normalized fluxes (loads) at five long-term monitoring sites in the western Lake Erie drainage basin show increases of 109 to 322% over the period 1995 to 2015. About one-third of the increase appears attributable to increasing discharge trends, while the remaining two-thirds appears to be driven by changes in concentration versus discharge relationships reflecting higher concentrations for any given discharge during recent years. Trends in total phosphorus and three nitrogen parameters (total nitrogen, nitrate-nitrite, and total Kjeldahl nitrogen) at the 10 sites analyzed were much less pronounced, and commonly show decreases in concentration-discharge relationships accompanied by increases in discharge, resulting in little net change in total flux. Trends in monthly SRP fluxes and discharge, dissolved versus particulate fractions of nutrients, and N:P flux ratios were also evaluated. The methods described here provide tools to more clearly discern the effectiveness of nutrient-control strategies and can serve as ongoing measures of progress, or lack of progress, towards nutrient-reduction goals.

Journal of Great Lakes Research↗
USGS2005

Hydrologic properties of coal beds in the Powder River Basin, Montana I. Geophysical log analysis

As part of a multidisciplinary investigation designed to assess the implications of coal-bed methane development on water resources for the Powder River Basin of southeastern Montana, six wells were drilled through Paleocene-age coal bedsalong a 31-km east–west transect within the Tongue River drainage basin. Analysis of geophysical logs obtained in these wells provides insight into the hydrostratigraphic characteristics of the coal and interbedded siliciclastic rocks and their possible interaction with the local stress field. Natural gamma and electrical resistivity logs were effective in distinguishing individual coal beds. Full-waveform sonic logs were used to determine elastic properties of the coal and an attendant estimate of aquifer storage is in reasonable agreement with that computed from a pumping test. Inspection of magnetically oriented images of the borehole walls generated from both acoustic and optical televiewers and comparison with coal cores infer a face cleat orientation of approximately N33°E, in close agreement with regional lineament patterns and the northeast trend of the nearby Tongue River. The local tectonic stress field in this physiographic province as inferred from a nearby 1984 earthquake denotes an oblique strike-slip faulting regime with dominant east–west compression and north–south extension. These stress directions are coincident with those of the primary fracture sets identified from the televiewer logs and also with the principle axes of the drawdown ellipse produced from a complementary aquifer test, but oblique to apparent cleat orientation. Consequently, examination of these geophysical logs within the context of local hydrologic characteristics indicates that transverse transmissivity anisotropy in these coals is predominantly controlled by bedding configuration and perhaps a mechanical response to the contemporary stress field rather than solely by cleat structure.

Montana↗
USGS2013

Trends in the suspended-sediment yields of coastal rivers of northern California, 1955–2010

Time-dependencies of suspended-sediment discharge from six coastal watersheds of northern California – Smith River, Klamath River, Trinity River, Redwood Creek, Mad River, and Eel River – were evaluated using monitoring data from 1955 to 2010. Suspended-sediment concentrations revealed time-dependent hysteresis and multi-year trends. The multi-year trends had two primary patterns relative to river discharge: (i) increases in concentration resulting from both land clearing from logging and the flood of record during December 1964 (water year 1965), and (ii) continual decreases in concentration during the decades following this flood. Data from the Eel River revealed that changes in suspended-sediment concentrations occurred for all grain-size fractions, but were most pronounced for the sand fraction. Because of these changes, the use of bulk discharge-concentration relationships (i.e., “sediment rating curves”) without time-dependencies in these relationships resulted in substantial errors in sediment load estimates, including 2.5-fold over-prediction of Eel River sediment loads since 1979. We conclude that sediment discharge and sediment discharge relationships (such as sediment rating curves) from these coastal rivers have varied substantially with time in response to land use and climate. Thus, the use of historical river sediment data and sediment rating curves without considerations for time-dependent trends may result in significant errors in sediment yield estimates from the globally-important steep, small watersheds.

California↗
USGS2004

Vestiges of an Iapetan rift basin in the New Jersey Highlands: Implfications for the Neoproterozoic Laurentian margin

Thin, discontinuous remnants of Neoproterozoic intracratonic rift-basin deposits of the Chestnut Hill Formation occur in the western New Jersey Highlands. These deposits form an important link between well-documented Iapetan rift-basins in both the northern and southern Appalachians. The close spatial relations of Chestnut Hill rocks to Paleozoic sedimentary rocks open the possibility that additional Iapetan rift-basins could be concealed beneath the rocks of the Valley and Ridge Province to the west indicating a much broader zone of rifting than has been previously proposed. The Chestnut Hill Formation is intermittently exposed along a 100 km-long band that extends northeast from Pennsylvania nearly to New York State. The lower part of the Chestnut Hill Formation is composed of interbedded lithic pebble- to boulder-conglomerate and feldspathic sandstone grading upward into interbedded phyllite, feldspathic and quartz sandstone, local paleosaprolite, quartz-pebble conglomerate, thin limestone lenses, volcanic, and volcaniclasic rocks, abundant bedded ironstone (hematite ore), and ultimately into diamictites that are interpreted as possible tilloids and containing rounded intra and extrabasinal clasts of the other lithologies. Extensive soft-sediment deformation, cross bedding, and clastic dikes are common in all but the lowest and upper facies. Banded hematite layers occur preferentially in fine-grained tuffs and tuffaceous sediments, but hematitization has affected most lithologies. Volcanic rocks consist of altered rhyolitic tuffs and lapilli tuffs that are interbedded with sediments. The Chestnut Hill Formation is interpreted to have been deposited in early alluvial, and later a complex of fluvial, lacustrine and deltaic environments. Provenance studies based upon petrographic and geochemical analysis of clastic rocks indicate that the sediments are predominantly immature and reflect derivation from local uplifted felsic basement sources in a rifted-margin tectonic setting. Low to moderate weathering of the source rocks is indicated by the geochemistry of most samples, as is the locally intense effect of hydrothermal alteration. Most occurrences of the Chestnut Hill Formation are associated with major faults that exhibit normal movement of apparent Neoproterozoic age. Rocks from the Morgan Hill fault near Easton, Pennsylvania display consistent normal shear sense and vary from low temperature S-C mylonites to breccia that contains deformed pieces of Chestnut Hill Formation. ?? 2004 Elsevier Ltd. All rights reserved.

Journal of Geodynamics↗
USGS2004

Geophysical modeling of the northern Appalachian Brompton-Cameron, Central Maine, and Avalon terranes under the New Jersey Coastal Plain

A regional terrane map of the New Jersey Coastal Plain basement was constructed using seismic, drilling, gravity and magnetic data. The Brompton-Cameron and Central Maine terranes were coalesced as one volcanic island arc terrane before obducting onto Laurentian, Grenville age, continental crust in the Taconian orogeny [Rankin, D.W., 1994. Continental margin of the eastern United States: past and present. In: Speed, R.C., (Ed.), Phanerozoic Evolution of North American Continent-Ocean Transitions. DNAG Continent-Ocean Transect Volume. Geological Society of America, Boulder, Colorado, pp. 129-218]. Volcanic island-arc rocks of the Avalon terrane are in contact with Central Maine terrane rocks in southern Connecticut where the latter are overthrust onto the Brompton-Cameron terrane, which is thrust over Laurentian basement. Similarities of these allochthonous island arc terranes (Brompton-Cameron, Central Maine, Avalon) in lithology, fauna and age suggest that they are faulted segments of the margin of one major late Precambrian to early Paleozoic, high latitude peri-Gondwana island arc designated as "Avalonia", which collided with Laurentia in the early to middle Paleozoic. The Brompton Cameron, Central Maine, and Avalon terranes are projected as the basement under the eastern New Jersey Coastal Plain based on drill core samples of metamorphic rocks of active margin/magmatic arc origin. A seismic reflection profile across the New York Bight traces the gentle dipping (approximately 20 degrees) Cameron's Line Taconian suture southeast beneath allochthonous Avalon and other terranes to a 4 sec TWTT depth (approximately 9 km) where the Avalonian rocks are over Laurentian crust. Gentle up-plunge (approximately 5 degrees) projections to the southwest bring the Laurentian Grenville age basement and the drift-stage early Paleozoic cover rocks to windows in Burlington Co. at approximately 1 km depth and Cape May Co. at approximately 2 km depths. The antiformal Shellburne Falls and Chester domes and Chain Lakes-Pelham dome-Bronson Hill structural trends, and the synformal Connecticut Valley-Gaspe structural trend can be traced southwest into the New Jersey Coastal Plain basement. A Mesozoic rift basin, the "Sandy Hook basin", and associated eastern boundary fault is identified, based upon gravity modeling, in the vicinity of Sandy Hook, New Jersey. The thickness of the rift-basin sedimentary rocks contained within the "Sandy Hook basin" is approximately 4.7 km, with the basin extending offshore to the east of the New Jersey coast. Gravity modeling indicates a deep rift basin and the magnetic data indicates a shallow magnetic basement caused by magnetic diabase sills and/or basalt flows contained within the rift-basin sedimentary rocks. The igneous sills and/or flows may be the eastward continuation of the Watchung and Palisades bodies. ?? 2004 Elsevier Ltd. All rights reserved.

Journal of Geodynamics↗
USGS2009

Evidence for an Alleghanian (Early Carboniferous to Late Permian) tectonothermal event in the New Jersey Coastal Plain basement from 40 Ar/ 39 Ar biotite data, geochemistry and gravity modeling

40 Ar/ 39 Ar dating of biotite from felsic orthogneiss recovered from the -3890-foot level of the Island Beach State Park (IBSP) well beneath the outer New Jersey Coastal Plain was accomplished using CO 2 laser incremental-heating techniques. Over 75% of the Ar released from the incremental-heating experiment form a well-behaved plateau with a calculated age of 243.98 ?? 0.10 Ma. The new 244 Ma biotite age reported here is a cooling age younger than the metamorphic event that crystallized or reheated the biotite. We consider reheating of older biotite to be unlikely because the concordant 40 Ar/ 39 Ar spectrum upon repeated incremental laser heating showed a well-developed plateau. Thus, biotites from the IBSP gneiss are interpreted as having crystallized during a single thermal event, followed by cooling to below 300 ??C. The IBSP well falls on a structural and geophysical anomaly trend that is along strike with rocks of the Bronson Hill anticlinorium to the north of the IBSP gneiss. Locally graphitic metasedimentary schists and gneisses recovered from New Jersey wells inboard of the IBSP well gneiss correlate to similar lithologies of the Connecticut Valley synclinorium west of the Hartford basin. Our reinterpretation of the IBSP gneiss as metamorphosed dacite or dacitic tuff is consistent with a correlation to some rocks of the Bronson Hill magmatic arc east of the Hartford basin. If correct, this would imply a Late Ordovician age for the protolith of the IBSP gneiss. Reported 40 Ar/ 39 Ar biotite ages of 235-253 Ma from southwestern Rhode Island, and of 238-247 Ma from southeastern Connecticut, are interpreted as cooling ages following a tectonothermal event associated with the Alleghanian orogeny (Early Carboniferous to Late Permian). Cooling ages of Alleghanian age (Early Carboniferous to Late Permian) are not recognized west of the Bronson Hill volcanic arc in either central Connecticut or in Massachusetts. Therefore, the 244 Ma cooling age presented here, and the geochemical affinity of the IBSP gneiss to some orthogneisses of the Bronson Hill arc, support an interpretation of the IBSP well as representing the southern continuation of the Bronson Hill arc into New Jersey. Moreover, it documents the presence of rocks beneath the outer New Jersey Coastal Plain that experienced a Permian Alleghanian metamorphism. ?? 2009 Elsevier Ltd. All rights reserved.

Journal of Geodynamics↗
USGS2020

The effects of fire on the thermal environment of sagebrush communities

Thermal heterogeneity provides options for organisms during extreme temperatures that can contribute to their fitness. Sagebrush ( Artemisia spp.) communities exhibit vegetation heterogeneity that creates thermal variation at fine spatial scales. However, fire can change vegetation and thereby variation within the thermal environment of sagebrush communities. To describe spatial and temporal thermal variation of sagebrush communities following wildfire, we measured black bulb temperature ( T bb ) at 144 random points dispersed within unburned and burned communities, for 24-h at each random point. We observed a wide thermal gradient in unburned (−7.3° to 63.3 °C) and burned (−4.6° to 64.8 °C) sagebrush communities. Moreover, unburned and burned sagebrush communities displayed high thermal heterogeneity relative to ambient temperature ( T air ). Notably, T bb varied by 47 °C in both unburned and burned communities when T air was 20 °C. However, fire greatly reduced the buffering capacity and thermal refuge of Wyoming big sagebrush ( A. tridentata wyomingensis ) communities during low and high T air . Furthermore, fire increased T bb in Wyoming big sagebrush and mountain big sagebrush ( A. t. vaseyana ) during the mid-day hours. These results demonstrate how fire changes the thermal environment of big sagebrush communities and the importance of shrub structure which can provide thermal refuge for organisms in burned communities during extreme low and high T air .

Journal of Thermal Biology↗
USGS2014

Towards understanding the puzzling lack of acid geothermal springs in Tibet (China): Insight from a comparison with Yellowstone (USA) and some active volcanic hydrothermal systems

Explanations for the lack of acid geothermal springs in Tibet are inferred from a comprehensive hydrochemical comparison of Tibetan geothermal waters with those discharged from Yellowstone (USA) and two active volcanic areas, Nevado del Ruiz (Colombia) and Miravalles (Costa Rica) where acid springs are widely distributed and diversified in terms of geochemical characteristic and origin. For the hydrothermal areas investigated in this study, there appears to be a relationship between the depths of magma chambers and the occurrence of acid, chloride-rich springs formed via direct magmatic fluid absorption. Nevado del Ruiz and Miravalles with magma at or very close to the surface (less than 1–2 km) exhibit very acidic waters containing HCl and H 2 SO 4 . In contrast, the Tibetan hydrothermal systems, represented by Yangbajain, usually have fairly deep-seated magma chambers so that the released acid fluids are much more likely to be fully neutralized during transport to the surface. The absence of steam-heated acid waters in Tibet, however, may be primarily due to the lack of a confining layer (like young impermeable lavas at Yellowstone) to separate geothermal steam from underlying neutral chloride waters and the possible scenario that the deep geothermal fluids below Tibet carry less H 2 S than those below Yellowstone.

Tibet, Yellowstone National Park↗
USGS2016

Multireaction equilibrium geothermometry: A sensitivity analysis using data from the Lower Geyser Basin, Yellowstone National Park, USA

A multireaction chemical equilibria geothermometry (MEG) model applicable to high-temperature geothermal systems has been developed over the past three decades. Given sufficient data, this model provides more constraint on calculated reservoir temperatures than classical chemical geothermometers that are based on either the concentration of silica (SiO 2 ), or the ratios of cation concentrations. A set of 23 chemical analyses from Ojo Caliente Spring and 22 analyses from other thermal features in the Lower Geyser Basin of Yellowstone National Park are used to examine the sensitivity of calculated reservoir temperatures using the GeoT MEG code (Spycher et al. 2013, 2014) to quantify the effects of solute concentrations, degassing, and mineral assemblages on calculated reservoir temperatures. Results of our analysis demonstrate that the MEG model can resolve reservoir temperatures within approximately ±15°C, and that natural variation in fluid compositions represents a greater source of variance in calculated reservoir temperatures than variations caused by analytical uncertainty (assuming ~5% for major elements). The analysis also suggests that MEG calculations are particularly sensitive to variations in silica concentration, the concentrations of the redox species Fe(II) and H 2 S, and that the parameters defining steam separation and CO 2 degassing from the liquid may be adequately determined by numerical optimization. Results from this study can provide guidance for future applications of MEG models, and thus provide more reliable information on geothermal energy resources during exploration.

Idaho, Montana, Wyoming↗