Northwestern Montana Research Articles

Evolution of eggshell structure during rapid range expansion in a passerine bird

Summary1. Environmental factors such as temperature, humidity and partial oxygen pressure can affect avian eggshell structure because gas exchange across the shell must allow sufficient water loss while preventing dehydration of the embryo. Studies of species with known chronology of colonization of novel environments provide a powerful insight into the relative importance of ecological factors shaping the evolution of eggshell structure.2. Here, we examined changes in eggshell structure that accompanied rapid range expansion of house finches (Carpodacus mexicanus) across North America. We analysed thickness and pore density in eggshells from three ecologically distinct populations: the native desert population in southwestern Arizona and two 30‐year‐old populations in the northwestern (north‐west Montana) and southeastern (south‐east Alabama) parts of the species’ range. We also conducted cross‐foster exchanges of freshly laid eggs within and between the northwestern and southeastern populations to examine consequences of population differences in eggshell structure on embryo development.3. Eggshell structure was most distinct in a recently established population inhabiting higher humidity environment (southeastern Alabama), where eggs were the largest, eggshells the thickest and pore density the lowest. Populations that experienced highly distinct ambient temperatures (southwestern Arizona and northwestern Montana) nevertheless had similar eggshell structure. These results were corroborated by experiments where humidity differences between cross‐fostered nests had twice the effect on embryo survival compared to the effect of change in ambient temperature. Correspondingly, experimental egg exchanges between southeastern. Alabama and northwestern Montana populations were associated with fourfold increase in embryo mortality compared to within‐population egg exchanges.4. We document rapid evolution of eggshell structure in response to colonization of novel environments and establish the relative importance of environmental factors on avian eggshells. We discuss these results in relation to population variation in incubation behaviour and its ability to shield eggshell structure from the selection exerted by novel environments.

Linking landscape characteristics to local grizzly bear abundance using multiple detection methods in a hierarchical model

Few studies link habitat to grizzly bear Ursus arctos abundance and these have not accounted for the variation in detection or spatial autocorrelation. We collected and genotyped bear hair in and around Glacier National Park in northwestern Montana during the summer of 2000. We developed a hierarchical Markov chain Monte Carlo model that extends the existing occupancy and count models by accounting for (1) spatially explicit variables that we hypothesized might influence abundance; (2) separate sub-models of detection probability for two distinct sampling methods (hair traps and rub trees) targeting different segments of the population; (3) covariates to explain variation in each sub-model of detection; (4) a conditional autoregressive term to account for spatial autocorrelation; (5) weights to identify most important variables. Road density and per cent mesic habitat best explained variation in female grizzly bear abundance; spatial autocorrelation was not supported. More female bears were predicted in places with lower road density and with more mesic habitat. Detection rates of females increased with rub tree sampling effort. Road density best explained variation in male grizzly bear abundance and spatial autocorrelation was supported. More male bears were predicted in areas of low road density. Detection rates of males increased with rub tree and hair trap sampling effort and decreased over the sampling period. We provide a new method to (1) incorporate multiple detection methods into hierarchical models of abundance; (2) determine whether spatial autocorrelation should be included in final models. Our results suggest that the influence of landscape variables is consistent between habitat selection and abundance in this system.

Observational Studies of Atmospheric Aerosols over Bozeman, Montana, Using a Two-Color Lidar, a Water Vapor DIAL, a Solar Radiometer, and a Ground-Based Nephelometer over a 24-h Period

Abstract Coordinated observational data of atmospheric aerosols were collected over a 24-h period between 2300 mountain daylight time (MDT) on 27 August 2009 and 2300 MDT on 28 August 2009 at Bozeman, Montana (45.66°N, 111.04°W, elevation 1530 m) using a collocated two-color lidar, a diode-laser-based water vapor differential absorption lidar (DIAL), a solar radiometer, and a ground-based nephelometer. The optical properties and spatial distribution of the atmospheric aerosols were inferred from the observational data collected using the collocated instruments as part of a closure experiment under dry conditions with a relative humidity below 60%. The aerosol lidar ratio and aerosol optical depth retrieved at 532 and 1064 nm using the two-color lidar and solar radiometer agreed with one another to within their individual uncertainties while the scattering component of the aerosol extinction measured using the nephelometer matched the scattering component of the aerosol extinction retrieved using the 532-nm channel of the two-color lidar and the single-scatter albedo retrieved using the solar radiometer. Using existing aerosol models developed with Aerosol Robotic Network (AERONET) data, a thin aerosol layer observed over Bozeman was most likely identified as smoke from forest fires burning in California; Washington; British Columbia, Canada; and northwestern Montana. The intrusion of the thin aerosol layer caused a change in the atmospheric radiative forcing by a factor of 1.8 ± 0.5 due to the aerosol direct effect.

