USGS Science⌕ Search

SEARCH · USGS Science

Results for “The R Journal”

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

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

At least 1,333 records · Page 74Linked to original sources

Factors associated with the deposition of Cladophora on Lake Michigan beaches in 2012

Deposition of the macroalgae Cladophora spp. was monitored on 18 beaches around Lake Michigan during 2012 at a high temporal frequency. We observed a high degree of spatial variability in Cladophora deposition among beaches on Lake Michigan, even within local regions, with no clear regional pattern in the intensity of Cladophora deposition. A strong seasonal pattern in Cladophora deposition was observed, with the heaviest deposition occurring during mid-summer. Several beaches exhibited high temporal variability in Cladophora deposition over short time scales, suggesting that drifting algal mats may be extremely dynamic in nearshore environments of the Great Lakes. Cladophora deposition on Lake Michigan beaches was primarily related to the presence of nearshore structures, local population density, and nearshore bathymetry. There was relatively little evidence that waves, winds, or currents were associated with Cladophora deposition on beaches, but this may be due to the relatively poor resolution of existing nearshore hydrodynamic data. Developing a predictive understanding of beach-cast Cladophora dynamics in Great Lakes environments may require both intensive Cladophora monitoring and fine-scale local hydrodynamic modeling efforts.

Illinois, Indiana, Michigan, Wisconsin↗

Use of an artificial stream to monitor avoidance behavior of larval sea lamprey in response to TFM and niclosamide

The lampricide 3-trifluoromethyl-4-nitrophenol (TFM) has been used in liquid form to control larval sea lamprey ( Petromyzon marinus ) in Great Lakes tributaries since the late 1950s. In the 1980s a dissolvable TFM bar was developed as a supplemental tool for application to small tributaries as a deterrent to larvae seeking water not activated with TFM. The size, mass, and number of bars needed in some streams, as well as the location of the streams, limit the utility of a TFM bar. The development and use of an alternative niclosamide bar has the potential to use fewer bars to achieve similar results. However, the use of a niclosamide bar is dependent upon its larval deterrent capability compared to the TFM bar. In this study, we developed a laboratory-scale, simulated stream fluvarium with several avoidance areas including two side channels and a seep. The objective was to evaluate the deterrent capabilities of TFM and niclosamide. We found similar behavioral responses, with TFM and niclosamide having similar capabilities to prevent sea lamprey from seeking refuge in side channels and seep avoidance areas. TFM-treated side channels and seep increased sea lamprey occupancy in the main channel 2.56 times more than the untreated-controls (95% CI 1.63–4.14) whereas niclosamide-treated side channels and seep increased sea lamprey occupancy of the main channel 2.68 times more than the untreated-controls (95% CI 1.72–4.32). These responses indicate a niclosamide bar would effectively prevent sea lamprey escapement into freshwater during a lampricide treatment at concentrations unlikely to harm aquatic organisms.

Journal of Great Lakes Research↗

Shining a light on Laurentian Great Lakes cisco (Coregonus artedi): How ice coverage may impact embryonic development

Changes in winter conditions, such as decreased ice coverage and duration, have been observed in the Laurentian Great Lakes for more than 20 years. Such changes have been hypothesized to be linked to low Coregonus spp. survival to age-1 as most cisco ( Coregonus artedi ) populations are autumn spawners whose embryos incubate under ice throughout the winter. The quantity of light during winter is regulated by ice coverage, and light affects embryo survival and development in some teleosts. We experimentally evaluated how cisco embryos from lakes Superior and Ontario respond to three light treatments that represented day-light intensity under 0–10, 40–60, and 90–100% ice coverage. Embryonic response measures included two developmental factors (embryo survival and incubation period) and two morphological traits (length-at-hatch and yolk-sac volume). Embryo survival was highest at the medium light treatment and decreased at high and low treatments for both populations, suggesting cisco may be adapted to withstand some light exposure from inter-annual variability in ice coverage. Light intensity had no overall effect on length of incubation. Increasing light intensity decreased length-at-hatch in Lake Superior but had no effect in Lake Ontario. Yolk-sac volume was positively correlated with increasing light in Lake Superior and negatively correlated in Lake Ontario. Contrasting responses in embryo development between lakes suggests differences in populations’ response to light is flexible. Our results provide a step towards better understanding the high variability observed in coregonine recruitment and may help predict what the future of this species may look like under current climate trends.

