Search
Search Results
-
This document is the 1998 annual progress report for studies of Pacific lampreys (Lampetra tridentata) conducted by the Confederated Tribes of the Umatilla Indian Reservation (CTUIR), Columbia River Inter-Tribal ...
Citation Citation
- Title:
- Pacific lamprey research and restoration project : annual report 1998
- Author:
- Close, David A.
This document is the 1998 annual progress report for studies of Pacific lampreys (Lampetra tridentata) conducted by the Confederated Tribes of the Umatilla Indian Reservation (CTUIR), Columbia River Inter-Tribal Fish Commission, and University of Minnesota (U of M). Bonneville Power Administration (BPA) funded activities through Project 94-026. The Pacific Lamprey Research and Restoration Project began after completion of a status report of Pacific lamprey in the Columbia River in 1995. The project started as a cooperative effort between the Confederated Tribes of the Umatilla Indian Reservation (CTUIR), Columbia River Inter-Tribal Fish Commission (CRITFC), and Oregon State University (OSU). Lamprey are a valuable subsistence food and cultural resource for Native Americans of the Pacific Northwest. The once abundant Pacific lampreys above Bonneville Dam are currently depressed (Close et al. 1995). Declines in Pacific lampreys have impacted treaty secured fishing opportunities by limiting tribal members catch and access to Pacific lampreys in the interior Columbia basin. Tribal members now harvest lampreys in lower Columbia River locations such as Willamette Falls near Oregon City, Oregon. Pacific lampreys are also an important part of the food web of North Pacific ecosystems, both as predator (Beamish 1980; Pike 1951; Roos and Gillohousin 1973), and prey (Semekula and Larkin 1968; Galbreath 1979; Roffe and Mate 1984; Merrell 1959; Wolf and Jones 1989) and as a vehicle for recruitment of marine nutrients. The decline of Pacific lampreys in the interior Columbia River basin has become a major concern. Effective recovery measures for Pacific lampreys can only be developed after we increase our knowledge of the biology and factors that are limiting the various life history stages. Prior to developing a restoration plan, we have carried out studies to review status, distribution, abundance, homing ability, and stock structure. These studies will culminate in the development and implementation of a restoration plan for the Umatilla River. Multiple pass electrofishing surveys to assess densities and distribution of lamprey larvae in the Umatilla River were conducted in 1998. Electrofishing surveys in the Umatilla River are useful for baseline comparison. Forty-two index sites were sampled from the mouth to river kilometer (RK) 124. Lamprey larvae were found in 4 of the 42 index plots. All sites with larvae were found at and below RK 9.3. Nine larvae were captured during the surveys. However, no larvae were caught on the second pass in each plot. Pacific lamprey larvae and adult lampreys were studied to determine their ability to produce and detect pheromones. Larval gall bladders were removed and gall bladder fluid was extracted and analyzed by high performance liquid chromatography (HPLC). Adult lampreys ability to detect pheromones were tested using electro-olfactogram (EOG) methods. Fifteen compounds including Petromyzonol sulfate (PS), a migratory pheromone found in sea lamprey larvae (Petromyzon marinus) (Li et al. 1995) were tested. Larval lampreys produced large amounts of (PS). Adult Pacific lamprey can detect PS and have an olfactory sensitivity to pheromones that is similar to sea lampreys. iv Pacific lamprey abundance, as indexed by fish ladder counts in 1998, was; Bonneville 37,478; The Dalles 7,665; John Day 12,579; McNary 3,393; Ice Harbor 763; Lower Monumental 69; Little Goose 90; Lower Granite 110; Rock Island 1,410; and Rock Reach 819 dams, respectively. Enumerating Pacific lamprey at counting stations remained extremely problematic, since excessive up- and downstream movement at the counting windows reduces the confidence in fish ladder passage estimates. This may be an indication of passage problems encountered by Pacific lampreys. In-season homing of Pacific lamprey was studied using radio telemetry. Pacific lampery were captured at Willamette Falls and Bonneville Dam, outfitted with radio transmitters and released approximately 26 km downstream of the Willamette River confluence. A total of 50 fish were instrumented. Results will be presented in next year’s report. Natal homing was also investigated using mtDNA analysis of fish captured at Bonneville Dam and from Willamette Falls. These results will also be presented next year. We collected lamprey tissues, from fish captured in several locations throughout the Columbia River Basin, to develop a genetic database for use in determining population structure. Additional samples for populations outside the Columbia River Basin were used to scale the results. Results from this investigation will be presented in next year’s annual report. Since the initiation of the CTUIR lamprey research and restoration project, additional lamprey studies have been proposed that have created uncertainties regarding the prioritization of projects and needs of lampreys. At the request of the Northwest Power Planning Council, a multi-agency Pacific lamprey technical workgroup (TWG) was established in 1996. Annual meetings are held to coordinate projects and prioritize research needs. The TWG identified critical uncertainties and needs to help in determining priorities of ongoing and proposed projects (Appendix A). Finally, an annotated bibliography of relevant lamprey literature was compiled (Appendix B).
-
2372. [Article] Students as pro-social bystanders : opportunities, past behaviors, and intentions to intervene in sexual assault risk situations
Sexual assault is a major public health concern in the U.S, and college students are particularly vulnerable to victimization. A health issue that affects nearly one in four women (Fisher, Cullen, & Turner, ...Citation Citation
- Title:
- Students as pro-social bystanders : opportunities, past behaviors, and intentions to intervene in sexual assault risk situations
- Author:
- Hoxmeier, Jill C.
Sexual assault is a major public health concern in the U.S, and college students are particularly vulnerable to victimization. A health issue that affects nearly one in four women (Fisher, Cullen, & Turner, 2000; Karjane, Cullen, & Turner, 2005) and that is associated with severe negative health outcomes, including depression substance abuse, suicide ideation, and risky sexual behaviors (CDC, 2012), warrants effective prevention programs. Moving away from traditional prevention efforts, which target females as potential victims in risk reduction programs and males as potential perpetrators in attitudinal-shifting programs, bystander engagement programs have become increasingly more widespread. These programs aim to engage all students on the college campus as potential bystanders who can intervene to prevent a sexual assault or reduce the harm of an assault that has already occurred (Banyard, Moynihan & Plante, 2007). Burn (2009) investigated potential barriers to pro-social bystander intervention using the Situational Model of Bystander Intervention, a model based on the original research of bystander behavior of Latanè and Darley (1970). The model outlines five barriers that influence students' intent to intervene as witnesses to sexual assault: failure to notice the situation, failure to identify the situation as high risk, failure to take intervention responsibility, failure to intervene due to skills deficit, and failure to intervene due to audience inhibition (Burn, 2009). She found that students' perception of barriers negatively correlated with intervention behaviors as bystanders to sexual assault (Burn, 2009). Although bystander engagement programs have shown initial promise in increasing students' intent to intervene, more needs to be known about the opportunities students have to intervene, their past intervention actions, and their intent to intervene in the future across the wide range of situations that encompass sexual assault risk. In addition, to develop effective programs that aim to increase pro-social behavior, understanding the salient influences of students' intent is critical. This study uses the Theory of Planned Behavior (TPB; Ajzen & Fishbein, 1991) to examine the influences of students’ intent to perform 12 different pro-social bystander behaviors. The TPB asserts that individuals' behavior is most proximally influenced by their behavioral intentions, and intentions are influences by their perceived behavioral control to perform the behavior, subjective norms that support performing the behavior, and attitudes toward the behavior (Ajzen & Fishbein, 1991). The four primary aims of this study were: 1) to examine the demographic correlates of students' opportunities, past intervention actions, and reported intent to intervene; 2) to examine any differences in students' intent to intervene based on the level of intervention (pre-, mid-, and post-assault) and type of intervention (with the potential or actual victim compared to the potential or actual perpetrator); 3) to examine the influences of perceived behavioral control, subjective norms, and attitudes on students' intent to intervene as bystanders; and 4) to compare the TPB-based model to the Situational Model of Bystander Intervention (Burn, 2009) in its ability to explain students' intent to intervene as bystanders. In the Fall of 2014, a sample of 815 undergraduate students at Oregon State University completed the Sexual Assault Bystander Behavior Questionnaire (SABB-Q), a tool comprised of items to measure students' opportunities, past behaviors, and future intent, in addition to measures assessing the influences of students' intent in line with the Theory of Planned Behavior and Burn's (2009) Situational Model of Bystander Intervention. Students who participate in Greek communities (fraternities and sororities) reported significantly greater odds of having the opportunity to perform four of the 12 intervention behaviors compared to non-Greek students, while student-athletes reported significantly greater odds of having the opportunity to perform two of the 12 intervention behaviors. Females reported significantly more past pro-social intervention behaviors (x̄ = 0.87) compared to males (x̄ = 0.79; p = 0.007). Regarding intent to intervene in the future, females reported significantly greater intent to intervene compared to males (x̄ = 6.07 vs. 5.68; p = 0.007). Students with friends who have been victims of sexual assault reported greater intent to intervene compared those without friends who have been victims (x̄ = 6.04 vs. 5.89; p = 0.02). Students with a personal history of victimization reported significantly greater intent compared to those without a personal history (x̄ =6.13 vs 5.93; p = 0.03). Students reported significantly greater intent to intervene with the potential or actual victim compared to the potential or actual perpetrator (x̄ = 6.19 vs. 5.74, p < 0.001). Females reported significantly greater intent to intervene with both the potential or actual victims and perpetrators (x̄ = 6.31 and 5.84, respectively) compared to males (x̄ = 5.88 and 5.49, respectively). Both males and females reported the greatest intent to perform post-assault intervention behavior (x̄ = 6.23), followed by pre-assault (x̄ = 6.08) and mid-assault behaviors (x̄ = 5.57). Females reported significantly greater intent to perform nine of the 12 pro-social intervention behaviors compared to males. A multiple regression analysis revealed that perceived behavioral control, subjective norms, and attitudes explained a significant proportion of the variance in intent to intervene (R² = 0.55, F(3, 771) = 315.68, p < 0.000). Perceived behavioral control was highly significant (β = 0.48, p < 0.001), as were subjective norms (β = 0.15, p < 0.001) and attitudes (β = 0.30, p < 0.001). Gender differences were also observed. For females, perceived behavioral control was highly significant (β = 0.49, p < 0.001), as were subjective norms (β = 0.15, p < 0.001) and attitudes (β = 0.29, p < 0.001). For males, perceived behavioral control was highly significant (β = 0.49, p < 0.001), as were attitudes (β = 0.29, p < 0.001). However, males' subjective norms were not significantly related (β = 0.07, p = 0.199) to their intent to intervene. Further analysis revealed a significant interaction between gender and subjective norms (β = -0.28; p = 0.039). The TPB-based model including this moderation effect explained a significant proportion of the variance in students' intent to intervene (R² = 0.57, F(6, 766) = 168.46, p < 0.000). Interveners reported significantly greater perceived behavioral control than non-interveners for seven of the 12 intervention behaviors; more supportive subjective norms than non-interveners for six of the 12 intervention behaviors; more positive attitudes than non-interveners for only one of the 12 intervention behaviors; and greater intent to intervene in the future for six of the 12 intervention behaviors. However, differences in the three TPB variables between interveners and non-interveners were not consistent for the 12 intervention behaviors. Regarding Burn's (2009) Situational Model of Bystander Intervention, a multiple regression analysis revealed two of the five barriers were significantly related to students’ intent to intervene: the failure to take intervention responsibility barrier (β = -0.29, p < 0.001) and the failure to intervene due to audience inhibition barrier (β = -0.22, p < 0.001). The model in whole explained a large proportion of the variance (R2 = 0.25, F(5, 768) = 50.14, p < 0.000). Gender differences were also observed. For females, failure to take intervention responsibility (β = -0.23; p < 0.000) and failure to intervene due to audience inhibition (β = -0.23; p < 0.001) both had a significant, negative influence on their intent to intervene. For males, failure to take intervention responsibility (β = -0.21; p < 0.014) had a significant, negative influence on intent to intervene. Additional analysis revealed no significant interactions between gender and any of the five barriers. The TPB-based model explained a greater proportion of the variance (R2 = 0.55) compared to Situational Model of Bystander Intervention (R² = 0.25) in the multiple regression analysis using all 12 intervention behaviors. All three variables in the TPB-based model were significantly related to students’ intent, whereas only two of the five barriers were significantly related. A final multiple regression analysis was conducted using all three significant TPB variables and the two significant barriers to explain students' intent to intervene. The combined model explained a significant proportion of variance in students' intent (R² = 0.58 F(5, 756) = 206.19, p < 0.000) and significantly improved upon the TPB-based model (Δ R² = 0.03; p < 0.000). The results of this study have several implications for future research and public health practice. First, it is important to ask students about their opportunities to intervene in addition to their actual intervention behaviors because this information helps paint a clearer picture of bystander engagement. This assessment could also help identify high-risk groups: students who have greater opportunities to intervene as bystanders and/or report fewer intervention behaviors compared to their reported opportunities. Second, students may conceptualize intervention behaviors differently depending on the phase of the assault and with whom the intervention behavior requires intervening. Accordingly, programs aimed at encouraging students to intervene should take these differences into consideration. Third, the Theory of Planned Behavior, used to explain and change other health-related behaviors, can effectively be applied to help uncover determinants of pro-social bystander behaviors. Perceived behavioral control, subjective norms, and attitudes appear to be salient influences in students' intent to intervene. Therefore, bystander engagement programs should incorporate activities to heighten students' skills to intervene, change social norms that support bystander intervention, and shift attitudes toward the benefits of intervening. This study demonstrates the importance of using an established, evidenced-based theoretical framework to explain behavioral influences and strengthens the argument for continued use of theory to identify, and potentially change, salient influences in behavioral performance. Students as pro-social bystanders have the potential to make a positive impact on the reduction of sexual assault on the college campus. Although the responsibility for sexual assault rests on those who perpetrate such acts, and primary prevention strategies aimed at those demonstrating a risk for perpetration are imperative, sexual assault is a public health issue that warrants a multi-pronged approach to reduce its incidence and migrate its associated harms. Programs that engage students as pro-social bystanders have the potential to make a positive impact on the reduction of sexual assault incidence in the absence of effective primary prevention strategies. The findings of this study make a contribution to the literature examining influences of students' pro-social bystander intervention to sexual assault situations and provide suggestions for strategies to increase bystander engagement.
-
2373. [Article] Evaluating tools used to estimate and manage browse available to wintering moose on the Copper River Delta, Alaska
Ungulates comprise some of the most well researched and intensely managed wildlife populations on earth. As such, they are recognized as ideal study subjects for developing and modifying management tools ...Citation Citation
- Title:
- Evaluating tools used to estimate and manage browse available to wintering moose on the Copper River Delta, Alaska
- Author:
- Smythe, Sharon E.
Ungulates comprise some of the most well researched and intensely managed wildlife populations on earth. As such, they are recognized as ideal study subjects for developing and modifying management tools or theories (Danell et al. 1994, Shipley 2010). An introduced moose (Alces alces gigas) population on the Copper River Delta (CRD; Delta) of south-central Alaska functions as a valuable resource for the residents of Cordova and an isolated research population on which to test managerial techniques. Since its introduction (1949-1958), the founding population of 23 moose has grown to over 830 in 2013, divided into two sub-populations. However, in 1964, the largest earthquake recorded in U.S. history (9.2 magnitude) uplifted the Delta by 1.0-4.0 m, initiating delta-wide changes in hydrology, vegetation distributions, and successional processes. The proportion of stands dominated by woody species, especially alder (Alnus viridis sinuata) and spruce (Picea sitchensis), increased visibly. Furthermore, previous research (1987-89) observed that 90% of the moose on the western region of the Delta wintered within 9-24% of the total land area, possibly restricting their available winter browse. Because moose diets on the CRD are dominated by willows (Salix spp.), managers were concerned that the combined effects of a restricted winter range and earthquake-initiated vegetation changes would negatively influence the population. Managers have responded to this concern by 1) supporting work to estimate the nutritional carrying capacity (NCC; i.e., the forage available to a population within a specified area and time) of the Delta, and 2) by exploring the feasibility of mechanical treatment as a means of stimulating browse production for the moose. Thus, the objectives of this thesis were to 1) explore the factors influencing NCC for moose on the west CRD while combining updated digital and field-collected data to estimate NCC, and 2) to assess the effects of mechanical treatment on the production of moose browse across stand types and over time. We collected field data and evaluated differences in the past (1988-89) and present (2012-13) biomass-predicting regression equations for two willow species (Barclay's and Hooker's willow, Salix barclayi and hookeriana) used within NCC models to determine 1) whether past and present models of Barclay's willow predicted similarly and 2) whether Hooker’s and Barclay's willow differed in average available biomass, nutritional quality, and utilization by moose. The linear coefficients for the current (2013) Barclay’s willow, Hooker's willow, and combined equations were 2.2x, 1.6x, and 1.9x larger, respectively, than that derived from the 1988 model for Barclay’s willow (which possibly included Hooker's willow data). Thus, willows on the CRD may now be supporting more biomass per stem than predicted by prior models. Hooker's and Barclay's willow did not differ in mean available biomass, nutritional values, or utilization rates. These results suggest a need to evaluate the accuracy of older allometric regressions, though separate identification of the visually-similar Barclay's and Hooker's willow may be unnecessary for future biomass-estimating efforts on the CRD. To further explore the factors influencing the biomass available to moose and their associated NCC estimates, we compared 5 NCC model types across 4 winter ranges and under 3 winter-severity scenarios for the western CRD moose population. We conducted a sensitivity analysis (Sx) of our final model to determine the relative influence of factors affecting NCC estimates. Lastly, we compared current (2012-13) browse available biomass, stand type areas (2011), and NCC results to those obtained by past research (1987-89, MacCracken et al. 1997 and 1959/1986, Stephenson et al. 2006) to determine changes over time, while evaluating the effects of models incorporating satellite-based estimates of stand areas and forage nutritional values on NCC estimates. Because recent aerial survey observations suggest expansion of the moose winter range, our final model estimated NCC between 2,198-3,471 moose depending on winter severity within a winter range encompassing the entire west Delta. These results suggest the current western moose population (approximately 600 in 2013) is below NCC. Model components with the largest and smallest Sx were snow depth and tannin- and lignin-caused reductions in forage nutritional quality, respectively. Changes from 1987-2013 in available biomass of forage species ranged from -66-493%, while changes from 1959-2011 in stand type areas ranged from -60-661%. Overall, NCC estimates only declined by 2% from 1959-2013, however inclusion of forage nutritional quality in models reduced NCC estimates by 60%. Lastly, we assessed the use of hydraulic-axing (i.e., hydro-axing) as tool for increasing the available willow biomass. We evaluated treatment effects on biomass, height, nutritional quality (crude protein, lignin, and tannin levels), utilization, and snow burial of the winter forage species within 3 winter-severity scenarios. Sites were treated in 4 winters (1990-92, 2008, 2010, and 2012) within 5 stand types in 20 locations varying from 0.86-63.40 ha in size. Results indicate few significant differences relative to controls, though treatment significantly increased the ratio of willow to alder. Our results may be limited by sample sizes (n = 1-9 per stand type or treatment year) as visual comparison suggests treatment via hydraulic-axing may be an effective method for increasing willow biomass without influencing nutritional quality. However, willows 20-23 years post-treatment are still significantly shorter than untreated willows (P = 0.03). Thus, treatment may result in decreased forage available to moose in severe winters. Management concerns regarding continued earthquake-initiated changes in vegetation distributions and successional processes prompted our investigation. However, studies on the vegetation dynamics of the CRD suggest the vegetation distribution of the Delta may be relatively stable (Boggs 2000, Thilenius 2008). If so, our current estimates of NCC suggest the west Delta can support a larger moose population than is currently present. Hydro-axing may not be necessary to ensure the future of the moose population, though it could be used to counter increases in alder (Alnus viridis sinuata) which are likely within certain successional sequences. However, together with earthquakes, the processes determining the future forage available on the CRD include complex, interactive forces such as glaciers, the Copper River, oceanic tides, and zoological- and human-caused influences. These forces and their effects on the vegetation create a dynamic ecosystem for the moose population, are difficult to predict, and may be further complicated by climate change. As a result, application of any managerial tool may be temporary and often difficult. However, this guarantees a constant need for further revision and redevelopment of the tools used to manage the moose population, ensuring that the moose of the CRD will remain an important resource for researching and refining ungulate management worldwide.
-
The forest alpine tundra ecotone (FTE, also known as alpine treeline or subalpine parkland), is a conspicuous feature of mountain landscapes throughout the world. Climate change-driven increases in temperature ...
Citation Citation
- Title:
- Patterns of tree establishment and vegetation composition in relation to climate and topography of a subalpine meadow landscape, Jefferson Park, Oregon, USA
- Author:
- Zald, Harold Samuel James
The forest alpine tundra ecotone (FTE, also known as alpine treeline or subalpine parkland), is a conspicuous feature of mountain landscapes throughout the world. Climate change-driven increases in temperature are believed to result in FTE movement and tree invasion of subalpine meadows, which have been documented throughout the Northern Hemisphere across a wide range of geographic locations, climatic regimes, forest types, land use histories, and disturbance regimes. Climate-driven FTE movement may have numerous ecological effects such as: positive temperature feedbacks, increased net primary productivity and carbon storage, and declines of plant populations and species. The magnitude of these ecological effects is highly uncertain, but will be largely determined by the rates and spatial extent of FTE movement and meadow invasion. FTE movement and meadow invasion are often considered at global or regional spatial scales in relation to climate, yet they are fundamentally driven by tree regeneration processes that are influenced by a variety of climatic and biophysical factors at micro site, landscape, and regional scales. Much of the research on the FTE has not taken a landscape approach incorporating multi-scale processes. For example, species distribution models used to project climate change effects on future species distributions and plant biodiversity in mountainous landscapes rely on species distribution data that is often sparse and incomplete across FTE landscapes. This dissertation attempts to overcome many of the limitations in FTE research by taking a landscape approach to develop a greater understanding of past spatiotemporal patterns of tree invasion, current spatial patterns of vegetation composition and structure, and potential future patterns of climate-driven tree invasion in the FTE. The setting for this research is Jefferson Park, a 260 ha subalpine parkland landscape in the Oregon High Cascades, USA. This study uses field plots, remotely sensed imagery, airborne Light Detection and Ranging (LiDAR), and simulation modeling to: 1) predictively map current fine-scale species distributions, vegetation structure, and tree ages; 2) reconstruct patterns of tree invasion over the last fifty years in subalpine meadows in relation to climatic conditions, landforms, microtopography, and seed dispersal limitations; and 3) develop a statistical model that projects future patterns of tree invasion into subalpine meadows under different climate scenarios in Jefferson Park. In chapter two, I generated fine-scale spatially-explicit predictions of current vegetation composition, structure, and tree ages in the Jefferson Park study area. Objectives of this chapter were threefold: 1) to characterize spatial patterns of tree ages, vegetation composition, and vegetation structure in a FTE landscape in the Oregon Cascades using predictive mapping; 2) determine how vegetation composition and structure were associated with gradients of environmental factors derived from multispectral satellite imagery and Light Detection and Ranging (LiDAR) data; and 3) determine if predictive mapping characterizations of tree age, vegetation composition, and vegetation structure were improved by the inclusion of LiDAR data. Predictive mapping of vegetation attributes was accomplished using gradient analysis with nearest neighbor imputation; integrating field plots, multispectral SPOT 5 satellite imagery, and LiDAR data. Vegetation composition was best described by SPOT 5 imagery and LiDAR-derived topography, while vegetation structure was best described by LiDAR-derived vegetation heights. Predictions of species occurrence were most accurate for tree species, moderate for shrub species and vegetation groups, and highly variable for graminoid species. Tree age, which was the most accurately predicted vegetation structure variable, indicates the study area was largely un-forested in 1600, gradually invaded by trees from 1600 to the 1920's, and rapidly invaded from the 1920's to 1980. Predictive mapping of vegetation structure variables such as basal area and stand density were subject to large amounts of error, possibly resulting from scale incompatibilities between vegetation patterns and plot size, and/or heterogeneous FTE landscapes where forest structure does not develop along consistent trajectories with stand age. This study suggests integrating multispectral satellite imagery, LiDAR data, and field plots can accurately predict fine-scale spatial characterizations of species distributions and tree invasion in the FTE. This study also indicates that sample design can influence spatial patterns of model uncertainty, which needs to be considered if predictive mapping of vegetation and sensitive ecosystems is a component of inventory and monitoring programs. In chapter three, I focused on quantifying spatiotemporal patterns of subalpine parkland tree invasion in Jefferson Park over the past five decades in relation multi-scale climatic and biophysical controls. LiDAR data provided previously unavailable fine-scale spatial characterizations of microtopography and vegetation structure. I utilized LiDAR, georeferenced field plots, and tree establishment reconstructions to quantify spatiotemporal patterns of tree invasion in relation to late season snow persistence, landform types, fine-scale topographic variability, distances from potential seed sources, and climate variation within 130 ha of the subalpine parkland landscape of Jefferson Park. Tree occurrence (i.e. tree presence in 2 m plots and grid cells) occurred in 7.75% of study area meadows in 1950 and increased to 34.7% in 2007. Landform types and finer-scale patterns of topography and vegetation structure influenced summer snow depth, which influenced temporal and spatial patterns of tree establishment. Tree invasion rates were higher on debris flow landforms, which had lower summer snow depth than glacial landforms, suggesting potentially rapid treeline responses to disturbance events. Tree invasion rates were strongly associated with reduced annual snow fall on glacial landforms, but not on debris flows. Tree establishment was spatially constrained to micro sites with high topographic positions and close proximity to overstory canopy, site conditions associated with low summer snow depth. Seed source limitations placed an additional species-specific spatial constraint on where trees invaded meadows. Climate and topography had an interactive effect, with trees establishing on higher topographic positions during both high snow/low temperature and low snow/high temperature periods, but had greater than expected establishment on lower topographic positions during low snow/high temperature periods. Within the context of larger landform types, topography and proximity to overstory trees constrained where trees established in the meadows, even during climate periods with higher temperatures and lower snowfall. Results of this study suggest large scale climate-driven models of vegetation change may overestimate treeline movement and meadow invasion, because they do not account for biophysical controls limiting tree establishment at multiple spatial scales. In chapter four, I used field data and analyses from chapter 3 to parameterize a spatially and temporally explicit statistical model of fine-scale tree invasion within 130 ha of the Jefferson Park study area. The model incorporated both the climatic and biophysical controls found in chapter 3 to influence tree invasion. The model was used in two ways: (1) to spatially project patterns of tree invasion from 1950 to 2007 in response to historical climate; and (2) to project future tree invasion of the study area from 2007 to 2064 under six different annual snowfall scenarios. Modeling addressed the following questions: (1) Can fine-scale (2 m pixel size) patterns of historical tree invasion be accurately predicted? (2) How sensitive is future tree invasion (and therefore meadow persistence) to different future snowfall scenarios? (3) Are non-climatic factors such as landforms and biotic interactions associated with different spatial patterns of tree invasion? From 1950 to 2007, simulated historical meadow area declined from 82% to 65% of the study area. Model outputs of historical area, spatial distributions, and spatial clustering of tree invasion generally agreed with independent validation, and suggest biotic interactions due to young tree establishment facilitation are important on glacial landforms but not debris flows. Simulations of future scenarios indicated meadow declined to 36 to 43% of the study area by 2064. Projected meadow area declined with reduced annual snow fall, but not under prolonged high and low snow fall periods. Meadows persisted under all future scenarios in 2064. This model suggests subalpine meadows may significantly decline under climate warming, but will still persist in 2064. Micro sites and recruitment limitation may be equally or more important factors than climate change in influencing subalpine landscape change, suggesting local high-elevation persistence of subalpine meadows under future climate warming.
-
The ~1 Myr history of the Purico-Chascon volcanic complex (PCVC) records significant changes in the production and storage of magmas in the crust. At ~1 Ma activity at the PCVC initiated with the eruption ...
Citation Citation
- Title:
- Crustal architecture and magma dynamics in a large continental magmatic system : a case study of the Purico-Chascon Volcanic Complex, Northern Chile
- Author:
- Burns, Dale H.
The ~1 Myr history of the Purico-Chascon volcanic complex (PCVC) records significant changes in the production and storage of magmas in the crust. At ~1 Ma activity at the PCVC initiated with the eruption of a large 80-100 km³ crystal-rich dacite ignimbrite with restricted whole rock ⁸⁷Sr/⁸⁶Sr isotope ratios between 0.7085-0.7090. In-situ analyses of plagioclase from the Purico ignimbrite have ⁸⁷Sr/⁸⁶Sr=0.7087-0.7090. The dacite magma accumulated and evolved at relatively low temperatures around 800-850 °C in the upper crust at 4-8 km depth. Minor andesite and rhyolite pumice late in the sequence have similar restricted whole rock ⁸⁷Sr/⁸⁶Sr=0.7089-0.7091. The radiogenic isotopes of this 0.98 Ma activity are consistent with all these compositions resulting from 50 to 70% crustal assimilation by parental Central Andean "baseline" magmas at depths between 15-30 km. The final eruptions at the PCVC occurred <0.18 Ma producing three small < 5 km³ crystal-rich dacite lava domes with whole rock ⁸⁷Sr/⁸⁶Sr ratios 0.7075 to 0.7081 containing abundant basaltic-andesite enclaves with whole rock ⁸⁷Sr/⁸⁶Sr ratios of 0.7057- 0.7061. Plagioclase and amphibole from samples from the largest of these domes, Cerro Chascon, record two distinct magmatic environments; an upper crustal environment identical to the Purico ignimbrite and a second deeper, ~15-20 km depth, higher temperature (~922-1001 °C) environment consistent with conditions recorded in the basaltic andesite enclaves. Accordingly, plagioclase cores in the host dacite lava and enclaves have enriched in-situ ⁸⁷Sr/⁸⁶Sr isotopic compositions of 0.7083 to 0.7095 while plagioclase rims and microphenocrysts in the enclaves have ⁸⁷Sr/⁸⁶Sr isotope ratios from 0.7057 to 0.7065 and 0.7062 to 0.7064 respectively. Lavas from Cerro Chascon also contain abundant Fo82 olivine with spinel and basaltic melt inclusions that crystallized in a deep crustal environment (>1250 °C) consistent with a lower crustal MASH zone. The high baseline isotopic ratios observed in bulk rock and plagioclase crystals from Cerro Chascon (0.7057-0.7065) are consistent with MASH processes. The evolution of the PCVC is a microcosm of the Andean arc in this region where, from 10 - 1 Ma, dominantly dacitic upper crustal magmatism of the Altiplano-Puna Volcanic Complex ignimbrite flare-up persisted until ~1 Ma, when smaller volume, more heterogeneous and less isotopically enriched basaltic andesite to dacite composite volcanoes signal a return to steady state arc volcanism. I suggest that the PCVC captures the transition of the Andean arc from flare-up to steady state. The temporal trend at the PCVC is consistent with a waning thermal flux. High magmatic fluxes during the flare-up would have resulted in elevated geothermal gradients and efficient crustal processing leading to a dominantly dacitic upper crust (0 to 35 km) that fed the large volume Purico ignimbrite. As magmatic flux and thermal energy wanes, crustal isotherms relax resulting in greater thermal contrast between parental magmas, crust and remnant upper crustal dacite magma. This manifests in more heterogeneity and the survival of less isotopically enriched magmas in the upper crust. These arc scale magma dynamics are recorded even at the intra-crystalline scale. Individual crystals from Cerro Chascon also record vital information on the crystallization and evolution of mantle-derived magmas in continental magmatic arcs. Fo₈₂ olivine, olivine hosted spinel, and basaltic melt inclusions record the crystallization of olivine at >1250 °C in conditions consistent with a lower crustal (~70 km depth) MASH zone. Another significant crystallization event appears to have occurred at ~20 km depth, characterized by the crystallization of high An plagioclase (An₇₂₋₈₄) at ~1100-1050 °C followed by high-Al amphibole (~12-15 wt.% Al₂O₃) at ~1000-950 °C. The appearance of amphibole on the liquidus appears to have resulted from a nearly 2-fold increase in melt water content following ~45% crystallization of high An plagioclase. Following this extensive crystallization the highly crystalline mafic magma ascended into the upper crust and interacted with the remnant crystal mush from the Purico ignimbrite magma reservoir. Low An plagioclase (An₃₉₋₅₅), low Al amphibole (~6-9 wt.% Al₂O₃), sanidine, and biotite retain the chemical composition of the Purico ignimbrite magma, whereas, olivine, high An plagioclase, and high Al amphibole record the mafic recharge magma. The textures and compositions observed in Cerro Chascon are common in both continental and oceanic magmatic arcs worldwide and I propose that multiple crystallization events and upper crustal assimilation are fundamental processes intrinsic to arc magmatism. I have also used in situ ⁸⁷Sr/⁸⁶Sr isotope ratios in plagioclase from andesite, dacite, and rhyolite pumice from the ~1 Ma Purico ignimbrite to determine the cause for compositional zoning in the Purico ignimbrite magma reservoir. Andesite pumice contains two texturally, compositionally, and isotopically distinct types of plagioclase, small (<500 μm) subhedral to euhedral crystals with high MgO (130-490 ppm) and low ⁸⁷Sr/⁸⁶r crystals (0.7076-0.7084) record a hot (>900 °C) andesite magma derived from an ~20 km deep magma reservoir. In contrast, the second type of plagioclase in the andesite appear to broken fragments of larger crystals and have significantly lower MgO (90-240 ppm), higher ⁸⁷Sr/⁸⁶Sr (0.7096-0.7114), and appears to be derived from the lower temperature (crystallized at ~800-900 °C), upper crustal (<10 km) plutonic basement. Dacite pumice also contains two texturally and compositionally distinct types of plagioclase. However, both types have very restricted MgO (b.d.l.-200 ppm) and ⁸⁷Sr/⁸⁶Sr (0.7085-0.7095) ratios and appear to have grown at ~850°C. These crystals are also significantly larger (>1000 μm) than plagioclase from the andesite pumice and have clear euhedral rims. Rhyolite pumice from the Purico ignimbrite also contains distinct types of plagioclase. Both types of plagioclase are similar in size (<500 μm) and appear to be fragments of larger crystals. One type is characterized by low MgO (b.d.l.-240 ppm) and restricted ⁸⁷Sr/⁸⁶Sr isotope ratios (0.7088-0.7095) similar to plagioclase in the dacite pumice, and the other has significantly higher ⁸⁷Sr/⁸⁶Sr ratios (0.7095-0.7103) consistent with the upper crustal ignimbrite basement. The compositional variations observed in plagioclase crystals from the Purico ignimbrites are consistent with the recharge of a previously emplaced upper crustal (4-8 km depth) dacite magma reservoir by a hotter, deeper (20 km deep) andesite. During ascent, the andesite incorporated crystals from the surrounding upper crustal plutonic bodies before pooling against the residence dacite magma and crystallizing. Crystallization of the andesite resulted in the expulsion of a rhyolite interstitial melt that ascended through the dacite reservoir and pooled at the top of the reservoir. The rhyolite melt incorporated crystals from the dacite magma during ascent as well as crystals from the roof rock, which in the case of the Purico ignimbrite represents the plutonic remnants from other large silicic magmatic systems associated with the APVC. Thus, the compositional variations observed in the Purico ignimbrite results from a combination of crustal assimilation, crystallization, and melt extraction all initiated by mafic recharge.
-
Geologic mapping of the Longview-Kelso area and the measurement and description of a composite 650-meter thick stratigraphic section of the Cowlitz Formation (Tc) in Coal Creek using bio-, magneto-, litho-, ...
Citation Citation
- Title:
- Stratigraphy and sedimentology of the middle eocene Cowlitz Formation and adjacent sedimentary and volcanic units in the Longview-Kelso area, southwest Washington
- Author:
- McCutcheon, Mark S.
Geologic mapping of the Longview-Kelso area and the measurement and description of a composite 650-meter thick stratigraphic section of the Cowlitz Formation (Tc) in Coal Creek using bio-, magneto-, litho-, and sequence stratigraphy reveals a complex interplay of Cowlitz micaceous, lithic arkosic shelf to tidal/estuarine to delta plain facies associations, and Grays River basalt lava flows and interbedded basalt volcaniclastics from nearby Grays River eruptive centers (e.g., Mt. Solo and Rocky Point). The lower 100 meters of the Coal Creek section (informal unit 1, Chron 18r) consists of micaceous, lithic arkosic sandstone and siltstone and minor coals, was deposited as part of a highstand system tract (HST) at the base of 3rd order cycle number 3. This unit consists of four dominantly tidal shoaling-upward arkosic sandstone parasequences reflecting upper shoreface to delta plain depositional environments. The overlying unit 2 (Chron 18n) is defined by abundant Grays River basalt volcaniclastic interbeds that intertongue with Cowlitz lithic arkoses. This unit represents the latter part of 3rd order cycle 3, and consists of mostly fining- and thinning-upward parasequences of middle shoreface to delta plain successions of an aggradational to transgressive parasequence set. Near the top of unit 2 is a maximum marine flooding surface depositing lower shoreface lithic arkosic sandstone to shelf siltstones over upper shoreface micaceous lithic arkose. Unit 3 comprises 3rd order cycle 4 (Chron 17r), a lowstand system tract, and consists of 6 mostly fining- and thinning-upward parasequences of lower shoreface to delta plain facies associations. A parasequence or erosional boundary at the base of unit 5 (Chron 17r) consists of submarine channel-fill scoured into underlying micaceous siltstones, produced during a lowstand system tract (LST) of 3rd order cycle 5. This deep marine channel-fill sequence is overlain by thinlybedded to laminated overbank distal turbidites and hemipelagic siltstones that define the top of the Coal Creek section. These 5 informal units in Coal Creek lithologically and chronologically correlate to 5 similar informal units defined by Payne (1998) in the type section of Cowlitz Formation in Olequa Creek near Vader -30 km to the north. Middle Eocene Grays River Volcanics of the study area are mapped as two separate units: a lower unit over 150 meters thick in places, consisting of subaerial basaltic flows and invasive flows (Tgvl), intrusions (Tgvis and Tgvid), and volcaniclastics (Tgvsl); and an upper unit consisting of commonly mollusk-bearing, shallow marine basaltic sedimentary interbeds that intertongue with the Cowlitz Formation (Tgvs2), particularly Cowlitz unit 2 of the Coal Creek section. These volcaniclastic deposits are intrabasinal, derived from volcanic highlands to the west and northwest, and local phreatomagmatic tuff cones. The lower Grays River volcaniclastic unit typically overlies Grays River flows in the study area and is divided into 5 informal facies. Geochemically, Grays River flows in the study area fall within normal parameters (3 to 4% TiO2 and high iron tholeiitic basalts). However, basalt flows and bedded scoriaceous breccias near Rocky Point are anomalously low in TiO2 and are considered in this study to be a separate volcanic subunit (Rocky Point Basalts), time equivalent to and interfingering with Grays River lavas, but may represent mixing with shallower western Cascade calc-alkaline magma. Over 60 younger Grays River dikes intrude the Cowlitz Formation in Coal Creek. A dike low in the Coal Creek section is dated at 40 ± 0.36 Ma, and an invasive flow at Mt. Solo is dated at 36.98 ±.78 Ma. Volcanics capping the hills east of the Cowlitz River are chemically distinct as slightly younger western Cascade basaltic andesite flows, and two dikes east of the river are chemically distinct as western Cascade andesite. Overlying Grays River Volcanics and Cowlitz Formation in much of the study area, are clayey and commonly tuffaceous siltstones and silty sandstones, possibly of the late Eocene-early Oligocene Toutle Formation, a new unit to this area. The Toutle Formation is a mixture of wave and stream reworked micaceous and arkosic Cowlitz Formation and fresh silicic pyroclastic ash and pumice from the active western Cascade arc. An angular unconformity separates the Paleogene Grays River Volcanics, Cowlitz Formation, and Toutle Formation from the early to middle Miocene Columbia River Basalt Group. Based on lithology, geochemistry, stratigraphic relationships, and magnetic polarity, 6 individual Columbia River Basalt flows have been mapped in this study. The three lower Grande Ronde flows are of normal polarity and Ortley low MgO chemical composition. The lowermost flow (N2 Ortley #1) is absent in the Columbia Heights area, low MgO, about 10 meters thick and consists of pillow-palagonite sequences in the upper quarry on Mt. Solo. Aphyric N2 Ortley flow #2 is over 35 meters thick with well-developed upper and lower colonnade, and of intermediate MgO. N2 Ortley flow #3 is pillow-palagonite in the Storedahl Quarry and low MgO. A -4-meter thick tuffaceous overbank siltstone and basalt conglomeratic channel interbed separates the three low MgO Ortley flows from the overlying high MgO N2 Grande Ronde Sentinel Bluffs flow. A single exposure of well-developed large colonnade with sparse 1 cm labradorite laths, and reddish oxidized soil, defines the N Sand Hollow flow of the Frenchman Springs Member of the Wanapum Formation. The overlying Pomona Member is mapped based on previous work by other authors. Pliocene gravels and arkosic sand of the Troutdale Formation form upland terrace deposits up to 100 meters thick in southern parts of the study area, and represent the uplifted paleo-thalweg and overbank flood deposits of the downcutting, antecedent ancestral Columbia River. Well-rounded clasts are a mixture of extrabasinal granitic and metamorphic quartzite, and intrabasinal porphyritic basaltic andesite, dacite, and basalt from the western Cascades and Columbia River Basalts. Troutdale terrace gravels grade northward into contemporaneous volcanic pebble and cobble gravel terrace deposits produced along the ancestral Cowlitz River that are dominantly composed of porphyritic andesite gravel and volcanic sand from the western Cascades. Lower terraces along the Cowlitz River were deposited by the late Pleistocene Missoula Floods. All of these unconsolidated to semiconsolidated gravels and sands are prone to landslides, and the Aldercrest-Banyon landslide, the second worst landslide disaster in American history, occurred in the Troutdale Formation gravels. After eruption of the Grays River Volcanics and deposition of the Cowlitz Formation, the forearc underwent a period of transtension in the late-middle Eocene related to magmatic upwelling and reorganization of the subducting Farallon Plate. This event produced a northwest-trending set of oblique slip normal faults, along which Grays River dikes intruded. Starting in the early Miocene the region underwent a transpressional event, reactivating many of the northwest-trending faults, and producing the Columbia Heights Anticline, Hazel Dell Syncline, the Coal Creek Fault, and the Kelso Fault Zone. The paleotopography resulting from this event was stream eroded to a nearly flat plain before emplacement of the Columbia River Basalts, which are nearly horizontal today. Continued offset along the northwest-trending fault set has also offset the Columbia River Basalts. Continued oblique slip post-Miocene broad arching of the Coast Range and downcutting by the Columbia and Cowlitz Rivers has resulted in Pliocene and Pleistocene terraces, and produced an east-west fault set that offsets all earlier structural features.
-
2377. [Article] Contemporary regional forest dynamics in the Pacific Northwest
Recent climatic warming trends and increases in the frequency and extent of wildfires have prompted much concern regarding the potential for rapid change in the structure and function of forested ecosystems ...Citation Citation
- Title:
- Contemporary regional forest dynamics in the Pacific Northwest
- Author:
- Reilly, Matthew J. (Matthew Justin), 1975-
Recent climatic warming trends and increases in the frequency and extent of wildfires have prompted much concern regarding the potential for rapid change in the structure and function of forested ecosystems around the world. Episodes of mortality in wildfires and insect outbreaks associated with drought have affected large areas and altered landscapes, but little is known about the cumulative effects of these disturbances at the regional scales. I used data from two different forest inventories in the Pacific Northwest to develop a framework for tracking regional forest dynamics and examine variation in tree mortality rates among vegetation zones that differ in biophysical setting as well as recent and historical disturbance regimes. In the second chapter I developed an empirically based framework for tracking regional forest dynamics using regional inventory data collected from 2001 to 2010. I characterized the major dimensions of forest structure and developed a classification incorporating multiple attributes of forest structure including biomass, size, and density of live trees, the distribution and abundance of dead wood, and the cover of understory vegetation. A single dimension related to live tree biomass accounted for almost half of the variation in a principal components analysis of structural attributes, but dimensions related to density and size of live trees, dead wood, and understory vegetation accounted for as much additional variation. Snags and biomass of dead and downed wood were related to multiple dimensions while understory vegetation acted independent of other dimensions. Results indicated that structural development is more complex than a monotonic accumulation of live biomass and that some components act independently or emerge at multiple stages of structural development. The hierarchical classification reduced the data into three “groups” based on live tree biomass, followed by eleven "classes" that varied in density and size of live trees, and finally twenty-five structural types that differed further in the abundance of dead wood and cover of understory vegetation. Most structural types were geographically widespread but varied in age of dominant trees by vegetation zone indicating that similar structural conditions developed in environments with different biophysical setting, climate, and disturbance/successional histories. Low live biomass structural types (<25 Mg/ha) differed in live tree density and the abundance of live and dead legacies, demonstrating that the variation in early developmental stages depends on the rate of tree establishment and the nature and severity of recent disturbance. Forests in early developmental stages made up less than 20% of most vegetation zones and diverse types with live or dead legacies associated with wildfires were rare. Moderate live biomass structural types (25-99 Mg/ha) represented multiple mid, mature, and late developmental stages, some of which lack analogs in existing conceptual models of structural development such as lower density woodlands with big trees. These structural types included two that have high densities of snags indicative of recent episodes of mortality; together these made up as much as 10% of some dry vegetation zones. Several high live biomass structural types (100->300 Mg/ha) were identified and substantiated the diversity and relative dominance of mature and later developmental stages, particularly in wet vegetation zones. The relative abundance and make up of structural types varied widely by vegetation zone. Most forests in wet vegetation zones had moderate to high live biomass and were in mid and mature developmental stages, while diverse early developmental stage stages were extremely rare. Dry forests had a far greater range of variation in the relative abundance of structural types which is partially attributable to the greater range of climatic conditions they included, but also to the occurrence of recent episodes of mortality associated with wildfires and insects. In the third chapter I examined variation in tree mortality rates using a different regional inventory that occurred from the mid-1990s to the mid-2000s. I compared the distribution of rates among stands in different vegetation zones and stages of structural developmental. I developed a simple framework based on changes in live tree density and mean tree size and examined trends in structural change associated with disturbances at different levels of mortality across all stages of structural development. Most plots were within the range of "background" mortality rates reported in other studies (<1.0 %/yr) and extremely high "stand-replacing" levels of mortality (>25%/yr) were rare. Approximately 30% of plot mortality rates occurred at intermediate levels (>1%/yr and <25%/yr) as result of insects and fire, highlighting the importance of conceptualizing mortality as a continuum as opposed to just “background” or “stand-replacement” to fully represent dynamics at a regional scale. The distributions of mortality differed among many vegetation zones. Levels of mortality were primarily <2.5%/yr in western hemlock, silver fir, and mountain hemlock vegetation zones where fires were rare and insects and pathogens occurred predominantly at endemic levels. Rates were highest in subalpine forests and higher elevation grand fir and Douglas-fir forests as a result of fire and insects. Mortality rates in ponderosa pine, the hottest driest forest vegetation zone, were surprisingly low, and there was little to no mortality in plots with no evidence of disturbance. Mortality rates varied among developmental stages in all vegetation zones but few consistent patterns emerged. Levels of mortality were often lowest in early developmental stages but varied in later stages where they were lowest in wet vegetation zones and highest in subalpine and dry vegetation zones. Application of a simple framework indicated that multiple trajectories of structural change were common at levels of mortality <2.5%/yr, but structural change at higher levels was predominantly associated with a “thinning” trajectory defined by decreases in density and increases in mean tree size. Results indicated that the rate and magnitude of mortality related change during the study period varies widely across the region. Rapid change has occurred in subalpine, grand fir/white fir, Douglas-fir, and ponderosa pine vegetation zones where disturbances such as insects and fire were widespread. However, these disturbances have potentially restored some aspects of historical structure by reducing overall density and increasing the dominance of bigger trees. In western hemlock, silver fir, and mountain hemlock vegetation zones where higher levels of mortality related to disturbances were rare, wildfires have increased landscape diversity by creating diverse early successional habitats and most change was more subtle but may be manifest oevr longer periods if current trends continue. This examination of short-period mortality rates and associated structural change across a broad geographic provides context for understanding trends from localized studies and potential ecological consequences of mortality, but there is still a great deal of uncertainty as to how the effects of a changing climate and disturbance regimes will manifest themselves over longer time scales. This dissertation is one of the first field based assessments of recent forest dynamics at a regional scale. The results of both chapters, each based on a different dataset, told a similar story. The abundance of structural types in various vegetation zones estimated during the mid-2000s was consistent with the cumulative effects of tree mortality during the preceding decade. It was evident that wildfire effects and recent mortality were small relative to the regional extent of the study and have contributed to structural diversity and restoration of historic structure in stands where fire exclusion and past logging has increased total stand density and decreased the dominance of big trees. However, the rate of change and cumulative effects of recent forest dynamics varied widely by geographic location and vegetation zone and there was greater variability and uncertainty regarding the effects of mortality at smaller landscape scales where individual events like large wildfires have the potential to rapidly alter the landscape structure and composition. Assessing this variability and the scales at which trade-offs (e.g. losses of old-growth and creation of diverse early developmental stages) occur will be an important next step in understanding the cumulative ecological effects of recent wildfires and tree mortality on Pacific Northwest forests.
-
2378. [Article] Distribution and movements of Chinook salmon, Oncorhynchus tshawytscha, returning to the Yukon River basin
Chinook salmon, Oncorhynchus tshawytscha, returning to the Yukon River basin and other large river systems in western Alaska have declined dramatically since the late 1990s. This continuing trend has ...Citation Citation
- Title:
- Distribution and movements of Chinook salmon, Oncorhynchus tshawytscha, returning to the Yukon River basin
- Author:
- Eiler, John H.
Chinook salmon, Oncorhynchus tshawytscha, returning to the Yukon River basin and other large river systems in western Alaska have declined dramatically since the late 1990s. This continuing trend has raised concerns over the future status of the returns, and severely impacted commercial and subsistence fisheries within the drainage. Management is further complicated by the mixed-stock composition of the run, the presence of other temporally similar salmon species, and the need to equitably allocate harvests between the numerous fisheries and user groups scattered throughout the basin. Detailed information is needed on Chinook salmon run characteristics to better understand and manage the returns, and facilitate conservation efforts. However, this goal is exacerbated by the massive size and remote nature of the basin, the large number of highly mobile fish, and the compressed timing of the run. To address these challenges, radio telemetry was used to determine the stock composition and spawning distribution of the returns, and the migratory characteristics of the fish. The migratory patterns exhibited by returning salmon provide a number of insights into the status of the run. Since the Yukon River is essentially free-flowing (i.e., not regulated), this study also presented an opportunity to document the distribution and upriver movements of large returns of wild Chinook salmon under natural conditions. During 2002-2004, returning adult Chinook salmon were captured in the lower Yukon River (approximately 300 km upriver from the river mouth), tagged with radio transmitters, and tracked upriver using remote tracking stations located on important migratory routes and major spawning tributaries. Aerial tracking surveys were used to locate fish in spawning areas and between stations. The fish responded well to the capture and handling procedures, with most (2,790, 98%) resuming upriver movements. Although the fish initially displayed a negative tagging response, with slower migration rates observed immediately after release, the duration of this response was relatively short (several days) and less severe as the fish moved upriver. Independent measures indicated that the swimming speeds and timing of the fish upriver from the tagging area were comparable to untagged fish, suggesting that the tagging methods used were relatively benign. Fish returned to spawning areas throughout the basin, ranging from several hundred to over 3,000 km from the tagging area. Distribution patterns were similar across years, suggesting that the principal components of the run were identified. Most spawning fish were clustered in a number of key tributaries, with smaller numbers of fish located in other spatially isolated areas. The fish typically returned to clear water tributaries that were relatively entrenched, had moderate gradients, and were associated with upland areas. Fish were largely absent in lowland reaches characterized by meandering, low gradient, highly alluvial channels often associated with main river floodplains. There was suggestive evidence of mainstem spawning in reaches of the Upper Yukon. The status of fish remaining in other mainstem areas was less certain, and may represent local spawning activity or fish that died while in-transit to upriver areas. Although Chinook salmon spawned throughout the basin, the run was dominated by two regional components (Tanana and Upper Yukon), which annually comprised over 70% of the return. Substantially fewer fish returned to other areas ranging from 2-9% of the return, although the collective contribution of these stocks was appreciable. Most regional returns consisted of several principal stocks and a number of small, spatially isolated populations. Regional and stock composition estimates were similar across years even though differences in run abundance were reported, suggesting that these abundance differences were not related to regional or stock-specific differences. Run timing was relatively compressed compared to rivers in the southern portion of the range, with most stocks passing through the lower river over a 6-week period, ranging from 16 to 38 d. Run timing was generally earlier for stocks traveling farther upriver, although exceptions were noted. Lower basin stocks were primarily later run fish. Pronounced differences were observed in the migration rates (km/d) exhibited by regional stocks. Substantially slower swimming speeds were observed for fish returning to terminal tributaries in the lower basin ranging from 28-40 km/d compared to 52-62 km/d for upper basin stocks. The migratory patterns (migration rates in sequential reaches) of the fish also showed distinct regional differences. Average migration rates through the lower river were remarkably similar for the different stocks, ranging from 57-62 km/d, with most stocks exhibiting a general decline as the fish moved farther upriver. Tanana River stocks displayed a pronounced reduction in swimming speed after leaving the Yukon River main stem, with migration rates declining to 24 km/d on average as the fish approached their terminal tributaries. Conversely, upper basin stocks exhibited a relatively gradual (but variable) overall decline in migration rate even though these fish were traveling substantially greater distances upriver. Average migration rates for upper basin stocks ranged from 43-61 km/d as the fish approached their terminal tributaries. There was substantial variation in the migratory patterns exhibited by individual fish, although these patterns tended to be similar to the patterns exhibited by the regional stocks, particularly as the fish moved farther upriver from the tagging area. The dominant source of variation among fish reflected the average migration rate, with individual fish traveling slower in the lower basin exhibiting consistently slower migration rates as they moved upriver compared to their faster moving counterparts. This migratory pattern was consistent across stocks, and on average explained 74% of the within-stock variation in migration rate represented by the multivariate data. The second source of variation in migration rate reflected a shift in the relative swimming speeds of the individual fish as they progressed upriver. Although movement rates declined for nearly all of the fish during the migration, differences were observed in the pattern of the decline. Fish with faster migration rates in the lower river exhibited a pronounced decline in swimming speed as they moved upriver, whereas fish moving slower in the lower river displayed a more gradual decline in migration rate. On average, this migratory pattern explained 22% of the within-stock variation in migration rate represented by the multivariate data. Most fish (98%) exhibited continuous upriver movements and strong fidelity to the rivers they entered. However a small number of fish (n = 66) deviated from this pattern. Some of these individuals initially passed their final destination and continued upriver for varying distances before reversing direction, swimming back downstream, and entering their terminal tributary. Although most of these excursions were relatively short (< 30 km), there were several instances where fish traveled hundreds of kilometers out of their way. Thirty-four fish tracked to terminal tributaries subsequently left these rivers, and traveled to other terminal tributaries within the basin (n = 31) or were harvested in upriver fisheries (n = 3). Although most of these incidents involved nearby tributaries, major diversions were also observed, with several fish traveling over 300 km to natal rivers after leaving the initial tributary. Chinook salmon returns to the Yukon River typically consisted of a series of distinct and sizable increases in the number fish entering the river over the course of the run, commonly referred to as pulses. A large number of fish (n = 251) were radio tagged over a 4-day period during a pulse in 2003 to provide information on the progression of the pulse as it moved upriver. The time taken by the pulse to move past subsequent upriver locations increased as the fish moved farther upriver from the tagging area, with the fish passing sites located 580 and 800 km upriver over a span of 14 and 21 d, respectively. Although not surprising considering the extensive variation in migration rates observed among individual fish, this finding does suggest that these pulses do not represent cohesive aggregates of fish moving upriver. Unlike the well established methods used to estimate other life history characteristics, the development of quantitative methods for analyzing and modeling fish movements has lagged noticeably behind, due in part to the complexity associated with movement data and (prior to the advent of telemetry) the difficulty of collecting this type of information on free-ranging individuals. Two fundamentally different analytical approaches, hierarchical linear regression models and multivariate ordination, were used during this study to evaluate factors thought to influence the upriver movements of the fish. In spite of the inherent differences, both methods provided strikingly similar results, indicating that the study findings were not dependent on the approach used, and suggesting that the results were plausible based on the information available and the weight of evidence. Both analytical methods had advantages, and provided complementary information. With hierarchical linear models, it was possible to simultaneously evaluate a wide range of explanatory variables (in our case, both biological and environmental), which provided standardized comparisons and simplified the interpretation of the results. Since both fixed and random effects were incorporated in the models, it was possible to account for sources of variation when insufficient information was available to identify the underlining factors – an important consideration since few field studies provide comprehensive data. With multivariate ordination, separate analyzes were needed to examine the relationships between the migration rates and the biotic and physical variables. In addition to being cumbersome, this limitation made it more difficult to compare the relative influence of the different factors and interactions between factors. However, ordination was very useful as an exploratory tool. Although compartmentalized by stock, across fish comparisons were simple and relatively straightforward. Because the explanatory variables were evaluated separately in relation to the ordination score assigned to the fish, it was possible to examine and compare highly correlated variables. Ordination was also able to identify overall patterns within the data and assess the relative importance. While this can be accomplished within the framework of linear regression using mixture models to determine whether multiple distributions exist within the data, the process is much simpler with ordination. The migratory patterns of the fish were influenced by a wide range of factors, with evidentiary support for complex, multi-faceted relationships. Physical features of the basin demonstrated stronger explanatory power, accounting for over 70% of the observed variation in migration rate compared to 18% for the biological characteristics of the fish. Parameter estimates associated with the steepness of the migratory route and remaining distance the fish had to travel to reach their natal rivers were most strongly correlated with migration rate, with consistent relationships observed across stocks. Migration rates were also noticeably slower in extensively braided reaches of the basin. The weaker relationships between migration rate and biotic factors may reflect stabilizing selection on long-distance migrants. Smaller fish exhibited minimally faster swimming speeds on average than larger individuals. This relationship was stronger in highly braided reaches. Run timing was positively related to migration rate for most stocks. Surprisingly, upper basin stocks traveling farther upriver displayed progressively negative relationships, suggesting that late-run fish were moving slower. Ancillary information suggests that this decline may relate to deteriorating fish condition later in the season.