SEARCH · USGS Science

Results for “Interpretation”

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

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8Linked to original sources

USGS2023

Approaching the upper boundary of driver-response relationships: Identifying factors using a novel framework integrating quantile regression with interpretable machine learning

The identification of factors that may be forcing ecological observations to approach the upper boundary provides insight into potential mechanisms affecting driver-response relationships, and can help inform ecosystem management, but has rarely been explored. In this study, we propose a novel framework integrating quantile regression with interpretable machine learning. In the first stage of the framework, we estimate the upper boundary of a driver-response relationship using quantile regression. Next, we calculate “potentials” of the response variable depending on the driver, which are defined as vertical distances from the estimated upper boundary of the relationship to observations in the driver-response variable scatter plot. Finally, we identify key factors impacting the potential using a machine learning model. We illustrate the necessary steps to implement the framework using the total phosphorus (TP)-Chlorophyll a (CHL) relationship in lakes across the continental US. We found that the nitrogen to phosphorus ratio (N:P), annual average precipitation, total nitrogen (TN), and summer average air temperature were key factors impacting the potential of CHL depending on TP. We further revealed important implications of our findings for lake eutrophication management. The important role of N:P and TN on the potential highlights the co-limitation of phosphorus and nitrogen and indicates the need for dual nutrient criteria. Future wetter and/or warmer climate scenarios can decrease the potential which may reduce the efficacy of lake eutrophication management. The novel framework advances the application of quantile regression to identify factors driving observations to approach the upper boundary of driver-response relationships.

Frontiers of Environmental Science & Engineering↗

USGS2022

From data to interpretable models: Machine learning for soil moisture forecasting

Soil moisture is critical to agricultural business, ecosystem health, and certain hydrologically driven natural disasters. Monitoring data, though, is prone to instrumental noise, wide ranging extrema, and nonstationary response to rainfall where ground conditions change. Furthermore, existing soil moisture models generally forecast poorly for time periods greater than a few hours. To improve such forecasts, we introduce two data-driven models, the Naive Accumulative Representation (NAR) and the Additive Exponential Accumulative Representation (AEAR). Both of these models are rooted in deterministic, physically based hydrology, and we study their capabilities in forecasting soilmoisture over time periods longer than a fewhours. Learned model parameters represent the physically based unsaturated hydrological redistribution processes of gravity and suction. We validate our models using soil moisture and rainfall time series data collected from a steep gradient, post-wildfire site in southern California. Data analysis is complicated by rapid landscape change observed in steep, burned hillslopes in response to even small to moderate rain events. The proposed NAR and AEAR models are, in forecasting experiments, shown to be competitive with several established and state-of-the-art baselines. The AEAR model fits the data well for three distinct soil textures at variable depths below the ground surface (5, 15, and 30 cm). Similar robust results are demonstrated in controlled, laboratory-based experiments. Our AEAR model includes readily interpretable hydrologic parameters and provides more accurate forecasts than existing models for time horizons of 10–24 h. Such extended periods of warning for natural disasters, such as floods and landslides, provide actionable knowledge to reduce loss of life and property.

International Journal of Data Science and Analytic↗

USGS1983

An interpretation of carbon and sulfur relationships in Black Sea sediments as indicators of environments of deposition

Syngenetic iron sulfides in sediments are formed from dissolved sulfide resulting from sulfate reduction and catabolism of organic matter by anaerobic bacteria. It has been shown that in recent marine sediments deposited below oxygenated waters there is a constant relationship between reduced sulfur and organic carbon which is generally independent of the environment of deposition. Reexamination of data from recent sediments from euxinic marine environments ( e.g. , the Black Sea) also shows a linear relationship between carbon and sulfur, but the slope is variable and the line intercepts the S axis at a value between 1 and 2 percent S. It is proposed that the positive S intercept is due to watercolumn microbial reduction of sulfate using metabolizable small organic molecules and the sulfide formed is precipitated and accumulates at the sediment-water interface. The variation in slope and intercept of the C to S plots for several cores and for different stratigraphic zones for the Black Sea can be interpreted in relation to thickness of the aqueous sulfide layer or thinness of the oxygen containing layer and to deposition rate, but also may be influenced by availability of iron, and perhaps the type of organic matter (Leventhal, 1979).

