INTEGRATING GEOGRAPHIC INFORMATION SYSTEMS AND MODFLOW FOR GROUND WATER RESOURCE ASSESSMENTS1

ABSTRACT: The Floridan Aquifer is the primary source of water in the coastal area of Santa Rosa County, Florida. In order to optimize well field design and analyze aquifer stress problems, the USGS MODFLOW code (McDonald and Harbaugh, 1988) is applied to develop a numerical computer model of the aquifer. The Geographical Information System (GIS) is the primary tool used in the development of the model grid, performance of the modeling procedure, and model analysis. The GIS is used in generating multiple grids in which to simulate both regional scale and local scale flow. The grid topology is recorded in geographic coordinates which facilitates geo-referencing and orientation of the grid to base maps and data coyerages. The GIS allows data transfer from various coverages to the nodes of the block centered grid where hydrogeologic information is stored as attributes to the grid coverage. From this grid coverage, pertinent information is queried within the GIS environment and used to generate the input files for the MODFLOW simulation. After MODFLOW execution, simulated heads and drawdown are imported into the grid coverage where residual error and recharge rates can be calculated. Contoured surfaces are then created for selected data sets including simulated heads, drawdown, residual error, and recharge rates. Model calibration is conducted utilizing the GIS to generate and process data sets associated with model simulations.     

RECYCLING INPUT DATA DURING ANALYTIC ELEMENT MODELING NEAR INDIANAPOLIS,INDIANA1

ABSTRACT: Overlapping and adjacent ground water investigations are common in areas where aquifers are threatened by industrial development. In the Indianapolis area in Marion County, Indiana, a patchwork of ground water flow models have been used during the past 20 years to evaluate ground water resources and to determine the effects of local contamination. In every case these ground water models were constructed from scratch. Site specific finite difference grids or finite element meshes inhibit the direct reuse of input data when the area of interest shifts. Because the aquifer is not discretized into a grid or mesh with analytic element models, there are unique opportunities for direct reuse of model input data. In two applications of this principle we illustrate how the newly emerging analytic element method allows a fairly straightforward reuse of model input data from previous models in the same general area. In analytic element models of Central Indiana, streams and their tributaries are represented in different resolutions. Input data items of several modeling studies are stored and cataloged on disk in such a manner that they can be selectively retrieved by a data management program PREPRO. In this manner, a new ground water model can be set up quickly with input data which have been previously defined and tested during model calibration.     

A RATING-CURVE METHOD FOR HYDRODYNAMIC SIMULATIONS IN CONTAMINANT TRANSPORT MODELING1

ABSTRACT: Accurate prediction of hydrodynamics is of great importance to modeling contaminant transport and water quality in a river. Flow conditions are needed in estimating potential exposure contamination levels and the recovery time for a no-action alternative in contaminated sediments remediation. Considering highly meandering characteristics of the Buffalo River, New York, a three-dimensional hydrodynamic model was selected to route upstream flows through the 8-km river section with limited existing information based on the model's fully predictive capability and process-oriented feature. The model was employed to simulate changes in water depth and flow velocity with space and time in response to variation in flow rate and/or water surface elevation at boundaries for given bottom morphometry and initial conditions. Flow conditions of the river reach where historical flow data are not available were computed. A rating-curve approach was developed to meet continuous and event contaminant modeling needs. Rating curves (depth-discharge and velocity-discharge relationships) were constructed at selected stations from the 3-D hydrodynamic simulations of individual flow events. The curves were obtained as steady solutions to an unsteady problem. The rating-curve approach serves to link flow information provided by the hydrodynamic model to a contaminant transport model. With the approach, the linking problem resulting from incompatible model dimensions and grid sizes can be solved. The curves will be used to simulate sediment movement and to predict contaminant fate and transport in the river.     

AN APPLICATION OF THE ANALYTIC ELEMENT METHOD IN MODELING FLORIDA EVERGLADES HYDROLOGY1

