Estimation of forest harvesting-induced stream temperature changes and bioenergetic consequences for cutthroat trout in a coastal stream in British Columbia, Canada

Data from a paired-catchment study in south coastal British Columbia, Canada, were analyzed to assess the thermal effects of clearcut harvesting with no riparian buffer on a fish-bearing headwater stream. The approach used time series of daily mean water temperatures for East Creek (control) and A Creek (treatment), both before and after harvest. Statistical models were developed to predict (a) what the temperatures would have been in the post-harvest period had harvesting not occurred, and (b) what temperatures would have been in the pre-harvest period had harvesting already occurred. The Wisconsin Bioenergetics Model was used to simulate growth of coastal cutthroat trout (Oncorhynchus clarki clarki) for the first year following fry emergence using the predicted and observed stream temperatures to generate scenarios representing with-harvest and no-harvest thermal regimes. A Monte Carlo approach was used to quantify the effects of uncertainty associated with the regression models on predicted stream temperature and trout growth. Summer daily mean temperatures in the with-harvest scenario were up to $5^{\circ}\hbox{C}$ higher than those for the no-harvest scenario. Harvesting-induced warming reduced growth rates during summer, but increased growth rates during autumn and spring. In the with-harvest scenario, trout were 0.2?C2.0?g (absolute weight) smaller throughout the winter period than in the no-harvest scenario. However, the bioenergetic simulations suggest that trout growth may be more sensitive to potential changes in food supply following harvesting than to direct impacts of stream temperature changes.     

Variability and comparison of hyporheic water temperatures and seepage fluxes in a small Atlantic salmon stream

Ground water discharge is often a significant factor in the quality of fish spawning and rearing habitat and for highly biologically productive streams. In the present study, water temperatures (stream and hyporheic) and seepage fluxes were used to characterize shallow ground water discharge and recharge within thestreambed of Catamaran Brook, a small Atlantic salmon (Salmo salar) stream in central New Brunswick, Canada. Three study sites were instrumented using a total of 10 temperature sensors and 18 seepage meters. Highly variable mean seepage fluxes, ranging from 1.7 x 10(-4) to 2.5 cm3 m(-2) sec(-1), and mean hyporheic water temperatures, ranging from 10.5 degrees to 18.0 degrees C, at depths of 20 to 30 cm in the streambed were dependent on streambed location (left versus right stream bank and site location) and time during the summer sampling season. Temperature data were usefulfor determining if an area of the streambed was under discharge (positive flux), recharge (negative flux), or parallel flow (no flux) conditions and seepage meters were used to directly measure the quantity of water flux. Hyporheic water temperature measurements and specific conductance measurements of the seepage meter sample water, mean values ranging from 68.8 to 157.9 microS/cm, provided additional data for determining flux sources. Three stream banks were consistently under discharge conditions, while the other three stream banks showed reversal from discharge to recharge conditions over the sampling season. Results indicate that the majority of the water collected in the seepage meters was composed of surface water. The data obtained suggests that even though a positive seepage flux is often interpreted as ground water discharge, this discharging water may be of stream water origin that has recently entered the hyporheic zone.The measurement of seepage flux in conjunction with hyporheic water temperature or other indicators of water origin should be considered when attempting to quantify the magnitude of exchange and the source of hyporheic water.     

Estimation of vertical water fluxes from temperature time series by the inverse numerical computer program FLUX‐BOT

The application of heat as a hydrological tracer has become a standard method for quantifying water fluxes between groundwater and surface water. The typical application is to estimate vertical water fluxes in the shallow subsurface beneath streams or lakes. For this purpose, time series of temperatures in the surface water and in the sediment are measured and evaluated by a vertical 1D representation of heat transport by advection and conduction. Several analytical solutions exist to calculate the vertical water flux from the measured temperatures. Although analytical solutions can be easily implemented, they are restricted to specific boundary conditions such as a sinusoidal upper temperature boundary. Numerical solutions offer higher flexibility in the selection of the boundary conditions. This, in turn, reduces the effort of data preprocessing, such as the extraction of the diurnal temperature variation from the raw data. Here, we present software to estimate water fluxes based on temperatures—FLUX‐BOT. FLUX‐BOT is a numerical code written in MATLAB that calculates vertical water fluxes in saturated sediments based on the inversion of measured temperature time series observed at multiple depths. FLUX‐BOT applies a centred Crank–Nicolson implicit finite difference scheme to solve the one‐dimensional heat advection–conduction equation. FLUX‐BOT includes functions for the inverse numerical routines, functions for visualizing the results, and a function for performing uncertainty analysis. We present applications of FLUX‐BOT to synthetic and to real temperature data to demonstrate its performance.     