Geochimica et Cosmochimica Acta↗

USGS1986

Europa: Characterization and interpretation of global spectral surface units

The Voyager global multispectral mosaic of the Galilean satellite Europa ( T. V. Johnson, L. A. Soderblom, J. A. Mosher, G. E. Danielson, A. F. Cook, and P. Kupferman, 1983, J. Geophys. Res. 88, 5789–5805 ) was analyzed to map surface units with similar optical properties ( T. B. McCord, M. L. Nelson, R. N. Clark, A. Meloy, W. Harrison, T. V. Johnson, D. L. Matson, J. A. Mosher, and L. Soderblom, 1982, Bull Amer. Astron. Soc. 14, 737 ). Color assignments in the unit map are indicative of the spectral nature of the unit. The unit maps make it possible to infer extensions of the geologic units mapped by B. K. Lucchitta and L. A. Soderblom (1982, in Satellites of Jupiter , pp. 521–555, Univ. of Arizona Press, Tucson) beyond the region covered in the high-resolution imagery. The most striking feature in the unit maps is a strong hemispheric asymmetry. It is seen most clearly in the ultraviolet/violet albedo ratio image, because the asymmetry becomes more intense as the wavelength decreases. It appears as if the surface has been darkened, most intensely in the center of the trailing hemisphere and decreasing gradually, essentially as the cosine of the angle from the antapex of motion, to a minimum in the center of the leading hemisphere. The cosine pattern suggests that the darkening is exogenic in origin and is interpreted as evidence of alteration of the surface by ion bombardment from the Jovian magnetosphere.

Icarus↗

USGS1980

Sediment waves on the Monterey fan levee: a preliminary physical interpretation.

Sediment waves on the deep ocean floor occur mostly on the lower continental rise on slopes of 1° or less. Previous studies show that their amplitude and wavelength vary greatly, but little is known about their shape in plan. A detailed survey of a 30-km 2 area of abyssal-depth sediment waves associated with the levee of the Monterey fan valley shows a pattern of sinuous crests and troughs with parallel, well-bedded internal structure. Material in the upper 1 m of sediment consists predominantly of bioturbated, muddy coccolith ooze. A single thin, silty horizon can be correlated between adjoining waves. On the basis of measured wave dimensions and an estimate of flow velocity we use a simple two-layer model for water movement to infer approximate flow parameters. The sediment waves are considered to be formed most likely by low-velocity (10 cm/s), low-concentration turbidity flows approximately 100–800 m thick. This interpretation emphasizes the role of low-speed, low-concentration turbidity currents in the downslope movement of fine-grained material. This type of transport—deposition regime explains the formation of sediment waves very well except for certain waves occurring on depositional ridges in the deep ocean.

Marine Geology↗

USGS1989

Heavy-mineral reconnaissance off the coast of the Apalachicola river delta, northwest Florida: A summary and new interpretations

Two-hundred-and-fifty sediment samples were collected for heavy-mineral and textural analysis along the northwest Florida coastline from approximately 24 km offshore of Apalachee Bay to the same distance offshore of Pensacola Bay. The heavy-mineral suite characterizing sediments within this region consists of opaque minerals, kyanite, staurolite, tourmaline, zircon and rutile. Minor constituents of this suite include epidote, sphene, amphibole, sillimanite, garnet and leucoxene. The average heavy-mineral concentration within these sediments is approximately 0.12 wt.%. Specifically, the 2 to 3∅ grain-size interval contains an average of 0.51 wt.%, whereas the 3 to 4∅ interval contains an average of 4.39 wt.% heavy minerals. Note that the 3 to 4∅ interval typically represents only 4% of the sample volume. There is a general westward increase in heavy-mineral concentrations throughout the study area. Superimposed on this regional trend, areas of maximum heavy-mineral concentration occur within sediments offshore of St. George and Santa Rosa islands. The primary source of sediments in the region is the crystalline rocks of the southern Appalachians. Granulometric analyses of these sediments reveal a westward increase in values of sample mean grain size, and decrease in standard deviation, and percent fines. It is postulated from these data, in addition to the interpretation of sample grain-size distributions, skewness, and kurtosis, that these inner continental shelf sediments are primarily fluvial in origin. These sediments have been transported to the shelf by the Apalachicola and surrounding major rivers during Pleistocene low sea-level stands. Data also indicate evidence of reworking by coastal or marine offshore wave processes.

