USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geological Research”

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 1,099 records · Page 61Linked to original sources
USGS2001

Geological characterization of remote field sites using visible and infrared spectroscopy: Results from the 1999 Marsokhod field test

Upcoming Mars Surveyor lander missions will include extensive spectroscopic capabilities designed to improve interpretations of the mineralogy and geology of landing sites on Mars. The 1999 Marsokhod Field Experiment (MFE) was a Mars rover simulation designed in part to investigate the utility of visible/near-infrared and thermal infrared field spectrometers to contribute to the remote geological exploration of a Mars analog field site in the California Mojave Desert. The experiment simultaneously investigated the abilities of an off-site science team to effectively analyze and acquire useful imaging and spectroscopic data and to communicate efficiently with rover engineers and an on-site field team to provide meaningful input to rover operations and traverse planning. Experiences gained during the MFE regarding effective communication between different mission operation teams will be useful to upcoming Mars mission teams. Field spectra acquired during the MFE mission exhibited features interpreted at the time as indicative of carbonates (both dolomitic and calcitic), mafic rocks and associated weathering products, and silicic rocks with desert varnish-like coatings. The visible/near-infrared spectra also suggested the presence of organic compounds, including chlorophyll in one rock. Postmission laboratory petrologic and spectral analyses of returned samples confirmed that all rocks identified as carbonates using field measurements alone were calc-silicates and that chlorophyll associated with endolithic organisms was present in the one rock for which it was predicted. Rocks classified from field spectra as silicics and weathered mafics were recognized in the laboratory as metamorphosed monzonites and diorite schists. This discrepancy was likely due to rock coatings sampled by the field spectrometers compared to fresh rock interiors analyzed petrographically, in addition to somewhat different surfaces analyzed by laboratory thermal spectroscopy compared to field spectra.

California↗
USGS2012

Mineralogical studies of sulfide samples and volatile concentrations of basalt glasses from the southern Juan de Fuca Ridge

Sulfide samples obtained from the U.S. Geological Survey's DSRV Alvin dives on the southern Juan de Fuca Ridge closely resemble those from the same area described by Koski et al. (1984). Major minerals include sphalerite, wurtzite, pyrite, marcasite, isocubanite, anhydrite, and chalcopyrite. Equilibrium, if attained at all, during deposition of most sulfides was a transient event over a few tens of micrometers at most and was perturbed by rapid temperature and compositional changes of the circulating fluid. Two new minerals were found: one, a hydrated Zn, Fe hydroxy-chlorosulfate, and the other, a (Mn, Mg, Fe) hydroxide or hydroxy-hydrate. Both were formed at relatively low temperatures. Lizardite, starkeyite, and anatase were found for the first time in such an environment. Sulfide geothermometry involving the system Cu-Fe-S indicates a vent temperature of <328°C for one sample. Fluid inclusion studies on crystals from the same vicinity of the same sample give pressure-corrected homogenization temperatures of 268° and 285°C. Ice-melting temperatures on inclusions from the same sample are about −2.8°C, indicating that the equivalent salinity of the trapped fluid is about 50% greater than that of seawater. Volatile concentrations from vesicle-free basalt glass from the vent field are about 0.013 wt % CO 2 and 0.16 wt % H 2 O. CO 2 contents in these samples yield an entrapment depth of 2200 m of seawater, which is the depth from which the samples were collected.

Journal of Geophysical Research Solid Earth↗
USGS2018

Enhancement of a parsimonious water balance model to simulate surface hydrology in a glacierized watershed

