SEARCH · USGS Science

Results for “The Explorers Journal”

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

SourceResearch connectionFrom yearThrough year

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

At least 523 records · Page 29Linked to original sources

USGS2021

Considerations of variability and power for long-term monitoring of stream fish assemblages

Little attention has been given to optimizing statistical power for monitoring stream fish assemblages. We explored the relationship between temporal variability and statistical power using 34 metrics from fish community data collected annually at six sites over 10 years via electrofishing. Metric variability differed by the life stage and group of species considered, use of abundance or mass data, and data standardization technique. Lower variability was associated with community data, abundance data, and time-based standardizations, while greater variability was associated with young-of-the-year data, mass data, and area-based standardizations. Simulation-based power analysis indicated metric choice, and to a lesser degree, monitoring design (annual, biennial, endpoints, or haphazard sampling) influenced power to detect change. Across a fixed number of surveys (N = 60), endpoints sampling performed best. The N needed to detect change was heavily dependent upon metric choice for all monitoring designs, with the most biologically specific metrics requiring greater N. Large savings in effort and resource expenditure can be obtained utilizing biologically relevant metrics that are robust to temporal noise within an appropriate sampling design.

Canadian Journal of Fisheries and Aquatic Sciences↗

USGS2024

Long-term communication of aftershock forecasts: The Canterbury earthquake sequence in New Zealand

On 14 February 2016, a magnitude (M)5.7 earthquake struck in Christchurch New Zealand (Aotearoa in the Maori language). The shaking caused damage to historic facades, power outages, cliff collapses, rock falls, and liquefaction but no reported injuries or fatalities. This Valentine’s Day earthquake was an aftershock in the Canterbury earthquake sequence (CES), which began on 4 September 2010 with the M7.1 Darfield Earthquake and included the destructive and fatal M6.2 Christchurch aftershock on 22 February 2011. This study, eight months after the Valentine’s Day earthquake and six years after the initiation of the CES, is the first to explore long-term aftershock forecast information and communication needs. The exploratory study also aimed to gather feedback on aftershock scenarios, an alternative form for communicating the forecast. The qualitative study involved workshops with emergency managers, public health officials, and members of the public in Christchurch. Key findings for long-term communication throughout an earthquake sequence include: 1. divergent earthquake experiences affect aftershock communication response and information needs; 2. understanding aftershock sequence behavior is foundational to sense-making when large aftershocks occur; 3. strategic earthquake sequence updates from the trusted science agency and local agencies could serve as important reminders for earthquake preparedness; 4. communication of aftershock forecast uncertainty could aid with both the credibility of the information and living with uncertainty, and 5. inclusion of impact information and preparedness advice into aftershock forecast scenarios could provide links to actionable information. The paper derives implications for research and practice of long-term communications during an aftershock sequence.

Christchurch↗

USGS2023

Effects of landscape cover and yard features on feral and free-roaming cat (Felis catus) distribution, abundance and activity patterns in a suburban area

Feral and free-roaming domestic cats ( Felis catus ) are invasive predators throughout the world. In some areas, cats occur in higher densities than native mammalian predators and can have severe effects upon prey populations. We set 48 wildlife game cameras in residential yards in Arkansas, USA, to evaluate which landscape and yard features influenced cat abundance occurring in yards. In addition, we quantified the daily activity patterns of free-roaming cats and explored how habitat features or predator activity influenced the timing of cat activity. We found that cats were present in 70.8% of yards with an average of three recognizable individuals per yard. Abundance of cats was higher than all native mesopredators except for raccoon ( Procyon lotor ) and Virginia opossum ( Didelphis virginiana ). Cat abundance and minimum population decreased when forest cover was high within 400 m of the camera. Cats were active at all times of the day but tended to be more diurnal in areas closer to city centers or in agricultural settings. Conversely, cats were more nocturnal later in the summer and in areas that had high levels of predator activity. Our results indicate that cats are widespread in this region and their relative abundance is driven more by landscape features than by yard features, possibly due to their large home ranges. Cats may alter their activity to better coexist with predators. Alteration in yard features is unlikely to be an effective deterrent for cats and more direct control measures may be necessary.

Arkansas↗

USGS2021

Contributed reports of widely felt earthquakes in California, United States: If they felt it, did they report it?

