David Rassam

Conceptualisation and application of models for groundwater-surface water interactions and nitrate attenuation potential in riparian zones

Riparian zones can provide a protective buffer between streams and adjacent land-based activities, by removing nitrate from shallow groundwater flowing through them. Catchment scale water quality models are useful tools for predicting catchment behaviour under various climatic conditions and land use scenarios, but most do not account for the effect of riparian buffer zones. In this paper, we present conceptual models for surface water-groundwater interactions and formulate analytical mathematical functions that describe nitrate removal in the riparian zone. We restrict nitrate attenuation capacity to potential denitrification only and present sample calculations based on the limited field data available. The models are classed into two types according to their applicability to either ephemeral or perennial streams. In ephemeral, low-order streams with the potential to form a perched water table, a simple bucket model is used. During events, stream water flows laterally into the riparian zone and may be denitrified while residing there before discharging back to the surface water system. In perennial middle-order streams, nitrate removal may occur either as base flow intercepts the root zone or when water is temporarily stored in stream banks during flood events. We incorporate these concepts within a GIS modelling framework and investigate the potential of riparian zones to reduce nitrate delivery to streams in the Maroochy catchment located in South East Queensland, Australia. The modelling results for the Maroochy catchment show that the optimum rooting depth is about 4m and that increasing the riparian buffer width beyond 10m yields little further reduction in nitrate. The potential nitrate removal capacity per unit length of riparian buffer in each sub-catchment is an attribute that can be used to help to prioritise riparian rehabilitation activities aimed at reducing stream nitrogen loads.

Soil–water and solute movement under precision irrigation: knowledge gaps for managing sustainable root zones

Precision irrigation involves the accurate and precise application of water to meet the specific requirements of individual plants or management units and minimize adverse environmental impact. Under precision irrigation applications, water and associated solute movement will vary spatially within the root zone and excess water application will not necessarily result in deep drainage and leaching of salt below the root zone. This paper estimates that 10% of the irrigated land area (producing as much as 40% of the total annual revenue from irrigated land) could be adversely affected by root zone salinity resulting from the adoption of precision irrigation within Australia. The cost of increases in root zone salinisation due to inappropriate irrigation management in the Murray and Murrumbidgee irrigation areas was estimated at AUD 245 million (in 2000/01) or 13.5% of the revenue from these cropping systems. A review of soil–water and solute movement under precision irrigation systems highlights the gaps in current knowledge including the mismatch between the data required by complex, process-based soil–water or solute simulation models and the data that is easily available from soil survey and routine soil analyses. Other major knowledge gaps identified include: (a) effect of root distribution, surface evaporation and plant transpiration on soil wetted patterns, (b) accuracy and adequacy of using simple mean values of root zone soil salinity levels to estimate the effect of salt on the plant, (c) fate of solutes during a single irrigation and during multiple irrigation cycles, and (d) effect of soil heterogeneity on the distribution of water and solutes in relation to placement of water. Opportunities for research investment are identified across a broad range of areas including: (a) requirements for soil characterisation, (b) irrigation management effects, (c) agronomic responses to variable water and salt distributions in the root zone, (d) potential to scale or evaluate impacts at various scales, (e) requirements for simplified soil–water and solute modelling tools, and (f) the need to build skills and capacity in soil–water and solute modelling.

Development and application of the Riparian Mapping Tool to identify priority rehabilitation areas for nitrogen removal in the Tully-Murray basin, Queensland, Australia

One feature of riparian zones is their ability to significantly reduce the nitrogen loads entering streams by removing nitrate from the groundwater. A novel GIS model was used to prioritise riparian rehabilitation in catchments. It is proposed that high-priority areas are those with a high potential for riparian denitrification and have nearby land uses that generate high nitrogen loads. For this purpose, we defined the Rehabilitation Index, which is the product of two other indices, the Nitrate Removal Index and the Nitrate Interception Index. The latter identifies the nitrate contamination potential for each raster cell in the riparian zone by examining the extent and proximity of agricultural urban land uses. The former is estimated using a conceptual model for surface-groundwater interactions in riparian zones associated with middle-order gaining perennial streams, where nitrate is removed via denitrification when the base flow interacts with the carbon-rich riparian sediments before discharging to the streams. Riparian zones that are relatively low in the landscape, have a flat topography, and have soils of medium hydraulic conductivity are most conducive to denitrification. In the present study, the model was implemented in the Tully-Murray basin, Queensland, Australia, to produce priority riparian rehabilitation area maps.

Predicting Aquifer Response Time for Application in Catchment Modeling

It is well established that changes in catchment land use can lead to significant impacts on water resources. Where land-use changes increase evapotranspiration there is a resultant decrease in groundwater recharge, which in turn decreases groundwater discharge to streams. The response time of changes in groundwater discharge to a change in recharge is a key aspect of predicting impacts of land-use change on catchment water yield. Predicting these impacts across the large catchments relevant to water resource planning can require the estimation of groundwater response times from hundreds of aquifers. At this scale, detailed site-specific measured data are often absent, and available spatial data are limited. While numerical models can be applied, there is little advantage if there are no detailed data to parameterize them. Simple analytical methods are useful in this situation, as they allow the variability in groundwater response to be incorporated into catchment hydrological models, with minimal modeling overhead. This paper describes an analytical model which has been developed to capture some of the features of real, sloping aquifer systems. The derived groundwater response timescale can be used to parameterize a groundwater discharge function, allowing groundwater response to be predicted in relation to different broad catchment characteristics at a level of complexity which matches the available data. The results from the analytical model are compared to published field data and numerical model results, and provide an approach with broad application to inform water resource planning in other large, data-scarce catchments.

