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,369 records · Page 76Linked to original sources
USGS2015

Context of ancient aqueous environments on Mars from in situ geologic mapping at Endeavour Crater

Using the Mars Exploration Rover Opportunity , we have compiled one of the first field geologic maps on Mars while traversing the Noachian terrain along the rim of the 22 km diameter Endeavour Crater (Latitude −2°16′33″, Longitude −5°10′51″). In situ mapping of the petrographic, elemental, structural, and stratigraphic characteristics of outcrops and rocks distinguishes four mappable bedrock lithologic units. Three of these rock units predate the surrounding Burns formation sulfate-rich sandstones and one, the Matijevic Formation, represents conditions on early Mars predating the formation of Endeavour Crater. The stratigraphy assembled from these observations includes several geologic unconformities. The differences in lithologic units across these unconformities record changes in the character and intensity of the Martian aqueous environment over geologic time. Water circulated through fractures in the oldest rocks over periods long enough that texturally and elementally significant alteration occurred in fracture walls. These oldest pre-Endeavour rocks and their network of mineralized and altered fractures were preserved by burial beneath impact ejecta and were subsequently exhumed and exposed. The alteration along joints in the oldest rocks and the mineralized veins and concentrations of trace metals in overlying lithologic units is direct evidence that copious volumes of mineralized and/or hydrothermal fluids circulated through the early Martian crust. The wide range in intensity of structural and chemical modification from outcrop to outcrop along the crater rim shows that the ejecta of large (>8 km in diameter) impact craters is complex. These results imply that geologic complexity is to be anticipated in other areas of Mars where cratering has been a fundamental process in the local and regional geology and mineralogy.

Journal of Geophysical Research E: Planets↗
USGS2016

Fault segmentation: New concepts from the Wasatch Fault Zone, Utah, USA

The question of whether structural segment boundaries along multisegment normal faults such as the Wasatch fault zone (WFZ) act as persistent barriers to rupture is critical to seismic hazard analyses. We synthesized late Holocene paleoseismic data from 20 trench sites along the central WFZ to evaluate earthquake rupture length and fault segmentation. For the youngest (<3 ka) and best-constrained earthquakes, differences in earthquake timing across prominent primary segment boundaries, especially for the most recent earthquakes on the north-central WFZ, are consistent with segment-controlled ruptures. However, broadly constrained earthquake times, dissimilar event times along the segments, the presence of smaller-scale (subsegment) boundaries, and areas of complex faulting permit partial-segment and multisegment (e.g., spillover) ruptures that are shorter (~20–40 km) or longer (~60–100 km) than the primary segment lengths (35–59 km). We report a segmented WFZ model that includes 24 earthquakes since ~7 ka and yields mean estimates of recurrence (1.1–1.3 kyr) and vertical slip rate (1.3–2.0 mm/yr) for the segments. However, additional rupture scenarios that include segment boundary spatial uncertainties, floating earthquakes, and multisegment ruptures are necessary to fully address epistemic uncertainties in rupture length. We compare the central WFZ to paleoseismic and historical surface ruptures in the Basin and Range Province and central Italian Apennines and conclude that displacement profiles have limited value for assessing the persistence of segment boundaries but can aid in interpreting prehistoric spillover ruptures. Our comparison also suggests that the probabilities of shorter and longer ruptures on the WFZ need to be investigated.

Utah↗
USGS2018

Direct measurements of mean Reynolds stress and ripple roughness in the presence of energetic forcing by surface waves

Direct covariance observations of the mean flow Reynolds stress and sonar images of the seafloor collected on a wave‐exposed inner continental shelf demonstrate that the drag exerted by the seabed on the overlying flow is consistent with boundary layer models for wave‐current interaction, provided that the orientation and anisotropy of the bed roughness are appropriately quantified. Large spatial and temporal variations in drag result from nonequilibrium ripple dynamics, ripple anisotropy, and the orientation of the ripples relative to the current. At a location in coarse sand characterized by large two‐dimensional orbital ripples, the observed drag shows a strong dependence on the relative orientation of the mean current to the ripple crests. At a contrasting location in fine sand, where more isotropic sub‐orbital ripples are observed, the sensitivity of the current to the orientation of the ripples is reduced. Further, at the coarse site under conditions when the currents are parallel to the ripple crests and the wave orbital diameter is smaller than the wavelength of the relic orbital ripples, the flow becomes hydraulically smooth. This transition is not observed at the fine site, where the observed wave orbital diameter is always greater than the wavelength of the observed sub‐orbital ripples. Paradoxically, the dominant along‐shelf flows often experience lower drag at the coarse site than at the fine site, despite the larger ripples, highlighting the complex dynamics controlling drag in wave‐exposed environments with heterogeneous roughness.

