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.     

Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

The influence of a non‐associated flow rule on the calculation of the factor of safety of soil slopes

This paper presents a method that incorporates a non‐associated flow rule into the limit analysis to investigate the influence of the dilatancy angle on the factor of safety for the slope stability analysis. The proposed method retain's the advantage of the upper bound method, which is simple and has no stress involvement in the calculation of the energy dissipation and the factor of safety. Copyright © 2001 John Wiley & Sons, Ltd.     

Late Holocene great earthquakes and relative sea‐level change at Kenai,southern Alaska

Kenai, located on the west coast of the Kenai Peninsula, Alaska, subsided during the great earthquake of AD 1964. Regional land subsidence is recorded within the estuarine stratigraphy as peat overlain by tidal silt and clay. Reconstructions using quantitative diatom transfer functions estimate co‐seismic subsidence (relative sea‐level rise) between 0.28±0.28 m and 0.70±0.28 m followed by rapid post‐seismic recovery. Stratigraphy records an earlier co‐seismic event as a second peat‐silt couplet, dated to ～1500–1400 cal. yr BP with 1.14±0.28 m subsidence. Two decimetre‐scale relative sea‐level rises are more likely the result of glacio‐isostatic responses to late Holocene and Little Ice Age glacier expansions rather than to co‐seismic subsidence during great earthquakes. Comparison with other sites around Cook Inlet, at Girdwood and Ocean View, helps in constructing regional patterns of land‐level change associated with three great earthquakes, AD 1964, ～950–850 cal. yr BP and ～1500–1400 cal. yr BP. Each earthquake has a different spatial pattern of co‐seismic subsidence which indicates that assessment of seismic hazard in southern Alaska requires an understanding of multiple great earthquakes, not only the most recent. All three earthquakes show a pre‐seismic phase of gradual land subsidence that marked the end of relative land uplift caused by inter‐seismic strain accumulation. Copyright © 2005 John Wiley & Sons, Ltd.     

A note on turbulent perturbations of salinity in a partially mixed estuary

Measurements of salinity perturbations in a partially mixed estuary have been used to evaluate the usefulness of an inductive salinometer and to determine some of the characteristics of the salinity perturbations. The salinometer performed satisfactorily under most conditions. Although internal wave like effects were present, the turbulence fluctuations were dominant. The salinity fluctuations and the turbulent fluxes $\text{sw}$ and $\text{su}$ were found to behave in a manner similar to the density fluctuations in a thermally stratified atmospheric boundary layer and a laboratory open channel flow. A quadrant analysis suggested that the contribution of each quadrant to the turbulent flux changed with Ri. The turbulence parameters ν and c_γ were found to decrease and increase respectively as Ri increases.     

Late Quaternary submarine bedforms and ice‐sheet flow in Gerlache Strait and on the adjacent continental shelf,Antarctic Peninsula

Geophysical data from Gerlache Strait, Croker Passage, Bismarck Strait and the adjacent continental shelf reveal streamlined subglacial bedforms that were produced at the bed of the Antarctic Peninsula Ice Sheet (APIS) during the last glaciation. The spatial arrangement and orientation of these bedforms record the former drainage pattern and flow dynamics of an APIS outlet up‐flow, and feeding into, a palaeo‐ice stream in the Western Bransfield Basin. Evidence suggests that together, they represent a single ice‐flow system that drained the APIS during the last glaciation. The ice‐sheet outlet flowed north/northeastwards through Gerlache Strait and Croker Passage and converged with a second, more easterly ice‐flow tributary on the middle shelf to form the main palaeo‐ice stream. The dominance of drumlins with low elongation ratios suggests that ice‐sheet outlet draining through Gerlache Strait was comparatively slower than the main palaeo‐ice stream in the Western Bransfield Basin, although the low elongation ratios may also partly reflect the lack of sediment. Progressive elongation of drumlins further down‐flow indicates that the ice sheet accelerated through Croker Passage and the western tributary trough, and fed into the main zone of streaming flow in the Western Bransfield Basin. Topography would have exerted a strong control on the development of the palaeo‐ice stream system but subglacial geology may also have been significant given the transition from crystalline bedrock to sedimentary strata on the inner–mid‐shelf. In the broader context, the APIS was drained by a number of major fast‐flowing outlets through cross‐shelf troughs to the outer continental shelf during the last glaciation. Copyright © 2004 John Wiley & Sons, Ltd.     

Low-cost volcano surveillance from space: case studies from Etna, Krafla, Cerro Negro, Fogo, Lascar and Erebus

 Satellite data offer a means of supplementing ground-based monitoring during volcanic eruptions, especially at times or locations where ground-based monitoring is difficult. Being directly and freely available several times a day, data from the advanced very high resolution radiometer (AVHRR) offers great potential for near real-time monitoring of all volcanoes across large (3000×3000 km) areas. Herein we describe techniques to detect and locate activity; estimate lava area, thermal flux, effusion rates and cumulative volume; and distinguish types of activity. Application is demonstrated using data for active lavas at Krafla, Etna, Fogo, Cerro Negro and Erebus; a pyroclastic flow at Lascar; and open vent systems at Etna and Stromboli. Automated near real-time analysis of AVHRR data could be achieved at existing, or cheap to install, receiving stations, offering a supplement to conventional monitoring methods. Received: 21 January 1997 / Accepted: 3 April 1997     