In a recent study, Hough and Martin (2021) considered the extent to which socioeconomic factors influence the numbers and distribution of contributed reports available to characterize the effects of both historical and recent large earthquakes. In this study I explore the question further, focusing on analysis of widely felt earthquakes near major population centers in northern and southern California since 2002. For most of these earthquakes there is a correlation between average household income in a postal ZIP code and the population-normalized rate of responses to the DYFI system. As past studies have demonstrated, there is also a strong correlation between DYFI participation and the severity of shaking. This first-order correlation can obscure correlations with other factors that influence participation. Focusing on five earthquakes between 2011 and 2021 that generated especially uniform shaking across the greater Los Angeles, California, region, response rate varies by two orders of magnitude across the region, with a clear correlation with demographics, and consistent spatial patterns in response rate for earthquakes 10 years apart. While there is no evidence that uneven DYFI participation in California impacts significantly the reliability of intensity data collected, the results reveal that DYFI participation is significantly higher in affluent parts of southern California compared to economically disadvantaged areas.

California↗

USGS2015

Host–parasite behavioral interactions in a recently introduced, whooping crane population

The whooping crane Grus americana has a long conservation history, but despite multiple attempts across North America, introduction success is lacking. Recently introduced, captively reared whooping cranes have had periods of poor reproductive performance in central Wisconsin that sometimes coincided with black fly (Diptera: Simuliidae) emergences. Sandhill crane Grus canadensis reproductive performance in central Wisconsin is approximately double that of whooping cranes. We used comfort behaviors as a measure of black fly harassment to infer whether behavioral differences existed between nesting sandhill cranes and nesting whooping cranes and between successful and unsuccessful whooping crane pairs. To further explore the interaction between black flies and incubating whooping cranes, we examined differences in behaviors between incubating birds and their off-nest mates. Compared to their off-nest mates, incubating whooping cranes exhibited elevated comfort behaviors, suggesting a bird at a nest may experience greater harassment from black flies. Sandhill cranes had elevated head-flicks over whooping cranes. Whooping cranes exhibited more head-rubs than sandhill cranes, and successful whooping crane pairs had elevated head-rubs over pairs that deserted their nests. Behavioral differences between sandhill cranes and whooping cranes as well as differences in reproductive performance, could be explained by exposure to local breeding conditions. Whereas sandhill cranes have nested in the area for hundreds, if not thousands, of years, whooping cranes were only recently introduced to the area. Behavioral differences between the species as well as those between successful and unsuccessful whooping crane pairs could also be explained by the effect of captive exposure, which could affect all whooping crane introductions.

Wisconsin↗

USGS2015

Stock-specific advection of larval walleye ( Sander vitreus ) in western Lake Erie: Implications for larval growth, mixing, and stock discrimination

Physical processes can generate spatiotemporal heterogeneity in habitat quality for fish and also influence the overlap of pre-recruit individuals (e.g., larvae) with high-quality habitat through hydrodynamic advection. In turn, individuals from different stocks that are produced in different spawning locations or at different times may experience dissimilar habitat conditions, which can underlie within- and among-stock variability in larval growth and survival. While such physically-mediated variation has been shown to be important in driving intra- and inter-annual patterns in recruitment in marine ecosystems, its role in governing larval advection, growth, survival, and recruitment has received less attention in large lake ecosystems such as the Laurentian Great Lakes. Herein, we used a hydrodynamic model linked to a larval walleye ( Sander vitreus ) individual-based model to explore how the timing and location of larval walleye emergence from several spawning sites in western Lake Erie (Maumee, Sandusky, and Detroit rivers; Ohio reef complex) can influence advection pathways and mixing among these local spawning populations (stocks), and how spatiotemporal variation in thermal habitat can influence stock-specific larval growth. While basin-wide advection patterns were fairly similar during 2011 and 2012, smaller scale advection patterns and the degree of stock mixing varied both within and between years. Additionally, differences in larval growth were evident among stocks and among cohorts within stocks which were attributed to spatiotemporal differences in water temperature. Using these findings, we discuss the value of linked physical–biological models for understanding the recruitment process and addressing fisheries management problems in the world's Great Lakes.

Journal of Great Lakes Research↗

USGS2006

Evidence of phyllosilicates in Wooly Patch, an altered rock encountered at West Spur, Columbia Hills, by the Spirit rover in Gusev crater, Mars