Marine Geology↗

USGS1989

High-resolution seismic-reflection interpretations of some sediment deposits, Antarctic continental margin: Focus on the western Ross Sea

High-resolution seismic-reflection data have been used to a varying degree by geoscientists to interpret the history of marine sediment accumulations around Antarctica. Reconnaissance analysis of 1-, 3.5-, and 12-kHz data collected by the U.S. Geological Survey in the western Ross Sea has led to the identification of eight echo-character facies and six microtopographic facies in the sediment deposits that overlie the Ross Sea unconformity. Three depositional facies regions, each characterized by a particular assemblage of echo-character type and microtopographic facies, have been identified on the continental shelf. These suites of acoustic facies are the result of specific depositional processes that control type and accumulation of sediment in a region. Evidence of glacial processes and products is uncommon in regions 1 and 2, but is abundant in region 3. McMurdo Sound, region 1, is characterized by a monospecific set of acoustic facies. This unique assemblage probably represents turbidity current deposition in the western part of the basin. Most of the seafloor in region 2, from about latitude 77°S to 75°S, is deeper than 600 m below sealevel. The microtopographic facies and echo-character facies observed on the lower slopes and basin floor there reflect the thin deposits of pelagic sediments that have accumulated in the low-energy conditions that are typical of deep-water environments. In shallower water near the boundary with region 3, the signature of the acoustic facies is different from that in deeper water and probably indicates higher energy conditions or, perhaps, ice-related processes. Thick deposits of tills emplaced by lodgement during the most recent advance of the West Antarctic Ice Sheet are common from latitude 75°S to the northern boundary of the study area just south of Coulman Island (region 3). The signature of microtopographic facies in this region reflects the relief of the base of the grounded ice sheet prior to decoupling from the seafloor. Current winnowing and scour of shallow parts of the seafloor inhibits sediment deposition and maintains the irregular, hummocky relief that characterizes much of the region. Seafloor relief of this type in other polar areas could indicate the former presence of grounded ice.

Marine Geology↗

USGS1978

Interpretation of a Landsat image of an unusual flood phenomenon in Australia

A Landsat image of part of the flooded area of Cooper Creek, Queensland, Australia, in February 1974, shows large dark areas within the flooded valley. The dark areas are believed to be wet, but unflooded, areas of dark alluvial soil. These striking features, which have not previously been identified on Landsat images, must be properly interpreted so as not to confuse them with clear water.

Queensland↗

USGS1980

Nearshore current pattern off south Texas: an interpretation from aerial photographs.

Current patterns in a 4-km-wide zone along the south Texas coast were interpreted from patterns of water turbidity visible in aerial photographs taken during a winter day of moderate northerly winds. Features of the turbidity pattern remained recognizable on photographs taken 25 min apart. Currents measured from the movements of these features were southward and nearly parallel to shore, increasing from about 17 cm/sec in an offshore zone to about 40 cm/sec at the line of breaking waves. - from Authors

Remote Sensing of Environment↗

USGS1982

Further interpretation of the relation of organochlorine residues in brown pelican eggs to reproductive success

This study was conducted to provide additional interpretation of the relation of organochlorine pollutants to reproductive success and population stability of brown pelicans Pelecanus occidentalis in the southeastern United States from 1969 to 1977. The sample egg technique was employed in South Carolina; it involved collecting an egg from each of 156 marked nests, analysing the eggs for residues and relating residue levels to nest success. All of the organochlorines appeared to induce adverse effects on reproductive success of brown pelicans, but the correlative evidence indicated that DDE was responsible for most of the pollutant-related nest failure. Close intercorrelation of each of the organochlorines presented a problem in the separation of the effects of each pollutant. The critical level of DDE in brown pelican eggs was 3 μg g −1 , this level was associated with substantially impaired reproductive success and 4 μg g −1 of DDE was associated with sotal reproductive failure. The critical level for dieldrin was not determined, but it appeared to exceed 1 μg g −1 . Dieldrin, DDT, DDD and the PCBs apparently exerted a minor effect on reproductive success of brown pelicans in South Carolina. Endrin seemed to induce reproductive impairment in brown pelicans in Louisiana in 1975. The critical level in eggs is only roughly estimated at <0·5 μg g −1 because it was not possible to use the sample egg technique in this small population. Significant declines in organochlorine residues in South Carolina during this study were associated with increases in eggshell thickness, reproductive success and the breeding population. Of those avian species studied, the brown pelican remains the most sensitive to organochlorine contaminants, particularly DDE and endrin; factors related to this sensitivity are poorly known.