Journal of Great Lakes Research↗

Genomics reveals identity, phenology and population demographics of larval ciscoes (Coregonus artedi, C. hoyi, and C. kiyi) in the Apostle Islands, Lake Superior

We demonstrate, for the first time, the ability to reliably assign an assemblage of larval coregonines [Salmonidae Coregoninae] to shallow and multiple deepwater species. Larval coregonines from the Apostle Islands, Lake Superior, were genotyped using restriction site-associated DNA sequencing (RADseq) and were assigned to species using reference genotypes from adult corgonines from the same region. Of the 193 genotyped larvae, 101 were assigned as Coregonus artedi (average assignment probability = 97.6%), 57 were assigned as C. kiyi (average assignment probability = 95.5%), and 28 were assigned as C. hoyi (average assignment probability = 89.0%). Coregonus artedi were collected earliest in the season, followed by C. kiyi and then C. hoyi . Estimates of genetic diversity within each species provide a baseline for future monitoring in the Apostle Islands. Our success with species assignment indicates the promise of leveraging genomic data for larval coregonine identification, which could enable assessing and evaluating early life history dynamics and recruitment processes at the species level to the benefit of ongoing coregonine restoration and management efforts.

Wisconsin↗

Effects of warming winter embryo incubation temperatures on larval cisco (Coregonus artedi) survival, growth, and critical thermal maximum

Freshwater whitefishes, Salmonidae Coregoninae, are cold stenothermic fishes of ecological and socio-economic importance in northern hemisphere lakes that are warming in response to climate change. To address the effect of warming waters on coregonine reproduction we experimentally evaluated different embryo incubation temperatures on post-hatching survival, growth, and critical thermal maximum of larval cisco ( Coregonus artedi ) sampled from lakes Superior and Ontario. Embryos were incubated at water temperatures of 2.0, 4.4, 6.9, and 8.9 °C to simulate present and increased winter temperatures, and hatched larvae were reared in a common environment. For the populations from both lakes, larval survival and critical thermal maximum were negatively related to incubation temperature, and larval growth was positively related to incubation temperature. The magnitude of change across incubation temperatures was greater in the population sampled from Lake Superior than Lake Ontario for all traits examined. The more rapid decrease in survival and critical thermal maximum across incubation temperatures for larval cisco in Lake Superior, compared to those from Lake Ontario, suggests that Lake Superior larvae may possess a more limited ability to acclimate to and cope with increasing winter water temperatures. However, the rapid increase in growth rates across incubation temperatures in Lake Superior larvae suggests they could recover better from hatching at a small length induced by warm winters, as compared to Lake Ontario larvae. Our results suggest propagation and restoration programs may want to consider integrating natural habitat preferences and maximizing phenotypic variability to ensure offspring are set up for success upon stocking.

Lake Ontario, Lake Superior↗

Indicators of fish population responses to avian predation with focus on double-crested cormorants

Double-crested cormorants (Nannopterum auritum) have been implicated as causes of fish population declines in many locations across their breeding range. Two challenges facing managers are identifying fisheries population metrics indicative of cormorant impacts and determining when this evidence becomes actionable. Building upon existing studies, we conducted a meta-analysis of eight data-rich systems across the Laurentian Great Lakes region of the United States for common fish population responses to changes in cormorant abundance. Specifically, we examined trends in mean total female length at age-3 (TL3), female mean length and age at 50 % maturity, and mean age evenness as indicated by Shannon’s Equitability Index. Annual observations for these metrics were independently regressed linearly against cormorant density by system for walleye (Sander vitreus), yellow perch (Perca flavescens), smallmouth bass (Micropterus dolomieu), and northern pike (Esox lucius) populations. TL3 was the most sensitive with 9 of the 14 datasets statistically significant (r2 range 0.29 to 0.86). Maturity metrics were moderately sensitive to trends in cormorant predation with mean total length at 50 % maturity significant in 4 out of 11 datasets (r2 range 0.27–0.41) and mean age at 50 % maturity significant in 3 out of 11 datasets (r2 range 0.12 – 0.51). Least sensitive was age evenness with the Shannon Index significant in 3 out of 12 datasets (r2 typically < 0.25). Of metrics tested, TL3 was the most reliable indicator of changes in cormorant effects despite varying system changes and management responses among locations.