On its traverse to Columbia Hills, the Mars Exploration Rover Spirit investigated an outcrop designated “Wooly Patch” that exhibited morphological, mineralogical, and geochemical characteristics at the extreme ends of ranges observed among rocks studied at West Spur, a westward projecting salient near the foot of the Columbia Hills, Gusev crater. The major‐element composition and Fe‐mineralogy, as determined by the Alpha‐Particle X‐ray Spectrometer and Mössbauer Spectrometer, are inconsistent with any reasonable assemblage of basaltic minerals in that there is an excess of Si and Al. The combined data are best explained by the presence of 14–17% phyllosilicate minerals. Phyllosilicates that account for the composition and cation ratios include members of the kaolinite, serpentine, chlorite, and septechlorite groups. The potential existence of kaolinite‐type Al‐rich phyllosilicates within the Wooly Patch outcrop suggests a mildly acidic environment (pH 4–6) in the past and an open hydrologic system with good drainage conditions in the environment where these rocks were altered.

Journal of Geophysical Research E: Planets↗

USGS2026

Nonstationary demographic state-space models using unreplicated counts for species undergoing environmental stressors

A fundamental task in ecological statistics is to estimate abundance and growth rate distributions from wildlife monitoring data to inform conservation management. Modeling time series of wildlife populations presents a number of challenges from both statistical and ecological perspectives, including discreteness; lack of replication; nonstationarity; and observation, demographic, and other phenomenological processes. Nonstationary dynamics are often exhibited by populations undergoing environmental stressors. Models must account for these characteristics to produce reliable estimates of abundance and trends, yet estimation can be challenging with unreplicated data. We propose nonstationary demographic state-space models using unreplicated counts for populations undergoing environmental stressors. A reduced growth rate model matches the complexity of the unreplicated count data, and a fecundity bound on growth rate distributions allows the separation of processes affecting growth rates like environmental stressors from those affecting abundance external to growth rates like migration. NDSSMs allow for the embedding of nonstationary model components, and we explore the use of changepoints, volatility clustering, and migration processes. We apply the proposed nonstationary models in case studies of herons affected by predator/competitor reestablishment and three bat species affected by a fungal pathogen causing white-nose syndrome. Nonstationary models outperform stationary models and generalized linear mixed effects models according to model scoring and visual inspection of predictions, and provide estimates more consistent with published values. Incorporating migration improves model fit universally, even with approximate one-way immigration, most likely because populations are extirpated, recolonized, and increase multiple-fold over the upper bound set by species fecundity. In addition, estimates of the timing and severity of the environmental stressor differed for models with migration. Including nonstationary and demographic components in a fecundity-bounded growth rate model improves inference and benefits interpretability of hyperparameters. In turn, this adjusts uncertainties in predictions of abundance and growth rates over time, providing the ingredients needed for informed conservation analysis and for directing future monitoring of at-risk species.

Journal of Agricultural, Biological and Environmen↗

USGS2012

Comparing seasonal dynamics of the Lake Huron zooplankton community between 1983-1984 and 2007 and revisiting the impact of Bythotrephes planktivory

Zooplankton community composition can be influenced by lake productivity as well as planktivory by fish or invertebrates. Previous analyses based on long-term Lake Huron zooplankton data from August reported a shift in community composition between the 1980s and 2000s: proportional biomass of calanoid copepods increased while that of cyclopoid copepods and herbivorous cladocerans decreased. Herein, we used seasonally collected data from Lake Huron in 1983–1984 and 2007 and reported similar shifts in proportional biomass. We also used a series of generalized additive models to explore differences in seasonal abundance by species and found that all three cyclopoid copepod species ( Diacyclops thomasi, Mesocylops edax, Tropocyclops prasinus mexicanus ) exhibited higher abundance in 1983–1984 than in 2007. Surprisingly, only one ( Epischura lacustris ) of seven calanoid species exhibited higher abundance in 2007. The results for cladocerans were also mixed with Bosmina spp. exhibiting higher abundance in 1983–1984, while Daphnia galeata mendotae reached a higher level of abundance in 2007. We used a subset of the 2007 data to estimate not only the vertical distribution of Bythotrephes longimanus and their prey, but also the consumption by Bythotrephes in the top 20 m of water. This epilimnetic layer was dominated by copepod copepodites and nauplii, and consumption either exceeded (Hammond Bay site) or equaled 65% (Detour site) of epilimnetic zooplankton production. The lack of spatial overlap between Bythotrephes and herbivorous cladocerans and cyclopoid copepod prey casts doubt on the hypothesis that Bythotrephes planktivory was the primary driver underlying the community composition changes in the 2000s.