Journal of Geophysical Research C: Oceans↗
USGS2007

Diagenesis of diatomite from the Kolubara Coal Basin, Baroševac, Serbia

Diatomite associated with the Kolubara Coal Basin was studied to better understand early stage silica diagenesis of shallow water deposits. The Kolubara Basin consists of Neogene siliciclastic rocks, diatomite, marlstone and rare carbonates. Palaeozoic metamorphic and Mesozoic sedimentary and igneous basement rocks are transgressively overlain by Upper Miocene sandstone, siltstone, shale and mudstone. This Upper Miocene section is transgressively overlain by the Pontian section, which contains diatomite and coal beds. White and grey diatomite forms beds 0.7-2.2 m thick that are continuous over an area of about 2 km 2 . Siliceous rocks vary in composition from diatomite (81-89 per cent SiO 2 ) to diatom-bearing shale (58-60 per cent SiO 2 ). Siliceous deposits are laminated in places, with the laminae defined by variations in clay minerals, organic matter and diatoms. Diatomite shows only incipient diagenesis characterized by the fragmentation of diatom frustules, the minor to moderate corrosion of frustules and the formation of minor amounts of opal-A' (X-ray amorphous inorganic opal) cement. The low degree of diagenesis results from the young age of the deposits, low burial temperatures and possibly also from the presence of abundant organic matter and the dissolution of kaolinite. The presence of only weak diagenesis is also reflected by the characteristically poor consolidation of the rocks and low rank of the associated coal.

Baroševac↗
USGS2025

Changing dynamic phosphorus forms from field to stream during surface runoff events

The risk of water quality impairment from agricultural runoff depends on nutrient source, transport, and bioavailability. Phosphorus (P) spirals between dissolved and particulate forms as it is transported with suspended sediment (SS) from agricultural fields, through the stream network, to receiving water bodies. This dynamic sorption-desorption influences bioavailability. We quantified P form and abundance in samples collected during surface-runoff events from a farm field in the East River Basin, Wisconsin and compared them to those in stream water collected from the East River. We sampled five events between late March 2022 and June 2023. During most events, P in surface runoff was mainly in dissolved form, with particulate P sorbed to fine clay, the most abundant particle fraction transported from the field, whereas P in stream water was mainly in particulate form and sorbed to silt, even though fine clay was the most abundant particle fraction in the stream during events. Overall capacity for P sorption to SS in the stream varied among events. Total P and SS concentrations were lower during summer baseflow conditions and smaller surface runoff events; however, what SS was present was more P enriched. This shift in P form from field to stream indicates a potential for sorbing dissolved P to SS during transport through the stream network, which changes the bioavailability of P exported downstream with less bioavailable P as dissolved P binds to SS.

Wisconsin↗
USGS2025

Assessing the topographic distribution of legacy soil phosphorus in agricultural fields of the Delmarva Peninsula, Mid-Atlantic Coastal Plain, USA

Phosphorus (P) management remains a challenge in agricultural watersheds. The Choptank River Conservation Effects Assessment Project watershed, located in Maryland and Delaware and draining to the Chesapeake Bay, contains legacy soil P from historical dairy and poultry manure applications. These practices elevated soil P beyond crop needs, contributing to persistent P export to aquatic ecosystems. We assessed spatial P distribution and analyzed GIS (Geographic Information Systems)-derived landscape features driving legacy P movement on a farm (47 ha). We hypothesized that P accumulates in drained lowlands and depressional areas due to gravity-driven processes that accelerate P-enriched water to receiving waters via overland flow. In collaboration with the US Department of Agriculture Legacy P project, we collected 105 soil samples (0- to 5-cm and 5- to 15-cm depths) and 14 ditch sediment samples across five topographic openness classes from a farm with >100 years of dairy manure application. Average Mehlich-III P concentrations were 218 and 179 mg kg −1 at 0- to 5-cm and 5- to 15-cm depths, respectively, with legacy areas defined by P content > 100 mg kg −1 . Soil P and clay particle size were positively correlated ( r = 0.42, p < 0.05), increased as landscape openness decreased, and were negatively correlated with topographic openness (ranging from −0.2 to −0.4, p < 0.05), indicating accumulation of P and clay in low-lying areas. These patterns suggest that historical field-level managements have primarily shaped P distribution, while hydrologic and landscape properties further influence its redistribution via transport pathways and drainage. These findings support the development of landscape models to map critical source areas in low-relief watersheds and guide targeted mitigation in high-risk P export zones.

