Geology topics

Research about Canada, United States

Source-linked reports with geographic coverage including Canada, United States.

At least 91 records · Page 5Linked to original sources

USGS2025

Structural controls on splay fault rupture dynamics during Cascadia megathrust earthquakes

Great subduction earthquakes ( M w ≥ 8.0) can generate devastating tsunamis by rapidly displacing the seafloor and overlying water column. These potentially tsunamigenic seafloor offsets result from coseismic fault slip and deformation beneath or within the accretionary wedge. The mechanics of these shallow rupture phenomena and their dependence on subduction zone properties remain unresolved, partly due to the sparsity of offshore observations of shallow megathrust earthquake deformation. Here, we analyze how offshore structure influences shallow rupture mechanics and slip partitioning using 3D dynamic earthquake simulations of the Cascadia subduction zone (CSZ) megathrust with and without variably dipping seaward- or landward-vergent splay faults in the wedge that sole into the megathrust. Resulting tradeoffs between splay and megathrust slip reveal structural controls on rupture partitioning, with greater splay slip leading to less shallow megathrust slip updip. Gently dipping and seaward-vergent splays host more slip than those with steeper, landward-vergent splays. To isolate the underlying mechanisms, we compare models with Andersonian and plunging principal stresses. Results suggest distinct static and dynamic processes control the dip- and vergence-dependence of splay rupture: static (mis)alignment relative to far-field tectonic loading favors slip on more optimally oriented, shallowly dipping splay faults. In contrast, dynamic stress interactions of an updip-propagating megathrust rupture front with the free surface and potential branch faults favor forward branching onto seaward-vergent splays and inhibit backward branching onto landward-vergent splays. Resulting seafloor displacements suggest splay fault structure may influence coseismic tsunami source processes, highlighting the importance of dynamically viable rupture scenarios in subduction hazard assessments.

Cascadia subduction zone↗

USGS2025

Revised length categories and standard weight equation for Northern Pikeminnow

Objective Length and weight indices (e.g., proportional size distribution, relative weight) provide standardized benchmarks that are useful for comparing groups of fish, identifying ecological interactions, and evaluating the effect of management actions. However, the current length categories and standard weight ( W s ) equation for Northern Pikeminnow Ptychocheilus oregonensis , a species of important management focus in the Pacific Northwest, were developed using limited data, unclear methods, and no validation. As such, we sought to revise the length categories and W s equation for Northern Pikeminnow. Methods We used the all-tackle world record to develop length categories for Northern Pikeminnow. We compiled data from 100,663 Northern Pikeminnow from 114 populations in Idaho, Montana, Oregon, Washington, and British Columbia to develop a revised W s equation. Most fish were measured in fork length (FL), so we converted total lengths (TLs) to FLs using the equation TL = −2.301 + 0.916(FL); r 2 = 0.998. We used the regression line percentile, linear empirical percentile (EmP), and quadratic EmP methods to develop 50th percentile and 75th percentile W s equations. We then assessed length-related biases in our W s equations and the previously published equation. Results We propose minimum FLs of 17 cm (7 inches; stock), 26 cm (10 inches; quality), 35 cm (14 inches; preferred), 41 cm (16 inches; memorable), and 51 cm (20 inches; trophy) for proportional size distribution calculations. The previously published W s equation exhibited substantial length-related biases according to the Willis and weighted empirical quartile method tests. The EmP 50th-percentile W s equation was the only equation that we evaluated that did not exhibit length-related bias. The EmP 50th-percentile W s equation was the best performing equation (i.e., no length-related bias detected). The equation is log 10 ( W s ) = −5.258 + 3.135(FL), where W s is in grams and FL is in millimeters. The equation is valid for Northern Pikeminnow 90–580 mm FL. Conclusion The length categories, length conversion, and W s equation presented here will aid fisheries professionals in the study and management of Northern Pikeminnow populations.

British Columbia, Idaho, Montana, Oregon Washingto↗

USGS2025

Rapid fault healing from cementation controls the dynamics of deep slow slip and tremor