Michigan, Minnesota, New York↗

Behavioral and reproductive effects of the lampricides TFM and TFM:1% Niclosamide on native freshwater mussels

The lampricides TFM (3-trifluoromethyl-4′-nitrophenol) and Niclosamide (NIC, 2′, 5-dichloro-4′-nitrosalicylanilide) are used to control sea lamprey populations in the Great Lakes and associated tributaries. Niclosamide is often used as an additive to TFM to reduce the amount of TFM required to control sea lamprey. Concern is growing over the risk that lampricide treatments pose to native freshwater mussels residing in streams. Our objectives were to determine the acute toxicity of TFM and TFM:NIC to free glochidia (removed from the marsupial gills), compare the relative toxicity of TFM and TFM:NIC between free glochidia and brooded glochidia (within the marsupial gills), determine if glochidia age influences toxicity, and assess if exposure of gravid mussels to TFM and TFM:NIC alters behavior and reproduction. Three acute toxicity tests (2:TFM, 1:TFM : NIC) were conducted with glochidia and adults of the plain pocketbook mussel ( Lampsilis cardium ). In tests with glochidia, viability did not differ across TFM and TFM : NIC concentrations that encompassed typical stream treatments. Glochidia age influenced toxicity as glochidia obtained later in the brooding season were less viable than glochidia obtained earlier in the brooding season. Exposure of adults to elevated concentrations of lampricides often resulted in behavioral effects, but rarely affected reproductive endpoints. Because mussels are long-lived (30 to 100 y), even intermittent and short duration exposures may cumulatively affect mussels over their lifetime. The risks posed by lampricide treatments in the Great Lakes would be further informed by research on the sublethal effects of lampricides, particularly effects on non-target organisms such as mussels.

Journal of Great Lakes Research↗

Visioning and conceptual framework for coordinating Great Lakes connecting waters research and monitoring

The Laurentian Great Lakes are connected via naturally occurring straits and rivers: St. Marys River, Straits of Mackinac, St. Clair-Detroit River System, Niagara River, and the St. Lawrence River. Despite the historical ecological and economic importance of these waters, international agreements (e.g., Great Lakes Water Quality Agreement) only recently explicitly named the Great Lakes Connecting Waters (GLCWs), requiring governments to address the challenges of adequate restoration and protection from historical use and degradation. Standardized research and monitoring activities are needed; however, there is no established mechanism for coordination across the GLCWs. A three-day summit in 2023 convened experts to form the initial framework for a GLCWs Collaborative to increase standardizations and knowledge transfer. Participants drafted a governance structure and priorities following the principles of collective impact but allowed for place-based specificity for local connecting water organizations. Priorities and suggestions for success included: 1) co-development of the collaborative with all rights holders, stakeholders, and surrounding communities; 2) investment in research and technology specific to GLCWs; 3) investment in information transfer and training; 4) increased communication; and 5) better integration into existing Great Lakes research, monitoring, and funding programs. Expanding participation in all local GLCWs organizations with principles of inclusivity was identified as a larger collaborative goal. Next steps in the development of a GLCWs Collaborative include increased communication and formation of working groups and obtaining funding for a dedicated organization to begin supporting activities (communication, facilitation, logistics). By using a deliberate process for establishment, the potentially slower time frame for establishment may result in increased participation and success.

Laurentian Great Lakes↗

