USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Research of the U.S. Geological Survey”

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 775 records · Page 43Linked to original sources
USGS2006

Metallogeny of the Great Basin: Crustal evolution, fluid flow, and ore deposits

The Great Basin physiographic province in the Western United States contains a diverse assortment of world-class ore deposits. It currently (2006) is the world’s second leading producer of gold, contains large silver and base metal (Cu, Zn, Pb, Mo, W) deposits, a variety of other important metallic (Fe, Ni, Be, REE’s, Hg, PGE) and industrial mineral (diatomite, barite, perlite, kaolinite, gallium) resources, as well as petroleum and geothermal energy resources. Ore deposits are most numerous and largest in size in linear mineral belts with complex geology. U.S. Geological Survey (USGS) scientists are in the final year of a research project initiated in the fall of 2001 to increase understanding of relations between crustal evolution, fluid flow, and ore deposits in the Great Basin. Because of its substantial past and current mineral production, this region has been the focus of numerous investigations over the past century and is the site of ongoing research by industry, academia, and state agencies. A variety of geoinformatic tools was used to organize, reinterpret, and display, in space and time, the large amounts of geologic, geophysical, geochemical, and hydrologic information deemed pertinent to this problem. This information, in combination with concentrated research on (1) critical aspects of the geologic history, (2) an area in northern Nevada that encompasses the major mineral belts, and (3) important mining districts and deposits, is producing new insights about the interplay between key tectonic events, hydrothermal fluid flow, and ore genesis in mineral belts. The results suggest that the Archean to Holocene history of the Great Basin was punctuated by several tectonic events that caused fluids of different origins (sea water, basinal brine, meteoric water, metamorphic water, magmatic water) to move through the crust. Basement faults reactivated during these events localized deformation, sedimentation, magmatism, and hydrothermal fluid flow in overlying rocks to form mineral belts that contain ore deposits of different types and ages that are locally superimposed (demonstrating inheritance). Fluid flow in these systems also was influenced by the distribution of permeable lithologies and paleotopographic highs and lows. Hydrothermal fluids evolved from their initial chemistries towards compositions that reflect the ƒ O 2 and ƒ S 2 buffering capacity of, and the ligands and metals present in, the rocks (±older mineralization) through which they moved. In northern Nevada, where gold deposits are relatively common, carbonaceous, pyritic strata buffered fluids of diverse origins to H 2 S-rich compositions so they could transport gold repeatedly over Paleozoic-Cenozoic time (convergent evolution). Ore formed where metal-laden fluids encountered effective physicochemical traps. Maps of Neogene basin fill and erosion surfaces identify areas where preexisting ore deposits have been progressively exposed or concealed. Comparisons with analogous terrains and deposit types in other parts of the world provide global context. The initial findings and some of the databases, geologic maps, sections, reconstructions, hydrogeologic models, topical syntheses, regional overviews, short courses, field guides, and deposit comparisons produced by project staff and associated managers, contractors, and collaborators have been presented in numerous abstracts, symposia, USGS publications, and professional journals over the last 5 years (see the extensive bibliography). Notable among these was the 2005 Geological Society of Nevada symposium in Reno, Nevada, and the 2005 Geological Society of America annual meeting in Salt Lake City, Utah, where project results were presented to audiences from around the nation and world. The final results of the project will be submitted for publication in 2007 to appropriate USGS and professional journals. A special issue of GEOSPHERE, scheduled for publication in 2007, will be devoted to the results of this project and related work. This special issue will reach an international audience and be available worldwide on the internet. Much of the research for this project has concentrated on areas that will receive the focused attention of the mining industry in the future. As such, the data and interpretations generated by this project have direct use for land-use managers in Federal, State, and local agencies. Improved hydrogeologic models developed by this project will considerably enhance ongoing and future water resource investigations in the region. The increased understanding of when, where, and how hydrothermal systems produce significant economic deposits has direct uses for mineral exploration and for future USGS mineral resource assessments in the Great Basin and other parts of the world.

Great Basin↗
USGS2007

A new towed platform for the unobtrusive surveying of benthic habitats and organisms

Maps of coral ecosystems are needed to support many conservation and management objectives, as well as research activities. Examples include ground-truthing aerial and satellite imagery, characterizing essential habitat, assessing changes, and monitoring the progress of restoration efforts. To address some of these needs, the U.S. Geological Survey developed the Along-Track Reef-Imaging System (ATRIS), a boat-based sensor package for mapping shallow-water benthic environments. ATRIS consists of a digital still camera, a video camera, and an acoustic depth sounder affixed to a moveable pole. This design, however, restricts its deployment to clear waters less than 10 m deep. To overcome this limitation, a towed version has been developed, referred to as Deep ATRIS. The system is based on a light-weight, computer-controlled, towed vehicle that is capable of following a programmed diving profile. The vehicle is 1.3 m long with a 63-cm wing span and can carry a wide variety of research instruments, including CTDs, fluorometers, transmissometers, and cameras. Deep ATRIS is currently equipped with a high-speed (20 frames · s-1) digital camera, custom-built light-emitting-diode lights, a compass, a 3-axis orientation sensor, and a nadir-looking altimeter. The vehicle dynamically adjusts its altitude to maintain a fixed height above the seafloor. The camera has a 29° x 22° field-of-view and captures color images that are 1360 x 1024 pixels in size. GPS coordinates are recorded for each image. A gigabit ethernet connection enables the images to be displayed and archived in real time on the surface computer. Deep ATRIS has a maximum tow speed of 2.6 m · s-1and a theoretical operating tow-depth limit of 27 m. With an improved tow cable, the operating depth can be extended to 90 m. Here, we present results from the initial sea trials in the Gulf of Mexico and Biscayne National Park, Florida, USA, and discuss the utility of Deep ATRIS for map-ping coral reef habitats. Several example mosaics illustrate the high-quality imagery that can be obtained with this system. The images also reveal the potential for unobtrusive animal observations; fish and sea turtles are unperturbed by the presence of Deep ATRIS