Despite its status as one of the most important discoveries in geophysics, the physical mechanism(s) responsible for slow slip events (SSEs) are not well understood. Here, we synthesize observations of deep SSEs in the Cascadia Subduction Zone and argue that rapid, cohesive fault strengthening may control the dynamics of deep SSEs. Cohesive strength is frequently ignored in constitutive laws used to describe fault rheology in numerical simulations of earthquakes and SSEs alike. To demonstrate its importance, we perform and analyze a suite of petrological experiments that simulate fault healing under representative pressure and temperature conditions. We show that significant cohesive strength recovery caused by dissolution-precipitation processes occurs on timescales of just a few hours. Together, our experimental and observational results support the idea that cohesion is a key component of fault strength under SSE conditions and highlight the need for its inclusion in both future experiments and numerical models of fault slip.

British Columbia, Washington↗

USGS2025

Population genetics of the endangered narrowly endemic Island Marble butterfly (Euchloe ausonides insulanus)

The Island Marble butterfly ( Euchloe ausonides insulanus ) is an endangered species endemic to the San Juan Islands off the coast of Washington State, United States, and British Columbia, Canada. The species was thought to be extinct for ~ 90 years before it was rediscovered at American Camp, San Juan Island National Historical Park in 1998. Here, we report the results of the first population genetic analyses for insulanus , using DNA collected non-invasively from individuals in the last known stronghold for the species. We used DNA extracted from meconium, larval exuviae, and natural mortalities to generate and test thirteen new microsatellite markers to estimate genetic diversity, population structure, and kinship. We assembled and annotated mitochondrial genomes, which were used alongside museum specimens of insulanus collected ~ 100 years ago from Vancouver Island, and other members of the E. ausonides species complex, to infer the evolutionary history of the species. The results indicated that insulanus experiences low heterozygosity, a small effective population size (N e ), and low allelic diversity. High levels of inbreeding were found in some individuals, but inbreeding was uneven across the population. No population structure or partitioning of genetic variation by host plant was detected. The mitogenomes of extant insulanus were all identical and modern samples showed a loss of allelic diversity compared to insulanus from museums. Extant insulanus formed a clade with museum specimens and we identified multiple putatively diagnostic alleles to differentiate insulanus from other subspecies. Based on these results, we outline considerations for species management and genetic monitoring.

British Columbia, Washington↗

USGS2025

Same view through a different lens: Comparing population trends for North American birds using eBird and the Breeding Bird Survey

Confidently estimating population trends is of vital importance for a wide range of ecological, conservation, and management applications. North America has 2 major data sources for estimating population trends of breeding birds—the North American Breeding Bird Survey (BBS) and the global participatory science project eBird. Because the surveys differ in protocols, coverage, and data analysis, their trend estimates are expected to vary in magnitude, direction, and/or precision for at least some species and regions. Here, we compare independently derived estimates of population change between 2012 and 2022 for every combination of species and bird conservation region (BCR) covered by both surveys ( n = 5,577 combinations) as well as aggregated across entire ranges or within the U.S. or Canada. Uncertainty was substantial for both surveys, though more prevalent for BBS (81% of credibility intervals for estimates included zero) than eBird (34% of confidence intervals overlapped zero). We found agreement of trend directions between the 2 surveys. Only 1.3% of estimated trends were significant in opposite directions between the 2 surveys for all species/BCR combinations, with the median difference in trend magnitude being –0.02% (BBS minus eBird trend). Correlations between the 2 were strongest for estimates that were graded as being high credibility compared to estimates judged to have medium or low credibility. Both surveys were subject to species, taxonomic, and regional effects that influenced agreement. Overall, we show where trend estimates derived from BBS and eBird agree, explore where they diverge, present several comparisons to assist in interpreting results from both surveys, and inform efforts to integrate information from each.

Ornithological Applications↗

USGS2025

Lake sturgeon behavioral diversity in the Laurentian Great Lakes: Migratory patterns across populations and habitats