Tracking legacy mercury in the Hackensack River Estuary using mercury stable isotopes

Spatial redistribution of legacy mercury (Hg) contamination in the Hackensack River estuary (New Jersey, USA) was evaluated using mercury stable isotopes. Total Hg varied from 0.06 to 3.8 µg g-1 in sediment from the tidal Hackensack River and from 15 to 154 µg g-1 near historically contaminated sites in upper Berry's Creek, a tributary of the Hackensack River. delta202Hg values for total Hg from the Berry's Creek and Hackensack River estuaries varied over a fairly narrow range (-0.44‰ to -0.21‰), but were highest for sediment from upper Berry's Creek. Isotope mixing plots show that residual legacy mercury from upper Berry's Creek is partially diluted by a low concentration and low delta202Hg pool of mercury associated with low organic matter content sediments similar to those in Newark Bay. Based on an isotope mixing model, we estimate that upper Berry’s Creek contributes 21% to 82% of the mercury in sediments in the Hackensack River estuary and its tidal tributaries, including upstream marsh habitats far from the primary source. Our results show that mercury stable isotopes can be used to track the redistribution of mercury in tidal ecosystems and highlight the potentially large areas which may be affected by legacy mercury contamination in estuaries.

New Jersey, New York↗

Combining particle-tracking and geochemical data to assess public supply well vulnerability to arsenic and uranium

Flow-model particle-tracking results and geochemical data from seven study areas across the United States were analyzed using three statistical methods to test the hypothesis that these variables can successfully be used to assess public supply well vulnerability to arsenic and uranium. Principal components analysis indicated that arsenic and uranium concentrations were associated with particle-tracking variables that simulate time of travel and water fluxes through aquifer systems and also through specific redox and pH zones within aquifers. Time-of-travel variables are important because many geochemical reactions are kinetically limited, and geochemical zonation can account for different modes of mobilization and fate. Spearman correlation analysis established statistical significance for correlations of arsenic and uranium concentrations with variables derived using the particle-tracking routines. Correlations between uranium concentrations and particle-tracking variables were generally strongest for variables computed for distinct redox zones. Classification tree analysis on arsenic concentrations yielded a quantitative categorical model using time-of-travel variables and solid-phase-arsenic concentrations. The classification tree model accuracy on the learning data subset was 70%, and on the testing data subset, 79%, demonstrating one application in which particle-tracking variables can be used predictively in a quantitative screening-level assessment of public supply well vulnerability. Ground-water management actions that are based on avoidance of young ground water, reflecting the premise that young ground water is more vulnerable to anthropogenic contaminants than is old ground water, may inadvertently lead to increased vulnerability to natural contaminants due to the tendency for concentrations of many natural contaminants to increase with increasing ground-water residence time.

Journal of Hydrology↗

Comparison of bottom-track to global positioning system referenced discharges measured using an acoustic Doppler current profiler

A negative bias in discharge measurements made with an acoustic Doppler current profiler (ADCP) can be caused by the movement of sediment on or near the streambed. The integration of a global positioning system (GPS) to track the movement of the ADCP can be used to avoid the systematic negative bias associated with a moving streambed. More than 500 discharge transects from 63 discharge measurements with GPS data were collected at sites throughout the US, Canada, and New Zealand with no moving bed to compare GPS and bottom-track-referenced discharges. Although the data indicated some statistical bias depending on site conditions and type of GPS data used, these biases were typically about 0.5% or less. An assessment of differential correction sources was limited by a lack of data collected in a range of different correction sources and different GPS receivers at the same sites. Despite this limitation, the data indicate that the use of Wide Area Augmentation System (WAAS) corrected positional data is acceptable for discharge measurements using GGA as the boat-velocity reference. The discharge data based on GPS-referenced boat velocities from the VTG data string, which does not require differential correction, were comparable to the discharges based on GPS-referenced boat velocities from the differentially-corrected GGA data string. Spatial variability of measure discharges referenced to GGA, VTG and bottom-tracking is higher near the channel banks. The spatial variability of VTG-referenced discharges is correlated with the spatial distribution of maximum Horizontal Dilution of Precision (HDOP) values and the spatial variability of GGA-referenced discharges is correlated with proximity to channel banks.