The U.S. Geological Survey monthly water balance model (MWBM) was enhanced with the capability to simulate glaciers in order to make it more suitable for simulating cold region hydrology. The new model, MWBMglacier, is demonstrated in the heavily glacierized and ecologically important Copper River watershed in Southcentral Alaska. Simulated water budget components compared well to satellite‐based observations and ground measurements of streamflow, evapotranspiration, snow extent, and total water storage, with differences ranging from 0.2% to 7% of the precipitation flux. Nash Sutcliffe efficiency for simulated and observed streamflow was greater than 0.8 for six of eight stream gages. Snow extent matched satellite‐based observations with Nash Sutcliffe efficiency values of greater than 0.89 in the four Copper River ecoregions represented. During the simulation period 1949 to 2009, glacier ice melt contributed 25% of total runoff, ranging from 12% to 45% in different tributaries, and glacierized area was reduced by 6%. Statistically significant ( p < 0.05) decreasing and increasing trends in annual glacier mass balance occurred during the multidecade cool and warm phases of the Pacific Decadal Oscillation, respectively, reinforcing the link between climate perturbations and glacier mass balance change. The simulations of glaciers and total runoff for a large, remote region of Alaska provide useful data to evaluate hydrologic, cryospheric, ecologic, and climatic trends. MWBM glacier is a valuable tool to understand when, and to what extent, streamflow may increase or decrease as glaciers respond to a changing climate.

Journal of Geophysical Research F: Earth Surface↗
USGS2022

Magnetotelluric investigations of the Kīlauea Volcano, Hawaii

In 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the U.S. Geological Survey (USGS) Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kīlauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kīlauea and the East and Southwest Rift zones, which has implications for understanding Kīlauea's plumbing system. An improved understanding of the rift zones has implications in understanding large-scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi-station robust processing techniques. In total, data were acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone, where electrical contact resistances are on the order of 100 kΩ. A three-dimensional (3D) MT model study was done to guide interpretation of the observed MT measurements. Synthetic modeling demonstrates that conductive bodies in the upper 3 km can be spatially resolved where MT station sampling is good. Resistivity anomalies in the 3D inversions have a high degree of spatial correlation with previously published seismic velocity anomalies beneath Kīlauea. Melt fractions between 0.096 and 0.117 are calculated for the Kīlauea and Puʻuʻōʻō low resistivity anomalies, respectively.

Hawai'i↗
USGS2022

The influence of drying on the aeolian transport of river-sourced sand

Transgression and regression of water levels (stages) have impacted the evolution of aeolian landforms and sedimentary deposits throughout geologic history. We studied this phenomenon over a five-day period of reduced flow on the Colorado River in Grand Canyon National Park, AZ, USA, in March 2021. These transient low flows exposed river-channel sand deposits to the air, causing progressive desiccation (drying) and thereby making these deposits susceptible to aeolian transport. We measured aeolian threshold friction velocities ( u *t ) for sand saltation and PM10 dust emissions, as well as other characteristics, on a subaerially exposed sandbar and downwind aeolian dunefield during each day of the low river flow. The sandbar transitioned from supply-limited to transport-limited aeolian sediment transport conditions during the regression in river water stage. A possible tipping point between the two transport conditions occurred approximately 48 hours after the drop in river flow. The empirically measured u *t decreased as the sandbar sediment dried with increased subaerial exposure time. Theoretical estimates and empirical measurements of u *t corresponded closely on the aeolian dunefield and on the sandbar when it was drier during the third and fourth day of the experiment. Eighty-seven percent of the variability in u *t was explained by empirical models that provide practical estimates of aeolian transport potential of subaerial river sediment deposits using monitoring data that are commonly available in this and other river systems. The work provides theoretical insight into the response of aeolian processes to sediment supply changes driven by periods of anthropogenic activity, drought, and climate change.

Arizona↗
USGS2000

Contaminant distribution and accumulation in the surface sediments of Long Island Sound

