HYDROLOGIC RESPONSE TO PUMPING AND CONTAMINANT ADVECTION IN A FRACTURED ROCK ENVIRONMENT1

ABSTRACT: Ground water flow and supply at the Whiteshell Research Area (WRA) in southeastern Manitoba and the advective movement of contaminants from a hypothetical nuclear fuel waste disposal vault to the adjacent biosphere and a nearby ground water supply well are simulated using finite-element modeling and numerical particle-tracking technique. The hypothetical vault is located at a depth of 500 m, below the water table, in low-permeability plutonic rock of the Canadian Shield. The rock mass is intersected by high-permeability fracture zones (aquifers), which also act as conduits for vault contaminants to migrate to the ground surface. The ground water resource is, therefore, limited in quantity and quality and should be explored with care. A 30 m deep well, which pumps water at a rate of 120 m³/yr from a low-dipping fracture zone, LD1, reduces natural discharge from the system to augment natural recharge. At this pumping rate, a 100 m or 200 m deep well neither reduces discharge nor induces recharge into the system. Thus, at the WRA, a 30 m deep domestic water supply well best meets the water requirements of a one-person household at the rate of 120 m³/yr. The 100 m and 200 m wells best meet the requirements of a family of six and a family of six with light irrigation, respectively, without capturing contaminants’pathways from the vault to the ground surface. By virtue of the proximity of the 200 m well intake to the hypothetical vault, this well performs best as a purge well at pumping rates of 0,000 m³/yr and greater. Finite-element modeling is useful in evaluating the water supply potential of a fractured rock environment in which a nuclear waste disposal vault is proposed to be sited.     

HEAD LOSS AT MANHOLES IN SURCHARGED SEWER SYSTEMS1

ABSTRACT: A laboratory investigation was conducted to determine the head losses at sewer pipe junctions (manholes) under surcharged conditions. A physical model of a manhole/pipeline system was constructed for head-loss measurements. Head-loss coefficients were determined for a variety of outlet-flow Reynolds numbers, surcharge levels, pipe sizes, flow configurations, and inlet-flow rates. Empirical formulas were also developed to estimate head-loss coefficients. The results indicate that head loss is insensitive to the amount of surcharge, but depends heavily on the flow configuration, relative flow rate, and the change of pipe diameter within the pipelines.     

PREDICTION OF MEAN ANNUAL FLOWS IN NORTH AND CENTRAL FLORIDA1

ABSTRACT: A frequency analysis approach for the prediction of flow characteristics at ungaged locations is applied to a region of high annual precipitation and low topography in north and central Florida. Stationary time series of annual flows are fitted with the lognormal distribution and estimated parameters of the distribution are fitted by third order trend surfaces. These explain 65 and 74 percent of the observed variances in the mean and standard deviation, respectively. Predictions of parameters are then made for several locations previously unused in the study and they are used to estimate the return periods of various flows from the lognormal distribution. Application of the Kolmogorov-Smirnov goodness-of-fit test suggests that only one of the five test stations can be considered significantly different from the observed data, confirming the applicability of this technique.     

SPATIAL MODELLING WITH GEOGRAPHIC INFORMATION SYSTEMS FOR DETERMINATION OF WATER RESOURCES VULNERABILITY APPLICATION TO AN AREA IN MASSIF CENTRAL (FRANCE)1

ABSTRACT: A method for water resources protection based on spatial variability of vulnerability is proposed. The vulnerability of a water resource is defined as the risk that the resource will become contaminated if a pollutant is placed on the surface at one point as compared to another. A spatial modelling method is defined in this paper to estimate a travel time between any point of a catchment and a resource (river or well). This method is based on spatial analysis tools integrated in Geographical Information Systems (GIS). The method is illustrated by an application to an area of Massif Central (France) where three different types of flow appear: surface flow, shallow subsurface flow, and permanent ground water flow (baseflow). The proposed method gives results similar to classical methods of estimation of travel time. The contribution of GIS is to improve the mapping of vulnerability by taking the spatial variability of physical phenomena into account.     

Mass conservation and atmospheric dynamics in the Regional Atmospheric Modeling System (RAMS)

This paper examines the spatio-temporal patterns of atmospheric carbon dioxide transport predicted by the Regional Atmospheric Modeling System (RAMS). Forty-eight hour simulations over northern New England incorporating a simple representation of the diurnal summertime surface carbon dioxide forcing arising from biological activity indicate that, in its native formulation, RAMS exhibits a significant degree of mass non-conservation. Domain-wide rates of non-physical mass gain and mass loss are as large as three percent per day which translates into approximately eleven parts per million per day for carbon dioxide — enough to rapidly dilute the signature of carbon dioxide fluxes arising from biological activity. Analysis shows that this is due to the approximation used by RAMS to compute the Exner function. Substitution of the exact, physically complete equation improves mass conservation by two orders of magnitude. In addition to greatly improving mass conservation, use of the complete Exner function equation has a substantial impact on the spatial pattern of carbon dioxide predicted by the model, yielding predictions differing from a conventional RAMS simulation by as much as forty parts per million. Such differences have important implications both for comparisons of modeled atmospheric carbon dioxide concentrations to observations and for carbon dioxide inversion studies, which use estimates of atmospheric transport of carbon dioxide in conjunction with measurements of atmospheric carbon dioxide concentrations to infer the spatio-temporal distribution of surface carbon dioxide fluxes. Furthermore, use of the complete Exner function equation affects the vertical velocity and water mixing ratio fields, causing significant changes in accumulated precipitation over the region.     

