SEARCH · USGS Science

Results for “Geophysics”

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 1,405 records · Page 78Linked to original sources

USGS2014

Sources of global climate data and visualization portals

Climate is integral to the geophysical foundation upon which ecosystems are structured. Knowledge about mechanistic linkages between the geophysical and biological environments is essential for understanding how global warming may reshape contemporary ecosystems and ecosystem services. Numerous global data sources spanning several decades are available that document key geophysical metrics such as temperature and precipitation, and metrics of primary biological production such as vegetation phenology and ocean phytoplankton. This paper provides an internet directory to portals for visualizing or servers for downloading many of the more commonly used global datasets, as well as a description of how to write simple computer code to efficiently retrieve these data. The data are broadly useful for quantifying relationships between climate, habitat availability, and lower-trophic-level habitat quality - especially in Arctic regions where strong seasonality is accompanied by intrinsically high year-to-year variability. If defensible linkages between the geophysical (climate) and the biological environment can be established, general circulation model (GCM) projections of future climate conditions can be used to infer future biological responses. Robustness of this approach is, however, complicated by the number of direct, indirect, or interacting linkages involved. For example, response of a predator species to climate change will be influenced by the responses of its prey and competitors, and so forth throughout a trophic web. The complexities of ecological systems warrant sensible and parsimonious approaches for assessing and establishing the role of natural climate variability in order to substantiate inferences about the potential effects of global warming.

Conference Paper↗

USGS2018

An on-campus well field for hydrogeophysics education and undergraduate research

The emerging subdiscipline of hydrogeophysics is underdeveloped in undergraduate curricula relative to its importance in professional engineering/environmental practice. In 2001, the Bucknell Department of Geology initiated efforts to refocus an undergraduate geophysics course on near-surface geophysical methods for hydrologic, environmental, and engineering problems. In addition to offering students practical experience, treatment of hydrogeophysics provides important pedagogical opportunities. Field-based hydrogeophysics labs challenge students to integrate concepts from other geology courses, as well as from physics, math, and chemistry. We faced two challenges in revising our geophysics course: (1) access to wells for field exercises on borehole geophysics; and (2) the costs of acquiring and maintaining equipment. We pursued two strategies to solve these problems. First, we established an on-campus well field, which serves as a field laboratory for downhole and cross-hole experiments. Second, we incorporated field demonstrations and lectures by professional geoscientists, including alumni, into our courses. By adding field exercises to our syllabi and promoting undergraduate research, we are building a cutting-edge dataset that includes televiewer and standard wellbore logging, cross-hole tomography, and hydraulic-test data. Student-led analysis of these data has already provided valuable insights into the control of fractures on aquifer properties, and these observations are being combined with outcrop studies to place our results in a more regional context.

Pennsylvania↗

USGS1993

Evaluation of historical scour at selected stream crossings in Indian

Geophysical data were collected by means of ground-penetrating radar and tuned transducer systems to estimate the historical scour at ten bridges in Indiana. These geophysical data were used to compare and evaluate the results of 13 published pier-scour equations. In order to make this comparison, it was assumed that the measured historical scour was associated with the peak historical discharge. Because the geophysical data were not sufficient to map the lateral extent of the refilled scour hole, local scour could not be isolated from concentration scour. For the evaluation, computed contraction scour and pier scour were used in combination with the existing channel geometry to determine a computed bed elevation. This computed bed elevation was compared to be minimum historic bed elevation estimated from the geophysical data. None of the selected pier-scour equations, when combined with the contraction-scour equation, accurately represented the historical scour at all of the study sites. On the basis of the limited data presented, the equations currently recommended by the Federal Highway Administration provided a combination of accuracy and safety, required by design equations, equal to or better than the other equations evaluated.

Conference Paper↗

USGS2018

2D and 3D potential field mapping and modelling at the Fallon FORGE site, Nevada, USA

