Western Montana Research Articles

Late Paleogene paleotopographic evolution of the northern Cordilleran orogenic front: Implications for demise of the orogen

The paleotopographic history of the North American Cordilleran orogen holds the key to understanding mechanisms of orogenesis and subsequent orogenic collapse. It has been suggested that the orogenic front in western Montana (USA) and Alberta (Canada) was more than 4 km high during Late Cretaceous−early Eocene contractional deformation and during the initial phase of extension in the middle Eocene; however, the late Eocene−Oligocene topographic evolution during continued extensional collapse remains poorly constrained. Here we extend the paleotopographic record in the Kishenehn Basin in northwestern Montana and southeastern British Columbia (Canada) to the late Oligocene by studying δ18O values of fossil mollusks and cement and paleosol carbonates. The molluscan taxa changed from three sympatric groups with preferred habitats ranging from tropical wet, semi-arid subtropical, and temperate during the middle and late Eocene, to mainly a single group associated with temperate environment during the Oligocene, reflecting a decline in molluscan biodiversity induced by climate cooling across the Eocene−Oligocene transition. Reconstructed δ18O values of alpine snowmelt and basinal precipitation decreased by 1.4‰ and 3.8‰, respectively, from the middle to late Eocene, reflecting climate cooling and ∼1 km surface uplift of the basin floor. The reconstructed alpine snowmelt δ18O values then increased by 2.9‰ in the Oligocene suggesting a ∼0.5 km drop in elevation of the orogenic front. Collectively, the results of our new and previously published δ18O data chronicle the paleotopographic response to the change from flat-slab subduction to slab rollback over a 45 m.y. period. These data suggest that the orogenic front was characterized by high elevation (>4 km) in the ancestral Lewis-Clark-Livingston ranges during latest Cretaceous−early Eocene (ca. 75−52 Ma) contraction. The initial phase of extension related to the Kishenehn Basin created a lowland basin with a surface elevation of only ∼1.5 km during the early middle Eocene (ca. 46−44 Ma) whereas the ranges remained >4 km high. The high range elevations were sustained for at least 12 m.y. in the middle to late Eocene concurrent with extension, while the basin floor elevation was uplifted to ∼2.5 km by the latest Eocene (ca. 36−34 Ma). Basin aggradation can explain at most half of the 1 km basin floor uplift. The remaining amount (at least 0.5 km) and sustained high range elevation suggest that range denudation and crustal extension was compensated by the isostatic and thermal effects of slab rollback and/or passage of a slab window and infusion of hot asthenosphere beneath the continent. The range elevation in the orogenic front decreased ∼0.5 km by the late Oligocene (ca. 28 Ma), associated with a decrease in rock uplift rate associated with extension. A post-Oligocene elevation drop of ∼1 km resulted in both the ranges and basin floor reaching modern topography in the Kishenehn Basin drainage, likely due to the regional effect of Neogene Basin and Range extension. This study, along with the previous investigation of the Kishenehn Basin by Fan et al. (2017), are the first studies that systematically investigate paleorelief of the orogenic belt by reconstructing paleoelevations of the mountains and the basin at the same time. The results highlight that the Cordilleran orogenic front of northern Montana and southern British Columbia sustained its high elevation edifice for at least 12 m.y. after the start of extension. We suggest that initial crustal extension did not result in orogenic demise because of concurrent thermal and isostatic uplift.

A Framework for Collaborative Wolverine Connectivity Conservation

In the Rocky Mountain West, maintaining connectivity between high-elevation public lands is important for wolverines and other species of conservation concern and depends on building collaborative partnerships with key landowners. However, prioritizing conservation areas that benefit wolverine connectivity and identifying key partners presents several challenges: 1) the wolverine metapopulation functions over a large spatial extent and 2) anthropogenic changes such as housing development and climate change may shift connective pathways for wolverines in the future. We prioritized areas important for wolverine connectivity under future climate conditions across western Montana using a systematic conservation planning framework. We solved our planning problem by optimizing 10, 15, 20, and 50% of habitat features that mitigated human threats and maintained ecological processes for wolverines using integer linear programming. We identified 369 privately owned areas in the 10% solution, 572 privately owned areas in the 15% solution, 822 privately owned areas in the 20% solution, and 3,996 privately owned areas in the 50% solution where voluntary landowner easements would greatly improve the long-term landscape functionality for wolverine connectivity. The median estimated easement value ranged from $8,762 to $12,220 across the four solutions, and the total costs ranged from $14,874,371 to $196,346,714 if all landowners agreed to voluntary conservation easements. This effort demonstrates the utility of optimization problems for conserving connectivity and provides a mutually beneficial and proactive tool to engage potential collaborators in wolverine connectivity conservation.