Differential bedload transport rates in a gravel-bed stream: A grain-size distribution approach

Two methods that were recently proposed for calculating fractional bedload transport rates in gravel-bed streams are examined closely. Both of them employ the Oak Creek bedload data. The Diplas (PD) method is guided by dimensional analysis and, therefore, can be used to predict bedload transport in different gravel rivers. The only requirement for using this method is the knowledge of the subsurface material size distribution of the stream of interest. The expression for the fractional bedload transport obtained by the Shih and Komar (SK) method for Oak Creek cannot be used for other streams. Its use for a given stream requires information that is rarely available. For the Oak Creek case both methods demonstrate similar predictive ability.     

Population dynamics of mayflies in a constant temperature spring stream in West Carpathians

We examined the life history and secondary production of four Ephemeroptera species (Baetis alpinus Pictet, 1843–1845, Baetis rhodani Pictet, 1843–1845, Rhithrogena carpatoalpina Klonowska, Olechowska, Sartoriet & Weichselbaumer, 1987 and Habroleptoides confusa Sartori & Jacob, 1986) in a temperature stable cold spring stream at Prosiek valley (Cho?ské vrchy Mts., West Carpathians, Slovakia). We have found asynchronous bivoltine life cycle for the most abundant species B. alpinus with growth rate positively correlated to photoperiod length. R. carpatoalpina have shown unusual asynchronous univoltine life cycle and B. rhodani have shown uncommonly low abundance in mayfly community. Total secondary production of mayfly community was very low, reaching 1654.8 mg DW m^?2 y^?1. We suggest that the observed asynchrony in growth could be related to the lack of temperature control.     

Estimation of Sq variation by means of multiresolution and principal component analyses

A wavelet based method of deconvoluting the solar quiet variation from the low and mid-latitude H-component records is proposed. The resulting daily variation is non-constant, and its day-to-day variability is quantified by functional principal component scores. The procedure removes the signature of an enhanced ring current by comparing the scores at different stations. The method is fully algorithmic and is implemented in the statistical software R.     

Characterization of volcanic thermal anomalies by means of sub-pixel temperature distribution analysis

The simultaneous solution of the Planck equation (involving the widely used “dual-band” technique) using two shortwave infrared (SWIR) bands allows for an estimate of the fractional area of the hottest part of an active lava flow (f _h) and the background temperature of the cooler crust (T _c). The use of a high spectral and spatial resolution imaging spectrometer with a wide dynamic range of 15 bits (DAIS 7915) in the wavelength range from 0.501 to 12.67 μm resulted in the identification of crustal temperature and fractional areas for an intra-crater hot spot at Mount Etna, Italy. This study indicates the existence of a relationship between these T _c and f _h extracted from DAIS and Landsat TM data. When the dual band equation system is performed on a lava flow, a logarithmic distribution is obtained from a plot of the fractional area of the hottest temperature vs. the temperature of the cooler crust. An entirely different distribution is obtained over active degassing vents, where increases in T _c occur without any increase in f _h. This result indicates that we can use scatter plots of T _c vs. fh to discriminate between different types of volcanic activity, in this case between degassing vents and lava flows, using satellite thermal data.     

The effect of a vertical temperature gradient in the lithosphere on seismic and tectonic waves

A linear analysis of the stability of the lithosphere considered as a viscoelastic layer with an equilibrium vertical gradient of temperature is carried out. The problem is solved with a complete system of linearized equations of a continuous medium represented in the dimensionless form and containing a set of dimensionless parameters that determine thermomechanical properties of the lithosphere. As a result of the stability analysis, decrements are found that give the time dependence of perturbations and correspond to high-frequency seismic waves and low-frequency tectonic waves. The frequency and velocity of seismic waves are determined by the elasticity and inertial properties of the lithosphere, and their attenuation, by viscous properties of the lithosphere. The temperature gradient existing in the lithosphere influences seismic waves very weakly. On the contrary, the pattern of tectonic waves is controlled by the temperature gradient and viscous properties, while the effect of elastic and inertial properties on these waves is negligibly small. The stability of a viscoelastic lithosphere is examined using such rheological models as the Maxwell, standard linear, and Andrade media (the frequency of tectonic waves is zero in the Maxwell medium).     