Accurate geological characterization of Fallon FORGE is important for preparing the site as an EGS laboratory. As part of this effort, a 3D geologic map was constructed previously from well logs, surface geologic mapping, 2D seismic profiles, interpreted gravity & magnetic maps, and a gravity-inferred basement surface. In this study, we have conducted both 2D and 3D modelling of high-resolution gravity and magnetic data (pre-existing and new) in an effort to further refine and test this 3D geologic map at Fallon. This effort enabled a direct comparison of the 2D and 3D model results. Potential field modelling was guided by rock-property measurements of samples from drill-core and outcrop from the Fallon area. In total, five 2D potential field model profiles, up to 30 km in length, were constructed that extend across the Fallon area. The 3D gravity model volume was 8 km (N-S) x 8 km (E-W) x 4 km (thick). The majority of the 3D gravity model volume had 100 m cubic cells; but cells near the land surface were 1 m thick to adequately capture topography. Overall, the 2D & 3D geophysical modelling largely confirmed the previously constructed 3D geologic map at Fallon for three reasons: 1) lithologic boundaries in the 2D & 3D density models mostly agree with those in the 3D geologic map, 2) the rock properties used in the models lie within the range of independent measurements made on representative rock samples from the region, and 3) the match between the observed and calculated anomalies are largely within the measurement error of the observed fields. In places where the geophysical and geologic models differ, geophysical model results have revealed subsurface structural features that have helped refine geologic interpretations which, in turn, lead to adjustments to the 3D geologic map. In this paper, we present the 2D & 3D geophysical model results and discuss how they were utilized to confirm and refine our 3D geologic understanding of the Fallon FORGE site.

Geothermal Resources Council Transactions↗

USGS2005

The inverse problem of refraction travel times, part II: Quantifying refraction nonuniqueness using a three-layer model

This paper is the second of a set of two papers in which we study the inverse refraction problem. The first paper, "Types of Geophysical Nonuniqueness through Minimization," studies and classifies the types of nonuniqueness that exist when solving inverse problems depending on the participation of a priori information required to obtain reliable solutions of inverse geophysical problems. In view of the classification developed, in this paper we study the type of nonuniqueness associated with the inverse refraction problem. An approach for obtaining a realistic solution to the inverse refraction problem is offered in a third paper that is in preparation. The nonuniqueness of the inverse refraction problem is examined by using a simple three-layer model. Like many other inverse geophysical problems, the inverse refraction problem does not have a unique solution. Conventionally, nonuniqueness is considered to be a result of insufficient data and/or error in the data, for any fixed number of model parameters. This study illustrates that even for overdetermined and error free data, nonlinear inverse refraction problems exhibit exact-data nonuniqueness, which further complicates the problem of nonuniqueness. By evaluating the nonuniqueness of the inverse refraction problem, this paper targets the improvement of refraction inversion algorithms, and as a result, the achievement of more realistic solutions. The nonuniqueness of the inverse refraction problem is examined initially by using a simple three-layer model. The observations and conclusions of the three-layer model nonuniqueness study are used to evaluate the nonuniqueness of more complicated n-layer models and multi-parameter cell models such as in refraction tomography. For any fixed number of model parameters, the inverse refraction problem exhibits continuous ranges of exact-data nonuniqueness. Such an unfavorable type of nonuniqueness can be uniquely solved only by providing abundant a priori information. Insufficient a priori information during the inversion is the reason why refraction methods often may not produce desired results or even fail. This work also demonstrates that the application of the smoothing constraints, typical when solving ill-posed inverse problems, has a dual and contradictory role when applied to the ill-posed inverse problem of refraction travel times. This observation indicates that smoothing constraints may play such a two-fold role when applied to other inverse problems. Other factors that contribute to inverse-refraction-problem nonuniqueness are also considered, including indeterminacy, statistical data-error distribution, numerical error and instability, finite data, and model parameters. ?? Birkha??user Verlag, Basel, 2005.

Pure and Applied Geophysics↗

USGS2011

Linking middle-school teachers to Earthscope

EarthScope is a multidisciplinary geophysical investigation of the structure and deformation of the North American continent. Components include a transportable array of digital seismometers (USArray) that will image the continental crust and underlying mantle, and Global Positioning System (GPS) receivers and strainmeters (Plate Boundary Observatory, PBO) that will measure tectonic and volcanic deformations Translating geophysics for kindergarten through 12th grade (K‐12) teachers and their students involves inviting them to explore new understandings of how Earth works. When they construct their understanding of Earth science through geophysical problem solving, these learners can engage in Earth‐Scope, which offers unprecedented public access to data from seismology and geodesy.

Eos, Transactions, American Geophysical Union↗

USGS2006

Combined time-series resistivity and geochemical tracer techniques to examine submarine groundwater discharge at Dor Beach, Israel

