Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

Chemical alteration of tephra in the depositional environment: theoretical stability modelling

The study of the chemical stability of vitreous material in aqueous media is well‐established. There has to date been little consideration of the implications of variations in the chemical durability of tephra in Quaternary tephrochronology. Chemical alteration can take the form of cationic leaching from the matrix, or complete destruction of the silica network, either of which could constrain the ability to chemically identify distal tephra. Here we apply established models of vitreous durability to the published chemical analyses of a large number of Icelandic tephras in order to predict their relative durabilities under equivalent conditions. This suggests that some important tephras have relatively poor chemical stability, and that rhyolitic tephras are, in general, more stable than basaltic. We conclude that tephras should be expected to show predictable differential chemical stability in the post‐depositional environment. Copyright © 2003 John Wiley & Sons, Ltd.     

Seasonal variations in prey selection by estuarine black-headed gulls (Larus ridibundus)

The number of black-headed gulls (Larus ridibundus) in the Clyde Estuary is large. In summer the average density has reached 1350 gulls km^?2 and in winter 180 gulls km^?2. This paper compares prey selection and feeding efficiency in gulls during summer and winter on tidal flats, and considers how seasonal differences may be adaptations to cope with seasonal changes in prey availability.Gross and net rates of energy intake were highest in summer because gulls captured more of the polychaete N. diversicolor than the amphipod C. volutator. In winter, gulls selected for C. volutator and therefore an energetically less profitable diet. Throughout the year gulls selected more C. volutator relative to N. diversicolor than expected on energetic grounds and so apparently did not maximize potential net rate of energy intake.Gulls used three techniques to capture prey and made most intensive use of the ‘crouch’ technique. Crouching gulls attained a much higher net rate of energy intake than ‘upright’ or ‘paddling’ gulls.A log-linear model showed that (a) season, water depth and gull density determined feeding technique and (b) feeding technique and season independently determined foraging success and prey selection. Thus gull density and water depth acted on prey selection through imposed variations in feeding technique.Reasons for gulls selecting energetically unprofitable C. volutator and for the use of several distinct feeding techniques are discussed.     

A MESO-β SCALE SIMULATION OF THE EFFECTS OF BOREAL FOREST ECOSYSTEM ON THE LOWER ATMOSPHERE

Based on the Intensive Field Campaign(IFC-1)data of Boreal Ecosystem-Atmosphere Study(BOREAS).a three-dimensional meso-β scale model is used to simulate the effect of boreal forests onthe lower atmosphere.A fine horizontal resolution of 2 km×2 km is used in order to distinguish thevegetative heterogeneity in the boreal region.A total of 20×25 grid points cover the entire sub-modeling area in BOREAS' South Study Area(SSA).The ecosystem types and their coverage ineach grid square are extracted from the North American Land Cover Characteristics Data Base(NALCCD)generated by the U.S.Geographical Survey(USGS)and the University of Nebraska-Lincoln(UNL).The topography of the study area is taken from the Digital Elevation Map(DEM)of USGS.The model outputs include the components of the energy balance budget within the canopyand at the ground.the turbulence parameters in the atmospheric boundary layer and the wind.temperature and humidity profiles extending up to a height of 1500 m.In addition to the fine timeand spatial step,the unique feature of the present model is the incorporation of both dynamic andbiological effects of the Boreal forest into the model parameterization scheme.The model resultscompare favorably with BOREAS' IFC-1 data in 1994 when the forest was in the luxuriant growingperiod.     

On salinity regimes and the vertical structure of residual flows in narrow tidal estuaries

An examination is made of the circulation in narrow estuaries subject to a predominant tidal forcing. Velocity structures are derived separately for residual flow components associated with (a) river flow, (b) wind stress, (c) a well-mixed longitudinal density gradient and (d) a fully stratified saline wedge. Dimensionless parameters are introduced to indicate the magnitude of each component and these parameters are evaluated for 9 major estuaries, thereby revealing their sensitivity to each component.For a channel of constant breadth and depth, formulae are deduced for the length of saline intrusion, L. Comparisons with observed data show that such formulae may be used with confidence to predict changes in L arising from variations in river flow, tidal range or channel depths.The level of stratification is shown to be related to a product of two parameters, one associated with velocity structure and a second involving the square of the ‘flow ratio’ $\text{u}\text{u}\text{?}$ (i.e. residual velocity/amplitude of the tidal velocity). This relationship provides a simple classification system for estuarine stratification which can be used to indicate the sensitivity of any particular estuary to changing conditions.     

Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

The hydrological influence of black locust plantations in the loess area of northwest China

Black locust (Robina pseudoacacia) has become one of the most important shelter species in the loess area of northwest China. This paper summarizes recent research concerning its hydrological influence, including canopy interception, litter absorption capacity, its effect on rainfall kinetic energy, infiltration rates, surface runoff, soil moisture, and evapotranspiration, and its role in soil conservation. Several predictive models are listed. on the basis of existing results, optimum characteristics for an effective plantation are defined, and problems requiring further research are identified.     

A new stereo‐analytical method for determination of removal blocks in discontinuous rock masses

The paper provides a new stereo‐analytical method, which is a combination of the stereographic method and analytical methods, to separate finite removable blocks from the infinite and tapered blocks in discontinuous rock masses. The methodology has applicability to both convex and concave blocks. Application of the methodology is illustrated through examples. Addition of this method to the existing block theory procedures available in the literature improves the capability of block theory in solving practical problems in rock engineering. Copyright © 2003 John Wiley & Sons, Ltd.     

Hysteresis‐based analysis of overland metal transport

Introducing a concept of equivalent mass depth of flow, this study describes the phenomenon of non‐point source pollutant (metal) transport for pavement (or overland) flow in analogy with wave propagation in wide open channels. Hysteretic and normal mass rating curves are developed for runoff rate and mass of 12 dissolved and particulate‐bound metal elements (pollutants) using the rainfall‐runoff and water quality data of the 15 × 20 m² instrumented pavement in Cincinnati, USA. Normal mass rating curves developed for easy computation of pollutant load are found to be of a form similar to Manning's, which is valid for open channel flows. Based on the hysteresis analysis, wave types for dissolution and mixing of particulate‐bound metals are identified. The analysis finds that the second‐order partial‐differential equation normally used for metal transport does not have the efficacy to describe fully the strong non‐linear phenomena such as is described for various metal elements by dynamic waves. In addition, the proportionality concept of the popular SCS‐CN concept is extended for determining the potential maximum metal mass M_p of all the 12 elements transported by a rain storm and related to the antecedent dry period (ADP). For the primary metal zinc element, M_p is found to increase with the ADP and vice versa. Copyright © 2003 John Wiley & Sons, Ltd.