South Carolina↗

USGS1985

The copper-nickel concentration log: A tool for stratigraphic interpretation within the ultramafic and basal zones of the stillwater complex, Montana

An analogue to the electric well log was devised for copper-nickel concentration drill-hole data from the Basal and lower part of the Ultramafic zones of the Stillwater Complex using automated data processing. The copper-nickel concentration logs graphically represent intensity (concentration) values that reflect the distribution of the elements in sulfide and silicate minerals. Four major patterns are recognized by their characteristic variations in copper and nickel intensity: (1) relatively flat, low-level copper-intensity signatures associated with arcuate nickel-intensity patterns that correlate with rocks in the Peridotite member of the Ultramafic zone; (2) arcuate or bulb-like patterns of copper and nickel intensity that correlate closely with the Basal bronzite cumulate member of the Basal zone; (3) complex patterns consisting of intervals of low-intensity copper and moderate-intensity nickel, spikes of high nickel and copper intensity, and high copper intensity associated with low nickel intensity that correlate respectively with cordierite-pyroxene hornfels, massive sulfide, norites and mineralized diabase dikes in the Basal norite member; and (4) large intervals of extremely low copper and nickel intensity that correlate with quartz-orthopyroxene hornfels. The recognition and interpretation of these patterns allow two- and three-dimensional stratigraphic and lithologic reconstructions to be done by means of concentration-log correlations instead of variable quality lithologic logging.

Journal of Geochemical Exploration↗

USGS1987

Mineralogical basis for the interpretation of multi-element (ICP-AES), oxalic acid, and aqua regia partial digestions of stream sediments for reconnaissance exploration geochemistry

We have applied partial digestion procedures, primarily oxalic acid and aqua regia leaches, to several regional geochemical reconnaissance studies carried out using Inductively Coupled Plasma-Atomic Emission Spectroscopy (ICP-AES) analytical methods. We have chosen to use these two acids because the oxalic acid primarily attacks those compounds formed during secondary geochemical processes, whereas aqua regia will digest the primary sulfide phases as well as secondary phases. Application of the partial digestion technique has proven superior to total digestion because the concentration of metals in hydromorphic compounds and the sulfides is enhanced relative to the metals bound in the unattacked silicate phases. The aqua regia digestion attacks and leaches metals from the mafic chain silicates and the phyllosilicates (coordination number of VI or more), yielding a characteristic geochemical signature, but does not leach appreciable metal from many other silicates. In order to interpret the results from these leach studies, we have initiated an investigation of a large suite of hand-picked mineral separates. The study includes analyses of about two hundred minerals representing the common rock-forming minerals as well as end-member compositions of various silicates, oxides, sulfides, carbonates, sulfates, and some vanadates, molybdates, tungstates, and phosphates. The objective of this study is to evaluate the effect of leaching by acids of particular lattice sites in specific mineral structures.

Journal of Geochemical Exploration↗

USGS2013

Using simulated maps to interpret the geochemistry, formation and quality of the Blue Gem Coal Bed, Kentucky, USA

This study presents geostatistical simulations of coal-quality parameters, major oxides and trace metals for an area covering roughly 812 km 2 of the Blue Gem coal bed in southeastern Kentucky, USA. The Blue Gem, characterized by low ash yield and low sulfur content, is an important economic resource. Past studies have characterized the Blue Gem's geochemistry, palynology and petrography and inferred a depositional setting of a planar peat deposit that transitioned to slightly domed later in its development. These studies have focused primarily on vertical geochemical trends within the coal bed. Simulated maps of chemical elements derived from 45 measured sample locations across the study area provide an opportunity to observe changes in the horizontal direction within the coal bed. As the Blue Gem coal bed shows significant vertical chemical trends, care was taken in this study to try to select samples from a single, middle portion of the coal. By revealing spatial distribution patterns of elements across the middle of the bed, associations between different components of the coal can be seen. The maps therefore help to provide a picture of the coal-forming peat bog at an instant in geologic time and allow interpretation of a depositional setting in the horizontal direction. Results from this middle portion of the coal suggest an association of SiO 2 with both K 2 O and TiO 2 in different parts of the study area. Further, a pocket in the southeast of the study area shows elevated concentrations of elements attributable to observed carbonate-phase minerals (MgO, CaO, Ba and Sr) as well as elements commonly associated with sulfide-phase minerals (Cu, Mo and Ni). Areas of relatively high ash yield are observed in the north and south of the mapped area, in contrast to the low ash yields seen towards the east. Additionally, we present joint probability maps where multiple coal-quality parameters are plotted simultaneously on one figure. This application allows researchers to investigate associations of more than two components in a straight-forward manner useful in guiding resource exploration.