A high-resolution, stationary geophysical and geochemical survey was conducted at Dor Beach, Israel, to examine the shallow coastal hydrogeology and its control on the exchange of submarine groundwater with the shallow Mediterranean Sea. Time-series resistivity profiles using a new 56 electrode (112-m long) marine cable produced detailed profiles of the fresh water/salt water interface and the subtle response of this interface to tidal excursions and other forcing factors. Such information, when ground-truthed with representative pore water salinities and formation resistivity values, can provide unique information of the extent and rates of submarine groundwater discharge (SGD). Time-series 222Rn measurements of the adjacent coastal water column complemented these geophysical techniques and were modeled to yield integrated advective flow rates across the sediment/water interface, which at Dor Beach ranged from about 0 to 30 cm day-1 (mean = 7.1 cm d-1), depending on the tidal range. Such results suggest that the underlying hydrogeologic framework at Dor is favorable for substantial SGD. Extrapolating these SGD estimates across a 100-m wide coastal zone implies that the Rn-derived SGD rate would equal ???7.1 m3 d-1 per m of shoreline, and that the source of this discharging groundwater is a complex mixture of fresh groundwater derived from the upland Kurkar deposits, as well as locally recycled seawater. Copyright 2006 by the American Geophysical Union.

Geophysical Research Letters↗

USGS1997

Introduction to special section: The Trans-Alaska Crustal Transect (TACT) across Arctic Alaska

This special section of the Journal of Geophysical Research addresses the composition and structural evolution of the lithosphere in northern Alaska. Investigations reported in this section were mainly undertaken as part of the Trans‐Alaska Crustal Transect (TACT), an integrated geological and geophysical transect of the entire Alaskan lithosphere along a north‐south corridor undertaken from 1984 to 1992 (Figure 1). The onshore segment of the transect approximately follows along the route of the trans‐Alaskan pipeline; the offshore segment extends across the continental margin in the Gulf of Alaska to the Pacific plate. The TACT line is unique in that it provides a coordinated onshore/offshore geological and geophysical traverse of the North American plate in Alaska from the active convergent Pacific margin to the passive Arctic margin of the continent.

Alaska↗

USGS2012

A deep research drill hole at the summit of an active volcano, Kilauea, Hawaii

Drilling and geophysical logging data for a 1,262 m‐deep bore hole in the area inferred to overlie the magma reservoir of Kilauea Volcano, Hawaii, support earlier interpretations based on surface geophysical surveys that a zone of brackish or saline water lies above the reservoir. Temperatures encountered within the hole are not sufficiently high to warrant commercial interest; the maximum temperature, 137°C, is at the hole bottom. However, the temperature gradient toward the bottom of the hole (approximately 160 m below sea level) increases sharply to about 370°C/km, perhaps partly reflecting the effect of decreased water circulation as suggested by the geophysical logging data. If this gradient persists or increases with depth, magmatic temperatures would be attained within 3 km from the hole bottom (i.e., approximately 4 km from ground surface)—a depth in accord with data from ground‐deformation and seismic studies.

Hawaii↗

USGS2012

Tectonics and conductivity structures in the Southern Washington Cascades

The tectonic setting of the southern Washington Cascades has been studied with the aid of magnetotelluric (MT) and other geophysical data. The main feature of interest in the geophysical data is a broad high-conductivity anomaly mapped with MT and geomagnetic variation (GMV) data. This anomaly is located roughly within the triangle formed by the volcanoes Mount Rainier, Mount St. Helens, and Mount Adams but exceeds beyond Mount Rainier to the northwest. We interpret the cause of the anomaly to be conductive rocks with resistivities of 1–4 ohm m and thicknesses possibly greater than 15 km. These conductive rocks are found 2–8 km beneath the overlying less conductive volcanic and sedimentary rocks at the surface. Two aeromagnetic lows follow the trend of the conductivity anomaly, and linear belts of strike-slip seismicity are coincident with both these magnetic lows. One of the aeromagnetic lows is coincident with the western margin of the conductivity anomaly. The geophysical data appear to outline a suture zone of probable Eocene age caused by accretion of a large seamount complex (Siletzia) and that may contain large thicknesses of compressed forearc basin and accretionary prism sedimentary rocks of Cretaceous to Eocene age. Part of the shallower conductive rocks may be associated with carbonaceous continental and transitional marine sedimentary rocks of the Puget Group. The contact between the hypothesized compressed basin and the accreted terranes to the west may localize the release of shear stresses in this region of oblique subduction. Several possible explanations for the conductivity anomaly are considered in addition to the compressed basin hypothesis.