Lake Huron↗

USGS2018

Using geologic structures to constrain constitutive laws not accessible in the laboratory

In this essay, we explore a central problem of structural geology today, and in the foreseeable future, which is the determination of constitutive laws governing rock deformation to produce geologic structures. Although laboratory experiments provide much needed data and insights about constitutive laws, these experiments cannot cover the range of conditions and compositions relevant to the formation of geologic structures. We advocate that structural geologists address this limitation by interpreting natural experiments, documented with field and microstructural data, using continuum mechanical models that enable the deduction of constitutive laws. To put this procedure into a historical context, we review the founding of structural geology by James Hutton in the late 18th century, and the seminal contributions to continuum mechanics from Newton to Cauchy that provide the tools to model geologic structures. The procedure is illustrated with two examples drawn from recent and on-going field investigations of crustal and mantle lithologies . We conclude by pointing to future research opportunities that will engage structural geologists in the pursuit of constitutive laws during the 21st century.

Journal of Structural Geology↗

USGS2026

Availability of dark daytime refuge may limit mysid abundance in the Laurentian Great Lakes

The zooplankton Mysis diluviana is a major component of the Laurentian Great Lakes food web and has recently declined in abundance in both lakes Michigan and Huron. Drivers of these declines are not well understood. Here, we explore the hypothesis that recent increases in water clarity have contributed to the decline of M. diluviana (mysids) by limiting the availability of daytime dark refuge from visual predators. Using Secchi depth data from 1996 to 2021, we estimate that dark refuge has decreased substantially in lakes Michigan and Huron where mysids have declined, but dark refuge has remained more stable in lakes Ontario and Superior where mysid populations did not decline. Results from a 2021 uncrewed surface vessel hydroacoustic survey and lake-wide, net-based sampling in lakes Michigan and Huron revealed significantly more mysids in areas with dark refuge, such as Lake Michigan’s northern basin. Conversely, Lake Huron contains sparse dark refuge consistent with low mysid densities in that lake. Higher water clarity leading to increased predation may be a primary driver of mysid declines in the Great Lakes.

Laurentian Great Lakes↗

USGS2018

Museum metabarcoding: a novel method revealing gut helminth communities of small mammals across space and time

Natural history collections spanning multiple decades provide fundamental historical baselines to measure and understand changing biodiversity. New technologies such as next generation DNA sequencing (NGS) have considerably increased the potential of museum specimens to address significant questions regarding the impact of environmental changes on host and parasite/pathogen dynamics. We developed a new technique to identify intestinal helminth parasites and applied it to shrews (Eulipotyphla: Soricidae) because they are ubiquitous, occupy diverse habitats, and host a diverse and abundant parasite fauna. Notably, we included museum specimens preserved in various ways to explore the efficacy of using metabarcoding analyses that may enable identification of helminth symbiont communities from historical archives. We successfully sequenced the parasite communities (using 12S mtDNA, 16S mtDNA, 28S rDNA) of 23 whole gastrointestinal (GI) tracts. All GI tracts were obtained from the Museum of Southwestern Biology (MSB), USA, and from recent field collections, varying both in time since fixation (ranging from 4 months to 16 years) and preservation method (70% or 95% ethanol stored at room temperature, or flash frozen in liquid nitrogen and stored at -80°C). Our proof of concept demonstrates the feasibility of applying NGS techniques to authoritatively identify the parasite/pathogen communities within whole GI tracts from museum specimens of varying age and fixation, and the value of future preservation of host-associated whole GI tracts in public research archives. This powerful approach facilitates future comparative examinations of the distributions and interactions among multiple associated groups of organisms through time and space.

International Journal for Parasitology↗

USGS1993

Pressure increases, the formation of chromite seams, and the development of the ultramafic series in the Stillwater Complex, Montana