Evolution of the Cordilleran foreland basin system in northwestern Montana, U.S.A.

New lithostratigraphic and chronostratigraphic, geochronologic, and sedimentary petrologic data illuminate the history of development of the North American Cordilleran foreland basin system and adjacent thrust belt from Middle Jurassic through Eocene time in northwestern Montana. The oldest deposits in the foreland basin system consist of relatively thin, regionally tabular deposits of the marine Ellis Group and fluvial-estuarine Morrison Formation, which accumulated during Bajocian to Kimmeridgian time. U-Pb ages of detrital zircons and sandstone modal petrographic data indicate that by ca. 170 Ma, miogeoclinal strata were being deformed and eroded in hinterland regions. Sandstones of the Swift and Morrison Formations contain detrital zircons derived from the Intermontane belt. The Jurassic deposits probably accumulated in the distal, back-bulge depozone of an early foreland basin, as suggested by the slow rates of tectonic subsidence and tabular geometry. A regional unconformity separates the Jurassic strata from late Barremian(?) foredeep deposits. This unconformity possibly resulted as a combined effect of forebulge migration, decreased dynamic subsidence, and eustatic sea-level fall. The late Barremian(?)-early Albian Kootenai Formation is the first unit that consistently thickens westward, as would be expected in a foredeep depozone. The subsidence curve at this time begins to show the convex-upward pattern characteristic of foredeeps. By Albian time, the fold-and-thrust belt had propagated to the east and incorporated Proterozoic rocks of the Belt Supergroup, as indicated by sandstone compositions, detrital zircon ages in the Blackleaf Formation, and by crosscutting relationships in thrust sheets involving Belt Supergroup rocks in the thrust belt. A major episode of marine inundation and black shale deposition (Marias River Shale) occurred between the Cenomanian and mid-Santonian, and was followed by a regressive succession represented by the Upper Santonian–mid-Campanian Telegraph Creek, Virgelle, and Two Medicine Formations. Provenance data do not resolve the timing of individual thrust displacements during Cenomanian–early Campanian time. The Upper Campanian Bearpaw Formation represents the last major marine inundation in the foreland basin. By latest Campanian time, a major episode of slip on the Lewis thrust system had commenced, as recorded in the foreland by the Willow Creek and St. Mary River Formations in the proximal foredeep depozone. The final stage in the evolution of the Cordilleran fold-and-thrust belt and foreland basin system is recorded by the Paleocene–early Eocene Fort Union and Wasatch Formations, which were preserved in the distal foreland region. Regional extensional faulting along the fold-and-thrust belt began during the middle Eocene. The results presented here enable the establishment of links between previous geological work in Canada and the better known parts of the Cordilleran foreland basin in the United States.

Visualizing Spatial Heterogeneity of Western U.S. Climate Variability

Abstract Monthly climatologies (1971–2000 monthly averages) for stations in the western United States, obtained from the NOAA/National Climatic Data Center (NCDC), are used to illustrate the spatial variations in the annual cycle of climate. Animated map sequences of temperature and precipitation, their average, intermonthly changes, and the local timing of annual maxima or minima provide a comprehensive spatiotemporal baseline of regional climate. The animated maps illustrate three scales of variation: 1) broadscale patterns related to the annual cycle of insolation and hemispheric-scale atmospheric circulation features; 2) mesoscale patterns related to location on the continent and the influence of specific regional circulation features like those associated with the North American monsoon; and 3) smaller-scale spatial variations, related to the mediation by local physiography of the influence of large-scale circulation. Although most western U.S. stations have temperature maxima in July, a delay occurs at stations along the West Coast and interior Washington, northern Idaho, and Montana. A seesaw pattern of precipitation maxima is evident between coastal areas (winter dominated) and the interior (summer dominated). Cluster analyses of the ratio of monthly-to-annual precipitation values for each station identify regions with similar annual cycles of precipitation. Regions of high spatial heterogeneity in the timing of when precipitation occurs include the northern Rocky Mountains, Utah, Arizona, and northwestern Montana. The superimposition of these three scales of spatial variability leads to steep gradients and, in some regions, considerable spatial heterogeneity in annual precipitation. The regional patterns of precipitation heterogeneity highlight vulnerability to drought, especially in regions of the interior west that do not have a dominant precipitation month or season.