Washington↗

USGS2012

A magmatic model of Medicine Lake Volcano, California

Medicine Lake volcano is a Pleistocene and Holocene shield volcano of the southern Cascade Range. It is located behind the main Cascade arc in an extensional tectonic setting where high-alumina basalt is the most commonly erupted lava. This basalt is parental to the higher-silica calc-alkaline and tholeiitic lavas that make up the bulk of the shield. The presence of late Holocene, chemically identical rhyolites on opposite sides of the volcano led to hypotheses of a large shallow silicic magma chamber and of a small, deep chamber that fed rhyolites to the surface via cone sheets. Subsequent geophysical work has been unable to identify a large silicic magma body, and instead a small one has apparently been recognized. Some geologic data support the geophysical results. Tectonic control of vent alignments and the dominance of mafic eruptions both in number of events and volume throughout the history of the volcano indicate that no large silicic magma reservoir exists. Instead, a model is proposed that includes numerous dikes, sills, and small magma bodies, most of which are too small to be recognized by present geophysical methods.

Journal of Geophysical Research Solid Earth↗

USGS2014

Resistivity‐studies of some salt‐water boundaries in the Hawaiian Islands

In the course of a systematic survey of the ground‐water resources of the Hawaiian Islands which is being made under the direction of H. T. Stearns of the United States Geological Survey in cooperation with the Territorial Government of Hawaii, it was found desirable to test the utility of geophysical methods in the solution of certain Hawaiian water‐supply problems. A cooperative geophysical survey for this purpose was undertaken by the Geophysical Section, then with the Bureau of Mines but later transferred to the Geological Survey, during the course of which a series of resistivity‐studies of certain salt‐water boundaries was begun. Since the location of such boundaries is of considerable importance in many areas it has been thought worth while to present the results which have thus far been obtained.

Eos, Transactions, American Geophysical Union↗

USGS2014

Results to be expected from resistivity‐measurements

The work described in this paper was all done in connection with dam‐site investigations and was not directly connected with hydrology. However, geophysics is coming to have a place in hydrologic investigations, and these results may throw some light on what can be accomplished by resistivity‐measurements. We have found that,for many questions not involving exact determinations of depth, resistivity‐ measurements give conclusive answers. Ordinarily a reliable answer can be expected to the question of the existence of a buried channel if the covering is composed of unconsolidated material with a resistivity differing from that of the rock. For example, topographic surveys were made at two alternative dam‐sites on a river about four miles apart. Examination of the surface‐geology indicated that a channel burled under glacial debris possibly existed at each site, but resistivity‐measurements proved that such a channel existed at one site and not at the other. On the other hand, at another site the geologist suspected there might be an old channel on a steep side hill. Geophysical measurements showed a depth of overburden of 46 feet and showed that if a deeper channel exists it must be narrow; but they did not show positively that no such channel exists. Probably a careful survey with a large number of lines would have given a more definite answer, but the rough topography interfered with the resistivity‐work, and time and money were not available for a detailed survey. After completion of the geophysical work, the geologist located some outcrops which led him to conclude that no old channel exists at this site. At two other dam‐sites in Oregon resistivity‐measurements showed that there were no burled channels.

Eos, Transactions, American Geophysical Union↗

USGS2023

Geologic map of Okmok Volcano