Journal of Hydrology↗

The timing of scour and fill in a gravel-bedded river measured with buried accelerometers

A device that measures the timing of streambed scour and the duration of sediment mobilization at specific depths of a streambed was developed using data-logging accelerometers placed within the gravel substrate of the Cedar River, Washington, USA. Each accelerometer recorded its orientation every 20 min and remained stable until the surrounding gravel matrix mobilized as sediment was transported downstream and scour reached the level of the accelerometer. The accelerometer scour monitors were deployed at 26 locations in salmon-spawning habitat during the 2010&ndash;2011 flood season to record when the streambed was scoured to the depth of typical egg-pocket deposition. Scour was recorded at one location during a moderate high-flow event (65 m 3 /s; 1.25&ndash;1.5-year recurrence interval) and at 17 locations during a larger high-flow event (159 m 3 /s; 7-year recurrence interval). Accelerometer scour monitors recorded periods of intermittent sediment mobilization and stability within a high-flow event providing insight into the duration of scour. Most scour was recorded during the rising limb and at the peak of a flood hydrograph, though some scour occurred during sustained high flows following the peak of the flood hydrograph.

Journal of Hydrology↗

Characterizing response of total suspended solids and total phosphorus loading to weather and watershed characteristics for rainfall and snowmelt events in agricultural watersheds

Understanding the response of total suspended solids (TSS) and total phosphorus (TP) to influential weather and watershed variables is critical in the development of sediment and nutrient reduction plans. In this study, rainfall and snowmelt event loadings of TSS and TP were analyzed for eight agricultural watersheds in Wisconsin, with areas ranging from 14 to 110 km2 and having four to twelve years of data available. The data showed that a small number of rainfall and snowmelt runoff events accounted for the majority of total event loading. The largest 10% of the loading events for each watershed accounted for 73–97% of the total TSS load and 64–88% of the total TP load. More than half of the total annual TSS load was transported during a single event for each watershed at least one of the monitored years. Rainfall and snowmelt events were both influential contributors of TSS and TP loading. TSS loading contributions were greater from rainfall events at five watersheds, from snowmelt events at two watersheds, and nearly equal at one watershed. The TP loading contributions were greater from rainfall events at three watersheds, from snowmelt events at two watersheds and nearly equal at three watersheds. Stepwise multivariate regression models for TSS and TP event loadings were developed separately for rainfall and snowmelt runoff events for each individual watershed and for all watersheds combined by using a suite of precipitation, melt, temperature, seasonality, and watershed characteristics as predictors. All individual models and the combined model for rainfall events resulted in two common predictors as most influential for TSS and TP. These included rainfall depth and the antecedent baseflow. Using these two predictors alone resulted in an R 2 greater than 0.7 in all but three individual models and 0.61 or greater for all individual models. The combined model yielded an R 2 of 0.66 for TSS and 0.59 for TP. Neither the individual nor the combined models were substantially improved by using additional predictors. Snowmelt event models were statistically significant for individual and combined watershed models, but the model fits were not all as good as those for rainfall events (R 2 between 0.19 and 0.87). Predictor selection varied from watershed to watershed, and the common variables that were selected were not always selected in the same order. Influential variables were commonly direct measures of moisture in the watershed such as snowmelt, rainfall + snowmelt, and antecedent baseflow, or measures of potential snowmelt volume in the watershed such as air temperature.

Wisconsin↗

Nitrogen speciation and trends, and prediction of denitrification extent, in shallow US groundwater