Florida↗
USGS2003

Inverse modeling of BTEX dissolution and biodegradation at the Bemidji, MN crude-oil spill site

The U.S. Geological Survey (USGS) solute transport and biodegradation code BIOMOC was used in conjunction with the USGS universal inverse modeling code UCODE to quantify field-scale hydrocarbon dissolution and biodegradation at the USGS Toxic Substances Hydrology Program crude-oil spill research site located near Bemidji, MN. This inverse modeling effort used the extensive historical data compiled at the Bemidji site from 1986 to 1997 and incorporated a multicomponent transport and biodegradation model. Inverse modeling was successful when coupled transport and degradation processes were incorporated into the model and a single dissolution rate coefficient was used for all BTEX components. Assuming a stationary oil body, we simulated benzene, toluene, ethylbenzene, m , p -xylene, and o -xylene (BTEX) concentrations in the oil and ground water, respectively, as well as dissolved oxygen. Dissolution from the oil phase and aerobic and anaerobic degradation processes were represented. The parameters estimated were the recharge rate, hydraulic conductivity, dissolution rate coefficient, individual first-order BTEX anaerobic degradation rates, and transverse dispersivity. Results were similar for simulations obtained using several alternative conceptual models of the hydrologic system and biodegradation processes. The dissolved BTEX concentration data were not sufficient to discriminate between these conceptual models. The calibrated simulations reproduced the general large-scale evolution of the plume, but did not reproduce the observed small-scale spatial and temporal variability in concentrations. The estimated anaerobic biodegradation rates for toluene and o -xylene were greater than the dissolution rate coefficient. However, the estimated anaerobic biodegradation rates for benzene, ethylbenzene, and m , p -xylene were less than the dissolution rate coefficient. The calibrated model was used to determine the BTEX mass balance in the oil body and groundwater plume. Dissolution from the oil body was greatest for compounds with large effective solubilities (benzene) and with large degradation rates (toluene and o -xylene). Anaerobic degradation removed 77% of the BTEX that dissolved into the water phase and aerobic degradation removed 17%. Although goodness-of-fit measures for the alternative conceptual models were not significantly different, predictions made with the models were quite variable.

Minnesota↗
USGS2012

Use of ASTER and MODIS thermal infrared data to quantify heat flow and hydrothermal change at Yellowstone National Park

The overarching aim of this study was to use satellite thermal infrared (TIR) remote sensing to monitor geothermal activity within the Yellowstone geothermal area to meet the missions of both the U.S. Geological Survey and the Yellowstone National Park Geology Program. Specific goals were to: 1) address the challenges of monitoring the surface thermal characteristics of the > 10,000 spatially and temporally dynamic thermal features in the Park (including hot springs, pools, geysers, fumaroles, and mud pots) that are spread out over ~ 5000 km 2 , by using satellite TIR remote sensing tools (e.g., ASTER and MODIS), 2) to estimate the radiant geothermal heat flux (GHF) for Yellowstone's thermal areas, and 3) to identify normal, background thermal changes so that significant, abnormal changes can be recognized, should they ever occur (e.g., changes related to tectonic, hydrothermal, impending volcanic processes, or human activities, such as nearby geothermal development). ASTER TIR data (90-m pixels) were used to estimate the radiant GHF from all of Yellowstone's thermal features and update maps of thermal areas. MODIS TIR data (1-km pixels) were used to record background thermal radiance variations from March 2000 through December 2010 and establish thermal change detection limits. A lower limit for the radiant GHF estimated from ASTER TIR temperature data was established at ~ 2.0 GW, which is ~ 30–45% of the heat flux estimated through geochemical thermometry. Also, about 5 km 2 of thermal areas was added to the geodatabase of mapped thermal areas. A decade-long time-series of MODIS TIR radiance data was dominated by seasonal cycles. A background subtraction technique was used in an attempt to isolate variations due to geothermal changes. Several statistically significant perturbations were noted in the time-series from Norris Geyser Basin, however many of these did not correspond to documented thermal disturbances. This study provides concrete examples of the strengths and limitations of current satellite TIR monitoring of geothermal areas, highlighting some specific areas that can be improved. This work provides a framework for future satellite-based thermal monitoring at Yellowstone and other volcanic and geothermal systems.

Idaho, Montana, Oregon, Wyoming↗
USGS1986

Structural deformation and sedimentation in an active Caldera, Rabaul, Papua New Guinea

