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.     

The Roche Limit

The role of tides in deforming and possibly disrupting a secondary body orbiting about a primary body has been known for a considerable time. This was first inspired by the observations of ocean tides on Earth and then seen as playing an important role in the formation and evolution of the Earth–Moon system. Finally, in the beginning of the 20th century it was generally thought to have a significant role in the formation of the solar system through the tidal disruption of the Sun. Here, an overview of the historical developments of the ideas concerned with tidal disruption of a secondary body that can lead to mass loss is given. Some discussion of possible extensions to consider more realistic situations where the secondary body may not be moving on a circular orbit and may not rotate so as to maintain the phase-on configuration to the primary body is also given.     

Upscaling of elastic properties for large scale geomechanical simulations

Large scale geomechanical simulations are being increasingly used to model the compaction of stress dependent reservoirs, predict the long term integrity of under‐ground radioactive waste disposals, and analyse the viability of hot‐dry rock geothermal sites. These large scale simulations require the definition of homogenous mechanical properties for each geomechanical cell whereas the rock properties are expected to vary at a smaller scale. Therefore, this paper proposes a new methodology that makes possible to define the equivalent mechanical properties of the geomechanical cells using the fine scale information given in the geological model. This methodology is implemented on a synthetic reservoir case and two upscaling procedures providing the effective elastic properties of the Hooke's law are tested. The first upscaling procedure is an analytical method for perfectly stratified rock mass, whereas the second procedure computes lower and upper bounds of the equivalent properties with no assumption on the small scale heterogeneity distribution. Both procedures are applied to one geomechanical cell extracted from the reservoir structure. The results show that the analytical and numerical upscaling procedures provide accurate estimations of the effective parameters. Furthermore, a large scale simulation using the homogenized properties of each geomechanical cell calculated with the analytical method demonstrates that the overall behaviour of the reservoir structure is well reproduced for two different loading cases. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Contribution of storms to groundwater recharge in the semi‐arid region of Karnataka,India

This paper describes the application of environmental isotopes and injected tracer techniques in estimating the contribution of storms as well as annual precipitation to groundwater recharge and its circulation, in the semi‐arid region of Bagepalli, Kolar district, Karnataka. Environmental isotopes ²H, ¹⁸O and ³H were used to study the effect of storms on the hydrological system, and an isotope balance was used to compute the contribution of a storm component to the groundwater. Some of the groundwater samples collected during the post‐storm periods were highly depleted in stable isotope content with higher deuterium excess relative to groundwater from the pre‐storm periods. Significant variation in deuterium excess in groundwater from the same area, collected in two different periods, indicates the different origin of air masses. The estimated recharge component of a storm event of 600 mm to the groundwater was found to be in the range of 117–165 mm. There was no significant variation in environmental tritium content of post‐storm and pre‐storm groundwater, indicating the fast circulation of groundwater in the system. After completion of the environmental isotope work, an injected radiotracer ³H technique was applied to estimate the direct recharge of total precipitation to the groundwater. The estimated recharge to the groundwater is 33 mm of the 550 mm annual precipitation during 1992. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Freshening of the Intermediate Waters in the Northern South China Sea over the PastSix Decades

The properties of salinity in the South China Sea (SCS), a significant marginal sea connecting the Pacific andIndian Oceans, are greatly influenced by the transport of fresh water flux between the two oceans. However, the long-termchanges in the intermediate water in the SCS have not been thoroughly studied due to limited data, particularly in relationto its thermodynamic variations. This study utilized reanalysis data products to identify a 60-year trend of freshening in theintermediate waters of the northern South China Sea (NSCS), accompanied by an expansion of low-salinity water. Thestudy also constructed salinity budget terms, including advection and entrainment processes, and conducted an analysis ofthe salinity budget to understand the impacts of external and internal dynamic processes on the freshening trend of theintermediate water in the NSCS. The analysis revealed that the freshening in the northwest Pacific Ocean and theintensification of intrusion through the Luzon Strait at intermediate levels are the primary drivers of the salinity changes inthe NSCS. Additionally, a weakened trend in the intensity of vertical entrainment also contributes to the freshening in theNSCS. This study offers new insights into the understanding of regional deep sea changes in response to variations in boththermodynamics and oceanic dynamic processes.     

Transport and Housing Security in the Klang Valley, Malaysia

This paper explores links between transport and housing security issues for the urban poor using the example of the Klang Valley in Malaysia. The interface between these issues is identified as a gap in the literature, including policy debates, on both housing and transport. A number of linkages are shown to be important and likely to be relevant in many cities of the South, especially those with rapid motorisation and large numbers of "squatters". A simple framework for understanding these linkages is presented. Key examples include displacement to make way for transport infrastructure and the impact on transport problems for the poor of policies affecting the location of urban poor housing, including relocation sites and transit accommodation. The case study of the Klang Valley is used to illustrate and test the relevance of a focus on this issue and the utility of the conceptual framework. Some policy implications of the investigation and case study are suggested.     

Long-term bed load transport rate based on aerial-photo and ground penetrating radar surveys of fan-delta growth, Coast Mountains, British Columbia

Sequential aerial photography, sonar bathymetry, ground-penetrating radar (GPR), and sediment sampling and analysis provide the basis for calculating the volumetric and mass rate of progradation of the delta of Fitzsimmons Creek, a steep, high-energy, debris-flow-dominated channel draining about 100 km² of the southern Coast Mountains of British Columbia. Fitzsimmons Creek is typical of small mountain rivers in the region. GPR imaging is used to define the pre-depositional morphology of the receiving basin, a technique that improves the accuracy of the volumetric survey. The 52-year record (1947–1999) of progradation yielded an average annual volumetric transport rate of 1.00±0.16×10⁴ m³ year⁻¹ for bed load, corresponding to a mass transport rate of 1.60±0.28×10⁴ Mg year⁻¹. Bed load yields are consistent with those obtained in hydrogeomorphically similar basins in the region and elsewhere. Decade-based annual rates, which vary from 0.64±0.11×10⁴ to 2.85±0.38×10⁴ Mg year⁻¹, provide poor estimates of the 52-year average. Indeed, the 52-year record may also not be long enough to fully integrate the significant fluctuations in the sediment efflux from Fitzsimmons Creek. The methodology proposed in this paper can be transferred to other comparable mountain environments worldwide.