Uncertainties surrounding nitrogen cycling complicate assessments of the environmental effects of nitrogen use and our understanding of the global carbon&ndash;nitrogen cycle. In this paper, we synthesize data from 877 ambient-monitoring wells across the US to frame broad patterns of nitrogen speciation and trends. At these sites, groundwater frequently contains substantial co-occurring NO 3 &minus; and XSN 2 (N 2 from denitrification), reflecting active/ongoing denitrification and/or a mixture of undenitrified and denitrified groundwater. NO 3 &minus; and NH 4 + essentially do not co-occur, indicating that the dominant source of NH 4 + at these sites likely is not dissimilatory reduction of NO 3 &minus; to NH 4 + . Positive correlations of NH 4 + with apparent age, CH 4 , dissolved organic carbon, and indicators of reduced conditions are consistent with NH 4 + mobilization from degradation of aquifer organic matter and contraindicate an anthropogenic source of NH 4 + for most sites. Glacial aquifers and eastern sand and gravel aquifers generally have lower proportions of NO 3 &minus; and greater proportions of XSN 2 than do fractured rock and karst aquifers and western sand and gravel aquifers. NO 3 &minus; dominates in the youngest groundwater, but XSN 2 increases as residence time increases. Temporal patterns of nitrogen speciation and concentration reflect (1) changing NO 3 &minus; loads over time, (2) groundwater residence-time controls on NH 4 + mobilization from solid phases, and (3) groundwater residence-time controls on denitrification. A simple classification tree using readily available variables (a national coverage of soil water depth, generalized geology) or variables reasonably estimated in many aquifers (residence time) identifies categorical denitrification extent (<10%, 10&ndash;50%, and >50%) with 79% accuracy in an independent testing set, demonstrating a predictive application based on the interconnected effects of redox, geology, and residence time.

Journal of Hydrology↗

Multi-scale preferential flow processes in an urban streambed under variable hydraulic conditions

Spatially preferential flow processes occur at nested scales at the sediment-water interface (SWI), due in part to sediment heterogeneities, which may be enhanced in flashy urban streams with heavy road sand influence. However, several factors, including the flow-rate dependence of preferential hyporheic flow and discrete groundwater discharge zones are commonly overlooked in reach-scale models of groundwater/surface water exchange. Using a series of controlled-head tracer-injection experiments coupled with cm-scale geophysics within the highly reactive upper 30 cm of the hyporheic zone of an urban stream, we quantified the flow dependence of local less-mobile porosity volume, mass-transfer rate coefficient, and the resulting local residence time in the less-mobile pore space at three controlled downward fluid fluxes (0.8, 2, and 3 m/d). Experiments were performed in two adjacent streambed locations, representing different sediment bulk vertical permeability. Less-mobile porosity parameters were generally substantial and similar between the two streambed locations; though a more competent, thin, organic layer at ∼15 cm depth in one location strongly impacted tracer loading, flushing dynamics, and local residence times. Increased downward flux led to (1) a decrease in less-mobile porosity residence time in all experiments, and (2) an increase in less-mobile porosity fraction for most experiments. Additionally, at the larger stream reach-scale, surface electrodes for electrical resistivity measurement were installed along 22 m of the wetted stream channel. These surface electrode measurements were collected during a natural storm flow event, which revealed widespread, short-term, flushing (e.g. <3 h) of the hyporheic zone with stream water, followed by longer-term (e.g. >60 h) flushing of the SWI with riparian zone groundwater. Flow dependence of preferential hyporheic zone flowpaths, like in the controlled tracer experiments, was also observed in these reach-scale electrical resistivity tomography measurements. Our findings reveal that the spatial and temporal dependence of preferential flow processes create highly dynamic SWI conditions that will affect the physical and coupled biogeochemical functions of the SWI in urbanized, sand-impacted streams.

Journal of Hydrology↗

Anatahan, Northern Mariana Islands: Reconnaissance geological observations during and after the volcanic crisis of spring 1990, and monitoring prior to the May 2003 eruption

