Hydrological response to timber harvest in northern Idaho: implications for channel scour and persistence of salmonids

The potential for forest harvest to increase snowmelt rates in maritime snow climates is well recognized. However, questions still exist about the magnitude of peak flow increases in basins larger than 10 km² and the geomorphic and biological consequences of these changes. In this study, we used observations from two nearly adjacent small basins (13 and 30 km²) in the Coeur d'Alene River basin, one with recent, relatively extensive, timber harvest, and the other with little disturbance in the last 50 years to explore changes in peak flows due to timber harvest and their potential effects on fish. Peak discharge was computed for a specific rain‐on‐snow event using a series of physical models that linked predicted values of snowmelt input to a runoff‐routing model. Predictions indicate that timber harvest caused a 25% increase in the peak flow of the modelled event and increased the frequency of events of this magnitude from a 9‐year recurrence interval to a 3·6‐year event. These changes in hydrologic regime, with larger discharges at shorter recurrence intervals, are predicted to increase the depth and frequency of streambed scour, causing up to 15% added mortality of bull trout (Salvelinus confluentus) embryos. Mortality from increased scour, although not catastrophic, may have contributed to the extirpation of this species from the Coeur d'Alene basin, given the widespread timber harvest that occurred in this region. Copyright © 2008 John Wiley & Sons, Ltd.     

On the role of spatial resolution on snow estimates using a process‐based snow model across a range of climatology and elevation

Hydrological processes in mountainous settings depend on snow distribution, whose prediction accuracy is a function of model spatial scale. Although model accuracy is expected to improve with finer spatial resolution, an increase in resolution comes with modelling costs related to increased computational time and greater input data and parameter information. This computational and data collection expense is still a limiting factor for many large watersheds. Thus, this work's main objective is to question which physical processes lead to loss in model accuracy with regard to input spatial resolution under different climatic conditions and elevation ranges. To address this objective, a spatially distributed snow model, iSnobal, was run with inputs distributed at 50‐m—our benchmark for comparison—and 100‐m resolutions and with aggregated (averaged from the fine to the large resolution) inputs from the 50‐m model to 100‐, 250‐, 500‐, and 750‐m resolution for wet, average, and dry years over the Upper Boise River Basin (6,963 km²), which spans four elevation bands: rain dominated, rain–snow transition, and snow dominated below treeline and above treeline. Residuals, defined as differences between values quantified with high resolution (>50 m) models minus the benchmark model (50 m), of simulated snow‐covered area (SCA) and snow water equivalent (SWE) were generally slight in the aggregated scenarios. This was due to transferring the effects of topography on meteorological variables from the 50‐m model to the coarser scales through aggregation. Residuals in SCA and SWE in the distributed 100‐m simulation were greater than those of the aggregated 750 m. Topographic features such as slope and aspect were simplified, and their gradient was reduced due to coarsening the topography from the 50‐ to 100‐m resolution. Therefore, solar radiation was overestimated, and snow drifting was modified and caused substantial SCA and SWE underestimation in the distributed 100‐m model relative to the 50‐m model. Large residuals were observed in the wet year and at the highest elevation band when and where snow mass was large. These results support that model accuracy is substantially reduced with model scales coarser than 50 m.     

Development of an inexpensive automated streamflow monitoring system

Estimating streamflow is time and labour intensive due to the necessity of developing a rating curve. The development of a rating curve involves acquiring in-field measurements of streamflow across a wide range of flow levels, which can be costly and impractical in remote regions with limited seasonal access. Here we showcase an automated system which estimates streamflow multiple times each day, greatly facilitating the development of rating curves for remote or seasonally inaccessible sites. The system uses an emerging technique referred to as particle image velocimetry (PIV) to track the movement of objects and flow structure features on the mobile water surface to generate velocity vector grids. Velocity grids were used to calculate streamflow and facilitate the development of a rating curve. This work contributes to our field-based knowledge for using automated PIV systems to estimate streamflow in small streams (<5 m wide) using automatically distributed particles.     

