An Excel macro for generating trilinear plots

This computer note describes a method for creating trilinear plots in Microsoft Excel. Macros have been created in MS Excel's internal language: Visual Basic for Applications (VBA). A simple form has been set up to allow the user to input data from an Excel worksheet. The VBA macro is used to convert the triangular data (which consist of three columns of percentage data) into X-Y data. The macro then generates the axes, labels, and grid for the trilinear plot. The X-Y data are plotted as scatter data in Excel. By providing this macro in Excel, users can create trilinear plots in a quick, inexpensive manner.     

Simulating Piecewise-Linear Surface Water and Ground Water Interactions with MODFLOW

The standard MODFLOW packages offer limited capabilities to model piecewise-linear boundary conditions to describe ground water–surface water interaction. Specifically, MODFLOW is incapable of representing a Cauchy-type boundary with different resistances for discharge or recharge conditions. Such a more sophisticated Cauchy boundary condition is needed to properly represent surface waters alternatively losing water through the bottom (high resistance) or gaining water mostly near the water surface (low resistance). One solution would be to create a new package for MODFLOW to accomplish this. However, it is also possible to combine multiple instances of standard packages in a single cell to the same effect. In this specific example, the general head boundary package is combined with the drain package to arrive at the desired piecewise-linear behavior. In doing so, the standard USGS MODFLOW version can be used without any modifications at the expense of a minor increase in preprocessing and postprocessing and computational effort. The extra preprocessing for creating the input and extra postprocessing to determine the water balance in terms of the physical entities from the MODFLOW cell fluxes per package can be taken care of by a user interface.     

A User Interface for the Kansas Geological Survey Slug Test Model

The Kansas Geological Survey (KGS) developed a semianalytical solution for slug tests that incorporates the effects of partial penetration, anisotropy, and the presence of variable conductivity well skins. The solution can simulate either confined or unconfined conditions. The original model, written in FORTRAN, has a text-based interface with rigid input requirements and limited output options. We re-created the main routine for the KGS model as a Visual Basic macro that runs in most versions of Microsoft Excel and built a simple-to-use Excel spreadsheet interface that automatically displays the graphical results of the test. A comparison of the output from the original FORTRAN code to that of the new Excel spreadsheet version for three cases produced identical results.     

Piecewise-continuous boundaries using the MODFLOW FHB and MT3DMS HSS packages

It is often necessary to simulate the influx to a groundwater model of water containing dissolved contaminants. Until fairly recently, users of MODFLOW and MT3DMS were restricted to varying the flux of water and contaminants on a stress-period basis: when a time-varying loading pattern required simulation, the modeler's only recourse was to discretize the model into many stress periods. From a practical standpoint this is cumbersome, while from a technical standpoint it requires that the modeler define a priori an appropriate time discretization that can accurately reproduce time-varying flow and mass loading. This is particularly undesirable when attempting to infer a time-varying flow or mass loading using inverse methods. The advent of the Flow and Head Boundary (FHB) package in the late 1990s greatly mitigated these limitations from the flow perspective. The recent release of the Hydrocarbon Spill Source (HSS) package for MT3DMS has essentially removed the limitation from the contaminant mass perspective. This Methods Note verifies the FHB and HSS packages by comparison with more commonly used boundary packages and highlights some benefits of their combined use, with reference to the reconstruction of historic flow and mass fluxes through inverse modeling. (Note: The Flow and Head Boundary and Hydrocarbon Spill Source packages are referred to throughout as "FHB" and "HSS", respectively,--that is, omitting version number suffixes--as the discussion presented should apply to all releases of each package.).     

Using MODFLOW 2000 to model ET and recharge for shallow ground water problems