Recent seismic and tectonic activity in Rabaul Caldera, Papua New Guinea, suggests that magma is accumulating at a shallow depth beneath this partially submerged structure and that a new volcano may be developing. Changes in onshore elevation since 1971 (as much as 2 m on south Matupit Island) indicate that rapid and large-scale uplifts have occurred on the seafloor near the center of the caldera. The frequency of seismic events within the caldera has also increased during this period. Earthquake locations define an elliptical ring surrounding the center of this uplift within the caldera. A marine geophysical survey in 1982 by the U.S. Geological Survey's R/V "S.P. Lee" in Rabaul Caldera shows the development of a bulge in the seafloor near the center of the caldera. High-resolution seismic reflection profiles show that this bulge consists of two domal uplifts bounded and separated by two major north-south-trending fault zones. Deformed sediments overlie these zones; a prominent slump flanks the area of the bulge. Five major acoustic units were identified in the seismic reflection profiles: an acoustic basement and four sedimentary units consisting of irregularly layered, cross-layered, contorted, and well-layered sequences. The acoustic basement is probably composed of crystalline volcanic rocks, and the layered acoustic units are probably sediments, primarily ash deposited in different environments. The cross-layered, irregularly layered, and contorted units appear to have been deposited in a dynamic environment subjected to strong currents, seismicity, and/or mass wasting, while the well-layered units were deposited in a low-energy environment. Locally, well-layered sequences interfinger with the other sedimentary units, indicating a transitional environment that alternated between high-energy and low-energy depositional processes. A submarine channel cuts most of the acoustic units and appears to be the conduit for sediment transport out of the caldera; it occupies an older buried channel north of the caldera that is presently being exhumed. In the south, active erosion of well-layered sediments is taking place. What are believed to be several young volcanic cones also disrupt the depositional layers. We conclude that the bulge in the seafloor and the associated fault zones are a result of emplacement of magma at a shallow depth. Contorted sediment and slumps adjacent to the bulge are probably the result of uplift and seismic activity. The pattern of seismicity appears to reflect increased magma pressure at depth beneath the caldera floor. This activity may eventually lead to an eruption. ?? 1986.

Journal of Volcanology and Geothermal Research↗
USGS1955

U.S. Geological Survey Trace Elements and related reports through 1954

This report supersedes TEI-400, "U.S. Geological Survey Trace Elements and related reports through 1953." Part I is a numerical list of U.S. Geological Survey Trace Elements Investigations and Memorandum reports. This part lists not only reports (followed by a date) that have been transmitted to the U.S> Atomic Energy Commission, but also reports in preparation (followed by an asterisk) for which tentative titles were available on December 31, 1954. If a report, or portion of a report, has been published, the fact is indicated by the abbreviation of the medium of publication. (See also Part II.) Part II is a reference guide to Trace Elements and related reports that are available to the public (this part supersedes Part I of Bulletin 1019-B). These reports are grouped according to the type of publication or release. If a report, or portion of the report, has been published in more than one form, appropriate cross-references are given to indicate the other forms of publication. Abstracts published in Nuclear Science Abstracts are not included in Part II; however, if the abstract of a Trace Elements report has been published in NSA, the fact is indicated in Part I. Publications in process on December 31, 1954, are designated by an asterisk. Part III is a finding list of states, areas, and subjects and supersedes Part II of Bulletin 1019-B. It is based on information derived mostly from the titles of reports and, where titles are of general nature, from a cursory review of the reports. This list is not a complete index of the information given in Trace Elements and related reports, but is designed to find subjects of major interest, Because of the numerous entries for Colorado and Utah, information has been listed by counties and, where possible, by subject under these states. Other states have county listings only if a county is included in the title of a report; otherwise, area may be listed separately under the state. Major subjects are listed separately in the index and also, where appropriate under states. Analytical methods and subjects related to analytical research are listed under Analytical Methods and Research, but not separately throughout the index. Most mineralogic studies are included under the heading Mineralogy, but are not necessarily listed according to location. Part IV is an author index. The words “with” and “and” are used to indicate seniority of authorship. For example, a listing of Smith and Brown indicated that Smith is the senior author. A listing of Smith with Brown indicates that Brown is the senior author. In both parts III and IV all Trace Elements reports are listed, as well as other related reports that have not been issued as Trace Elements reports. The following abbreviations have been used: TEI, Trace Elements Investigations report; TEM, Trace Elements Memorandum report; OF, open file; TIS, Technical Information Service release; C, Circular; B, Bulletin; P, Professional Paper; MC, Maps and Charts; J, Journal; and NSA, Nuclear Science Abstracts.

Trace Elements Investigations↗
USGS1997

Thermal regime of the San Andreas fault near Parkfield, California