The distribution of contaminants in surface sediments has been measured and mapped as part of a U.S. Geological Survey study of the sediment quality and dynamics of Long Island Sound. Surface samples from 219 stations were analyzed for trace (Ag, Ba, Cd, Cr, Cu, Hg, Ni, Pb, V, Zn and Zr) and major (Al, Fe, Mn, Ca, and Ti) elements, grain size, and Clostridium perfringens spores. Principal Components Analysis was used to identify metals that may covary as a function of common sources or geochemistry. The metallic elements generally have higher concentrations in fine-grained deposits, and their transport and depositional patterns mimic those of small particles. Fine-grained particles are remobilized and transported from areas of high bottom energy and deposited in less dynamic regions of the Sound. Metal concentrations in bottom sediments are high in the western part of the Sound and low in the bottom-scoured regions of the eastern Sound. The sediment chemistry was compared to model results (Signell et al., 1998) and maps of sedimentary environments (Knebel et al., 1999) to better understand the processes responsible for contaminant distribution across the Sound. Metal concentrations were normalized to grain-size and the resulting ratios are uniform in the depositional basins of the Sound and show residual signals in the eastern end as well as in some local areas. The preferential transport of fine-grained material from regions of high bottom stress is probably the dominant factor controlling the metal concentrations in different regions of Long Island Sound. This physical redistribution has implications for environmental management in the region.

Long Island Sound↗
USGS2011

Regional-scale stratigraphy of surface units in Tyrrhena and Iapygia Terrae, Mars: Insights into highland crustal evolution and alteration history

The compositional, thermophysical and geologic characteristics of surface units in Iapygia and Tyrrhena Terra (60°E–100°E, 0°–30°S) provide new insights into the compositional stratigraphy of the region. Intercrater plains are dominated by two surface units. The older unit (unit 1) is deficient in olivine and more degraded and likely consists of a mixture of impact, volcanic and sedimentary materials. The younger unit (unit 2) is enriched in olivine, exhibits a resistant morphology and higher thermal inertia, and likely represents volcanic infilling of plains. Units 1 and 2 bear a strong resemblance to those previously mapped in Mare Serpentis, a section of highlands crust located northwest of Hellas Basin. Thus, the two major intercrater plains units are even more widespread than previously thought and therefore likely constitute important components of Mars' highland stratigraphy. Many craters in the region contain high thermal inertia deposits (unit 3) that are compositionally identical to unit 2. These may have formed via volcanic infilling or may represent sedimentary materials that have been eroded from crater walls and lithified. Less common units include olivine and/or pyroxene‐rich massifs and crater central peaks. These are primarily found within Hellas Basin rim units and may represent mantle materials brought toward the surface during the Hellas impact. Putative chloride deposits are primarily associated with olivine‐deficient surfaces (unit 1) that may be heavily degraded occurrences of unit 2. The observations raise a variety of questions related to Martian crustal evolution and alteration that may have more widespread implications outside the study region.

Journal of Geophysical Research - Planets↗
USGS2011

A Bayesian network to predict coastal vulnerability to sea level rise

Sea level rise during the 21st century will have a wide range of effects on coastal environments, human development, and infrastructure in coastal areas. The broad range of complex factors influencing coastal systems contributes to large uncertainties in predicting long-term sea level rise impacts. Here we explore and demonstrate the capabilities of a Bayesian network (BN) to predict long-term shoreline change associated with sea level rise and make quantitative assessments of prediction uncertainty. A BN is used to define relationships between driving forces, geologic constraints, and coastal response for the U.S. Atlantic coast that include observations of local rates of relative sea level rise, wave height, tide range, geomorphic classification, coastal slope, and shoreline change rate. The BN is used to make probabilistic predictions of shoreline retreat in response to different future sea level rise rates. Results demonstrate that the probability of shoreline retreat increases with higher rates of sea level rise. Where more specific information is included, the probability of shoreline change increases in a number of cases, indicating more confident predictions. A hindcast evaluation of the BN indicates that the network correctly predicts 71% of the cases. Evaluation of the results using Brier skill and log likelihood ratio scores indicates that the network provides shoreline change predictions that are better than the prior probability. Shoreline change outcomes indicating stability (-1 < rate < 1 m/yr) or erosion (rate < -1 m/yr) tend to occur for two sets of input scenarios. Stable shoreline change rates occur mainly for low rates of relative sea level rise and occur in low-vulnerability geomorphic settings. Rates indicating erosion result for cases where the rate of relative sea level rise is high and moderate-to-high vulnerability geomorphic settings occur. In contrast, accretion (rate > 1 m/yr) was not well predicted. We find that BNs can assimilate important factors contributing to coastal change in response to sea level rise and can make quantitative, probabilistic predictions that can be applied to coastal management decisions. Copyright ?? 2011 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2020