High egg rejection and low aggressive responses of Mountain Bluebirds (Sialia currucoides) toward Brown-headed Cowbirds (Molothrus ater)

Hosts of brood parasitic Brown-headed Cowbirds (Molothrus ater) have acquired a suite of behavioral responses against them, including egg rejection and nest defense. Several cavity-nesting species have been shown to exhibit egg rejection, including the Mountain Bluebird (Sialia currucoides). Previous work reported egg rejection rates as low as 20%, under the assumption that they may instead be relying heavily on aggressive interactions to defend their nests from parasitism. In the present study, we evaluated the egg rejection and aggression behaviors of a population of Mountain Bluebirds in western Montana. Results indicate a higher rate of rejection than previously documented, with model cowbird eggs rejected in over 40% of trials (n = 37), compared with 5% rejection of mimetic eggs. In contrast, we found little aggressive response by Mountain Bluebirds toward both an adult cowbird model and a nonthreatening American Goldfinch (Carduelis tristis) model. Aggression patterns were consistent between sexes, with the most common behaviors including sitting and watching from within 5 m and flying past the models (n = 21). Hovering above the model was observed on a few occasions, and physical contact by bluebirds happened on only 3 occasions toward the goldfinch model. In summary, Mountain Bluebirds in western Montana showed relatively high rates of cowbird egg rejection but little aggressive response toward adult heterospecific nest intruders.

Ultrasound imaging identifies life history variation in resident Cutthroat Trout.

Human activities that fragment fish habitat have isolated inland salmonid populations. This isolation is associated with loss of migratory life histories and declines in population density and abundance. Isolated populations exhibiting only resident life histories may be more likely to persist if individuals can increase lifetime reproductive success by maturing at smaller sizes or earlier ages. Therefore, accurate estimates of age and size at maturity across resident salmonid populations would improve estimates of population viability. Commonly used methods for assessing maturity such as dissection, endoscopy and hormone analysis are invasive and may disturb vulnerable populations. Ultrasound imaging is a non-invasive method that has been used to measure reproductive status across fish taxa. However, little research has assessed the accuracy of ultrasound for determining maturation status of small-bodied fish, or reproductive potential early in a species' reproductive cycle. To address these knowledge gaps, we tested whether ultrasound imaging could be used to identify maturing female Westslope Cutthroat Trout (Oncorhynchus clarkii lewisi). Our methods were accurate at identifying maturing females reared in a hatchery setting up to eight months prior to spawning, with error rates ≤ 4.0%; accuracy was greater for larger fish. We also imaged fish in a field setting to examine variation in the size of maturing females among six wild, resident populations of Westslope Cutthroat Trout in western Montana. The median size of maturing females varied significantly across populations. We observed oocyte development in females as small as 109 mm, which is smaller than previously documented for this species. Methods tested in this study will allow researchers and managers to collect information on reproductive status of small-bodied salmonids without disrupting fish during the breeding season. This information can help elucidate life history traits that promote persistence of isolated salmonid populations.

On the timing and metallogenic implications of the sediment-hosted stratiform copper–silver mineralization in the Creston Formation (Belt-Purcell Supergroup), British Columbia, Canada

Sediment-hosted stratiform copper–silver mineralization is common in Middle Creston strata of southeastern British Columbia, a correlative lithostratigraphic unit of the Ravalli Group of western Montana. The Ravalli Group hosts the largest sediment-hosted stratiform copper–silver deposits of the Belt-Purcell Supergroup. Strata of the Belt-Purcell Supergroup accumulated between ~1540 and 1300 Ma in a pericratonic rift basin (Belt-Purcell basin) constructed over Archean and Paleoproterozoic blocks of Laurentia. The basin underwent several events of deformation and metamorphism at ~1380–1325 Ma (East Kootenay orogeny), ~1200–1000 Ma (Grenvillian), and 900–800 Ma (Goat River orogeny) prior to Meso-Cenozoic Cordilleran tectonism. U–Pb geochronology of detrital zircons from Middle Creston rocks define a maximum age of deposition of ~1470 Ma, in agreement with the established age range for the Belt-Purcell Supergroup in southeastern British Columbia. Molybdenite is locally part of the sediment-hosted stratiform copper–silver mineralization in the region, either as isolated grains within mineralized strata or in paragenetic sequence with the cupriferous sulfides. Re–Os geochronology on molybdenite reveals repetitive episodes of fluid flow through Middle Creston strata in the study area, with principal events at 1170–1140 and 1043–990 Ma. Paragenetically constrained molybdenite, in growth sequence with bornite and chalcocite, yields an age of 1043 ± 6 Ma. Hence, the fluid responsible for the sediment-hosted stratiform copper–silver mineralization in southeastern British Columbia permeated the host sequence ~440 m.yr. after its deposition, coincident with the Grenville-age metamorphic event at ~1200–1000 Ma. Such tectono-thermal activity extended along the western margin of Laurentia beyond the limits of the Belt-Purcell basin, and was a far-field response to orogen-scale collisions taking place off Laurentia during Rodinia assembly.