ABSTRACT: The large volumes of ground water that are discharged from the Everglades toward the Miami metropolitan area have historically posed a significant environmental water supply problem. In order to analyze the effects of seepage barriers on these subsurface outflows, the analytic element modeling code GFLOW was used to construct a ground water flow model of a region that includes a portion of the Everglades along with adjacent developed areas. The hydrology of this region can be characterized by a highly transmissive surficial aquifer in hydraulic contact with wetlands and canals. Calibration of the model to both wet and dry season conditions yielded satisfactory results, and it was concluded that the analytic element method is a suitable technique for modeling ground water flow in the Everglades environment. Finally, the model was used to evaluate the potential effectiveness of a subsurface barrier approximately two miles long for increasing water levels within the adjacent fringes of the Everglades National Park. It was found that the barrier had a negligible effect on water levels due to both its relatively short length and the high transmissivity of the surficial aquifer.     

Groundwater vulnerability assessment for the Banyas Catchment of the Syrian coastal area using GIS and the RISKE method

Vulnerability assessment to delineate areas that are more susceptible to contamination from anthropogenic sources has become an important element for sensible resource management and landuse planning. This contribution aims at estimating aquifer vulnerability by applying the RISKE model in Banyas Catchment Area (BCA), Tartous Prefecture, west Syria. An additional objective is to demonstrate the combined use of the RISKE model and a geographical information system (GIS) as an effective method for groundwater pollution risk assessment. The RISKE model uses five environmental parameters (Rock of aquifer media, Infiltration, Soil media, Karst, and Epikarst) to characterize the hydro-geological setting and evaluate aquifer vulnerability. The elevated eastern and low western part of the study area was dominated by high vulnerability classes, while the middle part was characterized by moderate vulnerability classes. Based on the vulnerability analysis, it was found that 2% and 39% of BCA is under low and high vulnerability to groundwater contamination, respectively, while more than 52% and 5% of the area of BCA can be designated as an area of moderate and very high vulnerability to groundwater contamination, respectively. The GIS technique has provided an efficient environment for analyses and high capabilities of handling a large amount of spatial data.     

AN INTERACTWE MODELING ENVIRONMENT FOR NON-POINT SOURCE POLLUTION CONTROL1

ABSTRACT: Non-point source pollution cuntinues to be an important environmental and water quality management problem. For the moat part, analysis of non-point source pollution in watersheds has depended on the use of distributed models to identify potential problem areas and to assess the effectiveness of alternative management practices. To effectively use these models for watershed water quality management, users depend on integrated geographic information systems (GIS)-based interfaces for input/output data management. However, existing interfaces are ad-hoc and the utility of GIS is limited to organization of input data and display of output data. A highly interactive water quality modeling interface that utilizes the functional components and analytical capability of GIS is highly desirable. This paper describes the tight coupling of the Agricultural Non-point Source (AGNPS) water quality model and ARC/INFO GIS software to provide an interactive hybrid modeling environment for evaluation of non-point source pollution in a watershed. The modeling environment is designed to generate AGNPS input parameters from user-specified GIS coverages, create AGNPS input data files, control AGNPS model simulations, and extract and organize AGNPS model output data for display. An example application involving the estimation of pesticide loading in a southern Iowa agricultural watershed demonstrates the capability of the modeling environment. Compared with traditional methods of watershed water quality modeling using the AGNPS model or other ad-hoc interfaces between a distributed model and GIS, the interactive modeling environment system is efficient and significantly reduces the task of watershed analysis using tightly coupled GIS databases and distributed models.     

Spatial and temporal modeling of microbial contaminants on grazing farmlands