Maryland↗
USGS2013

Heterogeneous rupture in the great Cascadia earthquake of 1700 inferred from coastal subsidence estimates

Past earthquake rupture models used to explain paleoseismic estimates of coastal subsidence during the great A.D. 1700 Cascadia earthquake have assumed a uniform slip distribution along the megathrust. Here we infer heterogeneous slip for the Cascadia margin in A.D. 1700 that is analogous to slip distributions during instrumentally recorded great subduction earthquakes worldwide. The assumption of uniform distribution in previous rupture models was due partly to the large uncertainties of then available paleoseismic data used to constrain the models. In this work, we use more precise estimates of subsidence in 1700 from detailed tidal microfossil studies. We develop a 3-D elastic dislocation model that allows the slip to vary both along strike and in the dip direction. Despite uncertainties in the updip and downdip slip extensions, the more precise subsidence estimates are best explained by a model with along-strike slip heterogeneity, with multiple patches of high-moment release separated by areas of low-moment release. For example, in A.D. 1700, there was very little slip near Alsea Bay, Oregon (~44.4°N), an area that coincides with a segment boundary previously suggested on the basis of gravity anomalies. A probable subducting seamount in this area may be responsible for impeding rupture during great earthquakes. Our results highlight the need for more precise, high-quality estimates of subsidence or uplift during prehistoric earthquakes from the coasts of southern British Columbia, northern Washington (north of 47°N), southernmost Oregon, and northern California (south of 43°N), where slip distributions of prehistoric earthquakes are poorly constrained.

British Columbia;Washington;Oregon;California↗
USGS2013

The role of viscous magma mush spreading in volcanic flank motion at Kīlauea Volcano, Hawai‘i

Multiple mechanisms have been suggested to explain seaward motion of the south flank of Kīlauea Volcano, Hawai‘i. The consistency of flank motion during both waxing and waning magmatic activity at Kīlauea suggests that a continuously acting force, like gravity body force, plays a substantial role. Using finite element models, we test whether gravity is the principal driver of long-term motion of Kīlauea's flank. We compare our model results to geodetic data from Global Positioning System and interferometric synthetic aperture radar during a time period with few magmatic and tectonic events (2000-2003), when deformation of Kīlauea was dominated by summit subsidence and seaward motion of the south flank. We find that gravity-only models can reproduce the horizontal surface velocities if we incorporate a regional décollement fault and a deep, low-viscosity magma mush zone. To obtain quasi steady state horizontal surface velocities that explain the long-term seaward motion of the flank, we find that an additional weak zone is needed, which is an extensional rift zone above the magma mush. The spreading rate in our model is mainly controlled by the magma mush viscosity, while its density plays a less significant role. We find that a viscosity of 2.5 × 1017–2.5 × 1019 Pa s for the magma mush provides an acceptable fit to the observed horizontal surface deformation. Using high magma mush viscosities, such as 2.5 × 1019 Pa s, the deformation rates remain more steady state over longer time scales. These models explain a significant amount of the observed subsidence at Kīlauea's summit. Some of the remaining subsidence is probably a result of magma withdrawal from subsurface reservoirs.

Hawai'i↗
USGS2013

Effects of riparian vegetation on topographic change during a large flood event, Rio Puerco, New Mexico, USA