Restoration of genetic connectivity among Northern Rockies wolf populations

Probably no conservation genetics issue is currently more controversial than the question of whether grey wolves (Canis lupus) in the Northern Rockies have recovered to genetically effective levels. Following the dispersal-based recolonization of Northwestern Montana from Canada, and reintroductions to Yellowstone and Central Idaho, wolves have vastly exceeded population recovery goals of 300 wolves distributed in at least 10 breeding pairs in each of Wyoming, Idaho and Montana. With >1700 wolves currently, efforts to delist wolves from endangered status have become mired in legal battles over the distinct population segment (DPS) clause of the Endangered Species Act (ESA), and whether subpopulations within the DPS were genetically isolated. An earlier study by vonHoldt etal. (2008) suggested Yellowstone National Park wolves were indeed isolated and was used against delisting in 2008. Since then, wolves were temporarily delisted, and a first controversial hunting season occurred in fall of 2009. Yet, concerns over the genetic recovery of wolves in the Northern Rockies remain, and upcoming District court rulings in the summer of 2010 will probably include consideration of gene flow between subpopulations. In this issue of Molecular Ecology, vonHoldt etal. (2010) conduct the largest analysis of gene flow and population structure of the Northern Rockies wolves to date. Using an impressive sampling design and novel analytic methods, vonHoldt etal. (2010) show substantial levels of gene flow between three identified subpopulations of wolves within the Northern Rockies, clarifying previous analyses and convincingly showing genetic recovery.

Death from anthropogenic causes is partially compensatory in recovering wolf populations

There is substantial interest in how mortality rates affect animal populations, but mechanisms explaining when and under what circumstances particular causes of death incur demographic responses are far from clear. In theory, small or expanding populations should experience additive mortality from anthropogenic causes of death, but whether such effects are homogenous across a population or expressed only in certain high-risk individuals is open for debate. We used competing risks models to analyze mortality patterns among radio-collared wolves ( Canis lupus, n = 711) from three populations in northwestern United States (1982–2004), and evaluated the degree to which anthropogenic mortality was additive vs. compensatory to natural demographic processes. Almost 80% ( n = 320) of wolves dying of known fates were killed by anthropogenic causes (legal control, illegal killing, harvest in Canada, vehicle collision), and additive effects of anthropogenic mortality were most pronounced in northwestern Montana where wolf exposure to humans and livestock was high compared to either the Greater Yellowstone Area or central Idaho, where anthropogenic risk was lower. In contrast, risk from natural hazards was lower in northwestern Montana than in the other areas, implying some degree of compensatory mortality from anthropogenic risk. Animals recruited to the study following human–wolf conflict had markedly higher anthropogenic risk than those recruited for standard monitoring purposes, and juvenile wolves as well as dispersers, succumbed to higher anthropogenic risk. Multivariate models revealed that increasing wolf population density promoted higher anthropogenic risk and reduced natural risk, indicating that partially-compensatory effects of anthropogenic mortality actually became increasingly additive with population density. The observed compensatory mortality and hazard heterogeneity in our study implies that demographic responses to mortality risk may be complex and more subtle than previously thought; the density-dependent effect of anthropogenic mortality portends a stabilizing influence of humans on recovering wolf populations. We conclude that future assessment of the role of anthropogenic mortality should include individual-based hazard estimation as a complement to traditional population-level approaches.

Mammals from the St. Mary River Formation (Upper Cretaceous), Montana