This paper introduces an integrated spatial and temporal modeling system developed mathematically for assessing microbial contaminants on animal-grazed farmlands. The model uses fecal coliform, specifically Escherichia coli, as an indicator of fecal contamination and describes the sources, sinks, transport processes, and fate of E. coli contaminants in catchments and associated streams. Spatial features include grazing location, land topography, distance to a nearby stream, and distance through the stream network to the outlet. Temporal features are population dynamics on the land surface, in flow, and on streambeds. The model applies the principles of conservation of mass balance on two different types of pools: grid cells on land surfaces and networked stream segments. The model aims to improve the prediction of the effects of different land management strategies on the fecal contamination of waterways. This is achieved by characterizing the movement of fecal contaminants from land to streams and in-stream mobilization. Processes of attenuation, diffusion, and transport govern the movement. Our study site is a hill land catchment with an area of 140 ha and is used exclusively for animal grazing. The model was calibrated with previous research results, and then tested using the data collected at the outlet of the catchment. The sensitivity of the model predictions was analyzed for different scenarios: effect of stock rate, attenuation rate, and flow volumes. The similar pattern between monitored and predicted E. coli concentration proved that the model captures the key features that control the population dynamics of fecal contaminants. Further experiments are required to expand the model's functionality for covering more mitigation options.     

DYNAMIC RESPONSE OF AQUIFER SYSTEMS TO LOCALIZED RECHARGE1

ABSTRACT. The response of stream-unconfined aquifer systems to localized recharge is investigated by means of a two-dimensional finite element model. A variational approach is used in conjunction with the finite element method to solve the ground water flow equation. Linear approximated triangular elements are used to calculate the hydraulic head distribution in the flow region. The Crank-Nicholson centered scheme of numerical integration is employed to approximate the time derivative in the flow equation. A computer program is developed to calculate the hydraulic head distribution in the flow region. Solutions provided by the finite element model should prove useful in the evaluation of quantitative and qualitative changes in aquifer systems due to natural or artificial recharge. In addition, they should prove useful in the study of irrigation and drainage problems.     

USING SENSITIVITY ANALYSIS TO ASSIST PARAMETER ZONATION IN GROUND WATER FLOW MODEL1

ABSTRACT: For numerical modeling of ground water movement in a real aquifer system, the aquifer is usually divided into hydrogeologically defined zones, each with its own parameter values. The responses of the system, such as head or drawdown, are often available only in some of the zones. The estimated parameters of all the zones are based on the measured response in these limited zones. However, the estimates for some of the zones may be very uncertain, and these zones are therefore not justified by the data. In this paper, an approach is presented to understand which zone may produce uncertain parameter values and should be lumped with its neighbor. This approach is demonstrated using a regional numerical model for pumping test analysis in the Nottinghamshire aquifer, UK. A step-by-step process is used in identifying the aquifer zones and estimating their parameters based on the principle of using the smallest possible numbers of zones and parameters for adequate representation of the drawdown response. After the parameters of each zone are estimated, the sensitivity features of these parameters are examined. The results show that the parameters in one zone can be estimated properly by the drawdown in another zone only when there is significant sensitivity. For transmissivity, sensitivity between zones occurs when there is significant flow between them. For storativity, sufficient sensitivity can occur without large flows between the zones, provided that one zone causes significant drawdown in the other. This idea can be extended to the flow model for a large aquifer system. If the aquifer is divided in such a way that aquifer responses are not sensitive to the parameters in some of the zones, the parameters in those zones cannot be estimated properly and should be lumped into their neighboring zones. In this way, a simple but more reasonable model can be built.     

SPATIALLY DISTRIBUTED MODELING OF STREAM FLOW DURING STORM EVENTS1

ABSTRACT: The purpose of this paper is to present a new approach for the spatially distributed modeling of water flow during storm events. Distributed modeling of flow during storm events is an important basis for any environmental modeling, including turbidity or sediment transport. During the initial phase of a rainstorm, surface runoff is the main contributor of flow. To provide the spatial components for distributed hydrological modeling a Geographic Information System (GIS) was used to map and visualize contributing areas around a stream channel. Stream segments were defined using the hydrologic response unit (HRU) concept. Lateral flows were derived from GIS output for each segment of the stream and at each time interval of the rain storm and were routed using the kinematic routing equation. This approach is new in hydrological modeling and can be used to enhance many existing simulations. The model is also unique in the fine time scale (i.e., intervals are on the order of minutes). Model results showed good correlation with measured discharge values; however, further studies of contributing area behavior, its relationship with soil types and slope categories, and the influence of watershed size are needed to improve model performance. This model will be used in the future as the basis to model turbidity in streams.     