High Spatial and Temporal Variability of Dry Deposition in a Coastal Region

A real meteorological situation characterized by strong spatial and temporal variability of the meteorological fields in a coastal region of eastern Denmark is examined in view of the transport of a passive tracer and dry deposition. Model simulations using a full mesoscale NWP model (COAMPS^TM) at different horizontal resolutions are performed. A realistic simulation showed that the differences in the amount of dry deposited matter can reach one order of magnitude and larger, during the period of one afternoon, depending on the model horizontal resolution. In addition, the results of an idealized experiment with straight coastline indicate that the horizontal model resolution alone is responsible for most of the differences. This study confirms the importance of high spatial and temporal resolution modelling for environmental applications.     

An Approach for the Aggregation of Aerodynamic Surface Parameters in Calculating the Turbulent Fluxes over Heterogeneous Surfaces in Atmospheric Models

Experimental evidence indicates that there is a significant departure of the wind profile above the underlying surface consisting of patches of solid and liquid parts, and plant communities with different morphological from that predicted by the logarithmic relationship, which gives the values larger than those observed. This situation can seriously affect the transfer of momentum, heat and water vapor from the surface fluxes into the atmosphere.The object of this paper is to generalize the calculation of the exchange of momentum between the atmosphere and a very heterogeneous surface, find a general equation for the wind speed profile in a roughness sublayer under neutral conditions, and, then, derive aggregated roughness length and displacement height over the grid cell. The suggested expression for the wind profile is compared with some earlier approaches, using a common parameterization of aerodynamic parameters over the grid cell, and the observations obtained at an experimental site in Philadelphia, PA.     

GIS and geostatistics: Essential partners for spatial analysis

Initially, geographical information systems (GIS) concentrated on two issues: automated map making, and facilitating the comparison of data on thematic maps. The first required high quality graphics, vector data models and powerful data bases, the second is based on grid cells that can be manipulated by suites of mathematical operators collectively termed map algebra. Both kinds of GIS are widely available and are taught in many universities and technical colleges. After more than 20 years of development, most standard GIS provide both kinds of functionality and good quality graphic display, but until recently they have not included the methods of statistics and geostatistics as tools for spatial analysis. Recently, standard statistical packages have been linked to GIS for both exploratory data analysis and statistical analysis and hypothesis testing. Standard statistical packages include methods for the analysis of random samples of cases or objects that are not necessarily co-located in space—if the results of statistical analysis display a spatial pattern then that is because the underlying data also share that pattern. Geostatistics addresses the need to make predictions of sampled attributes (i.e., maps) at unsampled locations from sparse, often expensive data. To make up for lack of hard data geostatistics has concentrated on the development of powerful methods based on stochastic theory. Though there have been recent moves to incorporate ancillary data in geostatistical analyses, insufficient attention has been paid to using modern methods of data display for the visualization of results. GIS can serve geostatistics by aiding geo-registration of data, facilitating spatial exploratory data analysis, providing a spatial context for interpolation and conditional simulation, as well as providing easy-to-use and effective tools for data display and visualization. The value of geostatistics for GIS lies in the provision of reliable interpolation methods with known errors, methods of upscaling and generalization, and for supplying multiple realizations of spatial patterns that can be used in environmental modeling. These stochastic methods are improving understanding of how errors in models of spatial processes accrue from errors in data or incompleteness in the structure of the models. New developments in GIS, based on ideas taken from map algebra, cellular automata and image analysis are providing high level programming languages for modeling dynamic processes such as erosion or the development of alluvial fans and deltas. Research has demonstrated that these models need stochastic inputs to yield realistic results. Non-stochastic tools such as fuzzy subsets have been shown to be useful for spatial analysis when probabilistic approaches are inappropriate or impossible. The conclusion is that in spite of differences in history and approach, the linkage of GIS, statistics and geostatistics provides a powerful, and complementary suite of tools for spatial analysis in the agricultural, earth and environmental sciences.     

Threats to Avifauna on Oceanic Islands

Abstract:  Results of the study by Blackburn et al. (2004 a ) of avifauna on oceanic islands suggest that distance from the mainland and time since European colonization have major influences on species extinctions and that island area is a significant but secondary contributing factor. After augmenting the data of the study on geographical properties for some of the islands they examined, we used a causal analysis approach with structural equation modeling to reexamine their conclusions. In our model geographical properties of islands, such as island area and isolation, were considered constraints on biological factors, such as the number of introduced mammalian predators and existing number of avifauna, that can directly or indirectly influence extinction. Of the variables we tested, island area had the greatest total influence on the threat of extinction due to its direct and indirect effects on the size of island avifauna. Larger islands had both a greater number of threatened bird species and more avifauna, increasing the number of species that could become threatened with extinction. Island isolation also had a significant, positive, and direct effect on threats to island avifauna because islands farther from the mainland had fewer current extant avifauna. Time since European colonization had a significant negative, but relatively weaker, influence on threats compared with the traditional biogeographic factors of island area and distance to the mainland. We also tested the hypothesis that the amount of threat is proportionally lower on islands that have had more extinctions (i.e., there is a "filter effect"). Because the proportion of bird extinctions potentially explained only 2.3% of the variation in the proportion of threatened species on islands, our results did not support this hypothesis. Causal modeling provided a powerful tool for examining threat of extinction patterns of known and hypothesized pathways of influence.