In environments with shallow ground water elevation, small changes in the water table can cause significant variations in recharge and evapotranspiration fluxes. Particularly, where ground water is close to the soil surface, both recharge and evapotranspiration are regulated by a thin unsaturated zone and, for accuracy, must be represented using nonconstant and often nonlinear relationships. The most commonly used ground water flow model today, MODFLOW, was originally designed with a modular structure with independent packages representing recharge and evaporation processes. Systems with shallow ground water, however, may be better represented using either a recharge function that varies with ground water depth or a continuous recharge and evapotranspiration function that is dependent on depth to water table. In situations where the boundaries between recharging and nonrecharging cells change with time, such as near a seepage zone, a continuous ground water flux relationship allows recharge rates to change with depth rather than having to calculate them at each stress period. This research article describes the modification of the MODFLOW 2000 recharge and segmented evapotranspiration packages into a continuous recharge-discharge function that allows ground water flux to be represented as a continuous process, dependent on head. The modifications were then used to model long-term recharge and evapotranspiration processes on a saline, semiarid floodplain in order to understand spatial patterns of salinization, and an overview of this process is given.     

A visual basic spreadsheet macro for geochemical background analysis

A Visual Basic macro entitled BACKGROUND calculates geochemical background values of chemical parameters and estimates threshold values separating background data from anomalies. The macro uses two statistical methods, the iterative 2-sigma technique and the calculated distribution function, and integrates these model-based objective methods into a widely accessible platform (i.e., MS Excel). The macro offers the possibility for automated processing of geochemical data and enables an automated generation of background range and threshold values for chemical parameters.     

A spreadsheet method of estimating best-fit hydraulic gradients using head data from multiple wells

Hydraulic gradients from planar water tables, or piezometric surfaces, and horizontal flow regimes can be quickly and conveniently calculated from data sets involving numerous wells. The matrix-solving functions of a modem spreadsheet program (Excel) were used to determine the equation of a water-table plane, Ax + By + Cz - D = 0, and the equation coefficients were then used to determine the magnitude of the hydraulic gradient, according to gradient = square root of A2 + B2/C2, and its direction, according to alpha = arctan B/A, where alpha is the angle measured from the x-axis. A pre-prepared Excel file constructed to handle data from up to 20 wells at once is available for free downloading at www.geo.ku.edu/hydro/KUHydro.html.     

Development and application of a stochastic optimization model for groundwater management: crop pattern and conjunctive use consideration

The need for irrigation water in arid and semi-arid regions is mostly supplied by groundwater. Furthermore, the agricultural development in these areas is not generally based on a comprehensive plan, which can cause aquifers depletion. On the other hand, to properly manage an aquifer and to have an optimal crop plan, the stochastic nature of the different parameters of a groundwater system such as groundwater recharge and water demands should be taken into consideration. In this paper, we develop an explicit stochastic optimization model for Firouzabad aquifer in Iran. This formulation is based on the first and second moment analysis for groundwater head which has been initially proposed for surface water resources management by Fletcher and Ponnambalam. We extend the model to create a new random withdrawal policy for conjunctive use setting in which the randomness in available precipitation is taken into account. The interesting point is that the model provides the respective probabilities of shortage and surplus without imposing the extra decision variables into the optimization model. A genetic-based algorithm is used to solve the stochastic nonlinear and non-convex formulation. The outcome results indicate that the current crop pattern should be changed, that is, the allocated areas of some crops have to be meaningfully reduced. Finally, to validate our model efficiency, we demonstrate that how much close the statistical characteristics obtained from the optimization model are to those estimated from the Monte Carlo simulation. Furthermore, the optimal benefits obtained using the proposed optimization model are as suitable as the benefits achieved using the corresponding Monte Carlo-based optimization model.     

Ground water rights, spatial variation, and transboundary conflicts

The goal for any property rights system is to achieve equity, efficiency, and certainty. Trying to achieve these goals for ground water is difficult because a ground water right is not exclusive. To make matters more complicated, ground water is often under the jurisdiction of more than one political unit. The result is transboundary conflicts. Two critical elements must be included in any system of ground water rights. The system must define how the ground water can be used and define the relationships that each user and each use has with the other users and uses in the system. Unfortunately, these relationships are seldom completely defined and are made more complex by the transboundary scales at which they operate. As ground water moves horizontally across boundaries, different users or different jurisdictions have sequential control, creating conflicts between the first users and subsequent ones. Other problems occur because of vertical relationships, with more than one person or entity having control over ground water at the same time. This simultaneous exercise of authority can create conflicts between an individual who possesses a right to use ground water and a state or federal agency that regulates the same water. Transboundary conflicts occur at different scales and include conflicts between neighboring property owners as well as conflicts between countries. Scale, the property rights structure, and the nature of the relationship between users influence the way transboundary ground water conflicts are resolved.     