Origin and properties of hydrothermal tremor at Lone Star Geyser, Yellowstone National Park, USA

Geysers are rare geologic features that intermittently discharge liquid water and steam driven by heating and decompression boiling. The cause of variability in eruptive styles and the associated seismic signals are not well understood. Data collected from five broadband seismometers at Lone Star Geyser, Yellowstone National Park are used to determine the properties, location, and temporal patterns of hydrothermal tremor. The tremor is harmonic at some stages of the eruption cycle and is caused by near‐periodic repetition of discrete seismic events. Using the polarization of ground motion, we identify the location of tremor sources throughout several eruption cycles. During preplay episodes (smaller eruptions preceding the more vigorous major eruption), tremor occurs at depths of 7–10 m and is laterally offset from the geyser's cone by ~5 m. At the onset of the main eruption, tremor sources migrate laterally and become shallower. As the eruption progresses, tremor sources migrate along the same path but in the opposite direction, ending where preplay tremor originates. The upward and then downward migration of tremor sources during eruptions are consistent with warming of the conduit followed by evacuation of water during the main eruption. We identify systematic relations among the two types of preplays, discharge, and the main eruption. A point‐source moment tensor fit to low‐frequency waveforms of an individual tremor event using half‐space velocity models indicates average V S ≳ 0.8 km/s, source depths ~4–20 m, and moment tensors with primarily positive isotropic and compensated linear vector dipole moments.

Wyoming↗
USGS2004

Influence of near-surface stratigraphy on coastal landslides at Sleeping Bear Dunes National Lakeshore, Lake Michigan, USA

Lake-level change and landslides are primary controls on the development of coastal environments along the coast of northeastern Lake Michigan. The late Quaternary geology of Sleeping Bear Dunes National Lakeshore was examined with high-resolution seismic reflection profiles, ground-penetrating radar (GPR), and boreholes. Based on sequence-stratigraphic principles, this study recognizes ten stratigraphic units and three major unconformities that were formed by late Pleistocene glaciation and postglacial lake-level changes. Locally high sediment supply, and reworking by two regressions and a transgression have produced a complex stratigraphy that is prone to episodic failure. In 1995, a large landslide deposited approximately 1 million m3 of sediment on the lake floor. The highly deformed landslide deposits, up to 18 m thick, extend 3-4 km offshore and unconformably overlie well-stratified glacial and lacustrine sediment. The landslide-prone bluff is underlain by channel-fill deposits that are oriented nearly perpendicular to the shoreline. The paleochannels are at least 10 m deep and 400 m wide and probably represent stream incision during a lake-level lowstand about 10.3 ka B.P. The channels filled with sediment during the subsequent transgression and lake-level highstand, which climaxed about 4.5 ka B.P. As lake level fell from the highstand, the formation of beach ridges and sand dunes sealed off the channel and isolated a small inland lake (Glen Lake), which lies 5 m above the level of Lake Michigan and may be a source of piped groundwater. Our hypothesis is that the paleochannels act as conduits for pore water flow, and thereby locally reduce soil strength and promote slope failure.

Michigan↗
USGS2012

Crop classification modelling using remote sensing and environmental data in the Greater Platte River Basin, USA