ABSTRACT The St. Mary River Formation (Fm) crops out in restricted parts of southwestern Alberta and northwestern Montana. Mammals are poorly known from the formation, but material collected from Alberta has played a role in recognizing an ‘Edmontonian’ Land Mammal Age between the better sampled Judithian (∼79–74 Ma) and Lancian (∼67–65 Ma) land mammal ages. New and well-preserved multituberculate and metatherian material collected from a single locality in the lower third of the St. Mary River Fm of Montana adds complexity to the interpretation of a discrete ‘Edmontonian’ Land Mammal Age. We report here three taxa of multituberculates (Paracimexomys propriscus, sp. nov., Nidimys occultus, gen. et sp. nov., and a primitive ptilodontoid of comparable size to the smallest species of MesodmaJepsen, 1940, and CimexomysSloan and Van Valen, 1965) and two metatherian taxa (Leptalestes toevsi, sp. nov., and Turgidodon russelli). Paracimexomys propriscus and Turgidodon russelli are probably conspecific with Judithian forms that did not survive to the Lancian. Leptalestes toevsi and Nidimys occultus are taxa unique to the ‘Edmontonian’ but may represent lineages that separated from their closest relatives prior to the Judithian. Lack of Lancian-aspect mammals contrasts sharply with other ‘Edmontonian’ faunas, including an older fauna recovered from the Williams Fork Fm of Colorado, suggesting that the replacement of Judithian mammals by Lancian forms was a complex transition. We tentatively suggest that the Lancian fauna may have expanded into the northern part of the Western Interior in response to the appearance of new terrestrial habitats as sea level fell.

Evaluation of Bear Rub Surveys to Monitor Grizzly Bear Population Trends

Abstract: Wildlife managers need reliable estimates of population size, trend, and distribution to make informed decisions about how to recover at‐risk populations, yet obtaining these estimates is costly and often imprecise. The grizzly bear (Ursus arctos) population in northwestern Montana, USA, has been managed for recovery since being listed under the United States Endangered Species Act in 1975, yet no rigorous data were available to evaluate the program's success. We used encounter data from 379 grizzly bears identified through bear rub surveys to parameterize a series of Pradel model simulations in Program MARK to assess the ability of noninvasive genetic sampling to estimate population growth rates. We evaluated model performance in terms of 1) power to detect gender‐specific and population‐wide declines in population abundance, 2) precision and relative bias of growth rate estimates, and 3) sampling effort required to achieve 80% power to detect a decline within 10 years. Simulations indicated that ecosystem‐wide, annual bear rub surveys would exceed 80% power to detect a 3% annual decline within 6 years. Robust‐design models with 2 simulated surveys per year provided precise and unbiased annual estimates of trend, abundance, and apparent survival. Designs incorporating one survey per year require less sampling effort but only yield trend and apparent survival estimates. Our results suggest that systematic, annual bear rub surveys may provide a viable complement or alternative to telemetry‐based methods for monitoring trends in grizzly bear populations.

Spatial analysis of feline immunodeficiency virus infection in cougars

The cougar ( Puma concolor) is a large predatory feline found widely in the Americas that is susceptible to feline immunodeficiency virus (FIV), a fast-evolving lentivirus found in wild feline species that is analogous to simian immunodeficiency viruses in wild primates and belongs to the same family of viruses as human immunodeficiency virus. FIV infection in cougars can lead to a weakened immune system that creates opportunities for other infecting agents. FIV prevalence and lineages have been studied previously in several areas in the western United States, but typically without spatially explicit statistical techniques. To describe the distribution of FIV in a sample of cougars located in the northern Rocky Mountain region of North America, we first used kernel density ratio estimation to map the log relative risk of FIV. The risk surface showed a significant cluster of FIV in northwestern Montana. We also used Bayesian cluster models for genetic data to investigate the spatial structure of the feline immunodeficiency virus with virus genetic sequence data. A result of the models was two spatially distinct FIV lineages that aligned considerably with an interstate highway in Montana. Our results suggest that the use of spatial information and models adds novel insight when investigating an infectious animal disease. The results also suggest that the influence of landscape features likely plays an important role in the spatiotemporal spread of an infectious disease within wildlife populations.