Background Characterizing the diversity of migration behaviors from the individual to the population level is essential for understanding how organisms respond to environmental variation and how these responses affect survival and habitat use. Lake sturgeon ( Acipenser fulvescens ) is a species of special concern in the Laurentian Great Lakes that are long-lived and generally classified as intermittent, adfluvial spawners. Observations of lake sturgeon movements at ecologically relevant spatiotemporal scales have shown that migration behavior often varies among individuals within the same population. However, studies on individual populations, particularly when focused only on a part of the life cycle (e.g., often spawning), provide just a partial understanding of the species’ full migratory scope and processes underlying expression of different migratory behaviors. To better understand lake sturgeon migratory diversity, we characterized and compared migratory behaviors of six lake sturgeon populations occupying environments that varied in habitat availability and connectivity in different Laurentian Great Lakes. Methods Sequence analysis combined with agglomerative hierarchical clustering and visual inspection of daily location data were used to identify distinct lake sturgeon migratory behaviors present in each population. Results Seven distinct migratory behaviors were identified based on differential patterns of lake and river use that encompass spawning and other seasonal periods. Behaviors were categorized as annual spring river, intermittent spring river, intermittent two-step, annual summer river, annual winter river, and annual interlake migrants along with river residents. The presence and frequency of migratory behaviors varied substantially among populations. Conclusions Our study demonstrated that migratory diversity is a general feature of lake sturgeon life history that may be partially shaped by habitat availability and connectivity. Given these results, we propose a conceptual model that links habitat availability and connectivity to migratory diversity and predict a positive association between them. This updated framework provides a cohesive basis for understanding lake sturgeon migratory behavior across variable ecological contexts in the Laurentian Great Lakes and will help promote future research to refute or refine the model.

Laurentian Great Lakes↗

USGS2025

Quantifying benthic flux of Mysis biomass through diel vertical migration at the ecosystem scale

Mysis diluviana is a macroinvertebrate that couples benthic and pelagic habitats on a daily timescale through diel vertical migration (DVM). However, quantifying how much Mysis biomass is exchanged between benthic and pelagic habitats at an ecosystem scale is difficult because of sampling limitations and variability in Mysis DVM behavior related to light and depth. Although Mysis are benthic-pelagic migrators, a portion remains pelagic during the day offshore in Lake Ontario, partially contradicting the assumption of population-level DVM over deep areas. To estimate the amount of biomass transferred from benthic to pelagic habitat via DVM in Lake Ontario, we estimated the portion of pelagic biomass at night originating from benthic habitat as the difference between night and day pelagic estimates from net tows along a bathymetric depth gradient. We then modeled the portion as a function of depth, extrapolated these depth-dependent estimates to an existing lake wide night-pelagic dataset, and summed amounts across depth strata. We estimated more biomass was transferred from benthic to pelagic habitat at intermediate lake depths (100–160 m) despite greater offshore (> 180 m) night-pelagic biomass. Our results suggest ways to improve estimates of Mysis habitat coupling and how to account for important factors such as depth and light for modeling Mysis DVM behavior at the population- and ecosystem-levels.

Lake Ontario↗

USGS2025

Non-native prey availability and over-compensatory density dependence drive population dynamics of a native fish predator

Understanding the factors that regulate population dynamics is crucial for conserving imperiled species. Bull trout ( Salvelinus confluentus ), a piscivorous salmonid and one of North America's most threatened cold-water species, has declined significantly due to habitat loss, overfishing, invasive species, and climate change. While recovery efforts have primarily targeted these threats, the role of prey availability in influencing bull trout population dynamics under multiple stressors remains poorly understood. Using a stage-based integrated population model, we quantified the effects of non-native prey availability (kokanee; Oncorhynchus nerka ), angling pressure, climatic variation, and density-dependent processes on bull trout population dynamics in Lake Koocanusa, a transboundary reservoir and river system (United States and Canada), over a 40-year period (1980–2023). Our results show that bull trout populations are regulated by density-dependent processes, including over-compensation in sub-adult recruitment and reduced adult survival at high densities. Increased kokanee biomass and restricted harvest significantly enhanced bull trout survival and abundance, whereas reduced water availability had a limited negative effect on sub-adult production. Model simulations indicate that as kokanee biomass availability increases, the number of bull trout that can be sustainably harvested also increases. In fact, a modest annual fishery (300 individuals) can be sustained, especially under moderate to high kokanee biomass conditions. These results underscore the importance of prey availability, including non-native species, in supporting bull trout populations. Effective management of threatened apex fish predators like bull trout requires addressing the complex interplay between environmental threats, prey dynamics, and density-dependent mechanisms across all life stages.

British Columbia, Montana↗

USGS2025

Application of a data-limited stock assessment model to understand the historical population of cisco in Lake Erie