The spatial distribution of riparian vegetation can strongly influence the geomorphic evolution of dryland rivers during large floods. We present the results of an airborne lidar differencing study that quantifies the topographic change that occurred along a 12 km reach of the Lower Rio Puerco, New Mexico, during an extreme event in 2006. Extensive erosion of the channel banks took place immediately upstream of the study area, where tamarisk and sandbar willow had been removed. Within the densely vegetated study reach, we measure a net volumetric change of 578,050 ± ∼ 490,000 m 3 , with 88.3% of the total aggradation occurring along the floodplain and channel and 76.7% of the erosion focusing on the vertical valley walls. The sediment derived from the devegetated reach deposited within the first 3.6 km of the study area, with depth decaying exponentially with distance downstream. Elsewhere, floodplain sediments were primarily sourced from the erosion of valley walls. Superimposed on this pattern are the effects of vegetation and valley morphology on sediment transport. Sediment thickness is seen to be uniform among sandbar willows and highly variable within tamarisk groves. These reach-scale patterns of sedimentation observed in the lidar differencing likely reflect complex interactions of vegetation, flow, and sediment at the scale of patches to individual plants.

New Mexico↗
USGS2016

Effects of lek count protocols on greater sage-grouse population trend estimates

Annual counts of males displaying at lek sites are an important tool for monitoring greater sage-grouse populations ( Centrocercus urophasianus ), but seasonal and diurnal variation in lek attendance may increase variance and bias of trend analyses. Recommendations for protocols to reduce observation error have called for restricting lek counts to within 30&thinsp;minutes of sunrise, but this may limit the number of lek counts available for analysis, particularly from years before monitoring was widely standardized. Reducing the temporal window for conducting lek counts also may constrain the ability of agencies to monitor leks efficiently. We used lek count data collected across Wyoming during 1995&minus;2014 to investigate the effect of lek counts conducted between 30&thinsp;minutes before and 30, 60, or 90&thinsp;minutes after sunrise on population trend estimates. We also evaluated trends across scales relevant to management, including statewide, within Working Group Areas and Core Areas, and for individual leks. To further evaluate accuracy and precision of trend estimates from lek count protocols, we used simulations based on a lek attendance model and compared simulated and estimated values of annual rate of change in population size (&lambda;) from scenarios of varying numbers of leks, lek count timing, and count frequency (counts/lek/year). We found that restricting analyses to counts conducted within 30&thinsp;minutes of sunrise generally did not improve precision of population trend estimates, although differences among timings increased as the number of leks and count frequency decreased. Lek attendance declined >30&thinsp;minutes after sunrise, but simulations indicated that including lek counts conducted up to 90&thinsp;minutes after sunrise can increase the number of leks monitored compared to trend estimates based on counts conducted within 30&thinsp;minutes of sunrise. This increase in leks monitored resulted in greater precision of estimates without reducing accuracy. Increasing count frequency also improved precision. These results suggest that the current distribution of count timings available in lek count databases such as that of Wyoming (conducted up to 90&thinsp;minutes after sunrise) can be used to estimate sage-grouse population trends without reducing precision or accuracy relative to trends from counts conducted within 30&thinsp;minutes of sunrise. However, only 10% of all Wyoming counts in our sample (1995&minus;2014) were conducted 61&minus;90&thinsp;minutes after sunrise, and further increasing this percentage may still bias trend estimates because of declining lek attendance.

Journal of Wildlife Management↗
USGS2011

Seasonal variation in nutritional characteristics of the diet of greater white-fronted geese

We studied diet and habitat use of greater white-fronted geese ( Anser albifrons ) from autumn through spring on their primary staging and wintering areas in the Pacific Flyway, 1979–1982. There have been few previous studies of resource use and forage quality of wintering greater white-fronted geese in North America, and as a consequence there has been little empirical support for management practices pertaining to habitat conservation of this broadly distributed species. Observations of >2,500 flocks of geese and collections of foraging birds revealed seasonal and geographic variation in resource use reflective of changes in habitat availability, selection, and fluctuating physiological demands. Autumn migrants from Alaska arrived first in the Klamath Basin of California and southern Oregon, where they fed on barley, oats, wheat, and potatoes. Geese migrated from the Klamath Basin into the Central Valley of California in late autumn where they exploited agricultural crops rich in soluble carbohydrates, with geese in the Sacramento Valley feeding almost exclusively on rice and birds on the Sacramento–San Joaquin Delta primarily utilizing corn. White-fronted geese began their northward migration in late winter, and by early spring most had returned to the Klamath Basin where 37% of flocks were found in fields of new growth cultivated and wild grasses. Cereal grains and potatoes ingested by geese were low in protein (7–14%) and high in soluble nutrients (17–47% neutral detergent fiber [NDF]), whereas grasses were low in available energy (47–49% NDF) but high in protein (26–42%). Greater white-fronted geese are generalist herbivores and can exploit a variety of carbohydrate-rich cultivated crops, likely making these geese less susceptible to winter food shortages than prior to the agriculturalization of the North American landscape. However, agricultural landscapes can be extremely dynamic and may be less predictable in the long-term than the historic environments to which geese are adapted. Thus far greater white-fronted geese have proved resilient to changes in land cover in the Pacific Flyway and by altering their migration regime have even been able to adapt to changes in the availability of suitable forage crops.