Kentucky↗

USGS2017

Minimizing effects of methodological decisions on interpretation and prediction in species distribution studies: An example with background selection

Evaluating the conditions where a species can persist is an important question in ecology both to understand tolerances of organisms and to predict distributions across landscapes. Presence data combined with background or pseudo-absence locations are commonly used with species distribution modeling to develop these relationships. However, there is not a standard method to generate background or pseudo-absence locations, and method choice affects model outcomes. We evaluated combinations of both model algorithms (simple and complex generalized linear models, multivariate adaptive regression splines, Maxent, boosted regression trees, and random forest) and background methods (random, minimum convex polygon, and continuous and binary kernel density estimator (KDE)) to assess the sensitivity of model outcomes to choices made. We evaluated six questions related to model results, including five beyond the common comparison of model accuracy assessment metrics (biological interpretability of response curves, cross-validation robustness, independent data accuracy and robustness, and prediction consistency). For our case study with cheatgrass in the western US, random forest was least sensitive to background choice and the binary KDE method was least sensitive to model algorithm choice. While this outcome may not hold for other locations or species, the methods we used can be implemented to help determine appropriate methodologies for particular research questions.

Ecological Modelling↗

USGS2026

Sackung at Bald Eagle ridge, central Colorado: An updated interpretation of ridge-spreading movement, structures, and mechanisms from 50 years (1975–2025) of U.S. Geological Survey research

Slow gravitational failures of mountain peaks and ridges are poorly understood. Herein, we report on 50 years of studies at a slowly spreading castellate ridge in the Sawatch Range in central Colorado. The orientations of geomorphic-structural features indicate that the fractured Precambrian granitic rock underlying the ridge has extended and spread northwestward toward the formerly glacier-covered Busk Creek valley. Results from surveying, field-based geomorphic-structural mapping using lidar, rock mass quality measurements, a passive seismic survey, and satellite radar provide a major update to research started by U.S. Geological Survey researchers in the 1970s and 1980s. New insights include a recognition that the entire ridge has slowly moved by concurrent sliding along an inferred northwest dipping, compound basal-slip surface (or zone), and through the formation of multiple grabens by normal faulting and flexural toppling along sets of pre-existing fractures that dip moderately (∼45°) to the southeast. We were unable to distinguish the presence of a sudden and strong contrast in seismic velocity across the inferred slip surface. Movement during the 50-year study period has been episodic and gradually decreasing, in correspondence with decreasing cumulative annual precipitation and increasing mean annual air temperatures. The fastest moving area, just upslope from the glacier trimline, had an average horizontal velocity of 3 – 4 mm/yr. Evidence suggests that movement started as a paraglacial response mechanism, but because of the site’s proximity to the Rio Grande Rift, we cannot exclude earthquake shaking as a mechanism for initiation or enhancement of slope movement. An estimate of longer-term horizontal movement from the exposed basal-slip surface at the uphill side of the ridgetop graben is ∼1.1 mm/yr for the 13 – 14 ky post-glacial period. Broad implications of our work are that: (1) long-term measurements (decades or longer) of slope movement can add insights into how sackungen form and evolve through time; (2) the identification of thrust faults and toes in zones of compression near valley bottoms can be crucial for interpreting sackung failure mechanisms, and (3) the use of passive-seismic techniques to identify the depth to a slip surface may not be successful in granitic terrain dominated by planar fractures and subtle changes in rock-mass characteristics.

Colorado↗

USGS2008

