Groundwater flow and solute transport at the Mourquong saline-water disposal basin, Murray Basin, southeastern Australia

Saline groundwater and drainage effluent from irrigation are commonly stored in some 200 natural and artificial saline-water disposal basins throughout the Murray-Darling Basin of Australia. Their impact on underlying aquifers and the River Murray, one of Australia's major water supplies, is of serious concern. In one such scheme, saline groundwater is pumped into Lake Mourquong, a natural groundwater discharge complex. The disposal basin is hydrodynamically restricted by low-permeability lacustrine clays, but there are vulnerable areas in the southeast where the clay is apparently missing. The extent of vertical and lateral leakage of basin brines and the processes controlling their migration are examined using (1) analyses of chloride and stable isotopes of water (²H/¹H and ¹⁸O/¹⁶O) to infer mixing between regional groundwater and lake water, and (2) the variable-density groundwater flow and solute-transport code SUTRA. Hydrochemical results indicate that evaporated disposal water has moved at least 100 m in an easterly direction and that there is negligible movement of brines in a southerly direction towards the River Murray. The model is used to consider various management scenarios. Salt-load movement to the River Murray was highest in a "worst-case" scenario with irrigation employed between the basin and the River Murray. Present-day operating conditions lead to little, if any, direct movement of brine from the basin into the river. Electronic Publication     

Analysis of A New Model for Unsaturated Flow in Porous Media Including Hysteresis and Dynamic Effects

A new model for unsaturated flow in porous media, including capillary hysteresis and dynamic capillary effects, is analyzed. Existence and uniqueness of solutions are established and qualitative and quantitative properties of (particular) solutions are analyzed. Some results of numerical computations are given. The model under consideration incorporates simple play-type hysteresis and a dynamic term (time-derivative with respect to water content) in the capillary relation. Given an initial water content distribution, the model determines which parts of the flow domain are in drainage and which parts are in imbibition. The governing equations can be recast into an elliptic problem for fluid pressure and an evolution equation for water content. Standard methods are used to obtain numerical results. A comparison is given between J.R. Philip's semi-explicit similarity solution for horizontal redistribution in an infinite one-dimensional domain and solutions of the new model.     

应用自然电场法研究地下水流场

讨论了不同水文地质条件下自然电场“8”字图特征,对评价输水管线等渗漏引起的环境问题具有指导意义。     

Estimation of Scalar Source/Sink Distributions in Plant Canopies Using Lagrangian Dispersion Analysis: Corrections for Atmospheric Stability and Comparison with a Multilayer Canopy Model

Source/sink distributions of heat, water vapour andCO₂ within a rice canopy were inferred using aninverse Lagrangian dispersion analysis and measuredmean profiles of temperature, specific humidity andCO₂ mixing ratio. Monin–Obukhov similarity theorywas used to account for the effects of atmosphericstability on _w(z), the standard deviation ofvertical velocity and _L(z), the Lagrangian timescale of the turbulence. Classical surface layer scaling was applied in the inertial sublayer (z > z_ruf)using the similarity parameter = (z - d)/L, where z is height above ground, d is the zero plane displacementheight for momentum, L is the Obukhov length,and z_ruf 2.3h_c, where h_c iscanopy height. A single length scale h_c, was usedfor the stability parameter 3 = h_c/L in the height range 0.25 < z/h_c < 2.5. This choice is justified by mixing layer theory, which shows that within the roughness sublayer there is one dominant turbulence length scaledetermined by the degree of inflection in the windprofile at the canopy top. In the absence of theoretical or experimental evidence for guidance,standard Monin–Obukhov similarity functions, with = h_c/L, were used to calculate the stabilitydependence of _w(z) and _L(z) in the roughness sublayer. For z/h_c < 0.25 the turbulence length and time scales are influenced by the presence of the lowersurface, and stability effects are minimal. With theseassumptions there was excellent agreement between eddycovariance flux measurements and deductions from theinverse Lagrangian analysis. Stability correctionswere particularly necessary for night time fluxes whenthe atmosphere was stably stratified.The inverse Lagrangian analysis provides a useful toolfor testing and refining multilayer canopy models usedto predict radiation absorption, energy partitioningand CO₂ exchanges within the canopy and at thesoil surface. Comparison of model predictions withsource strengths deduced from the inverse analysisgave good results. Observed discrepancies may be dueto incorrect specification of the turbulent timescales and vertical velocity fluctuations close to theground. Further investigation of turbulencecharacteristics within plant canopies is required toresolve these issues.     

Salt transport on islands in the Okavango Delta: Numerical investigations

This study uses a numerical model to investigate the groundwater flow and salt transport mechanisms below islands in the Okavango Delta. Continuous evapotranspiration on the islands results in accumulation of solutes and the formation of a saline boundary layer, which may eventually become unstable. A novel Lagrangian method is employed in this study and compared to other numerical methods. The numerical results support the geophysical observations of density fingering on Thata Island. However, the process is slow and it takes some hundreds of years until density fingering is triggered. The results are sensitive to changes of the hydraulic gradient and the evapotranspiration rate. Small changes may lead to different plume developments. Results further demonstrate that density effects may be entirely overridden by lateral flow on islands embedded in a sufficiently high regional hydraulic gradient.     