Carbon storage in old-growth and second growth fire-dependent western larch (Larix occidentalis Nutt.) forests of the Inland Northwest, USA

Abstract There is limited understanding of the carbon (C) storage capacity and overall ecological structure of old-growth forests of western Montana, leaving little ability to evaluate the role of old-growth forests in regional C cycles and ecosystem level C storage capacity. To investigate the difference in C storage between equivalent stands of contrasting age classes and management histories, we surveyed paired old-growth and second growth western larch (Larix occidentalis Nutt)–Douglas-fir (Pseudostuga menziesii var. glauca) stands in northwestern Montana. The specific objectives of this study were to: (1) estimate ecosystem C of old-growth and second growth western larch stands; (2) compare C storage of paired old-growth–second growth stands; and (3) assess differences in ecosystem function and structure between the two age classes, specifically measuring C associated with mineral soil, forest floor, coarse woody debris (CWD), understory, and overstory, as well as overall structure of vegetation. Stands were surveyed using a modified USFS FIA protocol, focusing on ecological components related to soil, forest floor, and overstory C. All downed wood, forest floor, and soil samples were then analyzed for total C and total nitrogen (N). Total ecosystem C in the old-growth forests was significantly greater than that in second growth forests, storing over 3 times the C. Average total mineral soil C was not significantly different in second growth stands compared to old-growth stands; however, total C of the forest floor was significantly greater in old-growth (23.8 Mg ha−1) compared to second growth stands (4.9 Mg ha−1). Overstory and coarse root biomass held the greatest differences in ecosystem C between the two stand types (old-growth, second growth), with nearly 7 times more C in old-growth trees than trees found on second growth stands (144.2 Mg ha−1 vs. 23.8 Mg ha−1). Total CWD on old-growth stands accounted for almost 19 times more C than CWD found in second growth stands. Soil bulk density was also significantly higher on second growth stands some 30+ years after harvest, demonstrating long-term impacts of harvest on soil. Results suggest ecological components specific to old-growth western larch forests, such as coarse root biomass, large amounts of CWD, and a thick forest floor layer are important contributors to long-term C storage within these ecosystems. This, combined with functional implications of contrasts in C distribution and dynamics, suggest that old-growth western larch/Douglas-fir forests are both functionally and structurally distinctive from their second growth counterparts.

Jurassic onset of foreland basin deposition in northwestern Montana, USA: Implications for along-strike synchroneity of Cordilleran orogenic activity

Stratigraphic, provenance, and subsidence analyses suggest that by the Middle to Late Jurassic a foreland basin system was active in northwestern Montana (United States). U-Pb ages of detrital zircons and detrital modes of sandstones indicate provenance from accreted terranes and deformed miogeoclinal rocks to the west. Subsidence commenced ca. 170 Ma and followed a sigmoidal pattern characteristic of foreland basin systems. Thin Jurassic deposits of the Ellis Group and Morrison Formation accumulated in a backbulge depozone. A regional unconformity and/or paleosol zone separates the Morrison from Early Cretaceous foredeep deposits of the Kootenai Formation. The model presented here is consistent with regional deformation events registered in hinterland regions, and challenges previous interpretations of a strongly diachronous onset of Cordilleran foreland basin deposition from northwestern Montana to southern Canada.

Evolution of the northern Cordilleran foreland basin during the middle Cretaceous

A detailed study of the middle–late Cenomanian Belle Fourche Formation in southern Alberta helps to elucidate the evolution of the Cordilleran foreland basin during the middle Cretaceous. By combining isopach data from the study area with previously published research on middle Cretaceous stratigraphy in adjacent areas of northern Alberta and Montana, it is possible to defi ne the position of the proximal foredeep, forebulge, and backbulge depozones from late Albian to Cenomanian time. The foredeep probably lay to the west of the present fold-thrust belt in western Montana from the late Albian to middle Cenomanian, with the forebulge extending northwestward through the northwest corner of Wyoming and curving northeastward through the present-day location of Calgary. Thick Dunvegan deposits (middle Cenomanian) of northwestern Alberta are interpreted as foredeep deposits. Palinspastic reconstruction of the southern Rocky Mountains of Alberta is consistent with a location of the foredeep within the area of the present fold-thrust belt. At the end of the middle Cenomanian the forebulge retreated westward within Alberta, defi ning a nearly straight northwest-southeast trend through the present position of the Foothills and Front Ranges of Alberta and northwestern Montana. The change in trend of the foreland basin at the end of the middle Cenomanian may refl ect the change of convergence vectors along the western margin of North America during the middle Cretaceous.