Pacific Flyway↗
USGS2016

Proximate influences on female dispersal in white-tailed deer

Ultimate causes of animal dispersal have been hypothesized to benefit the dispersing individual because dispersal reduces competition for local resources, potential for inbreeding, and competition for breeding partners. However, proximate cues influence important features of dispersal behavior, including when dispersal occurs, how long it lasts, and direction, straightness, and distance of the dispersal path. Therefore, proximate cues that affect dispersal influence ecological processes (e.g., population dynamics, disease transmission, gene flow). We captured and radio-marked 277 juvenile female white-tailed deer ( Odocoileus virginianus ), of which 27 dispersed, to evaluate dispersal behavior and to determine proximate cues that may influence dispersal behavior. Female dispersal largely occurred at 1 year of age and coincided with the fawning season. Dispersal paths varied but generally were non-linear and prolonged. Physical landscape features (i.e., roadways, rivers, residential areas) influenced dispersal path direction and where dispersal terminated. Additionally, forays outside of the natal range that did not result in dispersal occurred among 52% of global positioning system (GPS)-collared deer ( n = 25) during the dispersal period. Our results suggest intra-specific social interactions and physical landscape features influence dispersal behavior in female deer. Female dispersal behavior, particularly the lack of directionality, the semi-permeable nature of physical barriers, and the frequency of forays outside of the natal range, should be considered in regard to population management and controlling the spread of disease.

Pennsylvania↗
USGS2017

Shorebird stopover habitat decisions in a changing landscape