The geologic map and description of map units presented here cover approximately 880 km2 of northeastern Umnak Island, Aleutian Islands, Alaska. This report focuses on Okmok Volcano and its eruptive products and updates the mid-20th-century geologic map of Byers (1959). Mapped deposits reflect the state of the volcano just prior to the 2008 eruption. Published information about other portions of Umnak Island geology, including Mount Recheshnoi and Mount Vsevidof, can be found in Byers (1959). The 2008 eruption and its deposits are described in Larsen and others (2009, 2013, 2015). Okmok Volcano is one of 54 historically active volcanoes in the Alaska–Aleutian volcanic arc that stretches across southern mainland Alaska and the Aleutian Islands (fig. 1; Wood and Kienle, 1990; Miller and others, 1998; Cameron and others, 2020). The highest point of the modern Okmok Caldera is along the caldera’s northern rim, 967 m in elevation, and formally named “Mount Okmok” (U.S. Board on Geographic Names, www.usgs.gov/core-science-systems/ngp/boardon-geographic-names/domestic-names). Okmok Volcano dominates the northeastern portion of Umnak Island, which is 100 km southwest of Unalaska/Dutch Harbor and 1,400 km southwest of Anchorage (figs. 1, 2). The Port of Dutch Harbor on Unalaska Island produces the highest volume of seafood for any port in the United States (see fisheries.noaa.gov/resource/document/fisheries-united-states-2018-report). Unalaska city and the Port of Dutch Harbor have been impacted by ash fall and drifting ash clouds from Okmok Volcano’s explosive eruptions as recently as 2008. Holocene and late Pleistocene volcanic rocks and deposits of Okmok Volcano rest upon glaciated Tertiary volcanic and sedimentary rocks (Byers, 1959). The first geologic mapping expedition to Okmok Volcano was by the U.S. Geological Survey (USGS) after the 1945 eruption, largely in response to concerns about volcanic hazards to U.S. military activities in the Aleutians Islands (Byers and others, 1947, 1959; Byers and Brannock, 1949; Byers, 1955, 1959, 1961). The State of Alaska conducted further mapping and geochemical studies as part of its geothermal exploration program in the 1980s (Nye, 1983; Nye and Reid, 1986; Motyka and others, 1993). Additional modern geological work focused on Okmok Volcano and the rest of Umnak Island to address the geochemistry and origin of primary Aleutian arc magmas and subduction zone mass recycling (Marsh, 1982; Brophy and Marsh, 1986; Nye and Reid, 1986; Myers and Marsh, 1987; Miller and others, 1992; Fournelle and others, 1994; Kay and Kay, 1994). In 1998, the Alaska Volcano Observatory (AVO) began a multi-year effort to expand geophysical monitoring in the central Aleutians Islands, including at Okmok Volcano. As part of this effort, AVO geologists from the University of Alaska Fairbanks Geophysical Institute (UAF/GI), the Alaska Division of Geological & Geophysical Surveys (DGGS), and USGS also began a renewed effort to document Okmok Volcano’s recent eruption products. The project started with reconnaissance fieldwork to document and sample products from the 1997 eruption within Okmok Caldera. This evolved into an effort to produce an updated geologic map of Okmok Volcano and gather more information about its eruptive history and hazards. Three significant eruptions occurred at Okmok Volcano in 1958, 1997, and 2008—after fieldwork had been conducted for the original 1:63,360-scale geologic map produced by Byers (1959)—resulting in new volcanic deposits not previously described. Okmok Volcano is one of the most frequently active volcanoes in the Aleutian volcanic arc. Seismic and geodetic monitoring indicate ongoing unrest at Okmok Volcano since at least 1997. Geodetic observations of inflation before and after the 1997 and 2008 eruptions indicate a nearly continuous input of new magma from a depth consistent with frequent eruptions of basalt and basaltic andesite magmas over the past 200 years (Larsen and others, 2013; Lu and others, 2000, 2003, 2005; Mann, 2002; Mann and others, 2002). To better understand the likelihood and character of future eruptions from Okmok Volcano, it is necessary to understand its past behavior, including eruptions since the first geologic map was published by Byers (1959).

Alaska↗

USGS1998

The Holocene sea-level highstand in the equatorial Pacific: Analysis of the insular paleosea-level database

A review of the literature provides 92 estimates of the middle to late Holocene sea-level highstand on Pacific Islands. These data generally support geophysical model calculations that predict a +1 to 3 m relative sea-level highstand on oceanic islands due to the Earth’s rheological response to the melting of the last continental ice sheets and subsequent redistribution of meltwater. Both predictions and observations indicate sea level was higher than present in the equatorial Pacific between 5000 and 1500 y B.P. A non-linear relationship exists between the age and elevation of the highstand peak, suggesting that different rates of isostatic adjustment may occur in the Pacific, with the highest rates of sea-level fall following the highstand near the equator. It is important to resolve detailed sea-level histories from insular sites to test and refine models of climatic, oceanographic, and geophysical processes including hydroisostasy, equatorial ocean siphoning, and lithospheric flexure that are invoked as mechanisms affecting relative sea-level position. We use a select subset of the available database meeting specific criteria to examine model relationships of paleosea-surface topography. This new evaluated database of paleosea-level positions is also validated for testing and constraining geophysical model predictions of past and present sea-level variations.

Coral Reefs↗

USGS2004

Combined interpretation of radar, hydraulic, and tracer data from a fractured-rock aquifer near Mirror Lake, New Hampshire, USA