This paper explores the hypothesis that chromite seams in the Stillwater Complex formed in response to periodic increases in total pressure in the chamber. Total pressure increased because of the positive δ V of nucleation of CO 2 bubbles in the melt and their subsequent rise through the magma chamber, during which the bubbles increased in volume by a factor of 4–6. By analogy with the pressure changes in the summit chambers of Kilauea and Krafla volcanoes, the maximum variation was 0⋅2–0⋅25 kbar, or 5–10% of the total pressure in the Stillwater chamber. An evaluation of the likelihood of fountaining and mixing of a new, primitive liquid that entered the chamber with the somewhat more evolved liquid already in the chamber is based upon calculations using observed and inferred velocities and flow rates of basaltic magmas moving through volcanic fissures. The calculations indicate that hot, dense magma would have oozed, rather than fountained into the chamber, and early mixing of the new and residual magmas that could have resulted in chromite crystallizing alone did not take place.

Journal of Petrology↗

USGS2015

The role of the geophysical template and environmental regimes in controlling stream-living trout populations

The importance of multiple processes and instream factors to aquatic biota has been explored extensively, but questions remain about how local spatiotemporal variability of aquatic biota is tied to environmental regimes and the geophysical template of streams. We used an individual-based trout model to explore the relative role of the geophysical template versus environmental regimes on biomass of trout ( Oncorhynchus clarkii clarkii ). We parameterized the model with observed data from each of the four headwater streams (their local geophysical template and environmental regime) and then ran 12 simulations where we replaced environmental regimes (stream temperature, flow, turbidity) of a given stream with values from each neighboring stream while keeping the geophysical template fixed. We also performed single-parameter sensitivity analyses on the model results from each of the four streams. Although our modeled findings show that trout biomass is most responsive to changes in the geophysical template of streams, they also reveal that biomass is restricted by available habitat during seasonal low flow, which is a product of both the stream’s geophysical template and flow regime. Our modeled results suggest that differences in the geophysical template among streams render trout more or less sensitive to environmental change, emphasizing the importance of local fish–habitat relationships in streams.

Canadian Journal of Fisheries and Aquatic Sciences↗

USGS2022

Computed tomography for measuring body fat reserves in threatened Mohave desert tortoise (Gopherus agassizii)

Noninvasive methods for measuring fat reserves in both captive and free-ranging animals are important for monitoring individual and population health, but chelonian anatomy and physiology present challenges to accurate measurements. Standard field-based methods for assessing body condition in Mojave desert tortoises ( Gopherus agassizii ) involve the qualitative body condition score, which relies on the apparent height of the temporalis muscle relative to the sagittal crest (in addition to other characteristics) and quantitative body condition indices that measure relative mass at size. However, it is unclear how these metrics relate to body fat reserves in this species. The aims of this study were to (1) describe the use of noninvasive computed tomography in measuring body fat volume of Mojave desert tortoises, (2) describe the location of fat reserves, (3) investigate relationships between fat reserves and body condition score and body condition index, and (4) explore whether relative temporalis muscle depth, measured via computed tomography, correlates with body condition score. Body condition scores were assessed for eight captive Mojave desert tortoises prior to euthanasia, and computed tomography was performed postmortem to quantify fat volume and measure temporalis muscle depth. At necropsy, the distribution of fat was documented. Fat volume calculated by computed tomography ranged from 2.83 to 145.38 cm 3 (0.07–2.5% body volume). Neither qualitative body condition score nor quantitative body condition index was correlated with fat volume. Bladder content did not compromise body condition index. Body condition score was not correlated with relative temporalis muscle depth. Computed tomography is a noninvasive method for successfully identifying fat reserves and estimating total fat volume in Mojave desert tortoises. The lack of a relationship between computed tomography-determined metrics and commonly used body condition metrics indicates that computed tomography fills a critical gap in the health assessment tool kit for captive and free-ranging Mojave desert tortoises.

Journal of Zoo and Wildlife Medicine↗

USGS2017

Causes of distal volcano-tectonic seismicity inferred from hydrothermal modeling