Knowledge of the temperature variation with depth near the San Andreas fault is vital to understanding the physical processes that occur within the fault zone during earthquakes and creep events. Parkfield is near the southern end of the Coast Ranges segment of the San Andreas fault. This segment has higher mean heat flow than the Cape Mendocino segment to the northwest or the Mojave segment to the southeast. Boreholes were drilled specifically for the U.S. Geological Survey's Parkfield earthquake prediction experiment or converted from other uses at 25 sites within a few kilometers of the fault near Parkfield. These holes, which range in depth from 150 to over 1500m, were intended mainly for the deployment of volumetric strain meters, water-level recorders, and other downhole instruments. Temperature profiles were obtained from all the holes, and heat flow values were estimated from 17 of them. For a number of reasons, including a paucity of thermal conductivity data and rugged local topography, the accuracy of individual determinations was not sufficiently high to document local variations in heat flow. Values range from 54 to 92 mW m −2 , with mean and 95% confidence limits of 74±4 mW m −2 . This mean is slightly lower than the mean (83±3) of 39 previously published values from the central Coast Ranges, but it is consistent with the overall pattern of elevated heat flow in the Coast Ranges, and it is transitional to the mean of 68±2 mW m −2 that characterizes the Mojave segment of the San Andreas fault immediately to the south. The lack of a heat flow peak near the fault underscores the absence of a frictional thermal anomaly and provides additional support for a very small resolved shear stress parallel to the San Andreas fault and the nearly fault-normal maximum compressive stress observed in this region. Estimates of subsurface thermal conditions indicate that the seismicaseismic transition for the Parkfield segment corresponds to temperatures in the range of 350°–400°C. Increasing heat flow to the northwest of Parkfield corresponds to a transition from locked to creeping sections and to a shallowing of the base of seismicity and confirms the importance of temperature in controlling the thickness of the seismogenic crust. Lateral variations in heat flow do not appear to have any major role in determining the regularity of M 5.5–6 earthquakes at Parkfield.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Effects of tow transit on the efficacy of the Chicago Sanitary and Ship Canal Electric Dispersal Barrier System

In 2016, the U.S. Fish and Wildlife Service, U.S. Geological Survey, and U.S. Army Corps of Engineers undertook a field study in the Chicago Sanitary and Ship Canal near Romeoville, Illinois to determine the influence of tow transit on the efficacy of the Electric Dispersal Barrier System (EDBS) in preventing the passage of juvenile fish (total length < 100 millimeters (mm)). Dual-frequency identification sonar data showed that large schools of juvenile fish (mean school size of 120 fish; n = 19) moved upstream and crossed the electric field of an array in the EDBS concurrent with downstream-bound (downbound) loaded tows in 89.5% of trials. Smaller schools of juvenile fish (mean school size of 98 fish; n = 15) moved downstream and crossed the electric field of an array in the EDBS concurrent with upstream-bound (upbound) loaded tows in 73.3% of trials. Observed fish passages through the EDBS were always opposite to the direction of tow movement, and not associated with propeller wash. These schools were not observed to breach the EDBS in the absence of a tow and showed no signs of incapacitation in the barrier during tow passage. Loaded tows transiting the EDBS create a return current of water flowing between the tow and the canal wall that typically travels opposite the direction of tow movement, and cause a decrease in the voltage gradient of the barrier of up to 88%. Return currents and decreases in voltage gradients induced by tow passage likely contributed to the observed fish passage through the EDBS. The efficacy of the EDBS in preventing the passage of small, wild fish is compromised while tows are moving across the barrier system. In particular, downbound tows moving through the EDBS create a pathway for the upstream movement of small fish, and therefore may increase the risk of transfer of invasive fishes from the Mississippi River Basin to the Great Lakes Basin.

Illinois↗
USGS2012

The crustal structure of the Wrangellia Terrane along the East Glenn Highway, eastern‐southern Alaska

Recently acquired seismic refraction data from eastern‐southern Alaska provide new information on the structure and composition of the Wrangellia and adjacent terranes. The data comprise a 160‐km‐long refraction profile along the East Glenn (Tok‐Cutoff) Highway that was collected as part of the U.S. Geological Survey's multidisciplinary Trans‐Alaska Crustal Transect program. The upper 3 km of the Wrangellia terrane and associated rocks is characterized by low compressional wave velocities (V p = 1.9, 3.3, 4.6, 5.6 km s −1 ) and high‐velocity gradients common to most onshore seismic refraction profiles. There is also clear seismic expression of the West Fork fault system as a steep, down‐to‐the‐southwest fault that separates the Peninsular terrane to the southwest and the metamorphic complex of Gulkana River to the northeast. In contrast, no seismic expression occurs for the Paxson Lake fault, which separates the Wrangellia terrane from the metamorphic complex of Gulkana River. Adjacent to the Denali fault, within the Wrangellia terrane, two high‐velocity bodies (V p = 6.6 km s −1 ) occur in the upper crust. One of these extends to ∼10‐km depth and correlates with a late Paleozoic dioritic complex, suggesting that the Wrangellia terrane is at least 10 km thick in this part of Alaska. From 5 to 23 km depth, the crust appears seismically homogeneous, with velocity increasing from V p = 6.2 to V p = 6.6 km s −1 . Beneath this level, the crust is less well resolved, although evidence exists for a low‐velocity zone between 23 and 26 km and a possible southwest dipping interface at 35 km. No identifiable mantle refraction or reflection is observed, possibly indicating a crust as thick as 55 km. The relatively low seismic velocities in the upper 23 km of the crust compare favorably with laboratory‐measured velocities on pelitic schists and intermediate‐composition plutonic rocks (granites and granodiorites), both of which are recognized in Wrangellia. We interpret the seismic velocities to indicate that silicic‐to‐intermediate‐composition rocks are important constituents of the basement of this part of Wrangellia. Geologic evidence indicates that the Alaskan part of the Wrangellia terrane is a Paleozoic and Mesozoic island arc: our seismic evidence indicates it may have been built mostly on continental crust as opposed to the fragment of Wrangellia from Vancouver Island which was probably built on oceanic crust.

Alaska↗
USGS2012

Seismic constraints on the nature of lower crustal reflectors beneath the extending Southern Transition Zone of the Colorado Plateau, Arizona