Characterization of CYP3A pharmacogenetic variation in American Indian and Alaska Native communities, targeting CYP3A4*1G allele function.

The frequencies of genetic variants in the CYP3A4 and CYP3A5 genes differ greatly across global populations, leading to profound differences in the metabolic activity of these enzymes and resulting drug metabolism rates, with important consequences for therapeutic safety and efficacy. Yet, the impact of genetic variants on enzyme activity are incompletely described, particularly in American Indian and Alaska Native (AIAN) populations. To characterize genetic variation in CYP3A4 and CYP3A5 and its effect on enzyme activity, we partnered with AIAN people living in two regions of Alaska: Yup’ik Alaska Native people living in the Yukon‐Kuskokwim Delta region of rural southwest Alaska and AIAN people receiving care at the Southcentral Foundation in Anchorage, Alaska. We identified low frequencies of novel and known variation in CYP3A4 and CYP3A5, including low frequencies of the CYP3A4*1G and CYP3A5*1 variants, and linkage disequilibrium patterns that differed from those we previously identified in an American Indian population in western Montana. We also identified increased activity of the CYP3A4*1G allele in vitro and in vivo. We demonstrated that the CYP3A4*1G allele confers increased protein content in human lymphoblastoid cells and both increased protein content and increased activity in human liver microsomes. We confirmed enhanced CYP3A4‐mediated 4β‐vitamin D hydroxylation activity in Yup’ik people with the CYP3A4*1G allele. AIAN people in Alaska and Montana who carry the CYP3A4*1G allele—coupled with low frequency of the functional CYP3A5*1 variant—may metabolize CYP3A substrates more rapidly than people with the reference CYP3A4 allele.

Providing APPE pharmacy students rural health assessment experience following wildfire event in western Montana

Background and purposeWe describe a novel, interprofessional, experiential training involving pharmacy students in response to a health emergency in rural Montana (MT). Educational activity and settingFourth-year pharmacy students on clinical rotations were recruited to participate in screening events assessing effects of wildfire smoke in Seeley Lake, MT. Students were required to fulfill at least two hours of supplementary training in addition to education on human research guidelines. Students assisted with patient surveys (demographics, health, and respiratory), physiological testing with biomedical researchers, blood pressure and medication counseling, and spirometry specialists. FindingsAt least 20 pharmacy students have participated in this project in addition to nursing (n = 8), public health (n = 1), and social work (n = 1) students. In initial and subsequent screenings, students worked alongside a team of biomedical researchers and faculty from the University of Montana. An initial cohort of 95 patients was recruited. SummaryThis unique experiential training opportunity has afforded pharmacy students access to rural community patient interaction and exposure to and performance of a variety of tests in response to an environmental health emergency. Furthermore, it enabled health professionals and researchers to assess individual and overall community health following an extreme wildfire smoke event, providing the groundwork for utilization of pharmacy students in healthcare responses to public health emergencies.

Microseismic Evidence for Bookshelf Faulting in Western Montana