LOS ALAMOS FOREST FIRE IMPACT MODELING1

ABSTRACT: A controlled burn at Bandelier National Monument got out of control and burned about 43,000 acres (17,400 hectares) near Los Alamos, New Mexico, in May 2000. The wildfire caused dramatic changes in infiltration capacity and wettability of soils in many of the watersheds above Los Alamos National Laboratory (LANL) and destroyed the duff layer, dramatically reducing the interception and infiltration capacity of the formerly forested watersheds. These sudden changes in basin hydrology necessitated a rapid assessment of hydrology and hydraulics for the canyons running through LANL property to evaluate flood risk, plan emergency flood protection measures, and assess potential sediment and actinide transport. This paper presents the results of hydrologic and hydraulic modeling of Los Alamos Canyon following the wildfire. The large scale modeling effort, with over 13,000 cross sections for the hydraulic model (5,000 for Los Alamos Canyon, 8,000 for Guaje Canyon), relied heavily on a geographic information system (GIS) for model input and floodplain delineation. The HEC‐geoRAS model provided good integration between the hydraulic model (HEC‐RAS, Version 3.0.1) and the GIS (ArcView, v. 3.3). These modeling results are being used in drainage master planning efforts at LANL and in the development of sediment transport models using HEC‐6T. Sediment transport modeling results will be used to develop actinide transport models for the canyons at LANL.     

DEM AGGREGATION FOR WATERSHED MODELING1

ABSTRACT: The widely available USGS 7.5‐minute Digital Elevation Model (DEM) has a cell size of approximately 30 m × 30 m. This high resolution topographic information is impractical for many applications of distributed hydrologic and water quality models. In this study, cells were aggregated into coarse‐resolution areal units, termed grids, and a method to approximate flow direction for coarse‐resolution grids from 30 m DEM cells was developed. The method considers the flow path defined from the fine‐resolution DEM in determining a grid's flow direction and makes flow directions for grids closely follow the flow pattern suggested by the DEM. The aggregation method was applied to a DEM of Goodwater Creek, a nearly flat watershed that is located in central Missouri. The drainage networks derived for different levels of cell aggregations showed that grid aggregates of the Goodwater Creek watershed provided an adequate representation of the landscape topography.     

IDENTIFICATION OF AQUIFER DISPERSIVITIES: METHODS OF ANALYSIS AND PARAMETER UNCERTAINTY1

ABSTRACT: This paper presents a method for estimating aquifer dispersivities in solute transport models. Sensitivity equations are derived for the calculation of sensitivity coefficients. A modified Gauss-Newton algorithm is used to perform the least-squares minimization. A statistical procedure is outlined to assess reliability of the estimated parameters. The solute transport model is solved by the upstream weighted, multiple cell balance method which combines the concepts of local mass balance and finite element approximations. A one-dimensional solute transport problem in a vertical column system is first used to illustrate the inverse technique. A second example considers the parameter identification problem for three-dimensional solute transport with a unidirectional steady and uniform flow field. The third example solves the parameter identification problem in a three-dimensional, stream-aquifer, solute transport system with steady state flow. Numerical experiments are conducted to study data requirements for parameter identification.     

STORMFLOW SIMULATION USING A GEOGRAPHICAL INFORMATION SYSTEM WITH A DISTRIBUTED APPROACH1