Application of sedimentary-structure interpretation to geoarchaeological investigations in the Colorado River Corridor, Grand Canyon, Arizona, USA

We present a detailed geoarchaeological study of landscape processes that affected prehistoric formation and modern preservation of archaeological sites in three areas of the Colorado River corridor in Grand Canyon, Arizona, USA. The methods used in this case study can be applied to any locality containing unaltered, non-pedogenic sediments and, thus, are particularly relevant to geoarchaeology in arid regions. Resolving the interaction of fluvial, aeolian, and local runoff processes in an arid-land river corridor is important because the archaeological record in arid lands tends to be concentrated along river corridors. This study uses sedimentary structures and particle-size distributions to interpret landscape processes; these methods are commonplace in sedimentology but prove also to be valuable, though less utilized, in geoarchaeology and geomorphology. In this bedrock canyon, the proportion of fluvial sediment generally decreases with distance away from the river as aeolian, slope-wash, colluvial, and debris-flow sediments become more dominant. We describe a new facies consisting of 'flood couplets' that include a lower, fine-grained fluvial component and an upper, coarser, unit that reflects subaerial reworking at the land surface between flood events. Grain-size distributions of strata that lack original sedimentary structures are useful within this river corridor to distinguish aeolian deposits from finer-grained fluvial deposits that pre-date the influence of the upstream Glen Canyon Dam on the Colorado River. Identification of past geomorphic settings is critical for understanding the history and preservation of archaeologically significant areas, and for determining the sensitivity of archaeological sites to dam operations. Most archaeological sites in the areas studied were formed on fluvial deposits, with aeolian deposition acting as an important preservation agent during the past millennium. Therefore, the absence of sediment-rich floods in this regulated river, which formerly deposited large fluvial sandbars from which aeolian sediment was derived, has substantially altered processes by which the prehistoric, inhabited landscape formed, and has also reduced the preservation potential of many significant cultural sites.

Geomorphology↗

USGS2023

Large-scale variation in lakebed properties interpreted from single-beam sonar in two Laurentian Great Lakes

Acoustic seabed classification (ASC) is an important method for understanding landscape-level physical and biological patterns in the aquatic environment. Bottom habitats in the Laurentian Great Lakes are poorly mapped to date, and will require a variety of contributors and data sources to complete. We repurposed a long-term split-beam echosounder dataset gathered for purposes of fisheries assessment to estimate lakebed properties utilizing unsupervised classification of echo return data. We interpreted first echo properties representing lakebed hardness and roughness to define and map three statistically supported lakebed classes revealed through cluster analysis. Our results indicate coherent and logical class boundaries and suggest that the dataset has promise for expanded use in ASC. To improve inferences using repeated measures, future work should focus on collecting ground truth information for areas previously surveyed and on collecting concurrent ground truth information when sampling acoustic data moving forward.

Michigan, Ontario, Wisconsin↗

USGS2010

Snowmelt hydrograph interpretation: Revealing watershed scale hydrologic characteristics of the Yellowstone volcanic plateau

Snowmelt hydrograph analysis and groundwater age dates of cool water springs on the Yellowstone volcanic plateau provide evidence of high volumes of groundwater circulation in watersheds comprised of quaternary Yellowstone volcanics. Ratios of maximum to minimum mean daily discharge and average recession indices are calculated for watersheds within and surrounding the Yellowstone volcanic plateau. A model for snowmelt recession is used to separate groundwater discharge from overland runoff, and compare groundwater systems. Hydrograph signal interpretation is corroborated with chlorofluorocarbon (CFC) and tritium concentrations in cool water springs on the Yellowstone volcanic plateau. Hydrograph parameters show a spatial pattern correlated with watershed geology. Watersheds comprised dominantly of quaternary Yellowstone volcanics are characterized by slow streamflow recession, low maximum to minimum flow ratios. Cool springs sampled within the Park contain CFC's and tritium and have apparent CFC age dates that range from about 50 years to modern. Watersheds comprised of quaternary Yellowstone volcanics have a large volume of active groundwater circulation. A large, advecting groundwater field would be the dominant mechanism for mass and energy transport in the shallow crust of the Yellowstone volcanic plateau, and thus control the Yellowstone hydrothermal system. ?? 2009 Elsevier B.V.

Journal of Hydrology↗