Estimation of time dependent virus inactivation rates by geostatistical and resampling techniques: application to virus transport in porous media

A methodology is developed for estimating temporally variable virus inactivation rate coefficients from experimental virus inactivation data. The methodology consists of a technique for slope estimation of normalized virus inactivation data in conjunction with a resampling parameter estimation procedure. The slope estimation technique is based on a relatively flexible geostatistical method known as universal kriging. Drift coefficients are obtained by nonlinear fitting of bootstrap samples and the corresponding confidence intervals are obtained by bootstrap percentiles. The proposed methodology yields more accurate time dependent virus inactivation rate coefficients than those estimated by fitting virus inactivation data to a first-order inactivation model. The methodology is successfully applied to a set of poliovirus batch inactivation data. Furthermore, the importance of accurate inactivation rate coefficient determination on virus transport in water saturated porous media is demonstrated with model simulations.     

Estimation of groundwater velocity in localized fracture zones from well temperature profiles

In fractured rocks, well temperature logs often exhibit abrupt temperature changes over a small distance, typically surging or dropping several degrees over several tens of meters. The abrupt temperature changes usually occur in localized fracture zones or small faults. Away from these zones, temperature quickly restores itself to the background linear conductive profile. A theoretical model is presented herein to interpret the pattern of such abrupt temperature changes. I illustrate that this type of temperature profile originates from a different physical process than the physical processes found in some other typical temperature patterns in large scales. I also demonstrate that, using the field data, the theoretical model can be effectively used to estimate fluid flow velocity in fracture zones or local faults. Using the temperature profiles from two boreholes, fluid flow velocity in fracture zones is estimated to be as high as 1.1×10⁻⁶ m/s. Although fracture flow is a highly localized feature, it can reach three orders of magnitude higher than fluid flow in the hosting rock matrices.     

Estimation of vertical water and heat fluxes in the semi-confined aquifers in tokyo metropolitan area,japan

Groundwater circulation is known to be one of the agents responsible for the redistribution of geothermal energy by acting as a source or sink in the course of its movement through porous media. Heat transport in groundwater systems is considered to be a coupled process and the theory based on this was used to analyse temperature profiles of 30 thermally stable observation wells in a deep, semi-confined aquifer system in the Tokyo Metropolitan area. Vertical water fluxes in the semi-confined aquifers and the associated upward heat fluxes were estimated from a heat flux equation that describes convection and conduction processes of heat transport in one dimension. The vertical downward water fluxes in Shitamachi lowland, Musashino and Tachikawa terraces were 0.69.26.91 × 10^?9, 1.46-70.92 × 10^?9 and 2.61.2204 × 10^?9 m/s, respectively. A vertical upward water flux of 1.80-33.60 × 10^?9 m/s was estimated in Shitamachi lowland. The water flux generally decreased with increasing depth for observation wells which intercepted more than one semi-confining layer. The estimated upward heat fluxes for Shitamachi lowland, Musashino and Tachikawa terraces were 0.32-1.12, 0.49-1.21 and 1.00-11.62 W/m², respectively. The heat flux was highest in Tachikawa terrace where a major fault, the Tachikawa fault, is located. Generally, the estimated heat flux was higher in the semi-confining layers than in the aquifers. Areas with heat sources and sinks as well as groundwater flow patterns in the semi-confined aquifers were revealed by heat flux and temperature distributions in the study area.     

Development of an electronic seepage chamber for extended use in a river

Seepage chambers have been used to characterize the flux of water across the water-sediment interface in a variety of settings. In this work, an electronic seepage chamber was developed specifically for long-term use in a large river where hydraulic gradient reversals occur frequently with river-stage variations. A bidirectional electronic flowmeter coupled with a seepage chamber was used to measure temporal changes in the magnitude and direction of water flux across the water-sediment interface over an 8-week period. The specific discharge measured from the seepage chamber compared favorably with measurements of vertical hydraulic gradient and previous specific discharge calculations. This, as well as other supporting data, demonstrates the effectiveness of the electronic seepage chamber to accurately quantify water flux in two directions over a multimonth period in this setting. The ability to conduct multimonth measurements of water flux at a subhourly frequency in a river system is a critical capability for a seepage chamber in a system where hydraulic gradients change on a daily and seasonal basis.     