Anatahan island is 9.5 km east–west by 3.5 km north–south and truncated by an elongate caldera 5 km east–west by 2.5 km north–south. A steep-walled pit crater ∼1 km across and ∼200 m deep occupies the eastern part of the caldera. The island is the summit region of a mostly submarine stratovolcano. The oldest subaerial rocks (stage 1) are exposed low on the outer flanks and in the caldera walls. These include thick (∼10 m) and thin (2–3 m) lava flows, well-indurated tuffs, and scoria units that make up the bulk of the island. Rock compositions range from basaltic andesite to dacite, and most are plagioclase-phyric. On the steep north and south flanks of the volcano, these rocks are cut by numerous east–west-oriented, few-hundred-m-long lineaments of undetermined origin. Indurated breccias unconformably overlie scarps cut into stage 1 units low on the south flank. Intermediate-age eruptive units (stage 2) include caldera-filling lava flows and pyroclastic deposits and, on the outer flanks, vents and valley-filling lava flows. The youngest pre-2003 volcanic unit on Anatahan (stage 3) is a hydromagmatic surge and fall deposit rich in accretionary lapilli. Prior to 2003, this unit was found over almost the entire island, and in many places original depositional surfaces and outcrops could be found in high-energy environments along the coast, indicating a young (but undetermined) age. During reconnaissance visits in 1990, 1992, 1994, and 2001, geothermal activity (fumaroles as well as pits with boiling, sediment-laden pools) was observed in the southern part of the pit crater. In March and April 1990, increased local seismicity, a large regional earthquake, and reported increased fumarolic activity in the pit crater prompted evacuation of Anatahan village, at the west end of the island. Our first field investigation took place in late April 1990 to assess the level of volcanic unrest, conduct reconnaissance geological observations, collect rock and geothermal water samples, and set up a geophysical monitoring network. Results at this time were inconclusive with respect to determining whether the activity was anomalous. Water in some of the geothermal pits within the pit crater was boiling, and pH values as low as 0.7 were recorded in the field. An electronic distance measurement (EDM) network was installed, and over a ∼1-week period, up to 9 cm of extension occurred across some lines but not others. Seismicity was characterized by intermittent local earthquakes but no sustained swarms or tremor. A brief visit in June 1990 revealed that the shallow lake near the boiling pits was gone, but activity in the pits themselves was similar to that of April 1990. Only minor extension had occurred along a single EDM line since the previous visit, and no earthquakes >M2.5 occurred during the visit. Subsequent 1- to 2-day visits occurred in October 1990, May 1992, May 1994, and June 2001. Activity within the geothermal pits was relatively constant during every visit, although during this 11-year period the level of the water in each pit decreased. In June 2001, a ∼50-m-wide region of mud pots and steaming ground in the central part of the geothermal area had developed. No geologic evidence, however, suggested that an eruption would occur < 2 years afterward. Most of the EDM lines showed slight extension between late 1990 and 1992, followed by very gradual contraction from 1992 to 2001. A more extensive seismic-monitoring system was installed on the Northern Mariana Islands during these visits, and it recorded a small seismic swarm at Anatahan from May to July 1993. The telemetry component of the seismic equipment broke prior to 2001 and had not been repaired by the time of the May 2003 eruption, so no precursory seismic data were recorded to indicate pre-eruption unrest.

Northern Mariana Islands↗

Characteristics, extent and origin of hydrothermal alteration at Mount Rainier Volcano, Cascades Arc, USA: Implications for debris-flow hazards and mineral deposits