ABSTRACT: With the increasing availability of digital and remotely sensed data such as land use, soil texture, and digital elevation models (DEMs), geographic information systems (GIS) have become an indispensable tool in preprocessing data sets for watershed hydrologic modeling and post processing simulation results. However, model inputs and outputs must be transferred between the model and the GIS. These transfers can be greatly simplified by incorporating the model itself into the GIS environment. To this end, a simple hydrologic model, which incorporates the curve number method of rainfall‐runoff partitioning, the ground‐water base‐flow routine, and the Muskingum flow routing procedure, was implemented on the GIS. The model interfaces directly with stream network, flow direction, and watershed boundary data generated using standard GIS terrain analysis tools; and while the model is running, various data layers may be viewed at each time step using the full display capabilities. The terrain analysis tools were first used to delineate the drainage basins and stream networks for the Susquehanna River. Then the model was used to simulate the hydrologic response of the Upper West Branch of the Susquehanna to two different storms. The simulated streamflow hydrographs compare well with the observed hydrographs at the basin outlet.     

A Stochastic Modeling Approach to Address Hydrogeologic Uncertainties in Modeling Wellhead Protection Boundaries in Karst Aquifers1

Abstract: Development of any numerical ground‐water model is dependent on hydrogeologic data describing the subsurface. These data are obtained from geologic core analyses, stratigraphic analyses, aquifer performance tests, and geophysical studies. But typically in remote areas, these types of data are very sparse and site‐specific in terms of the aerial extent of the resource to be modeled. Uncertainties exist as to how well the available data from a few locations defines a heterogeneous surficial aquifer such as the Biscayne Aquifer in Miami‐Dade County, Florida. This is particularly the case when an exceptionally conductive horizontal flow zone is detected at one site due to specialized testing that was not historically conducted at the other at sites that provided data for the model. Not adequately accounting for the potential effect of the high flow zone in the aquifer within a ground‐water numerical model, even though the zone may be of very limited thickness, might underpredict the well field protection capture boundaries. Applied Stochastic ground‐water modeling in determining well field protection zones is steadily becoming important in addressing the uncertainty of the hydrogeologic subsurface parameters, specifically in karstic heterogeneous aquifers. This is particularly important in addressing the uncertainty of a 60‐day travel time capture zone in the Northwest Well Field, Miami‐Dade County, where a predominantly high flow zone controls much of the flow in the production wells. A stochastic ground‐water modeling application along with combination of pilot points and regularization technique is presented to further consolidate the uncertainty of the subsurface.     

The Transboundary Nature of the Allende–Piedras Negras Aquifer Using a Numerical Model Approach

The Allende–Piedras Negras (APN) aquifer is located between the states of Texas (United States [U.S.]) and Coahuila (Mexico). The Rio Grande crosses the aquifer, acting as a natural and political divide between the countries. However, it remains unclear whether the APN aquifer can be considered a truly transboundary aquifer flow system, which would potentially require joint management by two different administrative jurisdictions. The main purpose of this study was to evaluate the transboundary nature of this aquifer. This was achieved by developing a detailed hydrogeological model to analyze the direction of volumetric fluxes within the APN aquifer using Visual MODFLOW. The model simulated a spatially averaged cumulative drawdown of 0.76 m for the entire aquifer over an 18‐year modeling period (2000–2017). The flow convergence zone, previously located below the Rio Grande, has shifted to the U.S. side in most locations, driven by higher pumping rates of the wells located near the river. This shift of the convergence zone from one country to the other means that groundwater recharge from one side flows underneath the river to the other side. This qualifies the APN aquifer as a “transboundary groundwater flow system.” The procedure followed in this study may be applied to other aquifers that straddle the U.S.–Mexico border and may motivate future modeling studies on other poorly studied transboundary aquifers around the world and thereby enable bi‐national aquifer management.     

Regional scale flood modeling using NEXRAD rainfall, GIS, and HEC-HMS/RAS: a case study for the San Antonio River Basin Summer 2002 storm event