Historically, cisco (Coregonus artedi) supported one of the largest commercial fisheries in Lake Erie with over 25 million kilograms landed in a single year during the fisheries’ peak in 1912. Like many other coregonines in the Laurentian Great Lakes, the Lake Erie cisco population eventually collapsed and is now considered extirpated. Management interest in restoring coregonine abundance and diversity across the Great Lakes motivated the development of a multi-agency, science-based framework aimed at restoring or enhancing populations. As part of that framework, we applied a data-limited model (depletion-based stock reduction analysis) to a reconstructed catch series to estimate historical levels of abundance of cisco in Lake Erie. Estimated historical carrying capacity was 171 million kilograms (66.5 kg/ha) and estimated annual exploitation was as high as 70 %. Model projections suggested the population should have rebuilt following a fishery collapse in 1925; however, the fishery remained depressed and cisco disappeared from Lake Erie despite the end of targeted exploitation in 1965. Other factors, such as habitat degradation and invasive species, are believed to have played a part in impeding recovery of cisco. The results of this study could be used in setting potential restoration targets for Lake Erie cisco given current conditions. The likelihood of achieving recovery targets based on historical population size will be dependent on the suitability of current conditions for sustaining the population.

Lake Erie↗

USGS2025

Unique thermal mixing patterns in Lake Ontario revealed by novel year-round observations of thermal stratification

Year-round records of thermal stratification in the Great Lakes are rare, and there are few observations of thermal stratification during winter. In this paper, we analyze temperature data from 13 temperature logger chains and from over 130 benthic acoustic receivers that were deployed across Lake Ontario for 2 yr. The timing and duration of the fall overturn correlate with the local average water depth, and shallow sites (< 50 m depth) overturn up to a month before deep sites (> 100 m depths). Likewise, in spring, the shallow sites warm faster. Lake Ontario has partial ice cover, so wind-driven mixing stirs the water column throughout winter, and inverse thermal stratification is largely absent. The depth-averaged winter water temperatures vary between 0°C and 4°C, with the coldest temperatures (near 0.1°C) found in the shallow Kingston basin and warmest temperatures (near 4°C) at sites near the 244 m deep Rochester Basin. Lake Ontario appears to be a warm monomictic lake, rather than having a dimictic mixing pattern as previously described—there is no sustained ice cover or inverse stratification that inhibits vertical mixing in winter. Winter is a poorly understood season for many aquatic processes, including fish bioenergetics, fish distribution, biochemical processes, invertebrate distribution, and production. Moreover, the lack of knowledge of winter has hampered the use of correct initial conditions for running large lake hydrodynamic models.

Lake Ontario↗

USGS2025

Quantifying the relative importance of survival threats to a long-lived reptile using expert elicitation

Long-term survival of a conservation-reliant species requires understanding the impact of threats on population growth rate and the management actions that can help mitigate these threats. We used a threat assessment with expert-elicited estimates to determine the relative effect of each stage-specific threat on the population growth rate of the wood turtle Glyptemys insculpta . In addition, we offered potential management actions that could mitigate these threats and examined the relative cost and benefit of each. The experts responded that predators had the largest effect on hatchling and juvenile survival and that road mortality had the largest effect on adult survival. The population growth rate of the simulated turtle population increased the most when predators were removed from the system, though the population trajectory remained negative. Finally, we found that predator control had the lowest cost:benefit ratio of the proposed management actions. The process used in this analysis of expert elicitation combined with modeling that accounts for uncertainty proved to be a useful technique that is less expensive and labor intensive than empirical studies and quicker to implement, although it relies on sufficient empirical studies to inform expert responses. This process could be replicated for other species to inform species status assessments.

Endangered Species Research↗

USGS2025

Ice Age biogeography corresponds with current climate vulnerability of freshwater fishes

1. Both local environmental factors and historical biogeography shape ecological communities, but determining which historical biogeographical patterns correspond with contemporary climate vulnerability is an underused conservation method. The historical colonization patterns of freshwater fishes following the Pleistocene (“Ice Age”) glaciations offers an ideal model for comparing historical biogeography and climate-change vulnerability. 2. We used current thermal niches and future stream-temperature projections to estimate the climate vulnerability of 29 Great Plains and Rocky Mountain fishes that we classified as either early or late colonists of the region in the wake of glacial retreat (~19,000 years ago). 3. Ninety-three percent of the most vulnerable species were amongst the earliest colonists of the region and we consider them “postglacial-pioneer species”. Median predicted site loss (number of historically occupied sites predicted to become too warm by end-of-century) was 0% for late colonizing species and 33% for early colonizing species. 4. We provide empirical evidence that postglacial-pioneer fishes are uniquely vulnerable to climate change, and we suggest this may apply to many taxa from formerly glaciated regions. More broadly, we demonstrate that evaluating the relationship between current species-environment patterns and historical biogeography may be a fruitful avenue for future climate change and conservation research.