With an ever expanding population, potential climate variability and an increasing demand for agriculture-based alternative fuels, accurate agricultural land-cover classification for specific crops and their spatial distributions are becoming critical to researchers, policymakers, land managers and farmers. It is important to ensure the sustainability of these and other land uses and to quantify the net impacts that certain management practices have on the environment. Although other quality crop classification products are often available, temporal and spatial coverage gaps can create complications for certain regional or time-specific applications. Our goal was to develop a model capable of classifying major crops in the Greater Platte River Basin (GPRB) for the post-2000 era to supplement existing crop classification products. This study identifies annual spatial distributions and area totals of corn, soybeans, wheat and other crops across the GPRB from 2000 to 2009. We developed a regression tree classification model based on 2.5 million training data points derived from the National Agricultural Statistics Service (NASS) Cropland Data Layer (CDL) in relation to a variety of other relevant input environmental variables. The primary input variables included the weekly 250 m US Geological Survey Earth Observing System Moderate Resolution Imaging Spectroradiometer normalized differential vegetation index, average long-term growing season temperature, average long-term growing season precipitation and yearly start of growing season. An overall model accuracy rating of 78% was achieved for a test sample of roughly 215 000 independent points that were withheld from model training. Ten 250 m resolution annual crop classification maps were produced and evaluated for the GPRB region, one for each year from 2000 to 2009. In addition to the model accuracy assessment, our validation focused on spatial distribution and county-level crop area totals in comparison with the NASS CDL and county statistics from the US Department of Agriculture (USDA) Census of Agriculture. The results showed that our model produced crop classification maps that closely resembled the spatial distribution trends observed in the NASS CDL and exhibited a close linear agreement with county-by-county crop area totals from USDA census data ( R 2 = 0.90).

Greater Platte River Basin↗
USGS1999

Preliminary results on photometric properties of materials at the Sagan Memorial Station, Mars

Reflectance measurements of selected rocks and soils over a wide range of illumination geometries obtained by the Imager for Mars Pathfinder (IMP) camera provide constraints on interpretations of the physical and mineralogical nature of geologic materials at the landing site. The data sets consist of (1) three small “photometric spot” subframed scenes, covering phase angles from 20° to 150°; (2) two image strips composed of three subframed images each, located along the antisunrise and antisunset lines (photometric equator), covering phase angles from ∼0° to 155°; and (3) full‐image scenes of the rock “Yogi,” covering phase angles from 48° to 100°. Phase functions extracted from calibrated data exhibit a dominantly backscattering photometric function, consistent with the results from the Viking lander cameras. However, forward scattering behavior does appear at phase angles >140°, particularly for the darker gray rock surfaces. Preliminary efforts using a Hapke scattering model are useful in comparing surface properties of different rock and soil types but are not well constrained, possibly due to the incomplete phase angle availability, uncertainties related to the photometric function of the calibration targets, and/or the competing effects of diffuse and direct lighting. Preliminary interpretations of the derived Hapke parameters suggest that (1) red rocks can be modeled as a mixture of gray rocks with a coating of bright and dark soil or dust, and (2) gray rocks have macroscopically smoother surfaces composed of microscopically homogeneous, clear materials with little internal scattering, which may imply a glass‐like or varnished surface.

Journal of Geophysical Research E: Planets↗
USGS2012

Tectonic history of the north portion of the San Andreas fault system, California, inferred from gravity and magnetic anomalies

Geologic and geophysical data for the San Andreas fault system north of San Francisco suggest that the eastern boundary of the Pacific plate migrated eastward from its presumed original position at the base of the continental slope to its present position along the San Andreas transform fault by means of a series of eastward jumps of the Mendocino triple junction. These eastward jumps total a distance of about 150 km since 29 Ma. Correlation of right-laterally displaced gravity and magnetic anomalies that now have components at San Francisco and on the shelf north of Point Arena indicates that the presently active strand of the San Andreas fault north of the San Francisco peninsula formed recently at about 5 Ma when the triple junction jumped eastward a minimum of 100 km to its present location at the north end of the San Andreas fault. Prior to 5 Ma the triple junction was located at the north end of a proposed northwesterly extension of the Pilarcitos fault. This jump took place at what is now about latitude 38°20′N on the North American plate and may have occurred during part of the time when the relative motion between the North American and Pacific plates rotated 20° clockwise, resulting in transpression along the earlier transcurrent San Andreas fault system. The proposed 150-km eastward movement of the triple junction explains the submarine topography near Cape Mendocino where the continental shelf south of the Mendocino fault extends about 130 km farther west than does the shelf directly north of the fault.