A spatially and temporally adaptive solution of Richards’ equation

Efficient, robust simulation of groundwater flow in the unsaturated zone remains computationally expensive, especially for problems characterized by sharp fronts in both space and time. Standard approaches that employ uniform spatial and temporal discretizations for the numerical solution of these problems lead to inefficient and expensive simulations. In this work, we solve Richards’ equation using adaptive methods in both space and time. Spatial adaption is based upon a coarse grid solve and a gradient error indicator using a fixed-order approximation. Temporal adaption is accomplished using variable order, variable step size approximations based upon the backward difference formulas up to fifth order. Since the advantages of similar adaptive methods in time are now established, we evaluate our method by comparison with a uniform spatial discretization that is adaptive in time for four different one-dimensional test problems. The numerical results demonstrate that the proposed method provides a robust and efficient alternative to standard approaches for simulating variably saturated flow in one spatial dimension.     

Integrated surface–groundwater flow modeling: A free-surface overland flow boundary condition in a parallel groundwater flow model

Interactions between surface and groundwater are a key component of the hydrologic budget on the watershed scale. Models that honor these interactions are commonly based on the conductance concept that presumes a distinct interface at the land surface, separating the surface from the subsurface domain. These types of models link the subsurface and surface domains via an exchange flux that depends upon the magnitude and direction of the hydraulic gradient across the interface and a proportionality constant (a measure of the hydraulic connectivity). Because experimental evidence of such a distinct interface is often lacking in field systems, there is a need for a more general coupled modeling approach.     

Influence of uncertain boundary conditions and model structure on flood inundation predictions

In this study, the GLUE methodology is applied to establish the sensitivity of flood inundation predictions to uncertainty of the upstream boundary condition and bridges within the modelled region. An understanding of such uncertainties is essential to improve flood forecasting and floodplain mapping. The model has been evaluated on a large data set. This paper shows uncertainty of the upstream boundary can have significant impact on the model results, exceeding the importance of model parameter uncertainty in some areas. However, this depends on the hydraulic conditions in the reach e.g. internal boundary conditions and, for example, the amount of backwater within the modelled region. The type of bridge implementation can have local effects, which is strongly influenced by the bridge geometry (in this case the area of the culvert). However, the type of bridge will not merely influence the model performance within the region of the structure, but also other evaluation criteria such as the travel time. This also highlights the difficulties in establishing which parameters have to be more closely examined in order to achieve better fits. In this study no parameter set or model implementation that fulfils all evaluation criteria could be established. We propose four different approaches to this problem: closer investigation of anomalies; introduction of local parameters; increasing the size of acceptable error bounds; and resorting to local model evaluation. Moreover, we show that it can be advantageous to decouple the classification into behavioural and non-behavioural model data/parameter sets from the calculation of uncertainty bounds.     

Pore scale modeling of immiscible and miscible fluid flows using smoothed particle hydrodynamics

A numerical model based on smoothed particle hydrodynamics (SPH) was developed and used to simulate immiscible and miscible fluid flows in porous media and to study effects of pore scale heterogeneity and anisotropy on such flows.     

A bed‐load transport model for rough turbulent open‐channel flows on plane beds

Data from flume studies are used to develop a model for predicting bed‐load transport rates in rough turbulent two‐dimensional open‐channel flows moving well sorted non‐cohesive sediments over plane mobile beds. The object is not to predict transport rates in natural channel flows but rather to provide a standard against which measured bed‐load transport rates influenced by factors such as bed forms, bed armouring, or limited sediment availability may be compared in order to assess the impact of these factors on bed‐load transport rates. The model is based on a revised version of Bagnold's basic energy equation i_bs_b = e_bω, where i_b is the immersed bed‐load transport rate, ω is flow power per unit area, e_b is the efficiency coefficient, and s_b is the stress coefficient defined as the ratio of the tangential bed shear stress caused by grain collisions and fluid drag to the immersed weight of the bed load. Expressions are developed for s_b and e_b in terms of G, a normalized measure of sediment transport stage, and these expressions are substituted into the revised energy equation to obtain the bed‐load transport equation i_b = ω G ^3·4. This equation applies regardless of the mode of bed‐load transport (i.e. saltation or sheet flow) and reduces to i_b = ω where G approaches 1 in the sheet‐flow regime. That i_b = ω does not mean that all the available power is dissipated in transporting the bed load. Rather, it reflects the fact that i_b is a transport rate that must be multiplied by s_b to become a work rate before it can be compared with ω. It follows that the proportion of ω that is dissipated in the transport of bed load is i_bs_b/ω, which is approximately 0·6 when i_b = ω. It is suggested that this remarkably high transport efficiency is achieved in sheet flow (1) because the ratio of grain‐to‐grain to grain‐to‐bed collisions increases with bed shear stress, and (2) because on average much more momentum is lost in a grain‐to‐bed collision than in a grain‐to‐grain one. Copyright © 2006 John Wiley & Sons, Ltd.