Restoration of wadi aquifers by artificial recharge with treated waste water

Fresh water resources within the Kingdom of Saudi Arabia are a rare and precious commodity that must be managed within a context of integrated water management. Wadi aquifers contain a high percentage of the naturally occurring fresh groundwater in the Kingdom. This resource is currently overused and has become depleted or contaminated at many locations. One resource that could be used to restore or enhance the fresh water resources within wadi aquifers is treated municipal waste water (reclaimed water). Each year about 80 percent of the country's treated municipal waste water is discharged to waste without any beneficial use. These discharges not only represent a lost water resource, but also create a number of adverse environmental impacts, such as damage to sensitive nearshore marine environments and creation of high-salinity interior surface water areas. An investigation of the hydrogeology of wadi aquifers in Saudi Arabia revealed that these aquifers can be used to develop aquifer recharge and recovery (ARR) systems that will be able to treat the impaired-quality water, store it until needed, and allow recovery of the water for transmittal to areas in demand. Full-engineered ARR systems can be designed at high capacities within wadi aquifer systems that can operate in concert with the natural role of wadis, while providing the required functions of additional treatment, storage and recovery of reclaimed water, while reducing the need to develop additional, energy-intensive desalination to meet new water supply demands.     

A visual basic spreadsheet macro for recession curve analysis

A Visual Basic program for an Excel spreadsheet was written to construct a master recession curve (MRC), using the adapted matching strip method, for recession analysis of ground water level time series. The program uses five different linear/nonlinear regression models to adjust individual recession segments to their proper positions in the MRC. The program can also be used to analyze the recession segments of other time series, such as daily stream discharge or stage. Some examples of field data from Croatia are used to illustrate the usefulness of its application.     

Modeling Aeration Efficiency of Stepped Cascades by Using ANFIS

The physical process of oxygen transfer or oxygen absorption from the atmosphere acts to replenish the used oxygen, a process termed re‐aeration or aeration. Aeration enhancement by macro‐roughness is well‐known in water treatment and one form is the aeration cascade. The macro‐roughness of the steps significantly reduces the flow velocities and leads to flow aeration along the stepped cascade. In this paper, the aeration efficiency in stepped cascade aerators was modeled by using the Adaptive Network Based Fuzzy Inference System (ANFIS). The obtained model was tested with experimental data. Test results showed that ANFIS can be used to estimate the aeration efficiency in stepped cascade aerators.     

Isolating the impacts of anthropogenic water use within the hydrological regime of north India

The effects of anthropogenic water use play a significant role in determining the hydrological cycle of north India. This paper explores anthropogenic impacts within the region's hydrological regime by explicitly including observed human water use behaviour, irrigation infrastructure and the natural environment in the CHANSE (Coupled Human And Natural Systems Environment) socio-hydrological modelling framework. The model is constrained by observed qualitative and quantitative information collected in the study area, along with climate and socio-economic variables from additional sources. Four separate scenarios, including business as usual (BAU, representing observed irrigation practices), groundwater irrigation only (where the influence of the canal network is removed), canal irrigation only (where all irrigation water is supplied by diverted surface water) and rainfed only (where all human interventions are removed) are used. Under BAU conditions the modelling framework closely matched observed groundwater levels. Following the removal of the canal network, which forces farmers to rely completely on groundwater for irrigation, water levels decrease, while under a canal-only scenario flooding occurs. Under the rainfed-only scenario, groundwater levels similar to current business-as-usual conditions are observed, despite much larger volumes of recharge and discharge entering and leaving the system under BAU practices. While groundwater abstraction alone may lead to aquifer depletion, the conjunctive use of surface and groundwater resources, which includes unintended contributions of canal leakage, create conditions similar to those where no human interventions are present. Here, the importance of suitable water management practices, in maintaining sustainable water resources, is shown. This may include augmenting groundwater resources through managed aquifer recharge and reducing the impacts on aquifer resources through occasional canal water use where possible. The importance of optimal water management practices that highlight trade-offs between environmental impact and human wellbeing are shown, providing useful information for policy makers, water managers and users. © 2019 John Wiley & Sons, Ltd.     