AbstractOne of the most seismically active regions in the United States, located hundreds of kilometers inland from the nearest plate boundary, is the Intermountain Seismic Belt (ISB). The 6 July 2017 M 5.8 earthquake occurred 11 km southeast of Lincoln, Montana, within the ISB. This was the largest earthquake to rupture in the state of Montana since the 1959 M 7.3 Hebgen Lake earthquake. We use continuous seismic data from the University of Montana Seismic Network, the Montana Regional Seismic Network, and the U.S. Geological Survey to investigate the Lincoln aftershock sequence and to evaluate crustal stress conditions. We manually picked P- and S-wave arrival times, computed 4110 hypocenter locations and 2336 double-difference relocations, and generated focal mechanisms for 414 aftershocks (12+ polarities) in the 2 yr following the mainshock. Based on the alignment of aftershocks, we infer that the mainshock occurred on a north–northeast-trending left-lateral strike-slip fault. The orientation of the fault is unexpected, given that it strikes nearly perpendicular to the prominent Lewis and Clark line (LCL) faults in the area. Although most aftershocks concentrate near the mainshock, several distinct clusters of microseismic activity emerge along subparallel faults located primarily to the west of the mainshock. The subparallel faults also exhibit left-lateral strike-slip motion oblique to the LCL. We postulate that the aftershocks reveal the clockwise rotation of local-scale crustal blocks about vertical axes within a larger, right-lateral shear zone. The inferred block rotations are consistent with a bookshelf-faulting mechanism, which likely accommodates differential crustal motion to the north and south of the LCL region. The tension axes of well-constrained focal mechanisms indicate local northeast–southwest extension with a mean direction of N60°E.

Single-locus species delimitation and ecological niche modelling provide insights into the evolution, historical distribution and taxonomy of the Pacific chorus frogs

AbstractThe Pacific chorus frogs are a complex of three wide-ranging species (i.e. Hyliola hypochondriaca, Hyliola regilla, Hyliola sierra) whose current taxonomy remains unresolved. We conducted species delimitation analyses of these taxa using fragments of the cytochrome b and 12S–16S mtDNA genes to assess the species diversity. Importantly, we included samples from new locations throughout the range to better understand species distributions and identify potential contact zones among clades. Our analyses revealed three slightly parapatric but distinct species-level clades. Molecular dating revealed that these species began diverging in the Pleistocene c. 1.4 Mya with H. hypochondriaca and H. sierra diverging more recently c. 0.8 Mya. We found that populations from western Montana and Idaho originated recently from populations to the southwest that belong to H. sierra, rather than from H. regilla populations directly to the west. Population sizes of each species expanded c. 130–80 Kya with H. hypochondriaca exhibiting a more pronounced expansion beginning c. 100 Kya than the more gradual expansion of the other two species. The climatic niche models suggest that distributions of the three species were similar during the last interglacial (LIG) as they are today. During the Last Glacial Maximum (LGM), H. hypochondriaca and H. sierra occupied a larger range than they do today whereas H. regilla occupied a smaller refugium, shifted south from the current distribution. This study highlights the continued effectiveness of utilizing single-locus data sets for species delimitation and biogeographic analyses.

Widespread Lead Exposure in Golden Eagles Captured in Montana

ABSTRACTLead poisoning threatens many species of raptors, including golden eagles (Aquila chrysaetos). Much of this lead likely comes from bullet fragments that remain in the carcasses of animals killed by hunters. The likelihood of lead exposure may peak during fall hunting seasons and early winter until carcasses from hunting become scarce. From 2011 to 2018 in western Montana, USA, we captured 91 golden eagles in winter, tested their blood lead levels (BLL), and outfitted a subset of birds (n = 29) with global positioning system [GPS] transmitters. Nearly all golden eagles (94.5%) had elevated BLL (≥10 μg/dL), and 8 of them had BLL above clinical exposure (>60 μg/dL), where they may lose coordination and experience a host of other neurological and physiological disorders. Golden eagles caught late in winter tended to have lower BLL than those caught earlier. At least 69% of the golden eagles equipped with GPS transmitters migrated northward, spending the summer throughout Alaska, USA, and northwestern Canada. Blood lead levels did not differ between migratory and nonmigratory golden eagles. Overall, elevated BLL are widespread among golden eagles throughout winter in western Montana. Promoting nonlead hunting ammunition in areas with high densities of golden eagles will reduce the birds' lead exposure. © 2020 The Wildlife Society.

Trace-element and Sr and Nd isotopic geochemistry of Cretaceous bentonites in Wyoming and South Dakota tracks magmatic processes during eastward migration of Farallon arc plutons