We determine the reflection polarity and exploit variations in P and S wave reflectivity and P wave amplitude versus offset (AVO) to constrain the origin of lower crustal reflectivity observed on new three-component seismic data recorded across the structural transition of the Colorado Plateau. The near vertical incidence reflection data were collected by Stanford University in 1989 as part of the U.S. Geological Survey Pacific to Arizona Crustal Experiment that traversed the Arizona Transition Zone of the Colorado Plateau. The results of independent waveform modeling methods are consistent with much of the lower crustal reflectivity resulting from thin, high-impedance layers. The reflection polarity of the cleanest lower crustal events is positive, which implies that these reflections result from high-velocity contrasts, and the waveform character indicates that the reflectors are probably layers less than or approximately equal to 200 m thick. The lower crustal events are generally less reflective to incident S waves than to P waves, which agrees with the predicted behavior of high-velocity mafic layering. Analysis of the P wave AVO character of lower crustal reflections demonstrates that the events maintain a constant amplitude with offset, which is most consistent with a mafic-layering model. One exception is a high-amplitude (10 dB above background) event near the base of lower crustal reflectivity which abruptly decreases in amplitude at increasing offsets. The event has a pronounced S wave response, which along with its negative AVO trend is a possible indication of the presence of fluids in the lower crust. The Arizona Transition Zone is an active but weakly extended province, which causes us to discard models of lower crustal layering resulting from shearing because of the high degree of strain required to create such layers. Instead, we favor horizontal basaltic intrusions as the primary origin of high-impedance reflectors based on (1) The fact that most xenoliths in eruptive basalts of the Transition Zone are of mafic igneous composition, (2) indications that a pulse of magmatic activity crossed the Transition Zone in the late Tertiary period, and (3) the high regional heat flow observed in the Transition Zone. The apparent presence of fluids near the base of the reflective zone may indicate a partially molten intrusion. We present a mechanism by which magma can be trapped and be induced to intrude horizontally at rheologic contrasts in extending crust.

Arizona↗
USGS2002

U.S. Geological Survey spatial data access

The U.S. Geological Survey (USGS) has done a progress review on improving access to its spatial data holdings over the Web. The USGS EROS Data Center has created three major Web-based interfaces to deliver spatial data to the general public; they are Earth Explorer, the Seamless Data Distribution System (SDDS), and the USGS Web Mapping Portal. Lessons were learned in developing these systems, and various resources were needed for their implementation. The USGS serves as a fact-finding agency in the U.S. Government that collects, monitors, analyzes, and provides scientific information about natural resource conditions and issues. To carry out its mission, the USGS has created and managed spatial data since its inception. Originally relying on paper maps, the USGS now uses advanced technology to produce digital representations of the Earth’s features. The spatial products of the USGS include both source and derivative data. Derivative datasets include Digital Orthophoto Quadrangles (DOQ), Digital Elevation Models, Digital Line Graphs, land-cover Digital Raster Graphics, and the seamless National Elevation Dataset. These products, created with automated processes, use aerial photographs, satellite images, or other cartographic information such as scanned paper maps as source data. With Earth Explorer, users can search multiple inventories through metadata queries and can browse satellite and DOQ imagery. They can place orders and make payment through secure credit card transactions. Some USGS spatial data can be accessed with SDDS. The SDDS uses an ArcIMS map service interface to identify the user’s areas of interest and determine the output format; it allows the user to either download the actual spatial data directly for small areas or place orders for larger areas to be delivered on media. The USGS Web Mapping Portal provides views of national and international datasets through an ArcIMS map service interface. In addition, the map portal posts news about new map services available from the USGS, many simultaneously published on the Environmental Systems Research Institute Geography Network. These three information systems use new software tools and expanded hardware to meet the requirements of the users. The systems are designed to handle the required workload and are relatively easy to enhance and maintain. The software tools give users a high level of functionality and help the system conform to industry standards. The hardware and software architecture is designed to handle the large amounts of spatial data and Internet traffic required by the information systems. Last, customer support was needed to answer questions, monitor e-mail, and report customer problems.

Journal of Geospatial Engineering↗
USGS2010

Analysis of nonvolcanic tremor on the San Andreas Fault near Parkfield, CA using U.S. Geological Survey Parkfield Seismic Array

Reports by Nadeau and Dolenc (2005) that tremor had been detected near Cholame Valley spawned an effort to use UPSAR (U. S. Geological Survey Parkfield Seismic Array) to study characteristics of tremor. UPSAR was modified to record three channels of velocity at 40–50 sps continuously in January 2005 and ran for about 1 month, during which time we recorded numerous episodes of tremor. One tremor, on 21 January at 0728, was recorded with particularly high signal levels as well as another episode 3 days later. Both events were very emergent, had a frequency content between 2 and 8 Hz, and had numerous high-amplitude, short-duration arrivals within the tremor signal. Here using the first episode as an example, we discuss an analysis procedure, which yields azimuth and apparent velocity of the tremor at UPSAR. We then provide locations for both tremor episodes. The emphasis here is how the tremor episode evolves. Twelve stations were operating at the time of recording. Slowness of arrivals was determined using cross correlation of pairs of stations; the same method used in analyzing the main shock data from 28 September 2004. A feature of this analysis is that 20 s of the time series were used at a time to calculate correlation; the longer windows resulted in more consistent estimates of slowness, but lower peak correlations. These values of correlation (peaks of about 0.25), however, are similar to that obtained for the S wave of a microearthquake. Observed peaks in slowness were traced back to source locations assumed to lie on the San Andreas fault. Our inferred locations for the two tremor events cluster near the locations of previously observed tremor, south of the Cholame Valley. Tremor source depths are in the 14–24 km range, which is below the seismogenic brittle zone, but above the Moho. Estimates of error do not preclude locations below the Moho, however. The tremor signal is very emergent but contains packets that are several times larger than the background tremor signal and lasts about 5 s. These impulsive wavelets are similar to low-frequency earthquakes signals seen in Japan but appear to be broader band rather than just higher in low-frequency energy. They may be more appropriately called high-energy tremor (HET). HET signals at UPSAR correlate well with the record of this event from station GHIB of the HRSN borehole array at Parkfield and HETs typically have a higher cross-correlation coefficient than the rest of the tremor event. The amplitudes of a large HET are consistent with a magnitude of 0.1 when compared with a M2.3 event that had about the same epicenter. Polarizations of the tremor episode at UPSAR are mostly just north of east. Both linearity and azimuth evolve over time suggesting a change in tremor source location over time and linearity is typically higher at the HETs.