Demography and Genetic Structure of a Recovering Grizzly Bear Population

ABSTRACT Grizzly bears (brown bears; Ursus arctos) are imperiled in the southern extent of their range worldwide. The threatened population in northwestern Montana, USA, has been managed for recovery since 1975; yet, no rigorous data were available to monitor program success. We used data from a large noninvasive genetic sampling effort conducted in 2004 and 33 years of physical captures to assess abundance, distribution, and genetic health of this population. We combined data from our 3 sampling methods (hair trap, bear rub, and physical capture) to construct individual bear encounter histories for use in Huggins—Pledger closed mark—recapture models. Our population estimate, Ň = 765 (95% CI = 715–831) was more than double the existing estimate derived from sightings of females with young. Based on our results, the estimated known, human‐caused mortality rate in 2004 was 4.6% (95% CI = 4.2–4.9%), slightly above the 4% considered sustainable; however, the high proportion of female mortalities raises concern. We used location data from telemetry, confirmed sightings, and genetic sampling to estimate occupied habitat. We found that grizzly bears occupied 33,480 km2 in the Northern Continental Divide Ecosystem (NCDE) during 1994–2007, including 10,340 km2 beyond the Recovery Zone. We used factorial correspondence analysis to identify potential barriers to gene flow within this population. Our results suggested that genetic interchange recently increased in areas with low gene flow in the past; however, we also detected evidence of incipient fragmentation across the major transportation corridor in this ecosystem. Our results suggest that the NCDE population is faring better than previously thought, and they highlight the need for a more rigorous monitoring program.

Determinants of Survival and Fecundity Through a Population Bottleneck in Pronghorn (Antilocapra americana)

The long-term viability of small, isolated populations has been questioned in light of stochasticity and bottlenecks. Following a drought in northwestern Montana in 2003, the pronghorn (Antilocapra americana )o f the National Bison Range experienced a demographic bottleneck. We used information theoretic approaches to examine the influence of sex, age, genetic variation, summer lactation, mate-search effort by females, and mating effort by males on subsequent winter survival and spring fecundity. Survival of males was influenced by age, whereas survival of females was influenced by prior energy expenditure and genetic variation, implicating inbreeding depression as the mechanism. Fecundity of females also was influenced by prior energy expenditure and genetic variation, with heterosis as the apparent mechanism. Our results agree with those of other studies that have emphasized the need to maintain genetic variation and limit inbreeding in small, isolated populations, and to account for stochasticity in population viability assessments and long-term management planning.

Dynamics of mycorrhizae during development of riparian forests along an unregulated river

In this study, we explore two mycorrhizal groups during development of riparian soils along a freely‐flowing river. We provide the first documentation of a shift in abundance between arbuscular mycorrhizae and ectomycorrhizae during floodplain succession. We used a chronosequence spanning 0–70 yr along a river in northwestern Montana, USA, to test the hypothesis that abundance of arbuscular mycorrhizal fungi (AMF) is greatest in early stages of soil development, and abundance of ectomycorrhizal fungi (ECMF) is greatest later in floodplain succession. We also measured the AMF‐mediated process of formation of soil aggregates during site development. AMF colonization of the dominant tree (black cottonwood,Populus trichocarpa) remained low (<5%), while AMF colonization of understory species was high (45–90%), across the chronosequence. Mycorrhizal inoculum potential (MIP) and hyphal length of AMF in soil peaked within the first 13 yr of succession and then declined. No single variable significantly correlated with AMF abundance, but AMF tended to decline as litter and soil organic matter increased. Density of ectomycorrhizal root tips in soil increased linearly throughout the chronosequence, and ectomycorrhizal colonization of cottonwood roots increased rapidly in early stages of succession. These patterns suggest that ECMF are not limited by dispersal, but rather influenced by abundance of host plants. Formation of water stable aggregates increased rapidly during the first third of the chronosequence, which was the period of greatest AMF abundance in the soil. The peak in AMF infectivity and hyphal length during early succession suggests that regular flooding and establishment of new sites promotes AMF abundance in this ecosystem. Regulation of rivers that eliminates creation of new sites may reduce contributions of AMF to riparian areas.