Simulating estuarine nitrous oxide production by means of a dynamic model

We describe a dynamic model developed from a commercially available modeling package (ECoS-III) to simulate estuarine dissolved inorganic nitrogen (DIN) dynamics, and consequent N(2)O production and atmospheric flux on the timescale of tidal cycles. Simulated model state variables were NH(4)(+), NO(3)(-) and N(2)O concentrations, and salinity. Model outputs were evaluated through comparison with summer field data for the Tyne estuary, UK. The model adequately reproduced the observed axial profiles of NH(4)(+), NO(3)(-) and N(2)O concentrations. Nitrification was shown to be the dominant N(2)O source and estimates of the ratios nitrification to DIN load and N(2)O emission to DIN load are considerably lower than the corresponding values adopted in global scale models of estuarine N(2)O emissions based on DIN transformations. Hence our results are consistent with the requirement imposed by atmospheric N(2)O growth rate constraints that the amount of atmospheric N(2)O arising from agriculturally related sources, including estuarine transformations of N, be revised downward.     

Implications of bioturbation induced by Procambarus clarkii on seepage processes in channel levees

River levees are subject to bioturbation by various animals which can actively excavate into earthen structures producing an internal erosion that,during the passage of a flood,can grow in time making the levee unstable.This phenomenon can lead to river levee breaching and,as a consequence,collapse,even for relatively minor flood events.A well-known animal burrower is represented by the North American crayfish Procambarus clarkii(P.clarkii),an invasive species in Europe,mainly introduced for commercial purposes,causing a decline in biodiversity and profound habitat changes.The physical damages caused by P.clarkii on levees and banks,such as in rice fields,irrigation ditches,and small channels,have not been fully studied and behavioral components underlying this impact are mostly occasional.To under-stand the impact of burrowing activity on the seepage process,a field survey was done in a drainage channel in Tuscany,Italy,to evaluate the density and geometry of the internal burrows that were excavated by the crayfish.Based on these observations and some previous laboratory experiments,three dimensional(3D)numerical simulations of the seepage processes were done inside burrowed levees.Numerical results allowed the increase in the hydraulic vulnerability of levees to the process of internal seepage to be disclosed.In particular,for a given river water level,the reduction of the time scale for the phreatic line to reach the levee field side appears to be a function of a quantity here defined as the burrow hydraulic gradient.This quantity is here defined as the ratio between the hydraulic head inside the burrow and the horizontal distance from its end to the field side of the levee.Moreover,a comparison between the 3D with the analogous more common two dimensional(2D)numerical simulations illus-trated the schematization which is better suited for describing the seepage processes when animal burrows,not only by crayfish,are present.     

VHF Radar observations of a Kelvin-Helmholtz instability in a subtropical jet stream

Abstract

In connection with the large dish at the Arccibo Observatory (Puerto Rico), the mobile SOUSY VHF Radar has been used to carry out velocity measurements in the troposphere with height and time resolutions of a few hundreds of meters and some tens of seconds, respectively. During the passage of a strong subtropical jet stream, vertical velocity oscillations have been observed at heights of maximum wind shear. The height and time variations of the velocity and the background wind profile deduced from radiosonde data indicate that a Kelvin-Helmholtz instability (KHI), with a period of 340s, a wavelength of 12 km and a velocity amplitude of about 1 ms^?1 was the source of the oscillations. Model calculations show good agreement with the observed quantities.     

Comparison of industrial and high-resolution thermal data in a sedimentary basin

Temperature, temperature gradient and heat flow, derived from four wells in the Western Canada Sedimentary Basin have been compared with similar quantities derived from maps constructed from industrial data. Individual industrial temperature data may differ from the high-resolution temperature log by up to 30 K, but linear regression of the collected data agrees within 10 K at all points observed. Some evidence suggests that measured conductivities, using drill-cuttings, are biased toward average values. Derived heat flows show agreement of heat flow within 10% within the Mesozoic section. In the Paleozoic section differences are greater, and more varied, with insufficient data for general conclusions. Both styles of measurement provide opportunities for interpretation, each contributing to thermal analysis of sedimentary basins.