Abstract Bentonite beds, which are clay deposits produced by the submarine alteration of volcanic tephra, preserve millions of years of volcanic products linked to magmatic systems for which records are otherwise lost through erosion and alteration. Cretaceous strata from the Bighorn Basin, Wyoming, and southwestern South Dakota contain bentonites that originated from arc magmatism produced by subduction of the Farallon plate. We analyzed the bulk major- and trace-element geochemistry, and the 87Sr/86Sr (n = 87) and 143Nd/144Nd (n = 26) isotopic compositions of individual bentonite beds from these areas spanning 40 m.y. of volcanism to recover signals of magmatic processes and to attempt to trace bentonite geochemical and isotopic signatures to contemporaneous Cordilleran plutonic rocks. Using multiple immobile elements (e.g., Zr, TiO2, Nb, Ta, and rare earth elements), distinct temporal trends show variations in the effects of mineral fractionation and changes in crustal thickness. Bentonite Sr and Nd isotopic compositions allow ash beds to be correlated with specific batholithic complexes in Idaho and western Montana. With this data set, we observed the following: (1) The volcanic arc migrated across the 0.706 isopleth between 115 and 105 Ma; (2) between 105 and 95 Ma, magmatism stalled in central Idaho and was supported through significant MASH (mixing-assimilation-storage-homogenization) processing; (3) by 85 Ma, a shallowing subduction angle resulted in the eastward migration of the volcanic front into western Montana while volcanism in Idaho diminished; and (4) around 75 Ma, evidence of Idaho volcanism is lost. Montana plutonism continued with significant assimilation of radiogenic basement and regional centers of local magma emplacement (i.e., Pioneer batholith).

Determining the source of placer gold in the Anaconda metamorphic core complex supradetachment basin using detrital zircon U-Pb geochronology, western Montana, USA

Abstract Despite the widespread occurrence and economic significance of gold placer deposits, modern provenance studies of placer sediments remain largely qualitative. This study applies detrital zircon (DZ) geochronology to determine the source of zircon in placer deposits. We then evaluate the provenance of the zircon to assess whether the gold might have been derived from the same sources, thereby providing a case study of the use of DZ geochronology applied to placers. We present a new set of DZ U-Pb ages (n = 1058) and Lu-Hf (n = 61) isotopic data from four placer deposit samples collected from the Pioneer District of western Montana (USA). Each of the four samples yielded similar age spectra, with a range of U-Pb ages between 3000 and 25 Ma. We interpret that ≥250 Ma zircons were recycled from the Mesoproterozoic Belt Supergroup, Paleozoic–Mesozoic sedimentary rocks, and the Upper Cretaceous–Paleocene Beaverhead Group. Our 237 DZ U-Pb ages ≤250 Ma reveal two prominent age-probability peaks centered at ca. 69 Ma and ca. 26 Ma, which we interpret to record first-cycle derivation from the Royal stock and nearby Dillon Volcanics, respectively. We evaluate these data using an inverse Monte Carlo DZ unmixing model that calculates relative contributions from plausible source units, determining a 12% contribution from the Royal stock and a 43% contribution from the Beaverhead Group. A current absence of the Beaverhead Group in the hypothesized source region suggests complete erosion of the unit into the placer-bearing basin. Detrital zircon geochronology, Hf isotopic data, and the unmixing modeling results offer the first zircon-based support for previous interpretations that the Late Cretaceous Royal stock precipitated gold along its contact with overlying Proterozoic–Mesozoic sedimentary strata. Subsequent exhumation and erosion of the lode source led to gold deposition in the Anaconda metamorphic core complex supradetachment basin during the late Oligocene–late Miocene. The worldwide occurrence of gold placer deposits with unknown source areas provides abundant opportunity to apply these techniques elsewhere.

Effect of terrain, environment and infrastructure on potential CO2 pipeline corridors: a case study from North-Central USA

In this paper, a study has been undertaken with the objective to delineate the potential CO2 pipeline corridors through the north-central region of the USA including North Dakota, Montana, Wyoming, Utah and Colorado to enable the implementation of carbon capture storage and utilization projects. A combination of GIS along with analytical hierarchy process is used to identify regions with high potential for pipeline development. A total of 19 thematic information layers have been utilized to map the study area which reflect the effect of ecology, environment and existent infrastructure towards laying new CO2 pipeline networks. Weights are assigned to each class in the thematic maps based on their characteristics, capacity of building and maintaining pipelines and the potential environmental risk of CO2 pipelines from the literature. A tract suitability index is developed which identified Western North Dakota, central Wyoming and Western Montana as regions highly suitable for development of CO2 pipelines. The study reveals that 54.5% of the region of interest is suitable for pipeline construction while 15.39% of land is classified as either poor or infeasible candidate for locating CO2 pipelines. It is also revealed that pipeline right-of-way, water bodies and population centres have the most significant impact on future pipeline development. The results of the study are also varied to show the effect different factors could have on the potential CO2 corridors and are then correlated with existent pipeline in study area to check the validity of the resulting analysis.