Simulating methane emission from a Chinese rice field as influenced by fertilizer and water level

Natural and agricultural wetlands are considered to be the major sources of global atmospheric methane (CH₄). A one‐dimensional model was developed to simulate methane emission and used to examine the influence of various physical processes on the rate of methane emission. Three processes involved in the methane emission are implemented in the model: production, reoxidation and transport. Three transport pathways were considered: diffusion across water–air or soil–air interfaces, ebullition and diffusion through plants. These pathways are influenced by soil properties, plant growth, water‐table conditions, temperature and external inputs (e.g. fertilizer). The model was used to examine the seasonal variation of the methane emission at a rice field in Hunan, China, which was observed during a field experiment for consecutive (early and late) rice seasons in 1992. The observed seasonal variations of methane emission, and role of plants in transporting methane to the atmosphere, are captured by the model simulation. Further model applications were conducted to simulate effects of fertilizer and water‐level condition on the methane emission. The results indicate that unfermented organic fertilizer produces a higher methane emission rate than mineral fertilizer. The simulations with treatments of a deep‐water covering and constant moisture reduced the methane emission. The rice field study provides a framework for further development of the model towards simulations based on spatially distributed variables (e.g. water table, soil temperature and vegetation) at a regional scale. Copyright © 2003 John Wiley & Sons, Ltd.     

Hydrodynamic control of phytoplankton in low salinity waters of the James River estuary,Virginia, U.S.A.

Autotrophic biomass and productivity as well as nutrient distributions and phytoplankton cell populations in the James River estuary, Virginia, were quantified both spatially and temporally over a 17-month period. Emphasis was placed on the very low salinity region of the estuary in order to gain information on the fate of freshwater phytoplankters. Differing amounts of freshwater plant biomass are advected into the estuary as living material, DOC or POC and the demonstrated variability of this input must play an important role in marine biogeochemical cycling.Late summer and fall maxima in both chlorophyll a and the photosynthetic production of particulate organic carbon in very low salinity regions were inversely correlated with river discharge.During periods of low river discharge greater than 50% of the chlorophyll a biomass measured at 0‰ disappeared within a narrow range of salinity (0–2‰). Cell enumeration data suggest that species introduced from the freshwater end-member tend to comprise the bulk of the biomass removed. Confounding factors, which may contribute to the regulation of both the abundance and species of phytoplankters mid-river, include the flocculation of colloidal material with phytoplankton cells, the presence of the turbidity maximum and the growth of endemic phytoplankton populations.An inverse relationship exists between the phytoplankton abundance in very low salinity waters and the abundance of biomass measured in the lower portion of the river (estuary). Thus, autotrophic production in the fresh and very low salinity areas may indirectly regulate the onset on the spring bloom in the estuary by controlling the amount of nutrients available.     

The stable isotope signature of kilometre-scale fracturedominated metamorphic fluid pathways, Mary Kathleen, Australia

Abstract Large calcite veins and pods in the Proterozoic Corella Formation of the Mount Isa Inlier provide evidence for kilometre-scale fluid transport during amphibolite facies metamorphism. These 10- to 100-m-scale podiform veins and their surrounding alteration zones have similar oxygen and carbon isotopic ratios throughout the 200 × 10-km Mary Kathleen Fold Belt, despite the isotopic heterogeneity of the surrounding wallrocks. The fluids that formed the pods and veins were not in isotopic equilibrium with the immediately adjacent rocks. The pods have δ¹³C_calcite values of –2 to –7% and δ¹⁸O_calcite values of 10.5 to 12.5%. Away from the pods, metadolerite wallrocks have δ¹⁸O_whole-rock values of 3.5 to 7%. and unaltered banded calc-silicate and marble wallrocks have δ¹³C_calcite of –1.6 to –0.6%, and δ¹⁸O_calcite of 18 to 21%. In the alteration zones adjacent to the pods, the δ¹⁸O values of both metadolerite and calc-silicate rocks approach those of the pods. Large calcite pods hosted entirely in calc-silicates show little difference in isotopic composition from pods hosted entirely in metadolerite. Thus, 100- to 500-m-scale isotopic exchange with the surrounding metadolerites and calc-silicates does not explain the observation that the δ¹⁸O values of the pods are intermediate between these two rock types. Pods hosted in felsic metavolcanics and metasiltstones are also isotopically indistinguishable from those hosted in the dominant metadolerites and calc-silicates. These data suggest the veins are the product of infiltration of isotopically homogeneous fluids that were not derived from within the Corella Formation at the presently exposed crustal level, although some of the spread in the data may be due to a relatively small contribution from devolatilization reactions in the calc-silicates, or thermal fluctuations attending deformation and metamorphism. The overall L-shaped trend of the data on plots of δ¹³C vs. δ¹⁸O is most consistent with mixing of large volumes of externally derived fluids with small volumes of locally derived fluid produced by devolatilization of calc-silicate rocks. Localization of the vein systems in dilatant sites around metadolerite/calc-silicate boundaries indicates a strong structural control on fluid flow, and the stable isotope data suggest fluid migration must have occurred at scales greater than at least 1 km. The ultimate source for the external fluid is uncertain, but is probably fluid released from crystallizing melts derived from the lower crust or upper mantle. Intrusion of magmas below the exposed crustal level would also explain the high geothermal gradient calculated for the regional metamorphism.