California↗
USGS2012

Teleseismic tomography of the compressional wave velocity structure beneath the Long Valley region, California

In 1982 and 1984 the U.S. Geological Survey used several seismic networks, totaling over 90 stations, to record teleseismic P waves and measure travel time residuals in an area centered on the Long Valley caldera. We inverted the travel time residuals to obtain a three-dimensional image of the velocity structure with resolution of 5–6 km to depths of 70 km beneath the array. Direct inversion of these data indicates that the 2- to 4-km-thick low-velocity caldera fill contaminates the signal from any midcrustal velocity anomalies beneath the caldera. Thus two methods were used to strip the effects of the upper crust from the travel time residuals: (1) ray tracing through upper crustal velocity models provided by seismic refraction experiments and gravity surveys, and (2) an iterative stripping scheme using the inversion itself. The methods produce essentially identical results and adequately remove the effects of the shallowest crustal structures, including the caldera fill and hydrothermal alteration effects. The resulting “stripped” models show two well-resolved midcrustal low-velocity bodies in the Long Valley region. The first body is centered between 7 and 20 km depth beneath the resurgent dome of the Long Valley caldera and has a volume of 150–600 km 3 . The second, with a similar volume, is centered between 10 and 20 km depth beneath the Mono Craters, about 10 km north of Long Valley. Velocity contrasts in both of these bodies are about 6–10%, and the features are interpreted as silicic magma chambers. This experiment does not preclude the presence of additional pockets of magma smaller than 5 km across in the upper crust, particularly beneath the resurgent dome of the caldera (which would be removed with the stripping methods). It is likely that the midcrust anomaly beneath the resurgent dome is a remnant of the caldera forming magma chamber. Upper mantle velocities are lower than average beneath Mono Craters and higher than average beneath Long Valley. We tentatively interpret these anomalies as mantle partial melt beneath the Mono Craters and refractory upper mantle residuum beneath Long Valley. The high eruptive rate of the Mono Craters and these upper mantle structures suggest that the focus of volcanism is shifting north from Long Valley to the Mono Craters.

Journal of Geophysical Research Solid Earth↗
USGS2026

Petrologic and geochemical characterization of the first year of the 2024–ongoing episodic Halemaumau eruption of Kilauea (Episodes 1–39)

Kīlauea has been erupting episodically since December 23, 2024 (and is ongoing at the time of this writing in 2026) with two consistently active vents in Halema‘uma‘u producing lava fountains ranging in height from 40 to 540 m. In the first year of the eruption (through December 24, 2025, episode 39), the U.S. Geological Survey's Hawaiian Volcano Observatory collected 178 samples for eruption monitoring purposes and used a subset of 80 for near-real-time geochemical analyses. Tephra, lava, and recycled lava lake ooze-out samples have been studied on inter-episode timescales using diverse microanalytical techniques, where research-scale investigations have proven to be advantageous pseudo-near-real-time techniques. Here we present a geochemical and petrologic overview of episodes 1–39 to better understand eruption dynamics and identify changes at Kīlauea following the 2018 summit caldera collapse. There are no significant intra-episode geochemical changes that might be related to variable fountain heights or episode durations, suggesting that eruption dynamics are instead primarily controlled by viscosity and conduit geometry changes throughout a single episode and between episodes. There is variable but clear influence of the still-molten lava lake interior on compositions of this eruption, including in tephra samples from lava fountains, indicating that the vent conduits are connected in some way at depth with the still molten passive lava lake core. Tephra glasses (7.0–7.8 wt.% MgO) show cyclic changes in MgO that were decoupled from lava fountain heights and behavior in episodes 1–28 but positively correlated after episode 29. Additionally, whole-rock and glass compositions are enriched in minor elements such as K 2 O and TiO 2 compared to 2008–2023 summit eruptions, tracking the arrival of a new magma composition that may be related to the prolonged episodic eruptive behavior. Olivine-hosted melt inclusions from five episodes have ~0.35–0.43 wt.% H2O, < 300 ppm CO 2 and < 1500 ppm S, with no significant variations between episodes. Melt and fluid inclusion barometry show that magmas driving these episodes are stored and sourced from the shallow 1–2 km Halemaʻumaʻu reservoir, as were olivine from the 2023 summit eruption. The ongoing episodic eruption also has higher bulk-rock MgO compared to the five previous summit eruptions (2020−2023) and the forsterite (Fo) content of olivine has varied systematically from 2015 to 2025 with the equilibrium olivine composition of the Halema‘uma‘u reservoir increasing from Fo 80.7 in 2015 to Fo 82.7 by 2025. We argue that the Halemaʻumaʻu reservoir is heating up because the rate of magma supply is outpacing effusion and cooling following caldera collapse in 2018.

