Sampling Fauna in Stream Sediments as well as Groundwater Using One Net Sampler

Groundwater bores act as traps. Net samplers are regularly used for sampling this type of trap for fauna. To enable direct comparisons of faunal communities in groundwater bores and stream sediments, stream sediment tubes were built similar to groundwater bores and were sampled with net samplers for fauna. These stream sediment tubes consisted of a tube anchored in the stream sediment, also called interstitial space. To test the efficacy of this trap method in stream sediments, it was compared to another type of trap, Hahn's trap. Faunal communities sampled by a net in the stream sediment tubes did not differ hugely from fauna in Hahn's trap samples. Physical and chemical factors of sampled water in both the stream sediment tubes, the surrounding interstitial sediments and the second type of traps, Hahn's traps, showed that water in both the tubes and Hahn's traps was closely related to interstitial water. The net sampler is inexpensive and easy to handle. It is suggested that sampling stream tubes with nets may be an appropriate method for long‐term monitoring studies.     

Bayesian updating revisited

Bayesian updating methods provide an alternate philosophy to the characterization of the input variables of a stochastic mathematical model. Here, a priori values of statistical parameters are assumed on subjective grounds or by analysis of a data base from a geologically similar area. As measurements become available during site investigations, updated estimates of parameters characterizing spatial variability are generated. However, in solving the traditional updating equations, an updated covariance matrix may be generated that is not positive-definite, particularly when observed data errors are small. In addition, measurements may indicate that initial estimates of the statistical parameters are poor. The traditional procedure does not have a facility to revise the parameter estimates before the update is carried out. alternatively, Bayesian updating can be viewed as a linear inverse problem that minimizes a weighted combination of solution simplicity and data misfit. Depending on the weight given to the a priori information, a different solution is generated. A Bayesian updating procedure for log-conductivity interpolation that uses a singular value decomposition (SVD) is presented. An efficient and stable algorithm is outlined that computes the updated log-conductivity field and the a posteriori covariance of the estimated values (estimation errors). In addition, an information density matrix is constructed that indicates how well predicted data match observations. Analysis of this matrix indicates the relative importance of the observed data. The SVD updating procedure is used to interpolate the log-conductivity fields of a series of hypothetical aquifers to demonstrate pitfalls and possibilities of the method.     

Geostatistics and Bayesian updating for transmissivity estimation in a multiaquifer system in Manitoba,Canada

In ground water flow and transport modeling, the heterogeneous nature of porous media has a considerable effect on the resulting flow and solute transport. Some method of generating the heterogeneous field from a limited dataset of uncertain measurements is required. Bayesian updating is one method that interpolates from an uncertain dataset using the statistics of the underlying probability distribution function. In this paper, Bayesian updating was used to determine the heterogeneous natural log transmissivity field for a carbonate and a sandstone aquifer in southern Manitoba. It was determined that the transmissivity in m²/sec followed a natural log normal distribution for both aquifers with a mean of -7.2 and - 8.0 for the carbonate and sandstone aquifers, respectively. The variograms were calculated using an estimator developed by Li and Lake (1994). Fractal nature was not evident in the variogram from either aquifer. The Bayesian updating heterogeneous field provided good results even in cases where little data was available. A large transmissivity zone in the sandstone aquifer was created by the Bayesian procedure, which is not a reflection of any deterministic consideration, but is a natural outcome of updating a prior probability distribution function with observations. The statistical model returns a result that is very reasonable; that is homogeneous in regions where little or no information is available to alter an initial state. No long range correlation trends or fractal behavior of the log-transmissivity field was observed in either aquifer over a distance of about 300 km.     

Variations in sediment sources and yields in the Finger Lakes and Catskills regions of New York

The proportional contributions of stream bank and surface sources to fine sediment loads in watersheds in New York State were quantified with uncertainty analysis. Eroding streamside glacial drift, including glaciolacustrine deposits, were examined to help explain variations in the proportional contributions made by bank erosion. Sediment sources were quantified by comparing concentrations of the bomb‐derived radionuclide ¹³⁷Cs in fluvial sediment with sediment from potential source areas such as agricultural soils, forest soils and stream banks. To compare sediment sources in streams containing abundant deposits of fine‐grained glacial drift with watersheds that lacked moderate or extensive streamside deposits, samples were taken from 15 watersheds in the region. The mean contribution of bank erosion to sediment loads in the six streams with glaciolacustrine deposits was 60% (range 46–76%). The proportional contribution of bank erosion was also important in one stream lacking glaciolacustrine deposits (57%) but was less important in the remainder, with contributions ranging from 0 to 46%. Data from this study on the varying contributions of bank erosion and data from past studies of sediment yield in 15 watersheds of New York State suggest that eroding streamside glacial deposits dominate sediment yield in many watersheds. In other watersheds, past impacts to streams, such as channelization, have also resulted in high levels of bank erosion. Copyright © 2007 John Wiley & Sons, Ltd.     

Clay aquitards as archives of holocene paleoclimate: delta18O and thermal profiling