Simulating ground water-lake interactions: approaches and insights

Approaches for modeling lake-ground water interactions have evolved significantly from early simulations that used fixed lake stages specified as constant head to sophisticated LAK packages for MODFLOW. Although model input can be complex, the LAK package capabilities and output are superior to methods that rely on a fixed lake stage and compare well to other simple methods where lake stage can be calculated. Regardless of the approach, guidelines presented here for model grid size, location of three-dimensional flow, and extent of vertical capture can facilitate the construction of appropriately detailed models that simulate important lake-ground water interactions without adding unnecessary complexity. In addition to MODFLOW approaches, lake simulation has been formulated in terms of analytic elements. The analytic element lake package had acceptable agreement with a published LAKI problem, even though there were differences in the total lake conductance and number of layers used in the two models. The grid size used in the original LAKI problem, however, violated a grid size guideline presented in this paper. Grid sensitivity analyses demonstrated that an appreciable discrepancy in the distribution of stream and lake flux was related to the large grid size used in the original LAKI problem. This artifact is expected regardless of MODFLOW LAK package used. When the grid size was reduced, a finite-difference formulation approached the analytic element results. These insights and guidelines can help ensure that the proper lake simulation tool is being selected and applied.     

Trends in eutrophication research and control

Eutrophication is the natural ageing process of lakes. It is characterized by a geologically slow shift from in-lake biological production driven by allochthonous (external to the water body) loading of nutrients, to production driven by autochthonous (in-lake) processes. This shift typically is accompanied by changes in species and biotic community composition, as an aquatic ecosystem is ultimately transformed into a terrestrial biome. However, this typically slow process can be greatly accelerated by human intervention in the natural biogeochemical cycling of nutrients within a watershed; the resulting cultural eutrophication can create conditions inimical to the continued use of the water body for human-driven economic purposes. Excessive algal and rooted plant growth, degraded water quality, extensive deoxygenation of the bottom water layers and increased fish biomass accompanied by decreased harvest quality, are some features of this process. Following the Second World War, concern with cultural eutrophication achieved an intensity that spurred a significant research effort, culminating in the identification of phosphorus as the single most significant, and controllable, element involved in driving the eutrophication process. During the late 1960s and throughout the 1970s, much effort was devoted to reducing phosphorus in wastewater effluents, primarily in the developed countries of the temperate zone. These efforts generally resulted in the control of eutrophication in these countries, albeit with varying degrees of success. The present effort in the temperature zone, comprising mostly developed nations, has now shifted to the control of diffuse sources of a broader spectrum of contaminants that impact human water use. In the developing countries of the inter-tropical zone, however, rapidly expanding populations, a growing industrial economy and extensive urbanization have only recently reached an intensity at which cultural eutrophication can no longer be ignored. Further, initial attempts at applying temperate zone control measures in this region have been largely unsuccessful. Modification of the temperate zone eutrophication paradigm will be needed, especially to address the differing climatic and hydrological conditions, if cultural eutrophication is to be contained in this region, where eutrophication-related diseases continue to be a primary cause of human distress.     

Monitoring Continental Surface Waters by Satellite Altimetry