To examine how habitat use by sandpipers ( Calidris spp.; Baird's sandpipers, dunlin, least sandpipers, pectoral sandpipers, semipalmated sandpipers, stilt sandpipers, and white-rumped sandpipers) varies across a broad suite of environmental conditions, we conducted surveys at wetlands throughout the spring migratory period in 2013 and 2014 in 2 important stopover regions: the Rainwater Basin (RWB) in Nebraska, USA, and the Prairie Pothole Region (PPR) in South Dakota, USA. Because providing adequate energetic resources for migratory birds is a high priority for wetland management, we also measured invertebrate abundance at managed wetlands in the RWB to determine how food abundance influences the occupancy and abundance of sandpipers on wetlands throughout the migratory period. To quantify habitat use, we surveyed wetlands every 7–10 days in both regions and visually estimated wetland attributes. Our results indicate that invertebrate abundance predicted occupancy, but not abundance, of sandpipers at wetlands in the RWB. The wetland vegetation characteristics that predict sandpiper occupancy are similar in both regions, but wetlands in the PPR support a higher abundance of sandpipers than wetlands in the RWB. Our results suggest that sandpipers make stopover decisions that balance local and regional wetland conditions. Managers should maintain the cues (i.e., mudflat) and ecological conditions beyond invertebrate abundance that predict sandpiper habitat use to successfully provide resources for sandpipers during migratory stopover if that is a goal of wetland management. © 2017 The Wildlife Society.

Nebraska, South Dakota↗
USGS2018

Estimates of tidal-marsh bird densities using Bayesian networks

Conserving tidal-marsh bird communities requires strategies to address continuing pressures from human development to the effects of increasing rates of sea-level rise. Knowing tidal-marsh bird distributions and population sizes are important for developing these strategies. In the Northeast United States, where estimates of sea-level rise are 3 times higher than the global average, 5 bird species are tidal-marsh specialists: clapper rail (Rallus crepitans), willet (Tringa semipalmata), Nelson's sparrow (Ammospiza nelsoni), saltmarsh sparrow (A. caudacuta), and seaside sparrow (A. maritima). We used a regional marsh bird survey to develop Bayesian network models to identify factors that influence patch-scale species density and to estimate regional population sizes. We modeled species density as a function of habitat covariates at the patch, local, landscape, and regional spatial scales. Densities were most sensitive to patch location and dimension, patch geomorphic setting, indices of human development, and changes in mean sea level. We estimated 110,000 clapper rails (95% CI = 61,000–159,000), 111,000 willets (95% CI = 70,000–152,000), 7,000 Nelson's sparrows (95% CI = 4,000–10,000), 60,000 saltmarsh sparrows (95% CI = 40,000–80,000), and 234,000 seaside sparrows (95% CI = 112,000–356,000) from the United States–Canada border to, and including, the mouth of the Chesapeake Bay, Virginia, USA. Our abundance estimates can be used to identify priority conservation areas at multiple geographic scales and our models help identify key habitat and landscape components for tidal-marsh restoration and management to benefit tidal-marsh birds and can be modified for other species.

Maine, New Hampshire, Massachusetts, Rhode Island,↗
USGS2025

Weather, habitat area, connectivity, and number of patches influence breeding ecology of ring-necked pheasants

Understanding habitat selection is critical in habitat prioritization for species of conservation and management concern. Information on habitat selection is particularly important for grassland bird species whose populations have suffered steep declines over the last few decades. We assessed ring-necked pheasants' ( Phasianus colchicus ) habitat selection in a dynamic agricultural landscape. The population dynamics of pheasants are partially related to nest survival, which may be influenced by the quality of nesting habitat. Consequently, knowledge of vegetation composition and structural characteristics associated with the selection and survival of nests would help inform management decisions to improve nest success. We monitored nests from 103 radio-collared pheasants inhabiting an agricultural landscape in South Dakota, USA, from 2017–2019 to determine the effect of landscape composition and configuration on nest-site selection and nest survival. We explored nesting behavior at 2 orders of selection: resource selection within the home range (third order) and selection of specific resource items from a resource patch (fourth order). Proportion of row crop and connectivity of row crop was negatively associated with nest-site selection at the third order. At the fourth order, pheasants tended to select for taller vegetation and greater percent grass cover than at paired random sites. Pheasants also selected areas with more grasslands. A 1% increase in grass cover and proportion of grassland increased the odds of nest-site selection by 1% and 2%, respectively. Connectivity of row crop patches was negatively associated with daily nest survival. We also evaluated factors affecting pheasant brood-site selection. A 1-unit increase in grass cover and Hemiptera biomass increased the odds of brood-site selection by 4%. The probability of brood-site selection also increased with fewer row crop patches. Weather played a crucial role in driving nest survival. The consideration of local weather trends and regional variation in habitat can inform habitat management for pheasants. Pheasant populations may benefit from research that identifies thermal landscapes and land management techniques that promote cooler microclimates for nesting and brood-rearing activities.

South Dakota↗
USGS2025

Evaluation of the acute toxicity of the piscicide TFM to Burbot

Non-target animal sensitivity remains a concern when treating Laurentian Great Lakes streams with 4-nitro-3-(trifluoromethyl)phenol (TFM), the main pesticide used to control Sea Lamprey Petromyzon marinus as part of the bi-national Great Lakes Fishery Commission's Sea Lamprey Control Program. Populations of Burbot Lota lota , a historically and culturally important fish, inhabit some of the streams that are treated with TFM. While many species of fish inhabiting the Great Lakes streams have been assessed for sensitivity to TFM, we are not aware of previous research to assess the risk to Burbot. We assessed the sensitivity of Burbot to TFM using replicate 12-hour flow-through diluter toxicity tests. We found Burbot to have a median lethal concentration (LC 50 ) of 9.74 mg/L, while the minimum lethal concentration (LC 99.9 ) for Sea Lamprey was predicted to be 2.5 mg/L in similar waters. The resulting toxicity ratio (LC 50 of non-target organism/LC 99.9 of Sea Lamprey) of Burbot was 3.90, well above the toxicity ratios for known sensitive species. Our results suggest Burbot are not expected to be adversely affected during a typical TFM stream treatment.

Great Lakes↗