Paleoclimatic conditions in the Holocene were reconstructed from a detailed vertical profile of pore water delta(18)O and a series of downhole thermal profiles at a thick, hydrogeologically simple, aquitard research site in the Northern Great Plains of Saskatchewan. Reconstructions were obtained using the theory of one-dimensional diffusive transport and an empirical Bayesian inversion technique. Inversion of the delta(18)O profile shows that input signal consists of a sudden increase of +6 per thousand (corresponding to a temperature increase of about 18 degrees C) at about 12,000 years before present (BP), after which no measurable change in delta(18)O is apparent to present day. This research shows, at this location, that there is no evidence of large amplitude temperature shifts in the Holocene and supports the commonly used assumption in ground water studies of constant atmospheric delta(18)O values throughout the Holocene. Inversion of the temperature profiles suggests the ground surface temperature increased primarily in the last half of the 20th century, with a peak temperature (about 3 degrees C) circa 1990. For both profiles, the ability to resolve historical variability decays rapidly with time. For the temperature profiles, the decay in resolution precluded the development of reliable estimates of climatic conditions prior to about 1950 and, in the case of the delta(18)O profile, it prevented the precise definition of climate changes (e.g., Hypsithermal and Little Ice Age) in the Holocene.     

Observed thermal pollution and post-development simulations of low-temperature geothermal systems in Winnipeg, Canada

Two issues arise in the long-term use of groundwater for thermal purposes: (1) the sustainability of an individual system; and (2) the effect of neighbouring systems on each other. Both of these effects are observed in an area of the Carbonate Rock Aquifer beneath Winnipeg in Manitoba, Canada, where groundwater has been exploited in thermal applications since 1965. In this area, there are four systems that utilize groundwater for cooling purposes that are closely spaced. The current temperatures observed in this area of the Carbonate Rock Aquifer and the results of the numerical modeling conducted in this study confirm that in each system, temperatures at the production well have risen as a result of breakthrough of injected water. The results of numerical modeling also indicate that interference effects are present in three of the four systems examined in this study. The influence of these systems on each other implies that these systems have a spacing that is smaller than the optimum spacing for such systems, and indicates that there is a limit to the density of development that can occur in a given aquifer.Contaminación termal y simulaciones post-desarrollo de sistemas geotermales de baja temperatura en Winnipeg, Canadá     

Transmissivity estimation for highly heterogeneous aquifers: comparison of three methods applied to the Edwards Aquifer,Texas, USA

Obtaining reliable hydrological input parameters is a key challenge in groundwater modeling. Although many quantitative characterization techniques exist, experience applying these techniques to highly heterogeneous real-world aquifers is limited. Three geostatistical characterization techniques are applied to the Edwards Aquifer, a limestone aquifer in south-central Texas, USA, for the purposes of quantifying the performance in an 88,000-cell groundwater model. The first method is a simple kriging of existing hydraulic conductivity data developed primarily from single-well tests. The second method involves numerical upscaling to the grid-block scale, followed by cokriging the grid-block conductivity. In the third method, the results of the second method are used to establish the prior distribution for a Bayesian updating calculation. Results of kriging alone are biased towards low values and fail to reproduce hydraulic heads or spring flows. The upscaling/cokriging approach removes most of the systematic bias. The Bayesian update reduced the mean residual by more than a factor of 10, to 6 m, approximately 2.5% of the total head variation in the aquifer. This agreement demonstrates the utility of automatic calibration techniques based on formal statistical approaches and lends further support for using the Bayesian updating approach for highly heterogeneous aquifers.     

The effects of climatic variability on estimates of recharge from temperature profiles

Using heat as a tracer allows for estimation of ground water recharge rates based on subsurface temperature measurements. While possible in theory, it may be difficult in practice to discriminate the effects of climate from the effects of ground water advection. This study uses synthetic simulations to determine the influence of variability of ground surface temperature (GST) on the ability to estimate vertical specific discharge from temperature profiles. Results suggest that in cases where temperature measurements are sufficiently deep and specific discharge is sufficiently high, estimates of specific discharges will be reasonably accurate. Increasing the number of times temperatures are measured, or producing models that incorporate variations in GST, will increase the reliability of any studies using temperatures to estimate specific discharge. Furthermore, inversions of temperature measurements should be combined with other methods of estimating recharge rates to improve the reliability of recharge estimates.     

Improved Recharge Estimation from Portable,Low‐Cost Weather Stations

Groundwater recharge estimation is a critical quantity for sustainable groundwater management. The feasibility and robustness of recharge estimation was evaluated using physical‐based modeling procedures, and data from a low‐cost weather station with remote sensor techniques in Southern Abbotsford, British Columbia, Canada. Recharge was determined using the Richards‐based vadose zone hydrological model, HYDRUS‐1D. The required meteorological data were recorded with a HOBO^TM weather station for a short observation period (about 1 year) and an existing weather station (Abbotsford A) for long‐term study purpose (27 years). Undisturbed soil cores were taken at two locations in the vicinity of the HOBO^TM weather station. The derived soil hydraulic parameters were used to characterize the soil in the numerical model. Model performance was evaluated using observed soil moisture and soil temperature data obtained from subsurface remote sensors. A rigorous sensitivity analysis was used to test the robustness of the model. Recharge during the short observation period was estimated at 863 and 816 mm. The mean annual recharge was estimated at 848 and 859 mm/year based on a time series of 27 years. The relative ratio of annual recharge‐precipitation varied from 43% to 69%. From a monthly recharge perspective, the majority (80%) of recharge due to precipitation occurred during the hydrologic winter period. The comparison of the recharge estimates with other studies indicates a good agreement. Furthermore, this method is able to predict transient recharge estimates, and can provide a reasonable tool for estimates on nutrient leaching that is often controlled by strong precipitation events and rapid infiltration of water and nitrate into the soil.