Native and non-native understory vegetation responses to restoration treatments in a dry conifer forest over 23 years

Restoration of fire-prone forests is a common practice intended to increase resilience to wildfire, drought, and bark beetles. However, the long-term effects of restoration treatments on understory species, particularly non-native species, are poorly understood. We investigated long-term (23 years) effects of restoration treatments on native and non-native understory vegetation at the Lick Creek Demonstration-Research Forest, a ponderosa pine/Douglas-fir forest in western Montana, USA. Restoration treatments included shelterwood-with-reserves harvesting implemented alone or in combination with prescribed burning under either wet or dry conditions, plus a no-action control. Changes in understory vegetation cover and species richness (measured pre-treatment, 2, 3, 4, 5, 15, and 23 years post-treatment) were analyzed by growth form (forb, graminoid, and shrub) and origin (native vs. non-native) to test how understory dynamics respond to disturbance intensity. We also accounted for basal area of trees to consider how differential overstory conditions affected final treatment responses. Native vegetation cover responses followed predictable successional patterns: after initial declines by grasses and shrubs immediately after disturbance (43% and 40% reduction, respectively), all groups increased, peaked about five years after disturbance (74% greater than pre-treatment on average) and then generally declined to levels similar to pre-treatment. The magnitude of these changes was related to the disturbance intensity (i.e., cut-and-dry-burn > cut-and-wet-burn > cut-and-no-burn > control). For all native groups, variation in tree basal area at year 23 explained remaining differences among treatments. Non-native grasses and forbs followed similar overall trajectories, but cover of these groups increased immediately after treatment and rose more steeply to their peak five years post-treatment (12 times greater than pre-treatment on average), with differences persisting through year 23. While non-native cover responses in year 23 were negatively correlated with tree basal area, as seen for native taxa, overstory conditions alone did not fully account for differences between burned and unburned treatments. This suggests that effects of burning unrelated to tree regeneration favored non-natives. Nonetheless, in the last sampling year treatment-induced increases in non-native cover were modest, suggesting that after 23 years understory responses in this forest type were fairly resilient to restoration treatments. However, because forest structure and composition were reverting to pre-treatment conditions by year 23, retreatment may be necessary at intervals of <23-years in our system, which could exacerbate non-native responses. Results of our long-term experiment provide insight into how different forest restoration treatments and their timing may affect understory susceptibility to non-native plant invasion.

Groundwater Upwelling Regulates Thermal Hydrodynamics and Salmonid Movements during High‐Temperature Events at a Montane Tributary Confluence

AbstractThe Smith River is a popular recreational sport fishery in western Montana, but salmonid abundances there are thought to be artificially limited by riparian land‐use alterations, irrigation water withdrawals, and high summer water temperatures. We used integrated networks of temperature loggers, PIT tag antenna stations, and in situ temperature mapping to investigate the thermal hydrodynamics and associated movements of PIT‐tagged salmonids at the confluence of Tenderfoot Creek, a major, unaltered coldwater tributary of the Smith River. Contrary to expectations, Tenderfoot Creek itself was not used as a thermal refuge by salmonids during periods of high water temperatures in Smith River; rather, its cool outflow plume into the main stem was used instead. Mean daily outflow water temperatures averaged 2.9°C lower than those of the Smith River during summer and ranged from 0.5°C to 6.1°C lower. Moreover, measured and estimated temperatures in the outflow were cooler (by up to 2.8°C) than in Tenderfoot Creek itself at times as a result of groundwater upwelling at the confluence. Detections of PIT‐tagged fish in the thermal plume increased, especially at night, when daily mean water temperatures exceeded 20°C in the main‐stem Smith River; more than four times as many PIT‐tagged fish were detected in the plume (N = 52) than along the opposite bank (N = 12), which ostensibly afforded better cover. Coldwater tributary confluences may provide superior thermal refuges for salmonids—cooler than the tributaries themselves—when water temperatures in river main stems are stressful.

Forest restoration treatments in a ponderosa pine forest enhance physiological activity and growth under climatic stress.