Formation of cordierite-bearing lavas during anatexis in the lower crust beneath Lipari Island (Aeolian arc,Italy)

Cordierite-bearing lavas (CBL;~105 ka) erupted from the Mt. S. Angelo volcano at Lipari (Aeolian arc, Italy) are high-K andesites, displaying a range in the geochemical and isotopic compositions that reflect heterogeneity in the source and/or processes. CBL consist of megacrysts of Ca-plagioclase and clinopyroxene, euhedral crystals of cordierite and garnet, microphenocrysts of orthopyroxene and plagioclase, set in a heterogeneous rhyodacitic-rhyolitic groundmass containing abundant metamorphic and gabbroic xenoliths. New petrographic, chemical and isotopic data indicate formation of CBL by mixing of basaltic-andesitic magmas and high-K peraluminous rhyolitic magmas of anatectic origin and characterize partial melting processes in the lower continental crust of Lipari. Crustal anatectic melts generated through two main dehydration-melting peritectic reactions of metasedimentary rocks: (1) Biotite + Aluminosilicate + Quartz + Albite = Garnet + Cordierite + K-feldspar + Melt; (2) Biotite + Garnet + Quartz = Orthopyroxene + Cordierite + K-feldspar + Melt. Their position into the petrogenetic grid suggests that heating and consequent melting of metasedimentary rocks occurred at temperatures of 725 < T < 900°C and pressures of 0.4–0.45 GPa. Anatexis in the lower crust of Lipari was induced by protracted emplacement of basic magmas in the lower crust (~130 Ky). Crustal melting of the lower crust at 105 ka affected the volcano evolution, impeding frequent mafic-magma eruptions, and promoting magma stagnation and fractional crystallization processes.     

Far-field tectonics associated with a large impact basin: applications to Caloris on Mercury and Imbrium on the Moon

Lithospheric readjustments after the formation of a large impact crater have been computed. The models predict tectonic perturbations over a major portion of the surface of the planet. The weight of the ejecta and the mechanical perturbations in the crater area give rise to membrane stresses. If the planetary lithosphere is in compression, the direction of maximum compression after the impact is perpendicular to a meridian (in a referential where the basin would be at the North pole). On Mercury the preferential orientation of compressive scarps radial to Caloris is explained in this manner. For a planet in extension grabens perpendicular and parallel to a meridian can appear successively as the mechanical state of the crater area is modified. On the Moon this result is consistent with the large-scale tectonics associated with Imbrium.     

Stability of the oceanic lithosphere with variable viscosity: an initial-value approach

We have studied the problem concerning the onset of convective instabilities below the oceanic lithosphere. A system of linear partial differential equations, in which the background temperature field is time-dependent, is integrated in time to monitor the evolution of incipient disturbances. Two types of rheologies have been examined. One depends strongly on temperature. The other involves a viscosity which is both temperature- and pressure-dependent. The results from this initial-value approach, in which the viscosity profiles migrate downward with time, reveal the importance of considering temperature- and pressure-dependent rheology in issues regarding the development of local instabilities in upper mantle convection. For temperature-dependent viscosity, viscosities of 0(10²⁰P) are required to produce instabilities with growth-rates of 0(.1/Ma). In contrast, these same growth rates can be attained for a temperature- and pressure-dependent viscosity profile with a mean value close to 0(10²⁰P) in the upper mantle, owing to the presence of a low viscosity zone, 0(10²⁰P), existing right below the lithosphere. Unlike the results of temperature-dependent viscosity, whose growth-rates increase with time, the amplification of disturbances in a fluid medium with temperature- and pressure-dependent rheology reaches a maximum at an early age, < 50 Ma, and decreases thereafter with time. This suggests the potential importance played by initial disturbances in the evolution of the oceanic lithosphere.