Hawaii↗
USGS2012

Petrology and hydrothermal mineralogy of US Geological Survey Newberry: 2. Drill core from Newberry Caldera, Oregon

U.S. Geological Survey Newberry 2 was drilled to a depth of 932 m within Newberry caldera. The bottom-hole temperature of 265°C is the highest reported temperature of any drill hole in the Cascades region of the United States. The upper part of the stratigraphic section penetrated by Newberry 2 consists of caldera fill below which are increasingly more mafic lavas ranging from rhyodacite at 501 m to basalt at 932 m. Measured temperatures shallower than 300 m are less than 35°C, and rock alteration consists of hydration of glass and local palagonitization of basaltic tuffs. Incipient zeolitization and partial smectite replacement of ash and pumice occurred throughout the pumiceous lithic tuffs from 300 to 500 m. Higher-temperature alteration of the tuffs to chlorite and mordenite occurs adjacent to a rhyodacite sill at 460–470 m; alteration minerals within the sill consist of pyrrhotite, pyrite, quartz, calcite, and siderite. Below 697 m the rocks are progressively more altered with depth mainly because of increased temperature along a conductive gradient from 100°C at 697 m to 265°C at 930 m. Fluid inclusions in quartz and calcite indicate that temperatures in the past have been higher than at present, most likely due to local confining pressures between impermeable lava flows. Flow breccias are more altered than the adjacent dense massive lava flows, regardless of composition, because of their much higher permeability. Hydrothermal minerals in this zone are mainly mixed-layer chlorite-smectite, quartz, calcite, and pyrite. Chlorite becomes more abundant than mixed-layer clays near the bottom of the hole. In the lowest two lava flows, epidote, anhydrite, and scarce hematite occur locally. Alteration and leaching in the basal 2 m are unique and have led to the postulation of a localized two-phase fluid zone consisting mainly of steam and CO 2 . The hydrothermal system of Newberry 2 is a simple evolving system associated with the evolution of Newberry Volcano. Only a few localized highly altered intervals where fracturing controlled fluid access occur in the core. There are no crosscutting fractures to indicate multiple hydrothermal systems. Chemical analyses of altered rocks and equivalent fresh rocks indicate that little chemical migration has taken place.

Journal of Geophysical Research Solid Earth↗
USGS2023

Monitoring sediment transport pathways from an artificial nearshore berm, South Padre Island, Texas, USA, August 2018 to November 2019: Implications for coastal management

During August 2018 – November 2019, the transport pathways of dredge material from a specially constructed nearshore feeder berm were investigated as part of a collaborative study by the City of South Padre Island, U.S. Army Corps of Engineers–Galveston District, U.S. Geological Survey, Partrac GeoMarine Inc., and Texas A&M University, into the efficacy of beneficial use dredge material (BUDM) as a method of replenishing the beach profile and shoreface at South Padre Island, Texas with sediment. Dual-signature (fluorescent and ferrimagnetic) tracer particles, designed to be hydraulically equivalent to the dredge material, were placed on the berm, and a sampling program was initiated to monitor the spatiotemporal movement of tracer particles under the influence of the prevailing hydrodynamic regime. Wave and current data were collected and were used together with available metocean data to identify the forcing mechanisms of sediment transport; improved understanding garnered from the consideration of multiple datasets can be used to inform future coastal management decisions. Tracer analysis results indicated low magnitude, shoreward transport of dredge material from the berm and along shore transport in both northerly and southerly directions. Small amounts of tracer detected in beach face samples demonstrated connectivity between the constructed feeder berm and the shoreface. However, the generally low tracer concentration within beach face samples indicated low rates of sediment transport for berm sediments, and a low magnitude sediment transport pathway from the berm directly to the beach face, with possible storage of sediment in the nearshore bar system. Given that the berm was at or beyond the closure depth and considering the relatively weak prevailing near bottom hydrodynamic conditions, the potential for sediment to be mobilized and transported during the study period was generally low. Periods of higher energy waves, driven by north-northwest winds, were identified as a key driver of sediment transport but due to the infrequent nature of these events the sediment transport regime across the area of interest was temporally limited. Longshore movement of sediment was both north and south, mainly dependent on prevailing wind directions and resulting longshore currents. We discuss implications for coastal zone management and the use of feeder berms as an artificial beach nourishment mechanism to replenish sediments lost by coastal erosional processes. The insight into nearshore emplacement of material can be used by local governments and coastal managers to optimize the efficacy of berm emplacement and implement such schemes as a cost-effective nourishment option, or when onshore placement of material is not feasible.

Texas↗
USGS2003

Deformation near the Casa Diablo geothermal well field and related processes Long Valley caldera, Eastern California, 1993-2000