Effluent flux prediction in variably saturated soil zones within a septic tank—soil absorption trench

The treatment and hydraulic mechanisms in a septic tank–soil absorption system (SAS) are highly influenced by the clogging layer or biomat zone which develops on bottom and lower sidewall surfaces within the trench. Flow rates through the biomat and sub-biomat zones are governed largely by the biomat hydraulic properties (resistance and hydraulic conductivity) and the unsaturated hydraulic conductivity of the underlying soil. One- and 2-dimensional models were used to investigate the relative importance of sidewall and vertical flow rates and pathways in SAS. Results of 1-dimensional modelling show that several orders of magnitude variation in saturated hydraulic conductivity (Ks) reduce to a 1 order of magnitude variation in long-term flow rates. To increase the reliability of prediction of septic trench hydrology, HYDRUS-2D was used to model 2-dimensional flow. In the permeable soils, under high trench loading, effluent preferentially flowed in the upper region of the trench where no resistant biomat was present (the exfiltration zone). By comparison, flow was more evenly partitioned between the biomat zones and the exfiltration zones of the low permeability soil. An increase in effluent infiltration corresponded with a greater availability of exfiltration zone, rather than a lower resistance of biomat. Results of modelling simulations demonstrated the important role that a permeable A horizon may play in limiting surface surcharge of effluent under high trench hydraulic loading.

Probabilistic modelling and uncertainty analysis of flux and water balance changes in a regional aquifer system due to coal seam gas development

Large scale development of coal seam gas (CSG) is occurring in many sedimentary basins around the world including Australia, where commercial production of CSG has started in the Surat and Bowen basins. CSG development often involves extraction of large volumes of water that results in depressurising aquifers that overlie and/or underlie the coal seams thus perturbing their flow regimes. This can potentially impact regional aquifer systems that are used for many purposes such as irrigation, and stock and domestic water. In this study, we adopt a probabilistic approach to quantify the depressurisation of the Gunnedah coal seams and how this impacts fluxes to, and from the overlying Great Artesian Basin (GAB) Pilliga Sandstone aquifer. The proposed method is suitable when effects of a new resource development activity on the regional groundwater balance needs to be assessed and account for large scale uncertainties in the groundwater flow system and proposed activity. The results indicated that the extraction of water and gas from the coal seam could potentially induce additional fluxes from the Pilliga Sandstone to the deeper formations due to lowering pressure heads in the coal seams. The median value of the rise in the maximum flux from the Pilliga Sandstone to the deeper formations is estimated to be 85ML/year, which is considered insignificant as it forms only about 0.29% of the Long Term Annual Average Extraction Limit of 30GL/year from the groundwater management area. The probabilistic simulation of the water balance components indicates only small changes being induced by CSG development that influence interactions of the Pilliga Sandstone with the overlying and underlying formations and with the surface water courses. The current analyses that quantified the potential maximum impacts of resource developments and how they influences the regional water balance, would greatly underpin future management decisions.

Non-stationarity of low flows and their relevance to river modelling during drought periods

Changes in groundwater storage lead to a reduction in groundwater contribution to river flow and present as non-stationarity, especially during low-flow conditions. Conventional river models typically ignore this non-stationarity, and, hence, their predictions of declines in low flows during drought periods are likely to be compromised. The present study assesses non-stationarity and highlights its implications for river modelling. A quantile regression analysis showed non-stationarity of low flows in the Namoi catchment (Australia), with statistically significant downward trends in the 10th percentile of log-transformed baseflow (10-LTB). This highlighted the usefulness of the 10-LTB metric to identify non-stationarity and, hence, alert modellers to the importance of adopting models that explicitly account for groundwater processes when modelling such river systems.

Trajectory river modelling – a decision-support tool to help manage multiple risks associated with planning around variable water resources

The application of river-system models to inform water-resource planning and management is a growing global phenomenon. This requires models to be applied so that they are useful to water decision makers charged with setting targets that provide adequate water flows to sustain landholders and communities. This article examines why and how the innovative application of river-system models can facilitate interactions between water science and water management in Australia’s Murray–Darling Basin (the Basin). A trajectory river-modelling method was applied to run multiple short historical climate sequences through a river-system model to provide historical probabilities. These can allow better assessment of the risks and impacts associated with stream flow and water availability. This method allows known historical variability to be presented, and produces relevant results for a 10–15-year water-sharing plan lifetime. The benefits were demonstrated in the Basin’s Lachlan Catchment where modelled river-flow results demonstrated the increased variability between shorter 15-year sequences than for a single 114-year run. This approach highlighted the benefits of expressing modelling results as historical probabilities to inform short-term and strategic water-planning efforts.