The monitoring of continental water stages is a requirement for meeting human needs and assessing ongoing climatic changes. However, regular gauging networks fail to provide the information needed for spatial coverage and timely delivery. Although the space missions discussed here were not primarily dedicated to hydrology, 18 years of satellite altimetry have furnished complementary data that can be used to create hydrological products, such as time series of stages, estimated discharges of rivers or volume change of lakes, river altitude profiles or leveling of in situ stations. Raw data still suffer uncertainties of one to several decimeters. These require specific reprocessing such as waveform retracking or geophysical correction editing; much work still remains to be done. Besides, measuring the flow velocity appears feasible owing to SAR interferometer techniques. Inundated surfaces, and the time variations of their extent, are currently almost routinely computed using satellite imagery. Thus, the compilation of the continuous efforts of the scientific community in these various investigative directions, such as recording from space the discharges of rivers or the change in water volume stored in lakes, can be foreseen in the near future.     

IHMS—Integrated Hydrological Modelling System. Part 1. Hydrological processes and general structure

A newly Integrated Hydrological Modelling System (IHMS) has been developed to study the impact of changes in climate, land use and water management on groundwater and seawater intrusion (SWI) into coastal areas. The system represents the combination of three models, which can, if required, be run separately. It has been designed to assess the combined impact of climate, land use and groundwater abstraction changes on river, drainage and groundwater flows, groundwater levels and, where appropriate, SWI. The approach is interdisciplinary and reflects an integrated water management approach. The system comprises three packages: the Distributed Catchment Scale Model (DiCaSM), MODFLOW (96 and 2000) and SWI models. In addition to estimating all water balance components, DiCaSM, produces the recharge data that are used as input to the groundwater flow model of the US Geological Survey, MODFLOW. The latter subsequently generates the head distribution and groundwater flows that are used as input to the SWI model, SWI. Thus, any changes in land use, rainfall, water management, abstraction, etc. at the surface are first handled by DiCaSM, then by MODFLOW and finally by the SWI. The three models operate at different spatial and temporal scales and a facility (interface utilities between models) to aggregate/disaggregate input/output data to meet a desired spatial and temporal scale was developed allowing smooth and easy communication between the three models. As MODFLOW and SWI are published and in the public domain, this article focuses on DiCaSM, the newly developed unsaturated zone DiCaSM and equally important the interfacing utilities between the three models. DiCaSM simulates a number of hydrological processes: rainfall interception, evapotranspiration, surface runoff, infiltration, soil water movement in the root zone, plant water uptake, crop growth, stream flow and groundwater recharge. Input requirements include distributed data sets of rainfall, land use, soil types and digital terrain; climate data input can be either distributed or non‐distributed. The model produces distributed and time series output of all water balance components including potential evapotranspiration, actual evapotranspiration, rainfall interception, infiltration, plant water uptake, transpiration, soil water content, soil moisture (SM) deficit, groundwater recharge rate, stream flow and surface runoff. This article focuses on details of the hydrological processes and the various equations used in DiCaSM, as well as the nature of the interface to the MODFLOW and SWI models. Furthermore, the results of preliminary tests of DiCaSM are reported; these include tests related to the ability of the model to predict the SM content of surface and subsurface soil layers, as well as groundwater levels. The latter demonstrates how the groundwater recharge calculated from DiCaSM can be used as input into the groundwater model MODFLOW using aggregation and disaggregation algorithms (built into the interface utility). SWI has also been run successfully with hypothetical examples and was able to reproduce the results of some of the original examples of Bakker and Schaars ( 2005 ). In the subsequent articles, the results of applications to different catchments will be reported. Copyright © 2010 John Wiley & Sons, Ltd.     

Impacts of droughts on low flows and water quality near power stations

ABSTRACT

Power plants often use river waters for cooling purposes and can be sensitive to droughts and low flows. Water quality is also a concern, due to algal blooms and sediment loads that might clog filters. We assessed the impacts of droughts on river flow and water quality from the point of view of power plant operation. The INCA (INtegrated CAtchment) water quality model was coupled with a climate model to create a dataset of flow and water quality time series, using the River Trent (UK) as a case study. The result hints to a significant decrease in flows and an increase in phosphorus concentrations, potentially enhancing algal production. Power plants should expect more stress in the future based on the results of this study, due to reduced cooling water availability and decreasing upstream water quality. This issue might have serious consequences also on the whole national power network.