North America↗

USGS2025

Three-dimensional seismic velocity model for the Cascadia Subduction Zone with shallow soils and topography, version 1.7

The U.S. Geological Survey’s seismic velocity model for the Cascadia Subduction Zone provides P- and S-wave velocity ( V P and V S , respectively) information from 40.2° to 50.0° N. latitude and −129.0° to −121.0° W. longitude, and is used to support a variety of research topics, including three-dimensional (3D) earthquake simulations and seismic hazard assessment in the Pacific Northwest. This report describes an update to the previous version (v) 1.6 of the 3D seismic velocity model for the Cascadia Subduction Zone. This new model (herein referred to as v1.7) contains more detailed near-surface structure for improved earthquake ground motion modeling. Updated features include the addition of a new shallow soil velocity model in the top few hundred meters and the option of adding user-specified topography. Although v1.6 of the Cascadia seismic velocity model has a minimum V S of 600 meters per second (m/s), the new model (v1.7) has a minimum V S of approximately 40 m/s. Overall, this update will allow for more accurate ground motion estimates from 3D simulations of scenario earthquakes in the Cascadia Subduction Zone region.

British Columbia, California, Oregon, Washington↗

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↗

USGS2025

Population structure and movement dynamics of Redband Trout in the Kootenai River basin

Redband Trout Oncorhynchus mykiss gairdneri is a species of high conservation concern in the Kootenai River basin, United States and Canada. However, identifying the specific mechanisms influencing Redband Trout is difficult as its life history structure is largely unknown. Using otolith microchemistry analysis, we found that 18.5% ( n = 66) of the fish sampled ( n = 329) were a product of mainstem Kootenai River spawning. The remainder were fish ( n = 264) that originated from one of 31 tributaries throughout the basin. Redband Trout were captured between 0 and 153 km from their assigned natal tributaries. Most migratory fish originated from tributaries in the uppermost portion of the study area and had migratory mothers. However, a diversity of life history types (e.g., multi-year resident, migratory) was observed in the population. Redband Trout were estimated to outmigrate from natal habitats at a variety of ages (0–3 years) and lengths (12–430 mm). Our results highlight diversity in the population structure of Redband Trout and provide information valuable to conservation efforts in the Kootenai River system. Interestingly, several tributaries identified as natal habitats of Redband Trout contain substantial obstacles originally considered barriers to upstream fish passage. In addition, some streams in the lower portion of the Kootenai River basin with degraded habitat that were not expected to contribute Redband Trout had natal assignments. These observations indicate previously unknown sources of Redband Trout and support the notion that improved fish passage in disconnected and degraded streams could increase production of Redband Trout throughout the Kootenai River basin. Restoration efforts in the mainstem that include considerations for spawning and natal habitat may also aid in conservation efforts for Redband Trout. Our research demonstrates the value of microchemistry in planning and evaluating habitat restoration to recover fishes in large rivers.

British Columbia, Idaho, Montana↗

USGS2025

Field test of the self-preserving eDNA filter and the importance of calibration when changing methods during long-term monitoring

Filtering water is currently the primary field method used for collecting aquatic environmental DNA (eDNA). One of the drawbacks of filtering is the need to transfer the filter from the filter housing to a preservative-filled container in the field. New products are being developed to avoid this handling step, but comparative studies are needed to ensure that the results produced by new protocols are transferable within and across eDNA monitoring programs. To meet this need, we evaluated two filter preservation methods (self-preserving filter housing vs. ethanol) of the 5.0-μm polyethersulfone (PES) filter membrane in a field trial typical of stream fisheries eDNA sampling. We compared DNA detection and yield for free-swimming rainbow trout, Oncorhynchus mykiss (Walbaum, 1792), from streams in Washington, United States, and British Columbia, Canada, while accounting for the effects of two environmental covariates: stream discharge and water temperature. As these streams were part of an ongoing fisheries eDNA monitoring program, we also compared these methods to the original protocol, which used a 0.45-μm cellulose nitrate (CN) filter membrane and ethanol preservative. We found that the self-preserving filter housings collected and preserved eDNA well and provided similar results to identical filters preserved in ethanol. The 5.0-μm PES filters preserved in ethanol significantly outperformed the original protocol in terms of both DNA detection and yield, highlighting the importance of calibration of eDNA results when changing sampling methods during an ongoing monitoring program.