This paper develops a framework for regional scale flood modeling that integrates NEXRAD Level III rainfall, GIS, and a hydrological model (HEC-HMS/RAS). The San Antonio River Basin (about 4000 square miles, 10,000 km2) in Central Texas, USA, is the domain of the study because it is a region subject to frequent occurrences of severe flash flooding. A major flood in the summer of 2002 is chosen as a case to examine the modeling framework. The model consists of a rainfall-runoff model (HEC-HMS) that converts precipitation excess to overland flow and channel runoff, as well as a hydraulic model (HEC-RAS) that models unsteady state flow through the river channel network based on the HEC-HMS-derived hydrographs. HEC-HMS is run on a 4 x 4 km grid in the domain, a resolution consistent with the resolution of NEXRAD rainfall taken from the local river authority. Watershed parameters are calibrated manually to produce a good simulation of discharge at 12 subbasins. With the calibrated discharge, HEC-RAS is capable of producing floodplain polygons that are comparable to the satellite imagery. The modeling framework presented in this study incorporates a portion of the recently developed GIS tool named Map to Map that has been created on a local scale and extends it to a regional scale. The results of this research will benefit future modeling efforts by providing a tool for hydrological forecasts of flooding on a regional scale. While designed for the San Antonio River Basin, this regional scale model may be used as a prototype for model applications in other areas of the country.     

SWATMOD‐Prep: Graphical User Interface for Preparing Coupled SWAT‐MODFLOW Simulations

This article presents SWATMOD‐Prep, a graphical user interface that couples a SWAT watershed model with a MODFLOW groundwater flow model. The interface is based on a recently published SWAT‐MODFLOW code that couples the models via mapping schemes. The spatial layout of SWATMOD‐Prep guides the user through the process of importing shape files (sub‐basins, hydrologic response units [HRUs], river network) from an existing SWAT model, creating a grid, performing necessary geo‐processing operations to link the models, writing out SWAT‐MODFLOW files, and running the simulation. The option of creating a new single‐layer MODFLOW model for near‐surface alluvial aquifers is available, with the user prompted to provide groundwater surface elevation (through a digital elevation model), aquifer thickness, and necessary aquifer parameter values. The option of simulating nitrate transport in the aquifer also is available, using the reactive transport model RT3D. The interface is in the public domain. It is programmed in Python, with various software packages used for geo‐processing operations (e.g., selection, intersection of rasters) and inputting/outputting data, and is written for Windows. The use of SWATMOD‐Prep is demonstrated for the Little River Experimental Watershed, Georgia. SWATMOD‐Prep and SWAT‐MODFLOW executables are available with an accompanying user's manual at: http://swat.tamu.edu/software/swat-modflow/ . The user's manual also accompanies this article as Supporting Information.     

Support of sustainable management of nitrogen contamination due to septic systems using numerical modeling methods

Nitrogen contamination is a serious concern to sustainable environmental management, and one important source of nitrogen contaminant is due to wastewater treatment using onsite sewage treatment and disposal systems (OSTDS, a.k.a., septic systems). This paper presents a study in which numerical modeling is used to support sustainable decision-making and management of nitrogen contamination by utilizing a recently developed GIS-based software, VZMOD, a Vadose Zone MODel for simulating nitrogen transformation and transport in vadose zone between drainfield of septic systems and water table. VZMOD is based on a physical model of unsaturated flow and nitrogen transformation and transport, and the model is solved numerically using the finite element methods. This is the major difference between VZMOD and other GIS-based software of nitrogen modeling. Using GIS techniques, VZMOD considers spatial variability of a number of hydrogeologic parameters such as hydraulic conductivity and porosity. A unique feature of VZMOD is that VZMOD addresses spatial variability of water table by using VZMOD together with ArcNLET, an ArcGIS-based software developed to simulate groundwater flow and nitrate load from septic systems to surface water bodies. VZMOD is designed to execute in different modes to be compatible with different levels of data availability in various management projects of nitrogen contamination. This paper presents an application of VZMOD at a neighborhood with hundreds of septic systems and heterogeneous hydraulic conductivity, porosity, and water table depth. The modeling results indicate that using septic systems at the considered neighborhood is unsustainable and more management means are necessary.