Journal of Geophysical Research Solid Earth↗
USGS2005

An efficient algorithm for double-difference tomography and location in heterogeneous media, with an application to the Kilauea volcano

Improving our understanding of crustal processes requires a better knowledge of the geometry and the position of geological bodies. In this study we have designed a method based upon double-difference relocation and tomography to image, as accurately as possible, a heterogeneous medium containing seismogenic objects. Our approach consisted not only of incorporating double difference in tomography but also partly in revisiting tomographic schemes for choosing accurate and stable numerical strategies, adapted to the use of cross-spectral time delays. We used a finite difference solution to the eikonal equation for travel time computation and a Tarantola-Valette approach for both the classical and double-difference three-dimensional tomographic inversion to find accurate earthquake locations and seismic velocity estimates. We estimated efficiently the square root of the inverse model's covariance matrix in the case of a Gaussian correlation function. It allows the use of correlation length and a priori model variance criteria to determine the optimal solution. Double-difference relocation of similar earthquakes is performed in the optimal velocity model, making absolute and relative locations less biased by the velocity model. Double-difference tomography is achieved by using high-accuracy time delay measurements. These algorithms have been applied to earthquake data recorded in the vicinity of Kilauea and Mauna Loa volcanoes for imaging the volcanic structures. Stable and detailed velocity models are obtained: the regional tomography unambiguously highlights the structure of the island of Hawaii and the double-difference tomography shows a detailed image of the southern Kilauea caldera-upper east rift zone magmatic complex. Copyright 2005 by the American Geophysical Union.

Hawaii↗
USGS2013

Doppler weather radar observations of the 2009 eruption of Redoubt Volcano, Alaska

The U.S. Geological Survey (USGS) deployed a transportable Doppler C-band radar during the precursory stage of the 2009 eruption of Redoubt Volcano, Alaska that provided valuable information during subsequent explosive events. We describe the capabilities of this new monitoring tool and present data captured during the Redoubt eruption. The MiniMax 250-C (MM-250C) radar detected seventeen of the nineteen largest explosive events between March 23 and April 4, 2009. Sixteen of these events reached the stratosphere (above 10 km) within 2–5 min of explosion onset. High column and proximal cloud reflectivity values (50 to 60 dBZ) were observed from many of these events, and were likely due to the formation of mm-sized accretionary tephra-ice pellets. Reflectivity data suggest that these pellets formed within the first few minutes of explosion onset. Rapid sedimentation of the mm-sized pellets was observed as a decrease in maximum detection cloud height. The volcanic cloud from the April 4 explosive event showed lower reflectivity values, due to finer particle sizes (related to dome collapse and related pyroclastic flows) and lack of significant pellet formation. Eruption durations determined by the radar were within a factor of two compared to seismic and pressure-sensor derived estimates, and were not well correlated. Ash dispersion observed by the radar was primarily in the upper troposphere below 10 km, but satellite observations indicate the presence of volcanogenic clouds in the stratosphere. This study suggests that radar is a valuable complement to traditional seismic and satellite monitoring of explosive eruptions.

Alaska↗
USGS2019

Heat flow in the Western Arctic Ocean (Amerasian Basin)