Hydrothermal alteration at Mount Rainier waxed and waned over the 500,000-year episodic growth of the edifice. Hydrothermal minerals and their stable-isotope compositions in samples collected from outcrop and as clasts from Holocene debris-flow deposits identify three distinct hypogene argillic/advanced argillic hydrothermal environments: magmatic-hydrothermal, steam-heated, and magmatic steam (fumarolic), with minor superimposed supergene alteration. The 3.8 km 3 Osceola Mudflow (5600 y BP) and coeval phreatomagmatic F tephra contain the highest temperature and most deeply formed hydrothermal minerals. Relatively deeply formed magmatic-hydrothermal alteration minerals and associations in clasts include quartz (residual silica), quartz–alunite, quartz–topaz, quartz–pyrophyllite, quartz–dickite/kaolinite, and quartz–illite (all with pyrite). Clasts of smectite–pyrite and steam-heated opal–alunite–kaolinite are also common in the Osceola Mudflow. In contrast, the Paradise lahar, formed by collapse of the summit or near-summit of the edifice at about the same time, contains only smectite–pyrite and near-surface steam-heated and fumarolic alteration minerals. Younger debris-flow deposits on the west side of the volcano (Round Pass and distal Electron Mudflows) contain only low-temperature smectite–pyrite assemblages, whereas the proximal Electron Mudflow and a < 100 y BP rock avalanche on Tahoma Glacier also contain magmatic-hydrothermal alteration minerals that are exposed in the avalanche headwall of Sunset Amphitheater, reflecting progressive incision into deeper near-conduit alteration products that formed at higher temperatures. The pre-Osceola Mudflow alteration geometry is inferred to have consisted of a narrow feeder zone of intense magmatic-hydrothermal alteration limited to near the conduit of the volcano, which graded outward to more widely distributed, but weak, smectite–pyrite alteration within 1 km of the edifice axis, developed chiefly in porous breccias. The edifice was capped by a steam-heated alteration zone, most of which resulted from condensation of fumarolic vapor and oxidation of H 2 S in the unsaturated zone above the water table. Weakly developed smectite–pyrite alteration extended into the west and east flanks of the edifice, spatially associated with dikes that are localized in those sectors; other edifice flanks lack dikes and associated alteration. The Osceola collapse removed most of the altered core and upper east flank of the volcano, but intensely altered rocks remain on the uppermost west flank. Major conclusions of this study are that: (1) Hydrothermal–mineral assemblages and distributions at Mount Rainier can be understood in the framework of hydrothermal processes and environments developed from studies of ore deposits formed in analogous settings. (2) Frequent eruptions supplied sufficient hot magmatic fluid to alter the upper interior of the volcano hydrothermally, despite the consistently deep (≥ 8 km) magma reservoir which may have precluded formation of economic mineral deposits within or at shallow depths beneath Mount Rainier. The absence of indicator equilibrium alteration-mineral assemblages in the debris flows that effectively expose the volcano to a depth of 1–1.5 km also suggests a low potential for significant high-sulfidation epithermal or porphyry-type mineral deposits at depth. (3) Despite the long and complex history of the volcano, intensely altered collapse-prone rocks were spatially restricted to near the volcano's conduit system and summit, and short distances onto the upper east and west flanks, due to the necessary supply of reactive components carried by ascending magmatic fluids. (4) Intensely altered rocks were removed from the summit, east flank, and edifice interior by the Osceola collapse, but remain on the upper west flank in the Sunset Amphitheater area and present a continuing collapse hazard. (5) Visually conspicuous rocks on the lower east and mid-to-lower west flanks are not intensely altered and probably have not significantly weakened the rock, and thus do not present significant collapse hazards. (6) Alteration developed most intensely within breccia units, because of their high permeability and porosity. Volcanoes with abundant near-conduit upper-edifice breccias are prone to alteration increasing the possibility of collapse, whereas those that are breccia-poor (e.g., massive domes) are less prone to alteration.

Journal of Volcanology and Geothermal Research↗

High-resolution satellite and airborne thermal infrared imaging of precursory unrest and 2009 eruption of Redoubt Volcano, Alaska

A combination of satellite and airborne high-resolution visible and thermal infrared (TIR) image data detected and measured changes at Redoubt Volcano during the 2008&ndash;2009 unrest and eruption. The TIR sensors detected persistent elevated temperatures at summit ice-melt holes as seismicity and gas emissions increased in late 2008 to March 2009. A phreatic explosion on 15 March was followed by more than 19 magmatic explosive events from 23 March to 4 April that produced high-altitude ash clouds and large lahars. Two (or three) lava domes extruded and were destroyed between 23 March and 4 April. After 4 April, the eruption extruded a large lava dome that continued to grow until at least early July 2009.

Alaska↗