As the climate warms, drought will increasingly occur under elevated temperatures, placing forest ecosystems at growing risk of extensive dieback and mortality. In some cases, increases in tree density following early 20th-century fire suppression may exacerbate this risk. Treatments designed to restore historical stand structure and enhance resistance to high-severity fire might also alleviate drought stress by reducing competition, but the duration of these effects and the underlying mechanisms remain poorly understood. To elucidate these mechanisms, we evaluate tree growth, mortality, and tree-ring stable-carbon isotope responses to stand-density reduction treatments with and without prescribed fire in a ponderosa pine forest of western Montana. Moderate and heavier cutting experiments (basal area reductions of 35% and 56%, respectively) were initiated in 1992, followed by prescribed burning in a subset of the thinned units. All treatments led to a growth release that persisted to the time of resampling. The treatments had little effect on climate-growth relationships, but they markedly altered seasonal carbon isotope signals and their relationship to climate. In burned and unburned treatments, carbon isotope discrimination (Δ13 C) increased in the earlywood (EW) and decreased in the latewood (LW) relative to the control. The sensitivity of LW Δ13 C to late-summer climate also increased in all treatments, but not in the control. Such increased sensitivity indicates that the reduction in competition enabled trees to continue to fix carbon for new stem growth, even when the climate became sufficiently stressful to stop new assimilation in slower-growing trees in untreated units. These findings would have been masked had we not separated EW and LW. The importance of faster growth and enhanced carbon assimilation under late-summer climatic stress became evident in the second decade post-treatment, when mountain pine beetle activity increased locally, and tree mortality rates in the controls of both experiments increased to more than twice those in their respective treatments. These findings highlight that, when thinning is used to restore historical forest structure or increase resistance to high-severity fire, there will likely be additional benefits of enhanced growth and physiological activity under climatic stress, and the effects may persist for more than two decades.

Biochar alters nitrogen and phosphorus dynamics in a western rangeland ecosystem

Application of biochar to soils has been proposed as a novel approach to managing wood residuals, enhancing soil carbon (C) storage and improving soil fertility; however, the majority of biochar studies have been conducted in agricultural systems that rely on tillage and nutrient inputs associated with annual cropping schemes. Few studies have evaluated the influence of biochar on soil processes in semi-natural rangeland ecosystems that feature more complex plant communities, lack deep soil disturbance, and have relatively few external nutrient inputs. In August 2018, biochar produced using wood waste from a lumber mill in Columbia Falls, MT, USA was applied to surface soils in replicated plots at an experimental ranch in western Montana to test the impact of biochar on soil C storage and nutrient management. A series of soil biochemical properties including total soil C and nitrogen (N), microbial N functional genes, available phosphorus (P) and the net accumulation of nutrients below surface soil layer were evaluated over a one-year period following biochar addition with or without a poultry litter based organic fertilizer. Biochar used alone slightly reduced soil NH4+, significantly increased soil nitrification potential, increased the relative abundance of the bacterial amoA gene, and increased the soil nitrate (NO3−) pool size, while having no net effect on soil inorganic N accumulation below surface soil. By contrast, biochar charged with poultry litter (termed “charged biochar”) had no effect on NH4+ availability, but had a positive effect on amoA abundance. Charged biochar significantly reduced NH4+ accumulation below 25 cm depth compared to poultry litter alone. Biochar additions led to a shift towards a more fungal dominated community and a general increase in P availability. However, biochar alone also contributed to a greater amount of soluble P collected below surface soil, an effect slightly attenuated when biochar was applied with poultry litter. Soil pH increased from 5.7 to 6.9 in response to biochar addition and was one of the dominant factors governing the observed changes in soil processes. Charged biochar helped retain soil nutrients and promoted soil C storage in this semi-natural rangeland system over one growing season. Changes in these soil pools and fluxes may influence various trophic groups affecting ecosystem functioning over time.

Restoration Treatments Improve Overstory Tree Resistance Attributes and Growth in a Ponderosa Pine/Douglas-Fir Forest