british Columbia, Washington↗

USGS2025

Movements and habitat use vary across the Rocky Mountain Population of trumpeter swans

The Rocky Mountain Population (RMP) of trumpeter swans Cygnus buccinator (hereafter, swans) in North America includes breeders in the Greater Yellowstone Area (GYA) and other western states (together, United States segment) and western provinces of Canada (Canada segment). Conservation concern for the United States segment stems from its slow population growth and the resident nature of GYA swans, which intermingle with migrating Canada segment swans in wintering habitats. Thus, understanding variation in migratory behavior and habitat use by swans in the two population segments can inform how management actions may affect the RMP. We used telemetry data from 55 RMP swans captured in the western United States to understand their movements and habitat use. For 45 swans (60 swan-years) that spent the summer in the United States, distance traveled between breeding and wintering areas ranged from 0 km (i.e., no migration in 22% of swan-years) to 473 km, with an average of nonzero movements of 118 ± 95 km (SD). Swans traveled farther distances when maximum temperatures were lower. For 10 swans (16 swan-years) that spent the summer in Canada, five appeared to molt but not to nest, and four appeared to nest in one or more years. Migration timing was similar for molting and nesting swans. All five molting swans and one nesting swan spent at least one previous summer in the GYA. Migratory connectivity of all birds was weaker in years when more swans migrated to Canada for the summer. During the breeding season, Canada swans used low-elevation lakes, but United States swans used high-elevation lakes. Both groups of swans increased use of crop fields outside of the breeding season. Our study shows interchange between the United States and Canada segments, a finding that challenges the efficacy of existing population designations. Furthermore, variation in movement behavior of GYA swans suggests possible actions, such as restoring winter habitats to increase swan distribution and migration, to support swan conservation.

Rocky Mountains↗

USGS2025

Seasonal synchronicity and multi-decadal stability of headwater biogeochemistry in the northern temperate zone

Temporal patterns in chemistry of headwater streams reflect responses of water and elemental cycles to perturbations occurring at local to global scales. We evaluated multi-scale temporal patterns in up to 32 y of monthly observations of stream chemistry (ammonium, calcium, dissolved organic carbon, nitrate, total dissolved phosphorus, and sulfate) in 22 reference catchments within the northern temperate zone of North America. Multivariate autoregressive state-space (MARSS) models were applied to quantify patterns at multi-decadal, seasonal, and shorter intervals during a period that encompassed warming climate, seasonal changes in precipitation, and regional declines in atmospheric deposition. Significant long-term trends in solute concentrations within a subset of the catchments were consistent with recovery from atmospheric deposition (e.g., calcium, nitrate, sulfate) and increased precipitation (e.g., dissolved organic carbon). Lack of evidence for multi-decadal trends in most catchments suggests resilience of northern temperate ecosystems or that subtle net effects of simultaneous changes in climate and disturbance regimes do not result in directional trends. Synchronous seasonal oscillations of solute concentrations occurred across many catchments, reflecting shared climate and biotic drivers of seasonality within the northern temperate zone. Despite shared patterns among catchments at a seasonal scale, multi-scale temporal patterns were statistically distinct among even adjacent headwater catchments, implying that local attributes of headwater catchments modify the signals imparted by atmospheric phenomena and regional disturbances. To effectively characterize hydrologic and biogeochemical responses to changing climate and disturbance regimes, catchment monitoring programs could include multiple streams with contributing areas that encompass regional heterogeneity in vegetation, topography, and elevation. Overall, detection of long-term patterns and trends requires monitoring multiple catchments at a frequency that captures periodic variation (e.g., seasonality) and a duration encompassing the perturbations of interest.

Biogeochemistry↗