Effectiveness of aerial seeding and straw mulch for reducing post-wildfire erosion, north-western Montana, USA

Various methods are available to reduce post-wildfire erosion, but there is limited quantitative information on the relative effectiveness of these techniques. We used rainfall simulations to compare the erosion and runoff rates from adjacent 0.5-m2 plots treated with aerial grass seeding and straw mulch with untreated control plots following the July 2002 Fox Creek Fire in north-west Montana. In the first summer after the fire, plots seeded at a rate of 9 kg ha–1 had a mean of less than 5% ground cover and the seeding treatment had no effect on the rainsplash erosion rate. In contrast, straw mulch application at a rate of 2.24 Mg ha–1 resulted in ~100% ground cover and an 87% reduction in rainsplash erosion relative to the control (P = 0.001). Measurements on a subset of the plots in the second summer after the fire indicated that ground cover in the treatments and the control averaged 39%, and neither treatment provided a significant increase in ground cover or reduction in erosion relative to the control. These results add to the growing weight of evidence that straw mulch application is highly effective in reducing erosion in the first year after fire, whereas grass seeding is often ineffective because of the limited increase in ground cover that it produces.

Mapping forest alliances and associations using fuzzy systems and nearest neighbor classifiers

The study and management of biological communities depends on systems of classification and mapping for the organization and communication of resource information. Recent advances in remote sensing technology may enable the mapping of forest plant associations using image classification techniques. But few areas outside Europe have alliances and associations described in detail sufficient to support remote sensing-based modeling. Northwestern Montana has one of the few completed plant association classifications in the United States compliant with the recently established National Vegetation Classification system. This project examined the feasibility of mapping forest plant associations using Landsat Enhanced Thematic Mapper Plus data and advanced remote sensing technology and image classification techniques. Suitable reference data were selected from an extensive regional database of plot records. Fifteen percent of the plot samples were reserved for validation of map products, the remainder of plots designated as training data for map modeling. Key differentiae for image classification were identified from a suite of spectral and biophysical variables. Fuzzy rules were formulated for partitioning physiognomic classes in the upper levels of our image classification hierarchy. Nearest neighbor classifiers were developed for classification of lower levels (alliances and associations), where spectral and biophysical contrasts are less distinct. Maps were produced to reflect nine forest alliances and 24 associations across the study area. Error matrices were constructed for each map based on stratified random selections of map validation samples. Accuracy for the alliance map was estimated at 60%. Association classifiers provide between 54 and 86% accuracy within their respective alliances. Alternative techniques are proposed for aggregating classes and enhancing decision tree classifiers to model alliances and associations for interior forest types.

Invasive Pathogens At Alpine Treeline: Consequences for Treeline Dynamics

The potential impact of invasive forest pathogens on alpine treeline dynamics has not previously been considered. Whitebark pine (Pinus albicaulis), a foundation and keystone species of subalpine forests and major component of alpine treeline in the northern Rocky Mountains of the United States and southern Canada, is infected nearly range-wide by the exotic pathogen Cronartium ribicola, which causes white pine blister rust. A major component of treeline in the northern Rocky Mountains, whitebark pine initiates tree islands on the eastern slope in northwestern Montana more than any other conifer species. Blister rust infects whitebark pine throughout the region, and both infection and mortality are already evident at treeline. We discuss the cascading ecological effects of the loss of treeline whitebark pine and expected changes in landscape vegetation patterns. Potential implications of the loss of whitebark pine for northwestern Montana treelines are examined in the context of climate change within a conceptual model. We speculate that exotic pathogens could potentially confound predictions of treeline responses to global warming in many geographic regions.

Local Awareness of the Swift Fox Reintroduction Program in North Western Montana, USA

Local Awareness of the Swift Fox Reintroduction Program in North Western Montana, USA