Distal volcano-tectonic (dVT) seismicity typically precedes eruption at long-dormant volcanoes by days to years. Precursory dVT seismicity may reflect magma-induced fluid-pressure pulses that intersect critically stressed faults. We explored this hypothesis using an open-source magmatic-hydrothermal code that simulates multiphase fluid and heat transport over the temperature range 0 to 1200 °C. We calculated fluid-pressure changes caused by a small (0.04 km 3 ) intrusion and explored the effects of flow geometry (channelized vs. radial flow), magma devolatilization rates (0–15 kg/s), and intrusion depths (5 and 7.5 km, above and below the brittle-ductile transition). Magma and host-rock permeabilities were key controlling parameters and we tested a wide range of permeability ( k ) and permeability anisotropies ( k h / k v ), including k constant, k ( z ), k ( T ), and k ( z , T , P ) distributions, examining a total of ~ 1600 realizations to explore the relevant parameter space. Propagation of potentially causal pressure changes (Δ P ≥ 0.1 bars) to the mean dVT location (6 km lateral distance, 6 km depth) was favored by channelized fluid flow, high devolatilization rates, and permeabilities similar to those found in geothermal reservoirs ( k ~ 10 − 16 to 10 − 13 m 2 ). For channelized flow, magma-induced thermal pressurization alone can generate cases of ∆ P ≥ 0.1 bars for all permeabilities in the range 10 − 16 to 10 − 13 m 2 , whereas in radial flow regimes thermal pressurization causes ∆ P < 0.1 bars for all permeabilities. Changes in distal fluid pressure occurred before proximal pressure changes given modest anisotropies ( k h / k v ~ 10–100). Invoking k ( z , T , P ) and high, sustained devolatilization rates caused large dynamic fluctuations in k and P in the near-magma environment but had little effect on pressure changes at the distal dVT location. Intrusion below the brittle-ductile transition damps but does not prevent pressure transmission to the dVT site.

Journal of Volcanology and Geothermal Research↗

USGS2026

Can a physics-based fire behavior model predict burn severity and post-fire debris flow hazard?

Post-wildfire flooding and debris flows pose a threat to ecosystems and infrastructure. As the frequency and size of severe wildfires increase, there is a growing need for predictive debris flow hazard modeling in unburned but fire-prone landscapes. However, no crosswalks have been established between physics-based models of fire behavior and the remotely sensed burn severity inputs needed for debris flow modeling. Here, we explore this connection by creating a large ensemble of fire behavior simulations using QUIC-Fire – a fast-running three-dimensional fire model – in basins across five major-disaster wildfires in California and Washington. We then conducted generalized linear mixed modeling (GLMM) and Random Forest modeling to identify QUIC-Fire outputs that best predicted remotely sensed burn severity. We found that the selected linear predictors exhibited weak and unintuitive relationships with soil burn severity and differenced Normalized Burn Ratio. Moreover, model evaluations of GLMMs and Random Forests showed they were poorly predictive of both metrics of remotely sensed burn severity. This suggests that there are barriers to integrating QUIC-Fire into a wildfire-hazards modeling framework, which could include mismatches in scale between fire models and remotely sensing, or simplifications of weather data and fire initiation. Further research is needed to establish connections between the fire physics outputs of QUIC-Fire and post-fire effects to vegetation and soils that are essential for evaluating debris flow hazard.

California, Washington↗

USGS2016

Synthesis of common management concerns associated with dam removal

Managers make decisions regarding if and how to remove dams in spite of uncertainty surrounding physical and ecological responses, and stakeholders often raise concerns about certain negative effects, regardless of whether or not these concerns are warranted at a particular site. We used a dam-removal science database supplemented with other information sources to explore seven frequently-raised concerns, herein Common Management Concerns (CMCs). We investigate the occurrence of these concerns and the contributing biophysical controls. The CMCs addressed are: degree and rate of reservoir sediment erosion, excessive channel incision upstream of reservoirs, downstream sediment aggradation, elevated downstream turbidity, drawdown impacts on local water infrastructure, colonization of reservoir sediments by non-native plants, and expansion of invasive fish. Biophysical controls emerged for some of the concerns, providing managers with information to assess whether a given concern is likely to occur at a site. To fully assess CMC risk, managers should concurrently evaluate site conditions and identify the ecosystem or human uses that will be negatively affected if the biophysical phenomenon producing the CMC occurs. We show how many CMCs have one or more controls in common, facilitating the identification of multiple risks at a site, and demonstrate why CMC risks should be considered in the context of other factors like natural watershed variability and disturbance history.

Journal of the American Water Resources Associatio↗