An integrated interpretation of field experimental cross-hole radar, tracer, and hydraulic data demonstrates the value of combining time-lapse geophysical monitoring with conventional hydrologic measurements for improved characterization of a fractured-rock aquifer. Time-lapse difference-attenuation radar tomography was conducted during saline tracer experiments at the US Geological Survey Fractured Rock Hydrology Research Site near Mirror Lake, Grafton County, New Hampshire, USA. The presence of electrically conductive saline tracer effectively illuminates permeable fractures or pathways for geophysical imaging. The geophysical results guide the construction of three-dimensional numerical models of ground-water flow and solute transport. In an effort to explore alternative explanations for the tracer and tomographic data, a suite of conceptual models involving heterogeneous hydraulic conductivity fields and rate-limited mass transfer are considered. Calibration data include tracer concentrations, the arrival time of peak concentration at the outlet, and steady-state hydraulic head. Results from the coupled inversion procedure suggest that much of the tracer mass migrated outside the three tomographic image planes, and that solute is likely transported by two pathways through the system. This work provides basic and site-specific insights into the control of permeability heterogeneity on ground-water flow and solute transport in fractured rock.

New Hampshire↗

USGS2006

Airborne microorganisms in the African desert dust corridor over the mid-Atlantic ridge, Ocean Drilling Program, Leg 209

The objective of this study was to enhance our understanding of the fate and trans-Atlantic transport of dustborne microorganisms from Northern Africa to the Caribbean and Americas, and more specifically to determine if culturable populations could be detected at a mid-ocean site, closer to the source of dust relative to land-based Caribbean sites, during the early summer months of May and June. Between the dates of 22 May and 30 June 2003, daily air samples were collected and evaluated for the presence of culturable bacterial and fungal colony-forming units (CFU). Here we report a statistically significant correlation between daily atmospheric CFU counts at a mid-ocean research site (???15??N, 45??W) and daily desert dust concentrations as determined by the U.S. Navy's Naval Aerosol Analysis and Prediction System (NAAPS) Global Aerosol Model (Honrath et al. (2004). Journal of Geophysical Research, 109; Johnson et al. (2003). Global Biogeochemical Cycles, 17, 1063; Reid et al. (2004). Geophysical Research Letters, 31; Schollaert, Yoder, Westphal, & O'Reilly (2003). Journal of Geophysical Research, 108, 3191). ?? Springer Science+Business Media B.V. 2006.

Aerobiologia↗

USGS1984

A review of crust and upper mantle structure studies of the Snake River Plain-Yellowstone volcanic system: A major lithospheric anomaly in the western U.S.A.

The Snake River Plain-Yellowstone volcanic system is one of the largest, basaltic, volcanic field in the world. Here, there is clear evidence for northeasterly progression of rhyolitic volcanism with its present position in Yellowstone. Many theories have been advanced for the origin of the Snake River Plain-Yellowstone system. Yellowstone and Eastern Snake River Plain have been studied intensively using various geophysical techniques. Some sparse geophysical data are available for the Western Snake River Plain as well. Teleseismic data show the presence of a large anomalous body with low P- and S-wave velocities in the crust and upper mantle under the Yellowstone caldera. A similar body in which compressional wave velocity is lower than in the surrounding rock is present under the Eastern Snake River Plain. No data on upper mantle anomalies are available for the Western Snake River Plain. Detailed seismic refraction data for the Eastern Snake River Plain show strong lateral heterogeneities and suggest thinning of the granitic crust from below by mafic intrusion. Available data for the Western Snake River Plain also show similar thinning of the upper crust and its replacement by mafic material. The seismic refraction results in Yellowstone show no evidence of the low-velocity anomalies in the lower crust suggested by teleseismic P-delay data and interpreted as due to extensive partial melting. However, the seismic refraction models indicate lower-than-normal velocities and strong lateral inhomogeneities in the upper crust. Particularly obvious in the refraction data are two regions of very low seismic velocities near the Mallard Eake and Sour Creek resurgent domes in the Yellowstone caldera. The low-velocity body near the Sour Creek resurgent dome is interpreted as partially molten rock. Together with other geophysical and thermal data, the seismic results indicate that a sub-lithospheric thermal anomaly is responsible for the time-progressive volcanism along the Eastern Snake River Plain. However, the exact mechanism responsible for the volcanism and details of magma storage and migration are not yet fully understood.

California, Idaho, Montana, Nevada, Oregon, Washin↗