From 1963 to 1973 the U.S. Geological Survey (USGS) measured heat flow at 356 sites in the Amerasian Basin (Western Arctic Ocean) from a drifting ice island (T-3). The resulting measurements, which are unevenly distributed on Alpha-Mendeleev Ridge (AMR) and in Canada and Nautilus basins, greatly expand available heat flow data for the Arctic Ocean. Average T-3 heat flow is ~54.7 ± 11.3 mW m-2, and Nautilus Basin, including Mendeleev Plain, is the only well-surveyed area (~13% of data) with significantly higher average heat flow (63.8 mW m-2). Heat flow and bathymetry are not correlated at a large scale, and turbiditic surficial sediments (Canada and Nautilus basins) have higher heat flow than the sediments that blanket the AMR. Thermal gradients are mostly near-linear, implying that conductive heat transport dominates and that near-seafloor sediments are in thermal equilibrium with overlying bottom waters. Combining the heat flow data with modern seismic imagery suggests that some of the observed heat flow variability may be explained by local changes in sediment thickness or lithology or the presence of basement faults that channel circulating seawater. A thermal model that incorporates thermal conductivity variations along a profile from Canada Basin (thick sediment on mostly oceanic crust) to Alpha Ridge (thin sediment over thick magmatic units associated with the High Arctic Large Igneous Province) predicts heat flow lower than that observed on Alpha Ridge. This, along with other observations, implies that circulating fluids modulate conductive heat flow and contribute to high variability in the T-3 dataset. .

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Deep volcanic tremor and magma ascent mechanism under Kilauea, Hawaii

Deep harmonic tremor originating at depths around 40 km under Kilauea was studied using records accumulated since 1962 at the Hawaii Volcano Observatory of the U.S. Geological Survey. The deep source of the tremor was determined by onset times and confirmed by the relative amplitude across the island-wide network of seismometers. The period of tremor was conclusively shown to be determined by the source effect and not by the path or station site effect because the period would change considerably in time but maintained uniformity across the seismic net during the tremor episode. The tremor appeared to be primarily composed of P waves. We interpret the observed period and amplitude in terms of the stationary crack model of Aki et al. (1977) and find that the seismic moment rates for deep tremors are considerably larger than those for shallow-tremors suggesting more vigorous transport for the former. We propose a kinematic source model which may be more appropriate for deep tremor. According to this model, a measurable quantity called ‘reduced displacement’ is directly proportional to the rate of magma flow. A systematic search for deep tremor episodes was made for the period from 1962 through 1979, and the amplitude, period, and duration of the tremor were tabulated. We then constructed a cumulative reduced-displacement plot over the 18-year period. The result shows a generally steady process which does not seem to be significantly affected by major eruptions and large earthquakes near the surface. The total magma flow estimated from the reduced displacement is however, one order of magnitude smaller than that estimated by Swanson (1972). It may be that most channels transport magma aseismically, and only those with strong barriers generate tremor.

Hawaii↗
USGS2025

Changes in phosphorus concentration and flux from 2011 to 2023 in major U.S. tributaries to the Laurentian Great Lakes

Reducing phosphorus (P) flux to the Great Lakes is critical for improving water quality and controlling eutrophication. We used 13 water years (2011–2023) of U.S. Geological Survey data from 24 major U.S. tributaries (representing 47% of the U.S. Great Lakes watershed area) to evaluate temporal changes in orthophosphate (PO 4 -P) and total P (TP) using Weighted Regressions on Time, Discharge, and Season. We assessed actual and flow-normalized P concentrations and fluxes. Between 2011 and 2023, P concentrations and fluxes declined in many tributaries, although the extent and significance of these declines varied. Decreases were more common and statistically likely for TP than PO 4 -P, and several high-loading watersheds had modest or non-significant changes. Flow-normalized PO 4 -P:TP flux ratios increased in over half the tributaries, suggesting that even where P reductions occurred, reductions in the more bioavailable P fraction were proportionally smaller. Actual P fluxes were strongly correlated with streamflow, and year-to-year variability in actual fluxes was, on average, three times greater than variability related to trends in flow-normalized fluxes. This underscores the role of hydrology in modulating P export and highlights how changing precipitation and runoff patterns can obscure or counteract management progress. Spring accounted for the largest share of annual P flux in most tributaries, though many showed declining spring contributions. Our basin-wide analysis reveals that while management efforts may have yielded progress in reducing TP in many watersheds, additional strategies would be needed to address PO 4 -P reductions and account for changing hydrology, especially in high-contributing watersheds.

Journal of Great Lakes Research↗