Regional first-order leveling lines, which extend from Lee Vining, CA, to Tom's Place, CA, have been surveyed periodically since 1957 by the U.S. Geological Survey (USGS), the National Geodetic Survey (NGS), and Caltrans. Two of the regional survey lines, or leveling networks, intersect at the Casa Diablo geothermal well field. These leveling networks, referenced to a distant bench mark (C916) near Lee Vining, provide time-series vertical control data of land-surface deformation that began around 1980. These data are also useful for delineating localized subsidence at Casa Diablo related to reservoir pressure and temperature changes owing to geothermal development that began in 1985. A comparison of differences in bench-mark elevations for five time periods between 1983 and 1997 shows the development and expansion of a subsidence bowl at Casa Diablo. The subsidence coincides spatially with the geothermal well field and temporally with the increased production rates and the deepening of injection wells in 1991, which resulted in an increase in the rate of pressure decline. The subsidence, superimposed on a broad area of uplift, totaled about 310 mm by 1997. The USGS established orthogonal tilt arrays in 1983 to better monitor deformation across the caldera. One tilt array (DBR) was established near what would later become the Casa Diablo geothermal well field. This array responded to magmatic intrusions prior to geothermal development, tilting away from the well field. With the start of geothermal fluid extraction in 1985, tilt at the DBR array reversed direction and began tilting into the well field. In 1991, geothermal power production was increased by a factor of four, and reservoir pressures began a period of steep decline. These changes caused a temporary three-fold increase in the tilt rate. The tilt rate became stable in 1993 and was about 40% lower than that measured in 1991-1992, but still greater than the rates measured during 1985-1990. Data from the local leveling networks spanning the well field and the bounding graben were analyzed for several 2-year periods (1993-1995, 1995-1997, and 1997-1999). Annual rates of change across the normal faults bounding the graben have steadily decreased for each 2-year period between 1993 and 1999, reflecting the slowing decline in geothermal reservoir pressure. Horizontal control data from a two-color electronic distance meter (EDM) defined the lateral extent of subsidence at Casa Diablo. The EDM and leveling data elucidate the localized effect of the shallow source of subsidence and the broader effect of the deeper magmatic inflation source. Data from bench marks common to both the vertical and the horizontal control networks were used to assess the effect of subsidence on the EDM base station (CASA). Modeling of geodetic data collected during periods of little or no magmatic inflation indicated that the CASA two-color EDM station is being drawn toward the well field at a rate of 3-5 mm/yr. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2022

Tungsten skarn potential of the Yukon-Tanana Upland, eastern Alaska, USA—A mineral resource assessment

Tungsten (W) is used in a variety of industrial and technological applications and has been identified as a critical mineral for the United States, India, the European Union, and other countries. These countries rely on W imports mostly from China, which leaves them vulnerable to supply disruption. Consequently, the U.S. government has a current initiative to understand domestic resource potential. The eastern Alaska portion of the Yukon-Tanana Upland (YTU), is prospective for W skarn deposits, the major source of global W supply. The regional geology consists of juxtaposed Paleozoic lithotectonic packages that were reaccreted to North America in the Mesozoic . Multiple subsequent episodes of arc-related magmatism intruded the lithotectonic packages, accompanied by W skarn formation mostly associated with 100–90 Ma intrusions; major W skarn deposits in Canada are part of the same metallogenic event (e.g., Mactung, Cantung). In this paper, we present an assessment for undiscovered W skarn resources for parts of the lesser-explored western (Alaskan) portion of the YTU. We used GIS proximity analysis to map the intersection of pluton and carbonate-bearing rocks to define three permissive tracts for W skarn deposits. The permissive tracts were qualitatively assessed by mineral potential mapping using region-wide sediment geochemistry and mineral concentrate datasets. This analysis showed that much of the western YTU has high potential for undiscovered W skarn deposits, whereas the eastern and southern YTU had only isolated areas of medium to high potential. Historical production and the quality of the geochemistry data of the western YTU tract (ca. 9200 km 2 ) permitted a quantitative assessment of undiscovered W resources. Probabilistic estimates by a panel of 20 experts predicted a 70% chance of one to three undiscovered W skarn deposits in the western YTU tract. The rationale for favorability employed by the expert panel included favorable lithology, previous production, clustering of previously mined deposits, W placers in the area, lack of recent exploration, pan concentrates containing W minerals, and W geochemical anomalies. Estimates were combined with a global grade and tonnage model for W skarns in a Monte Carlo simulation and provided a median estimate of undiscovered resources of 94 kt WO 3 . If the undiscovered W skarn deposits are located close to infrastructure (e.g., near Fairbanks, or close to roads and/or power grid), application of an economic filter indicates that the median total economically recoverable WO 3 is 63 kt with a net present value (NPV) of $330 million USD (2008 dollars). Whereas if deposits are far from infrastructure, median recoverable WO 3 is only 30 kt and the NPV is $44 million. Our models for contained WO 3 resources and NPV estimates for the western YTU tract are considerably lower than the known resources in skarns in adjacent areas in Canada. Estimates for the western YTU are also lower than preliminary estimates for undiscovered W skarn deposits in areas of the western conterminous United States. We speculate that lower permeability and continuity of favorable carbonate rock horizons in the relatively higher-grade metamorphic country rocks in the Alaska portion of the YTU may explain some of the differences in prospectivity. More detailed geologic mapping, modern geochemistry, and geophysical surveys are needed to refine the resource potential of the whole YTU. Regardless, quantitative mineral resource assessment provides a useful tool for making first-order regional estimates of undiscovered resources, identifying target areas for new data acquisition, and guiding research on the fundamental controls of district-scale metallogenic endowments.

Alaska↗