Research Highlights: This study provides much needed insight into the development of resistance to disturbance and growth dynamics of overstory trees in response to restoration-based fuel reduction, and will be useful to scientists and managers attempting to better grasp the relative merits of restoration treatment types. Background and Objectives: Restoration-based fuel reduction treatments are common in dry, fire-prone forests of the western United States. The primary objective of such treatments is to immediately reduce a stand’s crown fire hazard. However, the impact of these treatments on residual trees is relevant to assess their longevity and resistance to future disturbances. In this study, we evaluate the effects of restoration on retained overstory ponderosa pine (Pinus ponderosa Lawson &amp; C. Lawson) and Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) trees in western Montana, where treatments were experimentally implemented 13 years prior as part of the national Fire and Fire Surrogate study. Materials and Methods: We examined tree attributes in response to the following replicated treatments: thin-only, burn-only, thin + burn, and a no-action control. We analyzed three different tree attributes that confer resistance to common disturbances: height-to-diameter ratio (resistance to wind), bark thickness (resistance to surface fire), and growth efficiency (resistance to bark beetles). Results: Our models suggest that thinning (with or without burning) alters tree attributes relative to the control in a manner that may increase tree resistance to wind and snow breakage, surface fire, and biotic agents such as bark beetles. Further analysis of annual growth of ponderosa pine and Douglas-fir varied by treatment type: thinning-based restoration (thin-only and thin + burn) increased diameter growth for both species, crown length and width in ponderosa pine, and crown length in Douglas-fir relative to unthinned treatments. Burning (burn-only and thin + burn) did not significantly affect tree growth relative to unburned treatments. Conclusions: While low-severity prescribed burning treatments are often used for restoration and have various ecosystem benefits, this study demonstrates that thinning (alone or in addition to burning) produces more measureable, beneficial results to overstory tree disturbance resistance metrics and growth.

Reconstructing summer upper-level flow in the northern Rocky Mountains using an alpine larch tree-ring chronology

Mid-latitude mesoscale weather during the climatological summer is strongly influenced by fluctuations in synoptic-scale circulation patterns. Previous research has linked Arctic amplification to alterations in summer synoptic climatology, leading to more extreme weather events in the mid-latitudes. In this study, seasonal (JJA) upper-level (500 hPa) atmospheric flow is reconstructed in the mid-latitudes using an alpine larch Larix lyallii Parl. tree-ring chronology sampled from western Montana. Significant relationships were found between alpine larch radial growth and upper-level flow patterns derived from the North American Regional Reanalysis dataset (1979-2015). Meridional and zonal flows that manifest in ridging are associated with enhanced radial growth of alpine larch (i.e. meridional flow west [r = 0.504, p = 0.001] and zonal flow north [r = 0.642, p &lt; 0.001] of the study site). Meridional and zonal flows associated with troughing result in decreased radial growth (i.e. meridional flow east [r = -0.497, p = 0.001] and zonal flow south [r = -0.584, p &lt; 0.001] of the study site). Using the leave-one-out method, a linear regression model was calibrated and verified between a principal component analysis score derived from measurements of upper-level flow in western North America and alpine larch tree growth. The 444 yr climate reconstruction of summer 500 hPa flow suggests that ridging is becoming more intense over the western United States and Canada since the 1980s.

Pyrogenic Carbon Generation From Fire and Forest Restoration Treatments

Pyrogenic carbon (PyC) is a chemically stable form of carbon (C) generated during fire events and is one of the few legacies of fire recorded in soil; however, the significance of this material as a form of C storage in forest ecosystems has received only limited scientific attention, and currently relatively little is known regarding the quantity of PyC generated during forest restoration efforts that include prescribed (Rx) fire. Quantifying the rate of PyC production during fire events is essential to estimating the potential for ecosystem C storage in forest systems where wildfire is a common natural disturbance and Rx fire is used as a fuel management practice. In this paper, we report on the average mass of PyC formation during wildfire and Rx fire events in forest ecosystems through synthesizing data across 12 published studies; and we report empirical data on the mass of PyC generated in a replicated Rx fire study in a ponderosa pine (Pinus ponderosa) forest ecosystem of the Rocky Mountain West. Our data synthesis showed that PyC (consisting of PyC from downed wood, O horizon, and mineral soil) was, on average, produced at a rate of 5.2% of biomass exposed to fire in forest ecosystems, and a single wildfire or Rx fire event generated about 2.2 Mg C ha-1 in the form of PyC. In our empirical study, we collected O horizon and surface mineral soil (0 – 10 cm) samples from the Fire and Fire Surrogates study plots in western Montana (treatments: control, thin, burn, thin&burn) that was burned in 2002 and analyzed for PyC content using a wet digestion method. The two Rx fire treatments had significantly higher PyC contents in the O horizon compared to the control. In the thin&burn treatment, PyC was formed at a rate of about 12% of total biomass consumed yielding approximately 2.24 Mg PyC C ha-1 in the O horizon. Mineral soil PyC contents were not affected by Rx fire. Wildfire and Rx fires generate a substantial amount of PyC which initially resides in the O horizon and serves as a stable form of ecosystem C.